Whole Genome Analysis of Human Rotaviruses Reveals Single Gene Reassortant Rotavirus Strains in Zambia

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.

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.

Detection and full-genotype determination of rare and reassortant rotavirus A strains in Nizhny Novgorod in the European part of Russia

Reassortant DS-1-like rotavirus A strains have been shown to circulate widely in many countries around the world. In Russia, the prevalence of such strains remains unclear due to the preferred use of the traditional binary classification system. In this work, we obtained partial sequence data from all 11 genome segments and determined the full-genotype constellations of rare and reassortant rotaviruses circulating in Nizhny Novgorod in 2016-2019. DS-1-like G3P[8] and G8P[8] strains were found, reflecting the global trend. Most likely, these strains were introduced into the territory of Russia from other countries but subsequently underwent further evolutionary changes locally. G3P[8], G9P[8], and G12P[8] Wa-like strains of subgenotypic lineages that are unusual for the territory of Russia were also identified. Reassortant G2P[8], G4P[4], and G9P[4] strains with one Wa-like gene (VP4 or VP7) on a DS-1-like backbone were found, and these apparently had a local origin. Feline-like G3P[9] and G6P[9] strains were found to be phylogenetically close to BA222 isolated from a cat in Italy but carried some traces of reassortment with human strains from Russia and other countries. Thus, full-genotype determination of rotavirus A strains in Nizhny Novgorod has clarified some questions related to their origin and evolution.

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.

The Evolution of Post-Vaccine G8P[4] Group a Rotavirus Strains in Rwanda; Notable Variance at the Neutralization Epitope Sites

Africa has a high level of genetic diversity of rotavirus strains, which is suggested to be a possible reason contributing to the suboptimal effectiveness of rotavirus vaccines in this region. One strain that contributes to this rotavirus diversity in Africa is the G8P[4]. This study aimed to elucidate the entire genome and evolution of Rwandan G8P[4] strains. Illumina sequencing was performed for twenty-one Rwandan G8P[4] rotavirus strains. Twenty of the Rwandan G8P[4] strains had a pure DS-1-like genotype constellation, and one strain had a reassortant genotype constellation. Notable radical amino acid differences were observed at the neutralization sites when compared with cognate regions in vaccine strains potentially playing a role in neutralization escape. Phylogenetic analysis revealed that the closest relationship was with East African human group A rotavirus (RVA) strains for five of the genome segments. Two genome sequences of the NSP4 genome segment were closely related to bovine members of the DS-1-like family. Fourteen VP1 and eleven VP3 sequences had the closest relationships with the RotaTeq™ vaccine WC3 bovine genes. These findings suggest that the evolution of VP1 and VP3 might have resulted from reassortment events with RotaTeq™ vaccine WC3 bovine genes. The close phylogenetic relationship with East African G8P[4] strains from Kenya and Uganda suggests co-circulation in these countries. These findings highlight the need for continued whole-genomic surveillance to elucidate the evolution of G8P[4] strains, especially after the introduction of rotavirus vaccination.

Genomic Analysis of G2P[4] Group A Rotaviruses in Zambia Reveals Positive Selection in Amino Acid Site 7 of Viral Protein 3

The G2P[4] genotype is among the rotavirus strains that circulate commonly in humans. Several countries have reported its immediate upsurge after the introduction of rotavirus vaccination, raising concern about sub-optimal vaccine effectiveness against this genotype in the long term. This study aimed to gain insight into the evolution of post-vaccine Zambian G2P[4] group A rotavirus (RVA) strains and their overall genetic make-up by analysis of sequence alignments at the amino acid (AA) level. Twenty-nine Zambian G2P[4] rotavirus strains were subjected to whole-genome sequencing using the Illumina MiSeq^® platform. All the strains exhibited the typical DS-1-like genotype constellation, and the nucleotide sequences of the 11 genome segments showed high nucleotide similarities (>97%). Phylogenetic analyses together with representative global G2P[4] RVA showed that Zambian strains clustered into human lineages IV (for VP2, VP4, VP7, NSP1, and NSP5), V (for VP1, VP3, VP6, NSP2, and NSP3), and XXIII (for NSP4). The AA differences between the lineages where the study strains clustered and lineages of global reference strains were identified and analyzed. Selection pressure analysis revealed that AA site seven in the Viral Protein 3 (VP3) genome segment was under positive selection. This site occurs in the region of intrinsic disorder in the VP3 protein, and Zambian G2P[4] strains could potentially be utilizing this intrinsically disordered region to survive immune pressure. The Zambian G2P[4] strains from 2012 to 2016 comprised the G2P[4] strains that have been circulating globally since the early 2000s, highlighting the epidemiological fitness of these contemporary G2P[4] strains. Continuous whole-genome surveillance of G2P[4] strains remains imperative to understand their evolution during the post-vaccination period.

Isolation and Characterization of Distinct Rotavirus A in Bat and Rodent Hosts

Rotavirus A (RVA) causes diarrheal disease in humans and various animals. Recent studies have identified bat and rodent RVAs with evidence of zoonotic transmission and genome reassortment. However, the virological properties of bat and rodent RVAs with currently identified genotypes still need to be better clarified. Here, we performed virus isolation-based screening for RVA in animal specimens and isolated RVAs (representative strains: 16-06 and MpR12) from Egyptian fruit bat and Natal multimammate mouse collected in Zambia. Whole-genome sequencing and phylogenetic analysis revealed that the genotypes of bat RVA 16-06 were identical to that of RVA BATp39 strain from the Kenyan fruit bat, which has not yet been characterized. Moreover, all segments of rodent RVA MpR12 were highly divergent and assigned to novel genotypes, but RVA MpR12 was phylogenetically closer to bat RVAs than to other rodent RVAs, indicating a unique evolutionary history. We further investigated the virological properties of the isolated RVAs. In brief, we found that 16-06 entered cells by binding to sialic acids on the cell surface, while MpR12 entered in a sialic acid-independent manner. Experimental inoculation of suckling mice with 16-06 and MpR12 revealed that these RVAs are causative agents of diarrhea. Moreover, 16-06 and MpR12 demonstrated an ability to infect and replicate in a 3D-reconstructed primary human intestinal epithelium with comparable efficiency to the human RVA. Taken together, our results detail the unique genetic and virological features of bat and rodent RVAs and demonstrate the need for further investigation of their zoonotic potential. IMPORTANCE Recent advances in nucleotide sequence detection methods have enabled the detection of RVA genomes from various animals. These studies have discovered multiple divergent RVAs and have resulted in proposals for the genetic classification of novel genotypes. However, most of these RVAs have been identified via dsRNA viral genomes and not from infectious viruses, and their virological properties, such as cell/host tropisms, transmissibility, and pathogenicity, are unclear and remain to be clarified. Here, we successfully isolated RVAs with novel genome constellations from three bats and one rodent in Zambia. In addition to whole-genome sequencing, the isolated RVAs were characterized by glycan-binding affinity, pathogenicity in mice, and infectivity to the human gut using a 3D culture of primary intestinal epithelium. Our study reveals the first virological properties of bat and rodent RVAs with high genetic diversity and unique evolutional history and provides basic knowledge to begin estimating the potential of zoonotic transmission.

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.