key: cord-0000545-u4649rtx authors: Shan, Tongling; Lan, Daoliang; Li, Linlin; Wang, Chunmei; Cui, Li; Zhang, Wen; Hua, Xiuguo; Zhu, Caixia; Zhao, Wei; Delwart, Eric title: Genomic Characterization and High Prevalence of Bocaviruses in Swine date: 2011-04-15 journal: PLoS One DOI: 10.1371/journal.pone.0017292 sha: 21caf0c49864b6af4c62e68e3ff7374e14c16537 doc_id: 545 cord_uid: u4649rtx Using random PCR amplification followed by plasmid subcloning and DNA sequencing, we detected bocavirus related sequences in 9 out of 17 porcine stool samples. Using primer walking, we sequenced the nearly complete genomes of two highly divergent bocaviruses we provisionally named porcine bocavirus 1 isolate H18 (PBoV1-H18) and porcine bocavirus 2 isolate A6 (PBoV2-A6) which differed by 51.8% in their NS1 protein. Phylogenetic analysis indicated that PBoV1-H18 was very closely related to a ∼2 Kb central region of a porcine bocavirus-like virus (PBo-LikeV) from Sweden described in 2009. PBoV2-A6 was very closely related to the porcine bocavirus genomes PBoV-1 and PBoV2 from China described in 2010. Among 340 fecal samples collected from different age, asymptomatic swine in five Chinese provinces, the prevalence of PBoV1-H18 and PBoV2-A6 related viruses were 45–75% and 55–70% respectively, with 30–47% of pigs co-infected. PBoV1-A6 related strains were highly conserved, while PBoV2-H18 related strains were more diverse, grouping into two genotypes corresponding to the previously described PBoV1 and PBoV2. Together with the recently described partial bocavirus genomes labeled V6 and V7, a total of three major porcine bocavirus clades have therefore been described to date. Further studies will be required to elucidate the possible pathogenic impact of these diverse bocaviruses either alone or in combination with other porcine viruses. Members of the Parvoviridae family are small non-lipid enveloped viruses with a diameter of 18-26 nm, icosahedral symmetry (T = 1), encoded by a single-stranded linear DNA genome of approximately 4,000 to 6,000 nucleotides (nt) [1] . The family includes two subfamilies: Densovirinae, Parvovirinae. The subfamily of Densovirinae contains four genera: Densovirus, Iteravirus, Brevidensovirus and Pefudensovirus, which infect only invertebrates [1] . The Parvovirinae subfamily is currently subdivided into five genera: Parvovirus, Erythrovirus, Dependovirus (adeno-associated virus), Amdovirus, and Bocavirus, infecting vertebrates [1] . The bocavirus genus was recently assigned by the International Committee on Taxonomy of Viruses (ICTV) [1] to parvovirus genomes containing a third ORF (labeled NP1) between the NS1 and VP1/VP2 genes [2] . Bocaviruses were first identified in bovine and canine [3, 4] , samples from which it derives its genus name [1, 5] . Presently, the bocavirus genus contains eight members: bovine parvovirus, canine minute virus (CnMV), human bocavirus 1-4 (HBoV1-4), a gorilla bocavirus and a partially sequenced chimpanzee bocavirus [1, 6, 7] . The first human bocavirus (HBoV) was found in the nasopharyngeal secretion of a child with respiratory problems using a methodology closely related to that used here [8] . HBoV has been associated with lower respiratory tract symptoms and possibly diarrhea [5, [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] , and shows a very low degree of genetic variability worldwide [6, 23, 24] . HBoV2 was first reported in the stool of Pakistani children with non-polio acute flaccid paralysis (AFP) [10] , and then in Australian children and Chinese children with diarrhea [9, 25] . HBoV3 was first reported in the stool of Australian children with diarrhea [9] and then in stool from Nigerian, Tunisian, Nepalese and US children [5] . HBoV4 was reported in the stool of children with AFP from Nigeria and Tunisia [5] . HBoV 1/2/3 were also detected in untreated sewage water from throughout the US [26] . Recently, a novel bocavirus was identified in the feces of captive gorillas with diarrhea [6] and from wild gorillas and chimpanzees [7] . In 2009, a porcine bocalike virus (PBo-likeV) was reported in swine feces with postweaning multisystemic wasting syndrome in Sweden and 1854 bp of its partial genome sequenced [27] . In 2010, the nearly complete genomes of distinct porcine bocaviruses provisionally named PBoV1 and PBoV2 were characterized from feces of swine in China [28] . Finally, partial genome sequences of 2.4 Kb from another clade of porcine bocaviruses labeled 6V and 7V where also identified yielding three major bocavirus groups in swine (PBo-likeV, PBoV1/PBoV2, and 6V/7V). Random amplification and sequencing has been used to discover novel virus in human and animal [8, 10, . In this study, we found in swine feces highly distinct bocavirus whose genome we tentatively named porcine bocavirus 1 (PBoV1-H18), and porcine bocavirus 2 (PBoV 2-A6). The nearly complete genomes of both viruses were acquired and are described here. PBoV1-H18 and PBoV2-A6 were also screened for in 340 stool samples of asymptomatic swine from five provinces of China. A total of 340 porcine stool samples from different aged swine were collected from 17 middle or large-scale porcine farms (200-2,000 sows each) in five provinces of China from April 2008 to October 2009, of which 80 were collected from four farms located in Shanghai, 60 from three farms located in Jiangsu province, 120 from six farms located in Anhui province, 20 from a farm located in Shandong province, and 60 from three farms located in Guizhou province and stored at 280uC (Table 1 ). One stool sample was randomly selected from each of the 17 farms. The samples were suspended in PBS (0.01 M phosphate, pH7.2-7.4, 0.15NaCl), vortexed, centrifuged at 15000 g, and filtered through a 0.22-mm filter to remove eukaryotic-and bacterial-cell-sized particles [45, 46, 54, 55] . The filtrates were then treated with benzonase, DNase and RNase to digest non-particleprotected nucleic acid as reported [54, 55] . Viral nucleic acids were then extracted using the TIANamp virus DNA/RNA Kit (TIANGEN BIOTECH, BEIJING, CO., Ltd.). Viral cDNA synthesis was performed by incubation of the extracted viral RNA/DNA with 100 pmol of primer K-8N [56] with a degenerate 39 end and the use of Superscript reverse transcriptase, and the opposite strand of the cDNA was generated after melting and reannealing and primer extension using Klenow DNA polymerase [45, 46, [54] [55] [56] . PCR of extension products was performed as reported previously using the K-8 primer (K-8N without the degenerate 39 end) [54] . This protocol amplifies both viral RNA and DNA genomes [54, 55] . Random RT-PCR DNA products ran as smears on agarose gel and were gel purified (Axygen, CA, USA), then subcloned into pMD-18T plasmid vector (TaKaRa, Japan) for sequencing. The sequences were then screened for sequence similarities using tBLASTx and BLASTn against the nr database in GenBank. Sequences of each PCR product were assembled using SeqMan II program (DNASTAR, Inc). The identification of open reading The near-full genomes of PBoV1-H18, PBoV2-A6 and the partial NS1 and VP1 sequences from the diagnostic nPCR have been deposited in GenBank under accession numbers HQ291308-HQ291309 and HQ291310-HQ291343. Seventeen porcine samples stool supernatants from 17 farms were analyzed using a generic viral particle-protected nucleic acid enrichment procedure followed by random amplification of extracted RNA and DNA (see materials and methods) [46, [54] [55] [56] . Amplified DNA was then subcloned and 1190 plasmid inserts were sequenced (70 for each of 17 samples). Nine samples (totaling 82 subclones) showed the presence of fragments whose virtual translation products were related to canine and human bocaviruses using BLASTx. Twenty-four clones from pig sample H18 and 26 clones from pig sample A6 showed significant similarity with bocaviruses. H18 derived sequences showed high (.99%) identity with the recently described porcine bocavirus-like virus (PBo-LikeV), the only porcine bocavirus reported at the time of these experiments (GenBank GU902971) [27] . Sequences from porcine sample A6 showed low identity with those of H18 and PBo-likeV. The H18 and A6 samples were selected for targeted viral genome amplification and sequencing. The 24 sequences from H18 were assembled to form a continuous sequence of approximately 2700 nucleotides that appeared to lack .1300 and 1000 nucleotides from the 59and 39 ends of its genome. PCR primers based on the available H18 sequences and regions highly conserved between HBoV2 (FJ170278) and canine bocavirus (FJ214110) were used to amplify the nearly complete bocavirus genome we provisionally called PBoV1-H18 (5267 nt). To confirm this genome sequence, this sequence was re-amplified using 6 sets of PCR primers generating overlapping fragments of the genome which were directly sequenced. Using the same method, the nearly complete bocavirus genome (5117 nt) from sample A6 was also sequenced and was provisionally labeled PBoV2-A6. Using an open reading frame (ORF) finder (http://www.ncbi. nlm.nih.gov/gorf/gorf.html), three ORF were found in both genomes (Figure 1 ). The ORFs of PBoV1-H18 were 636 aa for NS1, 219 aa for NP1 and 621 aa for VP1/VP2. The ORFs of PBoV2-A6 were 703 aa for NS1, 221 aa for NP1 and 704 aa for VP1/VP2. The possible splicing of bocavirus NS1 transcripts recently shown to extend the length of NS1 proteins was not investigated here [6, 57] . Nucleic acids were extracted from 340 porcine stool samples. PBoV1-H18 related sequences were screened for using nested PCR with primers amplifying a 530-bp fragment of the NP1/VP1 region. The electrophoretic bands of the expected size were subcloned and sequenced. The results showed that the prevalence of PBoV1-H18 related viruses was high in China with 215 out of 340 (63.2%) porcine samples positive ( To determine the genetic relationship of PBoV1-H18 and PBoV2-A6 with recently described porcine bocaciruses and bocaviruses from other host species, both nucleotide and amino acid alignments were generated and used for phylogenetic analyses. When the whole genomes were considered porcine bocaviruses as a group (except for the V6/V7 variants with only NP1 sequences available), were most closely related to the canine bocavirus CnMV (Figure 2A ). Phylogenetic analyses of the 3 ORFs -NS1, NP1 and VP1/VP2 -were also performed ( Figure 2B-G) . In all three regions, PBoV1-H18 was most closely related to the Chinese PBoV1 and PBoV2 recently reported by Cheng et al [28] . PBoV2-A6 was closely related in NP1 to the first reported partial porcine bocavirus sequence PBo-likeV from Sweden, the only region available for comparison (Fig. 2C, F ) [27] . Table 2 numerically shows the protein similarities between PBoV1-H18 and PBoV2-A6 and other porcine and non-porcine bocaviruses. The partial VP1 sequence of PBoV2-A6 related variants from different farms was also phylogentically analyzed and fell into two major clades we named PBoV2 genotype 1 and 2 (PBoV2-G1 and PBoV2-G2). The two previously described Chinese PBoV1 and PBoV2 ''species'' grouped within these two genotypes. PBo-likeV was originally found in swine with postweaning multisystemic wasting syndrome (PMWS) in 2009 when approximately 35% of its genome sequence was reported [27] . The nearly full genomes of two distantly related porcine bocaviruses labeled PBoV1 and PBoV2 as well as two partial genomes labeled V6 and V7 were then reported in 2010 [28] . These bocaviruses grouped into 3 phylogenetic clades containing PBo-likeV, PBoV1/PBoV2, and V6/V7. In the present study, we characterize two nearly complete bocavirus genomes, one of which (PBoV1-H18) grouped with the PBo-likeV clade while the second (PBoV2-A6) fell with the PBoV1/PBoV2 clade. We provisionally named the genome from the H18 sample PBoV1-H18 since its closest homologue, PBo-likeV, was the first reported porcine bocavirus [27] . The virus from sample A6 was provisionally named PBoV2-A6 since it phylogenetically groups with the second reported set of porcine bocaviruses containing both PBoV1 and PBoV2 [28] . No close homologues of the V6/V7 clade were identified in this study. Under this proposed classification scheme, the viruses labeled PBoV1 and PBoV2 by Cheng et al, therefore both belong to the PBoV2 clade, the second reported clade of porcine bocaviruses [28] . Under this proposed taxonomic classification, the V6 and V7 bocaviruses [28] belong to the still only partially characterized PBoV3 clade. Sequence analysis of the PBoV2 clade showed that their VP1/ VP2 genes were highly diverse and could be classified into two genotypes ( Figure 3 ). The partial NP1 and VP1 genes of PBoV1-H18 related viruses detected in this study were more highly conserved (99-100% identity), consistent with a recent report by Zhai et al reporting PBoV1 in Chinese pigs using partial VP1/ VP2 nested PCR and sequencing [58] . Zhai et al also found a high prevalence of PBoV1 (69% in weanling piglets) with a higher frequency of PBoV1 in animals also infected with PCV2, PRRSV, PTTV or CSFV and in pigs with respiratory symptoms versus healthy pigs [58] . Members of the bocavirus genus contain an 3 rd ORF of unknown function labeled NP1 gene. Recently another parvovirus (PPV4) was identified in porcine feces that also contained a central 3 rd ORF, although unrelated in sequence to the bocaviruses NP1 [59] . None of the ORFs of PPV4 clustered phylogenetically with the bocaviruses (data not shown) but instead clustered with members of the parvovirus genus [59] . PPV4 is therefore unrelated to the bocaviruses reported here. HBoV1 has been associated with respiratory symptoms while other HBoV may be associated with diarrhea and acute flaccid paralysis [5, [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] 25] . A gorilla bocavirus was detected in captive animals in the US experiencing severe diarrhea [6] . PBoV1 was found in pigs with PMWS in Sweden [27] and in pigs with respiratory tract symptoms in China [58] . In the present study, both PBoV1 and PBoV2 were highly prevalent in both asymptomatic swine from five provinces of China. Further studies are needed to examine possible associations between infections with these different porcine bocaviruses, the viral loads excreted, the presence of co-infections and various porcine diseases. Virus Taxonomy: The Eighth Report of the International Committee on Taxonomy of Viruses Animal bocaviruses: a brief review Complete nucleotide sequence and genome organization of bovine parvovirus Recovery and characterization of a minute virus of canines Human bocaviruses are highly diverse, dispersed, recombination prone, and prevalent in enteric infections Identification and characterization of a new bocavirus species in gorillas Widespread infection of chimpanzees and gorillas with homologues of human parvoviruses B19, PARV4 and human bocavirus in the wild Cloning of a human parvovirus by molecular screening of respiratory tract samples A novel bocavirus associated with acute gastroenteritis in Australian children A newly identified bocavirus species in human stool Human bocavirus Human bocavirus and acute wheezing in children Detection of human bocavirus in Canadian children in a 1-year study Human bocavirus infection in hospitalized children during winter Bocavirus infection in hospitalized children Human bocavirus: a novel parvovirus epidemiologically associated with pneumonia requiring hospitalization in Thailand Human bocavirus respiratory infections in children Clinical and molecular epidemiology of human bocavirus in respiratory and fecal samples from children in Hong Kong Human bocavirus infection in children with gastroenteritis Role of Human Bocavirus infections in outbreaks of gastroenteritis Human bocavirus-a novel parvovirus to infect humans No gastroenteric Bocavirus in high risk patients stool samples Human bocavirus infection in children with acute gastroenteritis and healthy controls Human bocavirus in an immunocompromised child presenting with severe diarrhea The first detection of human bocavirus 2 infections in China Frequent detection of highly diverse variants of cardiovirus, cosavirus, bocavirus, and circovirus in sewage samples collected in the United States Detection of a novel porcine boca-like virus in the background of porcine circovirus type 2 induced postweaning multisystemic wasting syndrome Identification and nearly full-length genome characterization of novel porcine bocaviruses Identification of a third human polyomavirus A virus discovery method incorporating DNase treatment and its application to the identification of two bovine parvovirus species Cardioviruses are genetically diverse and cause common enteric infections in South Asian children Genetic detection and characterization of Lujo virus, a new hemorrhagic feverassociated arenavirus from Southern Africa Identification of cardioviruses related to Theiler's murine encephalomyelitis virus in human infections A metagenomic survey of microbes in honey bee colony collapse disorder Viral metagenomics Clonal integration of a polyomavirus in human Merkel cell carcinoma Complete genome sequence of a highly divergent astrovirus isolated from a child with acute diarrhea Identification of a novel polyomavirus from patients with acute respiratory tract infections The complete genome of klassevirus-a novel picornavirus in pediatric stool Klassevirus 1, a previously undescribed member of the family Picornaviridae, is globally widespread Novel borna virus in psittacine birds with proventricular dilatation disease New adenovirus species found in a patient presenting with gastroenteritis New DNA viruses identified in patients with acute viral infection syndrome Discovery of a novel human picornavirus in a stool sample from a pediatric patient presenting with fever of unknown origin A highly divergent picornavirus in a marine mammal A highly prevalent and genetically diversified Picornaviridae genus in South Asian children Recovery of divergent avian bornaviruses from cases of proventricular dilatation disease: identification of a candidate etiologic agent Genomic characterization of novel human parechovirus type A new arenavirus in a cluster of fatal transplant-associated diseases Characterization of a novel coronavirus associated with severe acute respiratory syndrome Identification of a novel gammaretrovirus in prostate tumors of patients homozygous for R462Q RNASEL variant A newly discovered human pneumovirus isolated from young children with respiratory tract disease Rapid identification of known and new RNA viruses from animal tissues Metagenomic analyses of viruses in stool samples from children with acute flaccid paralysis Characterization of virus isolates by particle-associated nucleic acid PCR Characterization of the gene expression profile of human bocavirus High prevalence of a novel porcine bocavirus in weanling piglets with respiratory tract symptoms in China Identification and molecular cloning of a novel porcine parvovirus