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.

Draft genomes of Aeromonas caviae from patients with cholera-like illness during the 2022-2023 cholera outbreak in Malawi

Aeromonas caviae is an increasingly recognized etiological agent of acute gastroenteritis. Here, we report five draft genomes of A. caviae isolated from suspected cholera cases during the 2022-2023 cholera outbreak in Malawi.

Cheminformatics bioprospection of selected medicinal plants metabolites against trypsin cleaved VP4 (spike protein) of rotavirus A

Rotaviruses have continued to be the primary cause of acute dehydrating diarrhoea in children under five years of age despite the global introduction of four World Health Organization (WHO) prequalified oral vaccines in over 106 countries. Currently, no medication is approved by the Food and Drug Administration (FDA) specifically for treating rotavirus A-induced diarrhoea. Consequently, it is important to focus on developing prophylactic and curative therapeutics to combat rotaviral infections. For the first time, this study computationally screened and identified metabolites from Spondias mombin, Macaranga barteri and Dicerocaryum eriocarpum as potential novel inhibitors with broad-spectrum activity against VP5* and VP8* (spike protein) of rotavirus A (RVA). The initial top 20 metabolites identified through molecular docking were further filtered using drug-likeness and pharmacokinetics parameters. The molecular properties of the resulting top-ranked compounds were predicted by conducting density functional theory (DFT) calculations, while molecular dynamics (MD) simulation revealed their thermodynamic compatibility with a significant affinity towards VP8* than VP5*. Except for ellagic acid (-11.78 kcal/mol), the lead compounds had higher binding free energy than the reference standard (VP5* (-11.81 kcal/mol), VP8* (-14.12 kcal/mol)) with 2SG (-20.98 kcal/mol) and apigenin-4'-glucoside (-23.56 kcal/mol) having the highest affinity towards VP5* and VP8*, respectively. Of all the top-ranked compounds, better broad-spectrum affinities for both VP5* and VP8* than tizoxanide were observed in 2SG (VP5* (-20.98 kcal/mol), VP8* (-20.13 kcal/mol)) and sericetin (VP5* (-20.46 kcal/mol), VP8* (-18.31 kcal/mol)). While the identified leads could be regarded as potential modulators of the investigated therapeutic targets for effective management of rotaviral infection, additional in vitro and in vivo evaluation is strongly recommended, and efforts are on-going in this regard.Communicated by Ramaswamy H. Sarma.

Metabolites profiling and cheminformatics bioprospection of selected medicinal plants against the main protease and RNA-dependent RNA polymerase of SARS-CoV-2

