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

Abstract

Genotype 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 in Blantyre, Malawi, these strains were only detected from 1997 to 1999 and then disappeared and re-emerged in 2017, five years after the introduction of the Rotarix rotavirus vaccine. Here we analysed 27 whole genome sequences to understand how G3 strains re-emerged in Malawi. We randomly selected samples each month between November 2017 and August 2019 from stool samples of children hospitalised with acute diarrhoea at the Queen Elizabeth Hospital in Blantyre, Malawi. We found three genotypes namely G3P[4] (n=20), G3P[6] (n=1) and G3P[8] (n=6) associated with the re-emergence of G3 strains in Malawi post-Rotarix vaccine introduction. The identified genotypes co-circulated at different time points and were associated with three typical human G3 strains consisting of either a Wa-like or DS-1-like genetic constellation and reassortant strains possessing Wa-like and DS-1-like genetic backbones. Time-resolved phylogenetic trees demonstrated that the most recent common ancestor for each segment of the re-emerged G3 strains emerged 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 intragenogroup interspecies reassortment; and VP6, NSP1 and NSP4 segments through intragenogroup reassortment likely before importation into Malawi. Additionally, the re-emerged 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 shows that multiple rather than a single genotype have driven the re-emergence of G3 strains likely from other countries highlighting the role of human mobility and genome reassortment events in the 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

“…The 20 N2 genes of Rwandan G8P [4] strains that clustered into human lineage V had the highest level of similarity (99.6-100%) to a Ugandan strain (Table 4). The one strain with the N1 genotype had the closest nucleotide similarity (99.7%) to the NSP2 gene (N1 genotype) of a G9P [8] Ugandan strain (Figure 1H). Rwandan T2 gene sequences clustered into lineage V alongside other human T2 strains, with the highest nucleotide similarities (99.4-100%) observed with a Ugandan strain (Table 4; Figure 1I).…”

“…The G1, G3, G4, G9, and G12 genotypes, in combination with the P [8] genotype, are commonly linked with the Wa-like genogroup [7]. However, there have been reports of atypical strains among the G1P [8], G3P [8], and G9P [8] strains, bearing a DS-1-like genome constellation [8][9][10][11]. In contrast, the G2P [4] strains are typically linked with the DS-1-like genogroup [7].…”

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

“…The G1P[8], G3P[8], G9P[8], and G12P[8] typically belong to the Wa-like genogroup [ 9 ]. However, atypical DS-1-like G1P[8] [ 11 , 12 ], G3P[8] [ 13 , 14 ], and G9P[8] strains [ 15 ] have been reported. On the other hand, the G2P[4] strains usually belong to the DS-1-like genogroup [ 16 ].…”