Molecular evolution of the capsid gene in human norovirus genogroup II

Integrated in-silico design and in vivo validation of multi-epitope vaccines for norovirus

BACKGROUND

Norovirus (NoVs) is a foodborne pathogen that causes acute gastroenteritis. The diversity of its principal antigenic protein poses a significant challenge to vaccine development and the prevention of large-scale outbreaks globally. Currently, no licensed vaccines against norovirus have been approved.

METHODS

We developed a novel pipeline that integrates multiple bioinformatics tools to design broad-spectrum vaccines against NoVs. Specifically, broad-spectrum T-cell epitope vaccines were designed based on consensus sequences and optimized epitope screening, while broad-spectrum B-cell spatial epitope vaccines were constructed using high-throughput antigenicity calculations and epitope mapping.

RESULTS

This pipeline underwent rigorous validation at three levels: firstly, In silico validation: Analysis of properties and structures demonstrated the appropriateness of amino acid composition and the structural integrity of the vaccine sequences. Secondly, theoretical assessment: Evaluation of human leukocyte antigen (HLA) subtype and antigenicity coverage indicated a broad theoretical protective spectrum for the designed vaccine immunogens. Furthermore, in silico simulation confirmed their ability to elicit an immune response. Finally, animal-level validation: Experiments in mice showed that both vaccine immunogens stimulated high levels of IgG and IgA. Notably, Vac-B induced a strong IgG response against GII.2 and a robust IgA response against GII.17, comparable to the immune response elicited by the wild-type NoV non-replicating virus-like particle (VLP) protein group.

CONCLUSIONS

Both in silico and in vivo experimental findings suggest that the proposed pipeline and vaccine immunogens could serve as valuable theoretical guidance for the development of multi-epitope vaccines against NoVs.

Molecular Evolutionary Analyses of the RNA-Dependent RNA Polymerase (RdRp) Region and VP1 Gene in Sapovirus GI.1 and GI.2

Human sapovirus (HuSaV) is a significant cause of gastroenteritis. This study aims to analyze the evolutionary dynamics of the RNA-dependent RNA polymerase (RdRp) and capsid (VP1) genes of the HuSaV GI.1 and GI.2 genotypes between 1976 and 2020. Using bioinformatics tools such as the Bayesian phylogenetics software BEAST 2 package (v.2.7.6), we constructed time-scale evolutionary trees based on the gene sequences. Most of the recent common ancestors (MRCAs) of the RdRp region and VP1 gene in the present HuSaV GI.1 diverged around 1930 and 1933, respectively. The trees of the HuSaV GI.1 RdRp region and VP1 gene were divided into two clusters. Further, the MRCAs of the RdRp region and VP1 gene in HuSaV GI.2 diverged in 1960 and 1943, respectively. The evolutionary rates were higher for VP1 gene in HuSaV GI.1 than that in HuSaV GI.2, furthermore, were higher in GI.1 Cluster B than GI.1 Cluster A. In addition, a steep increase was observed in the time-scaled genome population size of the HuSaV GI.1 Cluster B. These results indicate that the HuSaV GI.1 Cluster B may be evolving more actively than other genotypes. The conformational B-cell epitopes were predicted with a higher probability in RdRp for GI.1 and in VP1 for GI.2, respectively. These results suggest that the RdRp region and VP1 gene in HuSaV GI.1 and GI.2 evolved uniquely. These findings suggest unique evolutionary patterns in the RdRp region and VP1 gene of HuSaV GI.1 and GI.2, emphasizing the need for a 'One Health' approach to better understand and combat this pathogen.

Establishment of a Nucleic Acid Detection Method for Norovirus GII.2 Genotype Based on RT-RPA and CRISPR/Cas12a-LFS

To establish a rapid detection method for norovirus GII.2 genotype, this study employed reverse transcription recombinase polymerase amplification (RT-RPA) combined with CRISPR/Cas12a and lateral flow strip (RT-RPA-Cas12a-LFS). Here, the genome of norovirus GII.2 genotype was compared to identify highly conserved sequences, facilitating the design of RT-RPA primers and crRNA specific to the conserved regions of norovirus GII.2. Subsequently, the reaction parameters of RT-RPA were optimized and evaluated using agar-gel electrophoresis and LFS. The results indicate that the conserved sequences of norovirus GII.2 were successfully amplified through RT-RPA at 37°C for 25 minutes. Additionally, CRISPR/Cas12a-mediated cleavage detection was achieved through LFS at 37°C within 10 minutes using the amplification products as templates. Including the isothermal amplification reaction time, the total time is 35 minutes. The established RT-RPA-Cas12a-LFS method demonstrated specific detection of norovirus GII.2, yielding negative results for other viral genomes, and exhibited an excellent detection limit of 10 copies/μl. The RT-RPA-Cas12a-LFS method was further compared with qRT-PCR by analyzing 60 food-contaminated samples. The positive conformity rate was 100%, the negative conformity rate was 95.45%, and the overall conformity rate reached 98.33%. This detection method for norovirus GII.2 genotype is cost-effective, highly sensitive, specific, and easy to operate, offering a promising technical solution for field-based detection of the norovirus GII.2 genotype.

Simultaneous detection of human norovirus GI, GII and SARS-CoV-2 by a quantitative one-step triplex RT-qPCR

Background

There are many similarities in the clinical manifestations of human norovirus and SARS-CoV-2 infections, and nucleic acid detection is the gold standard for diagnosing both diseases. In order to expedite the identification of norovirus and SARS-CoV-2, a quantitative one-step triplex reverse transcription PCR (RT-qPCR) method was designed in this paper.

Methods

A one-step triplex RT-qPCR assay was developed for simultaneous detection and differentiation of human norovirus GI (NoV-GI), GII (NoV-GII) and SARS-CoV-2 from fecal specimens.

Results

