Identification and nearly full-length genome characterization of novel porcine bocaviruses

The genus bocavirus includes bovine parvovirus (BPV), minute virus of canines (MVC), and a group of human bocaviruses (HBoV1-4). Using sequence-independent single primer amplification (SISPA), a novel bocavirus group was discovered with high prevalence (12.59%) in piglet stool samples. Two nearly full-length genome sequences were obtained, which were approximately 5,100 nucleotides in length. Multiple alignments revealed that they share 28.7-56.8% DNA sequence identity with other members of Parvovirinae. Phylogenetic analyses indicated their closest neighbors were members of the genus bocavirus. The new viruses had a putative non-structural NP1 protein, which was unique to bocaviruses. They were provisionally named porcine bocavirus 1 and 2 (PBoV1, PBoV2). PBoV1 and PBoV2 shared 94.2% nucleotide identity in NS1 gene sequence, suggesting that they represented two different bocavirus species. Two additional samples (6V, 7V) were amplified for 2,407 bp and 2,434 bp products, respectively, including a partial NP1 gene and the complete VP1 gene; Phylogenetic analysis indicated that 6Vand 7V grouped with PBoV1 and PBoV2 in the genus of bocavirus, but were in the separate clusters. Like other parvoviruses, PBoV1, PBoV2, 6Vand 7V also contained a putative secretory phospholipase A(2) (sPLA(2)) motif in the VP1 unique region, with a conserved HDXXY motif in the catalytic center. The conserved motif YXGXF of the Ca(2+)-binding loop of sPLA2 identified in human bocavirus was also found in porcine bocavirus, which differs from the YXGXG motif carried by most other parvoviruses. The observation of PBoV and potentially other new bocavirus genus members may aid in molecular and functional characterization of the genus bocavirus.     

猪博卡病毒环状附加体的发现与鉴定

为了证实博卡病毒可以环状附加体形式存在于宿主体内,本研究利用半巢式PCR方法在健康猪粪便标本中筛查出2株猪博卡病毒环状附加体PBoV_G2-episome和PBoV_G3-episome。通过反复测序和序列拼接得到其末端非编码区序列（405nt和511nt）。经过对其进行序列分析及二级结构的预测,发现PBoV_G2-episome与人博卡病毒3附加体（HBoV3-episome）结构相似,而PBoV_G3-episome与博卡病毒属其他成员的末端二级结构存在较大差别。猪博卡病毒环状附加体的发现证实了有些博卡病毒与其他细小病毒的复制方式存在一定差异,也为今后博卡病毒感染性克隆的构建提供了一条研究思路。     

Ungulate bocaparvovirus 4 and rodent bocavirus are different genotypes of the same species of virus

Bocaviruses are associated with many human infectious diseases, such as respiratory tract infections, gastroenteritis, and hepatitis. Rats are known to be reservoirs of bocaviruses, including rodent bocavirus and rat bocavirus. Recently, ungulate bocaparvovirus 4, a known porcine bocavirus, has also been found in rats. Thus, investigating bocaviruses in rats is important for determining the origin of the viruses and preventing and controlling their transmission. To the best of our knowledge, no study to date has investigated bocaviruses in the livers of rats. In this report, a total of 624 rats were trapped in southern China between 2014 and 2017. Liver and serum samples from rats were tested for the prevalence of bocaviruses using PCR. Sequences related to ungulate bocaparvovirus 4 and rodent bocavirus were detected in both liver and serum samples. Interestingly, the prevalence of ungulate bocaparvovirus 4 (reference strain:KJ622366.1) was higher than that of rodent bocavirus (reference strain:KY927868.1) in both liver (2.24% and 0.64%, respectively) and serum samples (2.19% and 0.44%, respectively). The NS1 regions of ungulate bocaparvovirus 4 and rodent bocavirus related sequences displayed over 84% and 88% identity at the nucleic acid and amino acid levels, respectively. Furthermore, these sequences had similar genomic structure, genomic features, and codon usage bias, and shared a common ancestor. These viruses also displayed greater adaptability to rats than pigs. Our results suggested that ungulate bocaparvovirus 4 and rodent bocavirus may originate from rats and may be different genotypes of the same bocavirus species.     

