Molecular Evolution of the VP1 Gene in Human Norovirus GII.4 Variants in 1974-2015

A DNA vaccine against GII.4 human norovirus VP1 induces blocking antibody production and T cell responses

Human noroviruses (HuNoVs) are highly contagious and a leading cause of epidemics of acute gastroenteritis worldwide. Among the various HuNoV genotypes, GII.4 is the most prevalent cause of outbreaks. However, no vaccines have been approved for HuNoVs to date. DNA vaccines are proposed to serve as an ideal platform against HuNoV since they can be easily produced and customized to express target proteins. In this study, we constructed a CMV/R vector expressing a major structural protein, VP1, of GII.4 HuNoV (CMV/R-GII.4 HuNoV VP1). Transfection of CMV/R-GII.4 HuNoV VP1 into human embryonic kidney 293T (HEK293T) cells resulted in successful expression of VP1 proteins in vitro. Intramuscular or intradermal immunization of mice with the CMV/R-GII.4 HuNoV VP1 construct elicited the production of blocking antibodies and activation of T cell responses against GII.4 HuNoV VP1. Our collective data support the utility of CMV/R-GII.4 HuNoV VP1 as a promising DNA vaccine candidate against GII.4 HuNoV.

Biological Specimen Banking as a Time Capsule to Explore the Temporal Dynamics of Norovirus Epidemiology

Norovirus is recognised as a major cause of epidemic and sporadic acute gastroenteritis (AGE) in all age groups. Information on the genetic diversity of the noroviruses circulating in the 1980s and 1990s, before the development and adoption of dedicated molecular assays, is limited compared with the last decades. Between 1986 and 2020, uninterrupted viral surveillance was conducted in symptomatic children hospitalized with AGE in Palermo, Italy, providing a unique time capsule for exploring the epidemiological and evolutionary dynamics of enteric viruses. A total of 8433 stool samples were tested using real-time RT-PCR. All samples were stored at -20 or -80 °C until processing. In this 35-year long time span, noroviruses of genogroup II (GII) were detected in 15.6% of AGE requiring hospitalization, whilst GI noroviruses were detected in 1.4% of AGE. Overall, the predominant norovirus capsid (Cap) genotype was GII.4 (60.8%), followed by GII.3 (13.3%) and GII.2 (12.4%). Temporal replacement of the GII.4 Cap variants associated with different polymerase (Pol) types were observed over the study period. The chronology of emergence and circulation of the different GII.4 variants were consistent with data available in the literature. Also, for GII.3 and GII.2 NoVs, the circulation of different lineages/strains, differing in either the Cap or Pol genes or in both, was observed. This long-term study revealed the ability of noroviruses to continuously and rapidly modify their genomic makeup and highlights the importance of surveillance activities in vaccine design.

Determination and analysis of whole genome sequence of recombinant GII.6[P7] norovirus in Ningxia, China

While the GII.4 norovirus was the predominant genotype, non-GII.4 genotype was increasingly focused since the non-GII.4 genotype caused regional epidemics. In this study, the detection rate was16.51% (183/1108) in Ningxia from January to December 2020. Among identified genotypes, GII.4[P31] and GII.4[P16] were the dominant genotypes (n = 20 and 18, respectively) while GII.6[P7] was the main type (n = 6) in non-GII.4 strains which was mainly detected in from May to July. The whole genome sequences of the norovirus diarrhea samples identified as GII.6 [P7] with Ct ≤ 30 collected in 2020 were determined. In this study, the complete genome sequences of norovirus strains PL20-044 and QTX20-071 were identified and analyzed phylogenetically. Phylogenetic analysis of the ORF1and ORF2 regions showed that these strains evolved from the GII·P7-GII.6 strains detected in recent years from different country. The results showed that PL20-044 had intra-type recombination with GII·P7-GII.6c and GII·P7-GII.6a, while QTX20-071 had intre-type recombination within GII·P7-GII.6a. The evolutionary rates of the RdRp gene region of the GII·P7 genotype and the VP1 gene region of the GII.6 genotype were 2.91 × 10-3 (95%HPDs2.32-3.51 × 10-3) and 2.61 × 10-3 (95%HPDs2.14-3.11 × 10-3) substitutions/site/year, respectively. Comparative analysis of the amino acid mutation sites in VP1 with the GII·P7-GII.6a strains before 1997, the later detected strains have changed in aa131 and aa354. Moreover, PL20-044 strains showed special mutations at aa316 and aa395. These results help to understand the norovirus genotype circulating in the human population in Ningxia, and discover the evolutionary characteristics of the GII·P7-GII.6 strain.

Molecular and Genetics-Based Systems for Tracing the Evolution and Exploring the Mechanisms of Human Norovirus Infections

Human noroviruses (HuNoV) are major causes of acute gastroenteritis around the world. The high mutation rate and recombination potential of noroviruses are significant challenges in studying the genetic diversity and evolution pattern of novel strains. In this review, we describe recent advances in the development of technologies for not only the detection but also the analysis of complete genome sequences of noroviruses and the future prospects of detection methods for tracing the evolution and genetic diversity of human noroviruses. The mechanisms of HuNoV infection and the development of antiviral drugs have been hampered by failure to develop the infectious virus in a cell model. However, recent studies have demonstrated the potential of reverse genetics for the recovery and generation of infectious viral particles, suggesting the utility of this genetics-based system as an alternative for studying the mechanisms of viral infection, such as cell entry and replication.

