Molecular characterization of a human G20P[28] rotavirus a strain with multiple genes related to bat rotaviruses

Role of Brazilian bats in the epidemiological cycle of potentially zoonotic pathogens

Bats (Chiroptera) are flying mammals of great biodiversity and habits. These characteristics contribute for them being natural reservoirs and part of the epidemiological cycle of several potentially zoonotic pathogens, such as viruses, protozoa, fungi and bacteria. Brazil hosts approximately 15% of the world's bat diversity, with 181 distinct species, 68 genera and 9 families. About 60% of infectious diseases in humans are of zoonotic origin and, in the last decades, the detection of zoonotic pathogens in bats and their environment has been reported, such as Rabies virus (RABV) and Histoplasma capsulatum. Thus, the aim of this work was to review the reports of zoonotic pathogens associated with bats in Brazil in the past ten years. We reviewed the main pathogenic microorganisms described and the species of bats most frequently involved in the epidemiological cycles of these zoonotic agents. The obtained data show an upward trend in the detection of zoonotic pathogens in Brazilian bats, such as RABV, Bartonella sp., Histoplasma capsulatum and Leishmania spp., with emphasis on the bat species Artibeus lituratus, Carollia perspicillata, Desmodus rotundus and Molossus molossus. These findings highlight the importance of monitoring bat-associated microrganisms to early identify pathogens that may threaten bat populations, including potentially zoonotic microrganisms, emphasizing the importance of the One Health approach to prevent and mitigate the risks of the emergence of zoonotic diseases.

Emerging viruses: Cross-species transmission of coronaviruses, filoviruses, henipaviruses, and rotaviruses from bats

Emerging infectious diseases, especially if caused by bat-borne viruses, significantly affect public health and the global economy. There is an urgent need to understand the mechanism of interspecies transmission, particularly to humans. Viral genetics; host factors, including polymorphisms in the receptors; and ecological, environmental, and population dynamics are major parameters to consider. Here, we describe the taxonomy, geographic distribution, and unique traits of bats associated with their importance as virus reservoirs. Then, we summarize the origin, intermediate hosts, and the current understanding of interspecies transmission of Middle East respiratory syndrome coronavirus (MERS-CoV), severe acute respiratory syndrome coronavirus (SARS-CoV), SARS-CoV-2, Nipah, Hendra, Ebola, Marburg virus, and rotaviruses. Finally, the molecular interactions of viral surface proteins with host cell receptors are examined, and a comparison of these interactions in humans, intermediate hosts, and bats is conducted. This uncovers adaptive mutations in virus spike protein that facilitate cross-species transmission and risk factors associated with the emergence of novel viruses from bats.

Whole genome analysis of rotavirus strains circulating in Benin before vaccine introduction, 2016-2018

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

At Least Seven Distinct Rotavirus Genotype Constellations in Bats with Evidence of Reassortment and Zoonotic Transmissions

The increased research on bat coronaviruses after severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) allowed the very rapid identification of SARS-CoV-2. This is an excellent example of the importance of knowing viruses harbored by wildlife in general, and bats in particular, for global preparedness against emerging viral pathogens.

Bats host many viruses pathogenic to humans, and increasing evidence suggests that rotavirus A (RVA) also belongs to this list. Rotaviruses cause diarrheal disease in many mammals and birds, and their segmented genomes allow them to reassort and increase their genetic diversity. Eighteen out of 2,142 bat fecal samples (0.8%) collected from Europe, Central America, and Africa were PCR-positive for RVA, and 11 of those were fully characterized using viral metagenomics. Upon contrasting their genomes with publicly available data, at least 7 distinct bat RVA genotype constellations (GCs) were identified, which included evidence of reassortments and 6 novel genotypes. Some of these constellations are spread across the world, whereas others appear to be geographically restricted. Our analyses also suggest that several unusual human and equine RVA strains might be of bat RVA origin, based on their phylogenetic clustering, despite various levels of nucleotide sequence identities between them. Although SA11 is one of the most widely used reference strains for RVA research and forms the backbone of a reverse genetics system, its origin remained enigmatic. Remarkably, the majority of the genotypes of SA11-like strains were shared with Gabonese bat RVAs, suggesting a potential common origin. Overall, our findings suggest an underexplored genetic diversity of RVAs in bats, which is likely only the tip of the iceberg. Increasing contact between humans and bat wildlife will further increase the zoonosis risk, which warrants closer attention to these viruses.