The triplex RT-qPCR assay had high detection reproducibility (CV < 1%) and sensitivity. The lower limits of detection (LLOD95) of the triplex RT-qPCR assay for each target site were 128.5-172.8 copies/mL, and LLOD95 of the singleplex RT-qPCR assay were 110.3-142.0 copies/mL. Meanwhile, among the detection of clinical oropharyngeal swabs and fecal specimens, the results of the singleplex and triplex RT-qPCR assay showed high agreement.

Conclusion

The triplex RT-qPCR assay for simultaneous detection of NoV-GI, NoV-GII and SARS-CoV-2 from fecal specimens has high clinical application value.

Molecular Evolutionary Analyses of the Spike Protein Gene and Spike Protein in the SARS-CoV-2 Omicron Subvariants

To better understand the evolution of the SARS-CoV-2 Omicron subvariants, we performed molecular evolutionary analyses of the spike (S) protein gene/S protein using advanced bioinformatics technologies. First, time-scaled phylogenetic analysis estimated that a common ancestor of the Wuhan, Alpha, Beta, Delta variants, and Omicron variants/subvariants diverged in May 2020. After that, a common ancestor of the Omicron variant generated various Omicron subvariants over one year. Furthermore, a chimeric virus between the BM.1.1.1 and BJ.1 subvariants, known as XBB, diverged in July 2021, leading to the emergence of the prevalent subvariants XBB.1.5 and XBB.1.16. Next, similarity plot (SimPlot) data estimated that the recombination point (breakpoint) corresponded to nucleotide position 1373. As a result, XBB.1.5 subvariants had the 5' nucleotide side from the breakpoint as a strain with a BJ.1 sequence and the 3' nucleotide side as a strain with a BM.1.1.1 sequence. Genome network data showed that Omicron subvariants were genetically linked with the common ancestors of the Wuhan and Delta variants, resulting in many amino acid mutations. Selective pressure analysis estimated that the prevalent subvariants, XBB.1.5 and XBB.1.16, had specific amino acid mutations, such as V445P, G446S, N460K, and F486P, located in the RBD when compared with the BA.4 and BA.5 subvariants. Moreover, some representative immunogenicity-associated amino acid mutations, including L452R, F486V, R493Q, and V490S, were also found in these subvariants. These substitutions were involved in the conformational epitopes, implying that these mutations affect immunogenicity and vaccine evasion. Furthermore, these mutations were identified as positive selection sites. These results suggest that the S gene/S protein Omicron subvariants rapidly evolved, and mutations observed in the conformational epitopes may reduce the effectiveness of the current vaccine, including bivalent vaccines such as mRNA vaccines containing the BA.4/BA.5 subvariants.

Predicting Dominant Genotypes in Norovirus Seasons in Japan

Human noroviruses are an etiological agent of acute gastroenteritis. Since multiple genotypes co-circulate every season changing their proportions, it may be desirable to develop multivalent vaccines by formulating genotype composition of seed strains to match that of dominant strains. Here, performances of the models for predicting dominant genotypes, defined as the two most prevalent genotypes, were evaluated using observed genotype frequencies in Japan and genomic sequences for GI and GII strains. In the null model, genotype proportions in the target season were predicted to be the same as those in the immediately preceding season. In the fitness model, genotype proportions were predicted taking into account the acquisition of novel P-types through recombination and genotype-specific proliferation efficiency, as well as herd immunity to VP1 assuming the duration (d) of 0-10 years. The null model performed better in GII than in GI, apparently because dominant genotypes were more stable in the former than in the latter. Performance of the fitness model was similar to that of the null model irrespective of the assumed value of d. However, performance was improved when dominant genotypes were predicted as the union of those predicted with d = 0-10, suggesting that d may vary among individuals.

Molecular Evolutionary Analyses of the RNA-Dependent RNA Polymerase (RdRp) Region and VP1 Gene in Human Norovirus Genotypes GII.P6-GII.6 and GII.P7-GII.6

To understand the evolution of GII.P6-GII.6 and GII.P7-GII.6 strains, the prevalent human norovirus genotypes, we analysed both the RdRp region and VP1 gene in globally collected strains using authentic bioinformatics technologies. A common ancestor of the P6- and P7-type RdRp region emerged approximately 50 years ago and a common ancestor of the P6- and P7-type VP1 gene emerged approximately 110 years ago. Subsequently, the RdRp region and VP1 gene evolved. Moreover, the evolutionary rates were significantly faster for the P6-type RdRp region and VP1 gene than for the P7-type RdRp region and VP1 genes. Large genetic divergence was observed in the P7-type RdRp region and VP1 gene compared with the P6-type RdRp region and VP1 gene. The phylodynamics of the RdRp region and VP1 gene fluctuated after the year 2000. Positive selection sites in VP1 proteins were located in the antigenicity-related protruding 2 domain, and these sites overlapped with conformational epitopes. These results suggest that the GII.6 VP1 gene and VP1 proteins evolved uniquely due to recombination between the P6- and P7-type RdRp regions in the HuNoV GII.P6-GII.6 and GII.P7-GII.6 virus strains.

Epidemiology and genetic diversity of norovirus GII genogroups among children in Hubei, China, 2017-2019

Norovirus (NoV) is an important cause of viral acute gastroenteritis (AGE). To gain insights into the epidemiological characteristics and genetic diversity of NoV among children in Hubei, 1216 stool samples from children (≤ 5 years) obtained under AGE surveillance from January 2017 to December 2019 were analyzed. The results showed that NoV was responsible for 14.64% of AGE cases, with the highest detection rate in children aged 7-12 months (19.76%). Statistically significant differences were found between male and female infection rates (χ2 ​= ​8.108, P ​= ​0.004). Genetic analysis of RdRp and VP1 sequences showed that NoV GII genotypes were GII.4 Sydney [P31] (34.35%), GII.3 [P12] (25.95%), GII.2 [P16] (22.90%), GII.4 Sydney [P16] (12.98%), GII.17 [P17] (2.29%), GII.6 [P7] and GII.3 [P16] (each at 0.76%). GII.17 [P17] variants were divided into the Kawasaki323-like lineage and the Kawasaki308-like lineage. A unique recombination event was detected between strains of GII.4 Sydney 2012 and GII.4 Sydney 2016. Significantly, all GII.P16 sequences associated with GII.4/GII.2 obtained in Hubei were correlated with novel GII.2 [P16] variants that re-emerged in Germany in 2016. Antigenic site analysis of complete VP1 sequences from all GII.4 variants from Hubei identified notable variable residues of antibody epitopes. Genotyping under continuous AGE surveillance and observation of the antigenic sites of VP1 are important monitoring strategies for emerging NoV strains.

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.

