Skip to main content

Comparative genetic analysis and pathological characteristics of goose parvovirus isolated in Heilongjiang, China



Goose parvovirus (GPV) causes acute enteritis, hepatitis, myocarditis and high morbidity and mortality in geese and ducks. GPV H strain was isolated from a Heilongjiang goose farm where the geese were showing signs of hemorrhage in the brain, liver, and intestinal tract. In this study, we explored the genetic diversity among waterfowl parvovirus isolates and the pathological characteristics of GPV H in Shaoxing ducklings.


The complete capsid protein (VP) and non-structural (NS) sequences of the isolated H strain were sequenced, and phylogenetic trees of VP and NS were constructed in MEGA version 5.05 using the neighbor-joining method. Three-day-old Shaoxing ducklings were inoculated with GPV and were euthanized at 1, 2, 4, 6, and 8 days post-inoculation (PI), and their organs were removed and collected. The organs of 6-day PI ducklings were fixed in formalin, embedded in paraffin, sectioned for histology, stained with HE and analyzed for pathological lesions. The distribution of the GPV H strain in the tissues of the inoculated ducklings was detected using the polymerase chain reaction (PCR) method.


Genetic analysis of the NS and VP genes indicated that the H strain was closely related to strains circulating in China during 1999–2014, and the nucleic acid identity of those strains was 98%–99%. Classical symptoms were observed in the inoculated ducklings. GPV remained in many tissues and replicated in a majority of the tissues, leading to histopathological lesions in four tissues.


We first reported the distribution and histopathological lesions of a Chinese strain of GPV in infected shaoxing ducklings. This H strain was moderate pathogenic for Shaoxing ducklings.


GPV causes watery diarrhea, weight loss, and stunting in ducklings and goslings, with high morbidity and mortality, and this disease is known as 3-week disease or Derzsy’s disease [1, 2]. GPV is a small virus particle, and its genome is composed of a linear, single-stranded DNA, approximately 5.1 kb in length [3, 4]. The genome contains two main open reading frames (ORFs), with the left ORF encoding the NS protein and the right ORF encoding the structural proteins VP1, VP2 and VP3 [5, 6].

GPV circulated in the Fujian, Jiangsu, Heilongjiang, Shanghai, Guangdong, Sichuan and Shandong Provinces in China during 1999–2014 [7, 8]. Many isolated GPV strains were pathogenic and lethal to goslings and ducklings. Since 2015, many Cherry Valley ducklings flocks appeared short beak and dwarfism syndrome in Northern China, including Jiangsu, Anhui, and Shanghai Provinces, and this disease was caused by a distinct GPV-related virus, namely novel goose parvovirus (N-GPV) [9, 10]. Many phylogenetic analysis have indicated that the GPV strains of mainland China can be divided into three distinct groups, PRC-I, PRC-II and PRC-III, and N-GPV was closely related to these GPV strains [11]. There are several studies regarding the distribution of virulent GPV and attenuated GPV in Muscovy ducklings, detected by using the polymerase chain reaction (PCR) [12, 13]. However, the distribution of recent prevalent Chinese strains of GPV in Shaoxing ducklings has not previously been reported.

Wang Bin and his colleagues isolated the H strain from Heilongjiang, sequenced and analyzed its VP3 sequence, and speculated that it was a virulent strain [14]. In this study, we sequenced and analyzed the NS and VP genes of the H strain. We also explored the pathogenicity of the H strain in 3-day-old Shaoxing ducklings by analyzing clinical symptoms, pathological changes, and the distribution of GPV in the tissues of ducklings inoculated with the virus.


Virus and cells

H strain GPV allantoic fluid was isolated by inoculation of goose embryos from goose farms of the Heilongjiang Province [14]. Duck embryo fibroblasts (DEF) cells were prepared according to the standard procedure and cultured in DMEM supplemented with 5% fetal bovine serum and 1% antibiotics (Gibco) at 5% CO2 37 °C. The primary monolayer DEFs was washed two times by PBS until up to 80–90% and the H strain allantoic fluid was inoculated into the primary monolayer DEFs for 1 h, then removed the allantoic fluid and cultured in the basic DMEM for 5 days.