Molecular Evolution of RNA-Dependent RNA Polymerase Region in Norovirus Genogroup I

Norovirus is the leading viral agent of gastroenteritis in humans. RNA-dependent RNA polymerase (RdRp) is essential in the replication of norovirus RNA. Here, we present a comprehensive evolutionary analysis of the norovirus GI RdRp gene. Our results show that the norovirus GI RdRp gene can be divided into three groups, and that the most recent common ancestor was 1484. The overall evolutionary rate of GI RdRp is 1.821 × 10-3 substitutions/site/year. Most of the amino acids of the GI RdRp gene were under negative selection, and only a few positively selected sites were recognized. Amino acid substitutions in the GI RdRp gene accumulated slowly over time. GI.P1, GI.P3 and GI.P6 owned the higher evolutionary rates. GI.P11 and GI.P13 had the faster accumulation rate of amino acid substitutions. GI.P2, GI.P3, GI.P4, GI.P6 and GI.P13 presented a strong linear evolution. These results reveal that the norovirus GI RdRp gene evolves conservatively, and that the molecular evolutionary characteristics of each P-genotype are diverse. Sequencing in RdRp and VP1 of norovirus should be advocated in the surveillance system to explore the effect of RdRp on norovirus activity.

Identification and genetic characterization of a minor norovirus genotype, GIX.1[GII.P15], from China

Background

Human noroviruses, single-stranded RNA viruses in the family Caliciviridae, are a leading cause of nonbacterial acute gastroenteritis in people of all ages worldwide. Despite three decades of genomic sequencing and epidemiological norovirus studies, full-length genome analyses of the non-epidemic or minor norovirus genotypes are rare and genomic regions other than ORF2 and 3′-end of ORF1 have been largely understudied, which hampers a better understanding of the evolutionary mechanisms of emergence of new strains. In this study, we detected a rare norovirus genotype, GIX.1[GII.P15], in a vomit sample of a 60 year old woman with acute gastroenteritis using Raji cells and sequenced the complete genome.

Results

Using electron microscopy, a morphology of spherical and lace-like appearance of norovirus virus particles with a diameter of approximately 30 nm were observed. Phylogenetic analysis of VP1 and the RdRp region indicated that the KMN1 strain could be genotyped as GIX.1[GII.P15]. In addition, the VP1 region of KMN1 strain had 94.15% ± 3.54% percent nucleotide identity (PNI) compared to 26 genomic sequences available in GenBank, indicating a higher degree similarity between KMN1 and other GIX.1[GII.P15] strains. Further analysis of the full genome sequence of KMN1 strain showed that a total of 96 nucleotide substitutions (63 in ORF1, 25 in ORF2, and 8 in ORF3) were found across the genome compared with the consensus sequence of GIX.1[GII.P15] genome, and 6 substitutions caused amino acid changes (4 in ORF1, 1 in ORF2, and 1 in ORF3). However, only one nucleotide substitution results in the amino acid change (P302S) in the VP1 protein and the site was located near one of the predicted conformational B epitopes on the dimer structure.

Conclusions

The genomic information of the new GIX.1[GII.P15] strain KMN1, which was identified in Kunming, China could provide helpful insights for the study of the genetic evolution of the virus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01066-6.

Clinical and virological characteristics of viral shedding in children with norovirus gastroenteritis

BACKGROUND

The correlation between the clinical manifestations and fecal viral load of norovirus (NoV) infection remains unknown.

METHODS

We established a SYBR® Green-based real-time quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) method to quantify NoV and then sequenced its genomes from the feces of patients admitted at the Chang Gung Memorial Hospital from 2017 to 2018.

RESULTS

NoV GII.4 Sydney (n = 21, 36.2%) and GII.P16-GII.2 (n = 19, 32.8%), the two predominant genotypes found among 58 isolates, were closely related to the Taiwan variant 2012a cluster in the VP1 region and genotypes of China strain. An increase in viral load could be observed on Day 3 following the onset of NoV infection. The viral load then declined rapidly from days 10-15 but remained high for >1 month in a severe combined immunodeficiency patient. Significantly longer shedding was found in patients with fever (p = 0.03) or infected by the GII.4 Sydney strain (p < 0.01).

CONCLUSION

The qRT-PCR-mediated method proposed in this work could quantify the viral load in patients with NoV infection. Significant viral shedding over a period of 2 weeks in children with acute gastroenteritis and >1 month in an immunodeficient patient was observed. Significantly longer shedding could be correlated with infection by the GII.4 Sydney strain and febrile patients.

Evolutionary Analyses of Emerging GII.2[P16] and GII.4 Sydney [P16] Noroviruses

