Identification and genomic analysis of two novel duck-origin GPV-related parvovirus in China
BMC Veterinary Research volume 15, Article number: 88 (2019)
Since early 2015, mule duck and Cherry Valley duck flocks have been suffering from short beak and dwarfism syndrome. This widely spreading infectious disease is characterized by growth retardation, smaller beak and tarsus with high morbidity and low mortality rate. For better understanding, we identified and characterized virus isolates named AH and GD from diseased Cherry Valley duck and mule duck flocks and investigated the damage caused by novel parvovirus-related virus (NGPV) to tissues and organs, including kidney, brain, pancreas, liver, spleen, bursa of fabricius and myocardial tissues.
AH and GD isolates shared high nucleotide identity with goose parvovirus (GPV). Alignment studies of AH and GD isolates showed 94.5–99.2% identity with novel parvovirus-related virus (NGPV), 98.7–91.5% identity with GPV and 79.9–83.7% with muscovy duck parvovirus (MDPV). Compared with other NGPV, classical GPV and MDPV sequences, a four 14-nucleotide-pair insertion in GD isolate was found in left open reading frame (ORF) (87–100 nt and 350–363 nt) and in right ORF (4847–4861 nt and 5122–5135 nt). However, in AH isolate, a five 14-nucleotide-pair deletions similar to other NGPV were found. The complete genome sequence comparison of eleven NGPV isolates from mule ducks and cherry valley ducks revealed no remarkable difference between them. Notably, the myocardium and bursa of fabricius of both disease and healthy animals are perfectly normal while other tissues have inflammatory cells exudation.
The AH and GD strains are novel parvovirus-related virus that isolates from mule ducks or cherry valley ducks which DNA sequence has no remarkable difference. The histopathology of tissues and organs such as kidney, brain etc. revealed non-significant changes in experimental and control animals. Overall, this study has contributed better understanding of molecular biology of NGPV strains and will help to develop the candidate strain for vaccine preparation to get better protection against these viral infections.
Muscovy duck parvovirus (MDPV) and Goose parvovirus (GPV) are members of the family Parvoviridae [1, 2]. GPV are etiological agents for Derzsy’s disease and results in huge economic losses to the duck industry . The mortality rate in goslings during first four weeks of life ranges from 70 to 100% . The disease in infected animals is characterized by lethargy, stunting, anorexia, locomotor dysfunction, watery diarrhea and death within 3 to 5 days after disease onset [1, 5, 6]. MDPV infection mainly affects three-week-old Muscovy ducklings and is characterized by signs almost similar to GPV. The complete nucleotide sequence of the two parvoviruses shares 79.7 to 85.0% homology.
Waterfowl parvoviruses, including MDPV and GPV, contain a linear, single-stranded DNA genome (approx. 5.0–5.1 kb in length) . At both ends of GPV and MDPV genome, there are two inverted terminal repeats (ITR) forming a hairpin structure which contains a stem region and a bubble region having length 456 and 442 nucleotides for MDPV and GPV respectively [8, 9]. There are two open reading frames (ORFs) in waterfowl parvoviruses [10, 11]. The left ORF encodes non-structural (NS) proteins NS1 and NS2 and right ORF encodes structural proteins namely VP1, VP2 and VP3 [12, 13].
During 1970s, Short beak and dwarfism syndrome (SBDS) was first reported in France in mule ducks, caused by a novel parvovirus-related virus (NGPV) [14, 15]. Since March 2015, there are numerous SBDS outbreaks in mule duck and Cherry Valley duck flocks in various regions of China [16, 17]. The morbidity rate of NGPV was found to be 10–30% while mortality rate was about 2% . The ducks infected with NGPV, showed clinical symptoms like swollen tongue, shorter tibia and stunted growth . In this way, the disease causes huge economic losses to duck industry by reduction in weight and size of animals. However, the adult ducks were found resistant to SBDS disease .
In order to have better understanding of SBDS disease, we characterized two NGPV namely AH and GD isolates, obtained from mule duck and Cherry Valley duck flocks in Anhui and Guangdong province, China. For detailed analysis, the complete genome of these isolates was sequencedand homology studies and phylogenetic comparison was carried out. The nucleotide sequences of the AH and GD genome and ITR were aligned for comparison with other GPV/MDPV/NGPV strains. The amino acid sequences of non-structural and structural proteins were also subjected to phylogenetic and homology analysis. The histopathological features were also studied in experimental and control groups.
In 2014–2015, two mule duck and Cherry Valley duck flocks located in Guangdong and Anhui province, China were found to be suffering from SBDS disease. Liver and spleen samples were collected from diseased birds and processed for viral isolation.
The AH and GD strains of variant GPV were isolated from mule duck and Cherry Valley duck flocks with typical signs of SBDS. Liver and spleen tissues were grinded and homogenized in physiological saline, freeze-thawed three times and centrifuged at 10,000 x g for 15 min at 4 °C. The supernatants were subjected to filtration using 0.22 μm sterile filters to remove contaminants and inoculated into chorioallantoic cavity of 10-day-old duck embryonated eggs (200 μl/egg). The allantoic fluid was harvested 5 days after inoculation and stored at − 80 °C for further testing.
The allantoic fulid of AH and GD isolates was centrifuged at 5000 xg for 10 min to remove cellular debris and proceeded for DNA extraction using AxyPrep DNA Extraction Kit (ANYGEN). The stock DNA was dissolved in TE buffer for future use.
Amplification of AH and GD complete genome
To amplify the complete genome of GD and AH isolates, six primer pairs were designed, based on the conserved regions in the genome sequence of NGPV (Table 1). The polymerase chain reaction conditions were as follows: 95 °C for 5 min, followed by 35 cycles of denaturation (95 °C for 30 s), annealing (52–60 °C for 30 s), extension (72 °Cfor 2 min) and a final extension at 72 °C for 10 min. The complete genome sequence was obtained by assembling all six DNA fragments.
Cloning and sequencing
The amplified PCR products were purified using gel extraction kit (AXYGEN), and all fragments were cloned into pMD18-T vector (TaKaRa). The plasmids were sequenced by Sangon Tech (Sangon Biotech, shanghai, China). For ensuring accuracy of genomic sequences, triplicate clones of fragments were sequenced.
Using BLAST sequence alignment, all six fragments were assembled to form a full–length complete genome. The sequence alignment and homology studies for whole-genome, NS (non–structural proteins) and VP (structural proteins) of GD and AH isolates was done by comparison with already reported sequences in the Genbank using MegAlign and MEGA 5.0 software programs  (Table 2).
Histopathological investigations were carried out in the diseased flocks by comparison with the healthy (Control) animals. At 14 days post infection (DPI), the animals were euthanized and necropsied. Intravenous anesthesia of high concentrations of pentobarbital sodium (150 mg/kg) was used in this study. Samples of kidney, brain, pancreas, liver, spleen, bursa of fabricius and myocardial tissues were fixed with 10% formalin; all the tissues were dehydrated, and immersed in transparent wax, cut into slices and stained with H & E dye. All of the animal protocols were performed in accordance with the ‘Guidelines for Experimental Animals’ of Ministry of Science and Technology (Beijing, China).
Two NGPV strains (GD and AH) were isolated from NGPV-positive mule duck and Cherry Valley duck samples respectively. These isolates were not able to cause death in duck embryonic eggs 5 days post infection (dpi). The main signs and symptoms in AH or GD-infected ducklings were retarded in growth with low mortality rate (Fig. 1).
Complete genome sequence analysis
The complete genome sequence of mule duck isolate “GD” and Cherry Valley duck isolate “AH” was 5106 bp and 5053 bp respectively. The sequences were submitted to GenBank and their accession numbers are indicated in Table 2. The NGPV genome contained two identical inverted terminal repeats (ITRs) and two open reading frames (ORFs). The length of ITRs of AH and GD genome was 379 and 407 nucleotides (nt) respectively.
Phylogenetic analysis of NGPV, GPV and MDPV sequences revealed that all NGPV isolates have more similarity with GPV than MDPV (Fig. 2). Homology among NGPV isolates ranges from 93.4 to 99.9% at whole genome level. All NGPV isolates share 89.7–96.7% identity with classical GPV and 79.6–83.9% identity with MDPV. These results revealed that there is no significant difference among NGPV isolated either from mule duck or Cherry Valley duck and also the NGPV isolates are more closely related to GPV virus. Sequence alignment of NGPV with classical GPV and MDPV using MEGA 5.0 software revealed three 14-nucleotide sequences deletion in most of NGPV isolates (Fig. 3). A putative recombination was also observed in GD strain between 87 nt – 100 nt (Fig. 3).
Phylogenetic analysis of NS genes
In order to understand the genetic diversity of Non-structural proteins (NS) in NGPV, GPV and MDPV isolates, phylogenetic tree was constructed and homology analysis was also carried out (Fig. 4). Through sequence alignment, it was found that the length of NS proteins of all NGPV, GPV and MDPV sequences is 1881 nt, and they encode 626 amino acids. Through comparison of NS gene, both NGPV isolates (AH and GD) have more similarly with GPV than MDPV, which is similar to complete genome sequence findings.
Homology of NS proteins among NGPV isolates ranged between 93.7–100% at gene level and 95.4–100% at amino acid level. However, they share 93.4–97.9% identity with classical GPV and 81.0–82.4% identity with MDPV. At amino acid level, all NGPV isolates shared 96.3–99.5% identity with classical GPV and 89.0–90.1% identity with MDPV. As compared to classical GPV and MDPV, there are twelve synchronous specific mutations in NS amino acids (Table 3).
Phylogenetic analysis of VP genes
The genetic diversity of VP gene (Capsid protein gene) was also studied by phylogeny and homology analysis (Fig. 5). Through sequence alignment, it was found that the length of all VP genes of NGPV, GPV and MDPV is 2199 nt and they encode 732 amino acids. The comparison of VP gene also showed more similarity of NGPV (AH and GD) sequences with GPV than MDPV sequences, which is similar to the findings of complete genome sequence and NS gene analysis.
Homology among NGPV isolates ranged between 94.0–100% at VP gene level and 96.9–100% at amino acid level. The VP gene of NGPV isolates has 90.9–97.5% identity with classical GPV and 80.9–91.5% identity with MDPV. At VP amino acid level, NGPV isolates shared 95.1–98.2% identity with classical GPV, and 88.0–92.6% identity with MDPV. As compared to classical GPV and MDPV, there are eight synchronous specific mutations in VP amino acids (Table 4).
To find out the damage caused by NGPV to tissues and organs, including kidney, brain, pancreas, liver, spleen, bursa of fabricius and myocardial tissues, histopathological investigations were carried out (Fig. 6). In both experimental and control animals, myocardium and bursa of fabricius are perfectly normal and the other tissues have inflammatory cells exudation. The results suggested that NGPV has no significant influence on these tissues and organs.
In this study, we studied the complete genome sequence of AH and GD, isolated from mule duck and Cherry Valley duck flocks reported with SBDS from Anhui and Guangdong province, China. SBDS in ducks was first reported in France in 1970s , followed by outbreaks in Taiwan (1989), Poland (1995) [2, 14], Hungary (2009)  and mainland China (2014) . Although the genome sequences of NGPV, GPV and MDPV strains have been reported previously, but this is the first study, which reported variations among various parvovirus duck strains on the basis of phylogeny and homology analysis [4, 14, 17, 21]. The sequence analysis indicated that all NGPV isolates have more homology among them as compared to classical GPV and MDPV. Several studies reported SBDS in ducks maybe caused by novel-MDPV is also reported, but the whole genome sequence analysis has shown that GPV clade also exists in this paper, which the rusults might be due to coinfection with NGPV and MDPV strains .
The complete genome sequence of AH and GD strains contain 5053 nt and 5106 nt respectively and they shared 93.4–99.9% identity with other NGPV, 89.7–96.7% identity with classical GPV and 79.6–83.9% identity with MDPV strains. The length of ITRs of AH and GD is 379 and 407 nt respectively. Through sequence alignment, it was found that the length of NS protein of all NGPV, GPV and MDPV strains is 1881 nt encoding 626 amino acids and the length of VP gene of all NGPV, GPV and MDPV strains is 2199 nt which encodes 732 amino acids. When comparison was carried out between mule duck isolates including GD and Cherry Valley duck isolates including AH, it was found that NGPV isolated from mule duck and those isolated from Cherry Valley duck have no significant differences among them.
VP protein plays an important role in the stimulation of protective response in the body against the virus as it contains viral antigenic sites [23,24,25]. Interestingly VP protein of NGPV and classical GPV shares 90.9 to 97.5% homology, therefore it is very likely that vaccine against GPV is effective against NGPV. However, in this case, it’s protective efficacy may slightly decrease. Further studies are required to explore the effectiveness of currently available vaccines against NGPV and to develop a new vaccine for NGPV. It will be a difficult task as virus titer is low in duck embryo and large scale culturing of NGPV would be challenging.
The complete genome sequence analysis also showed that most of NGPV isolates have several 14-nucleotide sequence deletions not found in classical GPV and MDPV sequences. It was also noticed and reported for the first time that the 14-nucleotide sequences deletion is identical to each other and the front and back of deleted position also have a same 14-nucleotide sequences. These results lead us to suspect that the 14-nucleotide sequences deletion play a major role in migration of host, if sequencing mistakes are excluded. The SimPlot analysis of GD strain with other strains revealed no clear indication of increase in genome at at 87 nt − 100 nt. However, as it is too short (14 nt), therefore it might be possible that SimPlot software could not able to use analysis recombination.
The histopathological examination revealed no remarkable difference in kidney, brain, pancreas, liver, spleen, bursa of fabricius, myocardium tissues between experimental group and control group. These findings suggested that NGPV causes non-significant pathological changes in these organs and tissues, which might be the reason of low mortality in infected animals (2%).
In conclusion, the novel parvovirus-related virus is responsible for SBDS disease outbreaks in eastern China since 2014 [21, 26]. The low mortality and smaller beak is the major difference between NGPV and GPV infected animals. NGPV is hard to culture in duck embryos and DEF (duck embryos fibroblast). Since classical GPV does not infect mule duck and Cherry Valley ducks, therefore NGPV is not completely adapted to cause severe disease in mule duck and Cherry Valley duck which might be the reason for low mortality caused by these strains. The NGPV also has little pathological impact on kidney, brain and other organs and tissues. However, further studies are required to evaluate the protective efficacy of classical GPV vaccine against SBDS and to develop a new vaccine for this disease.
AH and GD strains isolated from mule duck and Cherry Valley duck flocks are novel parvovirus-related virus. The comparison of complete genome sequence of eleven NGPV isolates from mule ducks and cherry valley ducks found no remarkable difference between them. Through histopathological examination, it was found that both experimental and control groups have perfectly normal myocardium and bursa of fabricius while other tissues have inflammatory cells exudation. The results suggested that NGPV has no significant influence on the tissues and organs. This study has contributed better understanding of molecular biology of NGPV strains and will help to develop the candidate strain for vaccine preparation in future.
Inverted terminal repeats
Muscovy duck parvovirus
Novel parvovirus-related virus
Open reading frame
Short beak and dwarfism syndrome
Capsid protein gene
Glávits R, Zolnai A, Szabó E, Ivanics E, Zarka P, Mató T, Palya V. Comparative pathological studies on domestic geese (Anser anser domestica) and Muscovy ducks (Cairina moschata) experimentally infected with parvovirus strains of goose and Muscovy duck origin. Acta Vet Hung. 2005;53(1):73.
Woźniakowski G, Kozdruń W, Samoreksalamonowicz E. Genetic variance of Derzsy's disease strains isolated in Poland. Journal of Molecular & Genetic Medicine An International Journal of Biomedical Research. 2009;3(2):210.
Derzsy D, Kisary J. Viral disease of goslings: World Veterinary Congress; 1975.
Yu K, Ma X, Sheng Z, Qi L, Liu C, Dan W, Bing H, Feng L, Song M. Identification of goose-origin parvovirus as a cause of newly emerging beak atrophy and dwarfism syndrome in ducklings. J Clin Microbiol. 2016;54(8):1999.
Kisary J: Some growth characteristics of goose parvovirus strain "B". Acta Vet Acad Sci Hung 1974, 24(3):329.
Takehara K, Hyakutake K, Imamura T, Mutoh KI, Yoshimura M. Isolation, identification, and plaque titration of parvovirus from Muscovy ducks in Japan. Avian Dis. 1994;38(4):810.
Le GRG, Jestin V. Biochemical and genomic characterization of Muscovy duck parvovirus. Arch Virol. 1994;139(1–2):121–31.
Wang J, Duan J, Meng X, Gong J, Jiang Z, Zhu G. Cloning of the genome of a goose parvovirus vaccine strain SYG61v and rescue of infectious virions from recombinant plasmid in embryonated goose eggs. J Virol Methods. 2014;200(10):41.
Wang J, Huang Y, Zhou M, Hardwidge PR, Zhu G. Construction and sequencing of an infectious clone of the goose embryo-adapted Muscovy duck parvovirus vaccine strain FZ91-30. Virol J. 2016;13(1):1–8.
Ji J, Xie QM, Chen CY, Bai SW, Zou LS, Zuo KJ, Cao YC, Xue CY, Ma JY, Bi YZ. Molecular detection of Muscovy duck parvovirus by loop-mediated isothermal amplification assay. Poult Sci. 2010;89(3):477–83.
Shen H, Zhang W, Wang H, Zhou Y, Shao S. Identification of recombination between Muscovy duck parvovirus and goose parvovirus structural protein genes. Arch Virol. 2015;160(10):1–5.
Zádori Z, Stefancsik R, Rauch T, Kisary J. Analysis of the complete nucleotide sequences of goose and Muscovy duck Pervoviruses indicates common ancestral origin with adeno-associated virus 2. Virology. 1995;212(2):562–73.
Poonia B, Dunn PA, Lu H, Jarosinski KW, Schat KA. Isolation and molecular characterization of a new Muscovy duck parvovirus from Muscovy ducks in the USA. Avian Pathology. 2006;35(6):435–41.
Palya V, Zolnai A, Benyeda Z, Kovács E, Kardi V, Mató T. Short beak and dwarfism syndrome of mule duck is caused by a distinct lineage of goose parvovirus. Avian Pathology. 2009;38(2):175.
Villate D, Caldier P. Manuel pratique des maladies des palmipedes; 1989.
Chen H, Dou Y, Tang Y, Zhang Z, Zheng X, Niu X, Yang J, Yu X, Diao Y. Isolation and genomic characterization of a duck-origin GPV-related parvovirus from Cherry Valley ducklings in China. PLoS One. 2015;10(10):e0140284.
Li P, Lin S, Zhang R, Chen J, Sun D, Lan J, Song S, Xie Z, Jiang S. Isolation and characterization of novel goose parvovirus-related virus reveal the evolution of waterfowl parvovirus. Transboundary & Emerging Diseases. 2017;65(2).
Fan W, Sun Z, Shen T, Xu D, Huang K, Zhou J, Song S, Yan L. Analysis of evolutionary processes of species jump in waterfowl parvovirus. Front Microbiol. 2017;8:e1003537.
Tamura K. MEGA5: molecular evolutionary genetics analysis version 6.0. Molecular Biology & Evolution. 2013;30(12):2725–9.
Chen S, Wang S, Cheng X, Xiao S, Zhu X, Lin F, Wu N, Wang J, Huang M, Zheng M. Isolation and characterization of a distinct duck-origin goose parvovirus causing an outbreak of duckling short beak and dwarfism syndrome in China. Arch Virol. 2016;161(9):2407–16.
Xiao S, Chen S, Cheng X, Lin F, Wang S, Zhu X, Yu B, Huang M, Wang J, Wu N. The newly emerging duck-origin goose parvovirus in China exhibits a wide range of pathogenicity to main domesticated waterfowl. Vet Microbiol. 2017;203:252–6.
Fu Q, Huang Y, Wan C, Fu G, Qi B, Cheng L, Shi S, Chen H, Liu R, Chen Z. Genomic and pathogenic analysis of a Muscovy duck parvovirus strain causing short beak and dwarfism syndrome without tongue protrusion. Res Vet Sci. 2017;115:393.
Berns KI. Parvovirus replication. Microbiol Rev. 1990;54(3):316–29.
Salganik M, Aydemir F, Nam HJ, Mckenna R, Agbandjemckenna M, Muzyczka N. Adeno-associated virus capsid proteins may play a role in transcription and second-Strand synthesis of recombinant genomes. J Virol. 2014;88(2):1071–9.
Hueffer K, Parker JSL, Weichert WS, Geisel RE, Sgro JY, Parrish CR. The natural host range shift and subsequent evolution of canine parvovirus resulted from virus-specific binding to the canine transferrin receptor. J Virol. 2003;77(3):1718–26.
Ning K, Wang M, Qu S, Lv J, Yang L, Zhang D. Pathogenicity of Pekin duck- and goose-origin parvoviruses in Pekin ducklings. Vet Microbiol. 2017;210:17–23.
This was supported by National Program on Key Reasearch Project of China (No. 2017YFD0500800), Guangdong Provincial Scientific and Technological Planning Project (No. 201704020083), Guangdong Provincial special Project for Science and Technology development (2016B020234006,2017B090904029), Guangdong Province Technology development and industrialization (2015B020203006).
Availability of data and materials
The genomic sequence has been deposited in GenBank with accession number MH444513 and MH444514. Other data supporting the findings is contained within the manuscript.
Ethics approval and consent to participate
All of the animal protocols were performed in accordance with the ‘Guidelines for Experimental Animals’ of Ministry of Science and Technology (Beijing, China) and the study was approved by the Research Ethics Committee of the Institute of Veterinary Medicine College of South China Agricultural University. No specific permits were required for these locations and activities.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Bian, G., Ma, H., Luo, M. et al. Identification and genomic analysis of two novel duck-origin GPV-related parvovirus in China. BMC Vet Res 15, 88 (2019). https://doi.org/10.1186/s12917-019-1833-9
- Novel parvovirus-related virus
- Phylogenetic analysis
- Mule duck
- Cherry valley duck