Viral DNA extraction and genome sequencing

Viruses were propagated in DEFs, and culture supernatants were harvested at 5 days post-inoculation. Viral DNA was extracted using viral DNA kit (Omega, USA) according to the manufacturer’s instructions. Briefly, added 250 μl cell supernatants, 10 μl OB protease, 250 μl Buffer BL and 4 μl Liner Acrylamide to a sterile microcentrifuge tube, mixed thoroughly and incubated at 65 °C for 10 min, then added 260 μl absolute ethanol to the above lysate, transfered the mixed lysate into an HIBind DNA Mini column which assembled in a 2 ml collection tube, and centrifuged at 8000 g for 1 min. Then washed the HIBind DNA Mini column by pipeting 500 μl buffer HB, 700 μl DNA wash buffer and repeated 700 μl DNA wash buffer, successively. Dryed the column by centrifuging at 15000 g for 2 min, placed the column into a sterile 1.5 ml microfuge tube, added 50 μl preheated (65 °C) elution buffer and sit for 5 min at room temperature. Centrifuged at 8000 g for 1 min and stored the elute DNA at − 20 °C. The primers complementary to NS (536–2419) and VP3 (3032–4636) encoding regions were designed on the basis of the complete sequences of the GPV E strain (GenBank accession No. KC184133.1). The primer sequences were as follows: NS-F, ATGGCACTTTCTATG CCTCTTCTGA; NS-R, TTATTGTTCATTTTCAGCATCATCA and VP3-F, ATGGCAGAGGGAGGAGGC GGAGCTT; VP3-R, TTACAGATTTTGAGTTAGATATCTG, respectively. The design of specific primers for the fragment 1106 bp (1955–3100) containing C terminal of NS (1955–2419) and N terminal of VP sequence (2438–3100) were based on the NS and VP3 sequences of H strain, the primers were as follows: GPV-2-F, TGGGCCAATGATAATCTAGTTCCTGTTG and GPV-2-R, CGAGGCATTACCC ACTCCATCGGCACC. The PCR conditions of three primers consisted one cycle of 2 min at 94 °C, 30 cycles of 30s at 98 °C, 30s at 55°Cand 2 min at 72 °C, followed by 10 min at 72 °C.The three PCR products were purified by DNA gel extraction kit (Axygen, America) and cloned into the pMD18T vector. Then we selected four positive clones of each product and sequenced (Bo Shi Biotechnology Company, China).

Sequence and phylogenetic analyses of the H strain

The GPV-2 and VP3 fragments were assembled into the capsid protein VP gene sequence. The NS and VP sequences were aligned by clustalw. MEGA5 software was used to analyze sequences, and construct phylogenetic trees using the neighbor-joining method based on 1000 bootstrap replications and to determine phylogenetic distances [15, 16]. I, II, III and IV groups were sorted by 95%–100% nucleotide identity, and Ia, Ib, Ic, Id, Ie, IIIf, IIIh, and IIIg group divisions were based on 98%–100% nucleotide identity. The VP sequences of 38 strains were redivided into groups (I and II), the cut off value for condensed tree was 50% and the pair-wise distance cut off for grouping them into separate groups was default value for MEGA 5.05.

Assessing the pathogenicity of the H strain in Shaoxing ducklings

In order to identify the appropriate infectious dose, 23 11-day-old specific pathogen-free duck embryo eggs were divided into A(5), B(5), C(5), D(5) and control groups(3). A, B, C, D and control groups were inoculated with the serial dilutions 1.2 × 102 plaque forming units (PFU)/0.2 mL, 1.2 × 103PFU, 1.2 × 104PFU, 1.2 × 105PFU, and 0.2 mL PBS per duck embryo egg, respectively. The symptoms of duck embryo eggs were observed for 7 days post inoculation. The 15 3-day-old Shaoxing ducklings were divided into two groups (I and II). Ten ducklings from Group I were injected intramuscularly with 1.2 × 103 PFU/0.2 mL GPV cell supernatant per duckling and were observed for 8 days PI. The five remaining ducklings were used as negative controls and were inoculated intramuscularly with 0.2 mL PBS. Two ducklings from Group I were sacrificed per time point at 1, 2, 4, 6, or 8 days PI, and their organs were removed and collected, respectively. At the same time, one duckling from the negative control group was euthanized and necropsied the same way. The organs of 6-day PI ducklings were fixed in formalin, embedded in paraffin, sectioned, stained with HE and analyzed histopathologically.

