Whole-genome characterisation of G12P[6] rotavirus strains possessing two distinct genotype constellations co-circulating in Blantyre, Malawi, 2008

Genetic diversity and evolution of G12P[6] DS-1-like and G12P[9] AU-1-like Rotavirus strains in Brazil

In the early 2000s, the global emergence of rotavirus (RVA) G12P[8] genotype was noted, while G12P[6] and G12P[9] combinations remained rare in humans. This study aimed to characterize and phylogenetically analyze three Brazilian G12P[9] and four G12P[6] RVA strains from 2011 to 2020, through RT-PCR and sequencing, in order to enhance our understanding of the genetic relationship between human and animal-origin RVA strains. G12P[6] strains displayed a DS-1-like backbone, showing a distinct genetic clustering. G12P[6] IAL-R52/2020, IAL-R95/2020 and IAL-R465/2019 strains clustered with 2019 Northeastern G12P[6] Brazilian strains and a 2018 Benin strain, whereas IAL-R86/2011 strain grouped with 2010 Northern G12P[6] Brazilian strains and G2P[4] strains from the United States and Belgium. These findings suggest an African genetic ancestry and reassortments with co-circulating American strains sharing the same DS-1-like constellation. No recent zoonotic reassortment was observed, and the DS-1-like constellation detected in Brazilian G12P[6] strains does not seem to be genetically linked to globally reported intergenogroup G1/G3/G9/G8P[8] DS-1-like human strains. G12P[9] strains exhibited an AU-1-like backbone with two different genotype-lineage constellations: IAL-R566/2011 and IAL-R1151/2012 belonged to a VP3/M3.V Lineage, and IAL-R870/2013 to a VP3/M3.II Lineage, suggesting two co-circulating strains in Brazil. This genetic diversity is not observed elsewhere, and the VP3/M3.II Lineage in G12P[9] strains seems to be exclusive to Brazil, indicating its evolution within the country. All three G12P[9] AU-1-like strains were closely relate to G12P[9] strains from Paraguay (2006-2007) and Brazil (2010). Phylogenetic analysis also highlighted that all South American G12P[9] AU-1-like strains had a common origin and supports the hypothesis of their importation from Asia, with no recent introduction from globally circulating G12P[9] strains or reassortments with local G12 strains P[8] or P[6]. Notably, certain genes in the Brazilian G12P[9] AU-1-like strains share ancestry with feline/canine RVAs (VP3/M3.II, NSP4/E3.IV and NSP2/N3.II), whereas NSP1/A3.VI likely originated from artiodactyls, suggesting a history of zoonotic transmission with human strains. This genomic data adds understanding to the molecular epidemiology of G12P[6] and G12P[9] RVA strains in Brazil, offering insights into their genetic diversity and evolution.

VP4 Mutation Boosts Replication of Recombinant Human/Simian Rotavirus in Cell Culture

Rotavirus A (RVA) is the leading cause of diarrhea requiring hospitalization in children and causes over 100,000 annual deaths in Sub-Saharan Africa. In order to generate next-generation vaccines against African RVA genotypes, a reverse genetics system based on a simian rotavirus strain was utilized here to exchange the antigenic capsid proteins VP4, VP7 and VP6 with those of African human rotavirus field strains. One VP4/VP7/VP6 (genotypes G9-P[6]-I2) triple-reassortant was successfully rescued, but it replicated poorly in the first cell culture passages. However, the viral titer was enhanced upon further passaging. Whole genome sequencing of the passaged virus revealed a single point mutation (A797G), resulting in an amino acid exchange (E263G) in VP4. After introducing this mutation into the VP4-encoding plasmid, a VP4 mono-reassortant as well as the VP4/VP7/VP6 triple-reassortant replicated to high titers already in the first cell culture passage. However, the introduction of the same mutation into the VP4 of other human RVA strains did not improve the rescue of those reassortants, indicating strain specificity. The results show that specific point mutations in VP4 can substantially improve the rescue and replication of recombinant RVA reassortants in cell culture, which may be useful for the development of novel vaccine strains.

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 Analysis of African G12P[6] and G12P[8] Rotaviruses Provides Evidence of Porcine-Human Reassortment at NSP2, NSP3, and NSP4

