Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mBio
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Ecological and Evolutionary Science

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

Ceren Simsek, Victor Max Corman, Hermann Ulrich Everling, Alexander N. Lukashev, Andrea Rasche, Gael Darren Maganga, Tabea Binger, Daan Jansen, Leen Beller, Ward Deboutte, Florian Gloza-Rausch, Antje Seebens-Hoyer, Stoian Yordanov, Augustina Sylverken, Samuel Oppong, Yaw Adu Sarkodie, Peter Vallo, Eric M. Leroy, Mathieu Bourgarel, Kwe Claude Yinda, Marc Van Ranst, Christian Drosten, Jan Felix Drexler, Jelle Matthijnssens
John T. Patton, Editor
Ceren Simsek
aKU Leuven-University of Leuven, Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Leuven, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ceren Simsek
Victor Max Corman
bCharité-Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humbolt-Universität zu Berlin and Berlin Institute of Health, Institute of Virology, Berlin, Germany
oGerman Centre for Infection Research (DZIF), associated partner Charité Universitätsmedizin Berlin, Berlin, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Victor Max Corman
Hermann Ulrich Everling
cInstitute of Virology, University of Bonn Medical Centre, Bonn, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alexander N. Lukashev
dSechenov University, Moscow, Russia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrea Rasche
bCharité-Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humbolt-Universität zu Berlin and Berlin Institute of Health, Institute of Virology, Berlin, Germany
oGerman Centre for Infection Research (DZIF), associated partner Charité Universitätsmedizin Berlin, Berlin, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gael Darren Maganga
eCentre International de Recherches Médicales de Franceville, Franceville, Gabon
fUniversité des Sciences et Technique de Masuku, Institut National d’Agronomie et de Biotechnologies, Franceville, Gabon
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tabea Binger
iKumasi Centre for Collaborative Research in Tropical Medicine (KCCR), Kumasi, Ghana
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daan Jansen
aKU Leuven-University of Leuven, Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Leuven, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Daan Jansen
Leen Beller
aKU Leuven-University of Leuven, Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Leuven, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ward Deboutte
aKU Leuven-University of Leuven, Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Leuven, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Florian Gloza-Rausch
gNoctalis, Centre for Bat Protection and Information, Bad Segeberg, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Antje Seebens-Hoyer
gNoctalis, Centre for Bat Protection and Information, Bad Segeberg, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stoian Yordanov
hForestry Board Directorate of Strandja Natural Park, Malko Tarnovo, Bulgaria
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Augustina Sylverken
iKumasi Centre for Collaborative Research in Tropical Medicine (KCCR), Kumasi, Ghana
jKwame Nkrumah University of Science and Technology, Kumasi, Ghana
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Samuel Oppong
jKwame Nkrumah University of Science and Technology, Kumasi, Ghana
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yaw Adu Sarkodie
jKwame Nkrumah University of Science and Technology, Kumasi, Ghana
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter Vallo
kInstitute of Vertebrate Biology, Academy of Sciences of the Czech Republic, v.v.i., Brno, Czech Republic
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eric M. Leroy
lInstitut de Recherche pour le Développement, UMR 224 (MIVEGEC), IRD/CNRS/Montpellier University, Montpellier, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mathieu Bourgarel
mCIRAD, UMR ASTRE, Harare, Zimbabwe
nASTRE, Montpellier University, CIRAD, INRA, Montpellier, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kwe Claude Yinda
aKU Leuven-University of Leuven, Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Leuven, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Kwe Claude Yinda
Marc Van Ranst
aKU Leuven-University of Leuven, Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Leuven, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christian Drosten
bCharité-Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humbolt-Universität zu Berlin and Berlin Institute of Health, Institute of Virology, Berlin, Germany
oGerman Centre for Infection Research (DZIF), associated partner Charité Universitätsmedizin Berlin, Berlin, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jan Felix Drexler
bCharité-Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humbolt-Universität zu Berlin and Berlin Institute of Health, Institute of Virology, Berlin, Germany
oGerman Centre for Infection Research (DZIF), associated partner Charité Universitätsmedizin Berlin, Berlin, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jan Felix Drexler
Jelle Matthijnssens
aKU Leuven-University of Leuven, Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Leuven, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jelle Matthijnssens
John T. Patton
Indiana University Bloomington
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mBio.02755-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

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.

IMPORTANCE 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. The current effort to characterize bat rotavirus strains from 3 continents sheds light on the vast genetic diversity of rotaviruses and also hints at a bat origin for several atypical rotaviruses in humans and animals, implying that zoonoses of bat rotaviruses might occur more frequently than currently realized.

INTRODUCTION

Rotaviruses are the leading cause of diarrheal disease in the young of mammals and birds. In humans, rotaviruses are responsible for 122,000 to 216,000 deaths in under-5-year-old infants on a yearly basis, mainly in developing countries (1). The Rotavirus genus belongs to the family Reoviridae and contains 9 species designated A to I (RVA to RVI). The rotavirus genome consists of 11 dsRNA segments encoding 6 structural viral proteins (VP1-6) and 6 nonstructural proteins (NSP1-6) (2).

The RVA outer capsid antigens, VP4 and VP7, are used for a dual classification system defining P-genotype (VP4 is protease sensitive) and G-genotype (VP7 is glycosylated), respectively (2). However, as gene reassortment is a common phenomenon after coinfection for viruses with segmented genomes, a more comprehensive classification approach became necessary to better account for the genome evolution and genetic diversity of RVAs. In 2008, a nucleotide sequence-based, complete genome classification system was developed for RVA, defining genotypes for each of the 11 gene segments. These genotypes allowed extending the dual classification to a full “genotype constellation” classification (3, 4). The gene assignments are reported as Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx, where “x” denotes the particular genotype. The Rotavirus Classification Working Group (RCWG) was formed in order to assign new genotypes to rotavirus genes that could not be assigned to an established genotype (5).

Accumulating whole-genome sequencing data demonstrated that there are typical genotype constellations (GCs) present in most animal species. Two of them, Wa-like and DS-1-like, are responsible for most of the human disease and designated I1-R1-C1-M1-A1-N1-T1-E1-H1 and I2-R2-C2-M2-A2-N2-T2-E2-H2, respectively, for the non-G/P genotypes (3). Furthermore, various animal species are known to have specific GCs such as I2-R2-C2-M2-A3/A11/A13-N2-T6-E2-H3 for cattle and other even-toed ungulates (6), I1/I5-R1-M1-A1/A8-N1-T1/T7-E1-H1 for swine (3, 7), I2/I6-R2-C2-M3-A10-N2-T3-E2/E12-H7 for horses (8), and I3-R3-C3-M3-A3/A9-N2-T3-E3-H3/H6 for cats and dogs (9). Partially shared genotype patterns between established GCs, such as Wa-like human RVA strains and porcine RVAs, as well as DS-1-like human RVA strains and bovine RVAs, suggest a common origin and important zoonotic transfer events in the past (3).

Bats belong to the Chiroptera order, which is the second largest order of mammals (10). They harbor a high diversity of viruses, among which are zoonotic viruses such as lyssavirus, Hendra and Nipah viruses, filovirus, and several coronaviruses (11–17). Given their great population densities, migration ability, and proximity to human habitats, bats are often screened for emerging and reemerging viral pathogens (18, 19). Such screenings have resulted in the sporadic identification of rotavirus strains in bats in the last decade. There are reports of RVH in South Korean bats (20) and Cameroonian bats (21), and a novel rotavirus species (tentatively named RVJ) was identified from Schreiber's bats in Serbia (22); however, RVA is the most commonly detected species and there are currently more than 20 bat RVA strains identified. In 2010, Esona and colleagues reported the first partially sequenced RVA strain (RVA/Bat-wt/KEN/KE4852/07/2007/G25P[6]) in a Kenyan Eidolon helvum (straw-colored fruit bat), and the majority of retrieved gene segments were only distantly related to known mammalian RVA strains (23). During the subsequent decade, sporadic and scattered reports have been published about RVA strains in bats collected from serum, gut, and fecal samples in insectivorous and fruit bats. Several of these reports came from Chinese studies (24–27), but bat RVAs were also detected and (partially) characterized from France (28), Brazil (29), Zambia (30, 31), Cameroon (32), Kenya (33), and Saudi Arabia (34). These studies investigated samples from a variety of bat families, such as Rhinolophidae (24, 30), Hipposideridae (25, 26), Vespertilionidae (26, 28), Molossidae, Phyllostomidae (29), Emballonuridae (26, 33), Pteropodidae (23, 31–34), and Rhinopomatidae (34). From some of these novel bat RVA strains, a few gene segments were sequenced, whereas other strains were sequenced completely, often resulting in one or multiple novel genotypes (23, 26, 29, 31, 32).

Even though RVAs are generally considered to have a relatively restricted host range, a number of unusual strains have been described in the literature, suggestive of interspecies transmissions involving bat RVA strains. One example is the RVA/Horse-wt/ARG/E3198/2008/G3P[3] strain that was isolated from a diarrheic foal in Argentina in 2008 (35). Although its GC was distantly related to feline/canine-like RVA strains at that time, two more recent publications showed a closer relationship with Chinese bat RVA strains in several gene segments (24, 25). A second example was the unusual human G3P[3] RVA strain RVA/Human-wt/JPN/12638/2014/G3P[3], isolated from a 4-year-old child with severe gastroenteric symptoms in Japan. Three of its eleven gene segments were closely related to a South African bat RVA strain, suggesting a reassortment involving a bat RVA strain (36). A third example are two unique G20 human RVA strains, RVA/Human-wt/ECU/Ecu534/2006/G20P[28] (37) and RVA/Human-wt/SUR/2014735512/2013/G20P[28] (38). The recent identification of the G20 genotype in a Brazilian bat RVA strain (RVA/Bat-wt/BRA/3081/2013/G20P[x]) also suggests a potential bat reservoir for these human strains (29).