Distribution of the GPV H strain in Shaoxing ducklings

The duodenum, jejunum, ileum, rectum, cardiac muscle, liver, lung, kidney, pancreas, bursa of Fabricius, skeletal muscles and brain of sacrificed ducklings were collected. 100 mg tissues and organs were suspended evenly in 1 mL PBS buffer using a tissue homogenizer instrument. Tissue and organ fluid was frozen and thawed three times before centrifuging at 12000×g for 10 min. DNA from 250 μL supernatant was extracted as described above. The primers GPV-F, ATGGCAGAGGGAGGAG GCGG and GPV-R, CGCATGGTGCCTTCCGT were designed and used to detect GPV DNA.


Phylogenetic analysis of the VP and NS sequences of the H strain

We obtained the VP sequence that was 2199 nucleotides in length and encoded 731 amino acids (GenBank accession No. KY411078). Alignment the VP sequence of the H strain with 37 other strains indicated that the VP of the H strain shares 90%–99% nucleotide identity with the 37 other GPV strains. These 38 GPV strains were divided into two groups (I and II), with group I containing five subgroups (Ia, Ib, Ic Id and Ie) (Fig. 1). Group Ia contained the H strain, 5 strains isolated from Taiwan, 1 strain isolated from Shanghai, 5 strains from Jiangsu, and 6 strains from Anhui. The VP sequences of these strains shared 98%–99% nucleotide identity; in particular, the H strain shared 99.77% nucleotide identity with the E strain, and there were four nucleotides differenerece in the VP sequence between the H strain and the E strain. Thirty-nine and 2166 nuclei acid were synonymous, however, nucleotides 149A and 1567A in the E strain were mutated into 149C and 1567G in the H strain, respectively. This corresponded to amino acids (aa) mutation T50 N and K523E in the H strain. Group Ib contained 1 strain from Shanghai, 3 vaccine strains, 1 strain from Guangdong, 1 strain from Chongqing and 1 strain from Heilongjiang (Fig. 1). The VP sequences of Ib strains shared 95%–96% nucleotide identity with the H strain, and 523aa was also an E. Group Ic contained 3 European strains and 1 strain from Northeastern China, and these strains shared 96%–97% nucleotide identity with the H strain. In particular, the H strain shared 96.63% identity with the B strain; 73 nucleotides in the H strain varied from the B strain, resulting in 17 aa differences among the H strain and the B strain (Table 1). 521aa, 529aa, 558aa, and 615aa in the VP sequence of the B strain were tryosine (Tyr, Y), leucine (Leu, L), glutamic acid (Glu, E), and glycine (Gly, G), respectively, however, the same positions in the H strain were histidine (His, H), isoleucine (Ile, I), aspartic acid (Asp, D), and tryptophan (Trp, W). 35aa, 41aa, 207aa and 210aa in the VP sequence of the B, VG32/1, QH15, AH1, M15, DB3, PT and D strains were serine (Ser, S), aspartic acid (Asp, D), methionine (Met, M), and Ser, respectively. However, these VP amino acids in other GPV strains were threonine (Thr, T), alanine (Ala, A), Leu and Ala (Fig. 1). Group II contained 3 strains isolated from ducks that shared 90%–92% identity with the H strain (Fig. 1).

Fig. 1

Phylogenetic tree based on the VP sequences of 38 GPV strains constructed in MEGA version 5 using the neighbor-joining method. The round shape represents the H strain

Table 1 Amino acid differences in the VP sequences of the H and B strains of GPV