GII.2[P16] and GII.4 Sydney [P16] are currently the two predominant norovirus genotypes. This study sought to clarify their evolutionary patterns by analyzing the major capsid VP1 and RNA-dependent RNA polymerase (RdRp) genes. Sequence diversities were analyzed at both nucleotide and amino acid levels. Selective pressures were evaluated with the Hyphy package in different models. Phylogenetic trees were constructed by the maximum likelihood method from full VP1 sequences, and evolutionary rates were estimated by the Bayesian Markov Chain Monte Carlo approach. The results showed that (1) several groups of tightly linked mutations between the RdRp and VP1 genes were detected in the GII.2[P16] and GII.4[P16] noroviruses, and most of these mutations were synonymous, which may lead to a better viral fitness to the host; (2) although the pattern of having new GII.4 variants every 2–4 years has been broken, both the pre- and the post-2015 Sydney VP1 had comparable evolutionary rates to previously epidemic GII.4 variants, and half of the major antigenic sites on GII.4 Sydney had residue substitutions and several caused obvious changes in the carbohydrate-binding surface that may potentially alter the property of the virus; and (3) GII.4 Sydney variants during 2018–21 showed geographical specificity in East Asia, South Asia, and North America; the antigenic sites of GII.2 are strictly conserved, but the GII.2 VP1 chronologically evolved into nine different sublineages over time, with sublineage IX being the most prevalent one since 2018. This study suggested that both VP1 and RdRp of the GII.2[P16] and GII.4 Sydney [P16] noroviruses exhibited different evolutionary directions. GII.4[P16] is likely to generate potential novel epidemic variants by accumulating mutations in the P2 domain, similar to previously epidemic GII.4 variants, while GII.2[P16] has conserved predicted antigenicity and may evolve by changing the properties of nonstructural proteins, such as polymerase replicational fidelity and efficiency. This study expands the understanding of the evolutionary dynamics of GII.2[P16] and GII.4[P16] noroviruses and may predict the emergence of new variants.

Norovirus GII.2[P16] strain in Shenzhen, China: a retrospective study

Background

Norovirus (NoV) is the main cause of non-bacterial acute gastroenteritis (AGE) outbreaks worldwide. From September 2015 through August 2018, 203 NoV outbreaks involving 2500 cases were reported to the Shenzhen Center for Disease Control and Prevention.

Methods

Faecal specimens for 203 outbreaks were collected and epidemiological data were obtained through the AGE outbreak surveillance system in Shenzhen. Genotypes were determined by sequencing analysis. To gain a better understanding of the evolutionary characteristics of NoV in Shenzhen, molecular evolution and mutations were evaluated based on time-scale evolutionary phylogeny and amino acid mutations.

Results

A total of nine districts reported NoV outbreaks and the reported NoV outbreaks peaked from November to March. Among the 203 NoV outbreaks, 150 were sequenced successfully. Most of these outbreaks were associated with the NoV GII.2[P16] strain (45.3%, 92/203) and occurred in school settings (91.6%, 186/203). The evolutionary rates of the RdRp region and the VP1 sequence were 2.1 × 10^–3 (95% HPD interval, 1.7 × 10^–3–2.5 × 10^–3) substitutions/site/year and 2.7 × 10^–3 (95% HPD interval, 2.4 × 10^–3–3.1 × 10^–3) substitutions/site/year, respectively. The common ancestors of the GII.2[P16] strain from Shenzhen and GII.4 Sydney 2012[P16] diverged from 2011 to 2012. The common ancestors of the GII.2[P16] strain from Shenzhen and previous GII.2[P16] (2010–2012) diverged from 2003 to 2004. The results of amino acid mutations showed 6 amino acid substitutions (*77E, R750K, P845Q, H1310Y, K1546Q, T1549A) were found only in GII.4 Sydney 2012[P16] and the GII.2[P16] recombinant strain.

Conclusions

This study illustrates the molecular epidemiological patterns in Shenzhen, China, from September 2015 to August 2018 and provides evidence that the epidemic trend of GII.2[P16] recombinant strain had weakened and the non-structural proteins of the recombinant strain might have played a more significant role than VP1.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-06746-9.

Prevalence and Evolution of Noroviruses between 1966 and 2019, Implications for Vaccine Design

Noroviruses (NoVs), a group of single-stranded RNA viruses causing epidemic acute gastroenteritis in humans, are highly diverse, consisting of multiple genogroups with >30 genotypes. Their continual evolutions make NoV vaccine design and development difficult. Here, we report a study of NoV sequences obtained from a population-based diarrhea surveillance in Zhengding County of Hebei Province spanning from 2001 to 2019 and those available in the GenBank database from 1966 to 2019. NoV genotypes and/or variants that may evade immunity were screened and identified based on primary and conformational structures for vaccine design. We selected 366, 301, 139, 74 and 495 complete VP1-coding nucleotide sequences representing the predominant genotypes of GII.4, GII.2, GII.3, GII.6 and GII.17, respectively. A total of 16 distinct GII.4 variants were identified, showing a typical linear evolutionary pattern of variant replacement, while only 1–4 variants of the other genotypes were found to co-circulate over the 40–50-year period without typical variant replacement. The vaccine strain GII.4c is close to variant Sydney_2012 (0.053) in their primary structure, but they are distinct at epitopes A and E in conformations. Our data suggested GII.4 variant Sydney_2012, GII.2 variant A, a GII.3 strain, GII.6 variants B and C and GII.17 variant D are primary candidate strains for NoV vaccine development.

