Molecular Evolution of the Capsid Gene in Norovirus Genogroup I

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.

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.

Illuminating Human Norovirus: A Perspective on Disinfection of Water and Surfaces Using UVC, Norovirus Model Organisms, and Radiation Safety Considerations

Human noroviruses (HuNoVs) are a major cause of gastroenteritis and are associated with high morbidity because of their ability to survive in the environment and small inoculum size required for infection. Norovirus is transmitted through water, food, high touch-surfaces, and human-to-human contact. Ultraviolet Subtype C (UVC) light-emitting diodes (LEDs) can disrupt the norovirus transmission chain for water, food, and surfaces. Here, we illuminate considerations to be adhered to when picking norovirus surrogates for disinfection studies and shine light on effective use of UVC for norovirus infection control in water and air and validation for such systems and explore the blind spot of radiation safety considerations when using UVC disinfection strategies. This perspective also discusses the promise of UVC for norovirus mitigation to save and ease life.

Evolutionary dynamics of influenza B strains detected from paediatric acute respiratory infections in central Vietnam

Influenza virus B belongs to the family Orthomyxoviridae with segmented negative-sense RNA genomes. Since 1970s, influenza B has diverged intoVictoria and Yamagata, which differs in antigenic and evolutionary characteristics. Yet, molecular-epidemiological information of influenza B from developing nations is limited. In central Vietnam, influenza A subtype-specific circulation pattern and clinical characteristics were previously described. However, molecular evolutionary characteristics of influenza B has not been discussed to date. We utilized the influenza B positives obtained from paediatric ARI surveillance during 2007-2013. Influenza B HA and NA genes were amplified, sequenced, and phylogenetic/molecular evolutionary analysis was performed using Maximum Likelihood and Bayesian MCMC. Phylodynamics analysis was performed with Bayesian Skyline Plot (BSP). Furthermore, we performed selection pressure analysis and estimated N-glycosylation sites. In the current study, overall positive rate for influenza B was 3.0%, and Victoria lineage immediately became predominant in post-A/H1N1pdm09 period. The noticeable shift in Victoria lineage WHO Group occurred. With respect to the evolutionary rate (substitutions/site/year), Victoria lineage HA gene was evolving faster than Yamagata lineage (2.43 × 10^-3 vs 2.00 × 10^-3). Furthermore, the evolutionary rate of Victoria Group 5 was greater than Group 1. BSP presented the rapid growth in Effective Population Size (EPS) of Victoria lineage occurred soon after the 1st A/H1N1pdm09 case was detected whereas the EPS of Yamagata lineage was stable for both genes. N-glycosylation pattern between lineages and among WHO Groups were slightly different, and HA gene had a total of 6 amino acid substitutions under positive section pressure (4 for Victoria and 2 for Yamagata). The current results highlight the importance of Victoria lineage in post-A/H1N1pdm09 period. Difference in evolutionary characteristics and phylodynamics may indicate lineage and WHO Group-specific evolutionary dynamics. It is necessary to further continue the molecular-epidemiological surveillance in local setting to gain a better understanding of local evolutionary characteristics of influenza B strains.

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.

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.

Molecular evolutionary analysis of type-1 human astroviruses identifies putative sites under selection pressure on the capsid protein

Human astroviruses (HAstV) are important enteric pathogens that can be classified into eight sero/genotypes (HAstV-1 to -8). Although the various HAstV types show global spread, type-1 strains tend to be predominant. Molecular analysis of the genomic region encoding the capsid protein (ORF2) has revealed discrete sequence variation, with different lineages within each HAstV type and at least three major lineages have been identified within HAstV-1. Longitudinal epidemiological surveillance has revealed temporal shift of the various HAstV-1 lineages. Metadata analysis of HAstV-1 sequences available in the databases also revealed temporal shifts of the circulation of HAstV-1 lineages, suggesting possible antigenic-related mechanisms of selection at the sub-genotype level. By comparison of HAstV-1 capsid sequences, lineage-defining residues under positive selection were identified. Structural analysis of HAstV-1 capsid allowed identifying at least six residues exposed on the virion surface. Two residues were located in the VP34 (shell region) whilst four residues were mapped in the VP25/27 (protruding region) of HAstV capsid protein, in proximity of the putative receptor binding S site. These findings suggest that mechanisms similar to those observed and/or hypothesized for other enteric viruses are also shaping the evolution of HAstVs, with intra-typic diversification being a possible mechanism to decrease the antigenic pressure to which these viruses are exposed.

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.

