Complete genome constellation of a dominant Bovine rotavirus genotype circulating in Bangladesh reveals NSP4 intragenic recombination with human strains

Rotavirus A is a leading cause of non-bacterial gastroenteritis in humans and domesticated animals. Despite the vast diversity of bovine Rotavirus A strains documented in South Asian countries, there are very few whole genomes available for phylogenetic study. A cross-sectional study identified a high prevalence of the G6P[11] genotype of bovine Rotavirus A circulating in the commercial cattle population in Bangladesh. Next-generation sequencing and downstream phylogenetic analysis unveiled all 11 complete gene segments of this strain (BD_ROTA_CVASU), classifying it under the genomic constellation G6P[11]-I2-R2-C2-M2-A13-N2-T6-E2-H3, which belongs to a classical DS-1-like genomic backbone. We found strong evidence of intragenic recombination between human and bovine strains in the Non-structural protein 4 (NSP4) gene, which encodes a multifunctional enterotoxin. Our analyses highlight frequent zoonotic transmissions of rotaviruses in diverse human-animal interfaces, which might have contributed to the evolution and pathogenesis of this dominant genotype circulating in the commercial cattle population in Bangladesh.

Cross-sectional study and genotyping of rotavirus-A infections in ruminants in Kuwait

Background

Group A rotaviruses (RVA) are zoonotic pathogens responsible for acute enteritis in human and neonatal ruminants. This research aimed to determine the prevalence of RVA in ruminants (cattle, sheep, and goats) and investigate the circulating RVA genotypes in these animals in Kuwait. We conducted a cross-sectional study to detect RVA in ruminants, using an immunochromatography test (IC), direct sandwich ELISA test, and real-time RT-PCR (RT-qPCR) assay using fecal samples.

Results

A total of 400 cattle, 334 sheep, and 222 goats were examined. The prevalence of RVA was 5.3, 1.2, and 2.3%, respectively, using IC. The ELISA test detected RVA from 4.3% of cattle, 0.9% of sheep, and 1.8% of goats. There was a significant association between the occurrence of diarrhea and the presence of RVA in bovine fecal samples (p-value = 0.0022), while no statistical association between diarrhea and the presence of RVA in fecal samples of sheep and goats was observed (p-value = 0.7250; p-value = 0.4499, respectively). Twenty-three of the IC-positive samples (17 from cattle, two from sheep, and four from goats) were tested using a RT-qPCR RVA detection assay targeting the NSP3 gene. The results showed that 21 of 23 IC-positive samples tested positive by RT-qPCR. Detection of RVA genotypes revealed that G10P[11] was the predominant strain in cattle (58.8%), followed by G8P[1] (11.7%). One sheep sample was genotyped as G8P[1]. In addition, G6P[1] and G6P[14] were detected in goat samples.

Conclusion

The present study revealed that the IC was more sensitive in detecting RVA antigen in fecal samples than the ELISA test. A higher occurrence of RVA infection was observed in cattle than in sheep and goats. This study suggests that RVA might be a risk factor of diarrhea in bovine calves less than 2 weeks old. This research also demonstrates the circulation of RVA in sheep and goat populations in Kuwait. Finally, the G10P[11] RVA genotype was the most prevalent genotype identified from cattle samples.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-021-02944-4.

Epidemiology and Molecular Characterization of Rotavirus A in Fruit Bats in Bangladesh

Rotavirus A (RVA) is the primary cause of acute dehydrating diarrhea in human and numerous animal species. Animal-to-human interspecies transmission is one of the evolutionary mechanisms driving rotavirus strain diversity in humans. We screened fresh feces from 416 bats (201 Pteropus medius, 165 Rousettus leschenaultii and 50 Taphozous melanopogon) for RVA using rRT-PCR. We detected a prevalence of 7% (95% CI 3.5-10.8) and 2% (95% CI 0.4-5.2) in P. medius and R. leschenaultii, respectively. We did not detect RVA in the insectivorous bat (T. melanopogon). We identified RVA strains similar to the human strains of G1 and G8 based on sequence-based genotyping, which underscores the importance of including wildlife species in surveillance for zoonotic pathogens to understand pathogen transmission and evolution better.

Report of the 1st African Enteric Viruses Genome Initiative (AEVGI) Data and Bioinformatics Workshop on whole-genome analysis of some African rotavirus strains held in Bloemfontein, South Africa

