key: cord-280188-gir0y1m1 authors: Wang, Yanqun; Li, Yamin; Liu, Jun; Zhao, Yanjie; Xie, Zhengde; Shen, Jun; Tan, Wenjie title: Genetic characterization of human bocavirus among children with severe acute respiratory infection in China date: 2016-06-13 journal: J Infect DOI: 10.1016/j.jinf.2016.05.014 sha: doc_id: 280188 cord_uid: gir0y1m1 OBJECTIVES: To investigate the genetic character of Human bocavirus (HBoV) among children with severe acute respiratory infection (SARI) in China. METHODS: We screened 993 respiratory samples for HBoV by PCR among hospitalized children with SARI between September 2007 and March 2014. Four of HBoV1 samples were selected for complete genomes analysis by next-generation sequencing. RESULTS: The results show that 200 (20.1%) out of 993 samples were HBoV-positive, most of these HBoV belong to HBoV1 subtype (n = 197), HBoV2 (n = 1) and HBoV3 (n = 2) were also detected. Fifty (5.04%) of 993 SARI patient were detected as HBoV-positive only. Four HBoV1 genomes in this study were conserved and showed no significant difference among the nucleotide diversity from different regions. Analyses of evolutionary rates showed that NS1 exhibited the highest degree of conservation while the VP1 gene exhibited the fastest rate of evolution at 4.20 × 10(−4) substitutions/site/year. The nucleotide deletions and substitutions occurred in NP1 and VP1 represented novel molecular signatures enabling subtype differentiation between HBoVs. CONCLUSIONS: We described some new characteristics in the epidemiology of HBoV among children with SARI, these data will significantly expand the current knowledge of HBoV epidemic and genomic characterization among children with SARI. Human bocavirus1e4 (HBoV) represents a novel pathogen associated with gastrointestinal and respiratory tract illnesses. 1e3 According to the latest ICTV classification of parvovirus, HBoV1 and HBoV3 belonged to Primate bocaparvovirus 1 species; HBoV2 and HBoV4 are part of P. bocaparvovirus 2 species. 4 The genome of HBoV is w5.3 kb in length, and is divided into four partially overlapping genes, namely NS1, NP1, VP1, and VP2, VP2 is totally included within VP1. 5, 6 While the prototype HBoV (HBoV1) were first discovered from nasopharyngeal samples, 7 three additional viruses (HBoV2e4) have since been discovered in stool specimens, classified based upon their close phylogenetic relationship with HBoV1. 8, 9 HBoV1 is the most commonly reported genotype and occurring primarily in pediatric patients with respiratory tract infection, but also gastrointestinal symptoms are often observed. 10, 11 In contrast, HBoV2 are preferentially detected in stool samples and appear to be more strongly associated with enteric disease, HBoV3 and HBoV4 are occasionally detected in faeces and too rare for any associations. 2 Since the discovery of HBoVs, numerous epidemiological surveillance efforts examining HBoV prevalence in children have been performed across multiple regions, comprising Thailand, United States, France, Jordan and Brazil. 10e14 Severe acute respiratory infection (SARI) is among the leading causes of morbidity and mortality among children globally. 15 HBoV infection in children has been reported associated with respiratory tract infection, a few cases reported that HBoV as the cause of SARI among Children. 16e18 However, prolonged shed periods of HBoV and high coinfection detection resulted in the on-going debate of the HBoV as the agent of SARI 19 ; In addition, the comprehensive research of HBoV genome among children with SARI were limited, especially in China. To better understand the molecular epidemiology and characterization of HBoV genome in children with SARI, we investigated the prevalence of HBoV among 993 inpatient children with SARI in China. Analyses of genome characterization were also performed. From Sep 2007 to Mar 2014, a total of 993 nasopharyngeal aspirates (NPAs) or induced sputum (IS) were randomly collected from hospitalized children with SARI in Beijing (n Z 259), Shanghai (n Z 441) and Zhejiang (n Z 293) area. The case of SARI was defined according to the World Health Organization case definition for all hospitalized children in whom the onset of illness occurred within seven days of admission. Most of the patients had received clinical diagnosis of respiratory tract infection, including pneumonia, acute bronchitis/bronchiolitis, asthma exacerbation and acute pharyngitis. The most common respiratory symptoms included fever (temperature !38 C), cough, sore throat, shortness of breath, vomit, dyspnea and so on. 15 In addition, none of the samples come from the patients in Pediatric intensive care unit (PICU) and most of the children have no co-morbidities, such as heart and liver diseases. All the samples were collected by medical professionals and placed in a tube containing of viral transportation medium and stored at À80 C. This project was approved by the Research Ethics Committee of Beijing Children's Hospital, Children's Hospital of Fudan University, Wenzhou Medical College, and the Institutional Review Board at the China CDC. Written informed consent was obtained on the participants' behalf from their parents or guardians. Viral nucleic acids were extracted from virus transport medium by the QIAamp MinElute Virus Spin Kit (QIAGEN, Germany), according to the instructions provided by the manufacturer. As previously described, 8 partial VP1/VP2 gene fragment was amplified by nested PCR to screen and type HBoV infection. The first round-PCR primers were F1 (5 0 -CGCCGTGGCTCCTGCTCT-3 0 ) and R1 (5 0 -TGTTCGCCATCA-CAAAAGATGTG-3 0 ) with 609 bp product, the second-round PCR primers were F2 (5 0 -GGCTCCTGCTCTAGGAAATAAA-GAG-3 0 ) and R2 (5 0 -CCTGCTGTTAGGTCGTTGTTGTATGT-3 0 ) with 576 bp PCR product. Positive products were cloned into pMDT-18 vector and sequenced by ABI 3730xl automated sequencer. HBoV co-infection with other respiratory viruses, including HRSV, HRV/EV, HAdV, HMPV, HPIV1e4, influenza A/B virus and HCoVs (-OC43, -229E, -NL63, -HKU1), was also screened as described previously. 20, 21 Sequencing and phylogenetic analysis of HBoV1 genome Four samples of HBoV1 infection only in this study were used for complete genomes sequencing by next-generation sequencing. Samples were pretreated as previously, 22 and the amplified DNA was used as a template for Illumina Hiseq 2500 sequencing, paired-end reads (2 Â 125 bp reads) were assembled into contigs by CLC genomic workbench. To analyze genetic variation of HBoV detected, nucleotide sequences were compared to strains available from GenBank. Nucleotide sequence alignment was conducted through MAFFT version 5. 23 Phylogenetic and molecular evolutionary analyses were constructed by Neighbor-Joining Method using MEGA 5.0 24 with the bootstrap value of 1000. Evolutionary rates were calculated using the Bayesian Markov Chain Monte Carlo approach employed by BEAST version 1.7. 25 Fragments of the NS1, NP1 and VP1 genes of HBoV were aligned separately and used to calculate the rate of nucleotide substitutions/site/year, under an uncorrelated lognormal-relaxed clock model of rate variation among lineages. The best-fit models were selected by jMo-delTest, with the following models: HKY þ I þ G for NP1, GTR þ G for NS1 and GTR þ I þ G for VP1. The evolutionary rates of individual genes were then compared to identify the differences in conservation throughout the genome. DNA SP5 software 26 was used to analyze the nucleic acid sequence diversity of the HBoV genomes. The model of HBoV1 VP1 (NC_007455) was based on the crystal structure of a capsid viral protein of Adenoassociated virus (PDB code: 4IOV) 27 which shared a close genetic relationship with HBoVs. Molecular modeling was performed manually using the COOT software under the guidance of Fo-Fc and 2Fo-Fc electron density maps. Consequently, we refined the initial rigid body, and performed a series of restrained TLS refinements using Refmac5. Additional rounds of refinements were performed using the phenix refine program implemented in the Phenix package with isotropic ADP refinement and bulk solvent modeling. We assessed the stereochemical quality of the final models with the program PROCHECK. Eventually, the molecular model was generated using PyMOL (http://www.pymol.org/). The nucleotide sequences generated in this study have been deposited in GenBank under the accession numbers KM378039wKM378094, KM877548wKM877592, KM8877594wKM877614 and KM464728wKM464730. Data were analyzed by the chi-squared test using SAS software version 9.2. P < 0.05 was considered statistically significant. All the respiratory samples were collected from hospitalized children with SARI in Beijing, Shanghai and Zhejiang from 2007 to 2014. The sample information was provided in Table 1 . The median age of the population in each of the three regions studied is 7 months, 7.5 months and 1 year, respectively. No significant differences were observed regarding gender and age distribution (P > 0.05) and the characteristics of three regions' population was matched well. In all, 200 (20.1%) of 993 SARI patient were positive for HBoV and the infection rates of HBoV among inpatient in Beijing, Zhejiang and Shanghai area were 21.6% (56/ 259), 22.5% (66/293) and 17.7% (78/441), respectively. No significant difference was observed regarding infection rate (P > 0.05). Phylogenetic analysis based on partial VP1/VP2 sequences indicated that HBoV1 was the most prevalent in China. The strains isolated from Beijing were identified as HBoV1 (n Z 54), HBoV2 (n Z 1) and HBoV3 (n Z 1). Similarly, of the 66 HBoV-positive samples collected from Zhejiang Province, 65 were HBoV1, one was HBoV3 (Fig. 1) . In addition, all the HBoV-positive samples from Shanghai city belonged to HBoV1. Combined with the three cohorts, the most frequently detected strains in China was HBoV1 (98.5%, [197/200] ), followed by HBoV3 (1%, 2/200) HBoV2 (Table S1 ). Through high-throughput sequencing, four complete genomes of HBoV were obtained (KM464728, KM464729, Figure 1 Phylogenetic analysis of HBoV detected among Children with SARI. The phylogenetic tree was constructed based on partial VP1/VP2 gene sequence. All the sequences presented here were indicated in black solid ball and Reference strains of HBoV were indicated in red solid triangle. Sequences were aligned by Neighbor-Joining method with 1000 bootstrap replicates using MEGA 5.0. KM464730 and KM464731), which all belong to HBoV1 (Fig. 2) . All four genomes of HBoV1 sequenced in this study were shown to be closely related to each other at the nucleotide (nt) level (99.47%e99.74% nt identity for each genome) and to be closely related (99.60%e99.87%) to prototypes HBoV1 (NC_007455.1). A comparative analysis of the four HBoV1 genomes presented here and all available other HBoV genomes identified a total of three deletions within NP1 and VP1 together, which could be used to distinguishing HBoV1 from HBoV2e4. These differences include a 3-bp deletion in VP1 of HBoV1 at nucleotides 419e421, a 12-bp deletion in VP1 of HBoV2e4 at nucleotides 993e996 and 1004e1011, and a 3-bp deletion in NP1 of HBoV2e4 at nucleotides 88e90 (Fig. 3) . To further investigate the potential implications of VP1 variations in HBoV1, relative to HBoV2e4, we generated a 3D model of HBoV1 VP1 based on the crystal structure of a capsid viral protein of Adeno-associated virus (PDB code: 4IOV) (Fig. 4) . According to the previous report, 27 the structure of VP1 contains two parts: the capsid exterior and capsid interior. The unique substitutions and insertions (residues 334e337) within HBoV1 compared to HBoV2e4 mapped primarily to the capsid exterior. Interestingly, among the 32 unique substitutions, 19 were related with serine or threonine (substitutions of S/T with other amino acids), which has a dramatic influence on the hydrophobicity of the protein. These consensus nucleotide deletions or substitutions can be used to distinguish HBoV1 from HBoV2e4 and may be associated with pathogenicity of HBoV. To determine the differences in gene-specific mutation rates, we calculated and compared the mutation rates of NS1, NP1 and VP1 of HBoV1 prevalent in China. The evolutionary rates of NS1, NP1 and VP1 of HBoV1 were 2.840 Â 10 À5 , 3.917 Â 10 À4 and 4.204 Â 10 À4 substitutions/ site/year with the ESS values >200, respectively ( Table 2 ). Comparative analysis showed that the evolutionary rate of VP1 is much faster than that of NS1 for HBoV1, and HBoV evolved relatively slowly in China. No genetic recombination event was found among HBoV1 from China. In addition, analysis of DNA polymorphisms within HBoV1 and HBoV2 was performed using the DNA SP5 software. The overall mean of diversities (Pi) were 0.00371(HBoV1 in China), 0.00395(HBoV1 out of China), 0.00377 (HBoV1 in this study) and 0.03728(HBoV2), respectively. There was no significant difference among the nucleotide diversity of HBoV1 from different regions and the nucleotide sequences diversity of VP1 gene was greater than others. In addition, the high degree of nucleotide sequence diversity in HBoV2 is greater than that of HBoV1 as expected (Fig. 5 ). HBoV is a novel parvovirus associated with respiratory tract infections in infants or children, including four genotypes (HBoV1-HBoV4) . 8, 28 Among the four recognized HBoV genotypes, HBoV1 have more often been associated with SARI. The prevalence of different HBoV genotype varies among the same region. Only limited data is available on the prevalence of HBoV infection in SARI children. Here, we performed a comprehensive molecular epidemiological study of HBoV in China between 2007 and 2014 and made genomic characterization analyses. Our study indicated that HBoV was frequently detected in SARI children from 2007 to 2014 in China; the frequency of HBoV1 was much higher than that of other HBoV subtypes. However, the infection rates of HBoV among hospitalized children with SARI in three different regions were almost consistent (21.6%, 22.5% and 17.7%). Previous reports demonstrate the prevalence of HBoV varies considerably between 1.5% and 19% in children with acute respiratory tract infections (ARTIs). 29, 30 However, the prevalence of HBoV among children with SARI was more common (200/993, 20.1%) in this study. Furthermore, fifty (5.04%) of 993 SARI patient were detected as HBoV positive only, about 25% among all 200 HBoV positive patients with SARI. These data indicated that the HBoV infection may play an important role among children with SARI, although co-infection of HBoV with other viruses was much higher than that of HBoV infection alone. In addition other factors may also affect the infection rate of HBoV, including assay sensitivity, specificity, sampling time and sampling locations, more respiratory samples were needed to provide more detailed information about the prevalence of HBoV among the children with SARI. We calculated the evolutionary rates of individual HBoV1 genes prevalent in China, identifying strong conservation in NS1 (2.840 Â 10 À5 substitutions/site/year), compared with VP1, which exhibited substantially higher mutation rates (4.204 Â 10 À4 substitutions/site/year). Furthermore, the evolutionary rate of HBoV1 VP1 is more than 10 time than that of NS1, indicating the strongest degree of conservation in non-structural protein NS1; while VP1, which encodes the capsid protein, mutated significantly more rapidly than other genes. However, we only calculated the evolutionary rates of HBoV1 genes, rather than HBoV2e4 genes, due to genomic recombinant 31, 32 and limited number of genome sequence. More extensive studies will be needed to address the evolution of HBoVs. Comparative analysis of the HBoV1 genome sequences presented here revealed consistent and reproducible nucleotide deletions and substitution. Three sets of deletions within NP1 and VP1 were shown to be diagnostic for HBoV1, clearly differentiating these from those of HBoV2e4. This difference is consistent with phenotypic analyses, which show HBoV1 to be the most commonly occurring in respiratory tract samples, while HBoV2e4 are detected mainly in gastrointestinal samples and are presumably enteric. Furthermore, the unique substitutions within HBoV1 relatively to HBoV2e4 mapped primarily to the capsid exterior. These otherwise uncommon substitutions and insertions within HBoV1 VP1 are likely to play a major role in the antigenicity of HBoV species and may account for differences in tissue tropism between HBoV1 and HBoV2e4. Further studies are necessary to validate these findings and to confirm the effects of these nucleotide differences on tissue tropism and pathogenicity. In summary, we first reported the circulating HBoV genotype and their genome character among children with SARI in China, providing more evidence for a causal role of HBoV1 in SARI. Novel molecular signatures were identified for distinguishing HBoV1 from other HBoV subtypes. This study complements and significantly expands upon the current knowledge of HBoV infection and SARI among children in China. Nil. Human bocavirus and acute wheezing in children A novel bocavirus associated with acute gastroenteritis in Australian children Human bocavirus, a respiratory and enteric virus The family Parvoviridae Genomic features of the human bocaviruses Human bocavirus capsid structure: insights into the structural repertoire of the parvoviridae Cloning of a human parvovirus by molecular screening of respiratory tract samples Human bocaviruses are highly diverse, dispersed, recombination prone, and prevalent in enteric infections A newly identified bocavirus species in human stool Circulating of human bocavirus 1, 2, 3, and 4 in pediatric patients with acute gastroenteritis in Thailand Human bocavirus infection in young children in the United States: molecular epidemiological profile and clinical characteristics of a newly emerging respiratory virus Human bocavirus in French children Human bocavirus infection among children Human bocavirus species 2 and 3 in Brazil Global burden of respiratory infections due to seasonal influenza in young children: a systematic review and meta-analysis Human bocavirus (HBoV) in Thailand: clinical manifestations in a hospitalized pediatric patient and molecular virus characterization Human bocavirus infection as a cause of severe acute respiratory tract infection in children Genetic variability of human metapneumo-and bocaviruses in children with respiratory tract infections Correlation of viral load of respiratory pathogens and co-infections with disease severity in children hospitalized for lower respiratory tract infection Viral etiology and clinical profiles of children with severe acute respiratory infections in China Prevalence of human parvovirus B19, bocavirus, and PARV4 in blood samples from the general population of China and lack of a correlation between parvovirus and hepatitis B Co-Infection Identification of diverse alphacoronaviruses and genomic characterization of a novel severe acute respiratory syndrome-like coronavirus from bats in china MAFFT multiple sequence alignment software version 7: improvements in performance and usability MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0 BEAST: Bayesian evolutionary analysis by sampling trees DnaSP v5: a software for comprehensive analysis of DNA polymorphism data The structure of AAVrh32. 33, a novel gene delivery vector Human bocavirus-1 primary infection and shedding in infants Surveillance and genome analysis of human bocavirus in patients with respiratory infection in Human bocavirus: prevalence and clinical spectrum at a children's hospital Recombination analysis based on the complete genome of bocavirus Prevalence analysis of different human bocavirus genotypes in pediatric patients revealed intra-genotype recombination Supplementary data related to this article can be found at http://dx.doi.org/10.1016/j.jinf.2016.05.014.