A Longitudinal Study on Enteric Virus Contamination in Bivalves along the Coast of Ibaraki Prefecture, Japan

ABSTRACT

During the 2014 to 2018 seasons, we conducted a longitudinal study involving enteric virus surveillance in bivalves, including natural oysters and clams harvested in Ibaraki Prefecture, Japan. Some norovirus (NoV) contaminations were detected in natural oysters, whereas no enteric virus was found in clams. NoVs detected in oysters were of the genotypes GII.4 and GII.6, both of which are closely related genetically to the NoV strains prevalent in humans. We found low level of enteric virus contamination in bivalves collected along the coast of Ibaraki Prefecture. The possibility of food poisoning caused by these viruses appears low, and few cases of infectious disease have been observed in the surrounding area. The harvest timing was more related to contamination quantity than the harvest area in many enteric viruses. Our results highlight that contamination of bivalves by enteric viruses may depend upon the prevalence of human diarrhea and illness.

HIGHLIGHTS

Evolutionary history and spatiotemporal dynamic of GIII norovirus: From emergence to classification in four genotypes

Noroviruses belong to a genetically diverse group of viruses infecting a wide range of mammalian host species, and those detected in cattle and sheep are classified within genogroup III (GIII). The current classification of norovirus in genogroups and genotypes is based on phylogenetic clustering and average distances within and between these phylogenetic clusters; however, the classification studies have been focused mainly on human norovirus, being GIII norovirus relegated. Due to the increasing number of studies on GIII norovirus, the need of an updated and extensive classification is evident. The aim of this study was to update the classification of norovirus within GIII, to describe the emergence of a circulating recombinant strain, and to reconstruct the evolutionary history of this genogroup. Two P-types (GIII.P1-2) and four genotypes (GIII.1-4) were described. For the genogroup GIII, the evolutionary rate estimated was 2.78E-3 s/s/y (95%HPD, 1.79E-3 s/s/y-3.78E-3 s/s/y), and the tMRCA was estimated around 1500 (95%HPD, 1247-1688). Despite the long history of this genogroup, the genotypes detected at present emerged in the last 100 years. Interestingly, most of the recombinant GIII.2P[1] strains detected worldwide were originated from a single recombination event and this recombinant strain was later dispersed through the world. Finally, our results indicate that a scenario of genotypes replacement through the time is highly probable.

Broad-range and effective detection of human noroviruses by colloidal gold immunochromatographic assay based on the shell domain of the major capsid protein

BACKGROUND

Human noroviruses (HuNoVs) are a major cause of nonbacterial gastroenteritis in all age groups worldwide. HuNoVs can be detected in vitro using molecular assays such as RT-PCR and RT-qPCR. However, these molecular-based techniques require special equipment, unique reagents, experienced personnel, and extended time to obtain results. Besides, the diversity of viral genotypes is high. Therefore, methods that are rapid, broad-range and effective in the detection of HuNoVs are desiderated for screening the feces or vomit of infected people during outbreaks.

RESULTS

In this study, a colloidal-gold-based immunochromatographic assay (ICA) was developed for effective detection of HuNoVs in clinical samples. Monoclonal antibodies (MAbs) against the shell (S) domain in the major capsid protein of HuNoVs were used in the ICA. The limitations of detection for HuNoVs in clinical samples were 1.2 × 106 genomic copies per gram of stool sample (gc/g) and 4.4 × 105 gc/g for genogroup I and II (GI and GII) HuNoVs, respectively. A total of 122 clinical samples were tested for HuNoVs by ICA and compared against RT-qPCR. The relative sensitivity, specificity and agreement of ICA was 84.2% (95% CI: 83.6-84.8%), 100.0% (95% CI: 98.5-100.0%) and 87.7% (95% CI: 85.6-89.8%), respectively. No cross-reaction with other common enteric viruses or bacteria was observed. The ICA detected a broad range of genotypes, including GI.1, GI.3, GI.4, GI.6, GI.14, GII.2, GII.3, GII.4, GII.6, GII.13, and GII.17 HuNoVs.

CONCLUSIONS

This study demonstrates that ICA targeting the S domain of VP1 is a promising candidate for effectively identifying the different genotypes of HuNoVs in clinical samples with high sensitivity and specificity.

Identification of GII.14[P7] norovirus and its genomic mutations from a case of long-term infection in a post-symptomatic individual