High prevalence of norovirus GII.4 Sydney among children with acute gastroenteritis in Bangladesh, 2018-2021

BACKGROUND

Active molecular surveillance and rapid diagnosis method to track an outbreak of norovirus in Bangladesh is lacking. This study aims to determine the genotypic diversity, molecular epidemiology and evaluate a rapid diagnosis method.

METHODS

A total of 404 fecal specimens were collected from children aged below 60 months from January 2018 to December 2021. All samples were analyzed by reverse transcriptase polymerase chain reaction molecular sequencing of partial VP1 nucleotide. Immunochromatography kit (IC, IP Rota/Noro) was evaluated against reference test method.

RESULTS

We found norovirus in 6.7 % (27 of 404) fecal specimens. A wide diversity of norovirus genotype including GII.3, GII.4, GII.5, GII.6, GII.7, and GII.9 were detected. Norovirus strain GII.4 Sydney-2012 was the most predominant (74 %, 20 of 27) followed by GII.7 (7.4 %), GII.9 (7.4 %), GII.3 (3.7 %), GII.5 (3.7 %) and GII.6 (3.7 %), respectively. Co-infection of rotavirus and norovirus (19 [4.7 %] of 404) was the most prevalent. We found higher odds of prolonged health impact [OR 1.93 (95 % CI 0.87-3.12) (p = .001)] among patients with co-infection. The incidence of norovirus was significant among the children below 24 months (p = 0.001). Significant relation of temperature with the cases of norovirus was detected (p = 0.001). The IC kit provided high specificity (99.3 %) and sensitivity (100 %) for the detection of norovirus.

CONCLUSIONS

This study will provide an integrated insight on the genotypic diversity and rapid identification method of norovirus in Bangladesh.

Patterns and Temporal Dynamics of Natural Recombination in Noroviruses

Noroviruses infect a wide range of mammals and are the major cause of gastroenteritis in humans. Recombination at the junction of ORF1 encoding nonstructural proteins and ORF2 encoding major capsid protein VP1 is a well-known feature of noroviruses. Using all available complete norovirus sequences, we systematically analyzed patterns of natural recombination in the genus Norovirus both throughout the genome and across the genogroups. Recombination events between nonstructural (ORF1) and structural genomic regions (ORF2 and ORF3) were found in all analyzed genogroups of noroviruses, although recombination was most prominent between members of GII, the most common genogroup that infects humans. The half-life times of recombinant forms (clades without evidence of recombination) of human GI and GII noroviruses were 10.4 and 8.4-11.3 years, respectively. There was evidence of many recent recombination events, and most noroviruses that differed by more than 18% of nucleotide sequence were recombinant relative to each other. However, there were no distinct recombination events between viruses that differed by over 42% in ORF2/3, consistent with the absence of systematic recombination between different genogroups. The few inter-genogroup recombination events most likely occurred between ancient viruses before they diverged into contemporary genogroups. The recombination events within ORF1 or between ORF2/3 were generally rare. Thus, noroviruses routinely exchange full structural and nonstructural blocks of the genome, providing a modular evolution.

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.

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.

Molecular epidemiology and genetic diversity of norovirus infection in children with acute gastroenteritis in Bangladesh, 2014-2019

Acute gastroenteritis (AGE) is one of the most common diseases in children, and it continues to be a significant cause of morbidity and mortality worldwide. Norovirus is one of the major enteropathogens associated with both sporadic diarrhea and outbreaks of gastroenteritis. This study aims to investigate genotype diversity and molecular epidemiology of norovirus in Bangladesh. A total of 466 fecal specimens were collected from January 2014 to January 2019 from children below 5 years old with AGE in Bangladesh. All samples were analyzed by reverse transcriptase polymerase chain reaction to detect norovirus, and sequence analysis was conducted if found positive. Norovirus was detected in 5.1% (24 of 466) fecal specimens. Norovirus genotype GII.7 was predominant (62.5%, 15 of 24), followed by GII.3 (37.5%, 9 of 24). Coinfection between rotavirus and norovirus was found in 7 of 24 positive cases. Diarrhea (93.7%) and dehydration (89%) were the most common symptoms in children with AGE. About 80% of the positive cases were detected in children aged under 24 months. One seasonal peak (87.5% infection) was detected in the winter. This study suggests that norovirus continues to be one of the major etiologies of children AGE in Bangladesh. This study will provide a guideline to assess the burden of norovirus infection in Bangladesh, which will assist to combat against AGE.

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.

Evolution of the GII.3[P12] Norovirus from 2010 to 2019 in Jiangsu, China

Objectives

Norovirus genotype GII.3[P12] strains have been an important pathogen for sporadic gastroenteritis infection. In previous studies of GII.3[P12], the number of specimens and time span are relatively small, which is difficult to truly reflect the infection and evolution of this type of norovirus. Here we report a molecular epidemiological study of the NoVs prevalent in Jiangsu between 2010 and 2019 to investigate the evolution of the GII.3[P12] strains in China.

Methods

In this study 60 GII.3[P12] norovirus strains were sequenced and analyzed for evolution, recombination, and selection pressure using bioanalysis software.

Results