The NS sequence from the H strain was 1884 nucleotides in length (GenBank accession No. KY411079), and the phylogenetic relationships among NS sequences is shown in Fig. 2. The NS sequences of 28 strains shared 83%–99% nucleotide identity. They were divided into 2 groups (III and IV), and Group III was further divided into three subgroups (IIIf, IIIg and IIIh) (Fig. 2). Group IIIf contained 1 standard virulent B strain and 1 European vaccine strain, 2 strains isolated from Taiwan, 4 strains from Jiangsu, 1 Heibei strain, 1 Shanghai strain, 1 Chongqing strain, 1 Heilongjiang strain and 3 strains from Anhui, and these strains shared 98%–99% nucleotide identity (Fig. 2). The variations in the NS sequences included 5 nucleotide and 3 amino acid differences between the H strain and the E strain. However, there were 16 nucleotide and 4 amino acid differences between the H strain and the B strain (Table 2). Amino acids 149H and 343I were only found in the H strain, while the other strains contained Q and T, respectively. Group IIIg also contained 3 strains isolated from ducks, 1 Taiwan strain, and 1 European strain. These strains shared 96%–98% nucleotide identity with the H strain. Group IIIfh contained 2 vaccine strains, 1 Heilongjiang strain, 1 Shanghai strain and 1 Liaoning strain (Fig. 2). Group IV contained 2 GPV strains isolated from Muscovy ducks, those strains nucleotide identity with the H strain was 83% (Fig. 2).

Fig. 2

Phylogenetic tree based on the NS sequences of 28 GPV strains constructed in MEGA version 5 using the neighbor-joining method. The round shape represents the H strain

Table 2 Amino acid differences in the NS sequences of the H and B strains of GPV

The pathogenicity of the H strain in duck embryo

No duck embryos of five groups died during 7 days observation period and were gross lesions at necropsy. Duck embryo muscle tissues were detected by PCR. The 400 bp fragment can be detected in duck embryo muscle of A, B, C and D groups and the sequence of 400 bp fragment was identical to the sequence of GPV. However, the 400 bp fragment in duck embryo muscle tissues inoculated 1.2 × 103 and 1.2 × 104PFU GPV are brighter than other dilutions (Additional file 1: Figure S2). The result suggested that the 1.2 × 103 and 1.2 × 104PFU GPV replicate rapidly in duck embryo, therefore 1.2 × 103PFU could be used as the appropriate infectious dose for investigating the pathogenicty of GPV H strain.

The pathogenicity of the H strain in ducklings

There was 90% morbidity in the Group I ducklings. No ducklings died during 1–8 days PI. The ducklings were depressed, had diarrhea and lost their appetite at 24 h PI. The clinical symptoms were most obvious at 4–5 days PI; however, these symptoms decreased at 6–8 days PI. Two ducklings from Group I were necropsied at 24, 48, 96, 144, and 192 h PI, and their organs were removed and collected. The gallbladders of inoculated ducks were swollen and full of black-green fluid at 24, 48, 96, and 144 h PI, respectively. The terminus of the cecum swelled in one of two necropsied ducklings at 24, 96, and 192 h PI, however, the entire cecum of other ducklings was swollen compared to the control ducklings (Fig. 3). The thymus of inoculated ducklings showed bleed spots at 48, 96, and 144 h PI (Fig. 3). The intestinal tracts of inoculated ducks were swollen at 48, 96, 144, and 192 h PI (Fig. 3). Livers were swollen and hemorrhagic at 24 h PI (data not shown). Kidneys showed obvious hemorrhage at 192 h PI. Clinical symptoms and gross lesions were not observed in ducklings of the control group.

Fig. 3

Gross lesions of GPV-infected ducks. a, g Thymus and cecum of control groups, respectively. b Bleeding spots on thymus at 1 days PI. c, d and e Petechial hemorrhage of thymus at 2, 4, and 6 days PI, respectively. f No gross lesion in thymus at 8 days PI. h and j Swelling and full of green feces at end of cecum at 1 and 4 days PI. i, k, and l Swelling in entire cecum at 2, 6, and 8 days PI

Histopathological analysis of tissues from ducklings 6 days PI

