Comparative genomic analysis of genogroup 1 (Wa-like) rotaviruses circulating in the USA, 2006-2009

Whole-genome analysis of human group A rotaviruses in 1980s Japan and evolutionary assessment of global Wa-like strains across half a century

Historically, the Wa-like strains of human group A rotavirus (RVA) have been major causes of gastroenteritis. However, since the 2010s, the circulation of non-Wa-like strains has been increasingly reported, indicating a shift in the molecular epidemiology of RVA. Although understanding RVA evolution requires the analysis of both current and historical strains, comprehensive pre-1980's sequencing data are scarce globally. We determined the whole-genome sequences of representative strains from six RVA gastroenteritis outbreaks observed at an infant home in Sapporo, Japan, between 1981 and 1989. These outbreaks were mainly caused by G1 or G3 Wa-like strains, resembling strains from the United States in the 1970s-1980s and from Malawi in the 1990s. Phylogenetic analysis of these infant home strains, together with Wa-like strains collected worldwide from the 1970s to 2020, revealed a notable trend: pre-2010 strains diverged into multiple lineages in many genomic segments, whereas post-2010 strains tended to converge into a single lineage. However, Bayesian skyline plot indicated near-constant effective population sizes from the 1970s to 2020, and selection pressure analysis identified positive selection only at amino acid 75 of NSP2. These results suggest that evidence supporting the influence of rotavirus vaccines, introduced globally since 2006, on Wa-like RVA molecular evolution is lacking at present, and phylogenetic analysis may simply reflect natural fluctuations in RVA molecular evolution. Evaluating the long-term impact of RV vaccines on the molecular evolution of RVA requires sustained surveillance.

Whole-genome characterization of common rotavirus strains circulating in Vellore, India from 2002 to 2017: emergence of non-classical genomic constellations

Rotaviruses (RVs) are the most common etiological agent of acute gastroenteritis among young children, even after vaccine introduction in low-income countries. A whole-genome classification representing the 11 RV genes, was introduced for surveillance and characterization of RVs. This study characterized the common circulating strains in Vellore, India from 2002 to 2017 to understand rotavirus strain diversity and evolution using Whole genome sequencing (WGS) carried out on Illumina MiSeq. The 89% (92% of Wa-like, 86% of DS-1-like) of strains had classical constellations, while reassortant constellations were seen in 11% (8% of Wa-like, 14% of DS-1-like) of the strains. The rare E6-NSP4 in combination with DS-1 like G1P[8] and the emergence of the OP-354 subtype of P[8] were identified. Phylogenetics of RV strains revealed multiple subtypes circulating in the past 15 years, with strong evidence of animal to human gene transmission among several strains.

Rotavirus Strain Trends in United States, 2009-2016: Results from the National Rotavirus Strain Surveillance System (NRSSS)

Before the introduction of vaccines, group A rotaviruses (RVA) were the leading cause of acute gastroenteritis in children worldwide. The National Rotavirus Strain Surveillance System (NRSSS) was established in 1996 by the Centers for Disease Control and Prevention (CDC) to perform passive RVA surveillance in the USA. We report the distribution of RVA genotypes collected through NRSSS during the 2009-2016 RVA seasons and retrospectively examine the genotypes detected through the NRSSS since 1996. During the 2009-2016 RVA seasons, 2134 RVA-positive fecal specimens were sent to the CDC for analysis of the VP7 and VP4 genes by RT-PCR genotyping assays and sequencing. During 2009-2011, RVA genotype G3P[8] dominated, while G12P[8] was the dominant genotype during 2012-2016. Vaccine strains were detected in 1.7% of specimens and uncommon/unusual strains, including equine-like G3P[8] strains, were found in 1.9%. Phylogenetic analyses showed limited VP7 and VP4 sequence variation within the common genotypes with 1-3 alleles/lineages identified per genotype. A review of 20 years of NRSSS surveillance showed two changes in genotype dominance, from G1P[8] to G3P[8] and then G3P[8] to G12P[8]. A better understanding of the long-term effects of vaccine use on epidemiological and evolutionary dynamics of circulating RVA strains requires continued surveillance.

Comparative genomic analysis of genogroup 1 and genogroup 2 rotaviruses circulating in seven US cities, 2014-2016

For over a decade, the New Vaccine Surveillance Network (NVSN) has conducted active rotavirus (RVA) strain surveillance in the USA. The evolution of RVA in the post-vaccine introduction era and the possible effects of vaccine pressure on contemporary circulating strains in the USA are still under investigation. Here, we report the whole-gene characterization (eleven ORFs) for 157 RVA strains collected at seven NVSN sites during the 2014 through 2016 seasons. The sequenced strains included 52 G1P[8], 47 G12P[8], 18 G9P[8], 24 G2P[4], 5 G3P[6], as well as 7 vaccine strains, a single mixed strain (G9G12P[8]), and 3 less common strains. The majority of the single and mixed strains possessed a Wa-like backbone with consensus genotype constellation of G1/G3/G9/G12-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1, while the G2P[4], G3P[6], and G2P[8] strains displayed a DS-1-like genetic backbone with consensus constellation of G2/G3-P[4]/P[6]/P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Two intergenogroup reassortant G1P[8] strains were detected that appear to be progenies of reassortment events between Wa-like G1P[8] and DS-1-like G2P[4] strains. Two Rotarix^® vaccine (RV1) and two RV5 derived (vd) reassortant strains were detected. Phylogenetic and similarity matrices analysis revealed 2-11 sub-genotypic allelic clusters among the genes of Wa- and DS-1-like strains. Most study strains clustered into previously defined alleles. Amino acid (AA) substitutions occurring in the neutralization epitopes of the VP7 and VP4 proteins characterized in this study were mostly neutral in nature, suggesting that these RVA proteins were possibly under strong negative or purifying selection in order to maintain competent and actual functionality, but fourteen radical (AA changes that occur between groups) AA substitutions were noted that may allow RVA strains to gain a selective advantage through immune escape. The tracking of RVA strains at the sub-genotypic allele constellation level will enhance our understanding of RVA evolution under vaccine pressure, help identify possible mechanisms of immune escape, and provide valuable information for formulation of future RVA vaccines.