Predicting genotype compositions in norovirus seasons in Japan

Noroviruses cause acute gastroenteritis. Since multiple genotypes of norovirus co-circulate in humans, changing the genotype composition and eluding host immunity, development of a polyvalent vaccine against norovirus in which the genotypes of vaccine strains match the major strains in circulation in the target season is desirable. However, this would require prediction of changes in the genotype composition of circulating strains. A fitness model that predicts the proportion of a strain in the next season from that in the current season has been developed for influenza A virus. Here, such a fitness model that takes into account the fitness effect of herd immunity was used to predict genotype compositions in norovirus seasons in Japan. In the current study, a model that assumes a decline in the magnitude of cross immunity between norovirus strains according to an increase in the divergence of the major antigenic protein VP1 was found to be appropriate for predicting genotype composition. Although it is difficult to predict the proportions of genotypes accurately, the model is effective in predicting the direction of change in the proportions of genotypes. The model predicted that GII.3 and GII.4 may contract, whereas GII.17 may expand and predominate in the 2015-2016 season. The procedure of predicting genotype compositions in norovirus seasons described in the present study has been implemented in the norovirus forecasting system (NOROCAST).

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.

Static and Evolving Norovirus Genotypes: Implications for Epidemiology and Immunity

Noroviruses are major pathogens associated with acute gastroenteritis worldwide. Their RNA genomes are diverse, with two major genogroups (GI and GII) comprised of at least 28 genotypes associated with human disease. To elucidate mechanisms underlying norovirus diversity and evolution, we used a large-scale genomics approach to analyze human norovirus sequences. Comparison of over 2000 nearly full-length ORF2 sequences representing most of the known GI and GII genotypes infecting humans showed a limited number (≤5) of distinct intra-genotypic variants within each genotype, with the exception of GII.4. The non-GII.4 genotypes were comprised of one or more intra-genotypic variants, with each variant containing strains that differed by only a few residues over several decades (remaining "static") and that have co-circulated with no clear epidemiologic pattern. In contrast, the GII.4 genotype presented the largest number of variants (>10) that have evolved over time with a clear pattern of periodic variant replacement. To expand our understanding of these two patterns of diversification ("static" versus "evolving"), we analyzed using NGS the nearly full-length norovirus genome in healthy individuals infected with GII.4, GII.6 or GII.17 viruses in different outbreak settings. The GII.4 viruses accumulated mutations rapidly within and between hosts, while the GII.6 and GII.17 viruses remained relatively stable, consistent with their diversification patterns. Further analysis of genetic relationships and natural history patterns identified groupings of certain genotypes into larger related clusters designated here as "immunotypes". We propose that "immunotypes" and their evolutionary patterns influence the prevalence of a particular norovirus genotype in the human population.

Molecular evolution of respiratory syncytial virus subgroup A genotype NA1 and ON1 attachment glycoprotein (G) gene in central Vietnam

We performed molecular evolutionary analyses of the G gene C-terminal 3rd hypervariable region of RSV-A genotypes NA1 and ON1 strains from the paediatric acute respiratory infection patients in central Vietnam during the 2010-2012 study period. Time-scaled phylogenetic analyses were performed using Bayesian Markov Chain Monte Carlo (MCMC) method, and pairwise distances (p-distances) were calculated. Bayesian Skyline Plot (BSP) was constructed to analyze the time-trend relative genetic diversity of central Vietnam RSV-A strains. We also estimated the N-glycosylation sites within G gene hypervariable region. Amino acid substitutions under positive and negative selection pressure were examined using Conservative Single Likelihood Ancestor Counting (SLAC), Fixed Effects Likelihood (FEL), Internal Fixed Effects Likelihood (IFEL) and Mixed Effects Model for Episodic Diversifying Selection (MEME) models. The majority of central Vietnam ON1 strains detected in 2012 were classified into lineage 1 with few positively selected substitutions. As for the Vietnamese NA1 strains, four lineages were circulating during the study period with a few positive selection sites. Shifting patterns of the predominantly circulating NA1 lineage were observed in each year during the investigation period. Median p-distance of central Vietnam NA1 strains was wider (p-distance=0.028) than that of ON1 (p-distance=0.012). The molecular evolutionary rate of central Vietnam ON1 strains was estimated to be 2.55×10^-2 (substitutions/site/year) and was faster than NA1 (7.12×10^-3 (substitutions/site/year)). Interestingly, the evolutionary rates of both genotypes ON1 and NA1 strains from central Vietnam were faster than the global strains respectively. Furthermore, the shifts of N-glycosylation pattern within the G gene 3rd hypervariable region of Vietnamese NA1 strains were observed in each year. BSP analysis indicated the rapid growth of RSV-A effective population size in early 2012. These results suggested that the molecular evolution of RSV-A G gene detected in central Vietnam was fast with unique evolutionary dynamics.