Group A rotaviruses (RVA) represent the most common cause of pediatric gastroenteritis in children <5 years, worldwide. There has been an increase in global detection and reported cases of acute gastroenteritis caused by RVA genotype G12 strains, particularly in Africa. This study sought to characterize the genomic relationship between African G12 strains and determine the possible origin of these strains. Whole genome sequencing of 34 RVA G12P[6] and G12P[8] strains detected from the continent including southern (South Africa, Zambia, Zimbabwe), eastern (Ethiopia, Uganda), central (Cameroon), and western (Togo) African regions, were sequenced using the Ion Torrent PGM method. The majority of the strains possessed a Wa-like backbone with consensus genotype constellation of G12-P[6]/P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1, while a single strain from Ethiopia displayed a DS-1-like genetic constellation of G12-P[6]-I2-R2-C2-M2-A2-N2-T2-E2-H2. In addition, three Ethiopian and one South African strains exhibited a genotype 2 reassortment of the NSP3 gene, with genetic constellation of G12-P[8]-I1-R1-C1-M1-A1-N1-T2-E1-H1. Overall, 10 gene segments (VP1–VP4, VP6, and NSP1–NSP5) of African G12 strains were determined to be genetically related to cognate gene sequences from globally circulating human Wa-like G12, G9, and G1 strains with nucleotide (amino acid) identities in the range of 94.1–99.9% (96.5–100%), 88.5–98.5% (93–99.1%), and 89.8–99.0% (88.7–100%), respectively. Phylogenetic analysis showed that the Ethiopian G12P[6] possessing a DS-1-like backbone consistently clustered with G2P[4] strains from Senegal and G3P[6] from Ethiopia with the VP1, VP2, VP6, and NSP1–NSP4 genes. Notably, the NSP2, NSP3, and NSP4 of most of the study strains exhibited the closest relationship with porcine strains suggesting the occurrence of reassortment between human and porcine strains. Our results add to the understanding of potential roles that interspecies transmission play in generating human rotavirus diversity through reassortment events and provide insights into the evolutionary dynamics of G12 strains spreading across selected sub-Saharan Africa regions.

The origins of G12P[6] rotavirus strains detected in Lebanon

The G12 rotaviruses are an increasingly important cause of severe diarrhoea in infants and young children worldwide. Seven human G12P[6] rotavirus strains were detected in stool samples from children hospitalized with gastroenteritis in Lebanon during a 2011-2013 surveillance study. Complete genomes of these strains were sequenced using VirCapSeq-VERT, a capture-based high-throughput viral-sequencing method, and further characterized based on phylogenetic analyses with global RVA and vaccine strains. Based on the complete genomic analysis, all Lebanese G12 strains were found to have Wa-like genetic backbone G12-P[6]-I1-R1-C1-M1-A1-N1-T1-E1-H1. Phylogenetically, these strains fell into two clusters where one of them might have emerged from Southeast Asian strains and the second one seems to have a mixed backbone between North American and Southeast Asian strains. Further analysis of these strains revealed high antigenic variability compared to available vaccine strains. To our knowledge, this is the first report on the complete genome-based characterization of G12P[6] emerging in Lebanon. Additional studies will provide important insights into the evolutionary dynamics of G12 rotaviruses spreading in Asia.

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.

Group A rotavirus surveillance before vaccine introduction in Italy, September 2014 to August 2017

Introduction

Group A rotaviruses (RVA) are the leading cause of acute gastroenteritis (AGE) in young children, causing ca 250,000 deaths worldwide, mainly in low-income countries. Two proteins, VP7 (glycoprotein, G genotype) and VP4 (protease-sensitive protein, P genotype), are the basis for the binary RVA nomenclature. Although 36 G types and 51 P types are presently known, most RVA infections in humans worldwide are related to five G/P combinations: G1P[8], G2P[4], G3P[8], G4P[8], G9P[8].

Aim

This study aimed to characterise the RVA strains circulating in Italy in the pre-vaccination era, to define the trends of circulation of genotypes in the Italian paediatric population.

Methods