Norovirus is a leading cause of acute gastroenteritis worldwide. Norovirus shedding typically lasts one week to one month after the onset of diarrhea in immunocompetent hosts. The occurrence of mutations in the genome during infection has contributed to the evolution of norovirus. It has been suggested that genomic mutations in the P2-domain of capsid protein VP1, the major antigenic site for virus clearance, are involved in the evasion of host immunity and prolonged shedding of norovirus. In our previous study, we found a case of long-term shedding of GII.14 norovirus in a post-symptomatic immunocompetent individual that lasted about three months. In this study, we characterized the genomic sequence of the GII.14 strain to gain insight into the context of long-term shedding. By sequencing a 4.8 kb region of the genome corresponding to half of ORF1 and the entire ORF2 and ORF3, which encode several non-structural proteins and the structural proteins VP1 and VP2, the GII.14 strain was found to be classified as recombinant GII.14[P7]. Six point-mutations occurred during the three-month period of infection in a time-dependent manner in the genomic regions encoding RNA-dependent RNA polymerase, VP1, and VP2. Three of the six mutations were sense mutations, but no amino acid substitution was identified in the P2-domain of VP1. These results suggest that there is a mechanism by which long-term shedding of norovirus occurs in immunocompetent individuals independent of P2-domain mutations.

Preadaptation of pandemic GII.4 noroviruses in unsampled virus reservoirs years before emergence

The control of re-occurring pandemic pathogens requires understanding the origins of new pandemic variants and the factors that drive their global spread. This is especially important for GII.4 norovirus, where vaccines under development offer promise to prevent hundreds of millions of annual gastroenteritis cases. Previous studies have hypothesized that new GII.4 pandemic viruses arise when previously circulating pandemic or pre-pandemic variants undergo substitutions in antigenic regions that enable evasion of host population immunity, as described by conventional models of antigenic drift. In contrast, we show here that the acquisition of new genetic and antigenic characteristics cannot be the proximal driver of new pandemics. Pandemic GII.4 viruses diversify and spread over wide geographical areas over several years prior to simultaneous pandemic emergence of multiple lineages, indicating that the necessary sequence changes must have occurred before diversification, years prior to pandemic emergence. We confirm this result through serological assays of reconstructed ancestral virus capsids, demonstrating that by 2003, the ancestral 2012 pandemic strain had already acquired the antigenic characteristics that allowed it to evade prevailing population immunity against the previous 2009 pandemic variant. These results provide strong evidence that viral genetic changes are necessary but not sufficient for GII.4 pandemic spread. Instead, we suggest that it is changes in host population immunity that enable pandemic spread of an antigenically preadapted GII.4 variant. These results indicate that predicting future GII.4 pandemic variants will require surveillance of currently unsampled reservoir populations. Furthermore, a broadly acting GII.4 vaccine will be critical to prevent future pandemics.

Molecular evolution of the capsid (VP1) region in human norovirus genogroup II genotype 3

Norovirus GII.3 has been suggested to be a prevalent genotype in patients with acute gastroenteritis. However, the genetic properties of the VP1 region encoding the major GII.3 antigen remain unclear. Here, we performed molecular evolutionary analyses of the GII.3 VP1 region detected in various countries. We performed time-scaled phylogenetic analyses, selective pressure analyses, phylogenetic distance analyses, and conformational epitope analyses. The time-scaled phylogenetic tree showed that an ancestor of the GII.3 VP1 region diverged from the common ancestors of the GII.6, GII.11, GII.18, and GII.19 approximately 70 years ago with relatively low divergence. The evolutionary rate of the GII.3 VP1 region was rapid (4.82 × 10⁻³ substitutions/site/year). Furthermore, one positive site and many negative selection sites were observed in the capsid protein. These results suggest that the GII.3 VP1 region rapidly evolved with antigenic variations.

Norovirus Replication in Human Intestinal Epithelial Cells Is Restricted by the Interferon-Induced JAK/STAT Signaling Pathway and RNA Polymerase II-Mediated Transcriptional Responses

Human noroviruses (HuNoV) are a leading cause of viral gastroenteritis worldwide and a significant cause of morbidity and mortality in all age groups. The recent finding that HuNoV can be propagated in B cells and mucosa-derived intestinal epithelial organoids (IEOs) has transformed our ability to dissect the life cycle of noroviruses. Using transcriptome sequencing (RNA-Seq) of HuNoV-infected intestinal epithelial cells (IECs), we have found that replication of HuNoV in IECs results in interferon (IFN)-induced transcriptional responses and that HuNoV replication in IECs is sensitive to IFN. This contrasts with previous studies that suggested that the innate immune response may play no role in the restriction of HuNoV replication in immortalized cells. We demonstrated that inhibition of Janus kinase 1 (JAK1)/JAK2 enhanced HuNoV replication in IECs. Surprisingly, targeted inhibition of cellular RNA polymerase II-mediated transcription was not detrimental to HuNoV replication but instead enhanced replication to a greater degree than blocking of JAK signaling directly. Furthermore, we demonstrated for the first time that IECs generated from genetically modified intestinal organoids, engineered to be deficient in the interferon response, were more permissive to HuNoV infection. Taking the results together, our work revealed that IFN-induced transcriptional responses restrict HuNoV replication in IECs and demonstrated that inhibition of these responses mediated by modifications of the culture conditions can greatly enhance the robustness of the norovirus culture system.IMPORTANCE Noroviruses are a major cause of gastroenteritis worldwide, and yet the challenges associated with their growth in culture have greatly hampered the development of therapeutic approaches and have limited our understanding of the cellular pathways that control infection. Here, we show that human intestinal epithelial cells, which represent the first point of entry of human noroviruses into the host, limit virus replication by induction of innate responses. Furthermore, we show that modulating the ability of intestinal epithelial cells to induce transcriptional responses to HuNoV infection can significantly enhance human norovirus replication in culture. Collectively, our findings provide new insights into the biological pathways that control norovirus infection but also identify mechanisms that enhance the robustness of norovirus culture.

