The first identification of rare human group A rotavirus strain G3P[10] with severe infantile diarrhea in eastern India

Diverse human and bat-like rotavirus G3 strains circulating in suburban Bangkok

Although rotavirus vaccines are available in many parts of the world and are effective in reducing the overall incidence of rotavirus infection, it remains a major cause of diarrhea in less-developed countries. Among various rotavirus group A (RVA) strains, the increasingly common genotype G3 (defined by the VP7 gene) has been identified in both humans and animals. Our previous epidemiological surveillance in Bangkok found several unusual non-vaccine-like G3 strains in patients with diarrhea. In this study, we sequenced and characterized the genomes of seven of these G3 strains, which formed combinations with genotypes P[4], P[6], P[9], and P[10] (defined by the VP4 gene). Interestingly, we identified a bat-like RVA strain with the genome constellation G3-P[10]-I3-R3-C3-M3-A9-N3-T3-E3-H6, which has not been previously reported in the literature. The amino acid residues deduced from the nucleotide sequences of our G3 strains differed at the antigenic epitopes to those of the VP7 capsid protein of the G3 strain in RotaTeq vaccine. Although it is not unusual for the segmented genomes of RVA to reassort and give rise to emerging novel strains, the atypical G3 strains identified in this study suggest possible animal-to-human RVA zoonotic spillover even in urban areas.

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.

High prevalence and genomic characteristics of G6P[1] Bovine Rotavirus A in yak in China

Yak is an iconic species of the Qinghai-Tibet Plateau, which is the world's highest plateau. Here, a total of 541 yak diarrhoeic samples were collected from 69 farms in four provinces in the Qinghai-Tibet Plateau from April 2015 to June 2018, and 73.6 % of samples were detected as Bovine Rotavirus A (BRVA) positive by RT-PCR assay. Two G genotypes (G6, G10) and two P genotypes (P[1], P[11]) were determined, in which G6P[1] BRVA was the predominant strain. Moreover, VP7 and VP4 of these G6P[1] strains showed unique amino acid mutations, such that they clustered into an independent branch in the phylogenetic tree. A strain of BRVA designated as RVA/Yak-tc/CHN/QH-1/2015/G6P[1] was isolated successfully using MA104 cells, and the virus titre was determined as 10^5.84 TCID₅₀ ml^-1. The genome of strain QH-1 had a G6-P[1]-I2-R2-C2-M2-A3-N3-T6-E2-H3 genotype constellation. QH-1 was identified as a reassortment strain of BRVA, human RVA and ovine RVA based on the nucleotide identity and phylogenetic tree of 11 gene segments, indicating its public health significance. To the best of our knowledge, this is the first report on the molecular prevalence and genome characteristics of BRVA in yak, contributing to further understanding of the epidemic and genetic evolution of BRVA.

Genomic analysis reveals G3P[13] porcine rotavirus A interspecific transmission to human from pigs in a swine farm with diarrhoea outbreak

Rotavirus A (RVA) is a major diarrhoea-causing pathogen in young animals and children. The zoonotic potential of RVA has received extensive attention in recent years. In May 2018, an outbreak of diarrhoea among piglets occurred on a swine farm in Sichuan province, PR China. RVA was detected in 95.7 % (22/23) of piglet samples, 60 % (9/15) of sow samples and 100 % (3/3) of pig-breeder faecal samples. The predominant RVA genotype on this swine farm was G3P[13], and G3P[13] RVA was also detected in the three breeder faecal samples. Three G3P[13] RVA strains were isolated from a piglet faecal sample, a sow faecal sample and a pig-breeder faecal sample, and were named SCLS-X1, SCLS-3 and SCLS-R3, respectively. The complete sequences of 11 gene segments of these three isolates were derived. Phylogenetic analysis showed that ten gene segments (VP7, VP4, VP1-VP3 and NSP1-NSP5) of pig-breeder isolate SCLS-R3 were closely related to pig isolates SCLS-X1 and SCLS-3 from this farm. Only the VP6 gene shared higher homology with human RVA strain I321. Therefore, a G3P[13] porcine RVA strain most likely infected pig breeders. These results provided the first complete epidemiological link demonstrating interspecies transmission of G3P[13] RVA from pigs to human. Our data contribute to an improved understanding of the genetic evolution and interspecies transmission of RVA.