All in all, slowly emerging data on bat RVA strains start to show that some unusual human and animal RVA strains might actually have been derived from bats. Therefore, the global surveillance of novel and reassorted RVA bat strains has to continue in order to better understand the genetic diversity of bat RVA strains, as well as to maintain both public and animal health. Here, we report identification of 11 bat RVA strains from Bulgaria, Gabon, Ghana, and Costa Rica, suggesting evidence of multiple reassortment and host switching events from bats to bats and to other mammals.

RESULTS AND DISCUSSION

Bats are known hosts of various human pathogens, including viruses such as rabies virus, henipaviruses, Marburg virus, SARS-CoV and MERS-CoV (11–17). In addition, there have been sporadic reports on several other RNA viruses in bats, such as paramyxoviruses, picornaviruses, orthoreoviruses, and astroviruses (39–42). Bat rotaviruses have also been sporadically reported during the last decade and it is rotavirus A (RVA) that has been the most frequently reported rotavirus species. This is not very surprising, given that RVA has been detected in a wide range of mammals and birds (43–45). Furthermore, there are plenty of examples in the literature of this enteric pathogen being capable of interspecies transmission, sometimes in combination with reassortments, between various mammalian species, including humans (46). On some occasions, such animal-derived gene segments (e.g., VP7 genotypes from cattle [G8] and pigs [G9 and presumably G12]) or complete GCs (AU-1 like strains from cats) have become established in the human population. This established circulation either happened in a limited geographical region (AU-1 like or G8) or worldwide, such as the epidemiologically important human-pathogenic G9 and G12 RVAs (47, 48).

In order to further investigate the potential of bat RVA strains to spill over between bats or toward other mammalian species, we investigated RVA strains from over 2,000 bats, spanning five countries on three continents. The bat fecal samples that were collected from Bulgaria, Romania, Germany, Gabon, Ghana, and Costa Rica were screened for RVA using a nested reverse transcriptase PCR (RT-PCR) targeting a short piece of the highly conserved polymerase gene (VP1, Table S1 in the supplemental material). This screening yielded 18 positive samples out of the 2,142 screened samples (0.8%) (Table S2). The RVA detection rate per species ranged from 0 to 1.1%, except for Macronycteris gigas (previously Hipposideros gigas) (14.9%). The reason for this higher detection rate is unknown, but could be due to: (i) better matching oligonucleotides used for the detection, (ii) an ongoing RVA outbreak in the sampled caves, or (iii) higher circulation of enteropathogens in M. gigas. RVA-positive samples were collected from five bat families: Pteropodidae, Rhinolophidae, Hipposideridae, Phyllostomidae, and Vespertilionidae, which originated from all sampling sites except Romania.

TABLE S1

RT-PCR oligonucleotides for the initial rotavirus screening against VP1. Download Table S1, DOCX file, 0.01 MB.
Copyright © 2021 Simsek et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S2

Taxonomical annotation, sampling time and location, and RVA PCR detection information of the bat samples. Download Table S2, DOCX file, 0.02 MB.
Copyright © 2021 Simsek et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Eleven near-complete bat RVA genomes, including six novel genotypes.From 16 of the RVA-positive samples, a sufficient amount of sample was available for complete viral genome sequencing using the NetoVIR protocol (Table S3). A total of 118.9 million paired-end (PE) reads (2 × 150 base pairs) and an average of 7 million PE reads/sample were generated by Illumina sequencing (Table 1). Four samples from Gabon and one sample from Germany did not yield any RVA contigs longer than 500 bp and were therefore not investigated further. From 11 samples, near-complete RVA genomes could be retrieved. These RVA samples belonged to 5 out of the 46 tested bat species (10.8%) and 4 out of the 10 (40%) tested families, as shown in the bat phylogenetic tree (Table S2, Fig. S1). The percentage of reads mapping to RVA in each sample ranged from 0 to 90% (Table 1).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Meta-data and NGS summary of the sequenced RVA-positive samples

TABLE S3

RVA-positive bat samples detected by targeted RT-PCR and having undergone viral metagenomics. Download Table S3, DOCX file, 0.01 MB.
Copyright © 2021 Simsek et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S1

RVA-positive bat families and species. Download FIG S1, DOCX file, 0.1 MB.
Copyright © 2021 Simsek et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

The GCs of the 11 bat RVA strains are shown in Table 2. The genotype assignments, including novel VP6 (I30), NSP1 (A36), NSP2 (N23), and NSP4 (E28) genotypes for some of the Gabonese strains and NSP1 (A32) and NSP3 (T23) genotypes for the strain from Costa Rica were made according to the guidelines determined by the RCWG (4). Although the NSP5 gene segment of RVA/Bat-wt/CRC/KCR10-93/2010/G20P[47] most likely also represents a novel genotype, we were not able to retrieve the complete open reading frame (ORF) (despite several attempts using RT-PCR and Sanger sequencing), which is required for the assignment of a novel genotype (5). Particular GCs were identified in different geographic locations (Table 2). Gabonese strains were similar to each other, with certain genotypes shared with the Bulgarian strains (G3, P[3], C3, M3, N3, T3, and E3). However, they do not cluster phylogenetically closely together (see above), indicating nonrecent reassortment events. KCR10-93 also possessed a unique GC, except for the VP4 genotype P[47], which was shared with the Ghanaian strain. Interestingly, these two VP4 genes were very closely related (see above), suggesting a recent reassortment event. Gabonese GKS-912, GKS-926, and GKS-934 appeared to have a coinfection, as multiple genotypes were identified in these samples for VP2, VP3, VP4, NSP2, NSP3, and NSP4. For GKS-934, two near-complete VP7 gene segments were identified, both belonging to the G3 genotype and yet having substantial nucleotide-level dissimilarity (19%, see below). This was also the case for K212, possessing two distinct M14 genotypes with 12% nucleotide sequence distance.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 2

Color-coded GCs of the bat RVA strains identified in this studya

At least seven distinct bat RVA genotype constellations.Despite that most animal species, including humans, have a limited number of typical RVA GCs, the RVAs harbored by bats show a great genetic diversity. Combining our data with previously published bat RVA genomes showed there are at least 7 distinct bat RVA GCs circulating in the bat population (Table 3), ranging from completely unique to partially overlapping with each other. The Bulgarian RVA/Bat-wt/BGR/BB89-15/2008/G3P[3] and RVA/Bat-wt/BGR/BR89-60/2008/G3P[3] strains were identical or very similar to MSLH14-like RVA strains from China and a partially sequenced strain from Brazil (orange GC in Table 3). Even though at least three of the samples from Gabon possessed more than one RVA strain, they possessed at least three distinct but related GCs (purple GC in Table 3) not previously identified in bats. The RVA/Bat-wt/GHA/K212/2009/G30P[47] strain (green GC in Table 3) was identical or very similar to several previously identified Cameroonian bat RVA strains (32), as well as some partially sequenced bat RVA strains from Zambia (31). RVA/Bat-wt/CRC/KCR10-93/2010/G20P[47] had a distinct GC (brown GC in Table 3), including at least two previously undescribed genotypes, and shared the G20 genotype with RVA/Bat-wt/BRA/3081/2013/G20P[x]. Of interest was the P[47] genotype, which was shared with two African strains from the green GC. The yellow GC in Table 3 was composed of two strains with identical genotypes, RVA/Bat-wt/CMR/BatLy03/2014/G25P[43] and RVA/Bat-wt/SAU/KSA402/2012/G25P[43], detected in Cameroon and Saudi Arabia, respectively, as well as the partially sequenced strain RVA/Bat-wt/KEN/KE4852/2007/G25P[6] from Kenya. Two GCs (blue and dark gray in Table 3) were only represented by a single bat strain from Kenya (RVA/Bat-wt/KEN/BATp39/2015/G36P[51]) and China (RVA/Bat-wt/CHN/GLRL1/2005/G33P[48]), respectively (Table 3).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 3

Color-coded GCs for the bat RVA strains identified in this study, previously published bat RVA strains, as well as a selection of RVA strains from other host species potentially related to batsa

Reassortments among bat RVA strains.Though the GCs are somewhat conserved, there are ample examples for the occurrence of reassortments. In the orange GC, there are some unusual genotypes such as P[10] for VP4, R20 for VP1, and A29 for NSP1 (Table S4a and S4b), which are most likely the results of reassortment events with currently unknown RVA strains (25, 26). Reassortment also takes place between different bat RVA GCs, albeit to a limited extent. For example, RVA/Bat-wt/GAB/GKS-897/2009/G3P[3] is the only strain from the purple GC with the I8 VP6 genotype, which is shared with several strains from the orange GC (RVA/Bat-tc/CHN/MSLH14/2012/G3P[3], RVA/Bat-wt/CHN/BSTM70/2015/G3P[3], RVA/Bat-tc/CHN/MYAS33/2013/G3P[10], and RVA/Bat-wt/CHN/YSSK5/2015/G3P[3]), suggesting a reassortment event. A second example is the shared P[47] VP4 genotype between RVA/Bat-wt/GHA/K212/2009/G30P[47] and RVA/Bat-wt/CMR/BatLy17/2014/G30P[47] (green GC) and RVA/Bat-wt/CRC/KCR10-93/2010/G20P[47] (brown GC) (Table 2). Interestingly, these last three strains were 97 to 100% identical to each other on the nucleotide level for VP4, suggesting a recent reassortment event. Finally, there are also a few bat RVA strains with unusual genotype compositions that do not clearly fall into the seven described GCs. RVA strains RVA/Bat-wt/ZMB/LUS12-14/2012/G3P[3] and RVA/Bat-wt/CHN/YSSK5/2015/G3P[3] possess several genotypes typical for the orange GC, in addition to several other genotypes of unknown origin (Table S4b). Finally, RVA/Bat-wt/KEN/322/Kwale/2015/G3P[10] possesses both genotypes typical to the orange and purple GCs, in addition to some atypical bat RVA genotypes.

TABLE S4

Examples of reassortments and unusual genotype constellations among bat RVA strains and distinct RVA genotype constellations in the same bat species. Download Table S4, DOCX file, 0.02 MB.
Copyright © 2021 Simsek et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

