Genomic Analysis of Rwandan G9P[8] Rotavirus Strains Pre- and Post-RotaTeq® Vaccine Reveals Significant Distinct Sub-Clustering in a Post-Vaccination Cohort.

Although the introduction of rotavirus vaccines has substantially contributed to the reduction in rotavirus morbidity and mortality, concerns persist about the re-emergence of variant strains that might alter vaccine effectiveness in the long term. The G9 strains re-emerged in Africa during the mid-1990s and have more recently become predominant in some countries, such as Ghana and Zambia. In Rwanda, during the 2011 to 2015 routine surveillance period, G9P[8] persisted during both the pre- and post-vaccine periods. The pre-vaccination cohort was based on the surveillance period of 2011 to 2012, and the post-vaccination cohort was based on the period of 2013 to 2015, excluding 2014. The RotaTeq® vaccine that was first introduced in Rwanda in 2012 is genotypically heterologous to Viral Protein 7 (VP7) G9. This study elucidated the whole genome of Rwandan G9P[8] rotavirus strains pre- and post-RotaTeq® vaccine introduction. Fecal samples from Rwandan children under the age of five years (pre-vaccine n = 23; post-vaccine n = 7), conventionally genotyped and identified as G9P[8], were included. Whole-genome sequencing was then performed using the Illumina® MiSeq platform. Phylogenetic analysis and pair-wise sequence analysis were performed using MEGA6 software. Distinct clustering of three post-vaccination study strains was observed in all 11 gene segments, compared to the other Rwandan G9P[8] study strains. Specific amino acid differences were identified across the gene segments of these three 2015 post-vaccine strains. Important amino acid differences were identified at position N242S in the VP7 genome segment of the three post-vaccine G9 strains compared to the other G9 strains. This substitution occurs at a neutralization epitope site and may slightly affect protein interaction at that position. These findings indicate that the Rwandan G9P[8] strains revealed a distinct sub-clustering pattern among post-vaccination study strains circulating in Rwanda, with changes at neutralization epitopes, which may play a role in neutralization escape from vaccine candidates. This emphasizes the need for continuous whole-genome surveillance to better understand the evolution and epidemiology of the G9P[8] strains post-vaccination.

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.

Evolutionary changes between pre- and post-vaccine South African group A G2P[4] rotavirus strains, 2003-2017.

The transient upsurge of G2P[4] group A rotavirus (RVA) after Rotarix vaccine introduction in several countries has been a matter of concern. To gain insight into the diversity and evolution of G2P[4] strains in South Africa pre- and post-RVA vaccination introduction, whole-genome sequencing was performed for RVA positive faecal specimens collected between 2003 and 2017 and samples previously sequenced were obtained from GenBank (n=103; 56 pre- and 47 post-vaccine). Pre-vaccine G2 sequences predominantly clustered within sub-lineage IVa-1. In contrast, post-vaccine G2 sequences clustered mainly within sub-lineage IVa-3, whereby a radical amino acid (AA) substitution, S15F, was observed between the two sub-lineages. Pre-vaccine P[4] sequences predominantly segregated within sub-lineage IVa while post-vaccine sequences clustered mostly within sub-lineage IVb, with a radical AA substitution R162G. Both S15F and R162G occurred outside recognised antigenic sites. The AA residue at position 15 is found within the signal sequence domain of Viral Protein 7 (VP7) involved in translocation of VP7 into endoplasmic reticulum during infection process. The 162 AA residue lies within the hemagglutination domain of Viral Protein 4 (VP4) engaged in interaction with sialic acid-containing structure during attachment to the target cell. Free energy change analysis on VP7 indicated accumulation of stable point mutations in both antigenic and non-antigenic regions. The segregation of South African G2P[4] strains into pre- and post-vaccination sub-lineages is likely due to erstwhile hypothesized stepwise lineage/sub-lineage evolution of G2P[4] strains rather than RVA vaccine introduction. Our findings reinforce the need for continuous whole-genome RVA surveillance and investigation of contribution of AA substitutions in understanding the dynamic G2P[4] epidemiology.

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