Rotavirus G9P[4], G9P[6] and G1P[6] strains isolated from children with acute gastroenteritis in Pune, western India, 2013-2015: evidence for recombination in genes encoding VP3, VP4 and NSP1

Species A rotaviruses (RVAs) are genetically diverse pathogens. These are the most evolutionarily adaptable organisms, with a multitude of mechanisms for evolutionary change. To date, full-genome classification has been proved to be an excellent tool for studying the evolution of unusual rotavirus strains. As limited data are available from Pune (Maharashtra), western India, the current study was undertaken with the aim of understanding the genetic diversity in three (G1P[6], G9P[4] and G9P[4]) unusual RVA strains circulating in Pune, India during 2013-2015. Full-genome analysis of these strains classified them as G1-P[6]-I1-R1-C1-M1-A1-N1-T1-E1-H1, G9-P[4]-I2-R2-C2-[M1-M2_R]-[A1-A2_R]-N2-T2-E6-H2 and G9-[P4-P6_R]-I1-R1-C1-M1-A1-N1-T1-E1-H1. Sequencing and phylogenetic analysis of the structural and non-structural genes of these unusual RVA strains showed nucleotide/amino acid identities of 82.3-98.5 %/77.3-99.8 % and 86.6-97.6 %/89.6-97.8 % between the strains of the study. Evidence of recombination events was found within the genes encoding VP3, VP4 and NSP1, which showed a combination of genetic information for genogroup 1 [M1/P[6]/A1] and genogroup 2 [M2/P[4]/A2] strains. This study will facilitate future investigations into the molecular pathogenesis of such RVAs as the exchange of whole or partial genetic material between rotaviruses through recombination contributes directly to their diversification, adaptation and evolution.

Glycan binding patterns of human rotavirus P[10] VP8* protein

Background

Rotaviruses (RVs) are a major cause of acute children gastroenteritis. The rotavirus P [10] belongs to P[I] genogroup of group A rotaviruses that mainly infect animals, while the rotavirus P [10] was mainly identified from human infection. The rotavirus P [10] is an unusual genotype and the recognition pattern of cellular receptors remains unclear.

Methods

We expressed and purified the RV P [10] VP8* protein and investigated the saliva and oligosaccharide binding profiles of the protein. A homology model of the P [10] VP8* core protein was built and the superimposition structural analysis of P [10] VP8* protein on P [19] VP8* in complex with mucin core 2 was performed to explore the possible docking structural basis of P [10] VP8* and mucin cores.

Results

Our data showed that rotavirus P [10] VP8* protein bound to all ABO secretor and non-secretor saliva. The rotavirus P [10] could bind strongly to mucin core 2 and weakly to mucin core 4. The homology modeling indicated that RV P [10] VP8* binds to mucin core 2 using a potential glycan binding site that is the same to P [19] VP8* belonging to P[II] genogroup.

Conclusion

Our results suggested an interaction of rotavirus P [10] VP8* protein with mucin core 2 and mucin core 4. These findings offer potential for elucidating the mechanism of RV A host specificity, evolution and epidemiology.

Upsurge and spread of G3 rotaviruses in Eastern India (2014-2016): Full genome analyses reveals heterogeneity within Wa-like genomic constellation

Advent of new strains and shift in predominantly circulating genotypes are characteristics of group- A rotavirus (RVA), one of the major causes of childhood gastroenteritis. During diarrheal disease surveillance at Kolkata, India (2014-2016), a shift in circulating RVA strains from G1P[8] to G3P[8] was seen. Stool samples from children (n = 3048) with acute gastroenteritis were tested of which 38.7% were RVA positive. G1 was the predominant strain (65.3%) in 2014-2015 whereas in late 2015 and 2016, G3 became the preponderant strain (44.6%). In the past decade G3 strains were not observed in this region, we conducted whole genome sequencing of representative strains to gain insight into the phenomenon of emergence and genetic constellation of these circulating human G3 strains. The analyses revealed intergenogroup reassortment in G3P[4] strains (among Wa and DS-1-like genogroup) whereas G3P[8] strains were authentic Wa-like. Phylogenetic analysis revealed Kolkata G3 strains as polymorphic and thus they formed two sub-clusters due to antigenic differences in their VP7 protein. One of the sub-clusters had the wild-type threonine at 87 amino acid position while another sub-cluster had an isoleucine mutation. Presence of additional N-linked glycosylation site at amino acid 283 of VP7 glycoprotein suggests that the major neutralizing epitope on the VP7 (G3) of RotaTeq vaccine differs from the currently circulating G3 strains. The study is important as efficiency of rotavirus vaccine depends on the circulating heterogeneous genotype constellations. Continuous monitoring of circulating RVA strains in endemic settings like India is therefore important in pre- and post-vaccination period to monitor the emergence of new reassortant genotypes in addition to assessing vaccine efficacy.

Structural basis of glycan specificity of P[19] VP8*: Implications for rotavirus zoonosis and evolution

Recognition of specific cell surface glycans, mediated by the VP8* domain of the spike protein VP4, is the essential first step in rotavirus (RV) infection. Due to lack of direct structural information of virus-ligand interactions, the molecular basis of ligand-controlled host ranges of the major human RVs (P[8] and P[4]) in P[II] genogroup remains unknown. Here, through characterization of a minor P[II] RV (P[19]) that can infect both animals (pigs) and humans, we made an important advance to fill this knowledge gap by solving the crystal structures of the P[19] VP8* in complex with its ligands. Our data showed that P[19] RVs use a novel binding site that differs from the known ones of other genotypes/genogroups. This binding site is capable of interacting with two types of glycans, the mucin core and type 1 histo-blood group antigens (HBGAs) with a common GlcNAc as the central binding saccharide. The binding site is apparently shared by other P[II] RVs and possibly two genotypes (P[10] and P[12]) in P[I] as shown by their highly conserved GlcNAc-interacting residues. These data provide strong evidence of evolutionary connections among these human and animal RVs, pointing to a common ancestor in P[I] with a possible animal host origin. While the binding properties to GlcNAc-containing saccharides are maintained, changes in binding to additional residues, such as those in the polymorphic type 1 HBGAs may occur in the course of RV evolution, explaining the complex P[II] genogroup that mainly causes diseases in humans but also in some animals.

Complete genotyping of unusual species A rotavirus G12P[11] and G10P[14] isolates and evidence of frequent in vivo reassortment among the rotaviruses detected in children with diarrhea in Kolkata, India, during 2014

Species A rotaviruses (RVA) are the most important cause of acute gastroenteritis in the young of humans and many animal species globally. G1P[8], G2P[4], G3P[8], G4P[8], G9P[6/8] and G12P[6/8] are the predominantly isolated genotypes throughout the world including India. Unusual genotypes from different host species such as G5, G6, G8, G10 and G11 have also been reported in humans with low frequency. In the present study, among >650 RVA positive stool samples collected from children with diarrhea in Kolkata, India, during 2014, two isolates each of the genotype G12P[11] and G10P[14] were obtained and their genomes completely sequenced. The full genotype constellations were G12-P[11]-I1-R1-C1-M2-A1-N1-T2-E1-H1 and G12-P[11]-I1-R1-C1-M1-A5-N1-T1-E1-H1 for G12P[11] viruses, suggesting several reassortments between Wa- and DS-1-like human RVA strains, including possible reassortment of a simian NSP1 gene. The G10P[14] viruses (G10-P[14]-I2-R2-C2-M2-A11-N2-T6-E2-H3) were found to contain multiple genes closely related to RVAs of artiodactyl origin, highlighting the role of inter-host species transmissions of RVAs. From the G/P constellation of all RVA isolates, it could be concluded that approximately one quarter had likely arisen from reassortment events in vivo among RVAs of 'usual' genotypes.

Genetic diversity of G3 rotavirus strains circulating in Argentina during 1998-2012 assessed by full genome analyses

Seasonal shifts in the predominant strains and the periodic emergence of new strains are epidemiological features of human rotaviruses. After the sporadic detection in two samples in 1998, G3P[8] strains reemerged as the predominant rotavirus during 2008-2009 in Argentina. Notably, in 2011 6.3% (37/587) of samples presented the G3P[6] genotypes, which coincided with the recent detection of G3P[6] and G2P[6] strains in South America and Europe. Analyses of the 11 gene segments of four G3P[8] and two G3P[6] strains revealed that G3P[8] strains detected a decade apart (1998 and 2009) presented minor differences, while the G3P[6] strains presented a complete different genomic constellation albeit showing a similar VP7 gene. This study provides insights in the dynamics and evolution of one of the genotypes with the wider range of hosts and inter-species transmission potential.

Review of global rotavirus strain prevalence data from six years post vaccine licensure surveillance: is there evidence of strain selection from vaccine pressure?

Comprehensive reviews of pre licensure rotavirus strain prevalence data indicated the global importance of six rotavirus genotypes, G1P[8], G2P[4], G3P[8], G4P[8], G9P[8] and G12P[8]. Since 2006, two vaccines, the monovalent Rotarix (RV1) and the pentavalent RotaTeq (RV5) have been available in over 100 countries worldwide. Of these, 60 countries have already introduced either RV1 or RV5 in their national immunization programs. Post licensure vaccine effectiveness is closely monitored worldwide. This review aimed at describing the global changes in rotavirus strain prevalence over time. The genotype distribution of the nearly 47,000 strains that were characterized during 2007-2012 showed similar picture to that seen in the preceding period. An intriguing finding was the transient predominance of heterotypic strains, mainly in countries using RV1. Unusual and novel antigen combinations continue to emerge, including some causing local outbreaks, even in vaccinated populations. In addition, vaccine strains have been found in both vaccinated infants and their contacts and there is evidence for genetic interaction between vaccine and wild-type strains. In conclusion, the post-vaccine introduction strain prevalence data do not show any consistent pattern indicative of selection pressure resulting from vaccine use, although the increased detection rate of heterotypic G2P[4] strains in some countries following RV1 vaccination is unusual and this issue requires further monitoring.

Analysis of human rotaviruses from a single location over an 18-year time span suggests that protein coadaption influences gene constellations

Rotaviruses (RVs) are 11-segmented, double-stranded RNA viruses that cause severe gastroenteritis in children. In addition to an error-prone genome replication mechanism, RVs can increase their genetic diversity by reassorting genes during host coinfection. Such exchanges allow RVs to acquire advantageous genes and adapt in the face of selective pressures. However, reassortment may also impose fitness costs if it unlinks genes/proteins that have accumulated compensatory, coadaptive mutations and that operate best when kept together. To better understand human RV evolutionary dynamics, we analyzed the genome sequences of 135 strains (genotype G1/G3/G4-P[8]-I1-C1-R1-A1-N1-T1-E1-H1) that were collected at a single location in Washington, DC, during the years 1974 to 1991. Intragenotypic phylogenetic trees were constructed for each viral gene using the nucleotide sequences, thereby defining novel allele level gene constellations (GCs) and illuminating putative reassortment events. The results showed that RVs with distinct GCs cocirculated during the vast majority of the collection years and that some of these GCs persisted in the community unchanged by reassortment. To investigate the influence of protein coadaptation on GC maintenance, we performed a mutual information-based analysis of the concatenated amino acid sequences and identified an extensive covariance network. Unexpectedly, amino acid covariation was highest between VP4 and VP2, which are structural components of the RV virion that are not thought to directly interact. These results suggest that GCs may be influenced by the selective constraints placed on functionally coadapted, albeit noninteracting, viral proteins. This work raises important questions about mutation-reassortment interplay and its impact on human RV evolution.

IMPORTANCE

Rotaviruses are devastating human pathogens that cause severe diarrhea and kill >450,000 children each year. The virus can evolve by accumulating mutations and by acquiring new genes from other strains via a process called reassortment. However, little is known about the relationship between mutation accumulation and gene reassortment for rotaviruses and how it impacts viral evolution. In this study, we analyzed the genome sequences of human strains found in clinical fecal specimens that were collected at a single hospital over an 18-year time span. We found that many rotaviruses did not reassort their genes but instead maintained them as specific sets (i.e., constellations). By analyzing the encoded proteins, we discovered concurrent amino acid changes among them, which suggests that they are functionally coadapted to operate best when kept together. This study increases our understanding of how rotaviruses evolve over time in the human population.

Rotaviruses: is their surveillance needed?

Rotaviruses, a major cause of gastroenteritis in children worldwide accounts for around 0.5 million deaths annually. Owing to their segmented genome and frequently evolving capability, these display a wide variation in their genotypes. In addition to commonly circulating genotypes (G1, G2, G3, G4, G9, P[4] and P[8]), a number of infrequent genotypes are being continuously reported to infect humans. These viral strains exhibit variation from one geographical setting to another in their distribution. Though the introduction of vaccines (RotaTeq and Rotarix) proved to be very effective in declining rotavirus associated morbidity and mortality, the number of infections remained same. Unusual genotypes significantly contribute to the rotavirus associated diarrhoeal burden, may reduce the efficacy of the vaccines in use and hence vaccinated individuals may not be benefited. Vaccine introduction may bring about a notable impact on the distribution and prevalence of these viruses due to selection pressure. Moreover, there is a sudden emergence of G2 and G3 in Brazil and United States, respectively, during the years 2006-2008 post-vaccination introduction; G9 and G12 became predominant during the years 1986 through 1998 before the vaccine introduction and now are commonly prevalent strains; and disparity in the predominance of strains after introduction of vaccines and their natural fluctuations poses a vital question on the impact of vaccines on rotavirus strain circulation. This interplay between vaccines and rotavirus strains is yet to be explored, but it certainly enforces the need to continuously monitor these changes in strains prevalence in a particular region. Furthermore, these fluctuations should be considered while administration or development of a vaccine, if rotavirus associated mortality is ever to be controlled.

First report of human rotavirus G8P[4] gastroenteritis in India: evidence of ruminants-to-human zoonotic transmission

Group A rotaviruses are the major cause of childhood gastroenteritis worldwide. Due to close proximity of human and cattle in rural areas of developing countries like India, interspecies transmission or zoonotic transmission is a major source of rapid generation of reassortants and genetic or antigenic variants. Previously, many human group A rotaviruses were found with porcine or bovine characteristics from eastern and north-eastern India. In this study, four unusual human G8P[4] strains were identified which had artiodactyl-like origins. During an ongoing community based surveillance for epidemiological profiling of diarrheal pathogens, these unusual human group A rotavirus G8P[4] strains were detected from the stool samples of 3-14 months old children with acute diarrhea in Sonarpur, eastern India. Analysis of eleven complete and/or partial gene segments of these unusual G8P[4] strains were done by reverse transcription-PCR (RT-PCR) followed by nucleotide sequencing and phylogenetic analysis. The VP7 nucleotide sequences revealed a close phylogenetic relationship to the G8P[7] porcine strain D-1 and bovine strain KJ59-2 from South Korea. Whereas the VP4 gene segments were also related closely to human rotavirus prototype strain DS-1. Other nine gene segments of these G8P[4] rotaviruses were related closely to either animal or animal-derived rotavirus members of the DS-1-like family. These results suggest that origin of these G8P[4] strains might have been resulted from multiple reassortment events between artiodactyls and ruminant-derived reassortant human rotaviruses. To date, this is the first report of G8P[4] rotavirus from India and the first genomic analysis of G8P[4] strains from Asian continent.