The histopathological analysis of brain, liver, pancreas, spleen, and kidney sections show no pathological changes at 6 days PI. GPV cannot cause pathological lesions in these tissues, this may be due to the low titer of inoculation compared to the IH and IHC strain inoculation [13, 17]. However, lymphocytes in bursa mildly decreased, small pulmonary veins wall degenerated, and edema and hemorrhage surrounding pulmonary venules. A few lymphocytes in trachea submucosa infiltration, bleeding in myocardiac fibers, epicardium hemorrhage, a large number of lymphocytes infiltrated in small cardiac veins, and thymic hemorrhage. No microscopic lesions were observed in the control group (Fig. 4). The results indicated that GPV could not cause pathological lesions in many tissues; however, serious lesions were found throughout the heart. This phenomenon indicated that GPV targeted the heart, bursa of Fabricius, and thymus.

Fig. 4

Histopathological changes in the GPV-inoculated Shaoxing ducklings at 6 days PI. Con and Ino represent the control and Inoculated groups, respectively. a, b, c, d, i, j, k, l, q, r, and s Brain, bursa of Fabricius, lung, spleen, trachea, kidney, thymus, pancreas, and different parts of the heart, respectively in the control group. e Brain. h Spleen. n Kidney. p Pancreas. f lymphocytes in bursa mildly decreased. g Small pulmonary veins wall degenerated. m A few lymphocytes in trachea submucosa infiltration. o thymic hemorrhage. t a large number of lymphocytes infiltrated in small cardiac veins. u bleeding in myocardiac fibers. v epicardium hemorrhage. Magnification is 200 ×

The distribution of viral DNA in Shaoxing ducklings

GPV distribution in the tissues of ducklings was assessed by amplifying an approximately 400-bp fragment of the GPV genome. This fragment bands could be detected in the liver, lung, jejunum, and ileum during the observation period (Table 3). The 400 bp DNA band could be always detected in cardiac muscle, bursa of Fabricius, duodenum, rectum, brain, kidney and skeletal muscle at 24, 48, 96, 144, 192 h PI (Table 3). The DNA bands in skeletal and cardiac muscle were especially bright. (Additional file 1: Figure S1). Viral DNA band could not be detected in the pancreas and large intestine (Additional file 1: Figure S1). These results suggested that GPVwas primarily remained in muscle.

Table 3 GPV DNA detected by PCR in tissues of Shaoxing ducklings inoculated with GPV


GPV was first isolated in Yangzhou, China in 1961. This disease was continuously reported in several provinces during 1999–2014, including Jiangsu, Shanghai, Guangdong, and Heilongjiang. The genomes of many GPV strains have been sequenced, and phylogenetic analyses have been reported [18, 19]. However, phylogenetic analyses of NS and VP sequences of GPV from the Heilongjiang Province have been rarely reported. We sequenced the NS and VP genes of the H strain isolated from Heilongjiang, China in 2010. The phylogenetic analysis of the H strain VP and NS genes showed that the H strain was closely related to Chinese mainland strains that were isolated over the past 20 years. However, H strain formed a distinct lineage from the HG5/82 strain which was also isolated from Heilongjiang [20]. While GPV strains isolated from ducks were closely related to each other [21]. This suggested that there were more than one GPV genetic groups circulating in Heilongjiang, China. However, the most prevalent GPV strains during the past 10 years arised from a common ancestor.

In our study, we established an experimental model of H strain GPV infection to investigate the pathogenicity of GPV in Shaoxing ducklings. Our results indicated that the incubation period of GPV was less than 24 h, with clinical symptoms becoming obvious at 4 days PI before gradually weakening at 6–8 days PI. The course of the disease took 8 days, and the development of disease had a close relationship with the viral DNA in tissues. GPV was a pantropic virus and replicated rapidly in many tissues of infected Shaoxing ducklings, especially in the heart, and trachea, causing serious pathological lesions. This is the first study of the tissue distribution and histopathological analysis of a Chinese GPV strain in Shaoxing ducklings.