RVA interspecies transmission in bats and potential host range restriction.As demonstrated by the orange GC, RVAs belonging to certain bat families might undergo multiple host switching events. The Bulgarian RVA strains were isolated from rhinolophid bats, whereas the Chinese MSLH14-like strains were found in bats from the Rhinolophidae, Hipposideridae, and Emballonuridae families (Table S4a).

In addition to RVAs potentially being able to infect multiple bat families, individual bat families could also harbor more than one GC, as is shown in Table S4c. Pteropodid bats harbor completely unique GCs (green and yellow), suggesting that the associated RVA strains have a high epidemiologic fitness in these populations. This further indicates that the Pteropodidae, which includes the straw-colored fruit bats, has been a substantial virus reservoir for a long time already, as also shown for Marburg virus, Hendra and Nipah viruses (12–14).

Wide geographic dispersal of bat RVA GCs.The global distribution of the bat RVA GCs revealed several patterns regarding RVA circulation in bats, as shown in Fig. 1. Bat RVAs belonging to the brown, purple, blue, and dark gray GCs have so far only been identified in Costa Rica (and perhaps Brazil), Gabon, Kenya, and China, respectively. On the other hand, the green and yellow GCs were confirmed to be further dispersed, from Cameroon to Saudi Arabia (G25P[43]), and from Ghana and Cameroon to Zambia, respectively, as was previously suggested by Sasaki et al. (31). However, highly similar RVA strains belonging to the orange MSLH14-like GCs span at least three different continents and subcontinents, e.g., Asia, Europe, and possibly Central America. Furthermore, it was also shown that RVA strains with distinct GCs could cocirculate in the same region, as is the case in Cameroon (green, yellow, and purple GCs) and China (orange and dark gray GCs) (Fig. 1).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Geographic distribution of the currently known bat RVA GCs. The colored dots on the map represent the circulating genotypes at the specified locations according to the GCs shown in Table 3.

With powered flight, migratory bats can travel long distances between summer and winter roosts, for foraging, and for searching for a mate (49). Among long-distance migratory bats, E. helvum can cover a range of 270 to 2,500 km (50), and vespertilionid “tree bats” and the subtropical/tropical molossid bats can fly over 1,000 km (51, 52). Global distribution and intercontinental bat virus transfers are also typical of other bat viruses (53). In addition to migration across vast distances, the fact that some distinct GCs seem to have overlapping geographical ranges (such as in China and West Africa in Fig. 1) suggest a fitness advantage for these particular genotypes occurring together. However, there is also ample evidence of gene reassortment events among established GCs (e.g., P[47] in the green and brown GCs, or I8 in the purple and orange GCs), or with RVA strains of currently unknown origin (e.g., A29, A15, or E27).

It is clear that more bats should be sampled in order to have a comprehensive understanding of the driving and restricting forces of bat RVA genetic diversity, or the lack thereof. The detection of P[47] reassortment between Ghanaian and Costa Rican bat RVAs, which are located more than 9,000 km apart, cannot only be explained by the flight ability of bats, but rather the lack of sampling between these two locations. We hypothesize that with the increasing bat RVA sequencing efforts, the geographical and host range of most GCs (such as the blue, dark gray, yellow, and brown) will be significantly expanded.

Potential of interspecies transmissions of bat RVA to mammalian hosts.We further investigated whether there is potential for unusual RVA strains detected in other mammals (including humans) to be a result of an interspecies transmission from bat strains identified in the current and other studies (Table 3).

(a) Likely transmission of bat RVA strains to a horse. In 2013, Miño and colleagues reported an unusual Argentinian equine G3P[3] RVA strain, RVA/Horse-wt/ARG/E3198/2008/G3P[3]. Based on the GC, it was speculated to have a common ancestor with both feline/canine RVA strains as well as the unusual rhesus RVA strain RRV. However, the nucleotide identities were below 90% for most of the genome segments, suggesting that the original host may not be identified yet (35). When more bat RVA genomes became available in subsequent years, Xia and colleagues, and later also Biao He and colleagues, suggested that E3198 might be of bat origin, based on the GCs and nucleotide similarities (25, 26). The close genetic relationship between E3198 and the Bulgarian strains presented here, across all 11 gene segments, might further suggest a bat origin of this unusual equine RVA strain (nucleotide similarities 87 to 97%) (Fig. 2, Fig. 3, Fig. 4, Fig. S2a).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Maximum likelihood trees of the VP1, VP2, VP3, and NSP1 genes of the identified bat RVA strains with known human, bat, and other mammal RVAs. Only bootstrap values above 70 are shown. The genotypes are listed on the right side of the trees. The bat RVA strains identified in this study are shown in bold and colored to their GC, while previously reported bat RVA strains are shown in bold in black. Non-bat RVA strains related to the bat RVA strains identified in this study are marked with silhouettes indicating their host.

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Maximum likelihood trees of the VP4, VP6, and VP7 genes of the identified bat RVA strains with known human, bat, and other mammal RVAs. Only bootstrap values above 70 are shown. The genotypes are listed on the right side of the trees. The bat RVA strains identified in this study are shown in bold and colored to their GC, while previously reported bat RVA strains are shown in bold in black. Non-bat RVA strains related to the bat RVA strains identified in this study are marked with silhouettes indicating their host.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Maximum likelihood trees of the NSP2, NSP3, NSP4, and NSP5 genes of the identified bat RVA strains with known human, bat and other mammal RVAs. Only bootstrap values above 70 are shown. The genotypes are listed on the right side of the trees. The bat RVA strains identified in this study are shown in bold and colored to their GC, while previously reported bat RVA strains are shown in bold in black. Non-bat RVA strains related to the bat RVA strains identified in this study are marked with silhouettes indicating their host.

FIG S2

Heatmap of pairwise nucleotide identities (NI) of the unusual RVA strains RVA/Horse-wt/ARG/E3198/2008/G3P[3] (a), RVA/Simian-tc/ZAF/SA11-H96/1958/G3P[2] (b), RVA/Human-tc/KEN/B10/1987/G3P[2] (c), and RVA/Human-wt/SUR/2014735512/2013/G20P[28] (d). The gray color indicates the nucleotide identities that are below 0.6 or have lack of sequence information for the compared strain. Download FIG S2, DOCX file, 0.3 MB.
Copyright © 2021 Simsek et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

(b) Unexpected high similarities between bat and simian RVA strains. RVA strain RVA/Simian-tc/ZAF/SA11-H96/1958/G3P[2] was isolated from an overtly healthy vervet monkey and has subsequently been used extensively as a laboratory strain in RVA growth, virulence, genome replication, and, in recent years, also reverse genetics research (54–56). However, its origin remained obscure, as related strains were never identified in vervet monkeys or other nonhuman primates. In 2011, Ghosh and colleagues identified an unusual RVA strain, RVA/Human-tc/KEN/B10/1987/G3P[2], from a child in Kenya, which shared 8 of 11 genotypes with SA11-H96. They speculated about a simian or other animal origin of this strain (57). Around the same time, a second human RVA strain, RVA/Human/CHN/ZTR-5/XXXX/G3P[2], nearly identical to SA11-H96 (Figure S2b), was deposited in GenBank as a potential vaccine candidate. However, the controversy about the origin of these SA11-like strains (SA11-H96, B10, and ZTR-5) remained. To our surprise, the purple GC described in this paper, containing only the bat RVA strains from Gabon, showed up to seven genotypes in common with these SA11-like strains (Table 3), with various degrees of nucleotide similarities (Fig. S2b and c). According to phylogenetic analyses, the bat RVAs from Gabon and Kenya clustered with B10 for the VP1, VP6, and NSP4 gene segments, and with all 3 strains (B10, SA11-H96, and ZTR-5) for VP2-4, NSP1, NSP3, and NSP5 (Fig. 2, Fig. 3, Fig. 4).

Not only for SA11-H96, but also for RVA/Simian-tc/USA/RRV/1975/G3P[3] and RVA/Rhesus-tc/USA/TUCH/2002/G3P[24], some close relationships with bat RVA strains were noted. The VP1, VP3, VP4, VP6, VP7, and NSP1-5 genes of RRV clustered closely with one or multiple bat and bat-related RVA strains (Fig. 2, Fig. 3, Fig. 4). For TUCH, the VP1, NSP1, and NSP5 gene segments also clustered close to bat RVA strains (Fig. 2, Fig. 3, Fig. 4).

The finding that the purple SA11-like GC was found in multiple bats in Gabon but only on a single occasion in vervet monkeys and in two unrelated human cases suggests bats as the prime suspect for being the major hosts of these viruses; this makes the monkey and humans strains putative examples of interspecies transmissions. It should, however, be noted that the phylogenetic clustering between these bat, simian, and human strains is still rather variable, and the nucleotide similarities are not as high as between bat RVA strains and RVA/Horse-wt/ARG/E3198/2008/G3P[3] (Figure S2a, Fig. 2, Fig. 3, Fig. 4), suggesting that more RVAs from yet-unsampled animal species will likely cluster in between. However, two other bat strains are of further interest: (i) the bat RVA strain RVA/Bat/KEN/322/Kwale/2015/G3P[10] (only available as a GenBank entry at this point) seems to have a mixed GC possessing characteristics of both the orange and purple GCs (Table 3), where the purple genotypes R8, M5, and A5 of 322/Kwale are of special interest as they are much more closely related to the SA11-like strains than the Gabon bat RVA strains (Fig. 2); and (ii) the bat RVA strain RVA/Bat-wt/KEN/BATp39/2015/G36P[51] (only available in GenBank) possesses a single purple genotype I16, and again this is more closely related to the SA11-like strain B10 compared to the Gabon RVA strains. Taken together, we speculate that with further RVA screenings in bat populations, more bat RVA strains that are closely related to the vervet monkey RVA strain SA11-H96 and human SA11-like RVA strains may be detected.