The GII.3[P12] strains were continuously detected during the study period, which showed a high constituent ratio in males, in winter and among children aged 0–11 months, respectively. A time-scaled evolutionary tree showed that both GII.P12 RdRp and GII.3 VP1 sequences were grouped into three major clusters (Cluster I–III). Most GII.3[P12] strains were mainly located in sub-cluster (SC) II of Cluster III. A SimPlot analysis identified GII.3[P12] strain to be as an ORF1-intragenic recombinant of GII.4[P12] and GII.3[P21]. The RdRp genes of the GII.3[P12] showed a higher mean substitution rate than those of all GII.P12, while the VP1 genes of the GII.3[P12] showed a lower mean substitution rate than those of all GII.3. Alignment of the GII.3 capsid sequences revealed that three HBGA binding sites of all known GII.3 strains remained conserved, while several amino acid mutations in the predicted antibody binding sites were detected. The mutation at 385 was within predicted antibody binding regions, close to host attachment factor binding sites. Positive and negative selection sites were estimated. Two common positively selected sites (sites 385 and 406) were located on the surface of the protruding domain. Moreover, an amino acid substitution (aa204) was estimated to be near the active site of the RdRp protein.

Conclusions

We conducted a comprehensive analysis on the epidemic and evolution of GII.3[P12] noroviruses and the results suggested that evolution was possibly driven by intergenic recombination and mutations in some key amino acid sites.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13099-021-00430-8.

Molecular Evolution of Human Norovirus GII.2 Clusters

Background

The human norovirus GII.2 outbreak during the 2016–2017 winter season was of unprecedented scale and geographic distribution.

Methods

We analyzed 519 complete VP1 gene sequences of the human norovirus GII.2 genotype sampled during the 2016–2017 winter season, as well as prior (dating back to 1976) from 7 countries. Phylodynamic analyses of these sequences were performed using maximum likelihood and Bayesian statistical frameworks in order to estimate viral evolutionary and population dynamics associated with the outbreak.

Results

Our results revealed an increase in the genetic diversity of human norovirus GII.2 during the recent Asian outbreak and diversification was characterized by at least eight distinct clusters. Bayesian estimation of viral population dynamics revealed a highly fluctuating effective population size, increasing in frequency during the past 15 years.

Conclusion

Despite an increasing viral diversity, we found no evidence of an elevated evolutionary rate or significant selection pressure in human norovirus GII.2, indicating viral evolutionary adaptation was not responsible for the volatility of or spread of the virus during this time.

Linear B-Cell Epitopes in Human Norovirus GII.4 Capsid Protein Elicit Blockade Antibodies

Human norovirus (HuNoV) is the leading cause of nonbacterial gastroenteritis worldwide with the GII.4 genotype accounting for over 80% of infections. The major capsid protein of GII.4 variants is evolving rapidly, resulting in new epidemic variants with altered antigenic potentials that must be considered for the development of an effective vaccine. In this study, we identify and characterize linear blockade B-cell epitopes in HuNoV GII.4. Five unique linear B-cell epitopes, namely P2A, P2B, P2C, P2D, and P2E, were predicted on the surface-exposed regions of the capsid protein. Evolving of the surface-exposed epitopes over time was found to correlate with the emergence of new GII.4 outbreak variants. Molecular dynamic simulation (MD) analysis and molecular docking revealed that amino acid substitutions in the putative epitopes P2B, P2C, and P2D could be associated with immune escape and the appearance of new GII.4 variants by affecting solvent accessibility and flexibility of the antigenic sites and histo-blood group antigens (HBAG) binding. Testing the synthetic peptides in wild-type mice, epitopes P2B (336-355), P2C (367-384), and P2D (390-400) were recognized as GII.4-specific linear blockade epitopes with the blocking rate of 68, 55 and 28%, respectively. Blocking rate was found to increase to 80% using the pooled serum of epitopes P2B and P2C. These data provide a strategy for expanding the broad blockade potential of vaccines for prevention of NoV infection.

Molecular epidemiology and genetic diversity of norovirus infection in children hospitalized with acute gastroenteritis in East Java, Indonesia in 2015-2019

Noroviruses are recognized as a leading cause of outbreaks and sporadic cases of acute gastroenteritis (AGE) among individuals of all ages worldwide, especially in children <5 years old. We investigated the epidemiology of noroviruses among hospitalized children at two hospitals in East Java, Indonesia. Stool samples were collected from 966 children with AGE during September 2015-July 2019. All samples were analyzed by reverse transcription-polymerase chain reaction (RT-PCR) for the amplification of both the RNA-dependent RNA polymerase (RdRp) and the capsid genes of noroviruses. The genotypes were determined by phylogenetic analyses. In 2015-2019, noroviruses were detected in 12.3% (119/966) of the samples. Children <2 years old showed a significantly higher prevalence than those ≥2 years old (P = 0.01). NoV infections were observed throughout the year, with the highest prevalence in December. Based on our genetic analyses of RdRp, GII.[P31] (43.7%, 31/71) was the most prevalent RdRp genotype, followed by GII.[P16] (36.6%, 26/71). GII.[P31] was a dominant genotype in 2016 and 2018, whereas GII.[P16] was a dominant genotype in 2015 and 2017. Among the capsid genotypes, the most predominant norovirus genotype from 2015 to 2018 was GII.4 Sydney_2012 (33.6%, 40/119). The most prevalent genotype in each year was GII.13 in 2015, GII.4 Sydney_2012 in 2016 and 2018, and GII.3 in 2017. Based on the genetic analyses of RdRp and capsid sequences, the strains were clustered into 13 RdRp/capsid genotypes; 12 of them were discordant, e.g., GII.4 Sydney[P31], GII.3[P16], and GII.13[P16]. The predominant genotype in each year was GII.13[P16] in 2015, GII.4 Sydney[P31] in 2016, GII.3[P16] in 2017, and GII.4 Sydney[P31] in 2018. Our results demonstrate high detection rates and genetic diversity of norovirus GII genotypes in pediatric AGE samples from Indonesia. These findings strengthen the importance of the continuous molecular surveillance of emerging norovirus strains.

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.

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.

Increasing Recombinant Strains Emerged in Norovirus Outbreaks in Jiangsu, China: 2015-2018