Our results indicated that viral DNA can be detected in 12 types of tissue but not in the pancreas or large intestine. In particular, an especially bright viral DNA band seen in skeletal muscle and cardiac muscle, which was in accordance with the pathological changes in the heart. The goslings inoculated with virulent strain B show myocardial, lympho-histiocytic infiltration and hepatic degeneration at postmortern examination. However, the livers of Shaoxing ducklings infected with strain H were normal and no gross lesions [22, 23]. This result suggested that H GPV mainly targeted muscle. This agrees with previous analyses of the IH and IHC GPV strains. Strain IH was the 20th-Muscovy egg passaged virus and the strain IHC developed from IH strain at the 15th –egg-passage and then passaged in Muscovy duck embryo fibroblasts (MDEFs) 50 times, the two strains DNA were detected constantly in cardiac muscle [13, 17]. The viral DNA was always detected in the trachea except in one duckling at 8 days PI, and these pathological changes in the trachea were obvious. This indicated that Shaoxing ducklings may excrete virus through the trachea for a long period of time, which posed a challenge for preventing and controlling this disease. The histopathological finding that many tissues were normal, including the kidney and brain, did not agree with the viral DNA detection results. The viral DNA detection results suggested that GPV was present and replicate in many tissues. The serious pathological changes in the immune tissues may cause immunosuppression, which would create favorable conditions for coinfections with other pathogens.

This animal experiment demonstrated that the H strain produced moderate pathology in Shaoxing ducklings, and the morbidity, mortality, clinical symptoms and pathological changes caused by the H strain in Shaoxing ducklings was similar to those caused by the N-GPV strain in Cherry Valley ducks, the N-GPV can cause atrophic beak, watery diarrhea, shorter tibia, slightly depauperate liver, and swelling and hemorrhage in the thymus in young mule duckling and Cherry Valley ducklings [9, 10]. This study also provided evidence for N-GPV being closely related to goose-origin parvovirus.


In this study, we replicated the GPV H strain in duck embryos and DEFs and sequenced and analyzed the NS and VP genes. Phylogenetic analyses indicated that there was more than one genetic group circulating in mainland China. The H strain was closely related to strains isolated during 1999–2014. Phylogenetic analysis of the molecular divergence among strains indicated that it is urgently needed to prevent and control goose plague. To our knowledge, this study is the first to explore the distribution and pathological changes of a GPV strain isolated from Northeastern China in infected Shaoxing ducklings. It is also the first study finding that the GPV H strain was moderately pathogenic to Shaoxing ducklings, and the trachea of infected ducklings may secrete the virus. These data provide new insights into outbreaks and the pathogenesis of GPV.



Duck embryo fibroblast


Goose parvovirus




Open reading frame


Polymerase chain reaction


Post inoculation