Complete genome sequence of a novel species of Porcine Bocavirus, PBoV5

Porcine bocavirus 5 is a novel porcine bocavirus species found in a pig with clinical diarrhea from a farm in China. Here, we report the complete genome sequence of strain PBoV5/JS677, which will help toward understanding the molecular and evolutionary characteristics of the porcine bocavirus.     

The fecal virome of pigs on a high-density farm

Swine are an important source of proteins worldwide but are subject to frequent viral outbreaks and numerous infections capable of infecting humans. Modern farming conditions may also increase viral transmission and potential zoonotic spread. We describe here the metagenomics-derived virome in the feces of 24 healthy and 12 diarrheic piglets on a high-density farm. An average of 4.2 different mammalian viruses were shed by healthy piglets, reflecting a high level of asymptomatic infections. Diarrheic pigs shed an average of 5.4 different mammalian viruses. Ninety-nine percent of the viral sequences were related to the RNA virus families Picornaviridae, Astroviridae, Coronaviridae, and Caliciviridae, while 1% were related to the small DNA virus families Circoviridae, and Parvoviridae. Porcine RNA viruses identified, in order of decreasing number of sequence reads, consisted of kobuviruses, astroviruses, enteroviruses, sapoviruses, sapeloviruses, coronaviruses, bocaviruses, and teschoviruses. The near-full genomes of multiple novel species of porcine astroviruses and bocaviruses were generated and phylogenetically analyzed. Multiple small circular DNA genomes encoding replicase proteins plus two highly divergent members of the Picornavirales order were also characterized. The possible origin of these viral genomes from pig-infecting protozoans and nematodes, based on closest sequence similarities, is discussed. In summary, an unbiased survey of viruses in the feces of intensely farmed animals revealed frequent coinfections with a highly diverse set of viruses providing favorable conditions for viral recombination. Viral surveys of animals can readily document the circulation of known and new viruses, facilitating the detection of emerging viruses and prospective evaluation of their pathogenic and zoonotic potentials.     

WLL-1株博卡病毒(Bocavirus)全基因组序列分析

儿童下呼吸道感染已成为当前儿科发病率最高的一种疾病,而病毒是小儿下呼吸道感染的重要原因。临床上有相当比例的小儿下呼吸道感染病因未能作出明确的实验室诊断,给临床诊断与治疗带来了较大的困难。小儿下呼吸道感染可以由不同病毒引起,其中包括呼吸道合胞病毒、流感病毒、腺     

Porcine Bocavirus: A 10-Year History since Its Discovery

Virologica Sinica - Porcine bocavirus (PBoV) is a single-stranded DNA virus, belongs to the genus Bocaparvovirus of family Parvoviridae. It was discovered along with porcine circovirus 2 (PCV 2)...     

Identification and Characterization of a New Bocavirus Species in Gorillas

A novel parvovirus, provisionally named Gorilla Bocavirus species 1 (GBoV1), was identified in four stool samples from Western gorillas (Gorilla gorilla) with acute enteritis. The complete genomic sequence of the new parvovirus revealed three open reading frames (ORFs) with an organization similar to that of known bocaviruses. Phylogenetic analysis using complete capsid and non structural (NS) gene sequence suggested that the new parvovirus is most closely related to human bocaviruses (HBoV). However, the NS ORF is more similar in length to the NS ORF found in canine minute virus and bovine parvovirus than in HBoV. Comparative genetic analysis using GBoV and HBoV genomes enabled characterization of unique splice donor and acceptor sites that appear to be highly conserved among all four HBoV species, and provided evidence for expression of two different NS proteins in all primate bocaviruses. GBoV is the first non-human primate bocavirus identified and provides new insights into the genetic diversity and evolution of this highly prevalent and recently discovered group of parvoviruses.     

Whole Genomic Analysis of Human G12P[6] and G12P[8] Rotavirus Strains that Have Emerged in Myanmar

