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

Emergence of a Novel G4P[6] Porcine Rotavirus with Unique Sequence Duplication in NSP5 Gene in China

Rotavirus is a major causative agent of diarrhoea in children, infants, and young animals around the world. The associated zoonotic risk necessitates the serious consideration of the complete genetic information of rotavirus. A segmented genome makes rotavirus prone to rearrangement and the formation of a new viral strain. Monitoring the molecular epidemiology of rotavirus is essential for its prevention and control. The quantitative RT-PCR targeting the NSP5 gene was used to detect rotavirus group A (RVA) in pig faecal samples, and two pairs of universal primers and protocols were used for amplifying the G and P genotype. The genotyping and phylogenetic analysis of 11 genes were performed by RT-PCR and a basic bioinformatics method. A unique G4P[6] rotavirus strain, designated S2CF (RVA/Pig-tc/CHN/S2CF/2023/G4P[6]), was identified in one faecal sample from a piglet with severe diarrhoea in Guangdong, China. Whole genome sequencing and analysis suggested that the 11 segments of the S2CF strain showed a unique Wa-like genotype constellation and a typical porcine RVA genomic configuration of G4-P[6]-I1-R1-C1-M1-A8-N1-T1-E1-H1. Notably, 4 of the 11 gene segments (VP4, VP6, VP2, and NSP5) clustered consistently with human-like RVAs, suggesting independent human-to-porcine interspecies transmission. Moreover, a unique 344-nt duplicated sequence was identified for the first time in the untranslated region of NSP5. This study further reveals the genetic diversity and potential inter-species transmission of porcine rotavirus.

Genomic characterization of the rotavirus G3P[8] strain in vaccinated children, reveals possible reassortment events between human and animal strains in Manhiça District, Mozambique

Mozambique introduced the rotavirus vaccine (Rotarix®; GlaxoSmithKline Biologicals, Rixensart, Belgium) in 2015, and since then, the Centro de Investigação em Saúde de Manhiça has been monitoring its impact on rotavirus-associated diarrhea and the trend of circulating strains, where G3P[8] was reported as the predominant strain after the vaccine introduction. Genotype G3 is among the most commonly detected Rotavirus strains in humans and animals, and herein, we report on the whole genome constellation of G3P[8] detected in two children (aged 18 months old) hospitalized with moderate-to-severe diarrhea at the Manhiça District Hospital. The two strains had a typical Wa-like genome constellation (I1-R1-C1-M1-A1-N1-T1-E1-H1) and shared 100% nucleotide (nt) and amino acid (aa) identities in 10 gene segments, except for VP6. Phylogenetic analysis demonstrated that genome segments encoding VP7, VP6, VP1, NSP3, and NSP4 of the two strains clustered most closely with porcine, bovine, and equine strains with identities ranging from 86.9-99.9% nt and 97.2-100% aa. Moreover, they consistently formed distinct clusters with some G1P[8], G3P[8], G9P[8], G12P[6], and G12P[8] strains circulating from 2012 to 2019 in Africa (Mozambique, Kenya, Rwanda, and Malawi) and Asia (Japan, China, and India) in genome segments encoding six proteins (VP2, VP3, NSP1-NSP2, NSP5/6). The identification of segments exhibiting the closest relationships with animal strains shows significant diversity of rotavirus and suggests the possible occurrence of reassortment events between human and animal strains. This demonstrates the importance of applying next-generation sequencing to monitor and understand the evolutionary changes of strains and evaluate the impact of vaccines on strain diversity.

Interspecies transmission of porcine-originated G4P[6] rotavirus A between pigs and humans: a synchronized spatiotemporal approach