From January 2015 to December 2018, 213 norovirus outbreaks with 3,951 patients were reported in Jiangsu, China. Based on viral RdRp and VP1 genes, eight genotypes, GII.2[P16] (144, 67.6%), GII.3[P12] (21, 9.9%), GII.6[P7] (5, 2.3%), GII.14[P7] (4, 1.9%), GII.4 Sydney[P31] (3, 1.4%), GII.1[P33] (1, 0.5%), GII.2[P2] (3, 1.4%), and GII.17[P17] (16, 7.5%) were identified throughout the study period. These genotypes were further regrouped as GII.R (Recombinant) and GII.Non-R (Non-recombinant) strains. In this report we showed that GII.R strains were responsible for at least 178 (83.6%) of 213 norovirus-positive outbreaks with a peak in 2017 and 2018. Most norovirus outbreaks occurred in primary schools and 94 of 109 (86.2%) outbreaks in primary schools were caused by GII.R, while GII.Non-R and GII.NT (not typed) strains accounted for 6 (5.5%) and 9 (8.3%) norovirus outbreaks, respectively. The SimPlot analysis showed recombination breakpoints near the ORF1/2 junction for all six recombinant strains. The recombination breakpoints were detected at positions varying from nucleotides 5009 to 5111, localized in the ORF1 region for four strains (GII.2[P16], GII.3[P12], GII.6[P7], and GII.14[P7]) and in the ORF2 region for the other (GII.4 Sydney[P31] and GII.1[P33]). We identified four clusters, Cluster I through IV, in the GII.P7 RdRp gene by phylogenetic analysis and the GII.14[P7] variants reported here belonged to Cluster IV in the RdRp tree. The HBGA binding site of all known GII.14 strains remained conserved with several point mutations found in the predicted conformational epitopes. In conclusion, gastroenteritis outbreaks caused by noroviruses increased rapidly in the last years and these viruses were classified into eight genotypes. Emerging recombinant noroviral strains have become a major concern and challenge to public health.

Molecular Evolution of Classic Human Astrovirus, as Revealed by the Analysis of the Capsid Protein Gene

Classic human astroviruses (HAstV) are major global viral agents for gastroenteritis, but the molecular characteristics of classic HAstVs are not well understood. Here, we presented the molecular evolution of all classic HAstV serotypes by the analysis of the capsid protein sequences. Our results show that classic HAstVs can be divided into four groups with the most recent common ancestor (TMRCA) of 749. The overall evolutionary rate of classic HAstVs on the capsid gene was 4.509 × 10⁻⁴ substitutions/site/year, and most of the serotypes present a clock-like evolution with an amino acid accumulation of mutations over time. The mean effective population size of classic HAstVs is in a downward trend, and some positive and more than 500 negative selection sites were determined. Taken together, these results reveal that classic HAstVs evolve at the intra-serotype level with high genetic heterogeneity and are driven by strong purifying selection. Long-term surveillance of classic HAstVs are needed to enrich the genomic data for further analysis.

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.

Quantitative risk assessment of norovirus and adenovirus for the use of reclaimed water to irrigate lettuce in Catalonia

Wastewater is an important resource in water-scarce regions of the world, and its use in agriculture requires the guarantee of acceptable public health risks. The use of fecal indicator bacteria to evaluate safety does not represent viruses, the main potential health hazards. Viral pathogens could complement the use of fecal indicator bacteria in the evaluation of water quality. In this study, we characterized the concentration and removal of human adenovirus (HAdV) and norovirus genogroup II (NoV GII), highly abundant and important viral pathogens found in wastewater, in two wastewater treatment plants (WWTPs) that use different tertiary treatments (constructed wetland vs conventional UV, chlorination and Actiflo® treatments) for a year in Catalonia. The main objective of this study was to develop a Quantitative Microbial Risk Assessment for viral gastroenteritis caused by norovirus GII and adenovirus, associated with the ingestion of lettuce irrigated with tertiary effluents from these WWTPs. The results show that the disease burden of NoV GII and HAdV for the consumption of lettuce irrigated with tertiary effluent from either WWTP was higher than the WHO recommendation of 10-6 DALYs for both viruses. The WWTP with constructed wetland showed a higher viral reduction on average (3.9 and 2.8 logs for NoV GII and HAdV, respectively) than conventional treatment (1.9 and 2.5 logs) but a higher variability than the conventional WWTP. Sensitivity analysis demonstrated that the input parameters used to estimate the viral reduction by treatment and viral concentrations accounted for much of the model output variability. The estimated reductions required to reach the WHO recommended levels in tertiary effluent are influenced by the characteristics of the treatments developed in the WWTPs, and additional average reductions are necessary (in WWTP with a constructed wetland: A total of 6.7 and 5.1 logs for NoV GII and HAdV, respectively; and in the more conventional treatment: 7 and 5.6 logs). This recommendation would be achieved with an average quantification of 0.5 genome copies per 100 mL in reclaimed water for both viruses. The results suggest that the analyzed reclaimed water would require additional treatments to achieve acceptable risk in the irrigation of vegetables with reclaimed water.

Genome characterization and temporal evolution analysis of a non-epidemic norovirus variant GII.8

Noroviruses are the primary cause of non-bacterial acute gastroenteritis worldwide, and GII.8 belongs to a non-epidemic genotype with a limited understanding currently. In this study, we assembled the first GII.8 norovirus genome from China and clarified the temporal evolutionary process of this non-epidemic variant. Using the "4+1+1" application strategy with newly designed primer sets, the genome of one GII.8 strain GZ2017-L601 from China was firstly sequenced that comprised 7476 nucleotides. The homology of the new genome and the previous only GII.8 genome reached 93.8% identity at the nucleotide level, but only 10, 6, 7 amino acid mutations occurred in three ORFs. When compared the new strain with other GII reference strains, p22 and P2 were calculated as the variable encoding regions, and NTPase, VPg, 3CL, RdRp and S were shown as the conserved ones. We then reconstructed the evolutionary process of the GII.8 genotype using other available sequences in GenBank. Based on the partial N/C region, all GII.8 strains could be subdivided chronologically into four clusters, which spans 1967-1994, 1997-2005, 2003-2009, and 2007-2017, respectively. Moreover, differences of capsid P proteins between GII.8 strains and the epidemic GII.4 strain VA387 were also compared. There existed 147/310 distinct amino acid sites in the alignment, including two insertion and three deletion mutations. Distribution of antigen epitopes of two GII.8 variants was comparable, but the numbers of antigenic sites of GII.8 strains were less than that of VA387. In summary, the first GII.8 genome from China was assembled and extensively characterized, and a time-order evolutionary process of this genotype was identified with a static pattern of antigenic variations.