Molecular evolution of the capsid gene in human norovirus genogroup II

Capsid protein of norovirus genogroup II (GII) plays crucial roles in host infection. Although studies on capsid gene evolution have been conducted for a few genotypes of norovirus, the molecular evolution of norovirus GII is not well understood. Here we report the molecular evolution of all GII genotypes, using various bioinformatics techniques. The time-scaled phylogenetic tree showed that the present GII strains diverged from GIV around 1630CE at a high evolutionary rate (around 10(-3) substitutions/site/year), resulting in three lineages. The GII capsid gene had large pairwise distances (maximum > 0.39). The effective population sizes of the present GII strains were large (>10(2)) for about 400 years. Positive (20) and negative (over 450) selection sites were estimated. Moreover, some linear and conformational B-cell epitopes were found in the deduced GII capsid protein. These results suggested that norovirus GII strains rapidly evolved with high divergence and adaptation to humans.

Molecular evolution of the fusion protein gene in human respiratory syncytial virus subgroup A

We studied the molecular evolution of the fusion protein (F) gene in the human respiratory syncytial virus subgroup A (HRSV-A). We performed time-scaled phylogenetic analyses using the Bayesian Markov chain Monte Carlo (MCMC) method. We also conducted genetic distance (p-distance), positive/negative selection, and Bayesian skyline plot analyses. Furthermore, we mapped the amino acid substitutions of the protein. The MCMC-constructed tree indicated that the HRSV F gene diverged from the bovine RSV (BRSV) gene approximately 550years ago and had a relatively low substitution rate (7.59×10(-4) substitutions/site/year). Moreover, a common ancestor of HRSV-A and -B diverged approximately 280years ago, which has since formed four distinct clusters. The present HRSV-A strains were assigned six genotypes based on F gene sequences and attachment glycoprotein gene sequences. The present strains exhibited high F gene sequence similarity values and low genetic divergence. No positive selection sites were identified; however, 50 negative selection sites were identified. F protein amino acid substitutions at 17 sites were distributed in the F protein. The effective population size of the gene has remained relatively constant, but the population size of the prevalent genotype (GA2) has increased in the last 10years. These results suggest that the HRSV-AF gene has evolved independently and formed some genotypes.

Human norovirus hyper-mutation revealed by ultra-deep sequencing

Human noroviruses (NoVs) are a major cause of gastroenteritis worldwide. It is thought that, similar to other RNA viruses, high mutation rates allow NoVs to evolve fast and to undergo rapid immune escape at the population level. However, the rate and spectrum of spontaneous mutations of human NoVs have not been quantified previously. Here, we analyzed the intra-patient diversity of the NoV capsid by carrying out RT-PCR and ultra-deep sequencing with 100,000-fold coverage of 16 stool samples from symptomatic patients. This revealed the presence of low-frequency sequences carrying large numbers of U-to-C or A-to-G base transitions, suggesting a role for hyper-mutation in NoV diversity. To more directly test for hyper-mutation, we performed transfection assays in which the production of mutations was restricted to a single cell infection cycle. This confirmed the presence of sequences with multiple U-to-C/A-to-G transitions, and suggested that hyper-mutation contributed a large fraction of the total NoV spontaneous mutation rate. The type of changes produced and their sequence context are compatible with ADAR-mediated editing of the viral RNA.

•

Norovirus U-to-C hyper-mutants are present in patient samples.

•

Analysis of hyper-mutants in cell culture suggests ADAR-mediated RNA edition.

•

Hyper-mutation may contribute to norovirus diversity and evolution.