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

Comparative analysis of the RVA VP7 and VP4 antigenic epitopes circulating in Iran and the Rotarix and RotaTeq vaccines

Analyzing the lineages and detecting antigenic variation in immunogenic motifs of Group A Rotavirus (RVA) variants is crucial because it can impact vaccine efficacy. This study investigated the circulating lineages of VP4 and VP7 proteins of human RVA isolates and their phylogeny in ≤24-month-old symptomatic, rotavirus-positive children with transudative diarrhea within 48 h of admission to Mofid Children's Hospital between December 2020 and March 2022 in Tehran, Iran. Antigen detection was performed by ELISA, RNA extraction, and semi-nested multiplex PCR for G/P genotypes, followed by sequencing and bioinformatic analysis using multiple sequence alignments in MEGA and phylogenetic analysis by Geneious Prime. The similarity of VP7 and VP4 amino acids with the RotaTeq and Rotarix vaccine strains for cytotoxic T cell and antigenic epitopes was evaluated using the UCSF Chimera Molecular Modeling System. Overall, 27.3 % of the samples were RVA positive, showing untypeable (2.5 %), single (76.9 %), and mixed (20.5 %) genotypic characteristics. The strains clustered in the G1/II, G2/IV, G3/I, G4/I, G9/III, P (Kachooei et al., 2023) [8]/III, P (Howley et al., 2020) [4]/V, and P (Wahyuni et al., 2021) [6]/I lineages. Comparative analysis of VP7 antigenic epitopes showed that the G1/II strains were completely conserved, while the G2/IV, G3/I, G4/I, G6, G9/III strains contained 2, 3-5, 2, 4 and 9 amino acid substitutions, respectively. The P (Kachooei et al., 2023) [8]/III genotypes differed by 3 amino acids, while the P (Wahyuni et al., 2021) [6]/I genotype had the most substitutions. CTL epitopes were completely conserved in G3/I strains, but other genotypes differed by 1-4 amino acids compared to the vaccine strains. Given the diversity of circulating RVA genotypes and the observed mutations in neutralizing and CTL epitopes, immune escape by some of the strains is likely in Iran. This finding underscores the importance of evaluating the effect of rotavirus vaccines on local genotypes and related lineages before implementing a vaccination program.

Rotavirus Genotypes in the Postvaccine Era: A Systematic Review and Meta-analysis of Global, Regional, and Temporal Trends by Rotavirus Vaccine Introduction

BACKGROUND

Even moderate differences in rotavirus vaccine effectiveness against nonvaccine genotypes may exert selective pressures on circulating rotaviruses. Whether this vaccine effect or natural temporal fluctuations underlie observed changes in genotype distributions is unclear.

METHODS

We systematically reviewed studies reporting rotavirus genotypes from children <5 years of age globally between 2005 and 2023. We compared rotavirus genotypes between vaccine-introducing and nonintroducing settings globally and by World Health Organization (WHO) region, calendar time, and time since vaccine introduction.

RESULTS

Crude pooling of genotype data from 361 studies indicated higher G2P[4], a nonvaccine genotype, prevalence in vaccine-introducing settings, both globally and by WHO region. This difference did not emerge when examining genotypes over time in the Americas, the only region with robust longitudinal data. Relative to nonintroducing settings, G2P[4] detections were more likely in settings with recent introduction (eg, 1-2 years postintroduction adjusted odds ratio [aOR], 4.39; 95% confidence interval [CI], 2.87-6.72) but were similarly likely in settings with more time elapsed since introduction, (eg, 7 or more years aOR, 1.62; 95% CI, .49-5.37).

CONCLUSIONS

When accounting for both regional and temporal trends, there was no substantial evidence of long-term vaccine-related selective pressures on circulating genotypes. Increased prevalence of G2P[4] may be transient after rotavirus vaccine introduction.

Genetic diversity of G9, G3, G8 and G1 rotavirus group A strains circulating among children with acute gastroenteritis in Vietnam from 2016 to 2021