Molecular evolution of GII.P17-GII.17 norovirus associated with sporadic acute gastroenteritis cases during 2013-2018 in Zhoushan Islands, China

In this study, we investigated the molecular characteristics and spatio-temporal dynamics of GII.P17-GII.17 norovirus in Zhoushan Islands during 2013-2018. We collected 1849 samples from sporadic acute gastroenteritis patients between January 2013 and August 2018 in Zhoushan Islands, China. Among the 1849 samples, 134 (7.24%) samples were positive for human norovirus (HuNoV). The complete sequence of GII.17 VP1 gene was amplified from 31 HuNoV-positive samples and sequenced. A phylogenetic tree was constructed based on the full-length sequence of the VP1 gene. Phylogenetic analysis revealed that the GII.17 genotype detected during 2014-2018 belongs to the new GII.17 Kawasaki variant. Divergence analysis revealed that the time of the most recent common ancestor (TMRCA) of GII.17 in Zhoushan Islands was estimated to be between 1997 and 1998. The evolutionary rate of the VP1 gene of the GII.17 genotype norovirus was 1.14 × 10^-3 (95% HPD: 0.62-1.73 × 10^-3) nucleotide substitutions/site/year. The spatio-temporal diffusion analysis of the GII.17 genotype identified Hong Kong as the epicenter for GII.17 dissemination. The VP1 gene sequence of Zhoushan Island isolates correlated with that of Hong Kong and Japan isolates.

Molecular Evolution of the Fusion Protein (F) Gene in Human Respirovirus 3

To elucidate the evolution of human respirovirus 3 (HRV3), we performed detailed genetic analyses of the F gene (full-length) detected from hundreds of HRV3 strains obtained from various geographic regions. First, we performed time-scaled evolutionary analyses using the Bayesian Markov chain Monte Carlo method. Then, we performed analyses of phylodynamics, similarity, phylogenetic distance, selective pressure, and conformational B-cell epitope with the F-protein structural analyses. Time-scaled phylogenetic tree showed that the common ancestor of HRV3 and bovine respirovirus 3 diverged over 300 years ago and subdivided it into three major clusters and four subclusters during the most recent 100 years. The overall evolutionary rate was approximately 10^-3 substitutions/site/year. Indigenous similarity was seen in the present strains, and the mean phylogenetic distance were 0.033. Many negative selection sites were seen in the ectodomain. The conformational epitopes did not correspond to the neutralizing antibody binding sites. These results suggest that the HRV3 F gene is relatively conserved and restricted in this diversity to preserve the protein function, although these strains form many branches on the phylogenetic tree. Furthermore, HRV3 reinfection may be responsible for discordances between the conformational epitopes and the neutralizing antibody binding sites of the F protein. These findings contribute to a better understanding of HRV3 virology.

Molecular Evolution of the Protease Region in Norovirus Genogroup II

Noroviruses are a major cause of viral epidemic gastroenteritis in humans worldwide. The protease (Pro) encoded in open reading frame 1 (ORF1) is an essential enzyme for proteolysis of the viral polyprotein. Although there are some reports regarding the evolutionary analysis of norovirus GII-encoding genes, there are few reports focused on the Pro region. We analyzed the molecular evolution of the Pro region of norovirus GII using bioinformatics approaches. A time-scaled phylogenetic tree of the Pro region constructed using a Bayesian Markov chain Monte Carlo method indicated that the common ancestor of GII diverged from GIV around 1680 CE [95% highest posterior density (HPD), 1607-1749]. The GII Pro region emerged around 1752 CE (95%HPD, 1707-1794), forming three further lineages. The evolutionary rate of GII Pro region was estimated at more than 10-3 substitutions/site/year. The distribution of the phylogenetic distances of each genotype differed, and showed genetic diversity. Mapping of the negative selection and substitution sites of the Pro structure showed that the substitution sites in the Pro protein were mostly produced under neutral selection in positions structurally adjacent to the active sites for proteolysis, whereas negative selection was observed in residues distant from the active sites. The phylodynamics of GII.P4, GII.P7, GII.P16, GII.P21, and GII.P31 indicated that their effective population sizes increased during the period from 2005 to 2016 and the increase in population size was almost consistent with the collection year of these genotypes. These results suggest that the Pro region of the norovirus GII evolved rapidly, but under no positive selection, with a high genetic divergence, similar to that of the RNA-dependent RNA polymerase (RdRp) region and the VP1 region of noroviruses.

Foodborne viral outbreaks associated with frozen produce