Despite the existence of some vaccines, SARS-CoV-2 (S-2) infections persist for various reasons relating to vaccine reluctance, rapid mutation rate, and an absence of specific treatments targeted to the infection. Due to their availability, low cost and low toxicity, research into potentially repurposing phytometabolites as therapeutic alternatives has gained attention. Therefore, this study explored the antiviral potential of metabolites of some medicinal plants [Spondias mombin, Macaranga barteri and Dicerocaryum eriocarpum (Sesame plant)] identified using liquid chromatography-mass spectrometry (LCMS) as possible inhibitory agents against the S-2 main protease (S-2 MP) and RNA-dependent RNA polymerase (RP) using computational approaches. Molecular docking was used to identify the compounds with the best affinities for the selected therapeutics targets. Afterwards, compounds with poor physicochemical characteristics, pharmacokinetics, and drug-likeness were screened out. The top-ranked compounds were further subjected to a 120-ns molecular dynamics (MD) simulation. Only quercetin 3-O-rhamnoside (-48.77 kcal/mol) had higher binding free energy than the reference standard (zafirlukast) (-44.99 kcal/mol) against S-2 MP. Conversely, all the top-ranked compounds (ellagic acid hexoside, spiraeoside, apigenin-4'-glucoside and chrysoeriol 7-glucuronide) except gnetin L (-24.24 kcal/mol) had higher binding free energy (-55.19 kcal/mol, -52.75 kcal/mol, -47.22 kcal/mol and -43.35 kcal/mol) respectively, against S-2 RP relative to the reference standard (-34.79 kcal/mol). The MD simulations study further revealed that the investigated inhibitors are thermodynamically stable and form structurally compatible complexes that impede the regular operation of the respective S-2 therapeutic targets. Although, these S-2 therapeutic candidates are promising, further in vitro and in vivo evaluation is required and highly recommended.Communicated by Ramaswamy H. Sarma.

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.

Metagenomics characterization of respiratory viral RNA pathogens in children under five years with severe acute respiratory infection in the Free State, South Africa

Prompt detection of viral respiratory pathogens is crucial in managing respiratory infection including severe acute respiratory infection (SARI). Metagenomics next-generation sequencing (mNGS) and bioinformatics analyses remain reliable strategies for diagnostic and surveillance purposes. This study evaluated the diagnostic utility of mNGS using multiple analysis tools compared with multiplex real-time PCR for the detection of viral respiratory pathogens in children under 5 years with SARI. Nasopharyngeal swabs collected in viral transport media from 84 children admitted with SARI as per the World Health Organization definition between December 2020 and August 2021 in the Free State Province, South Africa, were used in this study. The obtained specimens were subjected to mNGS using the Illumina MiSeq system, and bioinformatics analysis was performed using three web-based analysis tools; Genome Detective, One Codex and Twist Respiratory Viral Research Panel. With average reads of 211323, mNGS detected viral pathogens in 82 (97.6%) of the 84 patients. Viral aetiologies were established in nine previously undetected/missed cases with an additional bacterial aetiology (Neisseria meningitidis) detected in one patient. Furthermore, mNGS enabled the much needed viral genotypic and subtype differentiation and provided significant information on bacterial co-infection despite enrichment for RNA viruses. Sequences of nonhuman viruses, bacteriophages, and endogenous retrovirus K113 (constituting the respiratory virome) were also uncovered. Notably, mNGS had lower detectability rate for severe acute respiratory syndrome coronavirus 2 (missing 18/32 cases). This study suggests that mNGS, combined with multiple/improved bioinformatics tools, is practically feasible for increased viral and bacterial pathogen detection in SARI, especially in cases where no aetiological agent could be identified by available traditional methods.

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.

Metagenomic Analysis of Respiratory RNA Virome of Children with and without Severe Acute Respiratory Infection from the Free State, South Africa during COVID-19 Pandemic Reveals Higher Diversity and Abundance in Summer Compared with Winter Period

Viral respiratory infections contribute to significant morbidity and mortality in children. Currently, there are limited reports on the composition and abundance of the normal commensal respiratory virome in comparison to those in severe acute respiratory infections (SARIs) state. This study characterised the respiratory RNA virome in children ≤ 5 years with (n = 149) and without (n = 139) SARI during the summer and winter of 2020/2021 seasons in South Africa. Nasopharyngeal swabs were, collected, pooled, enriched for viral RNA detection, sequenced using Illumina MiSeq, and analysed using the Genome Detective bioinformatic tool. Overall, Picornaviridae, Paramoxyviridae, Pneumoviridae, Picobirnaviridae, Totiviridae, and Retroviridae families were the most abundant viral population in both groups across both seasons. Human rhinovirus and endogenous retrovirus K113 were detected in most pools, with exclusive detection of Pneumoviridae in SARI pools. Generally, higher viral diversity/abundance was seen in children with SARI and in the summer pools. Several plant/animal viruses, eukaryotic viruses with unclear pathogenicity including a distinct rhinovirus A type, were detected. This study provides remarkable data on the respiratory RNA virome in children with and without SARI with a degree of heterogeneity of known viruses colonizing their respiratory tract. The implication of the detected viruses in the dynamics/progression of SARI requires further investigations.

Plant and Soil Core Mycobiomes in a Two-Year Sorghum–Legume Intercropping System of Underutilized Crops in South Africa

Fungal communities form close beneficial (mutualists) or detrimental (pathogens) associations with their plant hosts. Their diversity and abundance can be affected by agricultural practices which include cropping systems such as rotations and intercropping. Despite the importance of cropping systems in increasing productivity, knowledge of the fungal mycobiome and the core inhabitants for under-utilised cereal and legume crops, particularly over a period, is still limited. The core mycobiomes in plant tissues and bulk soils of a cereal–legume intercrop were characterized over two years using high-throughput sequencing. The intercropping trial consisted of sorghum, Bambara groundnut, cowpea, dry bean, and soybean. A greater number of molecular operational taxonomic units (MOTUs) were found in plant tissues compared to those from the soils and between year one and year two. Principal coordinate analyses revealed that fungal communities for each year were relatively distinct, particularly for the soils. The core mycobiome was dominated by a Davidiellaceae sp. (Cladosporium), Didymellaceae sp. 1 (Phoma), Didymellaceae sp. 2 (Epicoccum), Fusarium sp. 2, Unidentified (Ascomycota), and Cryptococcus MOTUs that were present in all plant tissues and soils of year one and two. Other key MOTUs were only specific to a year, substrate, or crop. Although the mycobiome of sorghum were more distinct than the cores of the legumes, there were still MOTUs dominant across all of the crops. Characterization of this baseline core across two years provides insight into those fungi that are always present in these crops, and that could be utilized in improving crop performance and productivity.

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance

Investment in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequencing in Africa over the past year has led to a major increase in the number of sequences that have been generated and used to track the pandemic on the continent, a number that now exceeds 100,000 genomes. Our results show an increase in the number of African countries that are able to sequence domestically and highlight that local sequencing enables faster turnaround times and more-regular routine surveillance. Despite limitations of low testing proportions, findings from this genomic surveillance study underscore the heterogeneous nature of the pandemic and illuminate the distinct dispersal dynamics of variants of concern-particularly Alpha, Beta, Delta, and Omicron-on the continent. Sustained investment for diagnostics and genomic surveillance in Africa is needed as the virus continues to evolve while the continent faces many emerging and reemerging infectious disease threats. These investments are crucial for pandemic preparedness and response and will serve the health of the continent well into the 21st century.

Pathogen Profile of Children Hospitalised with Severe Acute Respiratory Infections during COVID-19 Pandemic in the Free State Province, South Africa

Severe acute respiratory infections (SARI) contribute to mortality in children ≤5 years. Their microbiological aetiologies are often unknown and may be exacerbated in light of coronavirus disease 19 (COVID-19). This study reports on respiratory pathogens in children ≤5 years (n = 84) admitted with SARI during and between the second and third waves of COVID-19 infection in South Africa. Nasopharyngeal/oropharyngeal swabs collected were subjected to viral detection using QIAstat-Dx^® Respiratory SARS-CoV-2 Panel. The results revealed viral positivity and negativity detection rates of 88% (74/84) and 12% (10/84), respectively. Of the 21 targeted pathogens, human rhinovirus/enterovirus (30%), respiratory syncytial virus (RSV; 26%), and severe acute respiratory syndrome coronavirus 2 (24%) were mostly detected, with other viruses being 20% and a co-infection rate of 64.2% (54/84). Generally, RSV-positive samples had lower Ct values, and fewer viruses were detected during the third wave. Changes in the circulation patterns of respiratory viruses with total absence of influenza virus could be attributed to measures against COVID-19 transmission, which may result in waned immunity, thereby increasing susceptibility to severe infections in the following season. High viral co-infection rate, as detected, may complicate diagnosis. Nonetheless, accurate identification of the pathogens may guide treatment decisions and infection control.

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.

Poultry gut health - microbiome functions, environmental impacts, microbiome engineering and advancements in characterization technologies

The gastrointestinal tract (GIT) health impacts animal productivity. The poultry microbiome has functions which range from protection against pathogens and nutrients production, to host immune system maturation. Fluctuations in the microbiome have also been linked to prevailing environmental conditions. Healthy poultry birds possess a natural resistance to infection. However, the exploration of environmental impacts and other relevant factors on poultry growth and health have been underplayed. Since good performance and growth rate are central to animal production, the host-microbiome relationship remains integral. Prior to the emergence of metagenomic techniques, conventional methods for poultry microbiome studies were used and were low-throughput and associated with insufficient genomic data and high cost of sequencing. Fortunately, the advent of high-throughput sequencing platforms have circumvented some of these shortfalls and paved the way for increased studies on the poultry gut microbiome diversity and functions. Here, we give an up-to-date review on the impact of varied environments on microbiome profile, as well as microbiome engineering and microbiome technology advancements. It is hoped that this paper will provide invaluable information that could guide and inspire further studies on the lingering pertinent questions about the poultry microbiome.

A year of genomic surveillance reveals how the SARS-CoV-2 pandemic unfolded in Africa

The progression of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic in Africa has so far been heterogeneous, and the full impact is not yet well understood. In this study, we describe the genomic epidemiology using a dataset of 8746 genomes from 33 African countries and two overseas territories. We show that the epidemics in most countries were initiated by importations predominantly from Europe, which diminished after the early introduction of international travel restrictions. As the pandemic progressed, ongoing transmission in many countries and increasing mobility led to the emergence and spread within the continent of many variants of concern and interest, such as B.1.351, B.1.525, A.23.1, and C.1.1. Although distorted by low sampling numbers and blind spots, the findings highlight that Africa must not be left behind in the global pandemic response, otherwise it could become a source for new variants.

Genetic Characterisation of South African and Mozambican Bovine Rotaviruses Reveals a Typical Bovine-like Artiodactyl Constellation Derived through Multiple Reassortment Events

This study presents whole genomes of seven bovine rotavirus strains from South Africa and Mozambique. Double-stranded RNA, extracted from stool samples without prior adaptation to cell culture, was used to synthesise cDNA using a self-annealing anchor primer ligated to dsRNA and random hexamers. The cDNA was subsequently sequenced using an Illumina MiSeq platform without prior genome amplification. All strains exhibited bovine-like artiodactyl genome constellations (G10/G6-P[11]/P[5]-I2-R2-C2-M2-A3/A11/A13-N2-T6-E2-H3). Phylogenetic analysis revealed relatively homogenous strains, which were mostly related to other South African animal strains or to each other. It appears that these study strains represent a specific bovine rotavirus population endemic to Southern Africa that was derived through multiple reassortment events. While one Mozambican strain, MPT307, was similar to the South African strains, the second strain, MPT93, was divergent from the other study strains, exhibiting evidence of interspecies transmission of the VP1 and NSP2 genes. The data presented in this study not only contribute to the knowledge of circulating African bovine rotavirus strains, but also emphasise the need for expanded surveillance of animal rotaviruses in African countries in order to improve our understanding of rotavirus strain diversity.

Genetic characterization of G12P[6] and G12P[8] rotavirus strains collected in six African countries between 2010 and 2014

BACKGROUND

G12 rotaviruses were first observed in sub-Saharan Africa in 2004 and since then have continued to emerge and spread across the continent and are reported as a significant human rotavirus genotype in several African countries, both prior to and after rotavirus vaccine introduction. This study investigated the genetic variability of 15 G12 rotavirus strains associated with either P[6] or P[8] identified between 2010 and 2014 from Ethiopia, Kenya, Rwanda, Tanzania, Togo and Zambia.

METHODS

The investigation was carried out by comparing partial VP7 and partial VP4 sequences of the African G12P[6] and G12P[8] strains with the available GenBank sequences and exploring the recognized neutralization epitopes of these strains. Additionally, Bayesian evolutionary analysis was carried out using Markov Chain Monte Carlo (MCMC) implemented in BEAST to estimate the time to the most recent ancestor and evolutionary rate for these G12 rotavirus strains.

RESULTS

The findings suggested that the VP7 and VP4 nucleotide and amino acid sequences of the G12 strains circulating in African countries are closely related, irrespective of country of origin and year of detection, with the exception of the Ethiopian strains that clustered distinctly. Neutralization epitope analysis revealed that rotavirus VP4 P[8] genes associated with G12 had amino acid sequences similar to those reported globally including the vaccine strains in RotaTeq and Rotarix. The estimated evolutionary rate of the G12 strains was 1.016 × 10^- 3 substitutions/site/year and was comparable to what has been previously reported. Three sub-clusters formed within the current circulating lineage III shows the diversification of G12 from three independent ancestries within a similar time frame in the late 1990s.

CONCLUSIONS

At present it appears to be unlikely that widespread vaccine use has driven the molecular evolution and sustainability of G12 strains in Africa. Continuous monitoring of rotavirus genotypes is recommended to assess the long-term impact of rotavirus vaccination on the dynamic nature of rotavirus evolution on the continent.

Whole-Genome Analyses Identifies Multiple Reassortant Rotavirus Strains in Rwanda Post-Vaccine Introduction

Children in low-and middle-income countries, including Rwanda, experience a greater burden of rotavirus disease relative to developed countries. Evolutionary mechanisms leading to multiple reassortant rotavirus strains have been documented over time which influence the diversity and evolutionary dynamics of novel rotaviruses. Comprehensive rotavirus whole-genome analysis was conducted on 158 rotavirus group A (RVA) samples collected pre- and post-vaccine introduction in children less than five years in Rwanda. Of these RVA positive samples, five strains with the genotype constellations G4P[4]-I1-R2-C2-M2-A2-N2-T1-E1-H2 (n = 1), G9P[4]-I1-R2-C2-M2-A1-N1-T1-E1-H1 (n = 1), G12P[8]-I1-R2-C2-M1-A1-N2-T1-E2-H3 (n = 2) and G12P[8]-I1-R1-C1-M1-A2-N2-T2-E1-H1 (n = 1), with double and triple gene reassortant rotavirus strains were identified. Phylogenetic analysis revealed a close relationship between the Rwandan strains and cognate human RVA strains as well as the RotaTeq^® vaccine strains in the VP1, VP2, NSP2, NSP4 and NSP5 gene segments. Pairwise analyses revealed multiple differences in amino acid residues of the VP7 and VP4 antigenic regions of the RotaTeq^® vaccine strain and representative Rwandan study strains. Although the impact of such amino acid changes on the effectiveness of rotavirus vaccines has not been fully explored, this analysis underlines the potential of rotavirus whole-genome analysis by enhancing knowledge and understanding of intergenogroup reassortant strains circulating in Rwanda post vaccine introduction.

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.

Rotavirus Genotypes in Hospitalized Children with Acute Gastroenteritis Before and After Rotavirus Vaccine Introduction in Blantyre, Malawi, 1997 - 2019

Background

Rotavirus vaccine (Rotarix [RV1]) has reduced diarrhea-associated hospitalizations and deaths in Malawi. We examined the trends in circulating rotavirus genotypes in Malawi over a 22-year period to assess the impact of RV1 introduction on strain distribution.

Methods

Data on rotavirus-positive stool specimens among children aged <5 years hospitalized with diarrhea in Blantyre, Malawi before (July 1997–October 2012, n = 1765) and after (November 2012–October 2019, n = 934) RV1 introduction were analyzed. Rotavirus G and P genotypes were assigned using reverse-transcription polymerase chain reaction.

Results

A rich rotavirus strain diversity circulated throughout the 22-year period; Shannon (H′) and Simpson diversity (D′) indices did not differ between the pre- and postvaccine periods (H′ P < .149; D′ P < .287). Overall, G1 (n = 268/924 [28.7%]), G2 (n = 308/924 [33.0%]), G3 (n = 72/924 [7.7%]), and G12 (n = 109/924 [11.8%]) were the most prevalent genotypes identified following RV1 introduction. The prevalence of G1P[8] and G2P[4] genotypes declined each successive year following RV1 introduction, and were not detected after 2018. Genotype G3 reemerged and became the predominant genotype from 2017 onward. No evidence of genotype selection was observed 7 years post–RV1 introduction.

Conclusions

Rotavirus strain diversity and genotype variation in Malawi are likely driven by natural mechanisms rather than vaccine pressure.

Report of the 1st African Enteric Viruses Genome Initiative (AEVGI) Data and Bioinformatics Workshop on whole-genome analysis of some African rotavirus strains held in Bloemfontein, South Africa

The University of the Free State - Next Generation Sequencing (NGS) Unit, Bloemfontein, South Africa, hosted a data and bioinformatics workshop from 19 to 22 June 2018. The workshop was coordinated by the African Enteric Viruses Genome Initiative (AEVGI) with support from the Bill & Melinda Gates Foundation. The event introduced technologies in NGS and data analysis with focus on the rotavirus (RV) genome. The workshop fostered interactions and networking between professionals, scientific experts, technicians and students. The courses provided an overview of RV diarrhoea and its burden in Africa, while highlighting the key resources and methodologies in NGS and advanced bioinformatics in deciphering vaccine impact. It was concluded that, despite the reported significant decline in RV associated-diarrhoea mortality and morbidity in Africa due to RV vaccine impact, the need for continuous surveillance and genomic characterization to better understand the ever-changing dynamics of RV strains is imperative.

Whole-genome characterization of G12 rotavirus strains detected in Mozambique reveals a co-infection with a GXP[14] strain of possible animal origin

A high prevalence of G12 rotavirus strains has previously been reported in southern Mozambique. In this study, the full genomes of five Mozambican G12 strains were determined directly from stool using an Illumina Miseq platform. One sample (0060) contained an intergenogroup co-infection of a G12P[8] Wa-like strain and a GXP[14] DS-1-like strain. The sequences of seven genome segments, detected for the GXP[14] strain, clustered with a diverse group of mostly animal strains, suggesting co-infection with a strain of possible animal origin. The stool samples contained G12P[6] rotavirus strains with Wa-like backbones. Phylogenetic analyses of the VP4 and VP7 encoding segments of the G12P[6] strains suggested that they were reassortants containing backbones that are similar to that of the G12P[8] strain. The study confirms previous observations of interspecies transmission and emphasizes the importance of whole-genome sequencing in order to evaluate rotavirus co-infections and reassortants.

Whole genome analyses of DS-1-like Rotavirus A strains detected in children with acute diarrhoea in southern Mozambique suggest several reassortment events

We report the first whole genome constellations of Mozambican rotavirus A strains detected between 2012 and 2013 in the Mavalane General Hospital in Maputo city and Manhiça District Hospital in the Manhiça district. Consensus sequences for ten DS-1-like strains (G2P[4] and G8P[4]) were identified with an Illumina Miseq platform using cDNA prepared from dsRNA extracted from stool samples, without genome amplification or prior adaptation to cell culture. Comparison of previously reported genotyping results and the consensus sequences described in this study, indicated that the genotype primers specific for G12 and P[4] might require revision. Phylogenetic analyses indicated diversity among the G2P[4] Mozambican strains and suggested reassortment between G2P[4] and G8P[4] Mozambican strains, as well as the intragenogroup reassortment of all the genome segments encoding VP1, 2, 3 and 6 for strain RVA/Human-wt/MOZ/0045/2012G8P[4]. These results highlight the necessity to determine whole genome constellations to confirm surveillance data in Africa and to monitor the growing diversity in DS-1-like strains.

Whole-genome sequencing and analyses identify high genetic heterogeneity, diversity and endemicity of rotavirus genotype P[6] strains circulating in Africa

Rotavirus A (RVA) exhibits a wide genotype diversity globally. Little is known about the genetic composition of genotype P[6] from Africa. This study investigated possible evolutionary mechanisms leading to genetic diversity of genotype P[6] VP4 sequences. Phylogenetic analyses on 167 P[6] VP4 full-length sequences were conducted, which included six porcine-origin sequences. Of the 167 sequences, 57 were newly acquired through whole genome sequencing as part of this study. The other 110 sequences were all publicly-available global P[6] VP4 full-length sequences downloaded from GenBank. The strength of association between the phenotypic features and the phylogeny was also determined. A number of reassortment and mixed infections of RVA genotype P[6] strains were observed in this study. Phylogenetic analyses demostrated the extensive genetic diversity that exists among human P[6] strains, porcine-like strains, their concomitant clades/subclades and estimated that P[6] VP4 gene has a higher substitution rate with the mean of 1.05E-3 substitutions/site/year. Further, the phylogenetic analyses indicated that genotype P[6] strains were endemic in Africa, characterised by an extensive genetic diversity and long-time local evolution of the viruses. This was also supported by phylogeographic clustering and G-genotype clustering of the P[6] strains when Bayesian Tip-association Significance testing (BaTS) was applied, clearly supporting that the viruses evolved locally in Africa instead of spatial mixing among different regions. Overall, the results demonstrated that multiple mechanisms such as reassortment events, various mutations and possibly interspecies transmission account for the enormous diversity of genotype P[6] strains in Africa. These findings highlight the need for continued global surveillance of rotavirus diversity.

Complete genome analyses of the first porcine rotavirus group H identified from a South African pig does not provide evidence for recent interspecies transmission events

Rotaviruses (RVs) are classified into eight species/groups (RVA-RVH) according to the migration patterns of their 11 genome segments, as well as by serological and molecular properties of Viral Protein 6 (VP6). In 1997 a new unclassified RV was reported infecting adults in Bangladesh and China. This virus was initially named novel adult diarrhoea rotavirus (ADRV-N), but later renamed as RVH. Since then, RVH has been detected in humans only very sporadically. However, RVH is increasingly being detected in pig populations in the USA, Brazil and Japan, but not yet in Africa. Unfortunately, whole genome sequence data of porcine RVH strains in GenBank is currently restricted to a single strain (SKA-1) from Japan. Porcine diarrhoeic samples were collected in South Africa and analysed for rotavirus using an RVA ELISA and electropherotyping by PAGE. One sample displayed a 4:2:1:1:1:1:1 migration pattern, typical for RVH. In order to further investigate this strain, sequence-independent amplification followed by random sequencing using the 454/Roche GS FLX Sequencer was performed, resulting in the second complete porcine RVH strain (MRC-DPRU1575) available in databases. Phylogenetically, all segments of MRC-DPRU1575 clustered closely with the SKA-1 strain and in some segments with known porcine RVH strains from Brazil and the USA. In contrast, the porcine RVH strains were only distantly related to human RVH strains from Asia and a partial RVH-like strain recently detected in bats from Cameroon. Overall, strain MRC-DPRU1575 is the first complete genome of a porcine RVH from Africa and allows for the development of improved RVH screening methods. Our analyses indicate that RVH strains cluster according to their host species, not suggesting any evidence of recent interspecies transmission events. However, more RVH genomes from a wider host range are needed to better understand their evolutionary pathways and zoonotic potential.

Whole genome detection of rotavirus mixed infections in human, porcine and bovine samples co-infected with various rotavirus strains collected from sub-Saharan Africa

Group A rotaviruses (RVA) are among the main global causes of severe diarrhea in children under the age of 5years. Strain diversity, mixed infections and untypeable RVA strains are frequently reported in Africa. We analysed rotavirus-positive human stool samples (n=13) obtained from hospitalised children under the age of 5years who presented with acute gastroenteritis at sentinel hospital sites in six African countries, as well as bovine and porcine stool samples (n=1 each), to gain insights into rotavirus diversity and evolution. Polyacrylamide gel electrophoresis (PAGE) analysis and genotyping with G-(VP7) and P-specific (VP4) typing primers suggested that 13 of the 15 samples contained more than 11 segments and/or mixed G/P genotypes. Full-length amplicons for each segment were generated using RVA-specific primers and sequenced using the Ion Torrent and/or Illumina MiSeq next-generation sequencing platforms. Sequencing detected at least one segment in each sample for which duplicate sequences, often having distinct genotypes, existed. This supported and extended the PAGE and RT-PCR genotyping findings that suggested these samples were collected from individuals that had mixed rotavirus infections. The study reports the first porcine (MRC-DPRU1567) and bovine (MRC-DPRU3010) mixed infections. We also report a unique genome segment 9 (VP7), whose G9 genotype belongs to lineage VI and clusters with porcine reference strains. Previously, African G9 strains have all been in lineage III. Furthermore, additional RVA segments isolated from humans have a clear evolutionary relationship with porcine, bovine and ovine rotavirus sequences, indicating relatively recent interspecies transmission and reassortment. Thus, multiple RVA strains from sub-Saharan Africa are infecting mammalian hosts with unpredictable variations in their gene segment combinations. Whole-genome sequence analyses of mixed RVA strains underscore the considerable diversity of rotavirus sequences and genome segment combinations that result from a complex evolutionary history involving multiple host species.

Whole-genome analyses of DS-1-like human G2P[4] and G8P[4] rotavirus strains from Eastern, Western and Southern Africa

Group A rotaviruses (RVAs) with distinct G and P genotype combinations have been reported globally. We report the genome composition and possible origin of seven G8P[4] and five G2P[4] human RVA strains based on the genetic evolution of all 11 genome segments at the nucleotide level. Twelve RVA ELISA positive stool samples collected in the representative countries of Eastern, Southern and West Africa during the 2007-2012 surveillance seasons were subjected to sequencing using the Ion Torrent PGM and Illumina MiSeq platforms. A reference-based assembly was performed using CLC Bio's clc_ref_assemble_long program, and full-genome consensus sequences were obtained. With the exception of the neutralising antigen, VP7, all study strains exhibited the DS-1-like genome constellation (P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2) and clustered phylogenetically with reference strains having a DS-1-like genetic backbone. Comparison of the nucleotide and amino acid sequences with selected global cognate genome segments revealed nucleotide and amino acid sequence identities of 81.7-100 % and 90.6-100 %, respectively, with NSP4 gene segment showing the most diversity among the strains. Bayesian analyses of all gene sequences to estimate the time of divergence of the lineage indicated that divergence times ranged from 16 to 44 years, except for the NSP4 gene where the lineage seemed to arise in the more distant past at an estimated 203 years ago. However, the long-term effects of changes found within the NSP4 genome segment should be further explored, and thus we recommend continued whole-genome analyses from larger sample sets to determine the evolutionary mechanisms of the DS-1-like strains collected in Africa.

Whole genome analyses of G1P[8] rotavirus strains from vaccinated and non-vaccinated South African children presenting with diarrhea

Group A rotaviruses (RVAs) are the leading cause of severe gastroenteritis and eventually death among infants and young children worldwide, and disease prevention and management through vaccination is a public health priority. In August 2009, Rotarix™ was introduced in the South African Expanded Programme on Immunisation. As a result, substantial reductions in RVA disease burden have been reported among children younger than 5 years old. Rotavirus strain surveillance post-vaccination is crucial to, inter alia, monitor and study the evolution of vaccine escape strains. Here, full-genome sequence data for the 11 gene segments from 11 South African G1P[8] rotavirus strains were generated, including 5 strains collected from non-vaccinated children during the 2004-2009 rotavirus seasons and 6 strains collected from vaccinated children during the 2010 rotavirus season. These data were analyzed to gain insights into the overall genetic makeup and evolution of South African G1P[8] rotavirus strains and to compare their genetic backbones with those of common human Wa-like RVAs from other countries, as well as with the Rotarix™ and RotaTeq™ G1P[8] vaccine components. All 11 South African G1P[8] strains revealed a complete Wa-like genotype constellation of G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1. On the basis of sequence similarities, the South African G1P[8] strains (with the exception of strain RVA/Human-wt/ZAF/1262/2004/G1P[8]) were closely related to each other (96-100% identity in all gene segments). Comparison to the Rotarix™ and RotaTeq™ G1P[8] vaccine components revealed a moderate nucleotide identity of 89-96% and 93-95%, respectively. The results indicated that none of the gene segments of these 11 South African G1P[8] strains were vaccine-derived. This study illustrates that large-scale next generation sequencing will provide crucial information on the influence of the vaccination program on evolution of rotavirus strains. This is the first report to describe full genomic analyses of G1P[8] RVA strains collected from both non-vaccinated and vaccinated children in South Africa.

Genetic diversity of rotavirus genome segment 6 (encoding VP6) in Pretoria, South Africa

BACKGROUND

Rotavirus viral protein 6 (VP6), encoded by genome segment (GS) 6, is the primary target for rotavirus diagnosis by serological and some molecular techniques. Selected full length nucleotide sequences of GS 6 of rotavirus strains from South Africa were sequenced and analysed to determine genetic diversity and variations within the circulating rotaviruses.

FINDINGS

The VP6 amplicons were sequenced using the Sanger ABI 3130xl. Phylogenetic and pairwise analysis revealed that the VP6 genes of the study strains belonged to two different VP6 [I] genotypes. Five sequences were assigned genotype I1 and seven as genotype I2. Comparison of the group specific antigenic regions of the South African strains to the reference strains, shows that the South African VP6 sequences belonging to the VP6 genotype I2 were highly conserved, with only two amino acids changes at positions 239 (T›N) and 261(I›V). On the other hand, South African VP6 sequences belonging to I1 genotypes revealed several amino acid variations mostly within the antigenic region III.

CONCLUSIONS

Rotavirus strains with I1 and I2 genotype are predominantly circulating within the South African communities of which the later seems to be more conserved within the antigenic regions. The observed genetic variations observed within GS 6 of rotaviruses analysed in the current study are unlikely to impact negatively on the performance of the current VP6-based detection methods. Nevertheless, investigators should continually consider this diversity and adapt the primer design for the detection and characterization of the VP6 gene accordingly.

Novel NSP1 genotype characterised in an African camel G8P[11] rotavirus strain

Animal-human interspecies transmission is thought to play a significant role in influencing rotavirus strain diversity in humans. Proving this concept requires a better understanding of the complete genetic constellation of rotaviruses circulating in various animal species. However, very few whole genomes of animal rotaviruses, especially in developing countries, are available. In this study, complete genetic configuration of the first African camel rotavirus strain (RVA/Camel-wt/SDN/MRC-DPRU447/2002/G8P[11]) was assigned a unique G8-P[11]-I2-R2-C2-M2-A18-N2-T6-E2-H3 genotype constellation that has not been reported in other ruminants. It contained a novel NSP1 genotype (genotype A18). The evolutionary dynamics of the genome segments of strain MRC-DPRU447 were rather complex compared to those found in other camelids. Its genome segments 1, 3, 7-10 were closely related (>93% nucleotide identity) to those of human-animal reassortant strains like RVA/Human-tc/ITA/PA169/1988/G6P[14] and RVA/Human-wt/HUN/Hun5/1997/G6P[14], segments 4, 6 and 11 shared common ancestry (>95% nucleotide identity) with bovine rotaviruses like strains RVA/Cow-wt/CHN/DQ-75/2008/G10P[11] and RVA/Cow-wt/KOR/KJ19-2/XXXX/G6P[7], whereas segment 2 was closely related (94% nucleotide identity) to guanaco rotavirus strain RVA/Guanaco-wt/ARG/Rio_Negro/1998/G8P[1]. Its genetic backbone consisted of DS-1-like, AU-1-like, artiodactyl-like and a novel A18 genotype. This suggests that strain MRC-DPRU447 potentially emerged through multiple reassortment events between several mammalian rotaviruses of at least two genogroups or simply strain MRC-DPRU447 display a unique progenitor genotypes. Close relationship between some of the genome segments of strain MRC-DPRU447 to human rotaviruses suggests previous occurrence of reassortment processes combined with interspecies transmission between humans and camels. The whole genome data for strain MRC-DPRU447 adds to the much needed animal rotavirus data from Africa which is limited at the moment.