Surveillance of norovirus contamination in commercial fresh/frozen berries from Heilongjiang Province, China, using a TaqMan real-time RT-PCR assay

Norovirus (NoV), a major food-borne virus, causes non-bacterial acute gastroenteritis in humans. Berries are generally harvested from low-growing bushes by hand and are minimally processed before being sold to consumers. Therefore, the consumption of berries has been linked to numerous food-borne gastroenteritis outbreaks caused by NoV in many countries. We performed a survey of NoV contamination in commercial fresh/frozen berry fruits collected from 2016 to 2017 in the Heilongjiang Province, the main berry-producing area in China, using a TaqMan-based real-time reverse transcription-PCR assay. Among 900 frozen and 900 fresh domestic retail berry samples, the prevalence of NoV was 9% (81/900) and 12.11% (109/900), including 35.80% (29/81) and 29.36% (32/109) of genotype GI alone, 54.32% (44/81) and 60.55% (66/109) of GII alone, and 9.88% (8/81) and 10.09% (11/109) of both GI and GII, respectively. No NoV was detected among the 677 frozen berry samples for export. Thus, the occurrence of NoV contamination was significantly higher in domestic berries than in exported berries and higher in fresh berries than in frozen berries. This study highlights the need for further risk surveillance for NoV contamination in berries produced in the Heilongjiang Province and recommends region-extended monitoring of retail berries for NoV.

Control of norovirus infection

PURPOSE OF REVIEW

The purpose of the review is to provide an update on control measures for norovirus (NoV), which is the most commonly implicated pathogen in acute gastroenteritis and outbreaks, causing major disruption in nurseries, schools, hospitals and care homes.

RECENT FINDINGS

Important developments include the discovery that virus particles, previously considered to be the infectious unit, also occur in clusters, which appear to be more virulent than individual virus particles; a working culture system using human stem-cell derived enteroids; promising results from early phase clinical trials of candidate NoV vaccines, which appear to be safe and immunogenic; chronic NoV affects patients with primary and secondary immune deficiencies. Although several treatments have been used none are supported by well designed clinical trials; infection control procedures are effective if properly implemented.

SUMMARY

NoV remains an important cause of morbidity and mortality. Although there are exciting developments on the vaccine front, the mainstay of control remains good hand hygiene, adherence to infection control procedures and limiting contamination of food, water and the wider environment. Once vaccines are available there will be important decisions to be made about how best to implement them.

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.

Emergence of norovirus GII.P16-GII.2 strains in patients with acute gastroenteritis in Huzhou, China, 2016-2017

Background

In late 2016, an uncommon recombinant NoV genotype called GII.P16-GII.2 caused a sharp increase in outbreaks of acute gastroenteritis in different countries of Asia and Europe, including China. However, we did not observe a drastic increase in sporadic norovirus cases in the winter of 2016 in Huzhou. Therefore, we investigate the prevalence and genetic diversity of NoVs in the sporadic acute gastroenteritis (AGE) cases from January 2016 to December 2017 in Huzhou City, Zhejiang, China.

Methods

From January 2016 to December 2017, a total of 1001 specimens collected from patients with AGE were screened for NoV by real-time RT-PCR. Partial sequences of the RNA-dependent RNA polymerase (RdRp) and capsid gene of the positive samples were amplified by RT-PCR and sequenced. Genotypes of NoV were confirmed by online NoV typing tool and phylogenetic analysis. Complete VP1 sequences of GII.P16-GII.2 strains detected in this study were further obtained and subjected into sequence analysis.

Results

In total, 204 (20.4%) specimens were identified as NoV-positive. GII genogroup accounted for most of the NoV-infected cases (98.0%, 200/204). NoV infection was found in all age groups tested (< 5, 5–15, 16–20, 21–30, 31–40, 41–50, 51–60, and >60 years), with the 5–15 year age group having the highest detection rate (17/49, 34.7%). Higher activity of NoV infection could be seen in winter-spring season. The predominant NoV genotypes have changed from GII.Pe-GII.4 Sydney2012 and GII.P17-GII.17 in 2016 to GII.P16-GII.2, GII.Pe-GII.4 Sydney2012 and GII.P17-GII.17 in 2017. Phylogenetic analyses revealed that 2016–2017 GII.P16-GII.2 strains were most closely related to Japan 2010–2012 cluster in VP1 region and no common mutations were found in the amino acids of the HBGA-binding sites and the predicted epitopes.

Conclusions

We report the emergence of GII.P16-GII.2 strains and characterize the molecular epidemiological patterns NoV infection between January 2016 and December 2017 in Huzhou. The predominant genotypes of NoV during our study period are diverse. VP1 amino acid sequences of 2016–2017 GII.P16-GII.2 strains remain static after one year of circulation.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3259-6) contains supplementary material, which is available to authorized users.

Intrahost Norovirus Evolution in Chronic Infection Over 5 Years of Shedding in a Kidney Transplant Recipient

