Human P[6] Rotaviruses From Sub-Saharan Africa and Southeast Asia Are Closely Related to Those of Human P[4] and P[8] Rotaviruses Circulating Worldwide

Full genome-based characterization of an Asian G3P[6] human rotavirus strain found in a diarrheic child in Japan: Evidence for porcine-to-human zoonotic transmission

Human rotavirus strains having the unconventional G3P[6] genotype have been sporadically detected in diarrheic patients in different parts of the world. However, the full genomes of only three human G3P[6] strains from Asian countries (China, Indonesia, and Vietnam) have been sequenced and characterized, and thus the exact origin and evolution of G3P[6] strains in Asia remain to be elucidated. Here, we sequenced and characterized the full genome of a G3P[6] strain (RVA/Human-wt/JPN/SO1199/2020/G3P[6]) found in a stool sample from a 3-month-old infant admitted with acute gastroenteritis in Japan. On full genomic analysis, strain SO1199 was revealed to have a unique Wa-like genogroup configuration: G3-P[6]-I5-R1-C1-M1-A8-N1-T1-E1-H1. VP6 genotype I5 and NSP1 genotype A8 are commonly found in porcine rotavirus strains. Furthermore, phylogenetic analysis demonstrated that all 11 genes of strain SO1199 were closely related to those of porcine and/or porcine-like human rotaviruses and thus appeared to be of porcine origin. Thus, strain SO1199 was shown to possess a porcine-like genomic backbone and thus is likely to be the result of interspecies transmission of a porcine rotavirus strain. Of note is that all 11 genes of strain SO1199 were phylogenetically located in clusters, distinct from those of the previously identified porcine-like human G3P[6] strains from around the world including Asia, suggesting the occurrence of independent porcine-to-human zoonotic transmission events. To our knowledge, this is the first report on full genome-based characterization of a human G3P[6] strain that has emerged in Japan. Our findings revealed the diversity of unconventional human G3P[6] strains in Asia, and provide important insights into the origin and evolution of G3P[6] strains.

Elucidating the conformational dynamics of histo-blood group antigens and their interactions with the rotavirus spike protein through computational lens

In the present study, we investigated the conformational dynamics of histo-blood group antigens (HBGAs) and their interactions with the VP8* domain of four rotavirus genotypes (P[4], P[6], P[19], and P[11]) utilizing all-atom molecular dynamics simulations in explicit water. Our study revealed distinct changes in the dynamic behavior of the same glycan due to linkage variations. We observed that LNFPI HBGA having a terminal β linkage shows two dominant conformations after complexation, whereas only one was obtained for LNFPI with a terminal α linkage. Interestingly, both variants displayed a single dominant structure in the free state. Similarly, LNT and LNnT show a shift in their dihedral linkage profile between their two terminal monosaccharides because of a change in the linkage from β(1-3) to β(1-4). The molecular mechanics generalized Born surface area (MM/GBSA) calculations yielded the highest binding affinity for LNFPI(β)/P[6] (-13.93 kcal/mol) due to the formation of numerous hydrogen bonds between VP8* and HBGAs. LNnT binds more strongly to P[11] (-12.88 kcal/mol) than LNT (-4.41 kcal/mol), suggesting a single change in the glycan linkage might impact its binding profile significantly. We have also identified critical amino acids and monosaccharides (Gal and GlcNAc) that contributed significantly to the protein-ligand binding through the per-residue decomposition of binding free energy. Moreover, we found that the interaction between the same glycan and different protein receptors within the same rotavirus genogroup influenced the micro-level dynamics of the glycan. Overall, our study helps a deeper understanding of the H-type HBGA and rotavirus spike protein interaction.Communicated by Ramaswamy H. Sarma.

Rotavirus A during the COVID-19 Pandemic in Brazil, 2020-2022: Emergence of G6P[8] Genotype