Rotavirus group A (RVA) is the most common cause of severe childhood diarrhea worldwide. The introduction of rotavirus vaccination programs has contributed to a reduction in hospitalizations and mortality caused by RVA. From 2016 to 2021, we conducted surveillance to monitor RVA prevalence and genotype distribution in Nam Dinh and Thua Thien Hue (TT Hue) provinces where a pilot Rotavin-M1 vaccine (Vietnam) implementation took place from 2017 to 2020. Out of 6626 stool samples, RVA was detected in 2164 (32.6%) by ELISA. RT-PCR using type-specific primers were used to determine the G and P genotypes of RVA-positive specimens. Whole genome sequences of a subset of 52 specimens randomly selected from 2016 to 2021 were mapped using next-generation sequencing. From 2016 to 2021, the G9, G3 and G8 strains dominated, with detected frequencies of 39%, 23%, and 19%, respectively; of which, the most common genotypes identified were G9P[8], G3P[8] and G8P[8]. G1 strains re-emerged in Nam Dinh and TT Hue (29.5% and 11.9%, respectively) from 2020 to 2021. G3 prevalence decreased from 74% to 20% in TT Hue and from 21% to 13% in Nam Dinh province between 2017 and 2021. The G3 strains consisted of 52% human typical G3 (hG3) and 47% equine-like G3 (eG3). Full genome analysis showed substantial diversity among the circulating G3 strains with different backgrounds relating to equine and feline viruses. G9 prevalence decreased sharply from 2016 to 2021 in both provinces. G8 strains peaked during 2019-2020 in Nam Dinh and TT Hue provinces (68% and 46%, respectively). Most G8 and G9 strains had no genetic differences over the surveillance period with very high nucleotide similarities of 99.2-99.9% and 99.1-99.7%, respectively. The G1 strains were not derived from the RVA vaccine. Changes in the genotype distribution and substantial diversity among circulating strains were detected throughout the surveillance period and differed between the two provinces. Determining vaccine effectiveness against circulating strains over time will be important to ensure that observed changes are due to natural secular variation and not from vaccine pressure.

Characterization of Rotavirus Strains Responsible for Breakthrough Diarrheal Diseases among Zambian Children Using Whole Genome Sequencing

The occurrence of rotavirus (RV) infection among vaccinated children in high-burden settings poses a threat to further disease burden reduction. Genetically altered viruses have the potential to evade both natural infection and vaccine-induced immune responses, leading to diarrheal diseases among vaccinated children. Studies characterizing RV strains responsible for breakthrough infections in resource-limited countries where RV-associated diarrheal diseases are endemic are limited. We aimed to characterize RV strains detected in fully vaccinated children residing in Zambia using next-generation sequencing. We conducted whole genome sequencing on Illumina MiSeq. Whole genome assembly was performed using Geneious Prime 2023.1.2. A total of 76 diarrheal stool specimens were screened for RV, and 4/76 (5.2%) were RV-positive. Whole genome analysis revealed RVA/Human-wt/ZMB/CIDRZ-RV2088/2020/G1P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2 and RVA/Human-wt/ZMB/CIDRZ-RV2106/2020/G12P[4]-I1-R2-C2-M2-A2-N1-T2-E1-H2 strains were mono and multiple reassortant (exchanged genes in bold) respectively, whilst RVA/Human-wt/ZMB/CIDRZ-RV2150/2020/G12P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 was a typical Wa-like strain. Comparison of VP7 and VP4 antigenic epitope of breakthrough strains and Rotarix strain revealed several amino acid differences. Variations in amino acids in antigenic epitope suggested they played a role in immune evasion of neutralizing antibodies elicited by vaccination. Findings from this study have the potential to inform national RV vaccination strategies and the design of highly efficacious universal RV vaccines.

Coding-complete genome sequences of rotavirus A reference strains EDIM, Ph158, and CC425

This study reports the coding-complete genome sequences of three rotavirus A (RVA) reference strains previously adapted in tissue culture: RVA/Mouse-tc/USA/EDIM/XXXX/G16P[16] with a G16-P[16]-I7-R7-C7-M8-A7-N7-T10-E7-H9 genotype constellation, RVA/Human-tc/USA/Ph158/1998/G9P[6] with a G9-P[6]-I2-R2-C2-M2-A2-N2-T2-E2-H2 genotype constellation, and RVA/Human-tc/USA/CC425/1998/G3P[9] with a G3-P[9]-I2-R2-C2-M2-A3-N2-T1-E2-H3 genotype constellation.

Full genome-based characterization of an Asian G3P[6] human rotavirus strain found in a diarrheic child in Japan: Evidence for porcine-to-human zoonotic transmission

Human rotavirus strains having the unconventional G3P[6] genotype have been sporadically detected in diarrheic patients in different parts of the world. However, the full genomes of only three human G3P[6] strains from Asian countries (China, Indonesia, and Vietnam) have been sequenced and characterized, and thus the exact origin and evolution of G3P[6] strains in Asia remain to be elucidated. Here, we sequenced and characterized the full genome of a G3P[6] strain (RVA/Human-wt/JPN/SO1199/2020/G3P[6]) found in a stool sample from a 3-month-old infant admitted with acute gastroenteritis in Japan. On full genomic analysis, strain SO1199 was revealed to have a unique Wa-like genogroup configuration: G3-P[6]-I5-R1-C1-M1-A8-N1-T1-E1-H1. VP6 genotype I5 and NSP1 genotype A8 are commonly found in porcine rotavirus strains. Furthermore, phylogenetic analysis demonstrated that all 11 genes of strain SO1199 were closely related to those of porcine and/or porcine-like human rotaviruses and thus appeared to be of porcine origin. Thus, strain SO1199 was shown to possess a porcine-like genomic backbone and thus is likely to be the result of interspecies transmission of a porcine rotavirus strain. Of note is that all 11 genes of strain SO1199 were phylogenetically located in clusters, distinct from those of the previously identified porcine-like human G3P[6] strains from around the world including Asia, suggesting the occurrence of independent porcine-to-human zoonotic transmission events. To our knowledge, this is the first report on full genome-based characterization of a human G3P[6] strain that has emerged in Japan. Our findings revealed the diversity of unconventional human G3P[6] strains in Asia, and provide important insights into the origin and evolution of G3P[6] strains.