G12 rotaviruses are emerging rotavirus strains causing severe diarrhea in infants and young children worldwide. However, the whole genomes of only a few G12 strains have been fully sequenced and analyzed. In this study, we sequenced and characterized the complete genomes of six G12 strains (RVA/Human-tc/MMR/A14/2011/G12P[8], RVA/Human-tc/MMR/A23/2011/G12P[6], RVA/Human-tc/MMR/A25/2011/G12P[8], RVA/Human-tc/MMR/P02/2011/G12P[8], RVA/Human-tc/MMR/P39/2011/G12P[8], and RVA/Human-tc/MMR/P43/2011/G12P[8]) detected in six stool samples from children with acute gastroenteritis in Myanmar. On whole genomic analysis, all six Myanmarese G12 strains were found to have a Wa-like genetic backbone: G12-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 for strains A14, A25, P02, P39, and P43, and G12-P[6]-I1-R1-C1-M1-A1-N1-T1-E1-H1 for strain A23. Phylogenetic analysis showed that most genes of the six strains examined in this study were genetically related to globally circulating human G1, G3, G9, and G12 strains. Of note is that the NSP4 gene of strain A23 exhibited the closest relationship with the cognate genes of human-like bovine strains as well as human strains, suggesting the occurrence of reassortment between human and bovine strains. Furthermore, strains A14, A25, P02, P39, and P43 were very closely related to one another in all the 11 gene segments, indicating derivation of the five strains from a common origin. On the other hand, strain A23 consistently formed distinct clusters as to all the 11 gene segments, indicating a distinct origin of strain A23 from that of strains A14, A25, P02, P39, and P43. To our knowledge, this is the first report on whole genome-based characterization of G12 strains that have emerged in Myanmar. Our observations will provide important insights into the evolutionary dynamics of spreading G12 rotaviruses in Asia.     

Epistatic Roles of E2 Glycoprotein Mutations in Adaption of Chikungunya Virus to Aedes Albopictus and Ae. Aegypti Mosquitoes

Between 2005 and 2007 Chikungunya virus (CHIKV) caused its largest outbreak/epidemic in documented history. An unusual feature of this epidemic is the involvement of Ae. albopictus as a principal vector. Previously we have demonstrated that a single mutation E1-A226V significantly changed the ability of the virus to infect and be transmitted by this vector when expressed in the background of well characterized CHIKV strains LR2006 OPY1 and 37997. However, in the current study we demonstrate that introduction of the E1-A226V mutation into the background of an infectious clone derived from the Ag41855 strain (isolated in Uganda in 1982) does not significantly increase infectivity for Ae. albopictus. In order to elucidate the genetic determinants that affect CHIKV sensitivity to the E1-A226V mutation in Ae. albopictus, the genomes of the LR2006 OPY1 and Ag41855 strains were used for construction of chimeric viruses and viruses with a specific combination of point mutations at selected positions. Based upon the midgut infection rates of the derived viruses in Ae. albopictus and Ae. aegypti mosquitoes, a critical role of the mutations at positions E2-60 and E2-211 on vector infection was revealed. The E2-G60D mutation was an important determinant of CHIKV infectivity for both Ae. albopictus and Ae. aegypti, but only moderately modulated the effect of the E1-A226V mutation in Ae. albopictus. However, the effect of the E2-I211T mutation with respect to mosquito infections was much more specific, strongly modifying the effect of the E1-A226V mutation in Ae. albopictus. In contrast, CHIKV infectivity for Ae. aegypti was not influenced by the E2-1211T mutation. The occurrence of the E2-60G and E2-211I residues among CHIKV isolates was analyzed, revealing a high prevalence of E2-211I among strains belonging to the Eastern/Central/South African (ECSA) clade. This suggests that the E2-211I might be important for adaptation of CHIKV to some particular conditions prevalent in areas occupied by ECSA stains. These newly described determinants of CHIKV mosquito infectivity for Ae. albopictus and Ae. aegypti are of particular importance for studies aimed at the investigation of the detailed mechanisms of CHIKV adaptations to its vector species.     

B-Cell Responses to Human Bocaviruses 1–4: New Insights from a Childhood Follow-Up Study

