Next Article in Journal
Library Preparation Based on Transposase Assisted RNA/DNA Hybrid Co-Tagmentation for Next-Generation Sequencing of Human Noroviruses
Next Article in Special Issue
Genetic and Antigenetic Characterization of the Novel Kotalahti Bat Lyssavirus (KBLV)
Previous Article in Journal
Jumbo Phages: A Comparative Genomic Overview of Core Functions and Adaptions for Biological Conflicts
Previous Article in Special Issue
Evaluation of Two Real-Time, TaqMan Reverse Transcription-PCR Assays for Detection of Rabies Virus in Circulating Variants from Argentina: Influence of Sequence Variation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Novel Rhabdoviruses in Chinese Bats

1
CAS Key Laboratory of Special Pathogens and Biosafety, Wuhan Institute of Virology, Chinese Academy of Sciences, Wuhan 430071, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Institut Pasteur, Lyssavirus Epidemiology and Neuropathology Unit, 75724 Paris, France
4
Hunan Provincial Key Laboratory of Medical Virology, College of Biology, Hunan University, Changsha 410082, China
*
Authors to whom correspondence should be addressed.
Submission received: 31 October 2020 / Revised: 24 December 2020 / Accepted: 30 December 2020 / Published: 5 January 2021
(This article belongs to the Special Issue Lyssaviruses and Other Bat Rhabdoviruses)

Abstract

:
Bats, the second largest order of mammals worldwide, harbor specific characteristics such as sustaining flight, a special immune system, unique habits, and ecological niches. In addition, they are the natural reservoirs of a variety of emerging or re-emerging zoonotic pathogens. Rhabdoviridae is one of the most diverse families of RNA viruses, which consists of 20 ecologically diverse genera, infecting plants, mammals, birds, reptiles, and fish. To date, three bat-related genera are described, named Lyssavirus, Vesiculovirus, and Ledantevirus. However, the prevalence and the distribution of these bat-related rhabdoviruses remain largely unknown, especially in China. To fill this gap, we performed a large molecular retrospective study based on the real-time reverse transcription polymerase chain reaction (RT-qPCR) detection of lyssavirus in bat samples (1044 brain and 3532 saliva samples, from 63 different bat species) originating from 21 provinces of China during 2006–2018. None of them were positive for lyssavirus, but six bat brains (0.6%) of Rhinolophus bat species, originating from Hubei and Hainan provinces, were positive for vesiculoviruses or ledanteviruses. Based on complete genomes, these viruses were phylogenetically classified into three putative new species, tentatively named Yinshui bat virus (YSBV), Taiyi bat virus (TYBV), and Qiongzhong bat virus (QZBV). These results indicate the novel rhabdoviruses circulated in different Chinese bat populations.

1. Introduction