Genotypic investigation of a rotavirus cluster at a quaternary-care pediatric hospital

Rotavirus (RV) was a common healthcare-associated infection prior to the introduction of the RV vaccine. Following widespread RV vaccination, healthcare-associated rotavirus cases are rare. We describe an investigation of a cluster of rotavirus infections in a pediatric hospital in which an uncommon genotype not typically circulating in the United States was detected.

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.

Increase in rotavirus prevalence with the emergence of genotype G9P[8] in replacement of genotype G12P[6] in Sabah, Malaysia

Rotaviruses are major causative agents of acute diarrhea in children under 5 years of age in Malaysia. However, a rotavirus vaccine has not been included in the national vaccination program. To date, only two studies have been carried out in the state of Sabah, Malaysia, although children in this state are at risk of diarrheal diseases. Previous studies showed that 16%-17% of cases of diarrhea were caused by rotaviruses and that equine-like G3 rotavirus strains are predominant. Because the prevalence of rotaviruses and their genotype distribution vary over time, this study was conducted at four government healthcare facilities from September 2019 through February 2020. Our study revealed that the proportion of rotavirus diarrhea increased significantly to 37.2% (51/137) after the emergence of the G9P[8] genotype in replacement of the G12P[8] genotype. Although equine-like G3P[8] strains remain the predominant rotaviruses circulating among children, the Sabahan G9P[8] strain belonged to lineage VI and was phylogenetically related to strains from other countries. A comparison of the Sabahan G9 strains with the G9 vaccine strains used in the RotaSiil and Rotavac vaccines revealed several mismatches in neutralizing epitopes, indicating that these vaccines might not be effective in Sabahan children. However, a vaccine trial may be necessary to understand the precise effects of vaccination.

Coding-Complete Genome Sequences of G6P[14] Rotavirus Strain Detected in a Human Stool Specimen within the United States

In this study, we report the detection of a G6P[14] rotavirus strain from a human stool sample within the United States. The full genotype constellation of the G6P[14] strain was identified as G6-P[14]-I2-R2-C2-M2-A11-N2-T6-E2-H3.

Full genotype constellations analysis of unusual DS-1-like G12P[6] and G6P[8] rotavirus strains detected in Brazil, 2019

Rotavirus A (RVA) is a major cause of acute gastroenteritis (AGE) in children worldwide. We report unusual RVA G12P[6] and G6P[8] strains isolated from fecal samples from Brazilian children hospitalized for AGE. The characterized RVA have genome segments backbone: G12-P[6]/ G6-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2 of DS-1-like genogroup. Our study describes the first identification of G6P[8], a DS-1-like genogroup strain. Nucleotide analysis of VP7 and VP4 genes revealed that all G12 Brazilian strains clustered into the sub-lineages IIIB, mostly associated with P[6] lineage I. Additionally, our G6 lineage I strains were closely related to German G6 genotypes, bound with P[8] lineage III, differing from both vaccine strains. The comparative sequence analysis of our strains with vaccine strains revealed amino acid substitutions located in immunodominant regions of VP7 and VP4 proteins. Continuous monitoring of RVA genotypes is essential to evaluate the impact of vaccination on the dynamic nature of RVA evolution.

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.

Differences in epidemiology of enteropathogens in children pre- and post-rotavirus vaccine introduction in Kilifi, coastal Kenya

BACKGROUND

Kenya introduced Rotarix® (GlaxoSmithKline Biologicals, Rixensart, Belgium) vaccination into its national immunization programme beginning July 2014. The impact of this vaccination program on the local epidemiology of various known enteropathogens is not fully understood.

METHODS

We used a custom TaqMan Array Card (TAC) to screen for 28 different enteropathogens in 718 stools from children aged less than 13 years admitted to Kilifi County Hospital, coastal Kenya, following presentation with diarrhea in 2013 (before vaccine introduction) and in 2016-2018 (after vaccine introduction). Pathogen positivity rate differences between pre- and post-Rotarix® vaccination introduction were examined using both univariate and multivariable logistic regression models.