Human bocaviruses (HBoVs) 1–4 are recently discovered, antigenically similar parvoviruses. We examined the hypothesis that the antigenic similarity of these viruses could give rise to clinically and diagnostically important immunological interactions. IgG and IgM EIAs as well as qPCR were used to study ~2000 sera collected from infancy to early adolescence at 3–6-month intervals from 109 children whose symptoms were recorded. We found that HBoV1-4-specific seroprevalences at age 6 years were 80%, 48%, 10%, and 0%, respectively. HBoV1 infections resulted in significantly weaker IgG responses among children who had pre-existing HBoV2 IgG, and vice versa. Furthermore, we documented a complete absence of virus type-specific immune responses in six viremic children who had pre-existing IgG for another bocavirus, indicating that not all HBoV infections can be diagnosed serologically. Our results strongly indicate that interactions between consecutive HBoV infections affect HBoV immunity via a phenomenon called “original antigenic sin”, cross-protection, or both; however, without evident clinical consequences but with important ramifications for the serodiagnosis of HBoV infections. Serological data is likely to underestimate human exposure to these viruses.     

Complete genome analysis of three novel picornaviruses from diverse bat species

Although bats are important reservoirs of diverse viruses that can cause human epidemics, little is known about the presence of picornaviruses in these flying mammals. Among 1,108 bats of 18 species studied, three novel picornaviruses (groups 1, 2, and 3) were identified from alimentary specimens of 12 bats from five species and four genera. Two complete genomes, each from the three picornaviruses, were sequenced. Phylogenetic analysis showed that they fell into three distinct clusters in the Picornaviridae family, with low homologies to known picornaviruses, especially in leader and 2A proteins. Moreover, group 1 and 2 viruses are more closely related to each other than to group 3 viruses, which exhibit genome features distinct from those of the former two virus groups. In particular, the group 3 virus genome contains the shortest leader protein within Picornaviridae, a putative type I internal ribosome entry site (IRES) in the 5'-untranslated region instead of the type IV IRES found in group 1 and 2 viruses, one instead of two GXCG motifs in 2A, an L→V substitution in the DDLXQ motif in 2C helicase, and a conserved GXH motif in 3C protease. Group 1 and 2 viruses are unique among picornaviruses in having AMH instead of the GXH motif in 3C(pro). These findings suggest that the three picornaviruses belong to two novel genera in the Picornaviridae family. This report describes the discovery and complete genome analysis of three picornaviruses in bats, and their presence in diverse bat genera/species suggests the ability to cross the species barrier.     

Genetic diversity within human erythroviruses: identification of three genotypes

B19 virus is a human virus belonging to the genus Erythrovirus: The genetic diversity among B19 virus isolates has been reported to be very low, with less than 2% nucleotide divergence in the whole genome sequence. We have previously reported the isolation of a human erythrovirus isolate, termed V9, whose sequence was markedly distinct (>11% nucleotide divergence) from that of B19 virus. To date, the V9 isolate remains the unique representative of a new variant in the genus Erythrovirus, and its taxonomic position is unclear. We report here the isolation of 11 V9-related viruses. A prospective study conducted in France between 1999 and 2001 indicates that V9-related viruses actually circulate at a significant frequency (11.4%) along with B19 viruses. Analysis of the nearly full-length genome sequence of one V9-related isolate (D91.1) indicates that the D91.1 sequence clusters together with but is notably distant from the V9 sequence (5.3% divergence) and is distantly related to B19 virus sequences (13.8 to 14.2% divergence). Additional phylogenetic analysis of partial sequences from the V9-related isolates combined with erythrovirus sequences available in GenBank indicates that the erythrovirus group is more diverse than thought previously and can be divided into three well-individualized genotypes, with B19 viruses corresponding to genotype 1 and V9-related viruses being distributed into genotypes 2 and 3.     

Genetically Diverse Low Pathogenicity Avian Influenza A Virus Subtypes Co-Circulate among Poultry in Bangladesh