Rotavirus A (RVA) remains a leading cause of acute gastroenteritis (AGE) hospitalizations in children worldwide. During the COVID-19 pandemic, a reduction in vaccination coverage in Brazil and elsewhere was observed, and some reports have demonstrated a reduction in AGE notifications during the pandemic. This study aims to investigate the diversity and prevalence of RVA genotypes in children and adults presenting with AGE symptoms in Brazil during the COVID-19 pandemic between 2020 and 2022. RVA was screened using RT-qPCR; then, G and P genotypes were characterized using one-step multiplex RT-PCR. A total of 2173 samples were investigated over the three-year period, and we detected RVA in 7.7% of samples (n = 167), being 15.5% in 2020, 0.5% in 2021, and 13.8% in 2022. Higher RVA prevalence was observed in the Northeastern region (19.3%) compared to the Southeastern (6.1%) and Southern regions (5.5%). The most affected age group was children aged between 0 and 6 months old; however, this was not statistically significant. Genotyping and phylogenetic analysis identified the emergence of G6P[8] during the period; moreover, it was detected in 10.6% of samples in 2020 and in 83.5% in 2022. In contrast, the prevalence of G3P[8], the previous dominant genotype, decreased from 72.3% in 2020 to 11.3% in 2022. We also identified unusual strains, such as G3P[9] and G9P[4], being sporadically detected during the period. This is the first report on the molecular epidemiology and surveillance of RVA during the COVID-19 pandemic period in Brazil. Our study provides evidence for the importance of maintaining high and sustainable levels of vaccine coverage to protect against RVA disease. Furthermore, it highlights the need to maintain nationwide surveillance in order to monitor future trends and changes in the epidemiology of RVA in Brazil.

Full genotype constellations analysis of unusual DS-1-like G12P[6] and G6P[8] rotavirus strains detected in Brazil, 2019

Rotavirus A (RVA) is a major cause of acute gastroenteritis (AGE) in children worldwide. We report unusual RVA G12P[6] and G6P[8] strains isolated from fecal samples from Brazilian children hospitalized for AGE. The characterized RVA have genome segments backbone: G12-P[6]/ G6-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2 of DS-1-like genogroup. Our study describes the first identification of G6P[8], a DS-1-like genogroup strain. Nucleotide analysis of VP7 and VP4 genes revealed that all G12 Brazilian strains clustered into the sub-lineages IIIB, mostly associated with P[6] lineage I. Additionally, our G6 lineage I strains were closely related to German G6 genotypes, bound with P[8] lineage III, differing from both vaccine strains. The comparative sequence analysis of our strains with vaccine strains revealed amino acid substitutions located in immunodominant regions of VP7 and VP4 proteins. Continuous monitoring of RVA genotypes is essential to evaluate the impact of vaccination on the dynamic nature of RVA evolution.

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.

Structural basis of P[II] rotavirus evolution and host ranges under selection of histo-blood group antigens

Group A rotaviruses cause severe gastroenteritis in infants and young children worldwide, with P[II] genogroup rotaviruses (RVs) responsible for >90% of global cases. RVs have diverse host ranges in different human and animal populations determined by host histo-blood group antigen (HBGA) receptor polymorphism, but details governing diversity, host ranges, and species barriers remain elusive. In this study, crystal structures of complexes of the major P[II] genogroup P[4] and P[8] genotype RV VP8* receptor-binding domains together with Lewis epitope-containing LNDFH I glycans in combination with VP8* receptor-glycan ligand affinity measurements based on NMR titration experiments revealed the structural basis for RV genotype-specific switching between ββ and βα HBGA receptor-binding sites that determine RV host ranges. The data support the hypothesis that P[II] RV evolution progressed from animals to humans under the selection of type 1 HBGAs guided by stepwise host synthesis of type 1 ABH and Lewis HBGAs. The results help explain disease burden, species barriers, epidemiology, and limited efficacy of current RV vaccines in developing countries. The structural data has the potential to impact the design of future vaccine strategies against RV gastroenteritis.

Relationship Between Rotavirus P[6] Infection in Korean Neonates and Histo-Blood Group Antigen: a Single-Center Study

Background

Rotaviruses are a major cause of pediatric gastroenteritis. The rotavirus P[6] genotype is the most prevalent genotype isolated from Korean neonates but has rarely been reported in other countries. Histo-blood group antigen (HBGA) is known to play an important role in rotavirus infection. We investigated the relationship between rotavirus genotype and HBGA-Lewis blood type in Korean children and explored the reasons for the predominance of rotavirus P[6] strain in Korean neonates.

Methods

Blood and stool samples were collected from 16 rotavirus-infected patients. Rotavirus G (VP7) and P (VP4) genotyping was performed using reverse transcription-PCR and sequencing. Lewis antigen phenotypes (Le^a/Le^b) were tested, and HBGA-Lewis genotype was determined by sequencing the secretor (FUT2) and Lewis (FUT3) genes. Deduced amino acid sequences and three-dimensional structures of the VP8* portion of the rotavirus VP4 protein were analyzed.

Results

All P[6] rotaviruses were isolated from neonates under one month of age, who were negative or weakly positive for the Le^b antigen. However, 10 of the 11 non-P[6] rotaviruses were isolated from older children who were Le^b antigen-positive. The VP8* amino acid sequences differed among P[6], P[4], and P[8] genotypes. Korean P[6] strains showed a unique VP8* sequence with amino acid substitutions, including Y169 > L169, which differed from the sequences of P[6] strains from other countries.

Conclusions

The predominance of the rotavirus P[6] genotype in Korean neonates may be related to the interaction between HBGA-Lewis antigen and the VP8* portion of the VP4 protein, and this information will be helpful in future neonatal vaccine development.

The Rotavirus Vaccine Landscape, an Update

Rotavirus is the leading cause of severe acute childhood gastroenteritis, responsible for more than 128,500 deaths per year, mainly in low-income countries. Although the mortality rate has dropped significantly since the introduction of the first vaccines around 2006, an estimated 83,158 deaths are still preventable. The two main vaccines currently deployed, Rotarix and RotaTeq, both live oral vaccines, have been shown to be less effective in developing countries. In addition, they have been associated with a slight risk of intussusception, and the need for cold chain maintenance limits the accessibility of these vaccines to certain areas, leaving 65% of children worldwide unvaccinated and therefore unprotected. Against this backdrop, here we review the main vaccines under development and the state of the art on potential alternatives.

Rotavirus breakthrough infections responsible for gastroenteritis in vaccinated infants who presented with acute diarrhoea at University Teaching Hospitals, Children's Hospital in 2016, in Lusaka Zambia

Background

In Zambia, before rotavirus vaccine introduction, the virus accounted for about 10 million episodes of diarrhoea, 63 000 hospitalisations and 15 000 deaths in 2015, making diarrhoea the third leading cause of death after pneumonia and malaria. In Zambia, despite the introduction of the vaccine acute diarrhoea due to rotaviruses has continued to affect children aged five years and below. This study aimed to characterise the rotavirus genotypes which were responsible for diarrhoeal infections in vaccinated infants aged 2 to 12 months and to determine the relationship between rotavirus strains and the severity of diarrhoea in 2016.

Methods

Stool samples from infants aged 2 to 12 months who presented to the hospital with acute diarrhoea of three or more episodes in 24 hours were tested for group A rotavirus. All positive specimens that had enough sample were genotyped using reverse transcriptase Polymerase Chain Reaction (RT-PCR). A 20-point Vesikari clinical score between 1–5 was considered as mild, 6–10 as moderate and greater or equal to 11 as severe.

Results

A total of 424 stool specimens were tested of which 153 (36%, 95% CI 31.5% to 40.9%) were positive for VP6 rotavirus antigen. The age-specific rotavirus infections decreased significantly (p = 0.041) from 2–4 months, 32.0% (49/118) followed by a 38.8% (70/181) infection rate in the 5–8 months’ category and subsequently dropped in the infants aged 9–12 months with a positivity rate of 27.2%. 38.5% of infants who received a single dose, 34.5% of those who received a complete dose and 45.2% (19/42) of the unvaccinated tested positive for rotavirus. The predominant rotavirus genotypes included G2P[6] 36%, G1P[8] 32%, mixed infections 19%, G2P[4] 6%, G1P[6] 4% and G9P[6] 3%.

Discussion and conclusion

Results suggest breakthrough infection of heterotypic strains (G2P[6] (36%), homotypic, G1P[8] (32%) and mixed infections (19%) raises concerns about the effects of the vaccination on the rotavirus diversity, considering the selective pressure that rotavirus vaccines could exert on viral populations. This data indicates that the rotavirus vaccine has generally reduced the severity of diarrhoea despite the detection of the virus strains.

FUT2 Secretor Status Influences Susceptibility to VP4 Strain-Specific Rotavirus Infections in South African Children

Gastroenteritis is a preventable cause of morbidity and mortality worldwide. Rotavirus vaccination has significantly reduced the disease burden, but the sub-optimal vaccine efficacy observed in low-income regions needs improvement. Rotavirus VP4 'spike' proteins interact with FUT2-defined, human histo-blood group antigens on mucosal surfaces, potentially influencing strain circulation and the efficacy of P[8]-based rotavirus vaccines. Secretor status was investigated in 500 children <5 years-old hospitalised with diarrhoea, including 250 previously genotyped rotavirus-positive cases (P[8] = 124, P[4] = 86, and P[6] = 40), and 250 rotavirus-negative controls. Secretor status genotyping detected the globally prevalent G428A single nucleotide polymorphism (SNP) and was confirmed by Sanger sequencing in 10% of participants. The proportions of secretors in rotavirus-positive cases (74%) were significantly higher than in the rotavirus-negative controls (58%; p < 0.001). The rotavirus genotypes P[8] and P[4] were observed at significantly higher proportions in secretors (78%) than in non-secretors (22%), contrasting with P[6] genotypes with similar proportions amongst secretors (53%) and non-secretors (47%; p = 0.001). This suggests that rotavirus interacts with secretors and non-secretors in a VP4 strain-specific manner; thus, secretor status may partially influence rotavirus VP4 wild-type circulation and P[8] rotavirus vaccine efficacy. The study detected a mutation (rs1800025) ~50 bp downstream of the G428A SNP that would overestimate non-secretors in African populations when using the TaqMan® SNP Genotyping Assay.

Formulation and preclinical studies with a trivalent rotavirus P2-VP8 subunit vaccine

More effective rotavirus vaccines are essential for preventing extensive diarrheal morbidity and mortality in children under five years of age in low-resource regions. Nonreplicating rotavirus vaccines (NRRV) administered parenterally provide an alternate vaccination method to the current licensed oral vaccine. Live attenuated vaccines and may generate increased efficacy in low-resource settings because the parenteral administration route bypasses some of the challenges associated with oral administration, including differences in intestinal environments. Work described here supports development of a trivalent NRRV vaccine for parenteral administration to avoid complications of the gastrointestinal route. Recombinant VP8* subunit proteins representing some of the most prevalent strains of rotavirus infecting humans - DS-1 (P[4]), 1076 (P[6]), and Wa (P[8]) - were combined with an aluminum adjuvant and the P2 epitope of tetanus toxoid to enhance the immune response to this NRRV antigen. Vaccine formulation development included selection of aluminum hydroxide (Alhydrogel®) as an appropriate adjuvant as well as an optimal buffer to maintain antigen stability and optimize antigen binding to the adjuvant. Characterization assays were used to select the lead vaccine formulation and monitor formulation stability. The NRRV liquid formulation was stable for one year at 2°C to 8°C and four weeks at 37°C. Immunogenicity of the NRRV formulation was evaluated using a guinea pig model, where we demonstrated that the adjuvant provided a 20-fold increase in neutralization titer against a homologous antigen and that the P2-fusion also enhanced the serum neutralizing antibody responses. This vaccine candidate is currently being evaluated in human clinical trials.

Molecular basis of P[II] major human rotavirus VP8* domain recognition of histo-blood group antigens

Initial cell attachment of rotavirus (RV) to specific cell surface glycan receptors, which is the essential first step in RV infection, is mediated by the VP8* domain of the spike protein VP4. Recently, human histo-blood group antigens (HBGAs) have been identified as receptors or attachment factors for human RV strains. RV strains in the P[4] and P[8] genotypes of the P[II] genogroup share common recognition of the Lewis b (Le^b) and H type 1 antigens, however, the molecular basis of receptor recognition by the major human P[8] RVs remains unknown due to lack of experimental structural information. Here, we used nuclear magnetic resonance (NMR) spectroscopy-based titration experiments and NMR-derived high ambiguity driven docking (HADDOCK) methods to elucidate the molecular basis for P[8] VP8* recognition of the Le^b (LNDFH I) and type 1 HBGAs. We also used X-ray crystallography to determine the molecular details underlying P[6] recognition of H type 1 HBGAs. Unlike P[6]/P[19] VP8*s that recognize H type 1 HBGAs in a binding surface composed of an α-helix and a β-sheet, referred as the “βα binding site”, the P[8] and P[4] VP8*s bind Le^b HBGAs in a previously undescribed pocket formed by the edges of two β-sheets, referred to as the “ββ binding site”. Importantly, the P[8] and P[4] VP8*s retain binding capability to non-Le^b type 1 HBGAs using the βα binding site. The presence of two distinct binding sites for Le^b and non-Le^b HBGA glycans in the P[8] and P[4] VP8* domains suggests host-pathogen co-evolution under structural and functional adaptation of RV pathogens to host glycan polymorphisms. Assessment and understanding of the precise impact of this co-evolutionary process in determining RV host ranges and cross-species RV transmission should facilitate improved RV vaccine development and prediction of future RV strain emergence and epidemics.

Rotaviruses (RV)s are the main cause of severe diarrhea in humans and animals. Significant advances in understanding RV diversity, evolution and epidemiology have been made after discovering that RVs recognize histo-blood group antigens (HBGAs) as host cell receptors or attachment factors. While different RV strains are known to have distinct binding preferences for HBGA receptor ligands, their molecular basis in controlling strain-specific host ranges remains unclear. In this study, we used solution nuclear magnetic resonance spectroscopy and X-ray crystallography to determine the molecular-level details for interactions of the human P[8] and P[6] RV VP8* domains with their HBGA receptors ligands. The distinct binding patterns observed between these major human RVs and their respective glycan ligands provide insight into the evolutionary relationships between different P[II] genotypes that ultimately determine host ranges, disease burden, zoonosis and epidemiology, which may impact future strategies for development of vaccines to protect against RV infections.

Evolution of P[8], P[4], and P[6] VP8* genes of human rotaviruses globally reported during 1974 and 2017: possible implications for rotavirus vaccines in development

Non-replicating parenteral rotavirus (RV) vaccine candidates are in development in an attempt to overcome the lower efficacy and effectiveness of oral RV vaccines in low-income countries. One of the leading candidates is a truncated recombinant VP8* protein, expressed in Escherichia coli from original sequences of the prototype RV genotypes P[8], P[4], or P[6] isolated before 1983. Since VP8* is highly variable, it was considered useful to examine the evolutionary changes of RV strains reported worldwide over time in relation to the three P2-VP8 vaccine strains. Here, we retrieved from the GenBank 6,366 RV VP8* gene sequences of P[8], P[4], or P[6] strains isolated between 1974 and 2017, in 77 countries, and compared them with those of the three P2-VP8 vaccine strains: Wa (USA, 1974, G1P[8]), DS-1 (USA, 1976, G2P[4]), and 1076 (Sweden, 1983, G2P[6]). Phylogenetic analysis showed that 94.9% (4,328/4,560), 99.8% (1,141/1,143), and 100% (663/663) of the P[8], P[4], and P[6] strains, respectively, reported globally between 1974 and 2018 belong to non-vaccine lineages. These P[8], P[4], and P[6] RV strains have a mean of 9%, 5%, and 6% amino acid difference from the corresponding vaccine strains. Additionally, in the USA, the mean percentage difference between all the P[8] RV strains and the original Wa strain increased over time: 4% (during 1974-1980), 5% (1988-1991), and 9% (2005-2013). Our analysis substantiated high evolutionary changes in VP8* of the P[8], P[4], and P[6] major RV strains and their increasing variations from the candidate subunit vaccine strains over time. These findings may have implications for the development of new RV vaccines.

Genomic characterization of uncommon human G3P[6] rotavirus strains that have emerged in Kenya after rotavirus vaccine introduction, and pre-vaccine human G8P[4] rotavirus strains

A monovalent rotavirus vaccine (RV1) was introduced to the national immunization program in Kenya in July 2014. There was increased detection of uncommon G3P[6] strains that coincided temporally with the timing of this vaccine introduction. Here, we sequenced and characterized the full genomes of two post-vaccine G3P[6] strains, RVA/Human-wt/KEN/KDH1951/2014/G3P[6] and RVA/Human-wt/KEN/KDH1968/2014/G3P[6], as representatives of these uncommon strains. On full-genomic analysis, both strains exhibited a DS-1-like genotype constellation: G3-P[6]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Phylogenetic analysis revealed that all 11 genes of strains KDH1951 and KDH1968 were very closely related to those of human G3P[6] strains isolated in Uganda in 2012-2013, indicating the derivation of these G3P[6] strains from a common ancestor. Because the uncommon G3P[6] strains that emerged in Kenya are fully heterotypic as to the introduced vaccine strain regarding the genotype constellation, vaccine effectiveness against these G3P[6] strains needs to be closely monitored.

A DS-1 like G9P[6] human strain CDC-6 as a new rotavirus vaccine candidate

Human rotavirus vaccine Rotarix® (G1P[8]) has shown broad cross protection against homotypic and heterotypic Wa-like human rotavirus strains among children worldwide. This vaccine, however, appears to induce slightly less or non-consistent protection against DS-1 like rotavirus P[4] strains in some settings. In addition, children who are secretor or Lewis-negative and are vaccinated with Rotarix® often experience breakthrough infection with P[6] strains. By contrast, P[6] strains infect all children, irrespective of their secretor or Lewis status. In the present study, we report successful adaptation of a DS-1 like human rotavirus G9P[6] strain (CDC-6) to high growth in Vero cells and identify sequence changes that may be critical for enhanced growth in vitro and attenuation in vivo. This human G9P[6] strain could serve as a promising new and potential low-cost vaccine candidate for global use, particularly in targeted population with secretor or Lewis-negative status and high prevalent DS-1 like P[6] strains.

Glycan recognition in globally dominant human rotaviruses

Rotaviruses (RVs) cause life-threatening diarrhea in infants and children worldwide. Recent biochemical and epidemiological studies underscore the importance of histo-blood group antigens (HBGA) as both cell attachment and susceptibility factors for the globally dominant P[4], P[6], and P[8] genotypes of human RVs. How these genotypes interact with HBGA is not known. Here, our crystal structures of P[4] and a neonate-specific P[6] VP8*s alone and in complex with H-type I HBGA reveal a unique glycan binding site that is conserved in the globally dominant genotypes and allows for the binding of ABH HBGAs, consistent with their prevalence. Remarkably, the VP8* of P[6] RVs isolated from neonates displays subtle structural changes in this binding site that may restrict its ability to bind branched glycans. This provides a structural basis for the age-restricted tropism of some P[6] RVs as developmentally regulated unbranched glycans are more abundant in the neonatal gut.

Human rotaviruses (RV) bind to histo-blood group antigens (HBGA) for attachment, but how different viral genotypes interact with HBGA isn’t known. Here, Hu et al. report crystal structures of a prevalent and a neonate-specific RV in complex with HBGA and provide insights into glycan recognition and age-restricted tropism of RVs.

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.