The University of the Free State - Next Generation Sequencing (NGS) Unit, Bloemfontein, South Africa, hosted a data and bioinformatics workshop from 19 to 22 June 2018. The workshop was coordinated by the African Enteric Viruses Genome Initiative (AEVGI) with support from the Bill & Melinda Gates Foundation. The event introduced technologies in NGS and data analysis with focus on the rotavirus (RV) genome. The workshop fostered interactions and networking between professionals, scientific experts, technicians and students. The courses provided an overview of RV diarrhoea and its burden in Africa, while highlighting the key resources and methodologies in NGS and advanced bioinformatics in deciphering vaccine impact. It was concluded that, despite the reported significant decline in RV associated-diarrhoea mortality and morbidity in Africa due to RV vaccine impact, the need for continuous surveillance and genomic characterization to better understand the ever-changing dynamics of RV strains is imperative.

Whole genomic analysis of G2P[4] human Rotaviruses in Mymensingh, north-central Bangladesh

Rotavirus A (RVA) is a dominant causative agent of acute gastroenteritis in children worldwide. G2P[4] is one of the most common genotypes among human rotavirus (HRV) strains, and has been persistently prevalent in South Asia including Bangladesh. In the present study, whole genome sequences of a total of 16 G2P[4] HRV strains (8 strains each in 2010 and 2013) detected in Mymensingh, north-central Bangladesh were determined. These strains had typical DS-1-like genotype constellation. Most of gene segments from DS-1 genogroup exhibited high level sequence identities to each other (>98%), while slight diversity was observed for VP1, VP3, and NSP4 genes. By phylogenetic analysis, individual RNA segments were classified into one (V) or two-three lineages (V–VI or V–VII). In terms of lineages (sublineages) of 11 gene segments, the 16 Bangladeshi strains could be further classified into four clades (A-D) containing 8 lineage constellations, revealing the presence of three clades (A-C) with three lineage constellations in 2010, and a single clade (D) with four constellations in 2013. Therefore, co-existence of multiple G2P[4] HRV strains with different lineage constellations, and change in clades for the study period were demonstrated. Although amino acids in the antigenic regions on VP7 and VP4 were mostly identical to those of global G2P[4] strains after 2000, VP4 of clade D RVAs in 2013 had alanine and proline at positions 88 and 114, respectively, which are novel substitutions compared with recent global G2P[4] strains. Replacement of lineage constellations associated with unique amino acid changes in the antigenic region in VP4 suggested continuous genetic evolutionary state for emerging new G2P[4] rotavirus strains in Bangladesh.

Genomic constellation and evolution of Ghanaian G2P[4] rotavirus strains from a global perspective

Understanding of the genetic diversity and evolution of Rotavirus A (RVA) strains, a common cause of severe diarrhoea in children, needs to be based on the analysis at the whole genome level in the vaccine era. This study sequenced the whole genomes of six representative G2P[4] strains detected in Ghana from 2008 to 2013, and analysed them phylogenetically with a global collection of G2P[4] strains and African non-G2P[4] DS-1-like strains. The genotype constellation of the study strains was G2-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Strains from the same season were highly identical across the whole genome while strains from different seasons were more divergent from each other. The VP7, VP4, VP2, NSP1, and NSP5 genes belonged to lineage IVa; the VP6, VP1, NSP2, and NSP3 genes belonged to lineage V, and all these genes evolved in the same fashion as the global strains. In the NSP4 gene, lineages V (2008) and X (2009) were replaced by VI (2012/2013) whereas in the VP3 gene, lineage V (2008/2009) was replaced by VII (2012/2013) and these replacements coincided with the vaccine introduction period (2012). The evolutionary rate of the NSP4 gene was 1.2×10^-3 substitutions/site/year and was rather comparable to that of the remaining 10 genes. The multiple NSP4 lineages were explained by intra-genotype reassortment with co-circulating African human DS-1-like strains bearing G2[6], G3P[6], G6[6] and G8. There was no explicit evidence of the contribution of animal RVA strains to the genome of the Ghanaian G2P[4] strains. In summary, this study revealed the dynamic evolution of the G2P[4] strains through intra-genotype reassortment events leading to African specific lineages such IX and X in the NSP4 gene. So far, there was no evidence of a recent direct involvement of animal RVA genes in the genome diversity of African G2P[4] strains.