Rotarix® vaccine was implemented nationwide in Zambia in 2013. In this study, four unusual strains collected in the post-vaccine period were subjected to whole genome sequencing and analysis. The four strains possessed atypical genotype constellations, with at least one reassortant genome segment within the constellation. One of the strains (UFS-NGS-MRC-DPRU4749) was genetically and phylogenetically distinct in the VP4 and VP1 gene segments. Pairwise analyses demonstrated several amino acid disparities in the VP4 antigenic sites of this strain compared to that of Rotarix®. Although the impact of these amino acid disparities remains to be determined, this study adds to our understanding of the whole genomes of reassortant strains circulating in Zambia following Rotarix® vaccine introduction.

Metagenomic Analysis of the Enteric RNA Virome of Infants from the Oukasie Clinic, North West Province, South Africa, Reveals Diverse Eukaryotic Viruses.

Establishing a diverse gut microbiota after birth is essential for preventing illnesses later in life. However, little knowledge exists about the total viral population (virome) present in the gut of infants during the early developmental stage, with RNA viruses being generally overlooked. Therefore, this small pilot longitudinal study investigated the diversity and changes in the enteric RNA virome in healthy infants from South Africa. Faecal samples (n = 12) were collected from four infants at three time points (on average at 8, 13, and 25 weeks), and then sequenced on an Illumina MiSeq platform. The genomic analysis revealed a diverse population of human enteric viruses from the infants' stools, and changes in the enteric virome composition were observed over time. The Reoviridae family, more specifically the Rotavirus genus, was the most common and could be linked to viral shedding due to the administration of live-attenuated oral vaccines in South Africa, followed by the Picornaviridae family including parechoviruses, echoviruses, coxsackieviruses, enteroviruses, and polioviruses. Polioviruses were also linked to vaccine-related shedding. Astroviridae (astroviruses) and Caliciviridae (noroviruses) were present at low abundance. It is evident that an infant's gut is colonized by distinct viral populations irrespective of their health state. Further characterization of the human virome (with a larger participant pool) is imperative to provide more conclusive insights into the viral community structure and diversity that has been shown in the current study, despite the smaller sample size.

Whole Genome In-Silico Analysis of South African G1P[8] Rotavirus Strains Before and After Vaccine Introduction Over A Period of 14 Years.

Rotavirus G1P[8] strains account for more than half of the group A rotavirus (RVA) infections in children under five years of age, globally. A total of 103 stool samples previously characterized as G1P[8] and collected seven years before and seven years after introducing the Rotarix® vaccine in South Africa were processed for whole-genome sequencing. All the strains analyzed had a Wa-like constellation (G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1). South African pre- and post-vaccine G1 strains were clustered in G1 lineage-I and II while the majority (84.2%) of the P[8] strains were grouped in P[8] lineage-III. Several amino acid sites across ten gene segments with the exception of VP7 were under positive selective pressure. Except for the N147D substitution in the antigenic site of eight post-vaccine G1 strains when compared to both Rotarix® and pre-vaccine strains, most of the amino acid substitutions in the antigenic regions of post-vaccine G1P[8] strains were already present during the pre-vaccine period. Therefore, Rotarix® did not appear to have an impact on the amino acid differences in the antigenic regions of South African post-vaccine G1P[8] strains. However, continued whole-genome surveillance of RVA strains to decipher genetic changes in the post-vaccine period remains imperative.

Whole genome and in-silico analyses of G1P[8] rotavirus strains from pre- and post-vaccination periods in Rwanda.