(c) Evidence of bat RVA strains transmitted to humans? The G3 genotype is usually associated with P[8] genotype in humans RVAs, and combinations such as G3P[3] and G3P[9] are only sporadically found in the human population (58). Nonetheless, in the 2000 to 2001 season, the VP4, VP7, VP6, and NSP4 genes were sequenced from a rare human strain RVA/Human-wt/THA/CMH222/2001/G3P[3], detected in a 2-year-old, severely diarrheic patient in Thailand (41). This strain was reported to have a VP7 gene closely related to RVA/Simian-tc/USA/RRV/1975/G3P[3] and a VP4 gene that was caprine-like. Subsequently, Xia and colleagues speculated that this strain is distinct from typical human RVA GCs and very likely shared a common ancestor with Asian bat RVAs (33). Our study provides further evidence for the bat origin of CMH222, as the VP6 I8 genotype of CMH222 is closely related to RVA/Bat-wt/GAB/GKS-897/2009/G3P[3] (Fig. 3).

Later, Wang and colleagues contributed to the list of unusual Southeast Asian human RVA strains possessing the G3P[9] genotypes; both the RVA/Human-tc/CHN/L621/2006/G3P[9] and RVA/Human-wt/CHN/E2451/2011/G3P[9] strains were isolated from a symptomatic adult and a symptomatic child, respectively (59). Complete genome analyses revealed a high genetic relatedness to strains of feline/canine origin for almost all 11 genes. L621 and E2451 also clustered near the aforementioned unusual RVA/Horse-wt/ARG/E3198/2008/G3P[3] for the VP3, VP6, NSP2, and NSP5 genes, and L621 also clustered with the E3198 NSP3 gene. Here, we observed that these atypical Asian human strains were also closely related to the Bulgarian bat RVA strains for VP3, VP6, NSP2, NSP4, and NSP5 and to Gabonese bat strains for NSP2, NSP3, and NSP4 of the orange GC (Fig. 2, Fig. 3, Fig. 4). These findings further substantiate, as well as complicate, the identification of the likely bat host from which the L621 and E2451 strains likely jumped to humans.

Following these potential zoonosis reports, Esona and colleagues also revealed remarkable findings in Latin America in 2018, where only limited bat RVA information is present to date (38). A human strain, RVA/Human-wt/SUR/2014735512/2013/G20P[28], was isolated in Suriname, and possessed a rare G20 genotype, which was also detected in an Ecuadorian human RVA strain (Ecu534) in 2006. Remarkably, 2014735512 showed high similarities with bat strain RVA/Bat-wt/BRA/3081/2013/G20P[x] for the VP7, NSP3, and NSP5 genes (Figure S2d) and it was speculated to be of bat origin, as these genotypes have not been detected in any other animal species so far. RVA/Bat-wt/CRC/KCR10-93/2010/G20P[47] showed nucleotide similarities for 9 out of 11 gene segments ranging from 82% to 92% with 2014735512, and they also phylogenetically clustered together, albeit not very closely (Fig. 2, Fig. 3, Fig. 4). Even though more evidence is needed, this finding might indicate a bat RVA origin for this rare human RVA strain.

Conclusions. Despite the limited number of bat species that have been screened for rotaviruses, a surprisingly large genetic diversity of RVA strains is presented in this study, including six novel genotypes. With increasing screening efforts, it is doubtless that this diversity will expand both genetically and geographically. We also presented multiple examples of close genetic relatedness of several mammalian and bat rotaviruses. The indicated zoonoses have—to the best of our knowledge—always been restricted to sporadic cases so far and have never resulted in major outbreaks in humans. However, it is believed that the rotavirus genotype constellations currently circulating in humans also have a common ancestor with animal rotaviruses, highlighting that interspecies transmissions allowed establishment in the human population and that this could happen again (3).

Another notable finding is that several gene segments of bat RVA strains and the simian SA11 RVA strain (the latter being used in global rotavirus research for decades) have a common origin. Furthermore, the SA11 strain has been recently used as the backbone of a RVA reverse genetics system, and is therefore likely to be used even more in the future. It would be intriguing to test whether or not SA11 grows well in bat cell lines, or in in vivo infection experiments.

MATERIALS AND METHODS

Sample collection.Fecal samples were collected from 2,142 bats from 10 bat families, representing 46 bat species (Table S2 in the supplemental material). Sample collection took place in Ghana, Gabon, Bulgaria, Romania, Germany, and Costa Rica during 2008 to 2010 as part of investigations of other viruses in bats, such as coronavirus, astrovirus, and picornavirus, as described previously (53, 60–63). Bat species were determined by trained field biologists. For European and Costa Rican studies, bats were caught with mist nets, put into cotton bags, and fecal pellets are collected. Ghanaian fecal droppings were collected with plastic foil from the trees in which E. helvum bats were roosting. The pellets were kept in RNAlater RNA stabilization solution (Qiagen, Hilden, Germany). Gabonese bats were also captured with mist nets just before twilight and were individually euthanized. Bat feces were collected with the corresponding permissions of the host countries in all of the studies.

RT-PCR rotavirus screening and viral metagenomics.Viral RNA was isolated from the fecal specimens as described previously (62). To screen the RVA presence in bats, conserved RVA-specific primer pairs targeting the VP1 gene were used (277-nucleotide-long PCR product) in a hemi-nested and single round reverse transcription PCR (RT-PCR) assay (Table S1). Among the 18 positive specimens (Tables S2 and S3), 16 fecal samples for which sufficient material remained were shipped to the Laboratory of Clinical and Epidemiological Virology, Leuven, Belgium on dry ice, for further complete genome analyses (Table 1).

The NetoVIR protocol was used for viral enrichment of the fecal suspensions, as described previously (64). Briefly, the fecal samples were suspended in Dulbecco’s phosphate-buffered saline (DPBS) and homogenized with a MINILYS homogenizer (Bertin Technologies) for 20 s at 3,000 rpm. The homogenates were centrifuged for 3 min at 17,000 × g and filtered with 0.8-μm PES filters (Sartorius). Filtrates were treated with benzonase (Novagen) and micrococcal nuclease (New England BioLabs) at 37°C for 2 h to remove the free-floating nucleic acids. Subsequently, samples were extracted using the QIAamp Viral RNA minikit (Qiagen) according to the manufacturer's instructions, without addition of carrier RNA to the lysis buffer. Reverse transcription and second strand synthesis were performed by an adjusted version of the whole-transcriptome amplification (WTA2) protocol (Sigma-Aldrich), as described previously (32). A sequencing library was constructed with the Nextera XT library preparation kit (Illumina). The size of the library was checked with Bioanalyzer (Agilent Technologies) with a high sensitivity DNA chip and the 2 nM pooled libraries were sequenced on an Illumina NextSeq 500 platform (2 × 150 bp paired ends).

Data analysis.Low-quality reads, ambiguous bases, and primer and adapter sequences were removed from the paired-end reads with Trimmomatic v0.36 with default parameters (65). Trimmed reads were de novo assembled with metaSPAdes from SPAdes software v3.11.1 using 21, 33, 55, and 77 k-mer lengths (66). The obtained contigs were annotated with DIAMOND v0.9.10 against a nonredundant protein database (67). The contigs annotated as “Rotavirus” were further investigated using the nucleotide BLAST against a nucleotide reference database to identify the gene segments (68). The incomplete contigs were completed in silico by mapping the trimmed reads of corresponding samples against the reference sequence determined by the highest BLASTn nucleotide similarity with the lowest E value using BWA software v0.5.9 (69) and SAMtools v1.6 (70). Open reading frames were determined by the web-based NCBI ORF Finder tool (71) (www.ncbi.nlm.nih.gov/orffinder).

Assignment of GCs and phylogenetic analyses.The genotypes were assigned using the RotaC tool (http://rotac.regatools.be). The sequences whose genotypes could not be determined were sent to the RCWG for assignment of novel genotypes.

Reference strains were downloaded from GenBank in order to represent all the relevant genotypes per gene segment. Codon-based, nucleotide-level multiple sequence alignments were done using MUSCLE (72) with default parameters in MEGA software v7.0.26 (73). Pairwise nucleotide distances were calculated using the number of identical residues in relation to the length of the alignment with the bio3d package in R (74). Alignments were trimmed with trimAL v1.2 with automated1 parameter (75). The optimized number of bootstrap replicates (100 to 1,000) was determined by the autoMRE option, and maximum-likelihood trees were reconstructed with RaxML-NG (76). The GTR+G+I nucleotide substitution model was used for trees of all segments, except for NSP4 and NSP5, as they did not converge after 1,000 bootstraps under the GTR+G+I model (TIM3+I+G and HKY+I+G, respectively). FigTree v1.4.3 from the BEAST package was used for phylogenetic tree visualization and manipulation (77). The GCs were illustrated on a world map using the maps package in R software (78).

Data availability.The data have been deposited with links to BioProject accession number PRJNA562472 in the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject/). The data are also deposited to GenBank under the following accession numbers: MN433617 to MN433627 (BB89-15), MN539284 to MN539294 (BR89-60), MN528116 to MN528126 (GKS-897), MN477236 to MN477246 (GKS-912), MN528101 to MN528115 (GKS-926), MN528075 to MN528085 (GKS-929), MN528086 to MN528100 (GKS-934), MN551587 to MN551597 (GKS-941), MN477225 to MN477235 (GKS-954), MN551598 to MN551608 (KCR10-93), and MN567261 to MN567272 (K212). Reference strains that were used to construct the multiple sequence alignments are listed in Table S5.

TABLE S5

Genbank accession numbers of the reference RVA strains used in the study. Download Table S5, DOCX file, 0.04 MB.
Copyright © 2021 Simsek et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Ethical statement.Bat capture and sampling were conducted with the following permissions: the Wildlife and Hunting Department of the Gabonese Ministry of Water and Forestry (N°003/MEFE-PA/SG/DGEF/DCF and N°0021/MEFE-PA/SG/DGEF/DCF); under clearance 314/5327.74.1.6 from the State Office of Energy and Agriculture, the Environment and Rural Areas Schleswig-Holstein (LANU); and clearances 133/24.03.2008 and 192/26.03.2009 from the Bulgarian Ministry of Environment and Water. For the Ghanaian bats, ethics approval was obtained from the Committee for Human Research, Publications and Ethics of Komfo Anokye Teaching Hospital and School of Medical Sciences, Kwame Nkrumah University of Science and Technology, Kumasi. Research samples were exported under a state agreement between the Republic of Ghana and the Federal Republic of Germany, represented by the City of Hamburg. Additional export permission was obtained from the Veterinary Services of the Ghana Ministry of Food and Agriculture.

