Multiple Introductions and Predominance of Rotavirus Group A Genotype G3P[8] in Kilifi, Coastal Kenya, 4 Years after Nationwide Vaccine Introduction

Generation of Recombinant Authentic Live Attenuated Human Rotavirus Vaccine Strain RIX4414 (Rotarix®) from Cloned cDNAs Using Reverse Genetics

The live attenuated human rotavirus vaccine strain RIX4414 (Rotarix®) is used worldwide to prevent severe rotavirus-induced diarrhea in infants. This strain was attenuated through the cell culture passaging of its predecessor, human strain 89-12, which resulted in multiple genomic mutations. However, the specific molecular reasons underlying its attenuation have remained elusive, primarily due to the absence of a suitable reverse genetics system enabling precise genetic manipulations. Therefore, we first completed the sequencing of its genome and then developed a reverse genetics system for the authentic RIX4414 virus. Our experimental results demonstrate that the rescued recombinant RIX4414 virus exhibits biological characteristics similar to those of the parental RIX4414 virus, both in vitro and in vivo. This novel reverse genetics system provides a powerful tool for investigating the molecular basis of RIX4414 attenuation and may facilitate the rational design of safer and more effective human rotavirus vaccines.

Genomic characterization of the rotavirus G3P[8] strain in vaccinated children, reveals possible reassortment events between human and animal strains in Manhiça District, Mozambique

Mozambique introduced the rotavirus vaccine (Rotarix®; GlaxoSmithKline Biologicals, Rixensart, Belgium) in 2015, and since then, the Centro de Investigação em Saúde de Manhiça has been monitoring its impact on rotavirus-associated diarrhea and the trend of circulating strains, where G3P[8] was reported as the predominant strain after the vaccine introduction. Genotype G3 is among the most commonly detected Rotavirus strains in humans and animals, and herein, we report on the whole genome constellation of G3P[8] detected in two children (aged 18 months old) hospitalized with moderate-to-severe diarrhea at the Manhiça District Hospital. The two strains had a typical Wa-like genome constellation (I1-R1-C1-M1-A1-N1-T1-E1-H1) and shared 100% nucleotide (nt) and amino acid (aa) identities in 10 gene segments, except for VP6. Phylogenetic analysis demonstrated that genome segments encoding VP7, VP6, VP1, NSP3, and NSP4 of the two strains clustered most closely with porcine, bovine, and equine strains with identities ranging from 86.9-99.9% nt and 97.2-100% aa. Moreover, they consistently formed distinct clusters with some G1P[8], G3P[8], G9P[8], G12P[6], and G12P[8] strains circulating from 2012 to 2019 in Africa (Mozambique, Kenya, Rwanda, and Malawi) and Asia (Japan, China, and India) in genome segments encoding six proteins (VP2, VP3, NSP1-NSP2, NSP5/6). The identification of segments exhibiting the closest relationships with animal strains shows significant diversity of rotavirus and suggests the possible occurrence of reassortment events between human and animal strains. This demonstrates the importance of applying next-generation sequencing to monitor and understand the evolutionary changes of strains and evaluate the impact of vaccines on strain diversity.

Comparative whole genome analysis reveals re-emergence of human Wa-like and DS-1-like G3 rotaviruses after Rotarix vaccine introduction in Malawi

G3 rotaviruses rank among the most common rotavirus strains worldwide in humans and animals. However, despite a robust long-term rotavirus surveillance system from 1997 at Queen Elizabeth Central Hospital in Blantyre, Malawi, these strains were only detected from 1997 to 1999 and then disappeared and re-emerged in 2017, 5 years after the introduction of the Rotarix rotavirus vaccine. Here, we analysed representative twenty-seven whole genome sequences (G3P[4], n = 20; G3P[6], n = 1; and G3P[8], n = 6) randomly selected each month between November 2017 and August 2019 to understand how G3 strains re-emerged in Malawi. We found four genotype constellations that were associated with the emergent G3 strains and co-circulated in Malawi post-Rotarix vaccine introduction: G3P[4] and G3P[6] strains with the DS-1-like genetic backbone genes (G3-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2 and G3-P[6]-I2-R2-C2-M2-A2-N2-T2-E2-H2), G3P[8] strains with the Wa-like genetic backbone genes (G3-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1), and reassortant G3P[4] strains consisting of the DS-1-like genetic backbone genes and a Wa-like NSP2 (N1) gene (G3-P[4]-I2-R2-C2-M2-A2-N1-T2-E2-H2). Time-resolved phylogenetic trees demonstrated that the most recent common ancestor for each ribonucleic acid (RNA) segment of the emergent G3 strains was between 1996 and 2012, possibly through introductions from outside the country due to the limited genetic similarity with G3 strains which circulated before their disappearance in the late 1990s. Further genomic analysis revealed that the reassortant DS-1-like G3P[4] strains acquired a Wa-like NSP2 genome segment (N1 genotype) through intergenogroup reassortment; an artiodactyl-like VP3 through intergenogroup interspecies reassortment; and VP6, NSP1, and NSP4 segments through intragenogroup reassortment likely before importation into Malawi. Additionally, the emergent G3 strains contain amino acid substitutions within the antigenic regions of the VP4 proteins which could potentially impact the binding of rotavirus vaccine-induced antibodies. Altogether, our findings show that multiple strains with either Wa-like or DS-1-like genotype constellations have driven the re-emergence of G3 strains. The findings also highlight the role of human mobility and genome reassortment events in the cross-border dissemination and evolution of rotavirus strains in Malawi necessitating the need for long-term genomic surveillance of rotavirus in high disease-burden settings to inform disease prevention and control.

