Skip to main content

Molecular epidemiology analysis of symptomatic and asymptomatic norovirus infections in Chinese infants

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.

Fig. 1
figure 1

Distribution of norovirus genotypes in symptomatic and asymptomatic norovirus infections from a cohort study in China, from 2021–2022

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.

Fig. 2
figure 2

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.

Fig. 3
figure 3

Site variability was calculated at nucleotide level using Shannon entropy for norovirus strains isolated from our study. Diversity nucleotide plots, were shown the difference site of the ORF2 (A) and RdRp (B) sequence. Sequence locus information is referenced to GII.4 Sydney 2012 genome (JX459908)

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).

Table 1 Mapping of immunogenic sites within VP1 region of different GII.4 strains in this study

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.

Fig. 4
figure 4

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).

Fig. 5
figure 5

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

  1. Katayama K. Caliciviridae. Nihon Rinsho. 2003;61(Suppl 3):468–74.

    PubMed  Google Scholar 

  2. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  3. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Jiang X, Wang M, Wang K, Estes MK. Sequence and genomic organization of Norwalk virus. Virology. 1993;195:51–61.

    Article  CAS  PubMed  Google Scholar 

  6. 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.

    Article  CAS  PubMed  Google Scholar 

  7. Tan M, Jiang X. Histo-blood group antigens: a common niche for norovirus and rotavirus. Expert Rev Mol Med. 2014;16: e5.

    Article  PubMed  Google Scholar 

  8. 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.

    Article  CAS  PubMed  Google Scholar 

  9. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 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.

    Article  CAS  PubMed  Google Scholar 

  12. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Robilotti E, Deresinski S, Pinsky BA. Norovirus. Clin Microbiol Rev. 2015;28:134–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  17. The World Health Organization. 2017. https://www.who.int/news-room/fact-sheets/detail/diarrhoeal-disease. Accessed 2 May 2017.

  18. 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.

    Article  PubMed  Google Scholar 

  19. Singh S. Review of epidata entry and analysis freewares. Indian J Community Med. 2009;34:76–7.

    Article  PubMed  PubMed Central  Google Scholar 

  20. 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.

    Article  PubMed  Google Scholar 

  21. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.

    CAS  PubMed  Google Scholar 

  23. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. 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.

    Article  CAS  PubMed  Google Scholar 

  25. Hall BG. Building phylogenetic trees from molecular data with MEGA. Mol Biol Evol. 2013;30:1229–35.

    Article  CAS  PubMed  Google Scholar 

  26. Nemzer LR. Shannon information entropy in the canonical genetic code. J Theor Biol. 2017;415:158–70.

    Article  CAS  PubMed  Google Scholar 

  27. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Larsson A. AliView: a fast and lightweight alignment viewer and editor for large datasets. Bioinformatics. 2014;30:3276–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. 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.

  32. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Janson G, Paiardini A. PyMod 3: a complete suite for structural bioinformatics in PyMOL. Bioinformatics. 2021;37:1471–2.

    Article  CAS  PubMed  Google Scholar 

  34. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Martin DP, Murrell B, Khoosal A, Muhire B. Detecting and analyzing genetic recombination using RDP4. Methods Mol Biol. 2017;1525:433–60.

    Article  CAS  PubMed  Google Scholar 

  36. Martin D, Rybicki E. RDP: detection of recombination amongst aligned sequences. Bioinformatics. 2000;16:562–3.

    Article  CAS  PubMed  Google Scholar 

  37. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. 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.

    Article  PubMed  Google Scholar 

  39. 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.

  40. Bucardo F. Understanding asymptomatic norovirus infections. EClinicalMedicine. 2018;2–3:7–8.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 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.

    Article  CAS  PubMed  Google Scholar 

  42. 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.

    Article  CAS  PubMed  Google Scholar 

  43. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  44. 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.

    Article  CAS  PubMed  Google Scholar 

  45. 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.

    Article  PubMed  Google Scholar 

  46. 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.

    Article  Google Scholar 

  47. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  48. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. 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.

    Article  CAS  PubMed  Google Scholar 

  50. 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.

    Article  PubMed  Google Scholar 

  51. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  52. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  54. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  56. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Parra GI. Emergence of norovirus strains: A tale of two genes. Virus Evol. 2019;5:vez048.

    Article  PubMed  PubMed Central  Google Scholar 

  58. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Eden JS, Tanaka MM, Boni MF, Rawlinson WD, White PA. Recombination within the pandemic norovirus GII.4 lineage. J Virol. 2013;87:6270–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. 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.

    Article  CAS  PubMed  Google Scholar 

  62. Bull RA, Tanaka MM, White PA. Norovirus recombination. J Gen Virol. 2007;88:3347–59.

    Article  CAS  PubMed  Google Scholar 

  63. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Smertina E, Urakova N, Strive T, Frese M. Calicivirus RNA-dependent RNA polymerases: evolution, structure, protein dynamics, and function. Front Microbiol. 2019;10:1280.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Mahar JE, Bok K, Green KY, Kirkwood CD. The importance of intergenic recombination in norovirus GII.3 evolution. J Virol. 2013;87:3687–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

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

Authors

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

Correspondence to Tian-Yi Qiu or Xuan-Yi Wang.

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-023-02024-z

Keywords