ACKNOWLEDGMENTS

This project has received funding from the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie agreement no. 721367 granted to J.F.D., J.M., and M.V.R. EU Horizon 2020 projects EVAg (grant agreement number 653316) and COMPARE (agreement number 643476) granted to C.D., and the Russian Science Foundation grant 19-15-00055 to A.N.L. also provided funding to the current study. H.U.E. had a personal scholarship from the BONFOR intramural program at the University of Bonn. German Federal Ministry of Education and Research (BMBF) (project code 01KIO16D), Deutsche Forschungsgemeinschaft (DFG DR 772-3/1), Deutsche Forschungsgemeinschaft within the Africa Infectious Diseases program gave grants to C.D. and Y.A.S. (DR 772/3-1) and to S.O. (KA1241/18-1) that were also among the funding contributions. A personal scholarship granted to A.R. from the German Academic Exchange Service (DAAD) supported field work in Costa Rica. D.J. was supported by the Fonds Wetenschappelijk Onderzoek (Research foundation Flanders) (1S78019N). L.B. was supported by the Fonds Wetenschappelijk Onderzoek (1S61618N). K.C.Y. was funded by the Interfaculty Council for Development Cooperation (IRO) from the KU Leuven. The computing power in this work was provided by the VSC (Flemish Supercomputer Centre), financed by the FWO and the Flemish government, department EWI.

C.D., J.F.D., J.M., and M.V.R. designed the research; V.M.C., H.U.E., A.N.L., A.R., G.D.M., T.B., F.G.-R., A.S.-H., S.Y., A.S., S.O., Y.A.S., P.V., M.B., and E.M.L. were involved in sample collection; V.M.C., H.U.E., A.N.L., and C.S. performed the research; C.S., D.J., L.B., W.D., H.U.E., V.M.C., and K.C.Y. contributed to data analysis; J.F.D., C.D., C.S., and J.M. drafted the paper; the final version was approved by all coauthors.

