A metagenomic investigation of the faecal RNA virome structure of asymptomatic chickens obtained from a commercial farm in Durban, KwaZulu-Natal province, South Africa.

BACKGROUND: Virome studies on birds, including chickens are relatively scarce, particularly from the African continent. Despite the continuous evolution of RNA viruses and severe losses recorded in poultry from seasonal viral outbreaks, the information on RNA virome composition is even scantier as a result of their highly unstable nature, genetic diversity, and difficulties associated with characterization. Also, information on factors that may modulate the occurrence of some viruses in birds is limited, particularly for domesticated birds. Viral metagenomics through advancements in sequencing technologies, has enabled the characterization of the entire virome of diverse host species using various samples. METHODS: The complex RNA viral constituents present in 27 faecal samples of asymptomatic chickens from a South African farm collected at 3-time points from two independent seasons were determined, and the impact of the chicken's age and collection season on viral abundance and diversity was further investigated. The study utilized the non-invasive faecal sampling method, mRNA viral targeted enrichment steps, a whole transcriptome amplification strategy, Illumina sequencing, and bioinformatics tools. RESULTS: The results obtained revealed a total of 48 viral species spanning across 11 orders, 15 families and 21 genera. Viral RNA families such as Coronaviridae, Picornaviridae, Reoviridae, Astroviridae, Caliciviridae, Picorbirnaviridae and Retroviridae were abundant, among which picornaviruses, demonstrated a 100% prevalence across the three age groups (2, 4 and 7 weeks) and two seasons (summer and winter) of the 27 faecal samples investigated. A further probe into the extent of variation between the different chicken groups investigated indicated that viral diversity and abundance were significantly influenced by age (P = 0.01099) and season (P = 0.00099) between chicken groups, while there was no effect on viral shedding within samples in a group (alpha diversity) for age (P = 0.146) and season (P = 0.242). CONCLUSION: The presence of an exceedingly varied chicken RNA virome, encompassing avian, mammalian, fungal, and dietary-associated viruses, underscores the complexities inherent in comprehending the causation, dynamics, and interspecies transmission of RNA viruses within the investigated chicken population. Hence, chickens, even in the absence of discernible symptoms, can harbour viruses that may exhibit opportunistic, commensal, or pathogenic characteristics.

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.

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.

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 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 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.

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.

High Detection Frequency of Vaccine-Associated Polioviruses and Non-Polio Enteroviruses in the Stools of Asymptomatic Infants from the Free State Province, South Africa.

Enterovirus (EV) infections are widespread and associated with a range of clinical conditions, from encephalitis to meningitis, gastroenteritis, and acute flaccid paralysis. Knowledge about the circulation of EVs in neonatal age and early infancy is scarce, especially in Africa. This study aimed to unveil the frequency and diversity of EVs circulating in apparently healthy newborns from the Free State Province, South Africa (SA). For this purpose, longitudinally collected faecal specimens (May 2021-February 2022) from a cohort of 17 asymptomatic infants were analysed using metagenomic next-generation sequencing. Overall, seven different non-polio EV (NPEV) subtypes belonging to EV-B and EV-C species were identified, while viruses classified under EV-A and EV-D species could not be characterised at the sub-species level. Additionally, under EV-C species, two vaccine-related poliovirus subtypes (PV1 and PV3) were identified. The most prevalent NPEV species was EV-B (16/17, 94.1%), followed by EV-A (3/17, 17.6%), and EV-D (4/17, 23.5%). Within EV-B, the commonly identified NPEV types included echoviruses 6, 13, 15, and 19 (E6, E13, E15, and E19), and coxsackievirus B2 (CVB2), whereas enterovirus C99 (EV-C99) and coxsackievirus A19 (CVA19) were the only two NPEVs identified under EV-C species. Sabin PV1 and PV3 strains were predominantly detected during the first week of birth and 6-8 week time points, respectively, corresponding with the OPV vaccination schedule in South Africa. A total of 11 complete/near-complete genomes were identified from seven NPEV subtypes, and phylogenetic analysis of the three EV-C99 identified revealed that our strains were closely related to other strains from Cameroon and Brazil, suggesting global distribution of these strains. This study provides an insight into the frequency and diversity of EVs circulating in asymptomatic infants from the Free State Province, with the predominance of subtypes from EV-B and EV-C species. This data will be helpful to researchers looking into strategies for the control and treatment of EV infection.

Genomic insights into the 2022-2023Vibrio cholerae outbreak in Malawi.

Malawi experienced its deadliest Vibrio cholerae (Vc) outbreak following devastating cyclones, with >58,000 cases and >1700 deaths reported between March 2022 and May 2023. Here, we use population genomics to investigate the attributes and origin of the Malawi 2022-2023 Vc outbreak isolates. Our results demonstrate the predominance of ST69 clone, also known as the seventh cholera pandemic El Tor (7PET) lineage, expressing O1 Ogawa (~ 80%) serotype followed by Inaba (~ 16%) and sporadic non-O1/non-7PET serogroups (~ 4%). Phylogenetic reconstruction revealed that the Malawi outbreak strains correspond to a recent importation from Asia into Africa (sublineage AFR15). These isolates harboured known antimicrobial resistance and virulence elements, notably the ICEGEN/ICEVchHai1/ICEVchind5 SXT/R391-like integrative conjugative elements and a CTXφ prophage with the ctxB7 genotype compared to historical Malawian Vc isolates. These data suggest that the devastating cyclones coupled with the recent importation of 7PET serogroup O1 strains, may explain the magnitude of 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.

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.

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 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.

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.

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.

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.

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.