Over the past decade, frozen fruits have been a major vehicle of foodborne illnesses mainly attributed to norovirus (NoV) and hepatitis A virus (HAV) infections. Fresh produce may acquire viral contamination by direct contact with contaminated surface, water or hands, and is then frozen without undergoing proper decontamination. Due to their structural integrity, foodborne viruses are able to withstand hostile conditions such as desiccation and freezing, and endure for a long period of time without losing their infectivity. Additionally, these foods are often consumed raw or undercooked, which increases the risk of infection. Herein, we searched published literature and databases of reported outbreaks as well as the databases of news articles for the viral outbreaks associated with the consumption of frozen produce between January 2008 and December 2018; recorded the worldwide distribution of these outbreaks; and analysed the implication of consumption of different types of contaminated frozen food. In addition, we have briefly discussed the factors that contribute to an increased risk of foodborne viral infection following the consumption of frozen produce. Our results revealed that frozen fruits, especially berries and pomegranate arils, contributed to the majority of the outbreaks, and that most outbreaks were reported in industrialised countries.

Evolutionary Analysis of the VP1 and RNA-Dependent RNA Polymerase Regions of Human Norovirus GII.P17-GII.17 in 2013-2017

Human norovirus (HuNoV) GII.P17-GII.17 (Kawasaki2014 variant) reportedly emerged in 2014 and caused gastroenteritis outbreaks worldwide. To clarify the evolution of both VP1 and RNA-dependent RNA polymerase (RdRp) regions of GII.P17-GII.17, we analyzed both global and novel Japanese strains detected during 2013-2017. Time-scaled phylogenetic trees revealed that the ancestral GII.17 VP1 region diverged around 1949, while the ancestral GII.P17 RdRp region diverged around 2010. The evolutionary rates of the VP1 and RdRp regions were estimated at ~2.7 × 10-3 and ~2.3 × 10-3 substitutions/site/year, respectively. The phylogenetic distances of the VP1 region exhibited no overlaps between intra-cluster and inter-cluster peaks in the GII.17 strains, whereas those of the RdRp region exhibited a unimodal distribution in the GII.P17 strains. Conformational epitope positions in the VP1 protein of the GII.P17-GII.17 strains were similar, although some substitutions, insertions and deletions had occurred. Strains belonging to the same cluster also harbored substitutions around the binding sites for the histo-blood group antigens of the VP1 protein. Moreover, some amino acid substitutions were estimated to be near the interface between monomers and the active site of the RdRp protein. These results suggest that the GII.P17-GII.17 virus has produced variants with the potential to alter viral antigenicity, host-binding capability, and replication property over the past 10 years.

Variation of human norovirus GII genotypes detected in Ibaraki, Japan, during 2012-2018

Background

Human norovirus (HuNoV) is the major cause of viral acute gastroenteritis for all age groups in various countries. HuNoV GII in particular accounted for the majority of norovirus outbreaks, among which GII.4 caused repeated outbreaks for a long time. Besides GII.4, other norovirus genotypes, GII.2, GII.6, and GII.17, have also been prevalent in various contexts in recent years, but few detailed epidemiological studies of them have been performed and are poorly understood. We thus conducted an epidemiological analysis of HuNoV GII in Ibaraki Prefecture, Japan, by performing surveillance in the six seasons from September 2012 to August 2018.

Results

HuNoV GI occurred almost sporadically for all genotypes; however, each genotype of GII exhibited its typical epidemiological characteristics. Although the number of outbreaks of GII.4 decreased season by season, it reemerged in 2017/2018 season. The timing of the epidemic peak in terms of number of cases for GII.17 differed from that for the other genotypes. The patients age with GII.2 and GII.6 were younger and outbreak of GII.17 occurred frequently as food poisoning. Namely, the primarily infected outbreak group differed for each genotype of HuNoV GII. Moreover, the viral load of patients differed according to the genotype.

Conclusions

Various HuNoV genotypes including GII.2, GII.4, GII.6, and GII.17 were shown to be associated with various types of outbreak sites (at childcare and educational facilities, involving cases of food poisoning, and at elderly nursing homes) in this study. These genotypes emerged in recent years, and their prevalence patterns differed from each other. Moreover, differences in outbreak sites and viral load of patients among the genotypes were identified.

The coordinating role of the human norovirus minor capsid protein VP2 is essential to functional change and nuclear localization of the major capsid protein VP1

Global outbreaks of norovirus (NOV) gastroenteritis are associated with the most prevalent genotype, GII.4. Mutations in the protruding domain 2 (P2 domain) of the norovirus major capsid protein (VP1) result in the emergence of various NOV variants, however, it is unclear whether the minor capsid protein (VP2) also affects the generation of VP1 variants. In this study, using a human 293T expression system, we investigated the interactions of VP1 and VP2 of three GII.4 strains, focusing on the changes in expression and cellular localization. We found that co-transfection with VP1 and VP2 leads to a significant increase in expression of both proteins compared to that in cells transfected with VP1 or VP2 alone. In contrast to VP1 expressed in the absence of VP2, which was dispersed throughout the cytosol, VP2 expressed in the absence of VP1 was found to be located in the nucleus. This could be attributed to a predicted specific nuclear localization signal found in this gene. When both proteins were expressed, VP1 was found together with VP2 in the nucleus. These results thus suggest that the VP2 of GII.4 NOVs affects the function and cellular location of VP1 and that, with the cooperation of VP2, VP1 could play a critical role in affecting cell functions by impairing the downstream transcriptional signaling and chromatin remodeling in the cell nuclei.