Noroviruses are the leading cause of acute gastroenteritis, and they can affect humans of all age groups. In immunocompromised patients, norovirus infections can develop into chronic diarrhea or show prolonged asymptomatic virus shedding. Chronic norovirus infections are frequently reported for solid organ transplant recipients, with rapid intrahost norovirus evolution seen. In this report, we describe a case of chronic norovirus infection in an immunocompromised patient who was followed up for over 5 years. The purpose of the study was to specify the norovirus evolution in a chronically infected immunocompromised host and identify possible selection sites in norovirus capsid protein. During the follow-up period, 25 sequential stool samples were collected and nine of them were selected to generate amplicons covering viral RNA-dependent RNA polymerase (RdRp) and viral capsid protein (VP1) genes. Amplicons were sequenced using next-generation sequencing. Single nucleotide polymorphisms were defined, which demonstrated a nearly 3-fold greater mutation rate in the VP1 genome region compared to the RdRp genome region (7.9 vs. 2.8 variable sites/100 nucleotides, respectively). This indicates that mutations in the virus genome were not accumulated randomly, but are rather the result of mutant selection during the infection cycle. Using ShoRAH software we were able to reconstruct haplotypes occurring in each of the nine selected samples. The deduced amino-acid haplotype sequences were aligned and the positions were analyzed for selective pressure using the Datamonkey program. Only 12 out of 25 positive selection sites were within the commonly described epitopes A, B, C, and D of the VP1 protein. New positive selection sites were determined that have not been described before and might reflect adaptation of the norovirus toward optimal histo-blood-group antigen binding, or modification of the norovirus antigenic properties. These data provide new insights into norovirus evolutionary dynamics and indicate new putative epitope “hot-spots” of modified and optimized norovirus–host interactions.

Detailed genetic analyses of the HN gene in human respirovirus 3 detected in children with acute respiratory illness in the Iwate Prefecture, Japan

We performed detailed genetic analyses of the partial hemagglutinin-neuraminidase (HN) gene in 34 human respirovirus 3 (HRV3) strains from children with acute respiratory illness during 2013-2015 in Iwate Prefecture, Japan. In addition, we performed analyses of the evolutionary timescale of the gene using the Bayesian Markov chain Monte Carlo (MCMC) method. Furthermore, we analyzed pairwise distances and performed selective pressure analyses followed by linear B-cell epitope mapping and N-glycosylation and phylodynamic analyses. A phylogenetic tree showed that the strains diversified at around 1939, and the rate of molecular evolution was 7.6 × 10^-4 substitutions/site/year. Although the pairwise distances were relatively short (0.03 ± 0.018 [mean ± standard deviation, SD]), two positive selection sites (Cys544Trp and Leu555Ser) and no amino acid substitutions were found in the active/catalytic sites. Six epitopes were estimated in this study, and three mouse monoclonal antibody binding sites (amino acid positions 278, 281, and 461) overlapped with two epitopes belonging to subcluster C3 strains. Bayesian skyline plot analyses indicated that subcluster C3 strains have been increasing from 2004, whereas subcluster C1 strains have declined from 2004. Based on these results, Iwate strains were divided into two subclusters and each subcluster evolved independently. Moreover, our results suggested that some predicted linear epitopes (epitopes 3 and 5) are candidates for an HRV3 vaccine motif. To better understand the details of the molecular evolution of HRV, further studies are needed.

Genetic Analysis of Human Norovirus Strains in Japan in 2016-2017

In the 2016/2017 winter season in Japan, HuNoV GII.P16-GII.2 strains (2016 strains) emerged and caused large outbreaks of acute gastroenteritis. To better understand the outbreaks, we examined the molecular evolution of the VP1 gene and RdRp region in 2016 strains from patients by studying their time-scale evolutionary phylogeny, positive/negative selection, conformational epitopes, and phylodynamics. The time-scale phylogeny suggested that the common ancestors of the 2016 strains VP1 gene and RdRp region diverged in 2006 and 1999, respectively, and that the 2016 strain was the progeny of a pre-2016 GII.2. The evolutionary rates of the VP1 gene and RdRp region were around 10^-3 substitutions/site/year. Amino acid substitutions (position 341) in an epitope in the P2 domain of 2016 strains were not found in pre-2016 GII.2 strains. Bayesian skyline plot analyses showed that the effective population size of the VP1 gene in GII.2 strains was almost constant for those 50 years, although the number of patients with NoV GII.2 increased in 2016. The 2016 strain may be involved in future outbreaks in Japan and elsewhere.

Outbreaks of acute gastroenteritis associated with a re-emerging GII.P16-GII.2 norovirus in the spring of 2017 in Jiangsu, China

A total of 64 acute gastroenteritis outbreaks with 2,953 patients starting in December of 2016 and occurring mostly in the late spring of 2017 were reported in Jiangsu, China. A recombinant GII.P16-GII.2 norovirus variant was associated with 47 outbreaks (73.4%) for the gastroenteritis epidemic, predominantly occurring in February and March of 2017. Sequence analysis of the RNA-dependent RNA polymerase (RdRp) and capsid protein of the viral isolates from these outbreaks confirmed that this GII.P16-GII.2 strain was the GII.P16-GII.2 variant with the intergenotypic recombination, identified in Taiwan, Hong Kong, and other cities in China in 2016. This GII.P16-GII.2 recombinant variant appeared to a re-emerging strain, firstly identified in 2011-2012 from Japan and USA but might be independently originated from other GII.P16-GII.2 variants for sporadic and outbreaks of gastroenteritis in Japan and China before 2016. Further identification of unique amino acid mutations in both VP1 and RdRp of NoV strain as shown in this report may provide insight in explaining its structural and antigenic changes, potentially critical for the variant recombinant to gain its predominance in causing regional and worldwide epidemics.

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