Between September 2014 and August 2017, after routine screening in hospital by commercial antigen detection kit, 2,202 rotavirus-positive samples were collected in Italy from children hospitalised with AGE; the viruses were genotyped following standard European protocols.

Results

This 3-year study revealed an overall predominance of the G12P[8] genotype (544 of 2,202 cases; 24.70%), followed by G9P[8] (535/2,202; 24.30%), G1P[8] (459/2,202; 20.84%) and G4P[8] (371/2,202; 16.85%). G2P[4] and G3P[8] genotypes were detected at low rates (3.32% and 3.09%, respectively). Mixed infections accounted for 6.49% of cases (143/2,202), uncommon RVA strains for 0.41% of cases (9/2,202).

Conclusions

The emergence of G12P[8] rotavirus in Italy, as in other countries, marks this genotype as the sixth most common human genotype. Continuous surveillance of RVA strains and monitoring of circulating genotypes are important for a better understanding of rotavirus evolution and genotype distribution, particularly regarding strains that may emerge from reassortment events.

Phylogeography and evolutionary analysis of African Rotavirus a genotype G12 reveals district genetic diversification within lineage III

Group A rotavirus (RVA) genotype G12 has spread globally and has become one of the most prevalent genotypes of rotavirus in Africa. To understand the drivers for its genetic diversity and rapid spread we investigated the Bayesian phylogeography, viral evolution and population demography of Rotavirus G12 in Africa. We downloaded and aligned VP7 gene sequences of Rotavirus genotype G12, from thirteen African countries (n = 96). Phylogenetic analysis, Evolutionary analysis and Bayesian Phylogeography was carried out, using MEGA Vs 6, BEAST, and SPREAD3. Phylogenetic analysis revealed that all the African sequences fell into lineage III diversifying into two major clades. The evolutionary rate of the African rotavirus G12 sequences was 1.678×10^-3, (95% HPD, 1.201×10^-3 - 2.198×10^-3) substitutions/site/year, with TMRC of 16.8 years. The Maximum clade credibility (MCC) tree clustered into three lineages (II, III, IV), African strains fell within lineage III, and diversified into three clusters. Phylogeography suggested that South Africa seemed to be the epicentre of dispersal of the genotype. The demographic history of the G12 viruses revealed a steady increase between the years1998-2007, followed by a sharp decrease in effective population size between the years 2008-2011. We have shown the potential for genetic diversification of Rotavirus genotype G12 in Africa. We recommend the adoption of Molecular surveillance across Africa to further control spread and diversification of Rotavirus.

Genetic analysis of Ghanaian G1P[8] and G9P[8] rotavirus A strains reveals the impact of P[8] VP4 gene polymorphism on P-genotyping

The World Health Organisation rotavirus surveillance networks have documented and shown eclectic geographic and temporal diversity in circulating G- and P- genotypes identified in children <5 years of age. To effectively monitor vaccine performance and effectiveness, robust molecular and phylogenetic techniques are essential to detect novel strain variants that might emerge due to vaccine pressure. This study inferred the phylogenetic history of the VP7 and VP4 genes of previously non-typeable strains and provided insight into the diversity of P[8] VP4 sequences which impacted the outcome of our routine VP4 genotyping method. Near-full-length VP7 gene and the VP8* fragment of the VP4 gene were obtained by Sanger sequencing and genotypes were determined using RotaC v2.0 web-based genotyping tool. The genotypes of the 57 rotavirus-positive samples with sufficient stool was determined. Forty-eight of the 57 (84.2%) had the P[8] specificity, of which 43 (89.6%) were characterized as P[8]a subtype and 5 (10.4%) as the rare OP354-like subtype. The VP7 gene of 27 samples were successfully sequenced and their G-genotypes confirmed as G1 (18/27) and G9 (9/27). Phylogenetic analysis of the P[8]a sequences placed them in subcluster IIIc within lineage III together with contemporary G1P[8], G3P[8], G8P[8], and G9P[8] strains detected globally from 2006-2016. The G1 VP7 sequences of the study strains formed a monophyletic cluster with African G1P[8] strains, previously detected in Ghana and Mali during the RotaTeq vaccine trial as well as Togo. The G9 VP7 sequences of the study strains formed a monophyletic cluster with contemporary African G9 sequences from neighbouring Burkina Faso within the major sub-cluster of lineage III. Mutations identified in the primer binding region of the VP8* sequence of the Ghanaian P[8]a strains may have resulted in the genotyping failure since the newly designed primer successfully genotyped the previously non-typeable P[8] strains. In summary, the G1, G9, and P[8]a sequences were highly similar to contemporary African strains at the lineage level. The study also resolved the methodological challenges of the standard genotyping techniques and highlighted the need for regular evaluation of the multiplex PCR-typing method especially in the post-vaccination era. The study further highlights the need for regions to start using sequencing data from local rotavirus strains to design and update genotyping primers.

Sub-genotype phylogeny of the non-G, non-P genes of genotype 2 Rotavirus A strains

Recent increase in the detection of unusual G1P[8], G3P[8], G8P[8], and G9P[4] Rotavirus A (RVA) strains bearing the DS-1-like constellation of the non-G, non-P genes (hereafter referred to as the genotype 2 backbone) requires better understanding of their evolutionary relationship. However, within a genotype, there is lack of a consensus lineage designation framework and a set of common sequences that can serve as references. Phylogenetic analyses were carried out on over 8,500 RVA genotype 2 genes systematically retrieved from the rotavirus database within the NCBI Virus Variation Resource. In line with previous designations, using pairwise comparison of cogent nucleotide sequences and stringent bootstrap support, reference lineages were defined. This study proposes a lineage framework and provides a dataset ranging from 34 to 145 sequences for each genotype 2 gene for orderly lineage designation of global genotype 2 genes of RVAs detected in human and animals. The framework identified five to 31 lineages depending on the gene. The least number of lineages (five to seven) were observed in genotypes A2 (NSP1), T2 (NSP3) and H2 (NSP5) which are limited to human RVA whereas the most number of lineages (31) was observed in genotype E2 (NSP4). Sharing of the same lineage constellations of the genotype 2 backbone genes between recently-emerging, unusual G1P[8], G3P[8], G8P[8] and G9P[4] reassortants and many contemporary G2P[4] strains provided strong support to the hypothesis that unusual genotype 2 strains originated primarily from reassortment events in the recent past involving contemporary G2P[4] strains as one parent and ordinary genotype 1 strains or animal RVA strains as the other. The lineage framework with selected reference sequences will help researchers to identify the lineage to which a given genotype 2 strain belongs, and trace the evolutionary history of common and unusual genotype 2 strains in circulation.

Emergence of G12P[6] rotavirus strains among hospitalised children with acute gastroenteritis in Belém, Northern Brazil, following introduction of a rotavirus vaccine

Species A rotavirus still remains a major cause of acute gastroenteritis in infants and young children. Globally, six genotypes (G1P[8], G2P[4], G3P[8], G4P[8], G9P[8] and G12P[8]) account for >90% of circulating strains; however, genotype G12 in combination with P[6] or P[9] has been detected at increasing rates. We sought to broaden our knowledge about the rotavirus strains circulating during the early post-vaccine-introduction period. Stool samples were obtained from children hospitalised for acute gastroenteritis in Belém, Northern Brazil, from May 2008 to May 2011 and examined by reverse transcription polymerase chain reaction and nucleotide sequencing. A total of 122 out of the original 1076 rotavirus strains were judged to be non-typeable in the first analysis and were therefore re-examined. G2P[4] was the most prevalent genotype (58.0%), followed by G1P[8] (16.9%), and G12P[6] (7.5%). G12P[6] strains were identified at similar rates during the first (2.5%) and second (3.9%) years, and the rate jumped to 15.6% in the third year. Analysis of VP7 sequences of the G12P[6] strains showed that they belonged to lineage III. In addition, co-circulating G12P[6] strains displaying long and short RNA patterns were found to belong to the Wa-like and DS-1-like constellation, respectively. Additional unusual circulating strains G12P[9] and G3P[9] were also identified. This hospital-based study showed a high prevalence of G12P[6] strains in the third year of surveillance. Our results highlight the need for continuous longitudinal monitoring of circulating rotavirus strains after introduction of rotavirus vaccines in Brazil and elsewhere.

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.