Comparative genomic analysis of genogroup 1 and genogroup 2 rotaviruses circulating in seven US cities, 2014-2016

For over a decade, the New Vaccine Surveillance Network (NVSN) has conducted active rotavirus (RVA) strain surveillance in the USA. The evolution of RVA in the post-vaccine introduction era and the possible effects of vaccine pressure on contemporary circulating strains in the USA are still under investigation. Here, we report the whole-gene characterization (eleven ORFs) for 157 RVA strains collected at seven NVSN sites during the 2014 through 2016 seasons. The sequenced strains included 52 G1P[8], 47 G12P[8], 18 G9P[8], 24 G2P[4], 5 G3P[6], as well as 7 vaccine strains, a single mixed strain (G9G12P[8]), and 3 less common strains. The majority of the single and mixed strains possessed a Wa-like backbone with consensus genotype constellation of G1/G3/G9/G12-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1, while the G2P[4], G3P[6], and G2P[8] strains displayed a DS-1-like genetic backbone with consensus constellation of G2/G3-P[4]/P[6]/P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Two intergenogroup reassortant G1P[8] strains were detected that appear to be progenies of reassortment events between Wa-like G1P[8] and DS-1-like G2P[4] strains. Two Rotarix^® vaccine (RV1) and two RV5 derived (vd) reassortant strains were detected. Phylogenetic and similarity matrices analysis revealed 2-11 sub-genotypic allelic clusters among the genes of Wa- and DS-1-like strains. Most study strains clustered into previously defined alleles. Amino acid (AA) substitutions occurring in the neutralization epitopes of the VP7 and VP4 proteins characterized in this study were mostly neutral in nature, suggesting that these RVA proteins were possibly under strong negative or purifying selection in order to maintain competent and actual functionality, but fourteen radical (AA changes that occur between groups) AA substitutions were noted that may allow RVA strains to gain a selective advantage through immune escape. The tracking of RVA strains at the sub-genotypic allele constellation level will enhance our understanding of RVA evolution under vaccine pressure, help identify possible mechanisms of immune escape, and provide valuable information for formulation of future RVA vaccines.

Whole Genome In-Silico Analysis of South African G1P[8] Rotavirus Strains Before and After Vaccine Introduction Over A Period of 14 Years

Rotavirus G1P[8] strains account for more than half of the group A rotavirus (RVA) infections in children under five years of age, globally. A total of 103 stool samples previously characterized as G1P[8] and collected seven years before and seven years after introducing the Rotarix® vaccine in South Africa were processed for whole-genome sequencing. All the strains analyzed had a Wa-like constellation (G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1). South African pre- and post-vaccine G1 strains were clustered in G1 lineage-I and II while the majority (84.2%) of the P[8] strains were grouped in P[8] lineage-III. Several amino acid sites across ten gene segments with the exception of VP7 were under positive selective pressure. Except for the N147D substitution in the antigenic site of eight post-vaccine G1 strains when compared to both Rotarix® and pre-vaccine strains, most of the amino acid substitutions in the antigenic regions of post-vaccine G1P[8] strains were already present during the pre-vaccine period. Therefore, Rotarix® did not appear to have an impact on the amino acid differences in the antigenic regions of South African post-vaccine G1P[8] strains. However, continued whole-genome surveillance of RVA strains to decipher genetic changes in the post-vaccine period remains imperative.

Whole genome and in-silico analyses of G1P[8] rotavirus strains from pre- and post-vaccination periods in Rwanda