Genetic diversity and Genomic analysis of G3P[6] and equine-like G3P[8] in Children Under-five from Southern Highlands and Eastern Tanzania

Rotavirus group A genomic characterization is critical for understanding the mechanisms of rotavirus diversity, such as reassortment events and possible interspecies transmission. However, little is known about the genetic diversity and genomic relationship of the rotavirus group A strains circulating in Tanzania. The genetic and genomic relationship of RVA genotypes was investigated in children under the age of five. A total of 169 faecal samples were collected from under-five with diarrhea in Mbeya, Iringa and Morogoro regions of Tanzania. The RVA were screened in children under five with diarrhea using reverse transcription PCR for VP7 and VP4, and the G and P genotypes were determined using Sanger dideoxynucleotide cycle sequencing. Whole-genome sequencing was performed on selected genotypes. The overall RVA rate was 4.7% (8/169). The G genotypes were G3 (7/8) and G6 (1/8) among the 8 RVA positives, while the P genotypes were P[6] (4/8) and P[8] (2), and the other two were untypeable. G3P[6] and G3P[8] were the identified genotype combinations. The genomic analysis reveals that the circulating G3P[8] and G3P[6] isolates from children under the age of five with diarrhea had a DS-1-like genome configuration (I2-R2-C2-M2-Ax-N2-T2-E2-H2). The phylogenic analysis revealed that all 11 segments of G3P[6] were closely related to human G3P[6] identified in neighboring countries such as Uganda, Kenya, and other African countries, implying that G3P[6] strains descended from a common ancestor. Whereas, G3P[8] were closely related to previously identified equine-like G3P[P8] from Kenya, Japan, Thailand, Brazil, and Taiwan, implying that this strain was introduced rather than reassortment events. We discovered amino acid differences at antigenic epitopes and N-linked glycosylation sites between the wild type G3 and P[8] compared to vaccine strains, implying that further research into the impact of these differences on vaccine effectiveness is warranted. The phylogenic analysis of VP7 also identified a bovine-like G6. For the first time in Tanzania, we report the emergence of novel equine-like G3 and bovine-like G6 RVA strains, highlighting the importance of rotavirus genotype monitoring and genomic analysis of representative genotypes.

Molecular Epidemiology of Rotavirus Strains in Symptomatic and Asymptomatic Children in Manhiça District, Southern Mozambique 2008-2019

Group A rotaviruses remain the leading cause of diarrhoea in children aged <5 years. Mozambique introduced rotavirus vaccine (Rotarix^®) in September 2015. We report rotavirus genotypes circulating among symptomatic and asymptomatic children in Manhiça District, Mozambique, pre- and post-vaccine introduction. Stool was collected from enrolled children and screened for rotavirus by enzyme-immuno-sorbent assay. Positive specimens were genotyped for VP7 (G genotypes) and VP4 (P genotypes) by the conventional reverse transcriptase polymerase chain reaction. The combination G12P[8] was more frequently observed in pre-vaccine than in post-vaccine introduction, in moderate to severe diarrhoea (34%, 61/177 vs. 0, p < 0.0001) and controls (23%, 26/113 vs. 0, p = 0.0013) and mixed genotypes (36%, 24/67 vs. 7% 4/58, p = 0.0003) in less severe diarrhoea. We observed changes in post-vaccine compared to pre-vaccine introduction, where G3P[4] and G3P[8] were prevalent in moderate to severe diarrhoea (10%, 5/49 vs. 0, p = 0.0002; and 14%, 7/49 vs. 1%, 1/177, p < 0.0001; respectively), and in less severe diarrhoea (21%, 12/58 vs. 0, p = 0.003; and 24%, 14/58 vs. 0, p < 0.0001; respectively). Our surveillance demonstrated the circulation of similar genotypes contemporaneously among cases and controls, as well as switching from pre- to post-vaccine introduction. Continuous surveillance is needed to evaluate the dynamics of the changes in genotypes following vaccine introduction.