FOOTNOTES

    • Received 20 October 2020
    • Accepted 19 November 2020
    • Published 19 January 2021
  • Copyright © 2021 Simsek et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Clark A,
    2. Black R,
    3. Tate J,
    4. Roose A,
    5. Kotloff K,
    6. Lam D,
    7. Blackwelder W,
    8. Parashar U,
    9. Lanata C,
    10. Kang G,
    11. Troeger C,
    12. Platts-Mills J,
    13. Mokdad A,
    14. Sanderson C,
    15. Lamberti L,
    16. Levine M,
    17. Santosham M,
    18. Steele D
    . 2017. Estimating global, regional and national rotavirus deaths in children aged <5 years: current approaches, new analyses and proposed improvements. PLoS One 12:e0183392. doi:10.1371/journal.pone.0183392.
    OpenUrlCrossRef
  2. 2.↵
    1. Estes MK,
    2. Kapikian AZ
    . 2007. Rotaviruses, p 1917–1974. In Knipe DM, Howley PM, Griffin DE, Lamb RA, Martin MA, Roizman B, Straus SE (ed), Fields virology. Kluwer Health/Lippincott, Williams and Wilkins, Philadelphia, PA.
  3. 3.↵
    1. Matthijnssens J,
    2. Ciarlet M,
    3. Heiman E,
    4. Arijs I,
    5. Delbeke T,
    6. McDonald SM,
    7. Palombo EA,
    8. Iturriza-Gómara M,
    9. Maes P,
    10. Patton JT,
    11. Rahman M,
    12. Van Ranst M
    . 2008. Full genome-based classification of rotaviruses reveals a common origin between human Wa-like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains. J Virol 82:3204–3219. doi:10.1128/JVI.02257-07.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Matthijnssens J,
    2. Ciarlet M,
    3. McDonald SM,
    4. Attoui H,
    5. Bányai K,
    6. Brister JR,
    7. Buesa J,
    8. Esona MD,
    9. Estes MK,
    10. Gentsch JR,
    11. Iturriza-Gómara M,
    12. Johne R,
    13. Kirkwood CD,
    14. Martella V,
    15. Mertens PPC,
    16. Nakagomi O,
    17. Parreño V,
    18. Rahman M,
    19. Ruggeri FM,
    20. Saif LJ,
    21. Santos N,
    22. Steyer A,
    23. Taniguchi K,
    24. Patton JT,
    25. Desselberger U,
    26. Van Ranst M
    . 2011. Uniformity of rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group (RCWG). Arch Virol 156:1397–1413. doi:10.1007/s00705-011-1006-z.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Matthijnssens J,
    2. Ciarlet M,
    3. Rahman M,
    4. Attoui H,
    5. Bányai K,
    6. Estes MK,
    7. Gentsch JR,
    8. Iturriza-Gómara M,
    9. Kirkwood CD,
    10. Martella V,
    11. Mertens PPC,
    12. Nakagomi O,
    13. Patton JT,
    14. Ruggeri FM,
    15. Saif LJ,
    16. Santos N,
    17. Steyer A,
    18. Taniguchi K,
    19. Desselberger U,
    20. Van Ranst M
    . 2008. Recommendations for the classification of group A rotaviruses using all 11 genomic RNA segments. Arch Virol 153:1621–1629. doi:10.1007/s00705-008-0155-1.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Matthijnssens J,
    2. Potgieter CA,
    3. Ciarlet M,
    4. Parreño V,
    5. Martella V,
    6. Bányai K,
    7. Garaicoechea L,
    8. Palombo EA,
    9. Novo L,
    10. Zeller M,
    11. Arista S,
    12. Gerna G,
    13. Rahman M,
    14. Van Ranst M
    . 2009. Are human P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla? J Virol 83:2917–2929. doi:10.1128/JVI.02246-08.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Kim H-H,
    2. Matthijnssens J,
    3. Kim H-J,
    4. Kwon H-J,
    5. Park J-G,
    6. Son K-Y,
    7. Ryu E-H,
    8. Kim D-S,
    9. Lee WS,
    10. Kang M-I,
    11. Yang D-K,
    12. Hyun B-H,
    13. Park S-I,
    14. Park S-J,
    15. Cho K-O
    . 2012. Full-length genomic analysis of porcine G9P[23] and G9P[7] rotavirus strains isolated from pigs with diarrhea in South Korea. Infect Genet Evol 12:1427–1435. doi:10.1016/j.meegid.2012.04.028.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Matthijnssens J,
    2. Miño S,
    3. Papp H,
    4. Potgieter C,
    5. Novo L,
    6. Heylen E,
    7. Zeller M,
    8. Garaicoechea L,
    9. Badaracco A,
    10. Lengyel G,
    11. Kisfali P,
    12. Cullinane A,
    13. Collins PJ,
    14. Ciarlet M,
    15. O'Shea H,
    16. Parreño V,
    17. Bányai K,
    18. Barrandeguy M,
    19. Van Ranst M
    . 2012. Complete molecular genome analyses of equine rotavirus A strains from different continents reveal several novel genotypes and a largely conserved genotype constellation. J Gen Virol 93:866–875. doi:10.1099/vir.0.039255-0.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Matthijnssens J,
    2. De Grazia S,
    3. Piessens J,
    4. Heylen E,
    5. Zeller M,
    6. Giammanco GM,
    7. Bányai K,
    8. Buonavoglia C,
    9. Ciarlet M,
    10. Martella V,
    11. Van Ranst M
    . 2011. Multiple reassortment and interspecies transmission events contribute to the diversity of feline, canine and feline/canine-like human group A rotavirus strains. Infect Genet Evol 11:1396–1406. doi:10.1016/j.meegid.2011.05.007.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Schmid R,
    2. Wilson DE,
    3. Reeder DM
    . 1993. Mammal species of the world: a taxonomic and geographic reference. Taxon 42:512. doi:10.2307/1223169.
    OpenUrlCrossRef
  11. 11.↵
    1. Badrane H,
    2. Tordo N
    . 2001. Host switching in lyssavirus history from the Chiroptera to the Carnivora orders. J Virol 75:8096–8104. doi:10.1128/jvi.75.17.8096-8104.2001.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Chua KB,
    2. Lek Koh C,
    3. Hooi PS,
    4. Wee KF,
    5. Khong JH,
    6. Chua BH,
    7. Chan YP,
    8. Lim ME,
    9. Lam SK
    . 2002. Isolation of Nipah virus from Malaysian Island flying-foxes. Microbes Infect 4:145–151. doi:10.1016/S1286-4579(01)01522-2.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    1. Halpin K,
    2. Young PL,
    3. Field HE,
    4. Mackenzie JS
    . 2000. Isolation of Hendra virus from pteropid bats: a natural reservoir of Hendra virus. J Gen Virol 81:1927–1932. doi:10.1099/0022-1317-81-8-1927.
    OpenUrlCrossRefPubMedWeb of Science
  14. 14.↵
    1. Towner JS,
    2. Amman BR,
    3. Sealy TK,
    4. Carroll SAR,
    5. Comer JA,
    6. Kemp A,
    7. Swanepoel R,
    8. Paddock CD,
    9. Balinandi S,
    10. Khristova ML,
    11. Formenty PBH,
    12. Albarino CG,
    13. Miller DM,
    14. Reed ZD,
    15. Kayiwa JT,
    16. Mills JN,
    17. Cannon DL,
    18. Greer PW,
    19. Byaruhanga E,
    20. Farnon EC,
    21. Atimnedi P,
    22. Okware S,
    23. Katongole-Mbidde E,
    24. Downing R,
    25. Tappero JW,
    26. Zaki SR,
    27. Ksiazek TG,
    28. Nichol ST,
    29. Rollin PE
    . 2009. Isolation of genetically diverse Marburg viruses from Egyptian fruit bats. PLoS Pathog 5:e1000536. doi:10.1371/journal.ppat.1000536.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Li W,
    2. Shi Z,
    3. Yu M,
    4. Ren W,
    5. Smith C,
    6. Epstein JH,
    7. Wang H,
    8. Crameri G,
    9. Hu Z,
    10. Zhang H,
    11. Zhang J,
    12. McEachern J,
    13. Field H,
    14. Daszak P,
    15. Eaton BT,
    16. Zhang S,
    17. Wang L-F
    . 2005. Bats are natural reservoirs of SARS-like coronaviruses. Science 310:676–679. doi:10.1126/science.1118391.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Wang Q,
    2. Qi J,
    3. Yuan Y,
    4. Xuan Y,
    5. Han P,
    6. Wan Y,
    7. Ji W,
    8. Li Y,
    9. Wu Y,
    10. Wang J,
    11. Iwamoto A,
    12. Woo PCY,
    13. Yuen K-Y,
    14. Yan J,
    15. Lu G,
    16. Gao GF
    . 2014. Bat origins of MERS-CoV supported by bat coronavirus HKU4 usage of human receptor CD26. Cell Host Microbe 16:328–337. doi:10.1016/j.chom.2014.08.009.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Zhou P,
    2. Yang X-L,
    3. Wang X-G,
    4. Hu B,
    5. Zhang L,
    6. Zhang W,
    7. Si H-R,
    8. Zhu Y,
    9. Li B,
    10. Huang C-L,
    11. Chen H-D,
    12. Chen J,
    13. Luo Y,
    14. Guo H,
    15. Jiang R-D,
    16. Liu M-Q,
    17. Chen Y,
    18. Shen X-R,
    19. Wang X,
    20. Zheng X-S,
    21. Zhao K,
    22. Chen Q-J,
    23. Deng F,
    24. Liu L-L,
    25. Yan B,
    26. Zhan F-X,
    27. Wang Y-Y,
    28. Xiao G-F,
    29. Shi Z-L
    . 2020. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 579:270–273. doi:10.1038/s41586-020-2012-7.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Calisher CH,
    2. Childs JE,
    3. Field HE,
    4. Holmes KV,
    5. Schountz T
    . 2006. Bats: important reservoir hosts of emerging viruses. Clin Microbiol Rev 19:531–545. doi:10.1128/CMR.00017-06.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Olival KJ,
    2. Hosseini PR,
    3. Zambrana-Torrelio C,
    4. Ross N,
    5. Bogich TL,
    6. Daszak P
    . 2017. Host and viral traits predict zoonotic spillover from mammals. Nature 546:646–650. doi:10.1038/nature22975.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Kim HK,
    2. Yoon S-W,
    3. Kim D-J,
    4. Koo B-S,
    5. Noh JY,
    6. Kim JH,
    7. Choi YG,
    8. Na W,
    9. Chang K-T,
    10. Song D,
    11. Jeong DG
    . 2016. Detection of severe acute respiratory syndrome-like, Middle East respiratory syndrome-like bat coronaviruses and group H rotavirus in faeces of Korean bats. Transbound Emerg Dis 63:365–372. doi:10.1111/tbed.12515.
    OpenUrlCrossRef
  21. 21.↵
    1. Yinda CK,
    2. Ghogomu SM,
    3. Conceição-Neto N,
    4. Beller L,
    5. Deboutte W,
    6. Vanhulle E,
    7. Maes P,
    8. Van Ranst M,
    9. Matthijnssens J
    . 2018. Cameroonian fruit bats harbor divergent viruses, including rotavirus H, bastroviruses, and picobirnaviruses using an alternative genetic code. Virus Evol 4:vey008. doi:10.1093/ve/vey008.
    OpenUrlCrossRef
  22. 22.↵
    1. Bányai K,
    2. Kemenesi G,
    3. Budinski I,
    4. Földes F,
    5. Zana B,
    6. Marton S,
    7. Varga-Kugler R,
    8. Oldal M,
    9. Kurucz K,
    10. Jakab F
    . 2017. Candidate new rotavirus species in Schreiber’s bats, Serbia. Infect Genet Evol 48:19–26. doi:10.1016/j.meegid.2016.12.002.
    OpenUrlCrossRef
  23. 23.↵
    1. Esona MD,
    2. Mijatovic-Rustempasic S,
    3. Conrardy C,
    4. Tong S,
    5. Kuzmin IV,
    6. Agwanda B,
    7. Breiman RF,
    8. Banyai K,
    9. Niezgoda M,
    10. Rupprecht CE,
    11. Gentsch JR,
    12. Bowen MD
    . 2010. Reassortant group A rotavirus from straw-colored fruit bat (Eidolon helvum). Emerg Infect Dis 16:1844–1852. doi:10.3201/eid1612.101089.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. He B,
    2. Yang F,
    3. Yang W,
    4. Zhang Y,
    5. Feng Y,
    6. Zhou J,
    7. Xie J,
    8. Feng Y,
    9. Bao X,
    10. Guo H,
    11. Li Y,
    12. Xia L,
    13. Li N,
    14. Matthijnssens J,
    15. Zhang H,
    16. Tu C
    . 2013. Characterization of a novel G3P[3] rotavirus isolated from a lesser horseshoe bat: a distant relative of feline/canine rotaviruses. J Virol 87:12357–12366. doi:10.1128/JVI.02013-13.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Xia L,
    2. Fan Q,
    3. He B,
    4. Xu L,
    5. Zhang F,
    6. Hu T,
    7. Wang Y,
    8. Li N,
    9. Qiu W,
    10. Zheng Y,
    11. Matthijnssens J,
    12. Tu C
    . 2014. The complete genome sequence of a G3P[10] Chinese bat rotavirus suggests multiple bat rotavirus inter-host species transmission events. Infect Genet Evol 28:1–4. doi:10.1016/j.meegid.2014.09.005.
    OpenUrlCrossRef
  26. 26.↵
    1. He B,
    2. Huang X,
    3. Zhang F,
    4. Tan W,
    5. Matthijnssens J,
    6. Qin S,
    7. Xu L,
    8. Zhao Z,
    9. Yang L,
    10. Wang Q,
    11. Hu T,
    12. Bao X,
    13. Wu J,
    14. Tu C
    . 2017. Group A rotaviruses in Chinese bats: genetic composition, serology, and evidence for bat-to-human transmission and reassortment. J Virol 91:e02493-16. doi:10.1128/JVI.02493-16.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Zheng X,
    2. Qiu M,
    3. Guan W,
    4. Li J,
    5. Chen S,
    6. Cheng M,
    7. Huo S,
    8. Chen Z,
    9. Wu Y,
    10. Jiang L,
    11. Chen Q
    . 2018. Viral metagenomics of six bat species in close contact with humans in southern China. Arch Virol 163:73–88. doi:10.1007/s00705-017-3570-3.
    OpenUrlCrossRef
  28. 28.↵
    1. Dacheux L,
    2. Cervantes-Gonzalez M,
    3. Guigon G,
    4. Thiberge J-M,
    5. Vandenbogaert M,
    6. Maufrais C,
    7. Caro V,
    8. Bourhy H
    . 2014. A preliminary study of viral metagenomics of French bat species in contact with humans: identification of new mammalian viruses. PLoS One 9:e87194. doi:10.1371/journal.pone.0087194.
    OpenUrlCrossRef
  29. 29.↵
    1. Asano KM,
    2. Gregori F,
    3. Hora AS,
    4. Scheffer KC,
    5. Fahl WO,
    6. Iamamoto K,
    7. Mori E,
    8. Silva FDF,
    9. Taniwaki SA,
    10. Brandão PE
    . 2016. Group A rotavirus in Brazilian bats: description of novel T15 and H15 genotypes. Arch Virol 161:3225–3230. doi:10.1007/s00705-016-3010-9.
    OpenUrlCrossRef
  30. 30.↵
    1. Sasaki M,
    2. Orba Y,
    3. Sasaki S,
    4. Gonzalez G,
    5. Ishii A,
    6. Hang'ombe BM,
    7. Mweene AS,
    8. Ito K,
    9. Sawa H
    . 2016. Multi-reassortant G3P[3] group A rotavirus in a horseshoe bat in Zambia. J Gen Virol 97:2488–2493. doi:10.1099/jgv.0.000591.
    OpenUrlCrossRef
  31. 31.↵
    1. Sasaki M,
    2. Kajihara M,
    3. Changula K,
    4. Mori-Kajihara A,
    5. Ogawa H,
    6. Hang'ombe BM,
    7. Mweene AS,
    8. Simuunza M,
    9. Yoshida R,
    10. Carr M,
    11. Orba Y,
    12. Takada A,
    13. Sawa H
    . 2018. Identification of group A rotaviruses from Zambian fruit bats provides evidence for long-distance dispersal events in Africa. Infect Genet Evol 63:104–109. doi:10.1016/j.meegid.2018.05.016.
    OpenUrlCrossRef
  32. 32.↵
    1. Yinda CK,
    2. Zeller M,
    3. Conceição-Neto N,
    4. Maes P,
    5. Deboutte W,
    6. Beller L,
    7. Heylen E,
    8. Ghogomu SM,
    9. Van Ranst M,
    10. Matthijnssens J
    . 2016. Novel highly divergent reassortant bat rotaviruses in Cameroon, without evidence of zoonosis. Sci Rep 6:34209. doi:10.1038/srep34209.
    OpenUrlCrossRef
  33. 33.↵
    1. Waruhiu C,
    2. Ommeh S,
    3. Obanda V,
    4. Agwanda B,
    5. Gakuya F,
    6. Ge X-Y,
    7. Yang X-L,
    8. Wu L-J,
    9. Zohaib A,
    10. Hu B,
    11. Shi Z-L
    . 2017. Molecular detection of viruses in Kenyan bats and discovery of novel astroviruses, caliciviruses and rotaviruses. Virol Sin 32:101–114. doi:10.1007/s12250-016-3930-2.
    OpenUrlCrossRef
  34. 34.↵
    1. Mishra N,
    2. Fagbo SF,
    3. Alagaili AN,
    4. Nitido A,
    5. Williams SH,
    6. Ng J,
    7. Lee B,
    8. Durosinlorun A,
    9. Garcia JA,
    10. Jain K,
    11. Kapoor V,
    12. Epstein JH,
    13. Briese T,
    14. Memish ZA,
    15. Olival KJ,
    16. Lipkin WI
    . 2019. A viral metagenomic survey identifies known and novel mammalian viruses in bats from Saudi Arabia. PLoS One 14:e0214227. doi:10.1371/journal.pone.0214227.
    OpenUrlCrossRef
  35. 35.↵
    1. Miño S,
    2. Matthijnssens J,
    3. Badaracco A,
    4. Garaicoechea L,
    5. Zeller M,
    6. Heylen E,
    7. Van Ranst M,
    8. Barrandeguy M,
    9. Parreño V
    . 2013. Equine G3P[3] rotavirus strain E3198 related to simian RRV and feline/canine-like rotaviruses based on complete genome analyses. Vet Microbiol 161:239–246. doi:10.1016/j.vetmic.2012.07.033.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Okitsu S,
    2. Hikita T,
    3. Thongprachum A,
    4. Khamrin P,
    5. Takanashi S,
    6. Hayakawa S,
    7. Maneekarn N,
    8. Ushijima H
    . 2018. Detection and molecular characterization of two rare G8P[14] and G3P[3] rotavirus strains collected from children with acute gastroenteritis in Japan. Infect Genet Evol 62:95–108. doi:10.1016/j.meegid.2018.04.011.
    OpenUrlCrossRef
  37. 37.↵
    1. Solberg OD,
    2. Hasing ME,
    3. Trueba G,
    4. Eisenberg JNS
    . 2009. Characterization of novel VP7, VP4, and VP6 genotypes of a previously untypeable group A rotavirus. Virology 385:58–67. doi:10.1016/j.virol.2008.11.026.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Esona MD,
    2. Roy S,
    3. Rungsrisuriyachai K,
    4. Gautam R,
    5. Hermelijn S,
    6. Rey-Benito G,
    7. Bowen MD
    . 2018. Molecular characterization of a human G20P[28] rotavirus a strain with multiple genes related to bat rotaviruses. Infect Genet Evol 57:166–170. doi:10.1016/j.meegid.2017.11.025.
    OpenUrlCrossRef
  39. 39.↵
    1. Drexler JF,
    2. Corman VM,
    3. Müller MA,
    4. Maganga GD,
    5. Vallo P,
    6. Binger T,
    7. Gloza-Rausch F,
    8. Cottontail VM,
    9. Rasche A,
    10. Yordanov S,
    11. Seebens A,
    12. Knörnschild M,
    13. Oppong S,
    14. Sarkodie YA,
    15. Pongombo C,
    16. Lukashev AN,
    17. Schmidt-Chanasit J,
    18. Stöcker A,
    19. Carneiro AJB,
    20. Erbar S,
    21. Maisner A,
    22. Fronhoffs F,
    23. Buettner R,
    24. Kalko EKV,
    25. Kruppa T,
    26. Franke CR,
    27. Kallies R,
    28. Yandoko ERN,
    29. Herrler G,
    30. Reusken C,
    31. Hassanin A,
    32. Krüger DH,
    33. Matthee S,
    34. Ulrich RG,
    35. Leroy EM,
    36. Drosten C
    . 2012. Bats host major mammalian paramyxoviruses. Nat Commun 3:796. doi:10.1038/ncomms1796.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Lau SKP,
    2. Woo PCY,
    3. Lai KKY,
    4. Huang Y,
    5. Yip CCY,
    6. Shek C-T,
    7. Lee P,
    8. Lam CSF,
    9. Chan K-H,
    10. Yuen K-Y
    . 2011. Complete genome analysis of three novel picornaviruses from diverse bat species. J Virol 85:8819–8828. doi:10.1128/JVI.02364-10.
    OpenUrlAbstract/FREE Full Text
  41. 41.↵
    1. Pritchard LI,
    2. Chua KB,
    3. Cummins D,
    4. Hyatt A,
    5. Crameri G,
    6. Eaton BT,
    7. Wang L-F
    . 2006. Pulau virus; a new member of the Nelson Bay orthoreovirus species isolated from fruit bats in Malaysia. Arch Virol 151:229–239. doi:10.1007/s00705-005-0644-4.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    1. Chu DKW,
    2. Poon LLM,
    3. Guan Y,
    4. Peiris JSM
    . 2008. Novel astroviruses in insectivorous bats. J Virol 82:9107–9114. doi:10.1128/JVI.00857-08.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. Holland RE
    . 1990. Some infectious causes of diarrhea in young farm animals. Clin Microbiol Rev 3:345–375. doi:10.1128/cmr.3.4.345.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Guy JS
    . 1998. Virus infections of the gastrointestinal tract of poultry. Poult Sci 77:1166–1175. doi:10.1093/ps/77.8.1166.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Dhama K,
    2. Chauhan RS,
    3. Mahendran M,
    4. Malik SVS
    . 2009. Rotavirus diarrhea in bovines and other domestic animals. Vet Res Commun 33:1–23. doi:10.1007/s11259-008-9070-x.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Cook N,
    2. Bridger J,
    3. Kendall K,
    4. Gomara MI,
    5. El-Attar L,
    6. Gray J
    . 2004. The zoonotic potential of rotavirus. J Infect 48:289–302. doi:10.1016/j.jinf.2004.01.018.
    OpenUrlCrossRefPubMedWeb of Science
  47. 47.↵
    1. Rahman M,
    2. Matthijnssens J,
    3. Goegebuer T,
    4. De Leener K,
    5. Vanderwegen L,
    6. van der Donck I,
    7. Van Hoovels L,
    8. De Vos S,
    9. Azim T,
    10. Van Ranst M
    . 2005. Predominance of rotavirus G9 genotype in children hospitalized for rotavirus gastroenteritis in Belgium during 1999–2003. J Clin Virol 33:1–6. doi:10.1016/j.jcv.2004.09.020.
    OpenUrlCrossRefPubMedWeb of Science
  48. 48.↵
    1. Rahman M,
    2. Matthijnssens J,
    3. Yang X,
    4. Delbeke T,
    5. Arijs I,
    6. Taniguchi K,
    7. Iturriza-Gómara M,
    8. Iftekharuddin N,
    9. Azim T,
    10. Van Ranst M
    . 2007. Evolutionary history and global spread of the emerging G12 human rotaviruses. J Virol 81:2382–2390. doi:10.1128/JVI.01622-06.
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    1. Voigt CC,
    2. Frick WF,
    3. Holderied MW,
    4. Holland R,
    5. Kerth G,
    6. Mello MAR,
    7. Plowright RK,
    8. Swartz S,
    9. Yovel Y
    . 2017. Principles and patterns of bat movements: from aerodynamics to ecology. Q Rev Biol 92:267–287. doi:10.1086/693847.
    OpenUrlCrossRef
  50. 50.↵
    1. Richter HV,
    2. Cumming GS
    . 2008. First application of satellite telemetry to track African straw-coloured fruit bat migration. J Zoology 275:172–176. doi:10.1111/j.1469-7998.2008.00425.x.
    OpenUrlCrossRefWeb of Science
  51. 51.↵
    1. Cryan PM,
    2. Stricker CA,
    3. Wunder MB
    . 2014. Continental-scale, seasonal movements of a heterothermic migratory tree bat. Ecol Appl 24:602–616. doi:10.1890/13-0752.1.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Fleming TH,
    2. Eby P
    . 2005. Ecology of bat migration, p 159–166. In Kunz TH, Fenton B (ed), Bat ecology. University of Chicago Press, Chicago, IL.
  53. 53.↵
    1. Lukashev AN,
    2. Corman VM,
    3. Schacht D,
    4. Gloza-Rausch F,
    5. Seebens-Hoyer A,
    6. Gmyl AP,
    7. Drosten C,
    8. Drexler JF
    . 2017. Close genetic relatedness of picornaviruses from European and Asian bats. J Gen Virol 98:955–961. doi:10.1099/jgv.0.000760.
    OpenUrlCrossRef
  54. 54.↵
    1. Malherbe H,
    2. Harwin R
    . 1963. The cytopathic effects of vervet monkey viruses. S Afr Med J 37:407–411.
    OpenUrlPubMed
  55. 55.↵
    1. Small C,
    2. Barro M,
    3. Brown TL,
    4. Patton JT
    . 2007. Genome heterogeneity of SA11 rotavirus due to reassortment with “O” agent. Virology 359:415–424. doi:10.1016/j.virol.2006.09.024.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Komoto S,
    2. Sasaki J,
    3. Taniguchi K
    . 2006. Reverse genetics system for introduction of site-specific mutations into the double-stranded RNA genome of infectious rotavirus. Proc Natl Acad Sci U S A 103:4646–4651. doi:10.1073/pnas.0509385103.
    OpenUrlAbstract/FREE Full Text
  57. 57.↵
    1. Ghosh S,
    2. Gatheru Z,
    3. Nyangao J,
    4. Adachi N,
    5. Urushibara N,
    6. Kobayashi N
    . 2011. Full genomic analysis of a simian SA11-like G3P[2] rotavirus strain isolated from an asymptomatic infant: identification of novel VP1, VP6 and NSP4 genotypes. Infect Genet Evol 11:57–63. doi:10.1016/j.meegid.2010.10.010.
    OpenUrlCrossRefPubMed
  58. 58.↵
    1. Dóró R,
    2. László B,
    3. Martella V,
    4. Leshem E,
    5. Gentsch J,
    6. Parashar U,
    7. Bányai K
    . 2014. Review of global rotavirus strain prevalence data from six years post vaccine licensure surveillance: is there evidence of strain selection from vaccine pressure? Infect Genet Evol 28:446–461. doi:10.1016/j.meegid.2014.08.017.
    OpenUrlCrossRef
  59. 59.↵
    1. Wang Y-H,
    2. Pang B-B,
    3. Zhou X,
    4. Ghosh S,
    5. Tang W-F,
    6. Peng J-S,
    7. Hu Q,
    8. Zhou D-J,
    9. Kobayashi N
    . 2013. Complex evolutionary patterns of two rare human G3P[9] rotavirus strains possessing a feline/canine-like H6 genotype on an AU-1-like genotype constellation. Infect Genet Evol 16:103–112. doi:10.1016/j.meegid.2013.01.016.
    OpenUrlCrossRefPubMed
  60. 60.↵
    1. Corman VM,
    2. Rasche A,
    3. Diallo TD,
    4. Cottontail VM,
    5. Stöcker A,
    6. Souza B. F d C D,
    7. Corrêa JI,
    8. Carneiro AJB,
    9. Franke CR,
    10. Nagy M,
    11. Metz M,
    12. Knörnschild M,
    13. Kalko EKV,
    14. Ghanem SJ,
    15. Morales KDS,
    16. Salsamendi E,
    17. Spínola M,
    18. Herrler G,
    19. Voigt CC,
    20. Tschapka M,
    21. Drosten C,
    22. Drexler JF
    . 2013. Highly diversified coronaviruses in neotropical bats. J Gen Virol 94:1984–1994. doi:10.1099/vir.0.054841-0.
    OpenUrlCrossRefPubMed
  61. 61.↵
    1. Pfefferle S,
    2. Oppong S,
    3. Drexler JF,
    4. Gloza-Rausch F,
    5. Ipsen A,
    6. Seebens A,
    7. Müller MA,
    8. Annan A,
    9. Vallo P,
    10. Adu-Sarkodie Y,
    11. Kruppa TF,
    12. Drosten C
    . 2009. Distant relatives of severe acute respiratory syndrome coronavirus and close relatives of human coronavirus 229E in bats. Emerg Infect Dis 15:1377–1384. doi:10.3201/eid1509.090224.
    OpenUrlCrossRefPubMed
  62. 62.↵
    1. Drexler JF,
    2. Gloza-Rausch F,
    3. Glende J,
    4. Corman VM,
    5. Muth D,
    6. Goettsche M,
    7. Seebens A,
    8. Niedrig M,
    9. Pfefferle S,
    10. Yordanov S,
    11. Zhelyazkov L,
    12. Hermanns U,
    13. Vallo P,
    14. Lukashev A,
    15. Müller MA,
    16. Deng H,
    17. Herrler G,
    18. Drosten C
    . 2010. Genomic characterization of severe acute respiratory syndrome-related coronavirus in European bats and classification of coronaviruses based on partial RNA-dependent RNA polymerase gene sequences. J Virol 84:11336–11349. doi:10.1128/JVI.00650-10.
    OpenUrlAbstract/FREE Full Text
  63. 63.↵
    1. Rougeron V,
    2. Suquet E,
    3. Maganga GD,
    4. Jiolle D,
    5. Mombo IM,
    6. Bourgarel M,
    7. Motsch P,
    8. Arnathau C,
    9. Durand P,
    10. Drexler F,
    11. Drosten C,
    12. Renaud F,
    13. Prugnolle F,
    14. Leroy EM
    . 2016. Characterization and phylogenetic analysis of new bat astroviruses detected in Gabon, Central Africa. Acta Virol 60:386–392. doi:10.4149/av_2016_04_386.
    OpenUrlCrossRef
  64. 64.↵
    1. Conceição-Neto N,
    2. Zeller M,
    3. Lefrère H,
    4. De Bruyn P,
    5. Beller L,
    6. Deboutte W,
    7. Yinda CK,
    8. Lavigne R,
    9. Maes P,
    10. Ranst M,
    11. Van Heylen E,
    12. Matthijnssens J
    . 2015. Modular approach to customise sample preparation procedures for viral metagenomics: a reproducible protocol for virome analysis. Sci Rep 5:16532. doi:10.1038/srep16532.
    OpenUrlCrossRefPubMed
  65. 65.↵
    1. Bolger AM,
    2. Lohse M,
    3. Usadel B
    . 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120. doi:10.1093/bioinformatics/btu170.
    OpenUrlCrossRefPubMedWeb of Science
  66. 66.↵
    1. Nurk S,
    2. Meleshko D,
    3. Korobeynikov A,
    4. Pevzner PA
    . 2017. MetaSPAdes: a new versatile metagenomic assembler. Genome Res 27:824–834. doi:10.1101/gr.213959.116.
    OpenUrlAbstract/FREE Full Text
  67. 67.↵
    1. Buchfink B,
    2. Xie C,
    3. Huson DH
    . 2015. Fast and sensitive protein alignment using DIAMOND. Nat Methods 12:59–60. doi:10.1038/nmeth.3176.
    OpenUrlCrossRefPubMed
  68. 68.↵
    1. Altschul SF,
    2. Gish W,
    3. Miller W,
    4. Myers EW,
    5. Lipman DJ
    . 1990. Basic local alignment search tool. J Mol Biol 215:403–410. doi:10.1016/S0022-2836(05)80360-2.
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.↵
    1. Li H,
    2. Durbin R
    . 2010. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26:589–595. doi:10.1093/bioinformatics/btp698.
    OpenUrlCrossRefPubMedWeb of Science
  70. 70.↵
    1. Li H,
    2. Handsaker B,
    3. Wysoker A,
    4. Fennell T,
    5. Ruan J,
    6. Homer N,
    7. Marth G,
    8. Abecasis G,
    9. Durbin R, 1000 Genome Project Data Processing Subgroup
    . 2009. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078–2079. doi:10.1093/bioinformatics/btp352.
    OpenUrlCrossRefPubMedWeb of Science
  71. 71.↵
    1. Wheeler DL,
    2. Church DM,
    3. Federhen S,
    4. Lash AE,
    5. Madden TL,
    6. Pontius JU,
    7. Schuler GD,
    8. Schriml LM,
    9. Sequeira E,
    10. Tatusova TA,
    11. Wagner L
    . 2003. Database resources of the National Center for Biotechnology. Nucleic Acids Res 31:28–33. doi:10.1093/nar/gkg033.
    OpenUrlCrossRefPubMedWeb of Science
  72. 72.↵
    1. Edgar RC
    . 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32:1792–1797. doi:10.1093/nar/gkh340.
    OpenUrlCrossRefPubMedWeb of Science
  73. 73.↵
    1. Kumar S,
    2. Stecher G,
    3. Tamura K
    . 2016. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol 33:1870–1874. doi:10.1093/molbev/msw054.
    OpenUrlCrossRefPubMed
  74. 74.↵
    1. Grant BJ,
    2. Rodrigues APC,
    3. ElSawy KM,
    4. McCammon JA,
    5. Caves LSD
    . 2006. Bio3d: an R package for the comparative analysis of protein structures. Bioinformatics 22:2695–2696. doi:10.1093/bioinformatics/btl461.
    OpenUrlCrossRefPubMedWeb of Science
  75. 75.↵
    1. Capella-Gutierrez S,
    2. Silla-Martinez JM,
    3. Gabaldon T
    . 2009. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25:1972–1973. doi:10.1093/bioinformatics/btp348.
    OpenUrlCrossRefPubMedWeb of Science
  76. 76.↵
    1. Kozlov AM,
    2. Darriba D,
    3. Flouri T,
    4. Morel B,
    5. Stamatakis A
    . 2019. RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 35:4453–4455. doi:10.1093/bioinformatics/btz305.
    OpenUrlCrossRefPubMed
  77. 77.↵
    1. Rambaut A
    . FigTree, a graphical viewer of phylogenetic trees. http://tree.bio.ed.ac.uk/software/figtree.
  78. 78.↵
    1. Becker RA,
    2. Wilks AR,
    3. Deckmyn A,
    4. Brownrigg R,
    5. Minka TP
    . 2018. Maps: draw geographical maps. https://CRAN.R-project.org/package=maps.