Rwanda was the first low-income African country to introduce RotaTeq vaccine into its Expanded Programme on Immunization in May 2012. To gain insights into the overall genetic make-up and evolution of Rwandan G1P[8] strains pre- and post-vaccine introduction, rotavirus positive fecal samples collected between 2011 and 2016 from children under the age of 5 years as part of ongoing surveillance were genotyped with conventional RT-PCR based methods and whole genome sequenced using the Illumina MiSeq platform. From a pool of samples sequenced (n = 158), 36 were identified as G1P[8] strains (10 pre-vaccine and 26 post-vaccine), of which 35 exhibited a typical Wa-like genome constellation. However, one post vaccine strain, RVA/Human-wt/RWA/UFS-NGS:MRC-DPRU442/2012/G1P[8], exhibited a RotaTeq vaccine strain constellation of G1-P[8]-I2-R2-C2-M2-A3-N2-T6-E2-H3, with most of the gene segments having a close relationship with a vaccine derived reassortant strain, previously reported in USA in 2010 and Australia in 2012. The study strains segregated into two lineages, each containing a paraphyletic pre- and post-vaccine introduction sub-lineages. In addition, the study strains demonstrated close relationship amongst each other when compared with globally selected group A rotavirus (RVA) G1P[8] reference strains. For VP7 neutralization epitopes, amino acid substitutions observed at positions T91A/V, S195D and M217T in relation to the RotaTeq vaccine were radical in nature and resulted in a change in polarity from a polar to non-polar molecule, while for the VP4, amino acid differences at position D195G was radical in nature and resulted in a change in polarity from a polar to non-polar molecule. The polarity change at position T91A/V of the neutralizing antigens might play a role in generating vaccine-escape mutants, while substitutions at positions S195D and M217T may be due to natural fluctuation of the RVA. Surveillance of RVA at whole genome level will enhance further assessment of vaccine impact on circulating strains, the frequency of reassortment events under natural conditions and epidemiological fitness generated by such events.

Uncovering the First Atypical DS-1-like G1P[8] Rotavirus Strains That Circulated during Pre-Rotavirus Vaccine Introduction Era in South Africa

Emergence of DS-1-like G1P[8] group A rotavirus (RVA) strains during post-rotavirus vaccination period has recently been reported in several countries. This study demonstrates, for the first time, rare atypical DS-1-like G1P[8] RVA strains that circulated in 2008 during pre-vaccine era in South Africa. Rotavirus positive samples were subjected to whole-genome sequencing. Two G1P[8] strains (RVA/Human-wt/ZAF/UFS-NGS-MRC-DPRU1971/2008/G1P[8] and RVA/Human-wt/ZAF/UFS-NGS-MRC-DPRU1973/2008/G1P[8]) possessed a DS-1-like genome constellation background (I2-R2-C2-M2-A2-N2-T2-E2-H2). The outer VP4 and VP7 capsid genes of the two South African G1P[8] strains had the highest nucleotide (amino acid) nt (aa) identities of 99.6–99.9% (99.1–100%) with the VP4 and the VP7 genes of a locally circulating South African strain, RVA/Human-wt/ZAF/MRC-DPRU1039/2008/G1P[8]. All the internal backbone genes (VP1–VP3, VP6, and NSP1-NSP5) had the highest nt (aa) identities with cognate internal genes of another locally circulating South African strain, RVA/Human-wt/ZAF/MRC-DPRU2344/2008/G2P[6]. The two study strains emerged through reassortment mechanism involving locally circulating South African strains, as they were distinctly unrelated to other reported atypical G1P[8] strains. The identification of these G1P[8] double-gene reassortants during the pre-vaccination period strongly supports natural RVA evolutionary mechanisms of the RVA genome. There is a need to maintain long-term whole-genome surveillance to monitor such atypical strains.

Whole-gene analysis of inter-genogroup reassortant rotaviruses from the Dominican Republic: Emergence of equine-like G3 strains and evidence of their reassortment with locally-circulating strains

Inter-genogroup reassortant group A rotavirus (RVA) strains possessing a G3 VP7 gene of putative equine origin (EQL-G3) have been detected in humans since 2013. Here we report detection of EQL-G3P[8] RVA strains from the Dominican Republic collected in 2014-16. Whole-gene analysis of RVA in stool specimens revealed 16 EQL-G3P[8] strains, 3 of which appear to have acquired an N1 NSP1 gene from locally-circulating G9P[8] strains and a novel G2P[8] reassortant possessing 7 EQL-G3-associated genes and 3 genes from a locally-circulating G2P[4] strain. Phylogenetic/genetic analyses of VP7 gene sequences revealed nine G3 lineages (I-IX) with newly-assigned lineage IX encompassing all reported human EQL-G3 strains along with the ancestral equine strain. VP1 and NSP2 gene phylogenies suggest that EQL-G3P[8] strains were introduced into the Dominican Republic from Thailand. The emergence of EQL-G3P[8] strains in the Dominican Republic and their reassortment with locally-circulating RVA could have implications for current vaccination strategies.