RESULTS

In 665 specimens (92.6%), one or more enteropathogen was detected, while in 323 specimens (48.6%) three or more enteropathogens were detected. The top six detected enteropathogens were: enteroaggregative Escherichia coli (EAggEC; 42.1%), enteropathogenic Escherichia coli (EPEC; 30.2%), enterovirus (26.9%), rotavirus group A (RVA; 24.8%), parechovirus (16.6%) and norovirus GI/GII (14.4%). Post-rotavirus vaccine introduction, there was a significant increase in the proportion of samples testing positive for EAggEC (35.7% vs. 45.3%, p = 0.014), cytomegalovirus (4.2% vs. 9.9%, p = 0.008), Vibrio cholerae (0.0% vs. 2.3%, p = 0.019), Strongyloides species (0.8% vs. 3.6%, p = 0.048) and Dientamoeba fragilis (2.1% vs. 7.8%, p = 0.004). Although not reaching statistical significance, the positivity rate of adenovirus 40/41 (5.8% vs. 7.3%, p = 0.444), norovirus GI/GII (11.2% vs. 15.9%, p = 0.089), Shigella species (8.7% vs. 13.0%, p = 0.092) and Cryptosporidium spp. (11.6% vs. 14.7%, p = 0.261) appeared to increase post-vaccine introduction. Conversely, the positivity rate of sapovirus decreased significantly post-vaccine introduction (7.8% vs. 4.0%, p = 0.030) while that of RVA appeared not to change (27.4% vs. 23.5%, p = 0.253). More enteropathogen coinfections were detected per child post-vaccine introduction compared to before (mean: 2.7 vs. 2.3; p = 0.0025).

CONCLUSIONS

In this rural Coastal Kenya setting, childhood enteropathogen infection burden was high both pre- and post-rotavirus vaccination introduction. Children who had diarrheal admissions post-vaccination showed an increase in coinfections and changes in specific enteropathogen positivity rates. This study highlights the utility of multipathogen detection platforms such as TAC in understanding etiology of childhood acute gastroenteritis in resource-limited regions.

Whole genome analysis of rotavirus strains circulating in Benin before vaccine introduction, 2016-2018

Species A Rotaviruses (RVA) still play a major role in causing acute diarrhea in children under five years old worldwide. Currently, an 11-gene classification system is used to designate the full genotypic constellations of circulating strains. Viral proteins and non-structural proteins in the order VP7-VP4-VP6-VP1-VP2-VP3-NSP1-NSP2-NSP3-NSP4-NSP5/6 are represented by the genotypes Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx, respectively. In Benin, ROTAVAC® vaccine was introduced into the Expanded Programme on Immunization in December 2019. To monitor circulating RVA strains for changes that may affect vaccine performance, in-depth analysis of strains prior to vaccine introduction are needed. Here we report, the whole-gene characterization (11 ORFs) for 72 randomly selected RVA strains of common and unusual genotypes collected in Benin from the 2016-2018 seasons. The sequenced strains were 15 G1P[8], 20 G2P[4], 5 G9P[8], 14 G12P[8], 9 G3P[6], 2 G1P[6], 3 G2P[6], 2 G9P[4], 1 G12P[6], and 1 G1G9P[8]/P[4]. The study strains exhibited two genetic constellations designed as Wa-like G1/G9/G12-P[6]/P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 and DS-1-like G2/G3/G12-P[4]/P[6]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Genotype G9P[4] strains possessed a DS-1-like genetic constellation with an E6 NSP4 gene, G9-P[4]-I2-R2-C2-M2-A2-N2-T2-E6-H2. The mixed genotype showed both Wa-like and DS-1-like profiles with a T6 NSP3 gene G1/G9P[8]/[4]-I1/I2-R1/R2-C1/C2-M1/M2-A1/A2-N1/N2-T1/T6-E1/E6-H1/H2. At the allelic level, the analysis of the Benin strains, reference strains (with known alleles), vaccine strains (with known alleles) identified 2-13 and 1-17 alleles for DS-1-like and Wa-like strains, respectively. Most of the study strains clustered into previously defined alleles, but we defined 3 new alleles for the VP7 (G3=1 new allele and G12=2 new alleles) and VP4 (P[4]=1 new allele and P[6]=2 new alleles) genes which formed the basis of the VP7 and VP4 gene clusters, respectively. For the remaining 9 genes, 0-6 new alleles were identified for both Wa-like and DS-1-like strains. This analysis of whole genome sequences of RVA strains circulating in Benin described genetic point mutations and reassortment events as well as novel alleles. Further detailed studies on these new alleles are needed and these data can also provide a baseline for studies on RVA in the post-vaccination period.