Human norovirus (HuNoV) is a leading cause of viral gastroenteritis worldwide, of which GII.4 is the most predominant genotype. Unlike other genotypes, GII.4 has created various variants that escaped from previously acquired immunity of the host and caused repeated epidemics. However, the molecular evolutionary differences among all GII.4 variants, including recently discovered strains, have not been elucidated. Thus, we conducted a series of bioinformatic analyses using numerous, globally collected, full-length GII.4 major capsid (VP1) gene sequences (466 strains) to compare the evolutionary patterns among GII.4 variants. The time-scaled phylogenetic tree constructed using the Bayesian Markov chain Monte Carlo (MCMC) method showed that the common ancestor of the GII.4 VP1 gene diverged from GII.20 in 1840. The GII.4 genotype emerged in 1932, and then formed seven clusters including 14 known variants after 1980. The evolutionary rate of GII.4 strains was estimated to be 7.68 × 10⁻³ substitutions/site/year. The evolutionary rates probably differed among variants as well as domains [protruding 1 (P1), shell, and P2 domains]. The Osaka 2007 variant strains probably contained more nucleotide substitutions than any other variant. Few conformational epitopes were located in the shell and P1 domains, although most were contained in the P2 domain, which, as previously established, is associated with attachment to host factors and antigenicity. We found that positive selection sites for the whole GII.4 genotype existed in the shell and P1 domains, while Den Haag 2006b, New Orleans 2009, and Sydney 2012 variants were under positive selection in the P2 domain. Amino acid substitutions overlapped with putative epitopes or were located around the epitopes in the P2 domain. The effective population sizes of the present strains increased stepwise for Den Haag 2006b, New Orleans 2009, and Sydney 2012 variants. These results suggest that HuNoV GII.4 rapidly evolved in a few decades, created various variants, and altered its evolutionary rate and antigenicity.

Molecular Evolution of the RNA-Dependent RNA Polymerase and Capsid Genes of Human Norovirus Genotype GII.2 in Japan during 2004-2015

The RNA-dependent RNA polymerase (RdRp) and capsid (VP1) genes of 51 GII.2 human norovirus (HuNoV) strains collected during the period of 2004-2015 in Japan were analyzed. Full-length analyses of the genes were performed using next-generation sequencing. Based on the gene sequences, we constructed the time-scale evolutionary trees by Bayesian Markov chain Monte Carlo methods. Time-scale phylogenies showed that the RdRp and VP1 genes evolved uniquely and independently. Four genotypes of GII.2 (major types: GII.P2-GII.2 and GII.P16-GII.2) were detected. A common ancestor of the GII.2 VP1 gene existed until about 1956. The evolutionary rates of the genes were high (over 10^-3 substitutions/site/year). Moreover, the VP1 gene evolution may depend on the RdRp gene. Based on these results, we hypothesized that transfer of the RdRp gene accelerated the VP1 gene evolution of HuNoV genotype GII.2. Consequently, recombination between ORF1 (polymerase) and ORF2 (capsid) might promote changes of GII.2 antigenicity.

Molecular evolution of the fusion protein (F) gene in human respiratory syncytial virus subgroup B

In this study, we examined the molecular evolution of the fusion protein (F) gene in human respiratory syncytial virus subgroup B (HRSV-B). First, we performed time-scale evolution analyses using the Bayesian Markov chain Monte Carlo (MCMC) method. Next, we performed genetic distance, linear B-cell epitope prediction, N-glycosylation, positive/negative selection site, and Bayesian skyline plot analyses. We also constructed a structural model of the F protein and mapped the amino acid substitutions and the predicted B-cell epitopes. The MCMC-constructed phylogenetic tree indicated that the HRSV F gene diverged from the bovine respiratory syncytial virus gene approximately 580years ago and had a relatively low evolutionary rate (7.14×10^-4substitutions/site/year). Furthermore, a common ancestor of HRSV-A and -B diverged approximately 290years ago, while HRSV-B diverged into three clusters for approximately 60years. The genetic similarity of the present strains was very high. Although a maximum of 11 amino acid substitutions were observed in the structural model of the F protein, only one strain possessed an amino acid substitution located within the palivizumab epitope. Four epitopes were predicted, although these did not correspond to the neutralization sites of the F protein including the palivizumab epitope. In addition, five N-glycosylation sites of the present HRSV-B strains were inferred. No positive selection sites were identified; however, many sites were found to be under negative selection. The effective population size of the gene has remained almost constant. On the basis of these results, it can be concluded that the HRSV-B F gene is highly conserved, as is the F protein of HRSV-A. Moreover, our prediction of B-cell epitopes does not show that the palivizumab reaction site may be recognized as an epitope during naturally occurring infections.

Genetic analyses of the fusion protein genes in human parainfluenza virus types 1 and 3 among patients with acute respiratory infections in Eastern Japan from 2011 to 2015

PURPOSE

To genetically explore the fusion protein gene (F) in human parainfluenza virus type 1 (HPIV1) and type 3 (HPIV3) strains, we analysed them in patients with acute respiratory infections in Eastern Japan from 2011 to 2015.

METHODOLOGY

We constructed phylogenetic trees based on the HPIV and HPIV3 F gene using the maximum likelihood method and conducted P-distance and selective pressure analyses. We also predicted the linear epitopes of the protein in the prototype strains. Furthermore, we mapped the amino acid substitutions of the proteins.

RESULTS

Nineteen strains of HPIV1 and 53 strains of HPIV3 were detected among the clinical acute respiratory infection cases. The phylogenetic trees indicated that the HPIV1 and HPIV3 strains were classified into clusters II and III and cluster C, respectively. The P-distance values of the HPIV1 and HPIV3 F genes were <0.03. Two positive selection sites were inferred in the HPIV1 (aa 8 and aa 10), and one positive selection site was inferred in the HPIV3 (aa 108), but over 10 negative selection sites were inferred. Four epitopes were predicted for the HPIV1 prototype strains, while five epitopes were predicted for the HPIV3 prototype strain. A positive selection site (aa 108) or the HPIV3 F protein was involved in the predicted epitope. Additionally, we found that an amino acid substitution (R73K) in the LC76627 HPIV3 strain presumably may affect the resistance to neutralization by antibodies.

CONCLUSION

The F gene of HPIV1 and HPIV3 was relatively well conserved in the eastern part of Japan during the investigation period.

Complete Genome Sequence of a Recombinant GII.P16-GII.4 Norovirus Detected in Kawasaki City, Japan, in 2016

A recombinant norovirus, GII.P16-GII.4_Sydney2012, was first detected from nine patients with gastroenteritis in Kawasaki City, Japan, in 2016. The viral genome showed nucleotide sequence identities of 95.1% and 97.2% to the closest strains in the regions of 5′ terminus to ORF1 and ORF2 to 3′ terminus, respectively.