Emergence and Characterization of Unusual DS-1-Like G1P[8] Rotavirus Strains in Children with Diarrhea in Thailand

The emergence and rapid spread of unusual DS-1-like G1P[8] rotaviruses in Japan have been recently reported. During rotavirus surveillance in Thailand, three DS-1-like G1P[8] strains (RVA/Human-wt/THA/PCB-180/2013/G1P[8], RVA/Human-wt/THA/SKT-109/2013/G1P[8], and RVA/Human-wt/THA/SSKT-41/2013/G1P[8]) were identified in stool specimens from hospitalized children with severe diarrhea. In this study, we sequenced and characterized the complete genomes of strains PCB-180, SKT-109, and SSKT-41. On whole genomic analysis, all three strains exhibited a unique genotype constellation including both genogroup 1 and 2 genes: G1-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2. This novel genotype constellation is shared with Japanese DS-1-like G1P[8] strains. Phylogenetic analysis revealed that the G/P genes of strains PCB-180, SKT-109, and SSKT-41 appeared to have originated from human Wa-like G1P[8] strains. On the other hand, the non-G/P genes of the three strains were assumed to have originated from human DS-1-like strains. Thus, strains PCB-180, SKT-109, and SSKT-41 appeared to be derived through reassortment event(s) between Wa-like G1P[8] and DS-1-like human rotaviruses. Furthermore, strains PCB-180, SKT-109, and SSKT-41 were found to have the 11-segment genome almost indistinguishable from one another in their nucleotide sequences and phylogenetic lineages, indicating the derivation of the three strains from a common origin. Moreover, all the 11 genes of the three strains were closely related to those of Japanese DS-1-like G1P[8] strains. Therefore, DS-1-like G1P[8] strains that have emerged in Thailand and Japan were assumed to have originated from a recent common ancestor. To our knowledge, this is the first report on whole genome-based characterization of DS-1-like G1P[8] strains that have emerged in an area other than Japan. Our observations will provide important insights into the evolutionary dynamics of emerging DS-1-like G1P[8] rotaviruses.

Molecular epidemiology of contemporary G2P[4] human rotaviruses cocirculating in a single U.S. community: footprints of a globally transitioning genotype

Group A rotaviruses (RVs) remain a leading cause of childhood gastroenteritis worldwide. Although the G/P types of locally circulating RVs can vary from year to year and differ depending upon geographical location, those with G1P[8], G2P[4], G3P[8], G4P[8], G9P[8], and G12P[8] specificities typically dominate. Little is known about the evolution and diversity of G2P[4] RVs and the possible role that widespread vaccine use has had on their increased frequency of detection. To address these issues, we analyzed the 12 G2P[4] RV isolates associated with a rise in RV gastroenteritis cases at Vanderbilt University Medical Center (VUMC) during the 2010-2011 winter season. Full-genome sequencing revealed that the isolates had genotype 2 constellations typical of DS-1-like viruses (G2P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2). Phylogenetic analyses showed that the genome segments of the isolates were comprised of two or three different subgenotype alleles; this enabled recognition of three distinct clades of G2P[4] viruses that caused disease at VUMC in the 2010-2011 season. Although the three clades cocirculated in the same community, there was no evidence of interclade reassortment. Bayesian analysis of 328 VP7 genes of G2 viruses isolated in the last 39 years indicate that existing G2 VP7 gene lineages continue to evolve and that novel lineages, as represented by the VUMC isolates, are constantly being formed. Moreover, G2 lineages are characteristically shaped by lineage turnover events that introduce new globally dominant strains every 7 years, on average. The ongoing evolution of G2 VP7 lineages may give rise to antigenic changes that undermine vaccine effectiveness in the long term.

IMPORTANCE

Little is known about the diversity of cocirculating G2 rotaviruses and how their evolution may undermine the effectiveness of rotavirus vaccines. To expand our understanding of the potential genetic range exhibited by rotaviruses circulating in postvaccine communities, we analyzed part of a collection of rotaviruses recovered from pediatric patients in the United States from 2010 to 2011. Examining the genetic makeup of these viruses revealed they represented three segregated groups that did not exchange genetic material. The distinction between these three groups may be explained by three separate introductions. By comparing a specific gene, namely, VP7, of the recent rotavirus isolates to those from a collection recovered from U.S. children between 1974 and 1991 and other globally circulating rotaviruses, we were able to reconstruct the timing of events that shaped their ancestry. This analysis indicates that G2 rotaviruses are continuously evolving, accumulating changes in their genetic material as they infect new patients.