PreviousNext
Back to top
Download PDF
Citation Tools
At Least Seven Distinct Rotavirus Genotype Constellations in Bats with Evidence of Reassortment and Zoonotic Transmissions
Ceren Simsek, Victor Max Corman, Hermann Ulrich Everling, Alexander N. Lukashev, Andrea Rasche, Gael Darren Maganga, Tabea Binger, Daan Jansen, Leen Beller, Ward Deboutte, Florian Gloza-Rausch, Antje Seebens-Hoyer, Stoian Yordanov, Augustina Sylverken, Samuel Oppong, Yaw Adu Sarkodie, Peter Vallo, Eric M. Leroy, Mathieu Bourgarel, Kwe Claude Yinda, Marc Van Ranst, Christian Drosten, Jan Felix Drexler, Jelle Matthijnssens
mBio Jan 2021, 12 (1) e02755-20; DOI: 10.1128/mBio.02755-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mBio article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
At Least Seven Distinct Rotavirus Genotype Constellations in Bats with Evidence of Reassortment and Zoonotic Transmissions
(Your Name) has forwarded a page to you from mBio
(Your Name) thought you would be interested in this article in mBio.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
At Least Seven Distinct Rotavirus Genotype Constellations in Bats with Evidence of Reassortment and Zoonotic Transmissions
Ceren Simsek, Victor Max Corman, Hermann Ulrich Everling, Alexander N. Lukashev, Andrea Rasche, Gael Darren Maganga, Tabea Binger, Daan Jansen, Leen Beller, Ward Deboutte, Florian Gloza-Rausch, Antje Seebens-Hoyer, Stoian Yordanov, Augustina Sylverken, Samuel Oppong, Yaw Adu Sarkodie, Peter Vallo, Eric M. Leroy, Mathieu Bourgarel, Kwe Claude Yinda, Marc Van Ranst, Christian Drosten, Jan Felix Drexler, Jelle Matthijnssens
mBio Jan 2021, 12 (1) e02755-20; DOI: 10.1128/mBio.02755-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS AND DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Viral metagenomics
bat rotavirus
rotavirus genetic diversity
SA11
zoonosis

Related Articles

Cited By...

About

  • About mBio
  • Editor in Chief
  • Board of Editors
  • AAM Fellows
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Article Types
  • Ethics
  • Contact Us

Follow #mBio

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2150-7511