Capsid protein


  1. 1.

    Glavits R, Zolnai A, Szabo E, Ivanics E, Zarka P, Mato 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:73–89.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Palya V, Zolnai A, Benyeda Z, Kovacs E, Kardi V, Mato T. Short beak and dwarfism syndrome of mule duck is caused by a distinct lineage of goose parvovirus. Avian Pathol. 2009;38:175–80.

    Article  PubMed  Google Scholar 

  3. 3.

    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:477–83.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Le Gall-Recule G, Jestin V. Biochemical and genomic characterization of muscovy duck parvovirus. Arch Virol. 1994;139:121–31.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Zadori Z, Erdei J, Nagy J, Kisary J. Characteristics of the genome of goose parvovirus. Avian Pathol. 1994;23:359–64.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Tatar-Kis T, Mato T, Markos B, Palya V. Phylogenetic analysis of Hungarian goose parvovirus isolates and vaccine strains. Avian Pathol. 2004;33:438–44.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Zhu Y, Zhou Z, Huang Y, Yu R, Dong S, Li Z, Zhang Y. Identification of a recombinant Muscovy duck parvovirus (MDPV) in Shanghai, China. Vet Microbiol. 2014;174:560–4.

    Article  PubMed  Google Scholar 

  8. 8.

    Liu HM, Wang H, Tian XJ, Zhang S, Zhou XH, Qi KZ, Pan L. Complete genome sequence of goose parvovirus Y strain isolated from Muscovy ducks in China. Virus Genes. 2014;48:199–202.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Chen S, Wang S, Cheng X, Xiao S, Zhu X, Lin F, Wu N, Wang J, Huang M, Zheng M, et al. 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:2407–16.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    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:e0140284.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Wang S, Cheng XX, Chen SY, Lin FQ, Chen SL, Zhu XL, Wang JX, Huang MQ, Zheng M. Phylogenetic analysis of VP1 gene sequences of waterfowl parvoviruses from the mainland of China revealed genetic diversity and recombination. Gene. 2016;578:124–31.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Zhang Y, Li Y, Liu M, Zhang D, Guo D, Liu C, Zhi H, Wang X, Li G, Li N, et al. Development and evaluation of a VP3-ELISA for the detection of goose and Muscovy duck parvovirus antibodies. J Virol Methods. 2010;163:405–9.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Limn CK, Yamada T, Takehar K. Detection of goose parvovirus genome by polymerase chain reaction: distribution of goose parvovirus in Muscovy ducklings. Virus Res. 1996;42:167–72.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Bin Wang DS, Chai H, Wang D, Pan X. Analyzing of the VP3 gene sequence of goose parvovirus H strain. Chin J Anim Vet Sci. 2010;37:155–8.

    Google Scholar 

  15. 15.

    Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.

    CAS  PubMed  Google Scholar 

  16. 16.

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Takehara K, Saitoh M, Kiyono M, Nakamura M. Distribution of attenuated goose parvoviruses in Muscovy ducklings. J Vet Med Sci. 1998;60:341–4.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Shien JH, Wang YS, Chen CH, Shieh HK, Hu CC, Chang PC. Identification of sequence changes in live attenuated goose parvovirus vaccine strains developed in Asia and Europe. Avian Pathol. 2008;37:499–505.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Wang S, Cheng XX, Chen SY, Zhu XL, Chen SL, Lin FQ, Li ZL. Genetic characterization of a potentially novel goose parvovirus circulating in Muscovy duck flocks in Fujian Province, China. J Vet Med Sci. 2013;75:1127–30.

    Article  PubMed  Google Scholar 

  20. 20.

    Kong X, Lui G, Liu S, Ran D. Analysis of molecular characteristics of goose parvovirus strain HG5/82. Virol Sin. 2005;20(1):28–32.

    Google Scholar 

  21. 21.

    Li C, Li Q, Chen Z, Liu G. Novel duck parvovirus identified in Cherry Valley ducks (Anas Platyrhynchos Domesticus), China. Infect Genet Evol. 2016;44:278–80.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Kisary J, Derzsy D, Meszaros J. Attenuation of the goose parvovirus strain B. Laboratory and field trials of the attenuated mutant for vaccination against Derzsy's disease. Avian Pathol. 1978;7:397–406.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Derzsy D, Dren C, Szedo M, Surjan J, Toth B, Iro E. Viral disease of goslings. 3. Isolation, properties and antigenic pattern of the virus strains. Acta Vet Acad Sci Hung. 1970;20:419–28.

    CAS  PubMed  Google Scholar 

Download references


This study was supported by Chinese Academy of Agricultural Sciences Fundamental Scientific Research Funds (Y2016PT41) and Key Program of Natural Science Foundation of Heilongjiang Province (ZD2016006).

Availability of data and materials

The datasets analyzed during the current study available from the corresponding author on reasonable request.

Author information




LZ, YN, and HC conceived and designed the study. JL and FY provided GPV H strain. YN and BL performed the experiments and analyzed the data. HY and HH drew the graphs. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Hongyan Chen.

Ethics declarations

Ethics approval and consent to participate

The study was approved by Harbin Veterinary Research Institute and performed in accordance with animal ethics guidelines and approved protocols. The animal Ethics Committee approval number is Heilongjiang-SYXK-2006-032.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Figure S1.