As a leading viral cause of acute gastroenteritis in both humans and pigs, rotavirus A (RVA) poses a potential public health concern. Although zoonotic spillover of porcine RVA strains to humans is sporadic, it has been detected worldwide. The origin of chimeric human-animal strains of RVA is closely linked to the crucial role of mixed genotypes in driving reassortment and homologous recombination, which play a major role in shaping the genetic diversity of RVA. To better understand how genetically intertwined porcine and zoonotic human-derived G4P[6] RVA strains are, the present study employed a spatiotemporal approach to whole-genome characterization of RVA strains collected during three consecutive RVA seasons in Croatia (2018-2021). Notably, sampled children under 2 years of age and weanling piglets with diarrhea were included in the study. In addition to samples tested by real-time RT-PCR, genotyping of VP7 and VP4 gene segments was conducted. The unusual genotype combinations detected in the initial screening, including three human and three porcine G4P[6] strains, were subjected to next-generation sequencing, followed by phylogenetic analysis of all gene segments, and intragenic recombination analysis. Results showed a porcine or porcine-like origin for each of the eleven gene segments in all six RVA strains. The G4P[6] RVA strains detected in children most likely resulted from porcine-to-human interspecies transmission. Furthermore, the genetic diversity of Croatian porcine and porcine-like human G4P[6] strains was propelled by reassortment events between porcine and porcine-like human G4P[6] RVA strains, along with homologous intragenotype and intergenotype recombinations in VP4, NSP1, and NSP3 segments. Described concurrent spatiotemporal approach in investigating autochthonous human and animal RVA strains is essential in drawing relevant conclusions about their phylogeographical relationship. Therefore, continuous surveillance of RVA, following the One Health principles, may provide relevant data for assessing the impact on the protectiveness of currently available vaccines.

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

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

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

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

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

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

Maintaining Momentum for Rotavirus Immunization in Africa during the COVID-19 Era: Report of the 13th African Rotavirus Symposium

The 13th African Rotavirus Symposium was held as a virtual event hosted by the University of Nairobi, Kenya and The Kenya Paediatric Association on 3rd and 4th November 2021. This biennial event organized under the auspices of the African Rotavirus Network shapes the agenda for rotavirus research and prevention on the continent, attracting key international and regional opinion leaders, researchers, and public health scientists. The African Rotavirus Network is a regional network of institutions initially established in 1999, and now encompassing much of the diarrheal disease and rotavirus related research in Africa, in collaboration with the World Health Organization African Regional Office (WHO-AFRO), Ministries of Health, and other partners. Surges in SARS-CoV2 variants and concomitant travel restrictions limited the meeting to a webinar platform with invited scientific presentations and scientific presentations from selected abstracts. The scientific program covered updates on burden of diarrheal diseases including rotavirus, the genomic characterization of rotavirus strains pre- and post-rotavirus vaccine introduction, and data from clinical evaluation of new rotavirus vaccines in Africa. Finally, 42 of the 54 African countries have fully introduced rotavirus vaccination at the time of the meeting, including the two recently WHO pre-qualified vaccines from India. Nonetheless, the full benefit of rotavirus vaccination is yet to be realized in Africa where approximately 80% of the global burden of rotavirus mortality exists.

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 rotavirus strains circulating in Benin before vaccine introduction, 2016-2018

Species A Rotaviruses (RVA) still play a major role in causing acute diarrhea in children under five years old worldwide. Currently, an 11-gene classification system is used to designate the full genotypic constellations of circulating strains. Viral proteins and non-structural proteins in the order VP7-VP4-VP6-VP1-VP2-VP3-NSP1-NSP2-NSP3-NSP4-NSP5/6 are represented by the genotypes Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx, respectively. In Benin, ROTAVAC® vaccine was introduced into the Expanded Programme on Immunization in December 2019. To monitor circulating RVA strains for changes that may affect vaccine performance, in-depth analysis of strains prior to vaccine introduction are needed. Here we report, the whole-gene characterization (11 ORFs) for 72 randomly selected RVA strains of common and unusual genotypes collected in Benin from the 2016-2018 seasons. The sequenced strains were 15 G1P[8], 20 G2P[4], 5 G9P[8], 14 G12P[8], 9 G3P[6], 2 G1P[6], 3 G2P[6], 2 G9P[4], 1 G12P[6], and 1 G1G9P[8]/P[4]. The study strains exhibited two genetic constellations designed as Wa-like G1/G9/G12-P[6]/P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 and DS-1-like G2/G3/G12-P[4]/P[6]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Genotype G9P[4] strains possessed a DS-1-like genetic constellation with an E6 NSP4 gene, G9-P[4]-I2-R2-C2-M2-A2-N2-T2-E6-H2. The mixed genotype showed both Wa-like and DS-1-like profiles with a T6 NSP3 gene G1/G9P[8]/[4]-I1/I2-R1/R2-C1/C2-M1/M2-A1/A2-N1/N2-T1/T6-E1/E6-H1/H2. At the allelic level, the analysis of the Benin strains, reference strains (with known alleles), vaccine strains (with known alleles) identified 2-13 and 1-17 alleles for DS-1-like and Wa-like strains, respectively. Most of the study strains clustered into previously defined alleles, but we defined 3 new alleles for the VP7 (G3=1 new allele and G12=2 new alleles) and VP4 (P[4]=1 new allele and P[6]=2 new alleles) genes which formed the basis of the VP7 and VP4 gene clusters, respectively. For the remaining 9 genes, 0-6 new alleles were identified for both Wa-like and DS-1-like strains. This analysis of whole genome sequences of RVA strains circulating in Benin described genetic point mutations and reassortment events as well as novel alleles. Further detailed studies on these new alleles are needed and these data can also provide a baseline for studies on RVA in the post-vaccination period.

