- Research
- Open Access
- Published:
Molecular epidemiology analysis of symptomatic and asymptomatic norovirus infections in Chinese infants
Virology Journal volume 20, Article number: 60 (2023)
Abstract
Background
Norovirus is a leading cause of acute gastroenteritis among children. Previous studies based on symptomatic infections indicated that mutations, rather than recombination drove the evolution of the norovirus ORF2. These characteristics were found in hospital-based symptomatic infections, whereas, asymptomatic infections are frequent and contribute significantly to transmission.
Methods
We conducted the first norovirus molecular epidemiology analysis covering both symptomatic and asymptomatic infections derived from a birth cohort study in the northern China.
Results
During the study, 14 symptomatic and 20 asymptomatic norovirus infections were detected in 32 infants. Out of the 14 strains that caused symptomatic infections, 12 strains were identified as GII.3[P12], and others were GII.4[P31]. Conversely, 17 asymptomatic infections were caused by GII.4[P31], two by GII.2[P16], and one by GII.4[P16]. Regardless of symptomatic and asymptomatic infections, the mutations were detected frequently in the ORF2 region, and almost all recombination were identified in the RdRp-ORF2 region. The majority of the mutations were located around the predefined epitope regions of P2 subdomain indicating a potential for immune evasion.
Conclusion
The role of symptomatic as well as asymptomatic infections in the evolution of norovirus needs to be evaluated continuously.
Background
Norovirus a the nonenveloped, single-stranded positive-sense RNA virus belonging to the Caliciviridae family [1], is a leading cause for acute gastroenteritis (AGE) in human hosts worldwide, especially among children and the elderly [2,3,4]. The whole genome of norovirus is 7.7 kb, which can be organized into three open reading frames (ORF 1–3) [5]. Of these, ORF2 encodes the major capsid protein (VP1). X-ray crystallography of VP1 showed two structural domains, including the N-terminal shell (S) domain and the C-terminal protrusion (P) domain [6]. The P domain is further divided into the P1 subdomain and the variable P2 subdomain which contains putative neutralization sites and interacts with histo-blood group antigens (HBGAs) [7].
Norovirus continues to evolve in humans through both mutation and recombination, in particular, recombination appears to be a major force driving virus evolution [8]. Novel norovirus strains have emerged, and some classified into new tentative genogroups and genotypes [9,10,11]. In the past years, we also had studied the evolution of norovirus open reading frame 2 (ORF2), utilizing strains derived from population-based surveillance, and detected far more mutations than recombination events in RNA-dependent RNA polymerase (RdRp) region of ORF1 [12]. These evolutionary characteristics of norovirus are mostly based on symptomatic infections. Yet the fecal excretion of norovirus infection is common in asymptomatic individuals, especially in children [13]. Asymptomatic infections, constitute a significant reservoir for infection in the community, and may act as a source of endemic and epidemic disease [14, 15]. Our study characterized norovirus strains in samples collected weekly from symptomatic and asymptomatic infants participating in a birth cohort study. A phylogenetic analysis sequence encoding the norovirus RdRp and ORF2 was performed to classify the strains co-circulating during the surveillance period. The mutation and recombination patterns of norovirus strains were explored by recombination analysis.
Materials and methods
Study population and study design
To understand the transmission dynamics of norovirus in Chinese infants, a birth cohort study was conducted in 11 villages of Zhengding County, Hebei Province, China, which is located 270 km south of Beijing. Considering the potential interference of passively transferred maternal antibodies [16], infants aged between 3-months (90-days) and 4-months (120-days) were recruited. Exclusion criteria included serious medical conditions such as any neonatal disorder, kidney, liver, lung and/or heart disease, congenital disorders, clinically diagnosed enteropathy, and/or birth weight < 1.5 kg. Recruited infants were followed weekly for occurrence of AGE by trained doctors working in the village clinics. Stool specimens were collected from each infant once every week, and tested for norovirus using RT-PCR, regardless of the clinical manifestations of AGE. An AGE episode was defined as 3 or more loose bowel movements and/or 2 or more episodes of vomiting during a 24-h period [17]. Recovery of AGE was defined as the absence of loose stools and/or vomiting for 3 consecutive days. A symptomatic norovirus infection was defined as norovirus test-positive by RT-PCR in the stool samples collected from infants with AGE within 7 days after the onset, and an asymptomatic infection was defined correspondingly as norovirus test-positive by RT-PCR in the stool samples collected from infants without AGE [18].
This study was approved by the Institutional Review Board (IRB) of the Institutes of Biomedical Sciences, Fudan University. Written inform consent was obtained from a parent/guardian of each child included in the study. All CRFs were double entered into a custom-made data entry program (EpiData program, version 3.1) [19], the data management program included error and consistency checks.
Norovirus RNA extraction and genotyping
Viral RNA was extracted from 10% of the stool supernatant using TianLong Stool DNA/RNA Extraction Kit (TianLong Science & Technology, China), as described previously [20]. The G (capsid type) and P (polymerase type) typing were carried out by amplifying parts of the ORF2 and ORF1 sequence, respectively. PCR assays was carried out using Qiagen One-Step RT-PCR kit (Qiagen, USA), with primer-pairs MON432/G1SKR for norovirus GI, or MON431/G2SKR for norovirus GII, generating 579-bp and 570-bp PCR products [21]. The assay was carried out in a 25 μl reaction mixture, containing 12.5 μl One-Step RT-PCR 2X Master Mix, 3 μl RNA template, 1.5 μl of 10 μM forward and reverse primers, and sterile, de-ionized water were mixed into a final volume of 25 μl. The cycling conditions included 30 min of reverse transcription at 42 °C, and then 15 min of initial denaturation at 95 °C. PCR amplification was performed for 40 cycles (95 °C, 50 °C, and 72 °C for 1 min each), and finally, reactions were completed with 10 min elongation at 72 °C followed by cooling to 4 °C. PCR products were analyzed on a 1 × TAE 2% agarose gel and visualized with Gel-Red staining. Dideoxy sanger sequencing of the products was performed by Tsingke (TsingkeBiotech, China). Genotypes were determined using the online Norovirus Genotyping Tool (https://www.rivm.nl/mpf/typingtool/norovirus/). Once the genotype was assigned, PCR assays were carried out using Qiagen One-Step RT-PCR kit (Qiagen, USA) with genotype-specific primers to amplify and sequence the full-length of ORF2 (Additional file 1: Table S1).
Phylogenetic analysis
Phylogenetic trees were constructed using the neighbor-joining (NJ) method for each of the norovirus polymerase and capsid genes [22]. Multiple sequence alignments for phylogenetic analysis were performed with the Molecular Evolutionary Genetic Analysis software (MEGA, version 7.0.21) [23]. The best substitution model was selected based on the corrected Akaike’s information criterion (AICc) values [24]. Kimura-2 parameter models were selected, and gamma distributions were used to calculate genetic distance. Reliability analyses were performed using a bootstrap method. Sampling was repeated 1000 times, with less than 70% considered meaningless [25]. The sequences of norovirus strains used as reference were downloaded from GenBank database. The sequence of strains isolated from our study were available in the GenBank (GenBank accession nos. OQ451905–OQ451938).
Nucleotide variation and protein mutation analysis
We evaluate the mutation and recombination events of the norovirus strains isolated in our study. The mutation was evaluated through Shannon entropy, which could represent the diversity of the mutation on individual site. In this study, the Shannon entropy was calculated by Shannon Entropy-One, which could apply phylogeny into Shannon entropy as a measure of variation in nucleotide sequence alignments (https://www.hiv.lanl.gov/content/sequence/ENTROPY/entropy_one.html) [26]. Entropy-one calculates the entropy at each position in the input sequence and can be used as a measure of the relative variation in different positions or regions of an aligned nucleotide sequence [27]. Entropy values for each position were plotted in GraphPad Prism v8 (San Diego, USA).
The ORF2 protein sequences of GII.4[P31] strains were aligned using MEGA and analyzed using AliView viewer and editing tool (version 1.28) [28]. Linear, discontinuous, and non-peptidic norovirus GII.4 specific T and B cell receptor epitopes mapping in VP1 region were retrieved from the Immune Epitope Database (IEDB) [29]. Mutations were defined as amino acid changes presenting between two or more VP1 sequences of GII.4 strains. The effect of the resulting amino acid substitutions on the structure and function of capsid proteins was prognosticated using the Polymorphism Phenotyping v2 (PolyPhen-2) software [30]. The prediction is based on a number of sequences, phylogenetic, and structural features characterizing the substitution. For a given amino acid substitution in a protein, PolyPhen-2 extracts various sequence and structure-based features of the substitution site and feeds them to a probabilistic classifier [31].
For structure modelling, the online server of SWISS-MODEL was used to construct the 3D structure of the VP1 protein through protein homology modeling with defaulted template [32]. Then, the potential deleterious positions can be visualized and annotated on the.pbd format files through PyMOL (version 2.3.4) [33].
Recombination analysis
The identification of potential parental sequences and the localization of putative recombination breakpoints of norovirus strains isolated from our study were determined by the Recombination Detection Program version 4 (RPD4, version 101) [34]. Seven recombination detection methods (Bootscan, Chimaera, Geneconv, MaxChi, RDP, SISCAN, and 3Seq) in RPD4 were used to predict recombination events [35]. Recombination events were accepted when potential recombination signals (p-value cutoff of 0.05) detected by at least three out of the seven detection methods [36].
Results
Study population and characteristics
As of the end of October, 100 infants were recruited from 11 villages. Between November 1, 2021 and March 1, 2022, which was the epidemic season, 1,600 fecal specimens were collected from the 100 infants. Of these, 14 symptomatic and 20 asymptomatic episodes of norovirus infections were detected in 32 infants, resulting in a cumulative infection rate of 34%. Compared with November (23.5%) and February (17.6%), norovirus infection was identified more frequently in December (26.5%) and January (32.4%). Re-infection was detected in 2 infants with different genotypes. One infant was identified as an asymptomatic GII.4[P16] infection, and subsequently experienced a symptomatic GII.3[P12] infection after 46 days. Another infant presented with diarrhea caused by GII.3[P12] infection, and then experienced an asymptomatic GII.2[P16] infection after 52 days. Out of 14 strains causing symptomatic infections, 12 strains were identified as GII.3[P12], and other strains were GII.4[P31]. Conversely, 17 asymptomatic infections were attributed to GII.4[P31], two as GII.2[P16], and one as GII.4[P16] (Fig. 1). The duration of viral shedding in symptomatic infections (1.43 ± 0.51 weeks) was slightly longer than that of asymptomatic infections (1.35 ± 0.49 weeks). This difference was not statistically significant.
Phylogenetic analysis
To detect the relations between the previously worldwide circulating strains and the strains monitored through our cohort study, we constructed phylogenetic trees based on the partial sequences of ORF2 (1591nt) and RdRp (565nt) (Fig. 2). The majority of norovirus strains isolated from both symptomatic and asymptomatic individuals in our study belonged to norovirus GII.4 Sydney2012 (n = 20), including 19 GII.4[P31] and one GII.4[P16] (Fig. 2A). All norovirus strains with GII.4 Sydney2012 capsid type were clustered with the previously reported GII.4 strains from China, Japan, Philippine, Thailand, Korea and USA, creating a phylogenetically distinct monophyletic clade. The nucleotide sequence identities among the GII.4 Sydney2012 strains detected in this study range from 93.36 to 98.51%. The previously detected GII.4[P31] strains circulated from 2017 to 2022 in both northern and southern of China, including Beijing, Zhengzhou, Shanghai and Guangzhou, which were genetically close to the GII.4[P31] strains isolated in our study. Similar, the GII.4[P16] strain clustered with strains from China, Thailand and USA reported from 2016 to 2020. This finding suggests that this clade was imported between 2016 and 2019.
Phylogenetic analyses of norovirus isolated from symptomatic and asymptomatic infections in our cohort study during 2021–2022. A Phylogenetic tree based on the ORF2 nucleotide sequence B Phylogenetic tree based on the RdRp nucleotide sequence. Reference strains were downloaded from GenBank and indicated by GenBank accession number/strain/name/year/country (JP, Japan; CHN, China; USA, United States; TH, Thailand; KOR, Korea; AU, Australia; VN, Vietnam; RUS, Russia; CA, Canada). Strains obtained from our study are labeled as strain/sample name/source (AC, asymptomatic individuals; SC, symptomatic individuals)/year. GII.4 labeled as red circle; GII.2 labeled as orange circle; GII.3 labeled as blue circle; GII.P31 labeled as red triangle; GII.P12 labeled as blue triangle; GII.P16 labeled as green triangle. Neighbor-joining (NJ) phylogenetic tree was constructed with MEGA 7 software and bootstrap tests (1000 replicates), based on the Kimura-2 parameter model. Scale bar indicates distances between sequence pairs
Besides GII.4 clades, we detected the second most-common genotype of GII.3 (n = 12) and the third most-common genotype of GII.2 (n = 2) in this study. The norovirus GII.3 strains were 97.3%-99.8% identical to each other and clustered together on the phylogenetic tree. The GII.3 strains detected in 2021 were clustered together with the HNZZ1036 GII.3 strain detected in Zhengzhou, China in 2017. The strains detected in 2022 were most closely related to the norovirus strains isolated in 2016 and 2017 in Russia. Two norovirus GII.2 strains were isolated from asymptomatic infected children in 2021 and 2022. The nucleotide sequence identity of the GII.2 strains isolated in our study was 99.4%. The strains were located in the same branch as norovirus GII.2 reference strains reported in 2010 and 2013 from China, Japan and Australia.
We detected three major clades of GII.P31, GII.P12 and GII.P16 based on the phylogenetic tree of the RdRp sequence (Fig. 2B). The majority of norovirus strains isolated in our study located in the GII.P31 clade (n = 19), which clustered with strains from China and Japan reported from 2017 to 2022. This branch also clustered with strains from Asia, Northern America and Europe isolated since 2009. Similarly, the second GII.3[P12] strains from SJZ were mostly located in one branch with high sequence identity (n = 12). Two other strains clustered with strains detected in the US and China, one GII.4[P16] strain and two highly similar GII.2[P16] strains.
Nucleotide variation and protein mutation analysis
We evaluated the variation of norovirus strains isolated from our study at the nucleotide level. The Shannon entropy were illustrated in Fig. 3, including 20 strains of the GII.4 cluster based on ORF2 region (Fig. 3A) and 19 strains of the GII.P31 cluster based on RdRp region (Fig. 3B). The mean value of Shannon entropy spanning the whole GII.4 ORF2 region was 0.031 (range: 0–0.693), and the mean value across RdRp was 0.015 (range: 0–0.0681), suggesting that the polymerases proteins were relatively conserved. Results of other regions can be found in the Additional file 1: Fig. S1. As shown in the Fig. 3A, most variation were observed within ORF2, particularly in the 5900–6300 nucleotide sites, which located on the P2 subdomain.
To further determine whether nucleotide substitution led to amino acid variation, we analyzed the complete amino acid sequence of VP1 proteins of GII.4[P31] strains isolated from both symptomatic and asymptomatic infections in our study. The positions involving more than two mutated residues in different strains were considered significant mutation sites. Eight important mutation sites were detected in the VP1 protein. As shown in Table 1, 25.0% (2/8) amino acid changes were identified in the Shell domain of VP1 (aa site: 46–221), and 62.5% (5/8) were observed in the P2 subdomain (aa site: 276–417). The mutation site of V8A (1/8, 12.5%) was located in the N-terminal (aa site:1–45). No mutations were located in the P1 subdomain (aa site: 222–275 & 418–540).
In the shell domain, the S174P site was identified in the previously reported CD4 + T cell restricted epitopes 161VRNNFYHYNQSNDST175 (IEDB ID: 985331), and the mutation site I145V was reported in the CD8 + T cell restricted epitope, 134SPSQVTMFPHIIVDVRQL151 (IEDB ID: 561756). We detected the H414P mutation in both the CD8 + epitopes 410GRNTHNVHL418 (IEDB ID: 984324) and antibody Epitopes E (407, 411–414) in the P2 subdomain. G295C and S309N, occurred in the B cell restrict epitopes of 288GDVTHITGSRNYTMNLASQNWSNY311 (IEDB ID: 1334831). N373H was identified in the antibody Epitopes A (294–298, 368, 372–373) and G393S was identified in the antibody Epitopes D (393–397). Three mutations, V8A, G295C, H414P and were characterized in the PolyPhen-2 as damaging, and the rest were considered benign (Table 1). We mapped 5 previously determined epitope regions, labeled Epitope A to Epitope E on P2 subdomain (Fig. 4). By mapping the mutation sites on the P2 subdomain, we found most of the mutations were not directly located on the epitope regions, but close to them.
3D structure modeling of VP1 region of norovirus GII.4[P31] strain. Protein site information is referenced to GII.4 Sydney 2012 VP1 (AFV08795). Mutations identified in this study are shown in red. Epitopes A (294–298, 368, 372–373), B (333, 389), C (340, 376), D (393–395), and E (407, 411–414) are highlighted in cyan, magenta, green, orange, and blue, respectively. The 3D structures were generated with SWISS-MODEL and visualized using PyMOL (version 2.3.4)
Recombination analysis
Recombination is an alternative way for the norovirus to produce variants. We investigated the pattern of norovirus recombination by inclusion of both the norovirus strains identified from symptomatic and asymptomatic infants in cohort study. The results of recombination analysis can be found in Fig. 5 and the Additional file 1: Table S2. Two strains with full sequence length (Fig. 5A and B) and two strains with partial but complete RdRp/ORF2 sequence (Fig. 5C and 5D) were used for illustration. Results showed that the predicted recombination position is around 5020 to 5105, overlapped with the RdRp (4044–5097) and partial of ORF2 (5085–6707). Among all 34 isolated strains, the recombination position of 16 strains were located on RdRp, followed by 13 on the overlap region of RdRp and ORF2, and 5 on ORF2 (Additional file 1: Table S2). The results illustrated four possible recombinants GII.4[P31], GII.4[P16], GII.2[P16] and, GII.3[P12]. For GII.4[P31], the parental strains were Osaka_2007 (AB541319) for ORF1 and Apeldoorn_2008 (AB541268) for ORF2-ORF3, which are the common parental strains for other GII.4[P31] subtypes (Fig. 5A and Additional file 1: Table S2). The ORF1 parental strain for the GII.P16 strain was AY772730 strain from Germany in 2000, which has different recombination points from the GII.4[P16] strain of SJZ7125 (Recombination breakpoint site 5090) and the GII.2[P16] strain of SJZ8231 (Recombination breakpoint site 5068) (Fig. 5B). The ORF2-ORF3 parental strain for the GII.4[P16] strain of SJZ7125 was predicted as US95_96 (KC013592), which suggested a different origin than the GII.2[P16] strain of SJZ8231 with ORF2 origin from Japan in 2004 (DQ456824) (Fig. 5C). We found that all 12 GII.3[P12] strains originated from parental Norwalk-like strains isolate from Japan in 2000 (AB039775) and Mexican in 1995 (U22498) (Fig. 5D).
Genome recombination analysis of norovirus strains isolated from China in 2021–2022. The full length of the GII.4[P31] (SJZ1007) and GII.4[P16] (SJZ7125) sequence and partial of GII.2[P16] (SJZ8231) and GII.3[P12] (SJZ10110) were used for recombination analysis. RDP4 results show recombination breakpoints and relationships of representative strains to the parental strains. All plots were obtained using a 200 bp sliding window, 284 steps of 20 bp. Sequence locus information is referenced to GII.4 Sydney 2012 genome (JX459908). The horizontal line and vertical lines show nucleotide positions and nucleotide pair- wise identity. A Analysis of GII.4[P31] consensus sequence against Osaka_2007 strain (GenBank: AB541319, purple) and Apeldoorn_2008 strain (GenBank: AB541268, teal), predicted recombination position at 5081nt. B Analysis of GII.4[P16] consensus sequence against Berlin 2000 strain (GenBank: AY772730, purple) and US95_96 strain (GenBank: KC013592, teal), predicted recombination position at 5090nt. C Analysis of GII.2[P16] consensus sequence against Berlin 2000 strain (GenBank: AY772730, teal) and Tokyo 2004 strain (GenBank: DQ456824, purple), predicted recombination position at 5068nt. D Analysis of GII.3[P12] consensus sequence against Saitama 2002 strain (GenBank: AB039775, teal) and Mexican 1995 strain (GenBank: U22498, purple), predicted recombination position at 5039nt
Discussion
Norovirus is a leading cause of AGE, especially among children in China [37]. For a comprehensive understanding of virus evolution and transmission dynamics in symptomatic and asymptomatic infections, a birth cohort study is desirable [3]. Our previous study estimated that norovirus infection led to an annual incidence rate of 6.0, 15.6 and 5.5 per 100 persons per year in the three age groups a) all ages, b) less than 5Â years and c) older than 60Â years respectively [38]. We have studied the evolution of the major capsid protein (VP1) encoded by the ORF2 based on sequences isolated from population-based diarrhea surveillance in Zhengding county spanning between 2001 and 2019, and concluded that antigenic variation was the major mechanism for the emergence of novel VP1s (not including RdRp region), rather than recombination [12]. However, all of these findings were derived from symptomatic infections caused by norovirus. This study is the first birth cohort study to explore the transmission dynamics of norovirus in Chinese infants. More asymptomatic infections (20%) than symptomatic infections (14%), were detected, which suggested that the asymptomatic infections played an important role in the persistence of transmission, and thus, should be accounted when exploring the evolution of norovirus [39, 40]. Several factors might be involved in the mechanism for the high occurrence of asymptomatic infection in the community. Firstly, the infants recruited in our study were around 6Â months of age when they experienced the epidemic season. Some studies have demonstrated that the maternal antibody decay rapidly in the first months of life, although the measurable maternal antibodies were expected to have waned at 6Â months [41,42,43,44], and thus could still provide a certain level of protection, and eventually attenuated the severity of infections [45, 46]. Secondly, the study site is located in rural areas, where infants are generally breastfed up to 12Â months of age. It was reported that norovirus-specific immunoglobulin A in breast milk might protect against norovirus associated diarrhea but not norovirus infection [47]. Similar to some published studies [48, 49], symptomatic infections were not associated significantly with duration of viral shedding.
In this study, 14 strains that caused symptomatic infections were identified as GII.3[P12] and GII.4[P31], while 20 strains that caused asymptomatic infections were involving GII.4[P31], GII.2[P16] and GII.4[P16]. Among them, GII.4 Sydney2012 strains were detected in both AGE patients and asymptomatic carriers, which is in agreement with previous studies [50,51,52]. Since first detected in USA in 2012, GII.4 Sydney[P31] strains have been the predominated strain resulting in global pandemics [53]. The phylogenetic tree suggested that GII.4[P31] strains isolated from our study have a close genetic distance with previously circulating strains in Shanghai, Beijing, Zhengzhou and Guangzhou in 2016–2020. This finding proves that GII.4[P31] strains have circulated in the mainland China before 2016 [54, 55]. The GII.4[P16] strain in our study was similar to the viruses circulating in Beijing in 2019–2020, and strains detected in Thailand in 2019–2020. It has been reported that the GII.4 Sydney[P16] strain emerged in 2015 and was soon replaced by the GII.4 Sydney[P31] strain as the primary cause of outbreaks in the United States [56]. Despite the fact that GII.4 Sydney[P16] strains were rarely detected in our study, it remains necessary to monitor the prevalence of this strain. GII.3[P12] strains in our study clustered with strains from China, Japan, Russia and the United States in 2015–2018, and GII.2[P16] strains were clustered with China, Thailand and Japan in 2011–2013, suggesting that the strains were spread across the globe. We found that the GII.3[P12] strains were only detected in symptomatic infections, and GII.2[P16] strains were only detected in asymptomatic infections. This result may provide new perspective for the investigation of symptomatic and asymptomatic infections, but the results in this study are not sufficient to provide that the infection patterns are specifically distinguishable for those two genotypes. In the future, with the increasing accumulation of molecular epidemiological data, it may be possible for us to verify whether symptomatic and asymptomatic infections are related to the specific genotypes of norovirus.
The emergence and spread of novel norovirus strains is associated with point mutations in ORF1 and ORF2 region, and recombination events that produce chimeric viruses [57]. Although all norovirus utilized both mechanisms, different genotypes may preferentially emerge and persist in populations. Our study found that the mutation frequency of ORF2 region was higher than RdRp region of the ORF1, especially in the P2 subdomain. The amino acid mutations of the VP1 protein identified in our study (amino acid positions 309, 373, 414) corresponded to the same characteristics identified by another Chinese group, and the remaining amino acid mutation sites were located at or near the amino acids identified by this group [55]. Most of these mutations are located at or near the epitope region of the P2 sub-domain. It has been reported that residue changes of the epitope region are likely to enable the norovirus to escape the pressure of population immunity and cause global epidemic [58, 59]. The predicted recombination breakpoints in norovirus strains identified from symptomatic and asymptomatic infants in our study around 5020 to 5105nt, overlapped with the RdRp (4044–5097) and partially overlapped with ORF2 (5085–6707). This is consistent with previously reported strains isolated during gastroenteritis outbreaks [60]. RdRp of norovirus is a key enzyme responsible for viral transcription and replication and was suggested to be a driving factor in norovirus recombination [61]. Norovirus recombination typically occurs at the RdRp-ORF2 junction [62, 63], which is also the transcription initiation site for viral sub-genomic RNA. It has been shown that exchanges in RdRp region caused by genomic recombination result in an increased mutation rate through acquisition of an RdRp with lower fidelity and/or increased replicative ability, improving viral fitness under certain selective pressures and hence drive virus evolution [64]. Therefore, recombination may be one of the reason for the diversity of norovirus, which not only causing the genomic exchange of RdRp region, but also increasing the mutation rate in other regions such as VP1 [65]. Meanwhile, mutations in regions other than RdRp, may be caused by antigenic drift or shift that result in point mutations to escape immune response [57]. Together, the recombination and the point mutation may lead to the generation of new immune escape strains. Thus, large molecular epidemiological cohort studies are required to further verify the causes of norovirus diversity.
Conclusions
In conclusion, during the study period, more episodes of asymptomatic infection were detected in comparison with that of symptomatic infection. Genetic diversity between isolated strains of symptomatic and asymptomatic infection was not observed. The predominant circulating strains were GII.4, followed by GII.3 and GII.2, which origin possibly for recombination include GII.4[P31], GII.4[P16], GII.2[P16], and GII.3[P12] circulating in China, and other Asian countries in recent years. The role of asymptomatic infection in the evolution and transmission of noroviruses needs to be evaluated continuously.
Availability of data and materials
All data were collected from publicly available literatures, and all data generated or analyzed during this study are included in this published article and its additional files.
Abbreviations
- AGE:
-
Acute gastroenteritis
- RT-PCR:
-
Reverse transcription-polymerase chain reaction
- RdRp:
-
RNA-dependent RNA polymerase
References
Katayama K. Caliciviridae. Nihon Rinsho. 2003;61(Suppl 3):468–74.
Ahmed SM, Hall AJ, Robinson AE, Verhoef L, Premkumar P, Parashar UD, Koopmans M, Lopman BA. Global prevalence of norovirus in cases of gastroenteritis: a systematic review and meta-analysis. Lancet Infect Dis. 2014;14:725–30.
Cannon JL, Lopman BA, Payne DC, Vinje J. Birth cohort studies assessing norovirus infection and immunity in young children: a review. Clin Infect Dis. 2019;69:357–65.
Harris JP, Edmunds WJ, Pebody R, Brown DW, Lopman BA. Deaths from norovirus among the elderly, England and Wales. Emerg Infect Dis. 2008;14:1546–52.
Jiang X, Wang M, Wang K, Estes MK. Sequence and genomic organization of Norwalk virus. Virology. 1993;195:51–61.
Prasad BV, Hardy ME, Dokland T, Bella J, Rossmann MG, Estes MK. X-ray crystallographic structure of the Norwalk virus capsid. Science. 1999;286:287–90.
Tan M, Jiang X. Histo-blood group antigens: a common niche for norovirus and rotavirus. Expert Rev Mol Med. 2014;16: e5.
Ludwig-Begall LF, Mauroy A, Thiry E. Norovirus recombinants: recurrent in the field, recalcitrant in the lab - a scoping review of recombination and recombinant types of noroviruses. J Gen Virol. 2018;99:970–88.
Huang Z, Yao D, Xiao S, Yang D, Ou X. Full-genome sequences of GII.13[P21] recombinant norovirus strains from an outbreak in Changsha, China. Arch Virol. 2020;165:1647–52.
Wang C, Ao Y, Yu J, Xie X, Deng H, Jin M, Liu L, Duan Z. Complete genome sequence of a novel recombinant GII.P16-GII.1 norovirus associated with a gastroenteritis outbreak in Shandong Province, China, in 2017. Genome Announc. 2018;6:e01483-e1517.
Medici MC, Tummolo F, Martella V, Giammanco GM, De Grazia S, Arcangeletti MC, De Conto F, Chezzi C, Calderaro A. Novel recombinant GII.P16_GII.13 and GII.P16_GII.3 norovirus strains in Italy. Virus Res. 2014;188:142–5.
Zhou HL, Chen LN, Wang SM, Tan M, Qiu C, Qiu TY, Wang XY. Prevalence and evolution of noroviruses between 1966 and 2019 implications for vaccine design. Pathogens. 2021;10:1012.
Robilotti E, Deresinski S, Pinsky BA. Norovirus. Clin Microbiol Rev. 2015;28:134–64.
Phattanawiboon B, Nonthabenjawan N, Boonyos P, Jetsukontorn C, Towayunanta W, Chuntrakool K, Ngaopravet K, Ruchusatsawat K, Uppapong B, Sangkitporn S, et al. Norovirus transmission mediated by asymptomatic family members in households. PLoS ONE. 2020;15: e0236502.
Monica B, Ramani S, Banerjee I, Primrose B, Iturriza-Gomara M, Gallimore CI, Brown DW, Moses PD, Gray JJ, Kang G. Human caliciviruses in symptomatic and asymptomatic infections in children in Vellore, South India. J Med Virol. 2007;79:544–51.
Vielot NA, Brinkman A, DeMaso C, Vilchez S, Lindesmith LC, Bucardo F, Reyes Y, Baric RS, Ryan EP, Braun R, Becker-Dreps S. Breadth and dynamics of human norovirus-specific antibodies in the first year of life. J Pediatric Infect Dis Soc. 2022;11(10):463–6.
The World Health Organization. 2017. https://www.who.int/news-room/fact-sheets/detail/diarrhoeal-disease. Accessed 2 May 2017.
Huynen P, Mauroy A, Martin C, Savadogo LG, Boreux R, Thiry E, Melin P, De Mol P. Molecular epidemiology of norovirus infections in symptomatic and asymptomatic children from Bobo Dioulasso, Burkina Faso. J Clin Virol. 2013;58:515–21.
Singh S. Review of epidata entry and analysis freewares. Indian J Community Med. 2009;34:76–7.
Qiao N, Wang SM, Wang JX, Kang B, Zhen SS, Zhang XJ, Hao ZY, Ma JC, Qiu C, Zhao YL, et al. Variation analysis of norovirus among children with diarrhea in rural Hebei Province, north of China. Infect Genet Evol. 2017;53:199–205.
Chhabra P, Browne H, Huynh T, Diez-Valcarce M, Barclay L, Kosek MN, Ahmed T, Lopez MR, Pan CY, Vinje J. Single-step RT-PCR assay for dual genotyping of GI and GII norovirus strains. J Clin Virol. 2021;134: 104689.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Yamaoka K, Nakagawa T, Uno T. Application of Akaike’s information criterion (AIC) in the evaluation of linear pharmacokinetic equations. J Pharmacokinet Biopharm. 1978;6:165–75.
Hall BG. Building phylogenetic trees from molecular data with MEGA. Mol Biol Evol. 2013;30:1229–35.
Nemzer LR. Shannon information entropy in the canonical genetic code. J Theor Biol. 2017;415:158–70.
Tohma K, Lepore CJ, Martinez M, Degiuseppe JI, Khamrin P, Saito M, Mayta H, Nwaba AUA, Ford-Siltz LA, Green KY, et al. Genome-wide analyses of human noroviruses provide insights on evolutionary dynamics and evidence of coexisting viral populations evolving under recombination constraints. PLoS Pathog. 2021;17: e1009744.
Larsson A. AliView: a fast and lightweight alignment viewer and editor for large datasets. Bioinformatics. 2014;30:3276–8.
Fleri W, Paul S, Dhanda SK, Mahajan S, Xu X, Peters B, Sette A. The immune epitope database and analysis resource in epitope discovery and synthetic vaccine design. Front Immunol. 2017;8:278.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–9.
Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. 2013;Chapter 7:Unit7 20.
Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, Heer FT, de Beer TAP, Rempfer C, Bordoli L, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46:W296–303.
Janson G, Paiardini A. PyMod 3: a complete suite for structural bioinformatics in PyMOL. Bioinformatics. 2021;37:1471–2.
Martin DP, Murrell B, Golden M, Khoosal A, Muhire B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015;1:vev003.
Martin DP, Murrell B, Khoosal A, Muhire B. Detecting and analyzing genetic recombination using RDP4. Methods Mol Biol. 2017;1525:433–60.
Martin D, Rybicki E. RDP: detection of recombination amongst aligned sequences. Bioinformatics. 2000;16:562–3.
Wang LP, Zhou SX, Wang X, Lu QB, Shi LS, Ren X, Zhang HY, Wang YF, Lin SH, Zhang CH, et al. Etiological, epidemiological, and clinical features of acute diarrhea in China. Nat Commun. 2021;12:2464.
Zhou HL, Zhen SS, Wang JX, Zhang CJ, Qiu C, Wang SM, Jiang X, Wang XY. Burden of acute gastroenteritis caused by norovirus in China: a systematic review. J Infect. 2017;75:216–24.
Utsumi T, Lusida MI, Dinana Z, Wahyuni RM, Yamani LN, Juniastuti, Soetjipto, Matsui C, Deng L, Abe T, et al. Occurrence of norovirus infection in an asymptomatic population in Indonesia. Infect Genet Evol. 2017;55:1–7.
Bucardo F. Understanding asymptomatic norovirus infections. EClinicalMedicine. 2018;2–3:7–8.
Parker SP, Cubitt WD, Jiang X. Enzyme immunoassay using baculovirus-expressed human calicivirus (Mexico) for the measurement of IgG responses and determining its seroprevalence in London, UK. J Med Virol. 1995;46:194–200.
Parker SP, Cubitt WD, Jiang XJ, Estes MK. Seroprevalence studies using a recombinant Norwalk virus protein enzyme immunoassay. J Med Virol. 1994;42:146–50.
Vielot NA, Brinkman A, DeMaso C, Vilchez S, Lindesmith LC, Bucardo F, Reyes Y, Baric RS, Ryan EP, Braun R, Becker-Dreps S. Breadth and dynamics of human norovirus-specific antibodies in the first year of life. J Pediatric Infect Dis Soc. 2022;11:463–6.
Jing Y, Qian Y, Huo Y, Wang LP, Jiang X. Seroprevalence against Norwalk-like human caliciviruses in Beijing, China. J Med Virol. 2000;60:97–101.
Zhang S, Chen TH, Wang J, Dong C, Pan J, Moe C, Chen W, Yang L, Wang X, Tang H, et al. Symptomatic and asymptomatic infections of rotavirus, norovirus, and adenovirus among hospitalized children in Xi’an, China. J Med Virol. 2011;83:1476–84.
Zweigart MR, Becker-Dreps S, Bucardo F, Gonzalez F, Baric RS, Lindesmith LC. Serological humoral immunity following natural infection of children with high burden gastrointestinal viruses. Viruses. 2021;13:2023.
Labayo HKM, Pajuelo MJ, Tohma K, Ford-Siltz LA, Gilman RH, Cabrera L, Mayta H, Sanchez GJ, Cornejo AT, Bern C, et al. Norovirus-specific immunoglobulin A in breast milk for protection against norovirus-associated diarrhea among infants. EClinicalMedicine. 2020;27: 100561.
Newman KL, Moe CL, Kirby AE, Flanders WD, Parkos CA, Leon JS. Norovirus in symptomatic and asymptomatic individuals: cytokines and viral shedding. Clin Exp Immunol. 2016;184:347–57.
Teunis PF, Sukhrie FH, Vennema H, Bogerman J, Beersma MF, Koopmans MP. Shedding of norovirus in symptomatic and asymptomatic infections. Epidemiol Infect. 2015;143:1710–7.
Jeong MH, Song YH, Ju SY, Kim SH, Kwak HS, An ES. Surveillance to prevent the spread of norovirus outbreak from asymptomatic food handlers during the PyeongChang 2018 Olympics. J Food Prot. 2021;84:1819–23.
Wu QS, Xuan ZL, Liu JY, Zhao XT, Chen YF, Wang CX, Shen XT, Wang YX, Wang L, Hu Y. Norovirus shedding among symptomatic and asymptomatic employees in outbreak settings in Shanghai, China. BMC Infect Dis. 2019;19:592.
Ozawa K, Oka T, Takeda N, Hansman GS. Norovirus infections in symptomatic and asymptomatic food handlers in Japan. J Clin Microbiol. 2007;45:3996–4005.
Allen DJ, Adams NL, Aladin F, Harris JP, Brown DW. Emergence of the GII-4 norovirus Sydney 2012 strain in England, winter 2012–2013. PLoS ONE. 2014;9: e88978.
Fu J, Ai J, Bao C, Zhang J, Wu Q, Zhu L, Hu J, Xing Z. Evolution of the GII.3[P12] norovirus from 2010 to 2019 in Jiangsu, China. Gut Pathog. 2021;13:34.
Zhu X, He Y, Wei X, Kong X, Zhang Q, Li J, Jin M, Duan Z. Molecular epidemiological characteristics of gastroenteritis outbreaks caused by norovirus GII4 Sydney [P31] strains - China, October 2016–December 2020. China CDC Wkly. 2021;3:1127–32.
Lun JH, Hewitt J, Yan GJH, Enosi Tuipulotu D, Rawlinson WD, White PA. Recombinant GII.P16/GII.4 Sydney 2012 was the dominant norovirus identified in Australia and New Zealand in 2017. Viruses. 2018;10:548.
Parra GI. Emergence of norovirus strains: A tale of two genes. Virus Evol. 2019;5:vez048.
Lindesmith LC, Costantini V, Swanstrom J, Debbink K, Donaldson EF, Vinje J, Baric RS. Emergence of a norovirus GII.4 strain correlates with changes in evolving blockade epitopes. J Virol. 2013;87:2803–13.
Donaldson EF, Lindesmith LC, Lobue AD, Baric RS. Viral shape-shifting: norovirus evasion of the human immune system. Nat Rev Microbiol. 2010;8:231–41.
Eden JS, Tanaka MM, Boni MF, Rawlinson WD, White PA. Recombination within the pandemic norovirus GII.4 lineage. J Virol. 2013;87:6270–82.
Deval J, Jin Z, Chuang YC, Kao CC. Structure(s), function(s), and inhibition of the RNA-dependent RNA polymerase of noroviruses. Virus Res. 2017;234:21–33.
Bull RA, Tanaka MM, White PA. Norovirus recombination. J Gen Virol. 2007;88:3347–59.
Bull RA, Hansman GS, Clancy LE, Tanaka MM, Rawlinson WD, White PA. Norovirus recombination in ORF1/ORF2 overlap. Emerg Infect Dis. 2005;11:1079–85.
Smertina E, Urakova N, Strive T, Frese M. Calicivirus RNA-dependent RNA polymerases: evolution, structure, protein dynamics, and function. Front Microbiol. 2019;10:1280.
Mahar JE, Bok K, Green KY, Kirkwood CD. The importance of intergenic recombination in norovirus GII.3 evolution. J Virol. 2013;87:3687–98.
Acknowledgements
We thank the children and the parents for participating our study. We acknowledge the work of research assistants involved in enrollment of study children and samples collection in the local CDCs.
Funding
This work was supported by the National Science and Technology Major Projects for Significant New Drug Development (2018ZX09739002-006).
Author information
Authors and Affiliations
Contributions
S-MW, T-YQ, and X-YW contributed to the study’s conception and design. Z-YH, H-SZ, X-JZ, Y-LZ recruited participants and collected specimens. L-NC, S-JW, X-LF, W-JZ and CQ carried-out the laboratory analyses on samples. L-NC, and T-YQ conducted the bioinformatics analysis. L-NC, T-YQ, and X-YW wrote the paper. LS provided contributions in the writing of the paper. All authors provided important input to the study methods, revised the manuscript, and approved the final version. All authors had full access to all the study data and had final responsibility for the decision to submit for publication. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was reviewed and approved by the Institutional Review Board (IRB) of the Institutes of Biomedical Sciences, Fudan University (registration number: 2020-001). Written inform consent was obtained from a parent/guardian of each child included in the study. All experiments were performed in accordance with Fudan University Ethics Committee guidelines.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Table S1.
Genotype-specific primers were used to amplify and sequence the full-length of ORF2. Table S2. Analysis of predicted recombinant breakpoints of norovirus strains isolated from a cohort study, between 2021 and 2022. Fig. S1. Site variability was calculated at nucleotide level using Shannon entropy for norovirus strains isolated from our study. Diversity nucleotide plots, were shown the difference sites of the ORF2 and RdRp of norovirus GII.3[P12] and GII.2[P16] sequences. Sequence locus information is referenced to GII.4 Sydney2012 genome (JX459908)
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Chen, LN., Wang, SJ., Wang, SM. et al. Molecular epidemiology analysis of symptomatic and asymptomatic norovirus infections in Chinese infants. Virol J 20, 60 (2023). https://doi.org/10.1186/s12985-023-02024-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12985-023-02024-z
Keywords
- Norovirus
- Acute gastroenteritis
- Birth cohort
- Mutation
- Recombination