GPV DNA amplification results. M: marker; Li: Liver; Lu: Lung; Ki: Kidney; Bu: Bursa of Fabricius; Pa: Pancreas; Du: Duodenum; Re: Rectum; Je: Jejunum; Il: Ileum; La: Large intestine; Sk: Skeletal muscle; Br: Brain; Tr: Trachea; 24 and 192: Un-inoculated 24 and 192 h; 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10: Duck ID number in Table 3. 1: 4Li; 2: 4Sk; 3: 4Br; 4: 4La; 5: 7 Pa; 6: 2Li; 7: 2 Pa; 8: 3Br; 9:1 Pa;10:3Tr; 11: 4Ca; 12: 3Bu; 13: 4Je; 14: 1Ca; 15: 4Tr; 16:1Li; 17:1Lu; 18: 3Du; 19: 2Ki; 20: 3 Pa; 21:1Tr; 22: 4Lu; 23: 4Ki; 24: 3Ki; 25: 2Ca; 26: 3Re; 27: 7Lu; 28: 2Sk; 29: 2Lu; 30: 7Li; 31: 3La; 32:3Sk; 33: 10Du; 34: 7Tr; 35: 6 Pa; 36: 3Il; 37: 5Je; 38: 7Re; 39: 5 Pa; 40: 3Je; 41: 2Tr; 42: 6Br; 43: 7La; 44: 1Je; 45: 7Du; 46: 7Ca; 47: 1Bu; 48: 7Il; 49: 4Du; 50: 10Je; 51:1La; 52: 9La; 53: 8Bu; 54: 1Br; 55: 3Lu; 56: 2Br; 57:5Br; 58: 7Br; 59: 8Br; 60: 9Bu; 61: 1Il; 62: 10Il; 63: 2La; 64: 4Re; 65:10Br; 66: 10Ca; 67: 9Lu; 68: 9Tr; 69: 7Ki; 70: 4 Pa; 71:9Re; 72:8La; 73: 5Ca; 74: 6Ca; 75: 8Je; 76: 7Bu; 77: 5Li; 78: 6Li; 79: 2Je; 80: 8Ca; 81: 10La; 82: 5Lu; 83: 6Lu; 84: 2Il; 85:8Lu; 86:1Ki; 87:9Br; 88:5Ki; 89:6Ki; 90: 8Ki; 91: 9Ki; 92: 10Ki; 93: 2Bu; 94: 4Bu; 95: 5Bu; 96: 8Li; 97: 6Bu; 98: 3Li; 99:2Du; 100: 5Du; 101: 3Ca; 102: 9Li; 103: 6Du; 104: 4Il; 105:8Du; 106: 9Du; 107: 2Re; 108:10Li; 109:10Lu; 110: 5Re; 111:10Bu; 112: 8 Pa; 113:9 Pa; 114: 6Re; 115: 8Re; 116: 6Je; 117: 9Je; 118: 5Il; 119: 6Il; 120: 1Sk; 121:5Sk; 122 6Sk; 123:10 Pa; 124: 8Sk; 125: 9Sk; 126:10Re; 127:1Du; 128: 10Sk; 129: 5Tr; 130:1Re; 131:7Je; 132:8Tr; 133: 6Tr; 134: 10Tr; 135:9Ca; 136: 5La; 137: 6La; 138: 7Sk; 139: 8Il; 140: 9Il; 141: 24Ca; 142: 24Li; 143: 24Lu; 144: 24Ki; 145: 24Bu; 146: 24 Pa; 147: 24Du; 148: 24Re; 149: 24Il; 150: 24Je; 151: 24La; 152: 24Sk; 153: 24Br; 154: 24Tr; 155:192Ca; 156:192Li; 157:192Lu; 158:192Ki; 159:192Bu; 160: 192 Pa; 161: 192 Du; 162:192Re; 163:192Il; 164:192Je; 165:192La; 166:192Sk; 167:192Br; 168:192Tr; 169:192con. Figure S2. Duck embryo muscle tissues inoculated with serial dilutions. Dentified by PCR. A 1.2 × 102PFU. B 1.2 × 103PFU. C 1.2 × 104PFU. D 1.2 × 105PFU. (PDF 289 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Niu, Y., Zhao, L., Liu, B. et al. Comparative genetic analysis and pathological characteristics of goose parvovirus isolated in Heilongjiang, China. Virol J 15, 27 (2018).

Download citation


  • Goose parvovirus
  • Phylogenetic analysis
  • Pathological changes