G8 rotaviruses with conserved genotype constellations detected in Malawi over 10 years (1997-2007) display frequent gene reassortment among strains co-circulating in humans

Rotavirus A, the most common cause of severe diarrhoea in children worldwide, occurs in five major VP7 (G) and VP4 (P) genotype combinations, comprising G1P[8], G2P[4], G3P[8], G4P[8] and G9P[8]. However, G8, a common bovine rotavirus genotype, has been reported frequently among children in African countries. Surveillance of rotavirus gastroenteritis conducted in a sentinel hospital in Blantyre, Malawi between 1997 and 2007 provided a rare opportunity to examine the whole genotype constellation of G8 strains and their evolution over time. A sample of 27 (9.0 %) of 299 G8 strains was selected to represent each surveillance year and a range of P genotypes, which shifted in predominance from P[6] to P[4] and P[8] during the study period. Following cell culture adaptation, whole genome sequencing demonstrated that the genetic background of 26 strains possessed the DS-1 genotype constellation. A single G8P[6] strain was a reassortant in which both NSP2 and NSP5 genes from strains with the Wa genotype constellation had been inserted into a strain with the DS-1 genotype background. Phylogenetic analysis suggested frequent reassortment among co-circulating strains with the DS-1 genotype constellation. Little evidence was identified to suggest the introduction of contemporary bovine rotavirus genes into any of the 27 G8 strains examined. In conclusion, Malawian G8 strains are closely related to other human strains with the DS-1 genotype constellation. They have evolved over the last decade through genetic reassortment with other human rotaviruses, changing their VP4 genotypes while maintaining a conserved genotype constellation for the remaining structural and non-structural proteins.

Genotype constellation and evolution of group A rotaviruses infecting humans

Numerous rotavirus group A (RVA) strains with distinct G-genotype and P-genotype combinations have been described infecting humans worldwide. However, the increasing amount of complete RVA genome data which have become available, suggest that only RVA strains with 2 discrete genotype constellations have been successful in sustaining infection of humans worldwide over longer periods of time. Those genotype constellations have been designated I1-R1-C1-M1-A1-N1-T1-E1-H1 and I2-R2-C2-M2-A2-N2-T2-E2-H2 and are also known as Wa-like and DS-1-like, respectively. RVAs of other genotype constellations which were able to spread to a limited extent in the human population are AU-1-related RVA strains (I3-R3-C3-M3-A3/A12-N3-T3-E3-H3/H6) in combination with G3P[9] or G12P[9], and neonatal G10P[11] RVA strains in India (bovine×human Wa-like reassortants). On the basis of the analysis of complete genomes, it is suggested that the overall genetic diversity of epidemiologically widespread human RVA strains is more limited than generally assumed. This conclusion has consequences for how we look at host range restriction and the criteria according to which the effectiveness of RVA universal mass vaccination programs is assessed.

Molecular characterization of the NSP4 gene of human group A rotavirus strains circulating in Tunisia from 2006 to 2008

Non-structural protein 4 (NSP4), encoded by group A rotavirus (RVA) genome segment 10, is a multifunctional protein and the first recognized virus-encoded enterotoxin. Recently, a new classification system for RVAs was proposed and a total of 14 NSP4 genotypes (E1-E14) are currently described. The most common NSP4 genotypes in humans are Wa-like E1 and DS-1-like E2. This report represents the first investigation on the genetic diversity of RVA NSP4 genes in Tunisia from 2006 to 2008. In the present study, the NSP4-encoding genes of human RVA strains with different G/P-genotype combinations were analyzed. NSP4 genes of 261 RVA-positive fecal samples were analyzed using a semi-nested reverse transcriptase-polymerase chain reaction and in addition the NSP4 gene of 46 representative RVA strains were sequenced. Phylogenetic analysis of the Tunisian NSP4 nucleotide sequences revealed the presence of two NSP4 genotypes. Genotype E1 was found to be associated with G1P[8], G3P[6], G3P[8], G4P[6] and G4P[8], whereas genotype E2 was associated with G2P[4], G2P[6] and G6P[9] samples. These results support the hypothesis that P[8] carrying RVA strains usually possess the E1 genotype, whereas P[4] carrying RVA strains usually possess the E2 genotype. P[6] carrying strains were found with both E1 and E2. The unusual G6P[9] strains possessed a E2 genotype with a possible animal origin. These results underline the need for further investigations to assess the validity of NSP4 as a suitable target for epidemiologic surveillance of RVA infections and vaccine development.