Bats belong to the second largest mammalian order, Chiroptera, which includes over 1350 extant species belonging to 21 bat families, within the two suborders, Yinpterochiroptera and Yangochiroptera [1,2]. These mammals have an extensive geographical distribution and present various unique characteristics such as the capability of sustaining flight and migration, a gregarious behavior, and an important longevity compared with their size and their lifestyle [3]. In addition, multiple lines of evidence indicate that bats are infected with or can host various pathogens, mainly viruses, without displaying clinical symptoms [4,5]. This ability to resist or tolerate viral infection may be directly related to their special immune system [6,7,8]. Over 200 viruses of 27 families were isolated or detected in bats, and numerous viruses, including zoonotic agents, originate from bats [9,10,11]. Indeed, in the past two decades, numerous emerging or re-emerging infectious diseases were demonstrated or suspected to be connected with bats [11], including infections associated with coronavirus with severe acute respiratory syndrome (SARS) [12,13]; Middle-East respiratory syndrome (MERS) [14]; and, more recently, Coronavirus Disease 2019 (COVID-19) [15], as well as acute fatal hemorrhagic diseases caused by filoviruses (Ebola and Marburg viruses) [16,17] and henipaviruses (Hendra and Nipah viruses) [18,19,20], and viral encephalitis with rabies caused by lyssavirus infection [21,22,23,24].
For the latter, 17 currently recognized species and one tentative species have been described so far (https://talk.ictvonline.org/), which are all suspected to induce rabies. Nearly all of them were found in bats, indicating that these animals are the original and natural reservoir of lyssaviruses [25,26,27]. These viruses belong to the genus Lyssavirus in the family Rhabdoviridae, with Rabies lyssavirus (including rabies virus—RABV) acting as the prototype species (https://talk.ictvonline.org/). Apart from lyssaviruses, rhabdoviruses of two other genera have been found in bats, with Vesiculovirus and Ledantevirus. Within the genus Ledantevirus, two bat-associated members have been reported in America [28] and China [29], and several different ledanteviruses species have been previously reported in Africa [30,31] and Europe [32]. As for all rhabdoviruses, these bat-related rhabdoviruses are typically enveloped virions with bullet-shaped or bacilliform morphology, and are reported to be in the range of 100–460 nm in length and 45–100 nm in diameter. They are characterized by a linear single stranded negative-sense RNA genome of approximately 10–16 kb with five canonical genes encoding the nucleoprotein (N), the phosphoprotein (P), the matrix protein (M), the glycoprotein (G), and the RNA-dependent RNA polymerase (L) [33,34]. However, many rhabdovirus genomes also encode multiple additional accessory proteins such as vesiculoviruses and lyssaviruses, which commonly encode small proteins in alternative open reading frames (ORFs) in the P gene [35,36], and abundant accessory ORFs were reported between the P and M gene of hapavirus genomes [33,37].
The family Rhabdoviridae exhibits a large ecological diversity with members infecting plants or animals including mammals, birds, reptiles, insects, or fish [33,38,39], and some genera also including many important animal and human viruses such as Lyssavirus (e.g., rabies virus) [21,22,24] and Vesiculovius (e.g., Chandipura virus) with acute encephalitis in humans [40], Tibrovirus (e.g., Bas-Congo virus) with acute hemorrhagic fever in humans [41], or Vesiculovirus and Ephemerovirus affecting cattle [42,43]. There is increasing evidence that bats are playing an important role in rhabdovirus spillover and highlight the importance of performing active surveillance of this animal reservoir, a source of viral genetic diversity and potential rhabdovirus emergences. Indeed, new members in Lyssavirus, Vesiculovirus, and Lendatevirus genera have recently been found in bats. In addition, the presence of a wide diversity of rhabdoviral sequences has been described in bat samples and in bat ectoparasites, especially in Europe [44,45].
However, studies to investigate the presence of lyssavirus and, to a less extent, bat-associated rhabdoviruses, are limited in Asia, especially in Chinese bats. To date, two lyssaviruses isolated in Jilin and Taiwan have been reported [46,47], and rabies virus seroprevalence studies in bats have been conducted in the South of China [48]. Within the genus Vesiculovirus, only one bat-associated virus has been reported in Yunnan in 2017 [29]. No other data are available regarding the circulation of virus members belonging to the genus Ledantevirus or to other genera. To further investigate the prevalence of lyssavirus, and more generally rhabdovirus, in bats in China, we conducted a large retrospective study based on molecular detection, based on the analysis of 4576 samples (1044 brain and 3532 saliva samples) from 63 different bat species, covering a period of 13 years (2006–2018). Based on this large screening, we were not able to detect the presence of any lyssavirus. Three putative novel bat rhabdoviruses were identified from Rhinolophus sinicus or Rhinolophus affinis in two different locations, with Taiyi bat virus (TYBV) belonging to the genus Ledantevirus, and with Yinshui bat virus (YSBV) and Qiongzhong bat virus (QZBV), both related to the genus Vesiculovirus.

2. Materials and Methods

2.1. Ethics Statement

Sample collection and animal experiments were performed according to standard procedures approved by the Animal Ethics Committee of the Wuhan Institute of Virology (WIVH05210201 and WIVA05201703).

2.2. Sample Collection and Storage

Bat samples (brain and saliva) were part of the biocollection housed in the Wuhan Institute of Virology, Chinese Academy of Sciences. They were originally collected from the field, during the period 2006–2018. Bat brain samples were transported from the field in liquid nitrogen, and saliva swab samples were immersed in 1 mL Virus Transport Medium (VTM, pH 7.4 Hank’s balanced salt solution containing bovine serum albumin 1%, penicillin 100 U/mL, and streptomycin 50 µg/mL), before long-term storage at −80 °C. The lists of samples are found in Tables S1–S3.

2.3. RNA Extraction

For each brain sample, a piece of approximate 0.3 cm3 was excised and placed into a 1.5 mL nuclease free plastic tube (Thermo Fisher Scientific, Shanghai, China) containing four or five glass beads with a diameter of 3.5 mm (Sigma, St. Louis, MO, USA). A volume of 1 mL TRIzol reagent (Invitrogen, Shanghai, China) was added and the sample was disrupted at 30 rpm/min for 2 min with TissueLyser II (Qiagen, Shanghai, China). After addition of 0.2 mL of chloroform, vortexing, and incubation at room temperature for 2–15 min, centrifugation at 12,000× g for 15 min at 4 °C was performed. The upper phase was collected, added with 0.5 mL of isopropanol, and incubated for 5–10 min at room temperature after mixing. RNA was then precipitated after centrifugation at 12,000× g for 8 min at 4 °C, and washed with 1 mL of 75% ethanol solution. After another centrifugation of 5 min at 7500× g at 4 °C, the pellet was air dried for 30 min to 1 h at room temperature, resuspended in 50 μL of nuclease free water (Ambion, Shanghai, China), and incubated for 10 min at 60 °C. For saliva swab samples, pools of five samples (40 μL each) were similarly extracted (without the disrupting step) using 0.8 mL of TRIzol LS (Invitrogen), adding 2 μL of glycogen (5 mg/mL) (Ambion) during the isopropanol precipitation step. All RNA samples were stored in −80 °C for subsequent analysis.

2.4. Identification of Bat Species

Bat species was determined on morphological features, and doubtful or unidentified species were identified after partial sequencing of the cytochrome b (CYTB) gene. Briefly, amplification of an approximate 1180 bp fragment of the CYTB gene was performed by polymerase chain reaction (PCR) using forward and reverse primers L14724ag (5′-ATGATATGAAAAACCATCGTTG-3′) and H15915ag (5′-TTTCCNTTTCTGGTTTACAAGAC-3′), respectively, as previously described [49]. PCR products were Sanger sequenced and contigs from assembled sequences were analyzed by BLASTn using the NCBI database. Identification of host species was done on the most closely related sequences with the highest query coverage and a minimum identity of 95%.

2.5. Combo RT-qPCR for Lyssavirus Screening

Extracted RNA from bat brain and saliva samples was screened for lyssavirus with a previously described dual combined real-time reverse transcription polymerase chain reaction (combo RT-qPCR) method [50]. This technique includes two complementary technologies: a probe-based (TaqMan) RT-qPCR for rabies virus (RABV) detection (pan-RABV RT-qPCR) and an intercalating dye (SYBR Green) based RT-qPCR to detect other lyssavirus species (pan-lyssa RT-qPCR). Both of the two techniques amplify an approximate 120 bp conserved region of the polymerase gene. TaqMan RT-qPCR and SYBR Green RT-qPCR were performed with Superscript III Platinum One-Step qRT-PCR Kit (Invitrogen) and Superscript III Platinum SYBR Green One-Step qRT-PCR Kit (Invitrogen), respectively, according to Dacheux et al. (2016) [50]. For each assay, the RABV CVS strain (RNA extracted from virus suspension) and the European bat lyssavirus 1 (EBLV-1) 8918FRA strain (in vitro synthesized RNA from a plasmid encoding the viral target) were used as positive controls for the TaqMan RT-qPCR and SYBR Green RT-qPCR, respectively. Negative control (nuclease free water) was also included.
In addition, the sensitivity (limit of detection) of the SYBR Green RT-qPCR was evaluated for one representative (isolate D170001) of Yinshui bat virus (YSBV). Briefly, the viral target that includes the region used for detection (partial L gene sequence, positions 1641–2050 bp) was amplified using primers D170001_pan-lyssa_target_F1 (5′-TCTAATCAAGACCCATTTTG-3′) and D170001_pan-lyssa_target_R1 (5′-GCAAATTGACAATGCT CCAT-3′). Amplification was performed using TaKaRa Ex Taq kit (Takara Bio, Saint-Germain-en-Laye, France) with 5 μL 10X Ex Taq Buffer, 4 μL dNTP Mixture (2.5 mM each), 1 μL of each primer (10 μM), 0.25 μL TaKaRa Ex Taq (5 units/μL), 36.75 μL nuclease-free water, and 2 μL of sample. The amplification was performed as follows: 98 °C for 2 min followed by 35 cycles consisting of 98 °C for 10 s, 47 °C for 30 s, and 72 °C for 30 s, ending with 72 °C for 7 min. The 410 bp length amplicon was cloned into pCRII-TOPO Vector (Invitrogen, Illkirch, France). Recombinant plasmids were amplified in E. coli TOP10 and purified with NucleoSpin Plasmid, Mini kit for plasmid (Macherey-Nagel, Hoerdt, France). After confirmation of the presence of the insert by Sanger sequencing, positive clones were quantified by Nanodrop (Thermo Fischer Scientific) to determine the number of target copy per mL. Finally, serial dilutions of the plasmids in nuclease free water were used as standard curves.

2.6. Confirmation of Doubtful Samples by Sanger Sequencing

For each doubtful sample, amplicons were cloned using pGEM®-T Easy Vector Systems (Promega, Madison, WI, USA) according to the manufacturer’s instructions. A sample was considered as doubtful after comparison of the value and the shape of the melting temperature (Tm) curve with those obtained with the controls, after SYBR Green RT-qPCR analysis. Between 5 and 10 clones for each sample were submitted for Sanger sequencing, and contigs were subsequently analyzed by BLASTn and BLASTx using the NCBI database.

2.7. Determination of Rhabdovirus Genome Sequences

The genome sequences of the six new rhabdoviruses (isolates 958, 1127, 1017, D170001, D170022, and D170190) were determined by next generation sequencing (NGS), as previously described [43,51]. Host rRNA depletion was carried out with the kit Terminator 5′-Phosphate-Dependent Exonuclease (Epicentre Biotechnologies, Madison, Wi, USA). Between 2 and 4 μg of extracted RNA (from bat brain) was incubated for 1 h at 30 °C with 2 μL of buffer A and 1 μL of Terminator enzyme, in addition to 0.5 μL of RNAsin Ribonuclease inhibitor (Promega, Charbonnières-les-Bains, France), and adjusted to 20 μL with nuclease free water. Depleted RNA was purified using Agencourt RNAclean XP beads (Beckman Coulter, Villepinte, France) at a ratio of 1.8, as recommended by the manufacturer, and eluted in 20 μL nuclease free water. Eight microliters of purified RNA were then reverse transcribed in complementary DNA (cDNA) using random hexamers (Invitrogen, Illkirch, France) and Superscript III reverse transcriptase (Invitrogen), according to manufacturer’s instructions. Afterward, double-stranded DNA (dsDNA) was synthesized at 16 °C for 2 h in 80 μL final reaction mixture, including 20 μL of cDNA, 8 μL of 10× Second-Strand Reaction Buffer (New England Biolabs, Evry, France), 3 μL of dNTP mix (10 mM) (Invitrogen), 1 μL (10 U) of E. coli DNA ligase (New England Biolabs), 4 μL (40 U) of E. coli DNA Polymerase I (New England Biolabs), 1 μL (5 U) of E. coli RNase H (New England Biolabs), and 43 μL of nuclease free water. Finally, dsDNA was purified using AMPure XP (Beckman Coulter) at a ratio of 1.8, as recommended by the manufacturer, and eluted in 20 μL nuclease free. The dsDNA libraries were constructed using the Nextera XT kit (Illumina, Evry, France) and sequenced using a 2 × 150 nucleotide paired-end strategy on the NextSeq500 platform.
NGS data were analyzed using de novo assembly and mapping, as previously described [43]. Briefly, pre-processing of reads, de novo assembly, and mapping (both using CLC Assembly Cell, Qiagen) were performed with a dedicated workflow implemented on the Pasteur Galaxy platform [52] and slightly modified from Troupin et al. (2016) [53] (workflow available on request). In parallel, another dedicated de novo assembly workflow (named Viroscreen), was performed according to Dacheux et al. (2019) [43], after adaptation from Dacheux et al. (2014) [54]. The viral sequence contigs were assembled and manually edited to produce the final sequences of the viral genomes using Sequencher 5.2.4 (Gene Codes Corporation, Ann Arbor, MI, USA). The quality and accuracy of the final genome sequences were checked after a final mapping step of the original cleaned reads and visualized using Tablet [55].
In the case of low covering regions or the presence of gaps, different primer datasets were designed, based either on the sequences obtained by NGS (specific primers) or on the most closely related sequences available in GenBank (Table S4). These primers were used by conventional PCR/nested PCR, and amplicons were submitted for Sanger sequencing. Contig sequences were then used to fulfil the gaps into the genome sequences.
The 5′ ends of the viral genomes were obtained or confirmed by rapid amplification of cDNA ends (RACE) with the kit SMARTer RACE 5′/3′ (Takara, Beijing, China). The 3′ race ends were connected using a phosphorylated end by T4 RNA ligase (Takara) and amplified by hemi-nested PCR, as previously described [56]. Sequences were assembled and manually edited to produce final sequences of the viral genome with Sequencher 5.2.4.
Complete genome sequences were deposited in GenBank under the accession numbers MN607592–MN607597.

2.8. Analysis of Sequences

Putative open reading frames (ORFs) were identified by Sequencher 5.2.4, and the presence of accessory genes was evaluated in silico after comparison with similar genome regions available in GenBank of other rhabdoviruses belonging to the same genus.
Phylogenetic analysis of the newly rhabdoviruses was conducted using reference genome sequences of other rhabdoviruses available and downloaded from GenBank (Table S5). For each virus, ORFs were concatenated and amino acid sequence alignments were performed using ClustalW 2.0 [57]. Phylogenetic tree was constructed with MEGA7.0 [58] and PhyML3.0 [59], using the maximum likelihood method with 1000 bootstraps. Amino acid sequence identity of the viral polymerase was determined by MEGA 7.0.

2.9. Specific RT-qPCR for Rhabdovirus RNA Quantification in Infected Bats

Rhabdovirus RNA quantification in the different organs of the positive bats was performed by RT-qPCR using HiScript II One Step qRT-PCR SYBR Green Kit (Vazyme, Nanjing, China) following the manufacturer’s recommendations. The amplification was performed as follows: 50 °C for 15 min, 95 °C for 5 min, followed by 40 cycles consisting of 95 °C for 10 s and 60 °C for 30 s. A melting curve analysis was done at 60 °C. Specific primers were designed using PrimerExpress3.0.1 software (Applied Biosystems, Shanghai, China) and amplified an approximate 60 bp region of the polymerase gene. Primers 958-L-660F (5′-CAAGAATGATGATAGAGCGACACA-3′) and 958-L-728R (5′-TTTCAGGTGGGAAATGTTTATTGA-3′) were designed for isolate 958 (putative TYBV species); primers DB170001-L-1180R (5′-GTTGCAGCAGGCGTGGTAT-3′) and DB170001-L-1240F (5′-TTTTAAACACTGCGCGAACAGT-3′) for isolates 1017, D170001, and D170022 (putative YSBV species); and primers 1127-L-1421R (5′-CTACCATTCCTCACCCCCTAGA-3′) and 1127-L-1484F (5′-CCGAGACCTATCCCTCCACTT-3′) for isolate 1127 (putative YSBV species). Absolute RNA quantification was performed using specific standard curves. Standards were produced after amplification by conventional PCR of a fragment of nearly 500 bp in the polymerase gene (containing the region used for quantification) from samples 958, 1127, and D170001, which was subsequently cloned into pGEM-T Easy Vector (Promega). Primers’ sequences used for the production of standards were 958-L-378F (5′-CAAGAATGATGATAGAGCGACACA-3′) and 958-L-1169R (5′-ATCTGGGAAACAAAGGTC-3′) for isolate 958, DB170001-981F (5′-AATGAGCCCATCCAACTT-3′) and DB170001-1467R (5′-TAATCCAATGCCTCCACT-3′) for isolate DB170001, and 1127-L-1043F (5′-TTCAGAGCCGACATACAA-3′) and 1127-L-1537R (5′-TAAGACCACCAACGCAAT-3′) for isolate 1127. Amplification was performed using the Platinum Taq DNA Polymerase kit (Invitrogen, Shanghai, China), with 2.5 μL 10× PCR Buffer, 1 μL MgCl2 (50 mM), 0.5 μL dNTP Mixture (10 mM), 1 μL of each primer (10 μM), 0.1 μL Platinum Taq DNA Polymerase (10 units/μL), 17.9 μL nuclease-free water, and 2 μL of sample. For all sets of primers, amplification was performed as follows: 94 °C for 2 min followed by 40 cycles consisting of 94 °C for 20 s, 48 °C for 30 s, and 72 °C for 1 min, ending with 72 °C for 5 min. After cloning, recombinant plasmids were amplified in E. coli DH10B and purified with EndoFree Mini Plasmid Kit II (Tiangen Biotech, Beijing, China). In vitro RNA synthesis was carried out by MAXscript Kit (Ambion) for each recombinant plasmid, and then freshly transcribed RNA was purified by TRIzol (Invitrogen) and quantified by Nanodrop (Thermo Fischer Scientific, Shanghai, China) to determine the number of target copy per mL. Finally, serial dilutions of RNA in nuclease free water were used as standard curves.

2.10. Virus Isolation

Different cell lines were used for virus isolation: C6-36 (ATCC CRL-1660), Vero (ATCC CCL-81), Vero E6 (ATCC CRL-1586), and Rhinolophus sinicus cell lines (brain, lung, liver, kidney, and heart) [60]. Cells’ monolayers were maintained for 80% confluency in DMEM medium (Gibco, Thermo Fisher, USA) supplemented with 10% fetal calf serum (FCS) (Gibco, Thermo Fisher) in 24-well plates at 37 °C with 5% CO2, and then washed with serum free DMEM (sfDMEM) after removing the medium. Approximately 0.3 cm3 of crushed brain samples (in 200 μL sfDMEM) was centrifuged at 10,000× g for 10 min at 4 °C. The supernatant was diluted at 1:10 in sfDMEM before being added to the cells. After 1 h incubation at 37 °C with 5% CO2, the inoculum was removed and replaced with fresh DMEM medium containing 2% FCS. The cells were incubated at 37 °C with 5% CO2 and checked daily for cytopathic effect. After three blind passages, the cells were harvested and total RNA extraction was performed using TRIzol reagent. Detection of rhabdovirus RNA was realized with the specific RT-qPCR developed for viral RNA quantification.
In parallel, virus isolation was attempted in vivo, using newborn BALB/c mice (3 days old) (Charles river, Beijing, China). Samples were diluted at 1:10 in sfDMEM, and then centrifuged at 10,000× g for 10 min at 4 °C. Approximately 5–10 μL of each sample was inoculated intracerebrally to 4–5 newborn mice. The injected animals were monitored daily and euthanized at day 5 and day 7 after inoculation. Brains were collected and total RNA was extracted by TRIzol reagent. Detection of rhabdovirus RNA was realized with the specific RT-qPCR developed as previously described.

3. Results

3.1. Bat Brain and Saliva Sample Collection

A total of 4576 samples, including bat brain (n = 1044) and saliva swab (n = 3532) samples, were collected between 2006 and 2018 in 21 provinces of China from 63 bat species representing 23 genera belonging to 8 families (Figure 1, Table 1 and Tables S1–S3).
The four most represented Chinese provinces were Guangdong (n = 1285; 28.1%), Yunnan (n = 1210; 26.4%), Guangxi (n = 893; 19.5%), and Hubei (n = 582; 12.7%), all together encompassing 3970 bat samples (72.5%) (Figure 1, Tables S1 and S2). Among the eight different bat families sampled, all excepted one (Pteropidae family, 5 species) were associated with insectivorous bats (58 species), with the highest number of samples associated with the Vespertilionidae family (n = 1388; 30.3%), followed by the Rhinolophidae (n = 1243; 27.2%) and Hipposideridae (n = 946; 20.7%) families (Table 1 and Table S3). At the bat species level, the most abundant samples were collected from Rhinolophus sinicus (Rhinolophidae family), with 547 samples (12%), and from Hipposideros armiger (Hipposideridae), with 399 samples (8.7%) (Tables S1 and S2). Only a limited part of the bats was sampled for both brain and saliva samples. The total number of paired samples was 144 (3.1%) (8 species including 2 frugivorous bats species), which were collected in Hubei (30 samples, 2 species) and Yunnan (114 samples, 6 species) provinces from 2016 to 2018 (Table S6).

3.2. Bat Sample Screening and Description of Positive Specimens

All the saliva swab samples (n = 3532) tested with the combo RT-qPCR were found to be negative (Table S2). Similarly, none of the 1044 bat brain samples were detected as positive with the pan-RABV RT-qPCR, but six samples (0.6%, with samples 958, 1017, 1127, D170001, D170022, and D170190) were classified as doubtful with the pan-lyssa RT-qPCR (Table S1). Each amplicon (nearly 120 bp) was cloned and Sanger sequenced, and BLASTn/BLASTx analysis indicated the presence of rhabdovirus polymerase sequences. The ability of the combo RT-qPCR, and more specifically of the pan-lyssa RT-qPCR (SYBR Green), to detect rhabdoviruses other than lyssaviruses is explained by the targeted viral region, corresponding to a relatively conserved region within the polymerase gene of this virus family (Figure S2) [50]. However, despite this relatively high level of conservation and depending on the sequence of the rhabdovirus considered, the presence of mutations in the primers regions probably impacts the sensitivity of detection, as exemplified with the high limit of detection (1 × 107 copies per assay) for isolate D170001 (Figure S3).
All six were detected in bats belonging to the Rhinolophus genus and originating from two provinces: Hainan and Hubei (Figure 1, Table 1). Sample 1127 was obtained from Rhinolophus affinis in Qiongzhong city of Hainan province in August 2007. The five other samples were collected from bats belonging to the species Rhinolophus sinicus collected in Xianning city in Hubei province. Samples 958, 1017, and D170190 were collected in the Taiyi cave in May 2007 (samples 958 and 1017) and March 2017 (sample D170190), whereas samples D170001 and D170022 were collected in the Yinshui cave, 34 km from the Taiyi cave, in March 2017.
The detection rate of rhabdovirus in these three locations was evaluated for both Rhinolophus species (sinicus and affinis) during the period of the study (2006–2018) (Table 2). This overall detection rate was 11.1% (1/9) in Hainan province for Rhinolophus affinis, and 2.8% (3/107) or 4% (2/50) for Rhinolophus sinicus in Hubei province in Taiyi or Yinshui caves, respectively. Over the period, the overall detection rate in Hubei province was 3.2% (5/157). Interestingly, the positivity rate for Taiyi cave, which has been sampled 10 times over the study period, was 8.3% (2/24) in 2007 and 7.1% in 2017 (1/14) (Table 2). Remarkably, positive samples from Yinshui cave were obtained during the unique collection done in 2017 (overall detection rate of 4.7% (3/64) in Hubei province). However, the number of positive samples or total samples tested remained low.

3.3. Genome Characterization of Chinese Bat Rhabdoviruses

Determination of the nucleotide genome sequences of the six bat rhabdoviruses detected by combo RT-qPCR was initially done by NGS, generating 10 million reads in average per sample. Remaining gaps and low coverage regions were resolved by specific PCR or nested-PCR and Sanger sequencing of the corresponding amplicons (directly or after cloning). Determination of genome extremities was performed using RACE approaches. Based on NGS and/or Sanger sequencing, nearly complete genome sequences were obtained for all six samples, ranging from 10,868 to 10,952 nt in length (Figure 2). The 5′ trailer sequence was obtained for all except for isolates 958, whereas the 3′ leader sequence was not successfully recovered, even after several attempts. BLASTn and BLASTx analyses performed on the genome sequences demonstrated that all these viruses were most closely related to the Vesiculovirus genus, except isolate 958, which was found to be more related to the Ledantevirus genus. This taxonomic classification was further confirmed by phylogenetic analysis (see Section 3.4).
All six genomes exhibited a typical rhabdovirus organization, consisting of the five canonical genes encoding, in the following order, the nucleoprotein (N), the phosphoprotein (P), the matrix protein (M), the glycoprotein (G), and the polymerase (L) (Figure 2).
The genomic characteristics of the five isolates (1127, 1017, D170001, D170022, and D170190) related to vesiculoviruses were highly similar. Indeed, each of the N, M, and L proteins were identical in size, with 423 aa (1269 nt), 209 aa (627 nt) and 2105 aa (6315 nt), respectively. Only isolate 1127 exhibited differences in the length of the P (241 aa, 723 nt) and G (521 aa, 1563 nt) proteins compared with the four other isolates, with 255 aa (765 nt) and 513 aa (1539 nt). The transcription initiation (TI) signal was highly conserved, with the TARCAGR sequence (R = A/G), whereas the consensus sequence for the transcription termination (TTP) was YMTGA7 (Y = C/T and M = A/C) (Figure S1). In addition, various potential ORFs representing putative accessory genes (U) (cut off > 140 nt) were identified (from 4 to 5, depending on the isolates) (Figure 2). The ORFs coding for the putative proteins U2 (98 aa, 294 nt), U3 (52 aa, 156 nt), and U4 (65 aa, 196 nt) were found in the P, M, and G genes of the genome for the isolates 1017, respectively, and similar putative proteins were found in the same regions of the genomes for the isolates D170001, D170022, and D170190. Relatively similar ORFs in the P and G genes were observed for isolate 1127, with U2 (50 aa, 150 nt) and U5 (85 aa, 258 nt), respectively. The sequence of the 5′ trailer region was 56 nt in length and identical for isolates 1127, 1017, D170001, D170022, and D170190.
The isolate 958, related to ledanteviruses, presented the shortest genome (10,933 nt) among this genus (excluding the leader sequence, which was not available). The N, P, M, G, and L proteins were 428 aa (1284 nt), 283 aa (849 nt), 194 aa (582 nt), 536 aa (1608 nt), and 2121 aa (6363 nt) in length, respectively. Interestingly, the M protein of this isolate was the shortest compared with other members of the Ledantevirus genus (208–217 aa, 624–651 nt), and the G protein was the second shortest after the G protein of Mount Elmont bat virus (MEBV) (535 aa, 1605 nt). The consensus sequences of the TI and TTP signals were AACGWGW (W = C/G) and CRTGA7 (R = A/G), respectively (Figure S1). Only one additional ORF (U1) was found in the N gene, which putatively encoded a 48 aa protein (144 nt). The length of the 5′ trailer region was uncompleted (32 nt missing).

3.4. Sequences Analysis of Chinese Bat Rhabdoviruses

Phylogenetic analysis of the six Chinese bat rhabdoviruses was conducted on both the complete amino-acid sequences of the L protein and the nucleotide genome sequences, with representative members of the Rhabdoviridae family.
Maximum likelihood phylogenetic analysis performed on the conserved and complete L proteins of these Chinese bat rhabdoviruses, with sequences of 155 other rhabdovirus members available in GenBank (Table S5), demonstrated that isolates 1127, 1017, D170001, D170022, and D170190 were putative new members of the Vesiculovirus genus, whereas the isolate 958 clustered within the Ledantevirus genus (Figure 3). These results, strongly supported by bootstrap values, confirm the preliminary classification based on BLAST analysis. Interestingly, all the new Chinese bat vesiculoviruses clustered with other vesiculoviruses previously found in bats, with American bat vesiculovirus (ABV) from North America [30] and the recently described Jinghong bat virus (JHBV) originating from the same bat species (Rhinolophus affinis) in Yunnan province, China [29]. Within this cluster, the four isolates 1017, D170001, D170022, and D170190 formed a monophyletic group and were considered as a single putative new rhabdovirus species, tentatively named Yinshui bat virus (YSBV) (Figure 3). The isolate 1127, tentatively identified as Qiongzhong bat virus (QZBV) species, did not fall into this cluster, but was clustered with the Chinese vesiculovirus JHBV. Lastly, the isolate 958, which was tentatively classified as Taiyi bat virus (TYBV) species, was also associated with other bat-related members within the Ledantevirus genus, and more specifically with Mount Elgon bat virus (MEBV) (Figure 3) [61].
These results were confirmed by phylogenetic analysis conducted on the nucleotide genome sequences (Figure 4). Indeed, similar topologies of the tree were found for YSBV (isolates 1017, D170001, D170022, and D170190) and QZBV (isolate 1127), which clustered within the Vesiculovirus genus (Figure 4a), and for TYBV (isolate 958), which was strongly related to Ledantevirus genus and more precisely to MEBV (Figure 4b).
The phylogenetic relation between the six novel bat rhabdoviruses and the other members within the same genus (i.e., Vesiculovirus and Ledantevirus) was supported by high amino acid identities of each canonical proteins (Tables S7 and S8). Within the genus Vesiculovirus, the N and L proteins were the most conserved for QZBV species (isolate 1127), with amino acid identities ranging from 43% to 100% and 52% to 91%, respectively, as well as for YSBV species (isolates 1017, D170001, D170022, and D170190), ranging from 44% to 90% and 52% to 76%, respectively (Table S7). Identically, both of these proteins were the most conserved for the genus Ledantevirus, with amino acid identities ranging from 38% to 87% for TYBV species (isolate 958) and from 48% to 76% for N and L proteins respectively (Table S8). In contrast, the P protein was the most variable for YSBV (amino acid identities from 9% to 52%), QZBV (11% to 87%), and TYBV (14% to 74%) (Tables S7 and S8).
Based on sequence comparison and phylogenetic analysis, the four isolates 1017, D170001, D170022, and D170190 exhibited a high homology in the genome sequences, suggesting that they belong to the same potential new species YSBV (Figure 3 and Figure 4a). Indeed, only 109 nucleotide mutations were found after genome sequence comparison, with 21 mutations associated with amino acid residue modification (Figure 5). For isolates D170001 and D170022, isolates that have been collected from Rhinolophus sinicus in Yinshui cave in 2017, a single amino acid (G gene) residue mutation was observed between them, in addition to three synonymous mutations in other parts of the genome. In parallel, only 7 amino acids mutations (four mutations in the P gene, two mutations in the G gene, and one mutation in the L gene) and 26 synonymous mutations were identified between isolates 1017 and D170190, both collected in Taiyi cave from the same bat species in 2007 and 2017, respectively (Figure 5). A total of 21 amino acids mutations were found between D170001/D170022 and 1017/D170190 in P, M, G, and L genes (Figure 5).

3.5. Organ Tropism of the Novel Chinese Bat Rhabdoviruses

Specific RT-qPCRs were designed to quantify the viral RNA concentration of the six novel rhabdoviruses in different organ samples from the original infected bats. Depending on the specimens, several tissues were tested, including brain, heart, liver, spleen, lung, kidney, and intestine (Table 3). For all rhabdoviruses, the viral concentration was high, between 1.12 × 108 copies/g and 1.62 × 1012 copies/g, regardless of the sample tested. The only exception was the concentration of viral RNA found in intestine for isolate 958 (4.68 × 105 copies/g). These results demonstrated that multiple organs can be infected, with a high viral load. This is especially true for YSBV species with isolates D170001 (seven different organs tested) and D170190 (seven organs tested) and, to a lesser extent, for isolates 958 and 1017 (two organs tested with brain and intestine) (Table 3).

3.6. Bat Sample Screening by Specific RT-qPCRs

The three RT-qPCRs specifically designed to detect and quantify TYBV (isolate 958), QZBV (isolate 1127) and YSBV (isolates 1017, D170001, D170022, and D170190) viruses were used on different bat samples available from Rhinolophus sinicus and Rhinolophus affinis, collected in the same location as the positive bat specimens (Yinshui and Taiyi caves in Hubei province, Qiongzhong in Hainan province). A total of 394 samples were tested, including samples already analyzed by the combo RT-qPCR and additional samples (Table S9). All of the 160 brain samples, 29 saliva and 205 feces were found to be negative.

3.7. Virus Isolation

Attempts to isolate the six different new bat rhabdoviruses were conducted with different organs from the original infected specimens, on cell culture with different cell lines, including C6-36 (ATCC CRL-1660), Vero (ATCC CCL-81), Vero E6 (ATCC CRL-1586), and Rhinolophus sinicus cell lines (brain, lung, liver, kidney, and heart), without any visible cytopathogenic effects nor positive RNA detection with specific RT-qPCR. Similarly, virus isolation after intracerebrally inoculation of crushed tissues in newborn mice did not lead to mortality or any clinical symptoms, or to the detection of viral RNA in inoculated brain after euthanasia.

4. Discussion

It is now clearly demonstrated that bats are major reservoirs of viruses, and could also act as vectors of emerging or re-emerging viral zoonosis [4,5,62,63]. Among them, rabies virus was one of the first identified in these animals [64]. Indeed, rabies virus and other members of the Lyssavirus genus have been described in several part of the world in bats, including the Americas (with the unique presence of Rabies lyssavirus species), Africa (with at least three different species of bat related-lyssaviruses—Duvenhage lyssavirus, Lagos bat lyssavirus, and Shimoni bat lyssavirus—and two potential ones—Ikoma lyssavirus and Mokola lyssavirus), and Europe (five different species described—Bokeloh bat lyssavirus, European bat 1 lyssavirus, European bat 2 lyssavirus, Lleida bat lyssavirus, and West Caucasian bat lyssavirus—and one tentative species with Kotalahti bat lyssavirus) [22]. Only one bat lyssavirus has been described so far in Australia and in the Indian subcontinent—Australia bat lyssavirus and Gannoruwa bat lyssavirus, respectively. Excepted for Middle Asian regions, where three different bat lyssavirus species were discovered (Aravan lyssavirus, Irkut lyssavirus, and Khujand lyssavirus), the other parts of Asia, and more especially China, remain largely unexplored. To date, only two bat lyssavirus species have been identified—Taiwan bat lyssavirus in two Japanese pipistrelles (Pipistrellus abramus) in Taiwan in 2017 [46] and with Irkut lyssavirus in a single greater tube-nosed bat (Murina leucogaster) bat in Tonghua county, Jilin province in 2012 [47]. In addition, a single rabies seroprevalence study has been conducted by indirect fluorescent antibody test in bats sampled in South of China during the period 2005–2006, and identified a low rabies seroconversion: 2.2% (15/685) [48].
In this context, the aim of our study was to fill this gap by improving the knowledge on the prevalence of lyssavirus infection in Chinese bat. For this purpose, we realized one of the largest molecular screenings of retrospective bat samples (n = 4576), with 1044 brains and 3532 saliva collected during the period 2006–2018 from 63 different bat species originating from 21 provinces of China. As at least 135 different bat species were reported in this country [65], our representative dataset covered nearly 50% of the bat species diversity found in 62% of the country (21 out of the 34 existing Chinese provinces). We applied a validated technique for pan-lyssavirus detection (combo RT-qPCR), based on two complementary technologies: a probe-based (TaqMan) RT-qPCR for rabies virus (RABV) detection (pan-RABV RT-qPCR) and an intercalating dye (SYBR Green) based RT-qPCR to detect other lyssavirus species (pan-lyssa RT-qPCR) [50]. None of the samples tested were found to be positive for rabies virus, nor other lyssaviruses. This result was not surprising because most of the samples were collected on apparently healthy and flying bats with nets when they left their roost caves at dusk, or directly in the colonies with landing nets [66]. Indeed, multiple lines of evidence demonstrated that bats are sensitive to lyssavirus infection, which can induce signs of rabies (including inability to fly, abnormal behavior, aggressiveness) and ultimately lead to death [22,25,67], even if we cannot exclude asymptomatic infection. All new bat lyssavirus species have been detected in sick or dead animals [22,46,68,69], and the proportion of positive samples has been demonstrated to be higher when targeting this category of animals within bat colonies [70]. On the other hand, in such a healthy bat population, previous studies demonstrated that the positive rate for lyssavirus detection was expected to be very low, around 0.1% in the brain or saliva sample [71,72,73]. However, although the number of brains tested (nearly 1000 samples) may explain the lack of positivity, we would have expected to identify around three positive saliva samples from our dataset (nearly 3500 samples), based on this positive detection rate of 0.1%. One reason could rely on the dilution factor associated with the storage medium used during saliva sample collection and the extraction step. Indeed, saliva swabs were immersed in 1 mL VTM after collection, and only 40 μL was used for RNA extraction. This dilution factor (1:25) could decrease the overall sensitivity of detection. Usually, similar studies are conducted on saliva samples collected in a smaller volume (nearly 200 μL of VTM, TRIzol, or other) and RNA was extracted from the entire volume [71,72]. Another limitation of this study is due to the undersigned sampling. Although our sampling periods cross more than 12 years, the tissue samples were randomly collected in some bat populations with high density and from bats of accidental death during sampling, explaining the sampling disparities at the geographic temporal and bat species level (Figure 1, Tables S1 and S2).
Despite the absence of lyssavirus detection, we were able to detect the presence of other rhabdoviruses with the pan-lyssa RT-qPCR in six bat brain samples (0.6%), with samples 958, 1017, D170001, D170022, and D170190. Among these six positive brain samples, all collected from bats belonging to Rhinolophus genus, we were able to characterize one putative new species (tentatively named Taiyi bat virus or TYBV) of the Ledantevirus genus with isolate 958, and clustered within two potential other species of the Vesiculovirus genus, with Yinshui bat virus (YSBV) and Qiongzhong bat virus (QZBV) species for isolates 1017, D170001, D170022, and D170190, and for isolate 1127, respectively. The extension of the spectrum of detection of this assay, initially dedicated to lyssavirus, rely on the use of a relatively conserved region on the L gene (positions 7272–7390 according to the reference sequence PV- GenBank accession number M13215) as a target (Figure S3) [50]. Previous studies based on molecular techniques targeting similar conserved regions on this gene reported similar broad detection among the Rhabdovirus family [44]. Despite this level of conservation, the presence of potential mutations in these regions, depending on the rhabdovirus species considered, could negatively impact the sensitivity of this technique, and lead to not detecting animals with low viral loads (Figures S2 and S3). It should, however, be noted that the use of RT-qPCRs specific for the new bat rhabdoviruses confirmed the results obtained on the same samples with the combo RT-qPCR dedicated for lyssavirus detection. Additional screening of samples for Rhinolophus sinicus and Rhinolophus affinis from the same location/cave (i.e., Taiyi and Yinshui caves, Qiongzhong location) with specific RT-qPCR remained negative (Table S9).
According to the species demarcation criteria of vesiculoviruses proposed by ICTV, viruses assigned to different species within the genus Vesiculovirus have minimum amino acid sequence divergence of 20%, 10%, and 15% in L, N, and G proteins, respectively. YSBV (isolates 1017, D170001, D170022, and D170190) meets these requirements, with amino acid sequence divergence of 24%, 10%, and 27% for L, N, and G proteins, respectively, when compared with the most closely related other vesiculoviruses (JHBV) (Table S7). QZBV (isolate 1127) is assigned to the same tentative species with the vesiculovirus JHBV, with amino acid sequence divergence of 9%, 0%, and 15% for L, N, and G proteins, respectively (Table S7).
TYBV (isolate 958) was identified in a Rhinolophus sinicus bat collected in 2007 in Hubei province. It represents the first member of the Ledantevirus genus described in China, and the second in Asia, with Oita virus (OITAV) described in bat (Rhinolophus cornutus) in 1972 in Japan [74]. The seven other bat ledanteviruses were mainly reported in Africa, with the exception of Vaprio virus (VAPV) reported in Italy (Europe), in 2016 in a Pipistrellus kuhlii bat species [32], and Kern Canyon virus (KCV) detected in 1956 in California (USA) in a Myotis yumanensis bat specimen [75] (Table S5). Based on phylogenetic analysis, TYBV was most closely related to the Mount Elgon bat virus (MEBV) in the subgroup B of ledanteviruses, which was isolated in Kenya in 1964 in a Rhinolophus hilderbrandtii eloquens bat species [61] (Figure 4b). This subgroup was initially suggested to be related to viruses found in bats belonging to Rhinolophoidea, whereas subgroup A could encompass viruses detected in bats from the family Vespertilionidae [75]. Based on ICTV criteria for species demarcation, TYBV can also be considered as a new species, exhibiting amino acid divergence of 24% and 36% for L and G proteins, respectively, compared with the most closely related MEBV ledantevirus. It significantly exceeds the proposed minimum divergence of 7% and 15% in L and G protein, respectively (Table S8). Additionally, their accessory ORFs differ from all known species in terms of lengths and locations (Figure 2).
Interestingly, the new tentative species TYBV also clustered with two insect-related viruses. The first one was Wuhan Louse Fly virus 5 (WLFV5), identified in Hippoboscidae sp. fly in China in 2013 [76], whereas Kanyawara virus (KYAV) was described in nycteribiid bat flies in Uganda on collared fruit bats (Myonycteris sp) in 2010 [31], and more recently on Angolan soft-furred fruit bats (Lissonycteris angolensis ruwenzorii) in 2017 [77]. TYBV was also closely related to the recently described Bughendera virus (BUGV), which was identified in one nycteribiid bat fly during the same study on Angolan soft-furred fruit bat in 2017 [77] (Figure 4b). As suggested by these previous studies and by our phylogenetic data, we cannot exclude that TYBV could act as an arbovirus, similar to other members of ledanteviruses belonging to subgroup b (Figure 4b).
QZBV (isolate 1127) and YSBV (isolates 1017, D170001, D170022, and D170190), all identified in Rhinolophus sinicus bat species, clustered with a bat-related subgroup among the Vesiculovirus genus (Figure 4a), which was divided according to the geographical origin of isolates. QZBV and YSBV were associated with the other Chinese vesiculovirus Jinghong bat virus (JHBV), whereas the American bat vesiculovirus (ABV) from North America was genetically distinct. Within the Chinese vesiculovirus cluster, QZBV was the closest to JHBV, isolated from a related bat species (Rhinolophus affinis) in Yunnan province in 2011. Considering that QZBV was detected in Yunnan province, which is more than 1000 km away from Hainnan province, we can assume that this subcluster of vesiculovirus is widely distributed over the country of China, at least in bat Rhinolophus genus.
On the other hand, the four isolates of the tentative species YSBV were all discovered in two caves from the same region in Hubei province (Figure 1), which thus represents a unique epidemiological context, both in time and space. Indeed, isolates 1017 and D170190 were found in the Taiyi cave in 2007 and 2017, respectively, suggesting that YSBV in this cave were present (continuously or not) for at least 10 years. In addition, the high level of genetic homology (only seven amino acid mutation) between the two isolates demonstrated a high stability of YSBV in the Taiyi cave over at least a 10-year period (Figure 5). Similarly, isolates D170001 and D170022 identified in Yinshui cave during the same sampling campaign were closely genetically related, with only four mutations (including one non-synonymous position), suggesting that the same virus was circulating in the Rhinolophus sinicus bat colony at this time (Figure 5 and Table 1). Most interestingly and despite the close distance between the two caves, 21 amino acid mutations were found between YSBV isolates from Taiyi and Yinshui caves, indicating the presence of two different virus populations. Such markers would represent potential indicators to investigate the circulation of both YSBV isolates and associated bats (here, Rhinolophus sinicus) between the two caves, although this hypothesis needs to be confirmed on a larger number of isolates, in addition to the investigation of the circulation of bat individuals between the two caves, notably through ringing and capture/recapture campaigns. Indeed, little is known about the ecological factors of this species of bat, especially regarding their fight distance, interactions between colonies, or potential migration [78].
Based either on complete nucleotide genome sequences (Figure 3) nor on complete amino acid polymerase sequences (Figure 4a), both QZBV and YSBV isolates clustered in a distinct group within the Vesiculovirus genus. This phylogroup also encompasses the closely related Chinese bat vesiculovirus JHBV and the more distant American bat vesiculovirus ABV. Even if the number of representative isolates or species remains limited, we can hypothesize that this phylogroup is strongly associated with bats, with a potential geographical clustering (at least between Asian and North American viruses). It will be of interest to investigate the presence of similar bat vesiculoviruses in other parts of the world, especially Europe and Africa. In addition, and similar to TYBV and other ledanteviruses, we cannot exclude at this time the role of arthropods in the life cycle of these bat viruses, in particular, mainly all other vesiculoviruses are considered as arboviruses. Further investigation in order to develop this hypothesis is necessary, in particular through the collection and screening of parasites infecting bats that are present in Taiyi and Yinshui caves.
Interestingly, we were able to demonstrate that different bat rhabdoviruses could circulate simultaneously in the same cave and same bat species. Indeed, we detected YSBV (isolate 1017) and TYBV (isolate 958) viruses in Taiyi cave in different Rhinolophus sinicus bats during the same sampling period (May 2007). This observation is not limited to only rhabdoviruses, because other virus families were found at this time in the same bat species, with SARS-like coronaviruses (unpublished data). This illustrates that virus populations can be large and diverse in the same bat colony. In our study, we did not find bat specimens coinfected with different rhabdoviruses (or even SARS-like viruses, data not shown). However, infection with different viruses has been already observed in bats [54,66].
All QZBV, YSBV, and TYBV isolates were initially detected in the brain of infected bats. Additional molecular screening and virus quantification in other available tissues revealed that the infection was multi-organs, such as intestine, lung liver, or heart, associated with a high viral load. The American vesiculovirus ABV was also detected in the lung and liver of infected bats [28], whereas Chinese bat vesiculovirus JHBV was found in the intestine [29]. This systemic infection raises questions about the impact of infected bats, because they did not exhibit specific clinical signs during their capture, with was done using nets to target bats in flight. It would be interesting to know if these viruses lead to asymptomatic infection, or if the animals collected did not yet have these symptoms (although high viral loads were found in various organs). Further investigations will be also necessary to determine how these viruses are transmitted between individuals and, more particularly, if virus shedding can be possible through urine or feces, although a preliminary screening on 205 feces was negative (Table S9). It would also be important to retest the presence of these viruses in saliva, using less diluted samples. These elements are important to know in order to assess the zoonotic risk of these new viruses. However, it seems unlikely that these bat rhabdoviruses present a high risk to spillover to other animals and humans, especially when considering the low rate of infection and the low opportunity of contact between bat and human. Identification of non or less invasive samples such as saliva, urine, or feces to further investigate these viruses is also a key element in terms of bat conservation.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/13/1/64/s1, Figure S1: Sequences of the transcription initiation (TI) and the transcription termination (TTP) signals of the 6 Chinese bat rhabdoviruses gene junctions, Figure S2: Alignement of a dataset of 161 nucleotide sequences of rhabdoviruses, corresponding to the regions (in grey) targeted by the primers Taq5long and Taq16revlong of the pan-lyssa RT-qPCR, Figure S3: Validation of the pan-lyssa RT-qPCR used for the screening of bat rhabdovirus, with isolate D170001 (YSBV virus) as an example, Table S1: Detection of rhabdoviruses in bat brain samples, Table S2: Detection of rhabdoviruses in bat saliva samples, Table S3: Repartition of the different bat species tested in this study, according to the type of sample analyzed, Table S4: Primers designed for rhabdovirus genome determination, Table S5: List of genome sequences of rhabdovirus used in this study (based on ICTV classification, 2018), Table S6: Bat species with paired samples (brain and saliva), Table S7: Amino acid identities of each canonical proteins (N, P, M, G and L) of isolates, 1017, D170001, D170022 and D170190 (putative species YSBV) and isolate 1127 (putative species TYBV) within the Vesiculovirus genus, Table S8: Amino acid identities of each canonical proteins (N, P, M, G and L) of isolate 958 (putative species QZBV) within the Ledantevirus genus, Table S9: Bat rhabdovirus screening by specific SYBR Green RT-qPCRs for different bat samples of Rhinolophus sinicus and Rhinolophus affinis originated from Hubei and Hainan provinces (All negative).

Author Contributions

Conceptualization, B.H., H.B., L.D., X.-Y.G., and Z.-L.S.; methodology, D.-S.L. and L.D.; formal analysis, D.-S.L. and L.D.; validation, L.D. and Z.-L.S.; sample collection, D.-S.L. and J.W.; viruses screening, D.-S.L.; NGS experiments: D.-S.L. and Y.Z.; PCR experiments, B.L., D.-S.L., and X.-R.S.; virus isolation, D.-S.L. and R.-D.J.; data curation, D.-S.L. and L.D.; NGS data analysis, D.-S.L., L.D., and Y.F.; writing—original draft preparation, D.-S.L. and L.D.; writing—review and editing, D.-S.L., H.B., and L.D.; supervision, L.D. and Z.-L.S.; project administration, B.H., H.B., L.D., and Z.-L.S.; funding acquisition, H.B., L.D., X.-Y.G., and Z.-L.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was jointly funded by Campus France and China Scholarship Council (D.-S.L.) through the PHC Cai Yuanpei 2016 program under grant number 36724VF. This study also received funding from the National Natural Science Foundation of China (81290341), the Scientific and technological basis special project (2013FY113500), and Institut Pasteur, Paris.

Institutional Review Board Statement

Sample collection and animal experiments were performed according to standard procedures approved by the Animal Ethics Committee of the Wuhan Institute of Virology (WIVH05210201 and WIVA05201703).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the present article and in supplementary material. Complete genome sequences were deposited in GenBank under the accession numbers MN607592–MN607597.

Acknowledgments

We are deeply grateful to all the collaborators who participated in the sampling activities: Yun-Zhi Zhang from Dali University; Li-Biao Zhang from Guangdong Institute of Applied Biological Resources; Guang-Jian Zhu from EcoHealth Alliance; Bing-Jie Hu from Hong Kong Queen Mary Hospital; Li-Jun Wu from Wuhan Academy of Agricultural Sciences; Yu-Ji Zhang from Nanjing Liding Medical technology company; Jun-Fa Yuan from Huazhong Agricultural University; Yan Li from Zhejiang University; and Jia-Lu Li, Wen Li, Chu-Ming Luo, Bing Tan, Bo Wang, Ning Wang, Yong-Quan Wu, Xing-Lou Yang, Lei-Ping Zeng, and Qian Zhang from Wuhan Institute of Virology. We thank Shuang Tang and Fei Deng from Microorganisms & Viruses Culture Collection Center, Wuhan Institute of Virology, Chinese Academy of Sciences for providing us the positive control used for lyssavirus screening (RABV, CVS strain). We thank the sequencing facilities of the “Plate-forme de microbiologie mutalisée (P2M)” of Institut Pasteur, Paris, France for technical assistance with NGS sequencing. We would also thank Corinne Maufrais for her continuous assistance in NGS analysis.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Teeling, E.C. A Molecular Phylogeny for Bats Illuminates Biogeography and the Fossil Record. Science 2005, 307, 580–584. [Google Scholar] [CrossRef] [PubMed]
  2. Jones, G.; Teeling, E.C. The Evolution of Echolocation in Bats. Trends Ecol. Evol. 2006, 21, 149–156. [Google Scholar] [CrossRef] [PubMed]
  3. Simmons, N.B. Order Chiroptera. In Mammal Species of the World: A Taxonomic and Geographic Reference, 3rd ed.; Wilson, D.E., Reeder, D.M., Eds.; Johns Hopkins University Press: Baltimore, MD, USA, 2005; pp. 312–529. [Google Scholar]
  4. Wang, L.-F.; Anderson, D.E. Viruses in Bats and Potential Spillover to Animals and Humans. Curr. Opin. Virol. 2019, 34, 79–89. [Google Scholar] [CrossRef]
  5. Hayman, D.T.S. Bats as Viral Reservoirs. Annu. Rev. Virol. 2016, 3, 77–99. [Google Scholar] [CrossRef]
  6. Xie, J.Z.; Li, Y.; Shen, X.R.; Goh, G.; Zhu, Y.; Cui, J.; Wang, L.F.; Shi, Z.L.; Zhou, P. Dampened STING-Dependent Interferon Activation in Bats. Cell Host Microbe 2018, 23, 297–301.e4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Zhou, P.; Tachedjian, M.; Wynne, J.W.; Boyd, V.; Cui, J.; Smith, I.; Cowled, C.; Ng, J.H.J.; Mok, L.; Michalski, W.P.; et al. Contraction of the Type I IFN Locus and Unusual Constitutive Expression of IFN-Alpha in Bats. Proc. Natl. Acad. Sci. USA 2016, 113, 2696–2701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Zhang, G.J.; Cowled, C.; Shi, Z.L.; Huang, Z.Y.; Bishop-Lilly, K.A.; Fang, X.D.; Wynne, J.W.; Xiong, Z.Q.; Baker, M.L.; Zhao, W.; et al. Comparative Analysis of Bat Genomes Provides Insight into the Evolution of Flight and Immunity. Science 2013, 339, 456–460. [Google Scholar] [CrossRef] [Green Version]
  9. Luis, A.D.; Hayman, D.T.S.; O’Shea, T.J.; Cryan, P.M.; Gilbert, A.T.; Pulliam, J.R.C.; Mills, J.N.; Timonin, M.E.; Willis, C.K.R.; Cunningham, A.A.; et al. A Comparison of Bats and Rodents as Reservoirs of Zoonotic Viruses: Are Bats Special? Proc. R. Soc. B Biol. Sci. 2013, 280, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Brook, C.E.; Dobson, A.P. Bats as “special” Reservoirs for Emerging Zoonotic Pathogens. Trends Microbiol. 2015, 23, 172–180. [Google Scholar] [CrossRef]
  11. Calisher, C.H.; Childs, J.E.; Field, H.E.; Holmes, K.V.; Schountz, T. Bats: Important Reservoir Hosts of Emerging Viruses. Clin. Microbiol. Rev. 2006, 19, 531–545. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Ge, X.Y.; Li, J.L.; Yang, X.L.; Chmura, A.A.; Zhu, G.J.; Epstein, J.H.; Mazet, J.K.; Hu, B.; Zhang, W.; Peng, C.; et al. Isolation and Characterization of a Bat SARS-like Coronavirus That Uses the ACE2 Receptor. Nature 2013, 503, 535–538. [Google Scholar] [CrossRef]
  13. Li, W.; Shi, Z.-L.; Yu, M.; Ren, W.; Smith, C.; Epstein, J.H.; Wang, H.; Crameri, G.; Hu, Z.; Zhang, H.; et al. Bats Are Natural Reservoirs of SARS-Like Coronaviruses. Science 2005, 310, 676–679. [Google Scholar] [CrossRef]
  14. Lau, S.K.P.; Li, K.S.M.; Tsang, A.K.L.; Lam, C.S.F.; Ahmed, S.; Chen, H.L.; Chan, K.H.; Woo, P.C.Y.; Yuen, K.Y. Genetic Characterization of Betacoronavirus Lineage C Viruses in Bats Reveals Marked Sequence Divergence in the Spike Protein of Pipistrellus Bat Coronavirus HKU5 in Japanese Pipistrelle: Implications for the Origin of the Novel Middle East Respiratory Syndrome Coronavirus. J. Virol. 2013, 87, 8638–8650. [Google Scholar] [CrossRef] [Green Version]
  15. Zhou, P.; Yang, X.-L.; Wang, X.-G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.-R.; Zhu, Y.; Li, B.; Huang, C.-L.; et al. A Pneumonia Outbreak Associated with a New Coronavirus of Probable Bat Origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [Green Version]
  16. Towner, J.S.; Amman, B.R.; Sealy, T.K.; Carroll, S.A.R.; Comer, J.A.; Kemp, A.; Swanepoel, R.; Paddock, C.D.; Balinandi, S.; Khristova, M.L.; et al. Isolation of Genetically Diverse Marburg Viruses from Egyptian Fruit Bats. PLoS Pathog. 2009, 5, e1000536. [Google Scholar] [CrossRef]
  17. Leroy, E.M.; Kumulungui, B.; Pourrut, X.; Rouquet, P.; Hassanin, A.; Yaba, P.; Delicat, A.; Paweska, J.T.; Gonzalez, J.P.; Swanepoel, R. Fruit Bats as Reservoirs of Ebola Virus. Nature 2005, 438, 575–576. [Google Scholar] [CrossRef]
  18. Mbu’u, C.M.; Mbacham, W.F.; Gontao, P.; Sado Kamdem, S.L.; Nlôga, A.M.N.; Groschup, M.H.; Wade, A.; Fischer, K.; Balkema-Buschmann, A. Henipaviruses at the Interface Between Bats, Livestock and Human Population in Africa. Vector Borne Zoonotic Dis. 2019, 19, 455–465. [Google Scholar] [CrossRef]
  19. Kessler, M.K.; Becker, D.J.; Peel, A.J.; Justice, N.V.; Lunn, T.; Crowley, D.E.; Jones, D.N.; Eby, P.; Sánchez, C.A.; Plowright, R.K. Changing Resource Landscapes and Spillover of Henipaviruses. Ann. N. Y. Acad. Sci. 2018, 1429, 78–99. [Google Scholar] [CrossRef]
  20. Ang, B.S.P.; Lim, T.C.C.; Wang, L. Nipah Virus Infection. J. Clin. Microbiol. 2018, 56. [Google Scholar] [CrossRef] [Green Version]
  21. Rupprecht, C.E.; Turmelle, A.; Kuzmin, I.V. A Perspective on Lyssavirus Emergence and Perpetuation. Curr. Opin. Virol. 2011, 1, 662–670. [Google Scholar] [CrossRef]
  22. Banyard, A.C.; Evans, J.S.; Luo, T.R.; Fooks, A.R. Lyssaviruses and Bats: Emergence and Zoonotic Threat. Viruses (Basel) 2014, 6, 2974–2990. [Google Scholar] [CrossRef] [Green Version]
  23. Johnson, N.; Vos, A.; Freuling, C.; Tordo, N.; Fooks, A.R.; Muller, T. Human Rabies Due to Lyssavirus Infection of Bat Origin. Vet. Microbiol. 2010, 142, 151–159. [Google Scholar] [CrossRef]
  24. Dacheux, L.; Delmas, O.; Bourhy, H. Human Rabies Encephalitis Prevention and Treatment: Progress Since Pasteurs Discovery. IDDT 2011, 11, 251–299. [Google Scholar] [CrossRef] [PubMed]
  25. Nel, L.H.; Markotter, W. Lyssaviruses. Crit. Rev. Microbiol. 2007, 33, 301–324. [Google Scholar] [CrossRef] [PubMed]
  26. Shipley, R.; Wright, E.; Selden, D.; Wu, G.; Aegerter, J.; Fooks, A.R.; Banyard, A.C. Bats and Viruses: Emergence of Novel Lyssaviruses and Association of Bats with Viral Zoonoses in the EU. Trop. Med. Infect. Dis. 2019, 4, 31. [Google Scholar] [CrossRef] [Green Version]
  27. Badrane, H.; Tordo, N. Host Switching in Lyssavirus History from the Chiroptera to the Carnivora Orders. J. Virol. 2001, 75, 8096–8104. [Google Scholar] [CrossRef] [Green Version]
  28. Ng, T.F.F.; Driscoll, C.; Carlos, M.P.; Prioleau, A.; Schmieder, R.; Dwivedi, B.; Wong, J.; Cha, Y.; Head, S.; Breitbart, M.; et al. Distinct Lineage of Vesiculovirus from Big Brown Bats, United States. Emerg. Infect. Dis. 2013, 19, 1978–1980. [Google Scholar] [CrossRef] [Green Version]
  29. Xu, L.; Wu, J.; Jiang, T.; Qin, S.; Xia, L.; Li, X.; He, B.; Tu, C. Molecular Detection and Sequence Characterization of Diverse Rhabdoviruses in Bats, China. Virus Res. 2018, 244, 208–212. [Google Scholar] [CrossRef]
  30. Binger, T.; Annan, A.; Drexler, J.F.; Müller, M.A.; Kallies, R.; Adankwah, E.; Wollny, R.; Kopp, A.; Heidemann, H.; Dei, D.; et al. A Novel Rhabdovirus Isolated from the Straw-Colored Fruit Bat Eidolon Helvum, with Signs of Antibodies in Swine and Humans. J. Virol. 2015, 89, 4588–4597. [Google Scholar] [CrossRef] [Green Version]
  31. Goldberg, T.L.; Bennett, A.J.; Kityo, R.; Kuhn, J.H.; Chapman, C.A. Kanyawara Virus: A Novel Rhabdovirus Infecting Newly Discovered Nycteribiid Bat Flies Infesting Previously Unknown Pteropodid Bats in Uganda. Sci. Rep. 2017, 7. [Google Scholar] [CrossRef] [Green Version]
  32. Lelli, D.; Prosperi, A.; Moreno, A.; Chiapponi, C.; Gibellini, A.M.; De Benedictis, P.; Leopardi, S.; Sozzi, E.; Lavazza, A. Isolation of a Novel Rhabdovirus from an Insectivorous Bat (Pipistrellus Kuhlii) in Italy. Virol. J. 2018, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Walker, P.J.; Firth, C.; Widen, S.G.; Blasdell, K.R.; Guzman, H.; Wood, T.G.; Paradkar, P.N.; Holmes, E.C.; Tesh, R.B.; Vasilakis, N. Evolution of Genome Size and Complexity in the Rhabdoviridae. PLoS Pathog. 2015, 11, e1004664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Walker, P.J.; Blasdell, K.R.; Calisher, C.H.; Dietzgen, R.G.; Kondo, H.; Kurath, G.; Longdon, B.; Stone, D.M.; Tesh, R.B.; Tordo, N.; et al. ICTV Virus Taxonomy Profile: Rhabdoviridae. J. Gen. Virol. 2018, 99, 447–448. [Google Scholar] [CrossRef]
  35. Spiropoulou, C.F.; Nichol, S.T. A Small Highly Basic Protein Is Encoded in Overlapping Frame within the P Gene of Vesicular Stomatitis Virus. J. Virol. 1993, 67, 3103–3110. [Google Scholar] [CrossRef] [Green Version]
  36. Chenik, M.; Chebli, K.; Blondel, D. Translation Initiation at Alternate In-Frame AUG Codons in the Rabies Virus Phosphoprotein MRNA Is Mediated by a Ribosomal Leaky Scanning Mechanism. J. Virol. 1995, 69, 707–712. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Walker, P.J.; Dietzgen, R.G.; Joubert, D.A.; Blasdell, K.R. Rhabdovirus Accessory Genes. Virus Res. 2011, 162, 110–125. [Google Scholar] [CrossRef] [PubMed]
  38. Fort, P.; Albertini, A.; Van-Hua, A.; Berthomieu, A.; Roche, S.; Delsuc, F.; Pasteur, N.; Capy, P.; Gaudin, Y.; Weill, M. Fossil Rhabdoviral Sequences Integrated into Arthropod Genomes: Ontogeny, Evolution, and Potential Functionality. Mol. Biol. Evol. 2012, 29, 381–390. [Google Scholar] [CrossRef]
  39. Kuzmin, I.V.; Novella, I.S.; Dietzgen, R.G.; Padhi, A.; Rupprecht, C.E. The Rhabdoviruses: Biodiversity, Phylogenetics, and Evolution. Infect. Genet. Evol. 2009, 9, 541–553. [Google Scholar] [CrossRef]
  40. Sapkal, G.N.; Sawant, P.M.; Mourya, D.T. Chandipura Viral Encephalitis: A Brief Review. TOVJ 2018, 12, 44–51. [Google Scholar] [CrossRef]
  41. Grard, G.; Fair, J.N.; Lee, D.; Slikas, E.; Steffen, I.; Muyembe, J.-J.; Sittler, T.; Veeraraghavan, N.; Ruby, J.G.; Wang, C.; et al. A Novel Rhabdovirus Associated with Acute Hemorrhagic Fever in Central Africa. PLoS Pathog. 2012, 8, e1002924. [Google Scholar] [CrossRef]
  42. Rodríguez, L.L. Emergence and Re-Emergence of Vesicular Stomatitis in the United States. Virus Res. 2002, 85, 211–219. [Google Scholar] [CrossRef]
  43. Dacheux, L.; Dommergues, L.; Chouanibou, Y.; Doméon, L.; Schuler, C.; Bonas, S.; Luo, D.; Maufrais, C.; Cetre-Sossah, C.; Cardinale, E.; et al. Co-circulation and Characterization of Novel African Arboviruses (Genus Ephemerovirus) in Cattle, Mayotte Island, Indian Ocean, 2017. Transbound. Emerg. Dis. 2019, 66, 2601–2604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Fischer, M.; Freuling, C.M.; Müller, T.; Schatz, J.; Rasmussen, T.B.; Chriel, M.; Balkema-Buschmann, A.; Beer, M.; Hoffmann, B. Identification of Rhabdoviral Sequences in Oropharyngeal Swabs from German and Danish Bats. Virol. J. 2014, 11, 196. [Google Scholar] [CrossRef] [Green Version]
  45. Aznar-Lopez, C.; Vazquez-Moron, S.; Marston, D.A.; Juste, J.; Ibanez, C.; Berciano, J.M.; Salsamendi, E.; Aihartza, J.; Banyard, A.C.; McElhinney, L.; et al. Detection of Rhabdovirus Viral RNA in Oropharyngeal Swabs and Ectoparasites of Spanish Bats. J. Gen. Virol. 2013, 94, 69–75. [Google Scholar] [CrossRef] [Green Version]
  46. Hu, S.-C.; Hsu, C.-L.; Lee, M.-S.; Tu, Y.-C.; Chang, J.-C.; Wu, C.-H.; Lee, S.-H.; Ting, L.-J.; Tsai, K.-R.; Cheng, M.-C.; et al. Lyssavirus in Japanese Pipistrelle, Taiwan. Emerg. Infect. Dis. 2018, 24, 782–785. [Google Scholar] [CrossRef] [Green Version]
  47. Liu, Y.; Zhang, S.; Zhao, J.; Zhang, F.; Hu, R. Isolation of Irkut Virus from a Murina Leucogaster Bat in China. PLoS Negl. Trop. Dis. 2013, 7, e2097. [Google Scholar] [CrossRef] [Green Version]
  48. Jiang, Y.; Wang, L.; Lu, Z.; Xuan, H.; Han, X.; Xia, X.; Zhao, F.; Tu, C. Seroprevalence of Rabies Virus Antibodies in Bats from Southern China. Vector-Borne Zoonotic Dis. 2010, 10, 177–181. [Google Scholar] [CrossRef]
  49. Guillén-Servent, A.; Francis, C.M. A New Species of Bat of the Hipposideros Bicolor Group (Chiroptera: Hipposideridae) from Central Laos, with Evidence of Convergent Evolution with Sundaic Taxa. Acta Chiropterologica 2006, 8, 39–61. [Google Scholar] [CrossRef]
  50. Dacheux, L.; Larrous, F.; Lavenir, R.; Lepelletier, A.; Faouzi, A.; Troupin, C.; Nourlil, J.; Buchy, P.; Bourhy, H. Dual Combined Real-Time Reverse Transcription Polymerase Chain Reaction Assay for the Diagnosis of Lyssavirus Infection. PLoS Negl. Trop. Dis. 2016, 10, e0004812. [Google Scholar] [CrossRef] [Green Version]
  51. Pallandre, L.; Luo, D.; Feuvrier, C.; Lieffrig, F.; Pozet, F.; Dacheux, L.; Bigarré, L. Revisiting the Classification of Percid Perhabdoviruses Using New Full-Length Genomes. Viruses 2020, 12, 649. [Google Scholar] [CrossRef]
  52. Mareuil, F.; Doppelt-Azeroual, O.; Ménager, H. A Public Galaxy Platform at Pasteur Used as an Execution Engine for Web Services. F1000Research 2017, 1030. [Google Scholar] [CrossRef]
  53. Troupin, C.; Dacheux, L.; Tanguy, M.; Sabeta, C.; Blanc, H.; Bouchier, C.; Vignuzzi, M.; Duchene, S.; Holmes, E.C.; Bourhy, H. Large-Scale Phylogenomic Analysis Reveals the Complex Evolutionary History of Rabies Virus in Multiple Carnivore Hosts. PLoS Pathog. 2016, 12, e1006041. [Google Scholar] [CrossRef] [PubMed]
  54. Dacheux, L.; Cervantes-Gonzalez, M.; Guigon, G.; Thiberge, J.-M.; Vandenbogaert, M.; Maufrais, C.; Caro, V.; Bourhy, H. A Preliminary Study of Viral Metagenomics of French Bat Species in Contact with Humans: Identification of New Mammalian Viruses. PLoS ONE 2014, 9, e0087194. [Google Scholar] [CrossRef] [PubMed]
  55. Milne, I.; Stephen, G.; Bayer, M.; Cock, P.J.; Pritchard, L.; Cardle, L.; Shaw, P.D.; Marshall, D. Using Tablet for Visual Exploration of Second-Generation Sequencing Data. Brief. Bioinform. 2013, 14, 193–202. [Google Scholar] [CrossRef] [PubMed]
  56. Delmas, O.; Holmes, E.C.; Talbi, C.; Larrous, F.; Dacheux, L.; Bouchier, C.; Bourhy, H. Genomic Diversity and Evolution of the Lyssaviruses. PLoS ONE 2008, 3, e0002057. [Google Scholar] [CrossRef]
  57. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X Version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [Green Version]
  58. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  59. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [Green Version]
  60. Guo, H.; Hu, B.-J.; Yang, X.-L.; Zeng, L.-P.; Li, B.; Ouyang, S.; Shi, Z.-L. Evolutionary Arms Race between Virus and Host Drives Genetic Diversity in Bat SARS Related Coronavirus Spike Genes. J. Virol. 2020, 94, e00902-20. [Google Scholar] [CrossRef]
  61. Metselaar, D.; Williams, M.C.; Simpson, D.I.; West, R.; Mutere, F.A. Mount Elgon Bat Virus: A Hitherto Underscribed Virus from Rhinolophus Hildebrandtii Eloquens K. Anderson. Arch. Gesamte Virusforsch. 1969, 26, 183–193. [Google Scholar] [CrossRef]
  62. Shi, Z.L. Emerging Infectious Diseases Associated with Bat Viruses. Sci. China-Life Sci. 2013, 56, 678–682. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Wibbelt, G.; Moore, M.S.; Schountz, T.; Voigt, C.C. Emerging Diseases in Chiroptera: Why Bats? Biol. Lett. 2010, 6, 438–440. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Velasco-Villa, A.; Mauldin, M.R.; Shi, M.; Escobar, L.E.; Gallardo-Romero, N.F.; Damon, I.; Olson, V.A.; Streicker, D.G.; Emerson, G. The History of Rabies in the Western Hemisphere. Antivir. Res. 2017, 146, 221–232. [Google Scholar] [CrossRef]
  65. Jiang, Z.G.; Liu, S.Y.; Wu, Y.; Jiang, X.L.; Zhou, K.Y. China’s Mammal Diversity. Biodivers. Sci. 2017, 25, 886–895. [Google Scholar] [CrossRef] [Green Version]
  66. Hu, B.; Zeng, L.P.; Yang, X.L.; Ge, X.Y.; Zhang, W.; Li, B.; Xie, J.Z.; Shen, X.R.; Zhang, Y.Z.; Wang, N.; et al. Discovery of a Rich Gene Pool of Bat SARS-Related Coronaviruses Provides New Insights into the Origin of SARS Coronavirus. PLoS Pathog. 2017, 13, 27. [Google Scholar] [CrossRef] [PubMed]
  67. Banyard, A.C.; Hayman, D.; Johnson, N.; McElhinney, L.; Fooks, A.R. Bats and Lyssaviruses. Adv. Virus Res. 2011, 79, 239–289. [Google Scholar] [CrossRef] [PubMed]
  68. Aréchiga Ceballos, N.; Vázquez Morón, S.; Berciano, J.M.; Nicolás, O.; Aznar López, C.; Juste, J.; Rodríguez Nevado, C.; Aguilar Setién, A.; Echevarría, J.E. Novel Lyssavirus in Bat, Spain. Emerg. Infect. Dis. 2013, 19, 793–795. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Nokireki, T.; Tammiranta, N.; Kokkonen, U.-M.; Kantala, T.; Gadd, T. Tentative Novel Lyssavirus in a Bat in Finland. Transbound. Emerg. Dis. 2018, 65, 593–596. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Gunawardena, P.S.; Marston, D.A.; Ellis, R.J.; Wise, E.L.; Karawita, A.C.; Breed, A.C.; McElhinney, L.M.; Johnson, N.; Banyard, A.C.; Fooks, A.R. Lyssavirus in Indian Flying Foxes, Sri Lanka. Emerg. Infect. Dis. 2016, 22, 1456–1459. [Google Scholar] [CrossRef]
  71. Serra-Cobo, J. European Bat Lyssavirus Infection in Spanish Bat Populations. Emerg. Infect. Dis. 2002, 8, 413–420. [Google Scholar] [CrossRef]
  72. de Thoisy, B.; Bourhy, H.; Delaval, M.; Pontier, D.; Dacheux, L.; Darcissac, E.; Donato, D.; Guidez, A.; Larrous, F.; Lavenir, R.; et al. Bioecological Drivers of Rabies Virus Circulation in a Neotropical Bat Community. PLoS Negl. Trop. Dis. 2016, 10, e0004378. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Schatz, J.; Ohlendorf, B.; Busse, P.; Pelz, G.; Dolch, D.; Teubner, J.; Encarnação, J.A.; Hoffmann, B.; Kwasnitschka, L.; Balkema-Buschmann, A.; et al. Twenty Years of Active Bat Rabies Surveillance in Germany: A Detailed Analysis and Future Perspectives. Epidemiol. Infect. 2014, 142, 1155–1166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Iwasaki, T.; Inoue, S.; Tanaka, K.; Sato, Y.; Morikawa, S.; Hayasaka, D.; Moriyama, M.; Ono, T.; Kanai, S.; Yamada, A.; et al. Characterization of Oita Virus 296/1972 of Rhabdoviridae Isolated from a Horseshoe Bat Bearing Characteristics of Both Lyssavirus and Vesiculovirus. Arch. Virol. 2004, 149, 1139–1154. [Google Scholar] [CrossRef] [PubMed]
  75. Blasdell, K.R.; Guzman, H.; Widen, S.G.; Firth, C.; Wood, T.G.; Holmes, E.C.; Tesh, R.B.; Vasilakis, N.; Walker, P.J. Ledantevirus: A Proposed New Genus in the Rhabdoviridae Has a Strong Ecological Association with Bats. Am. J. Trop. Med. Hyg. 2015, 92, 405–410. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Li, C.X.; Shi, M.; Tian, J.H.; Lin, X.D.; Kang, Y.J.; Chen, L.J.; Qin, X.C.; Xu, J.; Holmes, E.C.; Zhang, Y.Z. Unprecedented Genomic Diversity of RNA Viruses in Arthropods Reveals the Ancestry of Negative-Sense RNA Viruses. eLife 2015, 4. [Google Scholar] [CrossRef]
  77. Bennett, A.; Paskey, A.; Kuhn, J.; Bishop-Lilly, K.; Goldberg, T. Diversity, Transmission, and Cophylogeny of Ledanteviruses (Rhabdoviridae: Ledantevirus) and Nycteribiid Bat Flies Parasitizing Angolan Soft-Furred Fruit Bats in Bundibugyo District, Uganda. Microorganisms 2020, 8, 750. [Google Scholar] [CrossRef]
  78. Mao, X.; He, G.; Zhang, J.; Rossiter, S.J.; Zhang, S. Lineage Divergence and Historical Gene Flow in the Chinese Horseshoe Bat (Rhinolophus Sinicus). PLoS ONE 2013, 8, e56786. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Repartition of the bat samples analyzed in this study and collected in 21 provinces of China (in grey and back) during the period 2006–2018, with the total number of samples indicated per province. The two provinces in black (Hainan and Hubei) indicate the presence of positive samples. Magnifications of these provinces are presented, with the position and the name of the positive sample caves or sites indicated with black dots, as well as the different rhabdovirus identified (TYBV: Taiyi bat virus, YSBV: Yinshui bat virus, and QZBV: Qiongzhong bat virus).
Figure 1. Repartition of the bat samples analyzed in this study and collected in 21 provinces of China (in grey and back) during the period 2006–2018, with the total number of samples indicated per province. The two provinces in black (Hainan and Hubei) indicate the presence of positive samples. Magnifications of these provinces are presented, with the position and the name of the positive sample caves or sites indicated with black dots, as well as the different rhabdovirus identified (TYBV: Taiyi bat virus, YSBV: Yinshui bat virus, and QZBV: Qiongzhong bat virus).
Viruses 13 00064 g001
Figure 2. Schematic genomic organization of the six Chinese bat rhabdoviruses. The grey arrows represent the five canonical open reading frames (ORFs) (N, P, M, G, and L), and the black arrows indicate the position of the putative additional ORFs. The nucleotide and amino acid lengths of each ORF are indicated.
Figure 2. Schematic genomic organization of the six Chinese bat rhabdoviruses. The grey arrows represent the five canonical open reading frames (ORFs) (N, P, M, G, and L), and the black arrows indicate the position of the putative additional ORFs. The nucleotide and amino acid lengths of each ORF are indicated.
Viruses 13 00064 g002
Figure 3. Phylogenetic classification of the sic Chinese bat rhabdoviruses. A maximum likelihood phylogenetic tree was done using MEGA7.0 on the full amino-acid sequence of the L protein including 155 previously reported rhabdoviruses from GenBank, using the LG+G+I+F model and with 1000 bootstrap replicates. Bat-related rhabdoviruses are indicated by a dedicated symbol. Unclassified rhabdoviruses are indicated by asterisks. The Chinese bat rhabdoviruses described in this study are indicated in grey. All bootstrap proportion values (BSP) > 80% are specified. Scale bar indicates nucleotide substitutions per site.
Figure 3. Phylogenetic classification of the sic Chinese bat rhabdoviruses. A maximum likelihood phylogenetic tree was done using MEGA7.0 on the full amino-acid sequence of the L protein including 155 previously reported rhabdoviruses from GenBank, using the LG+G+I+F model and with 1000 bootstrap replicates. Bat-related rhabdoviruses are indicated by a dedicated symbol. Unclassified rhabdoviruses are indicated by asterisks. The Chinese bat rhabdoviruses described in this study are indicated in grey. All bootstrap proportion values (BSP) > 80% are specified. Scale bar indicates nucleotide substitutions per site.
Viruses 13 00064 g003
Figure 4. Phylogenetic classification of the six Chinese bat rhabdoviruses. A maximum likelihood phylogenetic tree was done with PhyML3.0 on the nucleotide complete genome sequences, including representative members of the Vesiculovirus genus (a) and of the Ledantevirus genus (b), using the GTR+G+I model with a bootstrap of 1000 replicates. The main animal reservoirs for each virus are indicated by specific cartoons and the Chinese bat rhabdoviruses described in this study are indicated in grey. All bootstrap proportion values (BSP) > 80% are specified. Scale bar indicates nucleotide substitutions per site.
Figure 4. Phylogenetic classification of the six Chinese bat rhabdoviruses. A maximum likelihood phylogenetic tree was done with PhyML3.0 on the nucleotide complete genome sequences, including representative members of the Vesiculovirus genus (a) and of the Ledantevirus genus (b), using the GTR+G+I model with a bootstrap of 1000 replicates. The main animal reservoirs for each virus are indicated by specific cartoons and the Chinese bat rhabdoviruses described in this study are indicated in grey. All bootstrap proportion values (BSP) > 80% are specified. Scale bar indicates nucleotide substitutions per site.
Viruses 13 00064 g004
Figure 5. Identification of the nucleotide and amino acid mutations between the four isolates (1017, D170001, D170022, and D170190) of putative YSBV species. The nucleotide mutations are indicated in bold and the amino acid mutations in grey.
Figure 5. Identification of the nucleotide and amino acid mutations between the four isolates (1017, D170001, D170022, and D170190) of putative YSBV species. The nucleotide mutations are indicated in bold and the amino acid mutations in grey.
Viruses 13 00064 g005
Table 1. Distribution of the different bat families tested in this study, according to the type of sample.
Table 1. Distribution of the different bat families tested in this study, according to the type of sample.
Bat FamilyAll SampleSaliva SampleBrain Sample
No. (%)Species No.No. (%)Species No.No. (%)Species No.
Vespertilionidae1388 (30.3)331145 (32.4)28243 (23.3)20
Rhinolophidae1243 (27.2)14926 (26.2)8317 (30.4)12
Hipposideridae946 (20.7)5689 (19.5)5257 (24.6)5
Pteropidae1345 (7.5)5188 (5.3)4157 (15)5
Miniopteridae362 (7.9)3293 (8.3)269 (6.6)3
Emballonuridae210 (4.6)1210 (6)100
Megadermatidae79 (1.7)178 (2.2)100
Molossidae3 (0.1)13 (0.1)11 (0.1)1
Total4576 (100)633532 (100)501044 (100)46
1 Bats belonging to this family are frugivorous; all other bat species tested in this study are insectivorous.
Table 2. Detection rate of rhabdovirus in brain samples of Rhinolophus sinicus and Rhinolophus affinis originating from Hubei and Hainan provinces.
Table 2. Detection rate of rhabdovirus in brain samples of Rhinolophus sinicus and Rhinolophus affinis originating from Hubei and Hainan provinces.
Rhinolophus affinis
(Hainan, Qiongzhong)
Rhinolophus sinicus
(Hubei, Taiyi)
Rhinolophus sinicus
(Hubei, Yinshui)
Total (%)
Date of collection
(month/year)
No. positive/No. tested (%)
08/2006 0/4 0/4 (0)
09/2006 0/11 0/11 (0)
05/2007 2/24 (8.3)
958, 1017 1
2/24 (8.3)
08/20071/9 (11.1)
1127
1/9 (11.1)
04/2009 0/3 0/3 (0)
05/2010 0/1 0/1 (0)
04/2011 0/2 0/2 (0)
03/2017 1/14 (7.1)
D170190
2/50 (4)
D170001, D170022
3/64 (4.7)
09/2017 0/10 0/10 (0)
10/2017 0/13 0/13 (0)
05/2018 0/25 0/25 (0)
Total (%)1/9 (11.1)3/107 (2.8)2/50 (4)6/166 (3.6)
1 Identification of positive samples is indicated in italic.
Table 3. Virus RNA quantification of the six Chinese bat rhabdoviruses in different organs from the original infected bats.
Table 3. Virus RNA quantification of the six Chinese bat rhabdoviruses in different organs from the original infected bats.
Log10 Viral RNA Copies per g of Tissue
SamplesBrainHeartLiverSpleenLungKidneyIntestine
9588.06NA 1NANANANA5.67
101710.47NANANANANA8.05
D17000110.7611.7811.2112.2111.7811.229.06
D17002210.31NANANANANANA
D1701908.729.208.87NA10.019.65NA
112710.22NANANANANANA
1 NA: not applicable.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Luo, D.-S.; Li, B.; Shen, X.-R.; Jiang, R.-D.; Zhu, Y.; Wu, J.; Fan, Y.; Bourhy, H.; Hu, B.; Ge, X.-Y.; et al. Characterization of Novel Rhabdoviruses in Chinese Bats. Viruses 2021, 13, 64. https://0-doi-org.brum.beds.ac.uk/10.3390/v13010064

AMA Style

Luo D-S, Li B, Shen X-R, Jiang R-D, Zhu Y, Wu J, Fan Y, Bourhy H, Hu B, Ge X-Y, et al. Characterization of Novel Rhabdoviruses in Chinese Bats. Viruses. 2021; 13(1):64. https://0-doi-org.brum.beds.ac.uk/10.3390/v13010064

Chicago/Turabian Style

Luo, Dong-Sheng, Bei Li, Xu-Rui Shen, Ren-Di Jiang, Yan Zhu, Jia Wu, Yi Fan, Hervé Bourhy, Ben Hu, Xing-Yi Ge, and et al. 2021. "Characterization of Novel Rhabdoviruses in Chinese Bats" Viruses 13, no. 1: 64. https://0-doi-org.brum.beds.ac.uk/10.3390/v13010064

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop