Open Access

Identification and genomic characterization of a novel porcine parvovirus (PPV6) in china

Virology Journal201411:203

https://doi.org/10.1186/s12985-014-0203-2

Received: 28 May 2014

Accepted: 14 November 2014

Published: 2 December 2014

Abstract

Background

Parvoviruses are classified into two subfamilies based on their host range: the Parvovirinae, which infect vertebrates, and the Densovirinae, which mainly infect insects and other arthropods. In recent years, a number of novel parvoviruses belonging to the subfamily Parvovirinae have been identified from various animal species and humans, including human parvovirus 4 (PARV4), porcine hokovirus, ovine partetravirus, porcine parvovirus 4 (PPV4), and porcine parvovirus 5 (PPV5).

Methods

Using sequence-independent single primer amplification (SISPA), a novel parvovirus within the subfamily Parvovirinae that was distinct from any known parvoviruses was identified and five full-length genome sequences were determined and analyzed.

Results

A novel porcine parvovirus, provisionally named PPV6, was initially identified from aborted pig fetuses in China. Retrospective studies revealed the prevalence of PPV6 in aborted pig fetuses and piglets(50% and 75%, respectively) was apparently higher than that in finishing pigs and sows (15.6% and 3.8% respectively). Furthermore, the prevalence of PPV6 in finishing pig was similar in affected and unaffected farms (i.e. 16.7% vs. 13.6%-21.7%). This finding indicates that animal age, perhaps due to increased innate immune resistance, strongly influences the level of PPV6 viremia. Complete genome sequencing and multiple alignments have shown that the nearly full-length genome sequences were approximately 6,100 nucleotides in length and shared 20.5%–42.6% DNA sequence identity with other members of the Parvovirinae subfamily. Phylogenetic analysis showed that PPV6 was significantly distinct from other known parvoviruses and was most closely related to PPV4.

Conclusion

Our findings and review of published parvovirus sequences suggested that a novel porcine parvovirus is currently circulating in China and might be classified into the novel genus Copiparvovirus within the subfamily Parvovirinae. However, the clinical manifestations of PPV6 are still unknown in that the prevalence of PPV6 was similar between healthy pigs and sick pigs in a retrospective epidemiological study. The identification of PPV6 within the subfamily Parvovirinae provides further insight into the viral and genetic diversity of parvoviruses.

Keywords

Porcine parvovirusRetrospective studyPhylogenetic analysisPairwise comparison

Background

The family Parvoviridae encompasses small non-enveloped and negative single-stranded DNA viruses, and includes many human and animal pathogens. Porcine parvoviruses (PPV) are important pathogens that cause reproductive failure in swine, resulting in enormous losses in the pig industry worldwide [1]. Epidemiological studies and diagnostic surveys have demonstrated that PPV was the major causative agent responsible for embryonic and fetal death in swine [2]. Parvoviruses are classified into two subfamilies based on their host range: the Parvovirinae, which infect vertebrates, and the Densovirinae, which mainly infect insects and other arthropods. The subfamily Parvovirinae is currently proposed to be divided into eight genera by the International Committee on the Taxonomy of Viruses (ICTV); i.e., Protoparvovirus, Amdoparvovirus, Aveparvovirus, Bocaparvovirus, Dependoparvovirus, Erythroparvovirus, Copiparvovirus, and Tetraparvovirus[3],[4].

A typical parvovirus genome is 4–6.3 kb in length and contains generally two open reading frames (ORFs), encoding for a non-structural protein (NSP) and capsid protein (VP) [3],[5]. An additional ORF; i.e., ORF3, encodes an accessory protein and has been identified in some parvoviruses, such as NP1, an NSP only found in bocavirus and PPV4 [6],[7]. In recent years, with the use of molecular assays and pathogen discovery tools, a number of novel parvoviruses belonging to the subfamily Parvovirinae have been identified in various animal species and humans, including human parvovirus 4 (PARV4), porcine hokovirus (PPV3), ovine partetravirus, PPV4, and PPV5 [6]-[13]. Among them, PPV3 and ovine partetravirus have been shown to have a close relationship to PARV4, and currently have been proposed to be classified into the genus Tetraparvovirus by the ICTV [3],[4]. In addition, PPV4 and PPV5 were both identified in clinical samples from swine herds, forming a distinct branch with bovine parvovirus 2 (BPV2) based on phylogenetic analysis, which led to consideration of these viruses being classified in a novel genus; i.e., Copiparvovirus[11],[12]. Within parvoviruses, six different phylogenetic groups of parvoviruses have been identified from pigs, including classic PPV, PPV2, PPV3, PPV4, PPV5, and porcine bocaviruses (PBoV) [7]-[12]. The reported overall prevalence of parvoviruses in pig herds has varied from 6.4% to 20% for PPV2, 9.7% to 12.4% for PPV3, 1.5% to 39.7% for PBoV, and 2.6% to 6.6% for PPV5 [14]-[17]. The prevalence of PPV and PPV4 among the clinical samples in Chinese swine herds was 17.22% and 2.09%, respectively [11],[18].

Timely identification of novel pathogens is of great significance in the diagnosis, control and prevention of emerging human and animal infectious diseases. The development of the sequence-independent single primer amplification (SISPA) method in recent years has allowed the rapid identification of new viruses [19]. With the advent of this method, some human and animal viruses, including bovine parvovirus and bungowannah virus have been discovered [19],[20]. Using this method, a novel PPV, which is distinct from any known parvovirus, was identified and five almost full-length viral genome sequences were assembled and analyzed in this study.

Results and discussion

Initially, three mixed tissues samples (i.e., spleen, kidney, and tonsils) collected from aborted pig fetuses in Beijing, China were tested for suspicious agents associated with reproductive failure including pseudorabies virus (PRV), porcine reproductive and respiratory syndrome virus (PRRSV), PPV, porcine circovirus-2 (PCV2), classical swine fever virus (CSFV), swine influenzavirus (SIV), Japanese encephalitis virus (JEV), and Brucella suis using RT-PCR or PCR. All samples were found to be negative for the pathogens mentioned above and a nonspecific etiologic agent associated with sow abortion was detected. To explore whether there were any unknown viruses present in the aborted pig fetuses, a SISPA experiment was conducted [21],[22]. The sequencing reads were obtained and subjected to BLAST homology searching with sequences deposited in GenBank (BLASTx and tBLASTx programs; http://blast.ncbi.nlm.nih.gov.blast.cgi). The obtained nucleotide sequences included gene fragments of swine, virus, bacteria, and unknown sequences. Of the unknown nucleotide acid sequences, six sequences ranging in size from 298 bp to 650 bp in length exhibited no significant similarity to database sequences. Furthermore, the deduced amino acid sequences of these sequences exhibited either no putative conserved domains (four of six) or contained putative conserved domains of parvoviruses (two of six; BLASTp, E scores < =10−9). These results indicated the presence of a possible novel parvovirus. On the basis of its homology to parvoviruses and its swine host, this DNA virus was tentatively named porcine parvovirus type 6 (PPV6) after the two recently identified PPVs; i.e., PPV4 and PPV5 (data not shown) [12],[21]. Attempts to stably passage the PPV in PK15 (swine kidney), Vero (African green monkey kidney), and Marc 145 (fetal rhesus monkey kidney) cells were unsuccessful (data not shown).

To investigate the prevalence of PPV6 in clinical samples and the clinical picture associated with PPV6 infection. A retrospective epidemiological study of PPV6 infection was performed by PCR. A total of 171 field samples (160 of sera and 11 of tissues) were collected from apparently healthy pigs and sick pigs with similar reproductive system symptoms from four provinces in China (Beijing, Jiangsu, Tianjin, and Sichuan) during the period of 2010–2013. Among these samples, 48 specimens collected from farms experiencing reproductive disease included 6 aborted fetuses, 4 piglets, 26 sows, and 12 finishing pigs. Taking into consideration that the infection of porcine parvovirus is asymptomatic in growing pigs and multiparous sows, to investigate the distribution of PPV6 in different age group, one half of sow samples and all samples of finishing pigs were collected from healthy pigs. The other 123 specimens collected from healthy farms were all of finishing pigs and were previously identified as negative for PRV, PRRSV, PCV2, PPV, and CSFV. The results showed the prevalence of PPV6 in aborted pig fetuses and piglets(50% and 75%, respectively) was apparently higher than that in finishing pigs and sows (15.6% and 3.8% respectively) no matter what their health status. Furthermore, the prevalence of PPV6 in finishing pig was similar in affected and unaffected farms (i.e. 16.7% vs. 13.6%-21.7%) (Table 1) [23]. This finding indicates that animal age, perhaps due to increased innate immune resistance, strongly influences the level of PPV6 viremia, which is similar to the study result of PRRSV [24].
Table 1

Detection of PPV6 from field samples using PCR

Age

Health status of sampled farms

Location

Type of sample

No. positive/no. tested samples

Prevalence (%)

Overall prevalence (%)

Aborted fetuses

Reproductive failure

Beijing

MTS*

3/6

50

50

piglet

Reproductive failure

Beijing

Serum or MTS

3/4

75

75

Finishing pigs

Reproductive failure

Beijing

Serum

2/12

16.7

15.6

Healthy

Beijing

Serum

6/42

14.2

Jiangsu

Serum

5/36

13.9

Tianjin

Serum

3/22

13.6

Sichuan

Serum

5/23

21.7

Sow

Reproductive failure

Beijing

Serum

1/26

3.8

3.8

MTS*, Mixed Tissues Sample.

It is interesting to note that, in samples from aborted pig fetuses collected in Beijing, only PPV6 was detected. Targeted studies are needed to investigate the role of this virus in the sows with reproductive failure [2]. However, the number of clinical samples detected in this study was limited and more extensive epidemiologic studies are needed to elucidate the precise role of PPV6 as a causal agent for reproductive failure. In recent years, a number of novel parvoviruses belonging to the subfamily Parvovirinae have been identified, including PARV, PPV4, and PPV5. Even though these new parvoviruses have been extensively studied, the clinical manifestations of these viruses are still unknown. Confirming a causal relationship between a virus and the observed symptoms is an important but difficult issue that normally requires multiple separate studies to be ultimately resolved [19].

Furthermore, PPV6 was detected in samples from all four provinces suggesting that PPV6 has been circulating in a wide area of China. In light of the importance of swine as a potential source of genetic diversity for parvoviruses, identifying possible counterparts of PPV6 in humans and other animals is important in understanding its epidemiology, evolution, and pathogenesis [1].

Based on SISPA-generated fragment sequences, diverging primers were designed to obtain the intervening portions of the genome. The terminal sequences were then acquired using a 5΄ and 3΄ RACE method [8]. Near full-length genome data were generated from 5 positive samples (i.e., PPV6-BJ, BJ2, SC, JS, and TJ) and deposited in GenBank under accession Nos. KF999681-KF999685. Of these sequenced PPV6 isolates, strain BJ was from aborted pig fetuses derived from herds affected with an epizootic reproductive failure in Beijing, whereas the other four strains (i.e., BJ2, SC, JS, and TJ) were from finishing pigs without any clinical signs, independently identified in Beijing, Sichuan, Jiangsu, and Tianjin.

The near full-length genome of PPV6 contained 6,136 bases (BJ, BJ2, and SC strains) or 6,148 bases (JS and TJ strains) with a G + C content of 46.7–47.1%. Their genome sizes were expected to be larger, and sequencing of the ends was hampered by hair structures. Putative ORFs were obtained using the ORF Finder tool at NCBI (http://www.ncbi.nlm.nih.gov/projects/gorf/) and then were identified by protein blast analysis in the NCBI RefSeq database. The genome organization of PPV6 was similar to that of other parvoviruses, with the characteristic gene order 5′UTR-ORF1-ORF2-3′UTR (Figure 1A) [5]. The ORF1 encodes a putative NSP of 662 amino acids and ORF2 encodes a putative VP of 1,189 amino acids. The predicted sizes of VP (132 kDa) were bigger than most other parvoviruses (most parvoviruses are <90 kDa, including PPV4 and except for BPV2, which was 105 kDa, and PPV5, which was 112 kDa). The sequence similarity between the tested PPV6 strains ranged from 97.1%–99.6% for nucleotides (whole genome) and 97.3%–99.9% for amino acids (NSP and VP) respectively (data not shown). The polymorphic diversity of the VP protein sequences was greater than for NSP protein sequences. The overall amino acid diversity of NSP and VP of five PPV6 strains was 0.0%–0.4% and 0.1%–2.7%, respectively (data not shown). These results suggest that all isolates were closely related to each other genetically.
Figure 1

Genome structure and phylogenetic analysis: bootstrapped neighbor-joining tree based on full-length genomes and NSP or VP amino acid sequences of PPV6 and the Parvoviridae . A. The solid line represents the DNA genome and the boxes represent the ORFs. The position and size of putative ORFs are indicated. B. A phylogenetic tree was constructed based on 5 near full-length genomes (6.1 kb) of PPV6 and 31 genomes from GenBank using amaximum-likelihood method (GTR + G + I) with bootstrap analysis of 1,000 replicates. C and D. Phylogenetic trees were constructed based on amino acid sequences of NSP and VP using a Poisson correction model with bootstrap analysis of 1,000 replicates. Phylogenetic analyses were conducted using the MEGA6.0 (http://www.megasoftware.net) program. Branches corresponding to partitions reproduced in <70% bootstrap replicates are collapsed. The evolutionary distance scale is shown at the bottom. The sequence of GM-densovirus (NC004286) is used as outgroup to root the tree. The asterisk indicates the new genus proposed by ICTV.

Five PPV6 sequences were aligned with 31 reference genomic sequences of viruses in the subfamilies of Parvovirinae and Densovirinae from GenBank using the ClustalX (Ver.1.81) program [25]. Phylogenetic analyses based on the full-length genomes with a maximum-likelihood method (GTR + G + I) showed that PPV6 was distinct from any known parvoviruses and represented a deeply rooted lineage between BPV2 and PPV4 and PPV5 (Figure 1B). The basal phylogenetic position of PPV6 suggested early divergence from other mammalian parvovirus species, and PPV6 may share a common ancestor with PPV4 and PPV5.

Evolutionary trees were also constructed separately for the putative protein sequence of NSP and VP with a Poisson correction model using 500 bootstrap replicates [25]. The topologies of these trees were similar to that of the full-length genome tree (Figures 1C and 1D), and PPV6 formed a distinct cluster within parvoviruses. They also differed from other parvoviruses by their relatively large predicted VP protein. These analyses indicated that PPV6 viruses are likely to belong to the novel genus Copiparvovirus within the subfamily Parvovirinae[12].

Pairwise comparisons were performed for the nucleotide sequences and predicted amino acids sequences of PPV6 with other parvoviruses. The results showed the genomes of the five strains of PPV6 shared 20.5%–42.6% DNA sequence identity with other members of Parvovirinae and are most closely related to PPV4. At the amino acid level, PPV6 exhibited the largest amino acid similarity in NSP with PPV4 (49.8%) and in VP with PPV5 (29.8%; Table 2). PPV6 was found to possess <30% amino acid similarity in the NSP to those of other genera, whereas it exhibited 28.1%–49.8% (BPV2 and PPV4, respectively) identity to that of the genus Copiparvovirus. Based on the new criteria for a genus within the subfamily Parvovirinae by ICTV (generally >30% amino acid identical to NS1 proteins within a genus but <30% identical to other genera), the PPV6 would be classified as a novel species of the genus Copiparvovirus[4].
Table 2

Nucleotide and amino acid identities of PPV6 with those of other parvoviruses

Virus

PPV6-BJ

PPV4-JS0918a(HM031134)

PPV5-ND564(JX896322)

BPV2(NC006259)

Genome(nt)

NSP(aa)

VP(aa)

Genome(nt)

NSP(aa)

VP(aa)

Genome(nt)

NSP(aa)

VP(aa)

Genome(nt)

NSP(aa)

VP(aa)

Tetraparvovirus

PARV4 AY622943

27.5

22.4

17.1

30.5

19.8

17.2

29.5

21.0

19.9

29.3

19.5

18.0

PPV3F4-7SM JQ868704

22.4

14.2

17.3

24.4

11.3

17.5

23.0

12.0

19.8

23.7

12.8

17.3

PPV2US-135 JX101461

27.3

24.1

17.9

30.1

19.3

15.4

27.9

19.7

18.3

27.5

19.0

18.0

Erythroparvovirus

B19-C39 FN598217

24.8

18.7

14.2

27.7

17.1

16.3

28.5

17.9

17.2

25.6

19.5

14.5

Dependoparvovirus

AAV1 AF063497

22.5

25.2

11.6

25.7

22.9

15.8

24.6

23.7

14.5

24.5

22.8

13.8

Bocaparvovirus

PBoV3 NC016031

23.2

15.9

8.9

27.1

15.9

11.7

28.0

16.0

10.8

26.5

16.3

11.3

BPV NC001540

24.1

15.0

9.9

29.5

14.3

11.8

30.5

13.6

11.6

26.7

13.2

11.0

Amdoparvovirus

Amdovirus X97629

20.5

10.7

6.6

22.8

11.9

7.6

22.0

11.9

9.7

22.5

12.1

8.3

Protoparvovirus

PPV China AY583318

21.4

15.4

7.9

23.2

15.1

107

22.5

15.0

9.6

23.8

12.8

10.4

Copiparvovirus

PPV4JS0918aHM031134

42.6

49.8

23.9

100

100

100

62.9

85

38.2

34.6

32.5

19.8

PPV5ND564 JX896322

41.3

49.1

29.8

62.9

85

38.2

100

100

100

33.1

31.9

22.8

BPV2 NC006259

31.7

28.1

18.0

34.6

32.5

19.8

33.1

31.9

22.8

100

100

100

PPV6-BJ

100

100

100

42.6

49.8

23.9

41.3

49.1

29.8

31.7

28.1

18.0

Densovirinae

GM-densovirusNC004286

24.9

8.5

4.4

24.9

8.4

5.0

24.5

8.8

7.1

25.0

10.2

7.8

All sequence identities were calculated using clustal X. Nt, nucleotide; aa, amino acid.

Even though the overall amino acid homology of PPV6 with other parvoviruses is low, the conserved sequence motifs important for the function of parvoviruses were observed. In the alignment, the conserved replication initiator motif (I and II), NTP-binding and helicase domain (A, B, and C) in the NSP (data not shown) identified within Parvovirinae was also found in PPV6 [12],[26]. Within the VP unique region, detailed characterization at the amino acid level revealed PPV6 possesses the conserved motifs of the Ca2+ binding loop (YXGXR) and the catalytic center (HDXXY) of the putative secretory phospholipase A2 (PLA2) motif (Figure 2), which are present in the capsid protein of PPV5 but are lacking in PPV4. Furthermore, the conserved motifs of the Ca2+ binding loop of PLA2 is the “YXGXR” motif in PPV6, rather than the “YXGXG” or “YXGXF” motif found in most parvoviruses [6]. The PLA2 activity of the alternate motif found in PPV6 needs further study. However, alignment analysis suggested PPV6 is not closely related to any other known human or animal parvovirus, and represents a potential novel species of the genus Copiparvovirus, corresponding to the tBLASTx results. Although recombination phenomena have been demonstrated within parvoviral species [27],[28], in this study, no recombination signal was found in PPV6 with other parvoviruses by SimPlot analysis (data not shown), which indicated PPV6 is not a recombination of other parvoviruses. Since parvoviruses utilize host DNA polymerase, they are generally considered to be relatively stable, and the increasing identification of these novel parvoviruses with extensive genetic diversity suggests that the evolution of parvoviruses is far more complicated [11],[12].
Figure 2

Sequence alignment of the putative phospholipase A2 motif of PPV6 with those of other parvoviruses. The conserved amino acids of the Ca2+ binding loop (YXGXR) and the catalytic residues (HD--Y) are indicated at the bottom of the alignment. The GenBank numbers of the sequences are indicated.

Conclusion

In summary, we described the identification and genome characterization of a novel parvovirus (PPV6) from pigs, the closest neighbors of which were PPV4 and PPV5. The full genome of PPV6 is approximately 6.1 kb in length, and the genomic organization of PPV6 is similar to PPV5 but not to PPV4, which contains an ORF3 in the middle of the viral genome. However, the genome of PPV6 was slightly distinguished from PPV5 by the larger VP gene and the larger genome size. Phylogenetic analysis demonstrated that PPV6, together with PPV4 and PPV5, form a distinct branch that is genetically different from viruses of previously defined genera in the subfamily Parvovirinae, and might be classified into the novel genus Copiparvovirus[4],[12]. The identification of a novel PPV; i.e., PPV6, within the subfamily Parvovirinae provides, further insight into the viral and genetic diversity of parvoviruses. Further study is needed to explore the exact roles of PPV6, especially regarding its host range, geographical distribution, and relatedness to disease [7],[9],[10],[12]. Although PPV4, PPV5, and PPV6 were discovered in recent years, the biological characteristics of these viruses and relatedness to diseases are still not fully understood. Future endeavors to culture PPV6 will help address these questions.

Materials and methods

Ethics statement

All samples utilized were originated from sick or healthy pig case submissions to the Veterinary Diagnostic Laboratory of the China Animal Disease Control Center (CADC-VDL) for diagnostic workups. The protocol for this study was approved by the Biosafety Committee of the CADC, Beijing, China.

Sample collection

All samples were obtained from Beijing, Jiangsu, Tianjin, and Sichuan in China with the assistance of local veterinary practitioners. The samples were transported to the CADC-VDL at low temperatures. The 123 samples from apparently healthy pigs were collected during the epidemiological investigation of PRV, PRRSV, PPV, PCV2, and CSFV in China from 2010 to 2013. All 123 samples used in this study were identified as negative for the five pig pathogens listed above. A total of 48 samples from sick pigs with reproductive disease were collected in Beijing, and these included 6 mixed tissue samples from aborted fetuses, 4 serum or mixed tissues samples from piglets, 12 serum or mixed tissues samples from finishing pigs, and 26 serum samples from sows, which were collected in 2013 using procedures described previously [29].

Nucleic acid extraction and routine detection

Total viral RNA and DNA were extracted directly from sera and tissue samples separately using a RNeasy Mini kit (Qiagen) and QIAamp viral DNA mini kit (Qiagen) following the manufacturer’s instructions. The RNA and DNA was re-dissolved in RNase- and DNase-free water and stored at −80°C until further processing. Conventional PCR or RT-PCR for PRV, PRRSV, PPV, PCV2, CSFV, swine influenza virus (SIV), Japanese encephalitis virus (JEV), and Brucella suis were carried out at the CADC.

SISPA

Random detection of unknown virus genomic DNA or RNA from serum samples was performed using SISPA as previously described with some modifications [6]. To detect viral RNA, first-strand cDNA synthesis was performed using SuperScript™ III reverse transcriptase (Invitrogen, Carlsbad, CA, USA) primed with 5 pmol random hexamer primer 42°C for 1 h. The second-strand cDNA was then synthesized by adding 1 U of RNase H and 30 U of DNA polymerase I (Invitrogen) at 12°C for 1 h followed by 1 h at 22°C. To detect viral DNA, a complementary strand was generated by the incubation of extracted DNA with Klenow fragment (exo-; New England Biolabs, NEB, USA) and 5 pmol of random hexamers at 37°C for 1 h. Double-stranded DNA samples were then digested with the restriction enzyme TaqI (Takara, Dalian, China). Adaptors composed of hybridized oligonucleotides Taq24 (5′ ATCCCTCGGATAGCACAGTTGCCT3′) and NTaq11 (5′CGAGGCAACTG3′; 80 pmol) were then ligated to the restricted DNA using T4 DNA ligase (New England Biolabs, NEB) (1). Two microliters of the ligation reaction were used to prime PCRs containing 50 pmol of Taq24 and 2.5 units of Taq polymerase (Promega Madison, Wisconsin, USA). Cycling conditions were as follows: 40 cycles, which consisting of 94°C for 1 min and 72°C for 3 min [21],[22].

The Amplified PCR products were then analyzed by agarose gel electrophoresis. The distinct DNA bands were gel-purified with an E.Z.N.A.TM Gel Extraction Kit (OMEGA, USA) and ligated with pGEM®-T Easy Vector System (Promega). The recombinant plasmids were transformed into Escherichia coli Trans1 T1 competent cells (Transgen, China). A total of 45 positive clones were sequenced in both directions using universal primers T7 and SP6 promoter-specific primers with an ABI Prism 3730 DNA sequencer (Applied Biosystems, USA) at Invitrogen (Beijing, China). The sequencing reads were trimmed, clustered, and subjected to BLAST homology searching with sequences deposited in GenBank (BLASTn and BLASTx; http://blast.ncbi.nlm.nih.gov.blast.cgi).

PCR detection of PPV6

According to the sequence obtained from SISPA, two primers (PPV6F: 5′-GTCAAAGTGGGAACCCAATTG-3′ and PPV6R: 5′-CCTGGACAGCAAGAAGAAATG-3′) were designed to amplify a 371 nucleotides region within the ORF2. All samples were screened for PPV6 by PCR. The thermal cycling conditions were 94°C for 3 min, followed by 35 cycles of 94°C for 30 s, 52°C for 45 s, 72°C for 45 s, and a final elongation step at 72°C for 10 min. Finally, the PCR products were analyzed on 1.5% agarose gel electrophoresis ultraviolet imaging. Positive samples were determined with 371 bp amplified products.

Whole genome sequencing

To dissect the phylogenetic position of PPV6, five representative viruses were selected for whole-genome sequencing. Four isolates (i.e., BJ2, SC, JS, and TJ) from healthy pig herds were selected to represent four geographical regions in China (i.e., Beijing, Sichuan, Tianjin, and Jiangsu). Also, the PPV6-BJ strain was selected to represent the sick pig herd with reproductive failure in Beijing. First, based on SISPA-generated fragment sequences, diverging primers were designed to obtain the intervening portions of the genome. Second, the terminal sequences were then acquired using the 5΄ and 3΄ RACE method (Invitrogen, Cat. No. 18373 and 18374) [8]. For 5′RACE, the PPV6 genome was first linearly amplified (60 cycles, with 1 cycle consisting of 94°C for 30 s, 50°C for 30 s, and 72°C for 2 min) using Taq polymerase and PPV6-specific primer p5RACE (5′-TGCGCTTATCTTCATTCAGAC-3′). Amplification products were then purified and a poly(C) tail was added to the 3′ end using deoxycytidine and terminal deoxynucleotidyl transferase. The 5′ region was then amplified with 5 units of Taq polymerase using an abridged anchor primer and Taq polymerase (5 U). For 3′RACE, gene specific primers p3RACE (5′-TATGGCCCATGTAAACGCATC-3′) were designed based on the available partial sequence and was used in combination with the Oligo-dT anchor primer. RACE-PCR reactions were performed using Taq DNA polymerase for 60 cycles (1 cycle consisting of 94°C for 45 s, 55°C for 45 s, and 72°C for 2 min). The products were analyzed by agarose gel electrophoresis and then sequenced according to standard protocols. Third, the acquired sequences were assembled into a complete genome with the aid of Vector NTI Site 9.0 software (Invitrogen). Finally, the five near full-length genome data from 5 positive samples (i.e., PPV6-BJ, BJ2, SC, JS, and TJ) were deposited in GenBank under accession Nos. KF999681-KF999685.

Phylogenetic analysis

Sequences alignments were performed using the ClustalX (Ver.1.81) program. A neighbor-joining (NJ) tree was constructed using MEGA version 6 software (http://www.megasoftware.net). Reliability of the NJ tree was calculated using 1,000 bootstrap replicates. In addition to the PPV6 viruses, the complete sequences of various other parvoviruses were obtained from GenBank.

Virus isolation

Virus isolation of PPV6 was attempted in PK15 (swine kidney), Vero (African green monkey kidney), and Marc 145 cells, as previously described with modifications [30]. Cells were seeded at 1.0 × 105 cells ml−1 and, after 1 h incubation, they were inoculated with sample. Cells were grown in DMEM supplemented with penicillin (100 U/ml) and streptomycin (100 μg/ml), and were incubated at 37 °C in a humid environment containing 5 % CO2. If no CPE was observed at 7 days post-inoculation, the plates were frozen and thawed once and the supernatants were inoculated on new cells for a second passage. Inoculated cells at each passage were also tested by PCR, as described above. If CPE and PCR were negative after four passages, the virus isolation result was considered negative.

Declarations

Acknowledgments

We thank the laboratory scientists at the porcine disease laboratory of the China Animal Disease Control Center for sample collection and result discussion. The authors also thank Dr. Yuliang Liu and Yingju Xia for kind help in editing the manuscript. This study was supported by a grant from Modern Agricultural Industry Technology System Programs (GWZJ-2009-05).

The GenBank accession numbers for the nucleotide sequences of the genomes of PPV6 determined in this study are KF999681-KF999685.

Authors’ Affiliations

(1)
China Animal Disease Control Center Veterinary Diagnostic Center
(2)
Beijing Entry-Exit Inspection and Quarantine Bureau

References

  1. Hueffer K, Parrish CR: Parvovirus host rang, cell tropism and evolution. Curr Top Microbiol Immunol 2003, 6: 392-398.Google Scholar
  2. Van-Leengoed LA, Vos J, Gruys E, Rondhuis P, Brand A: Porcine Parvovirus infection: review and diagnosis in a sow herd with reproductive failure. Vet Q 1983, 5: 131-141. 10.1080/01652176.1983.9693887PubMedView ArticleGoogle Scholar
  3. Knowles NJ, Hovi T, Hyypiä T, King AMQ, Lindberg AM, Pallansch MA, Palmenberg AC, Simmonds P, Skern T, Stanway G, Yamashita T, Zell R: Family Parvoviridae. In Virus Taxonomy 9th report of the International Committee on Taxonomy of Viruses. Academic Press, London; 2011:405-425.Google Scholar
  4. Cotmore SF, Agbandje-McKenna M, Chiorini JA, Mukha DV, Pintel DJ, Qiu J, Soderlund-Venermo M, Tattersall P, Tijssen P, Gatherer D, Davison AJ: The family Parvoviridae . Arch Virol 2014, 159: 1239-1247. 10.1007/s00705-013-1914-1PubMedPubMed CentralView ArticleGoogle Scholar
  5. López-Bueno A, Villarreal LP, Almendral JM: Parvovirus variation for disease: a difference with RNA viruses? Curr Top Microbiol Immunol 2006, 299: 349-370.PubMedGoogle Scholar
  6. Cheng WX, Li JS, Huang CP, Yao DP, Liu N, Cui SX, Jin Y, Duan ZJ: Identification and nearly full-length genome characterization of novel porcine bocaviruses. PLoS One 2010, 5: e13583.Google Scholar
  7. Cheung AK, Wu G, Wang D, Bayles DO, Lager KM, Vincent AL: Identificationand molecular cloning of a novel porcine parvovirus. Arch Virol 2010, 155: 801-806. 10.1007/s00705-010-0646-8PubMedView ArticleGoogle Scholar
  8. Jones MS, Kapoor A, Lukashov VV, Simmonds P, Hecht F, Delwart E: New DNA viruses identified in patients with acute viral infection syndrome. J Virol 2005, 79: 8230-8236. 10.1128/JVI.79.13.8230-8236.2005PubMedPubMed CentralView ArticleGoogle Scholar
  9. Lau SK, Woo PC, Tse H, Fu CT, Au WK, Chen XC, Tsoi HW, Tsang TH, Chan JS, Tsang DN, Li KS, Tse CW, Ng TK, Tsang OT, Zheng BJ, Tam S, Chan KH, Zhou B, Yuen KY: Identification of novel porcine and bovine parvoviruses closely related to human parvovirus 4. J Gen Virol 2008, 89: 1840-1848. 10.1099/vir.0.2008/000380-0PubMedView ArticleGoogle Scholar
  10. Tse H, Tsoi HW, Teng JL, Chen XC, Liu H, Zhou B, Zheng BJ, Woo PC, Lau SK, Yuen KY: Discovery and genomic characterization of a novel ovine partetravirus and a new genotype of bovine partetravirus. PLoS One 2011, 6: e25619.Google Scholar
  11. Huang L, Zhai SL, Cheung AK, Zhang HB, Long JX, Yuan SS: Detection of a novel porcine parvovirus, PPV4, in Chinese swine herds. Virol J 2010, 7: 333.Google Scholar
  12. Xiao CT, Gimenez-Lirol LG, Jiang YH, Halbur PG, Opriessnig T: Characterization of a Novel Porcine ParvovirusTentatively Designated PPV5. PLoS One 2013, 8: e65312.Google Scholar
  13. Blomstro¨m AL, Belak S, Fossum C, McKillen J, Allan G, Wallgren P, Berg M: Detection of a novel porcine boca-like virus in the background of porcine circovirus type 2 induced postweaning multisystemic wasting syndrome. Virus Res 2009, 146: 125-129. 10.1016/j.virusres.2009.09.006View ArticleGoogle Scholar
  14. Zhang HB, Huang L, Liu YJ, Lin T, Sun CQ, Deng Y, Wei ZZ, Cheung AK, Long JX, Yuan SS: Porcine bocaviruses: genetic analysis and prevalence in Chinese swine population. Epidemiol Infect 2011, 139: 1581-1586. 10.1017/S0950268811000847PubMedView ArticleGoogle Scholar
  15. Csa’gola A, Lorincz M, Cadar D, Tombacz K, Biksi I, Tuboly T: Detection, prevalence and analysis of emerging porcine parvovirus infections. Arch Virol 2012, 157: 1003-1010. 10.1007/s00705-012-1257-3View ArticleGoogle Scholar
  16. Xiao CT, Gimenez-Lirola LG, Halbur PG, Opriessnig T: Increasing porcine PARV4 prevalence with pig age in the U.S. pig population. Vet Microbiol 2012, 160: 290-296. 10.1016/j.vetmic.2012.05.038PubMedView ArticleGoogle Scholar
  17. Xiao CT, Gerber PF, Gimenez-Lirola LG, Halbur PG, Opriessnig T: Characterization of porcine parvovirus type 2 (PPV2) which is highly prevalent in the USA. Vet Microbiol 2013, 161: 325-330. 10.1016/j.vetmic.2012.07.038PubMedView ArticleGoogle Scholar
  18. Xu J, Guo WZ, Chen Y, Xu ZW, Wang XY: Epidemiology Survey of Porcine Parvovirus and Porcine Circovirus Type 2 in Sichuan Province. Anim Husbandry Feed Sci 2012, 4: 130-132.Google Scholar
  19. Allander T, Emerson SU, Engle RE, Purcell RH, Bukh J: A virus discovery method incorporating DNase treatment and its application to the identification of two bovine parvovirus species. Proc Natl Acad Sci U S A 2001, 98: 11609-11614. 10.1073/pnas.211424698PubMedPubMed CentralView ArticleGoogle Scholar
  20. Kirkland PD, Frost MJ, Finlaison DS, King KR, Ridpath JF, Gu X: Identification of a novel virus in pigs-Bungowannah virus: a possible new species of pestivirus. Virus Res 2007, 129: 26-34. 10.1016/j.virusres.2007.05.002PubMedView ArticleGoogle Scholar
  21. Sachsenröder J, Twardziok S, Hammerl JA, Janczyk P, Wrede P, Hertwig S, Johne R: Simultaneous identification of DNA and RNA viruses present in pig faeces using process-controlled deep sequencing. PLoS One 2012, 7: e34631.Google Scholar
  22. Wang F, Wei Y, Zhu C, Huang X, Xu Y, Yu L, Yu XP: Novel parvovirus sublineage in the family of Parvoviridae. Virus Genes 2010, 41: 305-308. 10.1007/s11262-010-0506-3PubMedView ArticleGoogle Scholar
  23. Simmonds P, Manning A, Kenneil R, Carnie FW, Bell JE: Parenteral transmission of the novel human parvovirus PARV4. Emerg Infect Dis 2007, 13: 1386-1388. 10.3201/eid1309.070428PubMedPubMed CentralView ArticleGoogle Scholar
  24. Klinge KL, Vaughn EM, Roof MB, Bautista EM, Murtaugh MP: Age-dependent resistance to Porcine reproductive and respiratory syndrome virus replication in swine. Virol J 2009, 6: 177.Google Scholar
  25. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S: MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol Biol Evol 2013, 30: 2725-2729. 10.1093/molbev/mst197PubMedPubMed CentralView ArticleGoogle Scholar
  26. Sukhumsirichart W, Attasart P, Boonsaeng V, Panyim S: Complete nucleotide sequence and genomic organization of hepatopancreatic parvovirus (HPV) of Penaeus monodon. Virology 2006, 346: 266-277. 10.1016/j.virol.2005.06.052PubMedView ArticleGoogle Scholar
  27. Kapoor A, Slikas E, Simmonds P, Chieochansin T, Naeem A, Shaukat S, Alam MM, Sharif S, Angez M, Zaidi S, Delwart E: A newly identified bocavirus species in human stool. J Infect Dis 2009, 199: 196-200. 10.1086/595831PubMedPubMed CentralView ArticleGoogle Scholar
  28. Lukashov VV, Goudsmit J: Evolutionary relationships among parvoviruses: virus-host coevolution among autonomous primate parvoviruses and links between adeno-associated and avian parvoviruses. J Virol 2001, 75: 2729-2740. 10.1128/JVI.75.6.2729-2740.2001PubMedPubMed CentralView ArticleGoogle Scholar
  29. Ni J, Yang S, Bounlom D, Yu X, Zhou Z, Song J, Khamphouth V, Vatthana T, Tian K: Emergence and pathogenicity of highly pathogenic Porcine reproductive and respiratory syndrome virus in Vientiane, Lao People’s Democratic Republic. J Vet Diagn Invest 2012, 24: 349-354. 10.1177/1040638711434111PubMedView ArticleGoogle Scholar
  30. Zeeuw EJ, Leinecker N, Herwig V, Selbitz HJ, Truyen U: Study of the virulence and cross-neutralization capability of recent porcine parvovirus field isolates and vaccine viruses in experimentally infected pregnant gilts. J Gen Virol 2007, 88: 420-427. 10.1099/vir.0.82302-0PubMedView ArticleGoogle Scholar

Copyright

© Ni et al.; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Advertisement