Influenza virus surveillance, poultry outbreak investigations and genomic sequencing were assessed to understand the ecology and evolution of low pathogenicity avian influenza (LPAI) A viruses in Bangladesh from 2007 to 2013. We analyzed 506 avian specimens collected from poultry in live bird markets and backyard flocks to identify influenza A viruses. Virus isolation-positive specimens (n = 50) were subtyped and their coding-complete genomes were sequenced. The most frequently identified subtypes among LPAI isolates were H9N2, H11N3, H4N6, and H1N1. Less frequently detected subtypes included H1N3, H2N4, H3N2, H3N6, H3N8, H4N2, H5N2, H6N1, H6N7, and H7N9. Gene sequences were compared to publicly available sequences using phylogenetic inference approaches. Among the 14 subtypes identified, the majority of viral gene segments were most closely related to poultry or wild bird viruses commonly found in Southeast Asia, Europe, and/or northern Africa. LPAI subtypes were distributed over several geographic locations in Bangladesh, and surface and internal protein gene segments clustered phylogenetically with a diverse number of viral subtypes suggesting extensive reassortment among these LPAI viruses. H9N2 subtype viruses differed from other LPAI subtypes because genes from these viruses consistently clustered together, indicating this subtype is enzootic in Bangladesh. The H9N2 strains identified in Bangladesh were phylogenetically and antigenically related to previous human-derived H9N2 viruses detected in Bangladesh representing a potential source for human infection. In contrast, the circulating LPAI H5N2 and H7N9 viruses were both phylogenetically and antigenically unrelated to H5 viruses identified previously in humans in Bangladesh and H7N9 strains isolated from humans in China. In Bangladesh, domestic poultry sold in live bird markets carried a wide range of LPAI virus subtypes and a high diversity of genotypes. These findings, combined with the seven year timeframe of sampling, indicate a continuous circulation of these viruses in the country.     

Heterologous Prime-Boost Vaccination with MF59-Adjuvanted H5 Vaccines Promotes Antibody Affinity Maturation towards the Hemagglutinin HA1 Domain and Broad H5N1 Cross-Clade Neutralization

In an open label clinical study (2007), MF59-adjuvanted hemagglutinin (HA) vaccine from H5N1-A/Vietnam/1194/2004 (clade 1) was administered to subjects previously vaccinated (primed) with clade 0 H5N3 (A/duck/Singapore/97) vaccine at least 6 years earlier (in 1999 or 2001). The primed individuals responded rapidly and generated high neutralizing antibody titers against the H5N1-Vietnam strain within 7 days of a single booster vaccination. Furthermore, significant cross-neutralization titers were measured against H5N1 clade 0, 1, and 2 viruses. In the current study, the impact of MF59 adjuvant during heterologous priming on the quality of humoral polyclonal immune response in different vaccine arms were further evaluated using real time kinetics assay by surface plasmon resonance (SPR). Total anti-H5N1 HA1 polyclonal sera antibody binding from the heterologous prime-boost groups after a single MF59-H5N1 boost was significantly higher compared with sera from unprimed individuals that received two MF59-H5N1 vaccinations. The antigen-antibody complex dissociation rates (surrogate for antibody affinity) of the polyclonal sera against HA1 of H5N1-A/Vietnam/1194/2004 from the MF59-H5N3 primed groups were significantly higher compared to sera from unadjuvanted primed groups or unprimed individuals that received two MF59-H5N1 vaccines. Furthermore, strong inverse correlations were observed between the antibody dissociation off-rates of the immune sera against HA1 (but not HA2) and the virus neutralization titers against H5 vaccine strains and heterologous H5N1 strains. These findings supports the use of oil-in-water-adjuvanted pandemic influenza vaccines to elicit long term memory B cells with high affinity BCR capable of responding to potential variant pandemic viruses likely to emerge and adapt to human transmissions.     

Reassortment of Human and Animal Rotavirus Gene Segments in Emerging DS-1-Like G1P[8] Rotavirus Strains