Genomic characterization of an African G4P[6] human rotavirus strain identified in a diarrheic child in Kenya: Evidence for porcine-to-human interspecies transmission and reassortment

Human rotavirus strains having the unconventional G4P[6] genotype have been sporadically identified in diarrheic patients in different parts of the world. However, the whole genome of only one human G4P[6] strain from Africa (central Africa) has been sequenced and analyzed, and thus the exact origin and evolutionary pattern of African G4P[6] strains remain to be elucidated. In this study, we characterized the full genome of an African G4P[6] strain (RVA/Human-wt/KEN/KCH148/2019/G4P[6]) identified in a stool specimen from a diarrheic child in Kenya. Full genome analysis of strain KCH148 revealed a unique Wa-like genogroup constellation: G4-P[6]-I1-R1-C1-M1-A1-N1-T7-E1-H1. NSP3 genotype T7 is commonly found in porcine rotavirus strains. Furthermore, phylogenetic analysis showed that 10 of the 11 genes of strain KCH148 (VP7, VP4, VP6, VP1-VP3, NSP1, and NSP3-NSP5) appeared to be of porcine origin, the remaining NSP2 gene appearing to be of human origin. Therefore, strain KCH148 was found to have a porcine rotavirus backbone and thus is likely to be of porcine origin. Furthermore, strain KCH148 is assumed to have been derived through interspecies transmission and reassortment events involving porcine and human rotavirus strains. To our knowledge, this is the first report on full genome-based characterization of a human G4P[6] strain from east Africa. Our observations demonstrated the diversity of human G4P[6] strains in Africa, and provide important insights into the origin and evolutionary pattern of zoonotic G4P[6] strains on the African continent.

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.

A decade of rotavirus vaccination in Africa - Saving lives and changing the face of diarrhoeal diseases: Report of the 12th African Rotavirus Symposium

The African Rotavirus Network organised the 12^th African Rotavirus Symposium (ARS) from 30 July to 1 August 2019 in Johannesburg, South Africa. The symposium theme "A decade of rotavirus vaccination in Africa - Saving lives and changing the face of diarrhoeal diseases", included sessions aimed at sharing ideas and expertise on prevention and control of diarrhoeal disease in Africa. Inter alia, the delegates reviewed global and regional epidemiological trends on rotavirus diarrhoea, progress and experiences on rotavirus vaccine introduction, including vaccine safety monitoring and impact in Africa, scientific advances in developing newer rotavirus vaccines, surveillance and research on other diarrhoeal pathogens, and providing an enabling environment for networking. Importantly, the 12^th ARS served to commemorate the 20^th anniversary of the African Rotavirus Network (AfrRN) coinciding with the 50th anniversary of the South African Medical Research Council. Four oral, live-attenuated rotavirus vaccines are currently prequalified by the WHO (Rotarix, RotaTeq, Rotavac and RotaSiil). African countries utilising rotavirus vaccines in routine national immunisation programmes are realising their effectiveness and impact on diarrhoeal disease morbidity. An ~40% reduction in hospitalisations of <5-year-olds with acute gastroenteritis following rotavirus vaccine introduction, was reported between 2006 and 2018 in 92,000 children from the WHO-coordinated African Rotavirus Surveillance Network (AfrRSN) comprising 33 Member States. This was corroborated by a meta-analysis of published data, sourced from January 2000 to August 2018 that reported substantial reductions in rotavirus hospitalisations in countries using rotavirus vaccines. However, it was highlighted that the transition of some countries from Gavi-eligibility and vaccine supply shortfalls present significant challenges to achieving the full impact of rotavirus immunization in Africa. The wide diversity of rotavirus genotypes continues in Africa, with variation observed both geographically and temporally. There is currently no evidence to suggest that the emergence of rotavirus strains not included in the current vaccines do escape vaccine-induced immunity.

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.