Molecular Evolutionary Analyses of the RNA-Dependent RNA Polymerase Region in Norovirus Genogroup II

Noroviruses are the leading cause of viral gastroenteritis in humans across the world. RNA-dependent RNA polymerase (RdRp) plays a critical role in the replication of the viral genome. Although there have been some reports on a limited number of genotypes with respect to the norovirus evolution of the RdRp region, no comprehensive molecular evolution examination of the norovirus GII genotype has yet been undertaken. Therefore, we conducted an evolutionary analysis of the 25 genotypes of the norovirus GII RdRp region (full-length), collected globally using different bioinformatics technologies. The time-scaled phylogenetic tree, generated using the Bayesian Markov Chain Monte Carlo (MCMC) method, indicated that the common ancestor of GII diverged from GIV around 1443 CE [95% highest posterior density (HPD), 1336–1542]. The GII RdRp region emerged around 1731 CE (95% HPD, 1703–1757), forming three lineages. The evolutionary rate of the RdRp region of the norovirus GII strains was estimated at over 10⁻³ substitutions/site/year. The evolutionary rates were significantly distinct in each genotype. The composition of the phylogenetic distances differed among the strains for each genotype. Furthermore, we mapped the negative selection sites on the RdRp protein and many of these were predicted in the GII.P4 RdRp proteins. The phylodynamics of GII.P4, GII.P12, GII.P16, and GII.Pe showed that their effective population sizes increased during the period from 2003 to 2014. Our results cumulatively suggest that the RdRp region of the norovirus GII rapidly and uniquely evolved with a high divergence similar to that of the norovirus VP1 gene.

Genetic characterization of norovirus GII.4 variants circulating in Canada using a metagenomic technique

BACKGROUND

Human norovirus is the leading cause of viral gastroenteritis globally, and the GII.4 has been the most predominant genotype for decades. This genotype has numerous variants that have caused repeated epidemics worldwide. However, the molecular evolutionary signatures among the GII.4 variants have not been elucidated throughout the viral genome.

METHOD

A metagenomic, next-generation sequencing method, based on Illumina RNA-Seq, was applied to determine norovirus sequences from clinical samples.

RESULTS

Herein, the obtained deep-sequencing data was employed to analyze full-genomic sequences from GII.4 variants prevailing in Canada from 2012 to 2016. Phylogenetic analysis demonstrated that the majority of these sequences belong to New Orleans 2009 and Sydney 2012 strains, and a recombinant sequence was also identified. Genome-wide similarity analyses implied that while the capsid gene is highly diverse among the isolates, the viral protease and polymerase genes remain relatively conserved. Numerous amino acid substitutions were observed at each putative antigenic epitope of the VP1 protein, whereas few amino acid changes were identified in the polymerase protein. Co-infection with other enteric RNA viruses was investigated and the astrovirus genome was identified in one of the samples.

CONCLUSIONS

Overall this study demonstrated the application of whole genome sequencing as an important tool in molecular characterization of noroviruses.

Molecular epidemiology and spatiotemporal dynamics of norovirus associated with sporadic acute gastroenteritis during 2013-2017, Zhoushan Islands, China

A total of 1 590 fecal swabs and stool samples from sporadic acute gastroenteritis patients of all ages were collected from January 2013 to March 2018 in the Zhoushan Islands, China, with 99 (6.23%) samples subsequently identified as human norovirus (HuNoV) positive. Phylogenetic analysis of partial RdRp and VP1 gene regions identified 10 genotypes of the GII genogroup and 3 genotypes of the GI genogroup. The predominant genotype was GII.P17-GII.17 (42.86%, 33/77), followed by GII.Pe-GII.4_Sydney 2012 (24.68%, 19/77) and GII.P16-GII.2 (12.96%, 10/77). However, the prevailing genotype in the Zhoushan Islands has shifted on three separate occasions. The GII.Pe-GII.4_Sydney_2012 strain was dominant in 2013-2014, the GII.P17-17 strain was dominant in 2015-2016, and the GII.P16-GII.2 strain was dominant in 2017. Divergence analysis showed that the re-emerging GII.P16-GII.2 strains clustered with the Japanese 2010-2012 GII.P16-GII.2 strains, and the time of the most recent common ancestor was estimated to have occurred in 2012 to 2013. The evolutionary rates of the RdRp gene region of the GII.P16 genotype and the VP1 gene region of the GII.2 genotype were 2.64 × 10(-3) (95% HPD interval, 2.17-3.08 × 10(-3)) and 3.36 × 10(-3) (95% HPD interval, 2.66-4.04 × 10(-3)) substitutions/site/year, respectively. The migration pattern of the HuNoV GII.2 genotype in China demonstrated that the re-emerging GII.P16-GII.2 strains were first introduced into Hong Kong from Japan, and then spread from Hong Kong to other coastal areas. Our results also showed that the GII.P16-GII.2 strains in the Zhoushan Islands were likely introduced from Jiangsu Province, China, in 2016.