The emergence and rapid spread of novel DS-1-like G1P[8] human rotaviruses in Japan were recently reported. More recently, such intergenogroup reassortant strains were identified in Thailand, implying the ongoing spread of unusual rotavirus strains in Asia. During rotavirus surveillance in Thailand, three DS-1-like intergenogroup reassortant strains having G3P[8] (RVA/Human-wt/THA/SKT-281/2013/G3P[8] and RVA/Human-wt/THA/SKT-289/2013/G3P[8]) and G2P[8] (RVA/Human-wt/THA/LS-04/2013/G2P[8]) genotypes were identified in fecal samples from hospitalized children with acute gastroenteritis. In this study, we sequenced and characterized the complete genomes of strains SKT-281, SKT-289, and LS-04. On whole genomic analysis, all three strains exhibited unique genotype constellations including both genogroup 1 and 2 genes: G3-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2 for strains SKT-281 and SKT-289, and G2-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2 for strain LS-04. Except for the G genotype, the unique genotype constellation of the three strains (P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2) is commonly shared with DS-1-like G1P[8] strains. On phylogenetic analysis, nine of the 11 genes of strains SKT-281 and SKT-289 (VP4, VP6, VP1-3, NSP1-3, and NSP5) appeared to have originated from DS-1-like G1P[8] strains, while the remaining VP7 and NSP4 genes appeared to be of equine and bovine origin, respectively. Thus, strains SKT-281 and SKT-289 appeared to be reassortant strains as to DS-1-like G1P[8], animal-derived human, and/or animal rotaviruses. On the other hand, seven of the 11 genes of strain LS-04 (VP7, VP6, VP1, VP3, and NSP3-5) appeared to have originated from locally circulating DS-1-like G2P[4] human rotaviruses, while three genes (VP4, VP2, and NSP1) were assumed to be derived from DS-1-like G1P[8] strains. Notably, the remaining NSP2 gene of strain LS-04 appeared to be of bovine origin. Thus, strain LS-04 was assumed to be a multiple reassortment strain as to DS-1-like G1P[8], locally circulating DS-1-like G2P[4], bovine-like human, and/or bovine rotaviruses. Overall, the great genomic diversity among the DS-1-like G1P[8] strains seemed to have been generated through reassortment involving human and animal strains. To our knowledge, this is the first report on whole genome-based characterization of DS-1-like intergenogroup reassortant strains having G3P[8] and G2P[8] genotypes that have emerged in Thailand. Our observations will provide important insights into the evolutionary dynamics of emerging DS-1-like G1P[8] strains and related reassortant ones.     

New Genotype of Dengue Type 3 Virus Circulating in Brazil and Colombia Showed a Close Relationship to Old Asian Viruses

Dengue type 3 genotype V viruses have been recently detected in Brazil and Colombia. In this study, we described another Brazilian isolate belonging to this genotype. Phylogenetic analysis including dengue type 3 viruses isolated worldwide showed that Brazilian and Colombian viruses were closely related to viruses isolated in Asia more than two decades ago. The characteristic evolutionary pattern of dengue type 3 virus cannot explain the close similarity of new circulating viruses with old viruses. Further studies are needed to confirm the origin of the new dengue type III genotype circulating in Brazil and Colombia.     

Bocavirus episome in infected human tissue contains non-identical termini

Human bocaviruses (HBoV) are highly prevalent human infections whose pathogenic potential remains unknown. Recent identification of the first non-human primate bocavirus [1] in captive gorillas raised the possibility of the persistent nature of bocavirus infection. To characterize bocavirus infection in humans, we tested intestinal biopsies from 22 children with gastrointestinal disease for the presence of HBoV DNA. Four HBoV-positive tissue samples were analyzed to determine whether viral DNA was present in the linear genomic, the episomal closed circular or the host genome-integrated form. Whereas one tissue sample positive for HBoV3 contained the episomal form (HBoV3-E1), none had the genome-integrated form. The complete genome sequence of HBoV3-E1 contains 5319 nucleotides of which 513 represent the non-coding terminal sequence. The secondary structure of HBoV3-E1 termini suggests several conserved and variable features among human and animal bocaviruses. Our observation that HBoV genome exists as head-to-tail monomer in infected tissue either reflects the likely evolution of alternative replication mechanism in primate bocaviruses or a mechanism of viral persistence in their host. Moreover, we identified the HBoV genomic terminal sequences that will be helpful in developing reverse genetic systems for these widely prevalent parvoviruses. Significance: HBoV have been found in healthy human controls as well as individuals with respiratory or gastrointestinal disease. Our findings suggest that HBoV DNA can exist as episomes in infected human tissues and therefore can likely establish persistent infection in the host. Previous efforts to grow HBoV in cell culture and to develop reverse genetic systems have been unsuccessful. Complete genomic sequence of the HBoV3 episome and its genomic termini will improve our understanding of HBoV replication mechanism and its pathogenesis.