Amino Acid Substitutions in Positions 385 and 393 of the Hydrophobic Region of VP4 May Be Associated with Rotavirus Attenuation and Cell Culture Adaptation

Rotaviruses (RVs) are the leading cause of the acute viral gastroenteritis in young children and livestock animals worldwide. Although live attenuated vaccines have been applied to control RV infection for many years, the underlying mechanisms of RV attenuation following cell culture adaption are unknown. To study these mechanisms at the genomic level, we have sequenced and conducted a comparative analysis of two virulent human (Wa, G1P[8] and M, G3P[8]) and two virulent porcine (Gottfried, G4P[6] and OSU, G5P[7]) RV strains maintained in gnotobiotic piglets for 22, 11, 12 and 9 serial passages, respectively, with their attenuated counterparts serially passaged in MA-104 cell cultures for 25, 43, 54 and 43 passages, respectively. We showed that most of the mutations were clustered in the VP4 gene, with a relatively high nonsynonymous substitution rate (81.2%). Moreover, two amino acid substitutions observed in the VP4 gene were conserved between two or more strain pairs. D385N substitution was found in M, Wa and Gottfried strains, and another one, S471H/L was present in Wa and Gottfried strains. Importantly, D385 was reported previously in another study and may be involved in regulation of virus entry. Of interest, although no 385 substitution was found in OSU strains, the attenuated OSU strain contained a unique D393H substitution within the same VP4 hydrophobic domain. Collectively, our data suggest that the VP4 hydrophobic region may play an important role in RV attenuation and aa385 and aa393 may represent potential targets for RV vaccine development using reverse genetics and site-specific mutagenesis.

Molecular characterisation of wild-type G1P[8] and G3P[8] rotaviruses isolated in Vietnam 2008 during a vaccine trial

Rotavirus vaccines work better in developed countries than in developing countries, leading to the question of whether the circulating strains are different in these two settings. In 2008, a clinical trial of the pentavalent rotavirus vaccine was performed in Nha Trang, Vietnam, in which the efficacy was reported to be 64 %. Although samples were collected independently from the clinical trial, we examined faecal specimens from children hospitalised for rotavirus diarrhoea and found that G3P[8] and G1P[8] were co-dominant at the time of the clinical trial. The aim of this study was to explore whether they were divergent from the strains circulating in the developed countries where the vaccine efficacy is high. Two G3P[8] and two G1P[8] strains that were regarded as representatives based on their electropherotypes were selected for full-genome sequencing. The genotype constellation was G1/G3-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1. All but the VP4 genes, one of which belonged to the emerging P[8]b genotype (OP354-like VP4), clustered into one or more lineages/alleles with the strains circulating in developed countries, with ≥97.5 % nucleotide sequence identity. Additionally, 10 G1 and 12 G3 VP7 sequences as well as 31 VP4 sequences were determined. No amino acid differences were observed between the Vietnamese strains and strains in the developed countries that were likely to have affected the neutralisation specificity of their VP7 and VP4. In conclusion, apart from prevalent P[8]b VP4, virtually no differences were observed between the predominant strains circulating in Vietnam at the time of the clinical trial and the strains in the developed countries; hence, the lower vaccine efficacy was more likely to be due to factors other than strain divergence.

Genome-Wide Evolutionary Analyses of G1P[8] Strains Isolated Before and After Rotavirus Vaccine Introduction

Rotaviruses are the most important etiological agent of acute gastroenteritis in young children worldwide. Among the first countries to introduce rotavirus vaccines into their national immunization programs were Belgium (November 2006) and Australia (July 2007). Surveillance programs in Belgium (since 1999) and Australia (since 1989) offer the opportunity to perform a detailed comparison of rotavirus strains circulating pre- and postvaccine introduction. G1P[8] rotaviruses are the most prominent genotype in humans, and a total of 157 G1P[8] rotaviruses isolated between 1999 and 2011 were selected from Belgium and Australia and their complete genomes were sequenced. Phylogenetic analysis showed evidence of frequent reassortment among Belgian and Australian G1P[8] rotaviruses. Although many different phylogenetic subclusters were present before and after vaccine introduction, some unique clusters were only identified after vaccine introduction, which could be due to natural fluctuation or the first signs of vaccine-driven evolution. The times to the most recent common ancestors for the Belgian and Australian G1P[8] rotaviruses ranged from 1846 to 1955 depending on the gene segment, with VP7 and NSP4 resulting in the most recent estimates. We found no evidence that rotavirus population size was affected after vaccine introduction and only six amino acid sites in VP2, VP3, VP7, and NSP1 were identified to be under positive selective pressure. Continued surveillance of G1P[8] strains is needed to determine long-term effects of vaccine introductions, particularly now rotavirus vaccines are implemented in the national immunization programs of an increasing number of countries worldwide.