Whole-genomic analysis of rotavirus strains: current status and future prospects

Studies on genetic diversity of rotaviruses have been primarily based on the genes encoding the antigenically significant VP7 and VP4 proteins. Since the rotavirus genome has 11 segments of RNA that are vulnerable to reassortment events, analyses of the VP7 and VP4 genes may not be sufficient to obtain conclusive data on the overall genetic diversity, or true origin of strains. In the last few years following the advent of the whole-genome-based genotype classification system, the whole genomes of at least 167 human group A rotavirus strains have been analyzed, providing a plethora of new and important information on the complex origin of strains, inter- and intra-genogroup reassortment events, animal–human reassortment events, zoonosis, and genetic linkages involving different group A rotavirus gene segments. In addition, the whole genomes of a limited number of human group B, C and novel group rotavirus strains have been analyzed. This article briefly reviews the available data on whole-genomic analysis of human rotavirus strains. The significance and future prospects of whole-genome-based studies are also discussed.

Whole-genome analysis reveals the complex evolutionary dynamics of Kenyan G2P[4] human rotavirus strains

Although G2P[4] rotaviruses are common causes of acute childhood diarrhoea in Africa, to date there are no reports on whole genomic analysis of African G2P[4] strains. In this study, the nearly complete genome sequences of two Kenyan G2P[4] strains, AK26 and D205, detected in 1982 and 1989, respectively, were analysed. Strain D205 exhibited a DS-1-like genotype constellation, whilst strain AK26 appeared to be an intergenogroup reassortant with a Wa-like NSP2 genotype on the DS-1-like genotype constellation. The VP2-4, VP6-7, NSP1, NSP3 and NSP5 genes of strain AK26 and the VP2, VP4, VP7 and NSP1–5 genes of strain D205 were closely related to those of the prototype or other human G2P[4] strains. In contrast, their remaining genes were distantly related, and, except for NSP2 of AK26, appeared to originate from or share a common origin with rotavirus genes of artiodactyl (ruminant and camelid) origin. These observations highlight the complex evolutionary dynamics of African G2P[4] rotaviruses.

Uniformity of rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group (RCWG)

In April 2008, a nucleotide-sequence-based, complete genome classification system was developed for group A rotaviruses (RVs). This system assigns a specific genotype to each of the 11 genome segments of a particular RV strain according to established nucleotide percent cutoff values. Using this approach, the genome of individual RV strains are given the complete descriptor of Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx. The Rotavirus Classification Working Group (RCWG) was formed by scientists in the field to maintain, evaluate and develop the RV genotype classification system, in particular to aid in the designation of new genotypes. Since its conception, the group has ratified 51 new genotypes: as of April 2011, new genotypes for VP7 (G20-G27), VP4 (P[28]-P[35]), VP6 (I12-I16), VP1 (R5-R9), VP2 (C6-C9), VP3 (M7-M8), NSP1 (A15-A16), NSP2 (N6-N9), NSP3 (T8-T12), NSP4 (E12-E14) and NSP5/6 (H7-H11) have been defined for RV strains recovered from humans, cows, pigs, horses, mice, South American camelids (guanaco), chickens, turkeys, pheasants, bats and a sugar glider. With increasing numbers of complete RV genome sequences becoming available, a standardized RV strain nomenclature system is needed, and the RCWG proposes that individual RV strains are named as follows: RV group/species of origin/country of identification/common name/year of identification/G- and P-type. In collaboration with the National Center for Biotechnology Information (NCBI), the RCWG is also working on developing a RV-specific resource for the deposition of nucleotide sequences. This resource will provide useful information regarding RV strains, including, but not limited to, the individual gene genotypes and epidemiological and clinical information. Together, the proposed nomenclature system and the NCBI RV resource will offer highly useful tools for investigators to search for, retrieve, and analyze the ever-growing volume of RV genomic data.