Group A Rotavirus VP1 Polymerase and VP2 Core Shell Proteins: Intergenotypic Sequence Variation and In Vitro Functional Compatibility

Group A rotaviruses (RVAs) are classified according to a nucleotide sequence-based system that assigns a genotype to each of the 11 double-stranded RNA (dsRNA) genome segments. For the segment encoding the VP1 polymerase, 22 genotypes (R1 to R22) are defined with an 83% nucleotide identity cutoff value. For the segment encoding the VP2 core shell protein, which is a functional VP1-binding partner, 20 genotypes (C1 to C20) are defined with an 84% nucleotide identity cutoff value. However, the extent to which the VP1 and VP2 proteins encoded by these genotypes differ in their sequences or interactions has not been described. Here, we sought to (i) delineate the relationships and sites of variation for VP1 and VP2 proteins belonging to the known RVA genotypes and (ii) correlate intergenotypic sequence diversity with functional VP1-VP2 interaction(s) during dsRNA synthesis. Using bioinformatic approaches, we revealed which VP1 and VP2 genotypes encode divergent proteins and identified the positional locations of amino acid changes in the context of known structural domains/subdomains. We then employed an in vitro dsRNA synthesis assay to test whether genotype R1, R2, R4, and R7 VP1 polymerases could be enzymatically activated by genotype C1, C2, C4, C5, and C7 VP2 core shell proteins. Genotype combinations that were incompatible informed the rational design and in vitro testing of chimeric mutant VP1 and VP2 proteins. The results of this study connect VP1 and VP2 nucleotide-level diversity to protein-level diversity for the first time, and they provide new insights into regions/residues critical for VP1-VP2 interaction(s) during viral genome replication.IMPORTANCE Group A rotaviruses (RVAs) are widespread in nature, infecting numerous mammalian and avian hosts and causing severe gastroenteritis in human children. RVAs are classified using a system that assigns a genotype to each viral gene according to its nucleotide sequence. To date, 22 genotypes have been described for the gene encoding the viral polymerase (VP1), and 20 genotypes have been described for the gene encoding the core shell protein (VP2). Here, we analyzed if/how the VP1 and VP2 proteins encoded by the known RVA genotypes differ from each other in their sequences. We also used a biochemical approach to test whether the intergenotypic sequence differences influenced how VP1 and VP2 functionally engage each other to mediate RNA synthesis in a test tube. This work is important because it increases our understanding of RVA protein-level diversity and raises new ideas about the VP1-VP2 binding interface(s) that is important for viral replication.

Identification of group A rotaviruses from Zambian fruit bats provides evidence for long-distance dispersal events in Africa

Group A rotavirus (RVA) is a major cause of diarrhea in children worldwide. Although RVA infects many animals, little is known about RVA in bats. The present study investigated the genetic diversity of RVA in Zambian bats. We identified RVA from two straw-colored fruit bats (Eidolon helvum) and an Egyptian fruit bat (Rousettus aegyptiacus), and analyzed the genome sequences of these strains. Genome segments of the RVA strains from Zambian E. helvum showed 97%–99% nucleotide sequence identity with those of other RVA strains from E. helvum in Cameroon, which is 2800 km from the sampling locations. These findings suggest that migratory straw-colored fruit bat species, distributed across sub-Saharan Africa, have the potential to disseminate RVA across long distances. By contrast, the RVA strain from Zambian R. aegyptiacus carried highly divergent NSP2 and NSP4 genes, leading us to propose novel genotypes N21 and E27, respectively. Notably, this RVA strain also shared the same genotype for VP6 and NSP3 with the RVA strains from Zambian E. helvum, suggesting interspecies transmission and genetic reassortment may have occurred between these two bat species in the past. Our study has important implications for RVA dispersal in bat populations, and expands our knowledge of the ecology, diversity and evolutionary relationships of RVA.

•

Detection of group A rotavirus from Zambian fruit bats.

•

Some viral genes were almost identical to those of rotavirus from Cameroonian bats.

•

The findings provide evidence for long-distance dispersal events of rotavirus.

•

First report of novel N21 and E27 genotypes.