Rwanda was the first low-income African country to introduce RotaTeq vaccine into its Expanded Programme on Immunization in May 2012. To gain insights into the overall genetic make-up and evolution of Rwandan G1P[8] strains pre- and post-vaccine introduction, rotavirus positive fecal samples collected between 2011 and 2016 from children under the age of 5 years as part of ongoing surveillance were genotyped with conventional RT-PCR based methods and whole genome sequenced using the Illumina MiSeq platform. From a pool of samples sequenced (n = 158), 36 were identified as G1P[8] strains (10 pre-vaccine and 26 post-vaccine), of which 35 exhibited a typical Wa-like genome constellation. However, one post vaccine strain, RVA/Human-wt/RWA/UFS-NGS:MRC-DPRU442/2012/G1P[8], exhibited a RotaTeq vaccine strain constellation of G1-P[8]-I2-R2-C2-M2-A3-N2-T6-E2-H3, with most of the gene segments having a close relationship with a vaccine derived reassortant strain, previously reported in USA in 2010 and Australia in 2012. The study strains segregated into two lineages, each containing a paraphyletic pre- and post-vaccine introduction sub-lineages. In addition, the study strains demonstrated close relationship amongst each other when compared with globally selected group A rotavirus (RVA) G1P[8] reference strains. For VP7 neutralization epitopes, amino acid substitutions observed at positions T91A/V, S195D and M217T in relation to the RotaTeq vaccine were radical in nature and resulted in a change in polarity from a polar to non-polar molecule, while for the VP4, amino acid differences at position D195G was radical in nature and resulted in a change in polarity from a polar to non-polar molecule. The polarity change at position T91A/V of the neutralizing antigens might play a role in generating vaccine-escape mutants, while substitutions at positions S195D and M217T may be due to natural fluctuation of the RVA. Surveillance of RVA at whole genome level will enhance further assessment of vaccine impact on circulating strains, the frequency of reassortment events under natural conditions and epidemiological fitness generated by such events.

Molecular Characterisation of a Rare Reassortant Porcine-Like G5P[6] Rotavirus Strain Detected in an Unvaccinated Child in Kasama, Zambia.

A human-porcine reassortant strain, RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6], was identified in a sample collected in 2014 from an unvaccinated 12 month old male hospitalised for gastroenteritis in Zambia. We sequenced and characterised the complete genome of this strain which presented the constellation: G5-P[6]-I1-R1-C1-M1-A8-N1-T1-E1-H1. The genotype A8 is often observed in porcine strains. Phylogenetic analyses showed that VP6, VP7, NSP2, NSP4, and NSP5 genes were closely related to cognate gene sequences of porcine strains (e.g., RVA/Pig-wt/CHN/DZ-2/2013/G5P[X] for VP7) from the NCBI database, while VP1, VP3, VP4, and NSP3 were closely related to porcine-like human strains (e.g., RVA/Human-wt/CHN/E931/2008/G4P[6] for VP1, and VP3). On the other hand, the origin of the VP2 was not clear from our analyses, as it was not only close to both porcine (e.g., RVA/Pig-tc/CHN/SWU-1C/2018/G9P[13]) and porcine-like human strains (e.g., RVA/Human-wt/LKA/R1207/2009/G4P[6]) but also to three human strains (e.g., RVA/Human-wt/USA/1476/1974/G1P[8]). The VP7 gene was located in lineage II that comprised only porcine strains, which suggests the occurrence of independent porcine-to-human reassortment events. The study strain may have collectively been derived through interspecies transmission, or through reassortment event(s) involving strains of porcine and porcine-like human origin. The results of this study underline the importance of whole-genome characterisation of rotavirus strains and provide insights into interspecies transmissions from porcine to humans.

Uncovering the First Atypical DS-1-like G1P[8] Rotavirus Strains That Circulated during Pre-Rotavirus Vaccine Introduction Era in South Africa.

Emergence of DS-1-like G1P[8] group A rotavirus (RVA) strains during post-rotavirus vaccination period has recently been reported in several countries. This study demonstrates, for the first time, rare atypical DS-1-like G1P[8] RVA strains that circulated in 2008 during pre-vaccine era in South Africa. Rotavirus positive samples were subjected to whole-genome sequencing. Two G1P[8] strains (RVA/Human-wt/ZAF/UFS-NGS-MRC-DPRU1971/2008/G1P[8] and RVA/Human-wt/ZAF/UFS-NGS-MRC-DPRU1973/2008/G1P[8]) possessed a DS-1-like genome constellation background (I2-R2-C2-M2-A2-N2-T2-E2-H2). The outer VP4 and VP7 capsid genes of the two South African G1P[8] strains had the highest nucleotide (amino acid) nt (aa) identities of 99.6-99.9% (99.1-100%) with the VP4 and the VP7 genes of a locally circulating South African strain, RVA/Human-wt/ZAF/MRC-DPRU1039/2008/G1P[8]. All the internal backbone genes (VP1-VP3, VP6, and NSP1-NSP5) had the highest nt (aa) identities with cognate internal genes of another locally circulating South African strain, RVA/Human-wt/ZAF/MRC-DPRU2344/2008/G2P[6]. The two study strains emerged through reassortment mechanism involving locally circulating South African strains, as they were distinctly unrelated to other reported atypical G1P[8] strains. The identification of these G1P[8] double-gene reassortants during the pre-vaccination period strongly supports natural RVA evolutionary mechanisms of the RVA genome. There is a need to maintain long-term whole-genome surveillance to monitor such atypical strains.

Positive interactions between nitrogen-fixing legumes and four different neighbouring species in a biodiversity experiment.

The importance of facilitative processes due to the presence of nitrogen-fixing legumes in temperate grasslands is a contentious issue in biodiversity experiments. Despite a multitude of studies of fertilization effects of legumes on associated nonfixers in agricultural systems, we know little about the dynamics in more diverse systems. We hypothesised that the identity of target plant species (phytometers) and the diversity of neighbouring plant species would affect the magnitude of such positive species interactions. We therefore sampled aboveground tissues of phytometers planted into all plots of a grassland biodiversity-ecosystem functioning experiment and analysed their N concentrations, delta15N values and biomasses. The four phytometer species (Festuca pratensis, Plantago lanceolata, Knautia arvensis and Trifolium pratensis) each belonged to one of the four plant functional groups used in the experiment and allowed the effects of diversity on N dynamics in individual species to be assessed. We found significantly lower delta15N values and higher N concentrations and N contents (amount of N per plant) in phytometer species growing with legumes, indicating a facilitative role for legumes in these grassland ecosystems. Our data suggest that the main driving force behind these facilitative interactions in plots containing legumes was reduced competition for soil nitrate ("nitrate sparing"), with apparent N transfer playing a secondary role. Interestingly, species richness (and to a lesser extent functional group number) significantly decreased delta15N values, N concentrations and N content irrespective of any legume effect. Possible mechanisms behind this effect, such as increased N mineralisation and nitrate uptake in more diverse plots, now need further investigation. The magnitude of the positive interactions depended on the identity of the phytometer species. Evidence for increased N uptake in communities containing legumes was found in all three nonlegume phytometer species, with a subsequent strong increase in biomass in the grass F. pratensis across all diversity levels, and a lesser biomass gain in P. lanceolata and K. arvensis. In contrast, the legume phytometer species T. pratense was negatively affected when other legumes were present in their host communities across all diversity levels.

Effects of plant diversity on invertebrate herbivory in experimental grassland.

The rate at which a plant species is attacked by invertebrate herbivores has been hypothesized to depend on plant species richness, yet empirical evidence is scarce. Current theory predicts higher herbivore damage in monocultures than in species-rich mixtures. We quantified herbivore damage by insects and molluscs to plants in experimental plots established in 2002 from a species pool of 60 species of Central European Arrhenatherum grasslands. Plots differed in plant species richness (1, 2, 4, 8, 16, 60 species), number of functional groups (1, 2, 3, 4), functional group and species composition. We estimated herbivore damage by insects and molluscs at the level of transplanted plant individuals ("phytometer" species Plantago lanceolata, Trifolium pratense, Rumex acetosa) and of the entire plant community during 2003 and 2004. In contrast to previous studies, our design allows specific predictions about the relative contributions of functional diversity, plant functional identity, and species richness in relation to herbivory. Additionally, the phytometer approach is new to biodiversity-herbivory studies, allowing estimates of species-specific herbivory rates within the larger biodiversity-ecosystem functioning context. Herbivory in phytometers and experimental communities tended to increase with plant species richness and the number of plant functional groups, but the effects were rarely significant. Herbivory in phytometers was in some cases positively correlated with community biomass or leaf area index. The most important factor influencing invertebrate herbivory was the presence of particular plant functional groups. Legume (grass) presence strongly increased (decreased) herbivory at the community level. The opposite pattern was found for herbivory in T. pratense phytometers. We conclude that (1) plant species richness is much less important than previously thought and (2) plant functional identity is a much better predictor of invertebrate herbivory in temperate grassland ecosystems.