Next Article in Journal
Membrane Retention of West Nile Virus NS5 Depends on NS1 or NS3 for Enzymatic Activity
Previous Article in Journal
Novel Pestiviruses Detected in Cattle Interfere with Bovine Viral Diarrhea Virus Diagnostics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bat Rabies in the Americas: Is Myotis the Main Ancestral Spreader?

by
Diego A. Caraballo
1,*,
María Lorena Vico
2,
María Guadalupe Piccirilli
3,
Stella Maris Hirmas Riade
3,
Susana Russo
4,
Gustavo Martínez
3,
Fernando J. Beltrán
5 and
Daniel M. Cisterna
2,*
1
Instituto de Ecología, Genética y Evolución de Buenos Aires, CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Autónoma de Buenos Aires C1428EHA, Argentina
2
Departamento de Zoonosis Urbanas, Ministerio de Salud de la Provincia de Buenos Aires, Buenos Aires B1870, Argentina
3
Servicio de Neurovirosis, Instituto Nacional de Enfermedades Infecciosas, Administración Nacional de Laboratorios e Institutos de Salud (ANLIS), “Dr. Carlos G. Malbrán”, Ciudad Autónoma de Buenos Aires C1282AFF, Argentina
4
Dirección General de Laboratorio y Control Técnico (DILAB), Servicio Nacional de Sanidad y Calidad Agroalimentaria (SENASA), Buenos Aires B1640CZT, Argentina
5
Instituto de Zoonosis “Luis Pasteur”, Ciudad Autónoma de Buenos Aires C1405DCD, Argentina
*
Authors to whom correspondence should be addressed.
Viruses 2024, 16(8), 1302; https://doi.org/10.3390/v16081302
Submission received: 8 July 2024 / Revised: 5 August 2024 / Accepted: 15 August 2024 / Published: 16 August 2024
(This article belongs to the Section Animal Viruses)

Abstract

:
The rabies virus (RABV) is the exclusive lyssavirus affecting both wild and domestic mammalian hosts in the Americas, including humans. Additionally, the Americas stand out as the sole region where bat rabies occurs. While carnivore rabies is being increasingly managed across the region, bats are emerging as significant reservoirs of RABV infection for humans and domestic animals. Knowledge of the bat species maintaining rabies and comprehending cross-species transmission (CST) and host shift processes are pivotal for directing surveillance as well as ecological research involving wildlife reservoir hosts. Prior research indicates that bat RABV CST is influenced by host genetic similarity and geographic overlap, reflecting host adaptation. In this study, we compiled and analyzed a comprehensive nucleoprotein gene dataset representing bat-borne RABV diversity in Argentina and the broader Americas using Bayesian phylogenetics. We examined the association between host genus and geography, finding both factors shaping the global phylogenetic structure. Utilizing a phylogeographic approach, we inferred CST and identified key bat hosts driving transmission. Consistent with CST determinants, we observed monophyletic/paraphyletic clustering of most bat genera in the RABV phylogeny, with stronger CST evidence between host genera of the same family. We further discuss Myotis as a potential ancestral spreader of much of RABV diversity.

1. Introduction

Rabies (RABV) is the most widespread and oldest known lyssavirus infecting mammals, both domestic and wild, including humans. RABV (genus: Lyssavirus, family: Rhabdoviridae) causes acute, highly lethal encephalomyelitis and is responsible for numerous local epidemics, resulting in approximately 59,000 human deaths annually, primarily in Africa and Asia [1]. While Rhabdoviruses can infect a very wide range of hosts, including plants, invertebrates, and vertebrates via vectors, as well as vertebrates without vectors, the lyssaviruses are unique in that they are the only group known to exclusively infect mammals, without exception [2]. Bats are the primary hosts of most lyssaviruses, but rabies is unique in that it circulates not only among mammals of the order Chiroptera but also among a diverse set of Carnivora [3]. Occasional spillover infections occur in other animal species, including humans and domestic animals.
Phylogenetically, RABV can be divided into two primary groups: (i) one comprising viruses sampled from bats and carnivores in the Americas, referred to as the ‘bat-related’ virus group, and (ii) another primarily consisting of ‘dog-related’ viruses sampled from carnivores, excluding bats [4]. The ‘bat-related’ RABV group is limited to viruses circulating among various New World bat species, as well as some mesocarnivores such as skunks and raccoons, forming a distinct phylogenetic cluster within this group [3,4]. In Latin America, significant attention has been directed towards vampire bats as hosts, particularly due to their role as reservoir species for RABV, capable of infecting other mammals, including humans and livestock [5,6]. Nevertheless, there is a significant occurrence of spillovers, including the establishment of variants, from insectivorous bats to other mammals. This has resulted in the conceivable establishment of RABV variants in various species, not only in carnivores such as skunks, raccoons, foxes, coatis, and kinkajous but also in marmosets and presumably in capuchin monkeys [7,8,9,10]. These spillover events also carry significant public health implications, as demonstrated by a case where the virus was transmitted from Tadarida brasiliensis to a cat, resulting in a fatal human infection in Argentina through secondary transmission in 2021 [11].
The main determinants of bat RABV cross-species transmissions (CST) have been shown to be host genetic similarity and, to a lesser extent, geographic overlap between species from different genera [12,13]. In these studies, based on North American bat RABV, higher host genetic similarity correlated with a notably increased occurrence of both spillovers (in external nodes) and host shifts (limited to internal nodes). This is attributable to the assumption that highly related hosts will require less viral adaptation due to finding a similar physiological and immunological environment [14].
An expectable outcome of host adaptation is that virus occurrence will be constrained by host phylogeny. This was demonstrated by the non-random clustering of 34 RABV host species in the carnivore phylogeny [15]. This would imply that, in a phylogenetic reconstruction of RABV variants, host taxa will form monophyletic/paraphyletic clusters.
In this study, we generated a dataset representing the diversity of bat-borne RABV in Argentina to contextualize it within the broader spectrum of bat RABV diversity in the Americas through a Bayesian phylogenetic approach. We further analyze the association between host genus and geography, revealing that both factors contribute to the global phylogenetic structure. Employing a phylogeographic approach, we aimed to model CST and determine which bat hosts act as main spreaders. Based on the background on CST determinants and adaptation, we predict that (i) most bat genera will be monophyletic/paraphyletic in the RABV phylogeny, indicative of viral adaptation and (ii) evidence for CST will be higher between host genera of the same family. We further analyze the role of Myotis as the potential ancestral spreader of the majority of RABV diversity observed at the present.

2. Methods

2.1. Samples

A total of 62 brain samples from RABV-positive insectivorous bats were obtained from 2001 to 2023 from several Argentinian provinces (Figure 1, Table S1): Buenos Aires (50), Chubut (1), Córdoba (7), La Pampa (1), Río Negro (2), and Tierra del Fuego (1). These samples belong to the following bat species: Tadarida brasiliensis (50), Myotis dinelli (3), Myotis levis (3), Myotis sp. (1), Eptesicus furinalis (2), Dasypterus ega (1), Aeorestes cinereus (1), and Eumops bonariensis (1). These specimens were gathered as part of rabies virus surveillance conducted by the Avellaneda Zoonosis Center in Buenos Aires province (BAI), the Luis Pasteur Zoonosis Institute in the Autonomous City of Buenos Aires (CABA), and the National Agri-Food Health and Quality Service (SENASA) in Buenos Aires province.

2.2. RNA Extraction, PCR Amplification, and Sequencing

Viral RNA was obtained using the QIAmp Viral RNA Mini kit, Qiagen, following the manufacturer’s instructions. Reverse transcription and PCR amplification were performed using random primers for retrotranscription and primers Lys001, 550B, JW12, 550F, and 304R to obtain the complete RABV nucleoprotein gene, using PCR conditions as described in [16]. The product obtained was purified using the ExoSAP-IT reagent, Usb; Nucleotide sequences were obtained using the BigDye Terminator 3.1 reagent (Applied Biosystems, Waltham, MA, USA), purified with the XTerminator reagent (Applied Biosystems), and sequenced using the ABI 3500 genetic analyzer (Applied Biosystems). The sequences generated were deposited in GenBank under the accession numbers PP661584-PP661645 (Table S1).

2.3. Identification of Bat RABV Main Lineages

An initial analysis was conducted to identify all major bat-borne RABV (rabies virus) lineages. Subsequently, a subsample was selected to encompass all primary host lineages and sublineages, as well as the countries where these lineages were sampled. This set the stage for a secondary analysis aimed at testing the phylogenetic positioning of Argentinian samples and conducting host-switching analysis.
Reference sequences for the initial phylogenetic analysis were retrieved from the NCBI virus database (https://www.ncbi.nlm.nih.gov/labs/virus/vssi/, accessed on 1 April 2023) [17], using the following filters: Virus: Lyssavirus rabies (taxid: 11292), Host: Chiroptera (bats) (taxid: 9397), Protein: Nucleoprotein, and filtered for samples from the Americas. The search yielded 1864 hits. These sequences were aligned with Clustal Omega and subjected to an initial analysis on IQ-TREE [18], using the GTR + F + R5 model selected with the corrected Akaike Information Criterion (AICc). This preliminary analysis enabled the distinction of 22 RABV main host lineages (Figure S1). A subset of complete nucleoprotein sequences that represent all main lineages/sublineages (main host species) and all countries was selected for testing the sequences obtained in this study (Table S1). The validity of host species and genus names was verified using the Integrated Taxonomic Information System (ITIS) online database, www.itis.gov, CC0 (https://doi.org/10.5066/F7KH0KBK, accessed on 1 February 2024) and updated to reflect current taxonomy when necessary.

2.4. Bayesian Phylogeny and Cross-Species Transmission Dynamics

The complete matrix consisted of 132 sequences (62 generated in this study plus 70 N sequences from bats representing all main lineages, see File S1). A Bayesian phylogeny was run using BEAST 1.10.4 [19]. The molecular substitution model was selected with MrModeltest v2 using the second-order Akaike information criterion (AICc) [20], separating each codon position (GTR + I + G, HKY + I + G, and GTR + G, for 1st, 2nd, and 3rd codon positions, respectively). A genus-level host-switching dynamics analysis was performed using an asymmetric model, allowing for unequal transition rates. This method enables the reconstruction of ancestral host posterior probabilities at each internal node of the tree and estimates the transition rates between states (hosts). A tip-dated strict molecular clock was utilized, alongside a strict clock for the host genus. A Bayesian stochastic search variable selection (BSSVS) method was applied to assess the significance of transitions between trait states, using Bayes factors (BFs) as a measure of statistical significance. BFs were computed using SpreaD3 [21], and their interpretation followed Jeffreys’ guidelines [22]: (BF > 3: substantial support, BF > 10: strong support, BF > 30: very strong support, BF > 100: decisive support). Only transitions strongly supported by BFs are shown in figures, consistent with methodologies employed in prior studies. A coalescent with a constant-size tree prior was selected. A total of 1E8 Markov chain Monte Carlo (MCMC) generations were run, sampling every 1E4 generations, discarding the first 10% of the run (burn-in phase). Convergence of the MCMC runs was confirmed using Tracer. Trees were visualized in FigTree v1.4.4 [23] and plotted with ggtree v3.6.2 [24] and ggplot2 version 3.4.2 [25], using R version 4.2.3 [26] and RStudio version 2023.03.0 Build 386 [27].
The decision to base our analysis at the genus level has multiple reasons. Primarily, it allows us to reflect profound host shifting events, particularly those occurring between genera. This approach allows for a comprehensive analysis of the contribution of each bat genus, rather than focusing solely on species. Analyzing at the species level may lead to inflated rates of CST by potentially overlooking common variants shared among closely related species. This approach could also diminish the significance of a donor bat genus’s contribution to a recipient one. Additionally, our choice is grounded in the frequent challenges associated with species-level identification in bats. The prevalence of cryptic species and the potential for human error in identification underscore the practicality and reliability of analyzing at the genus level [28].

2.5. Association Tests

Statistical tests were conducted to examine trait clustering among the tips of the phylogeny. Bayesian Tip Association Significance Testing (BaTS version 0.9 beta) software was utilized for this purpose [29]. We assessed the clustering patterns of the bat genera and countries sampled. To evaluate the significance of trait clustering, we compared the association index (AI), parsimony score (PS), and monophyletic clade (MC) statistics derived from posterior samples of each run (with a burn-in of 10%) against null distributions created from 100 randomizations of traits assigned to tips across the sampled trees. A significance level of p < 0.05 was chosen for this analysis.

2.6. Recombination Analysis

To identify candidate recombinants, the alignment including 132 bat-related nucleoprotein sequences was tested for recombination using the Recombination Detection Program RDP version 4 (RDP4) [30]. Default settings for all algorithms were applied (RDP, GENECONV, Maxchi, Chimaera, BootScan, SiSscan, and 3Seq), incorporating SiSscan “primary scans” (the default option is that this method is only used to examine sequences in which recombination signals are detectable by other “primary scanning” methods). For each inferred event, a phylogenetic analysis was run in IQ-TREE using the recombinant segments to test for topological incongruence, estimating node support using an approximate likelihood-ratio test (1000 replicates).
An additional analysis to detect recombination was conducted using the GARD method [31], on the Datamonkey 2.0 web server [32]. The analysis employed Beta-Gamma site-to-site rate variation and utilized four classes to model rate heterogeneity.

2.7. Genetic Distances

Uncorrected pairwise genetic p-distances were calculated with MEGA X [33]. Distances were calculated by considering all sequences belonging to a specific lineage, defined as members of either a monophyletic or paraphyletic group within a given host genus. Some genera comprised multiple lineages, each treated separately. Furthermore, distances between all sequences from all lineages within the genera Eptesicus, Myotis, and Tadarida were calculated to evaluate overall variation within these widely distributed genera represented by multiple lineages in the tree. In the distance calculation of Myotis yumanensis, two sequences were excluded (KY203184 and KY203142), since these were considered in the calculations of the lineage “Myotis NA” (see below). Similarly, one sequence of Myotis austroriparius (AY039225) was excluded from the intralineage distance, because it pertains to a non-related lineage with respect to the other two sequences.

3. Results

3.1. Recombination Analysis

Detecting recombinants prior to phylogenetic analysis ensures the accuracy of evolutionary reconstructions by preventing the inclusion of sequences that may introduce conflicting signals and distort inferred relationships. Three putative recombination events were detected with RDP4. Event 1 is supported by four out of seven methods and involves two sequences from Tadarida (major parent and recombinant) and one from Eumops in the Aeorestes clade. The recombinant fragment spans nucleotides 718 to 828 (111 bp in length). It appears to be a true recombinant because the sequence is identical to the minor parent in the putative recombinant region (15 substitutions in 111 sites), as confirmed by a quick phylogenetic test conducted with IQ-TREE. However, this sequence was not generated in our analysis (GenBank accession number: JF826130), and thus, we cannot confirm if it was produced in error (contamination, assembly error, etc.). Even if it is a true recombinant, the involvement of such a short fragment causes this sequence to remain in the parental clade when the complete alignment is included. The second event involves two sequences from the Myotis SA clade (major parent and recombinant) and a sequence from Tadarida SA. The breakpoints of this event are not unambiguously identified, and in the recombinant fragment, there is no clear affinity to the minor parent (it shares multiple interspersed substitutions, comparable in number, with both parental sequences). This potential event is more likely to represent a divergent sequence rather than a recombinant. The third event lacks a clear signal of recombination and was detected only by one method (SiScan); thus, it is not considered a true recombination event.
In addition, GARD found no evidence of recombination. A total of 570 potential breakpoints were identified, of which only one was confirmed as a breakpoint. However, in the topological analysis, the multiple tree model could not be preferred over the single tree model by an evidence ratio of 100 or greater, indicating the inferred breakpoint may reflect rate variation instead of topological incongruence.
Taken together, the lack of evidence of recombination in the dataset allowed us to proceed with the phylogenetic inference and downstream analyses.

3.2. Molecular Phylogeny

A total of 62 N gene whole sequences were generated in this study. No premature stop codons were found in the dataset. Although there is low support near the base of the tree, there are four well-supported groups that encompass the diversity of bat RABV (Figure 2). File S2 contains the annotated tree, allowing exploration of additional data such as sample country annotations and support values for internal nodes. It is interesting to note that except for Myotis, Eptesicus, and Dasypterus, all bat genera are lineage-specific, occurring in a single monophyletic or paraphyletic cluster in the tree. The most-represented lineage is that of Tadarida brasiliensis from South America (Tadarida SA), which comprises 54 closely related sequences, including samples from Brazil, Argentina, Chile, and Uruguay. As confirmed in the preliminary analysis (Figure S1), Tadarida is paraphyletic with respect to Desmodus/Artibeus, being the North American variant (Tadarida NA) sister to this group. Similarly, the Lasiurus/Dasypterus/Aeorestes clade is paraphyletic to a clade containing variants circulating in Perimyotis and Lasionycteris. It is noteworthy that there are sequences from both subcontinents in this clade, suggesting that these bats are connected through migration. However, it is remarkable that the structure of this clade, comprising three closely related genera, is more influenced by geographical factors (country) than by the host species (File S2). Another well-supported clade is that of Histiotus (and Eptesicus) from Brazil and Argentina, which is paraphyletic to a variant found in Nyctinomops. This clade is related to four lineages found in Corynorhinus townsendii, Myotis volans, M. yumanensis, and Eptesicus fuscus, from the USA and Canada. Closely related Eptesicus and Myotis variants are also found in Brazil and Argentina, including RABV sequences obtained from Myotis levis and Myotis dinelli, and those obtained from E. furinalis. There is an additional Myotis clade from Canada and the USA which is paraphyletic with respect to Parastrellus and Antrozous. Finally, there is a third Eptesicus clade from the USA (Eptesicus NA2).

3.3. Association Tests

In the phylogeny, there is a clear hierarchical structure that responds to two levels: host and geography. Both global tests, PS and AI, indicate that the association exists at both levels. This suggests that related RABV sequences are more likely to share a host or geographic location by ancestry than they do by chance. At the host level, the genera that contribute to stronger host–virus associations (producing larger monophyletic clades) are Tadarida and Myotis, although there are several other genera with significant associations (Table S2). At the geographic level, Argentina, the USA, and Canada are characterized by the strongest virus–country associations (Table S3).

3.4. Cross-Species Transmission Dynamics

Two types of CST can be distinguished: those that took place in the recent past, and likely represent dead-ends in the transmission chain (spillovers), and those that can be inferred from ancestral states reconstruction, and that reflect the adaptation of the virus to a new host (Figure 2). Among the sequences generated in this study, an AeorestesEumops spillover was confirmed alongside two MyotisEptesicus events in independent clades of the tree (Figure 2). One occurred in a clade sister to the Histiotus lineage, while the other took place in the clade Myotis SA. Remarkably, both spillovers involved M. dinelli and E. furinalis, which may be indicative of an ecological association between these species.
The ancestral host reconstruction poses Aeorestes/Dasypterus as the most probable common ancestor of bat rabies. However, it should be taken into account that the relationships between the main clades are not strongly supported, and this result should be then taken cautiously.
The Bayesian stochastic search variable selection (BSSVS) procedure was employed to identify viral host switches between bat genera along the branches of the MCC tree, further calculating Bayes factors (BF) to estimate the significance of these switches. Fourteen highly supported (BF < 10) inter-genus host switches were identified (Figure 3). The highest supported switches (BF > 100) occurred from Myotis to Eptesicus, and from Desmodus to Artibeus. Myotis is the most frequent donor, with significant host switches to five other genera. This approach enables the identification of historical transmission chains such as TadaridaDesmodusArtibeus, and Myotis/EptesicusHistiotusNyctinomops. In addition, host shifts in both directions were found between Aeorestes/Dasypterus and Lasionycteris/Perimyotis. This could be indicative of recent host divergence (Aeorestes/Dasypterus) or ecological association, in the case of less-related taxa.

3.5. Genetic Distances

The estimates of evolutionary divergence between sequences varied significantly across different lineages (Figure 4, Table S4). Among the 22 lineages analyzed, 14 exhibited low levels of divergence (<0.025 substitutions per site). Certain lineages, such as Desmodus, Eptesicus NA2, and Myotis NA, displayed intermediate levels of genetic divergence. Although these values and their distribution may be influenced by the limited number of samples in some lineages, it is noteworthy that Histiotus and Parastrellus, despite their low representation in the dataset, exhibited high levels of genetic divergence. This suggests that these lineages experienced elevated substitution rates and/or reduced homogenization. In contrast, the lineage Tadarida SA, which possesses the largest representation in the dataset (comprising 54 sequences from four countries), displayed low levels of genetic differentiation, indicating significant homogenization among these bats.
Furthermore, lineages such as Eptesicus NA1 and Myotis SA displayed significant dispersion, where closely related sequences formed sublineages; however, these sublineages exhibited considerable divergence from each other. The highest levels of divergence, accompanied by significant dispersion, were observed in Dasypterus, with distance values exceeding 0.1 substitutions per site. This reflects both the lack of monophyly within this lineage and its extensive geographical distribution, spanning from South America to North America.
In the comparison of genetic distances at the genus level, the three genera analyzed showed contrasting patterns (Figure 5). Eptesicus has high intra-genus distances (mean 0.1 substitutions per site), reflecting the phylogenetic distance between its four non-related clades. Eptesicus NA2 originated from a host shift from Myotis and is more related to Eptesicus SA than to the other Eptesicus lineage from North America (even when Eptesicus NA1 and NA2 both circulate among E. fuscus, while Eptesicus SA circulates in E. furinalis). Myotis also depicts high intra-genus distances, in concordance with the phylogenetic diversity observed in RABV variants circulating in this genus. On the contrary, Tadarida exhibits a distinct bimodal pattern, indicative of the separation between North American and South American clades. The majority of distance values correspond to intra-clade differences (Tadarida SA), resulting in an exceptionally low mean of 0.02 substitutions per site.

4. Discussion

Surveying virus diversity is crucial for understanding viral ecology and developing effective strategies for disease prevention and control. In this study, we obtained and analyzed full nucleoprotein sequences of bat-related RABV in Argentina. Previous studies involved short fragments (191–320 bp) that, although allowing for variant distinction, were less resolutive at internal nodes of the phylogeny [34,35,36]. In addition, this study focuses on all available bat RABV sequences, encompassing the entirety of known virus diversity to date, and explicitly conducts a CST analysis. Research on CST has predominantly centered on lineages originating from North America, neglecting other geographical regions [12,13,37,38]. The single prior study examining the diversity across both subcontinents conducted a phylogenetic analysis utilizing partial N gene sequences (582 bp) [39]. This research was limited in its representation of several lineages, which were identified in the comprehensive phylogenetic analysis conducted in the present study. As a result, these additional lineages were included in both the phylogeny and the CST analysis. Therefore, the present study integrates more robust information (full N gene and all known bat RABV lineages) from both subcontinents, providing insights, particularly into widespread host taxa.
The phylogeny revealed a clear hierarchical structure that responds to two levels: host and geography. This result was further supported by the tip-association tests. The majority of bat genera are lineage-specific, confirming our first prediction. Interestingly, not only geographically restricted genera (e.g., Histiotus, Antrozous) but also three widespread genera found across the three American subcontinents, Tadarida, Desmodus, and Artibeus, represent this situation. Tadarida is divided into two related subclades, South America and North America, while Desmodus and Artibeus share RABV sequences that represent a unique monophyletic variant. This result was also confirmed by a phylogeographic study of vampire bat-associated RABV that revealed that the viruses found in bats from both North and South America originate from a single common ancestor, indicating that they constitute a singular, well-defined epizootic event [40].
In contrast to strict host specificity, the presence of multiple lineages within a host genus reflects dynamics dominated by host switching. We observe that most widespread host genera cluster in different non-related RABV lineages, suggesting that when species/populations are geographically isolated, there is more chance of variant extinction and recolonization by a less-adapted variant. The presence of multiple viral lineages within the widespread host genera Eptesicus, Dasypterus, and Myotis a priori (but see below for Myotis) suggests that their species acquired RABV infections independently, leading to simultaneous epizootics in different regions. Eptesicus is divided into three clades, two of which represent the North American E. fuscus, while the other corresponds to the South American E. furinalis. As occurs in Myotis, in North America, Eptesicus has two non-related cycles (both are related to Myotis). In the case of Dasypterus, it is interesting to note that its variants are interspersed among other closely related genera (Aeorestes and Lasiurus), and this may be indicative of ecological and phylogenetic proximity.
The presence of four distinct Myotis lineages in North America (Figure 2) might initially suggest multiple independent acquisitions of RABV through host shifts, as proposed in prior studies [13,37,39]. However, upon considering the ancestral host reconstruction and CST dynamics inferred from the BSSVS analysis, a different interpretation emerges. The major bat RABV clade, which encompasses all Myotis lineages alongside most host lineages, except for Dasypterus/Lasiurus/Aeorestes (plus Lasionycteris and Perimyotis), indicates Myotis as the common ancestor of these lineages. In the host shifting analysis (Figure 3), Myotis emerges as the predominant donor and never as a (significant) recipient of RABV. Taken together, these results suggest that Myotis did not acquire RABV variants from other bat genera; rather, its numerous species and populations harbored most of the observed RABV diversity in bats today and transferred RABV to other bat genera. The occurrence of multiple RABV lineages in Myotis aligns with its status as an ancient (>33 Mya) and globally distributed genus, encompassing over 120 described species [41]. Indeed, these bats host a wide variety of viruses, such as coronaviruses, astroviruses, and picornaviruses, and they have been identified as key players in the diversification of these pathogens [42,43,44]. The diversification of the genus Myotis resulted from a rapid speciation process [41], which was shaped by the isolation and reconnection of populations during the climatic fluctuations of the Quaternary, ultimately making it the most speciose mammal genus globally [45]. Throughout this evolutionary journey, accompanying pathogens probably co-diversified with Myotis, undergoing similar isolation and reconnection dynamics that elucidate the diversity of RABV observed in present times. It is likely that Myotis acquired RABV after the divergence of the New World and Old World clades, about 19 MYA [46] or that it was lost in the Old World lineage. Both possibilities are plausible for explaining the absence of bat RABV in the Old World.
Two contrasting patterns of genetic distances used as a proxy of phylogenetic diversity were found (Figure 4 and Figure 5). Some lineages, such as Histiotus and Parastrellus, depicted high levels of genetic divergence, suggesting that they could be experiencing high rates of substitution rates or reduced homogenization. However, these values are based on a few sequences, and a more detailed analysis should be carried out. In contrast, the lineage Tadarida SA, which possesses the largest representation in the dataset (comprising 54 sequences from four countries), displayed low levels of genetic differentiation, indicating significant homogenization among these bats, probably due to high levels of gene flow that characterizes this species [47].
Previous studies of North American bats across various genera revealed that the occurrence of sustained transmission and spillover of RABV increased with higher host phylogenetic relatedness [13,48], likely due to the lower number of adaptive changes necessary for viral exchange. However, between closely related species (within a genus), no significant relationship between host species relatedness and the frequency of CST events was found [37]. At this level, where all host species are highly related, ecological factors, such as the level of overlap in geographic ranges, gain preponderance. In the present study, we confirmed a predominance of host shifts among members of the same bat family, reflecting the existence of evolutionary constraints at this taxonomic level.
The two highest supported CST were MyotisEptesicus, and DesmodusArtibeus. It is worth noting that both pairs of CST were unilateral, where the inferred donor hosts are Myotis and Desmodus, and there are multiple host shifts to the recipient hosts Eptesicus and Artibeus, supported also by the broader phylogenetic analysis (Figure S1). A significant association between Myotis and Eptesicus is evident, as all Myotis lineages show an associated Eptesicus lineage [49]. While previous claims suggested that Desmodus and Artibeus maintained independent RABV cycles, clustered in reciprocally monophyletic lineages [50,51], the approach employed in this study suggests a more intricate scenario. The ancestral host is likely Desmodus which received RABV from North American Tadarida. Artibeus might have acquired the virus via CST from Desmodus twice, but numerous spillbacks ArtibeusDesmodus were also observed (see Figure S1). Since we aimed to include representatives of main host lineages and not spillovers to capture host shifts, in our CST approach we included only one member of each Artibeus lineage (Figure 2), so we could not quantify the level of spillbacks to Desmodus.
The results presented in this study confirm the general paradigm that CST of viruses is more likely to occur among genetically similar taxa, suggesting that genetic relatedness of hosts may imply similar immune and physiological environments, reducing barriers to adaptation to a new host [14]. However, there is considerable evidence demonstrating that bat-related RABVs have repeatedly jumped to several terrestrial and arboreal New World mammal species, including a wide range of mesocarnivores (cats, dogs, foxes, skunks, coatis, and kinkajous) but also primates (marmosets and capuchin monkeys) [7,8,9,52]. Some of these events led to the establishing of independent transmission cycles, in which these mammalian species also became RABV reservoirs. The raccoon (Procyon lotor) strain of the rabies virus has a shared ancestor with a skunk strain, and both are most closely related to rabies viruses found in bats [53]. In white-nosed coatis, a variant derived from Tadarida brasiliensis, has also been demonstrated to represent a stable host shift [54], while in marmosets, RABV underwent a stable host shift likely from bats of the genus Lasiurus [55]. Even in cases where stable host shifts have not been confirmed, it is common to observe repeated host shifts from the same donor to the same recipient host species, involving the same RABV variant. This phenomenon, together with the finding of pre-shift convergent mutations led to the hypothesis that significant pre-adaptation of bat viruses should take place before effective host shifting to terrestrial mammals [9]. The factors determining the success in bat-to-carnivore spillovers and stable host shifting are far from being elucidated but may include social behavior, lifespan, reproductive cycle duration, population density, and home range, among others [56]. The potential of wild and domestic mesocarnivores to be successful RABV reservoirs, in addition to the host shifts and spillovers cited above, should be an alert for the need of rapid diagnostics and typing in RABV surveillance for effective disease prevention and mitigation. Considering that bat rabies is non eradicable, in countries where canine rabies has been completely (or almost) completely eliminated, these bat-to-mesocarnivore spillovers may be the main cause of risk of spillover to humans.
In this study, we investigated the diversity of bat RABV from Argentina, placing it in the context of the total diversity known to date. We confirmed our first hypothesis that most bat hosts are clustered in single RABV lineages. Among widespread genera, we identified two patterns: those that acquired independent RABV variants multiple times (e.g., Eptesicus and Dasypterus) and Myotis, proposed as the ancestral host genus of the major bat RABV clade and the primary donor to other genera. In addition, we corroborated our second prediction that the majority of CST would occur between related bat genera (same family) and underscored the strong association between Myotis/Eptesicus and Desmodus/Artibeus, supported by both host shifts and spillovers. Although the panorama is far from complete—for example, there is a large gap in molecular data from surveillance in Central America—the identification of main hosts and the characterization of CST can be useful for better understanding RABV epizootiology and preparedness for viral emergence.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/v16081302/s1, Table S1. Sample list. Samples, host species, geographic origin, and collection dates of the sequences used in this study, as well as those retrieved from Genbank; Table S2. Tip association statistics for hosts genus; Table S3. Tip association statistics for hosts and geographic location; Table S4. Pairwise genetic distances; Figure S1. Maximum likelihood phylogenetic tree of rabies virus (RABV) sequences from GenBank. The maximum likelihood phylogenetic tree was constructed using all available sequences of the rabies virus (RABV) obtained from GenBank. The tree was generated to assess the overall diversity of RABV and to identify main lineages for subsequent phylogenetic and host-switching analyses, including test samples. The branch lengths are proportional to genetic distances, and the bootstrap values greater than 90% are indicated at the branch nodes to indicate the robustness of the tree topology; File S1. Alignment of the 1353 bp encompassing the full nucleoprotein gene; File S2. Annotated phylogenetic tree in nexus format.

Author Contributions

Conceptualization, D.A.C. and D.M.C.; methodology, D.A.C., M.L.V., M.G.P., S.M.H.R., S.R., G.M., F.J.B. and D.M.C.; formal analysis, D.A.C.; investigation, D.A.C. and D.M.C.; resources, D.A.C. and D.M.C.; data curation, D.A.C.; writing—original draft preparation, D.A.C.; writing—review and editing, D.A.C. and D.M.C.; visualization, D.A.C.; supervision, D.A.C. and D.M.C.; project administration, D.A.C. and D.M.C.; funding acquisition, D.A.C. and D.M.C. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by grant PIP-0587 awarded by the National Scientific and Technical Research Council of Argentina to Diego Caraballo and the Instituto Nacional de Enfermedades Infecciosas (INEI-ANLIS), “Dr. Carlos G. Malbrán” to Daniel Cisterna.

Institutional Review Board Statement

Ethical review and approval were waived for this study, due to the fact that all the samples belong to third-party euthanized animals which were submitted to zoonotic surveillance laboratories.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within this article or Supplementary Materials.

Acknowledgments

This work used computational resources from CCAD—Universidad Nacional de Córdoba (https://ccad.unc.edu.ar/, accessed on 1 April 2023), which are part of SNCAD—MinCyT, República Argentina.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hampson, K.; Coudeville, L.; Lembo, T.; Sambo, M.; Kieffer, A.; Attlan, M.; Barrat, J.; Blanton, J.D.; Briggs, D.J.; Cleaveland, S.; et al. Estimating the Global Burden of Endemic Canine Rabies. Carvalho MS, editor. PLoS Negl. Trop. Dis. 2015, 9, e0003709. [Google Scholar] [CrossRef]
  2. Holmes, E.C.; Harvey, E.H. The Diversity, Evolution and Emergence of Rabies Virus in the Americas. In History of Rabies in the Americas: From the Pre-Columbian to the Present, Volume I; Rupprecht, C.E., Ed.; Springer International Publishing: Cham, Switzeraland, 2023; pp. 43–59. [Google Scholar] [CrossRef]
  3. Fisher, C.R.; Streicker, D.G.; Schnell, M.J. The spread and evolution of rabies virus: Conquering new frontiers. Nat. Rev. Microbiol. 2018, 16, 241–255. [Google Scholar] [CrossRef] [PubMed]
  4. 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]
  5. Benavides, J.A.; Rojas Paniagua, E.; Hampson, K.; Valderrama, W.; Streicker, D.G. Quantifying the burden of vampire bat rabies in Peruvian livestock. PLoS Negl. Trop. Dis. 2017, 11, e0006105. [Google Scholar] [CrossRef] [PubMed]
  6. Johnson, N.; Aréchiga-Ceballos, N.; Aguilar-Setien, A. Vampire Bat Rabies: Ecology, Epidemiology and Control. Viruses 2014, 6, 1911–1928. [Google Scholar] [CrossRef]
  7. Aréchiga-Ceballos, N.; Velasco-Villa, A.; Shi, M.; Flores-Chávez, S.; Barrón, B.; Cuevas-Domínguez, E.; González-Origel, A.; Aguilar-Setién, A. New rabies virus variant found during an epizootic in white-nosed coatis from the Yucatan Peninsula. Epidemiol. Infect. 2010, 138, 1586–1589. [Google Scholar] [CrossRef]
  8. Dell’Armelina Rocha, P.R.; Velasco-Villa, A.; de Lima, E.M.; Salomoni, A.; Fusaro, A.; da Conceição Souza, E.; Negreiros, R.L.; Zafino, V.L.; Zamperin, G.; Leopardi, S.; et al. Unexpected rabies variant identified in kinkajou (Potos flavus), Mato Grosso, Brazil. Emerg. Microbes Infect. 2020, 9, 851–854. [Google Scholar] [CrossRef]
  9. Kuzmin, I.V.; Shi, M.; Orciari, L.A.; Yager, P.A.; Velasco-Villa, A.; Kuzmina, N.A.; Streicker, D.G.; Bergman, D.L.; Rupprecht, C.E. Molecular Inferences Suggest Multiple Host Shifts of Rabies Viruses from Bats to Mesocarnivores in Arizona during 2001–2009. PLoS Pathog. 2012, 8, e1002786. [Google Scholar] [CrossRef] [PubMed]
  10. 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]
  11. Firpo, S.; Piccirilli, M.G.; Urizar, R.; Vitta, N.; Hirmas Riade, S.M.; Leguizamón, C.; Vico, M.L.; Martínez, G.; Beltrán, F.J.; Cisterna, D.M. Human Rabies by Secondary Transmission in Argentina, 2021. Diseases 2022, 10, 17. [Google Scholar] [CrossRef]
  12. Faria, N.R.; Suchard, M.A.; Rambaut, A.; Streicker, D.G.; Lemey, P. Simultaneously reconstructing viral cross-species transmission history and identifying the underlying constraints. Philos. Trans. R. Soc. B Biol. Sci. 2013, 368, 20120196. [Google Scholar] [CrossRef] [PubMed]
  13. Streicker, D.G.; Turmelle, A.S.; Vonhof, M.J.; Kuzmin, I.V.; McCracken, G.F.; Rupprecht, C.E. Host Phylogeny Constrains Cross-Species Emergence and Establishment of Rabies Virus in Bats. Science 2010, 329, 676–679. [Google Scholar] [CrossRef]
  14. Mollentze, N.; Biek, R.; Streicker, D.G. The role of viral evolution in rabies host shifts and emergence. Curr. Opin. Virol. 2014, 8, 68–72. [Google Scholar] [CrossRef] [PubMed]
  15. Huang, S.; Bininda-Emonds, O.R.P.; Stephens, P.R.; Gittleman, J.L.; Altizer, S. Phylogenetically related and ecologically similar carnivores harbour similar parasite assemblages. J. Anim. Ecol. 2014, 83, 671–680. [Google Scholar] [CrossRef]
  16. Torres, C.; Lema, C.; Dohmen, F.G.; Beltran, F.; Novaro, L.; Russo, S.; Freire, M.C.; Velasco-Villa, A.; Mbayed, V.A.; Cisterna, D.M. Phylodynamics of vampire bat-transmitted rabies in Argentina. Mol. Ecol. 2014, 23, 2340–2352. [Google Scholar] [CrossRef]
  17. NCBI Virus; National Library of Medicine (US), National Center for Biotechnology Information: Bethesda, MD, USA, 2004.
  18. Nguyen, L.-T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef] [PubMed]
  19. Suchard, M.A.; Lemey, P.; Baele, G.; Ayres, D.L.; Drummond, A.J.; Rambaut, A. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 2018, 4, 1–5. [Google Scholar] [CrossRef]
  20. Nylander, J. Mr Modeltest 2.2. Program Distributed by the Author; Evolutionary Biology Centre, Uppsala University: Uppsala, Sweden, 2004. [Google Scholar]
  21. Bielejec, F.; Baele, G.; Vrancken, B.; Suchard, M.A.; Rambaut, A.; Lemey, P. SpreaD3: Interactive Visualization of Spatiotemporal History and Trait Evolutionary Processes. Mol. Biol. Evol. 2016, 33, 2167–2169. [Google Scholar] [CrossRef]
  22. Jeffreys, H. Theory of Probability; Clarendon Press: Oxford, UK, 1961. [Google Scholar]
  23. Rambaut, A. FigTree-v1. 4.2. 2014. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 1 April 2023).
  24. Yu, G. Using ggtree to Visualize Data on Tree-Like Structures. Curr. Protoc. Bioinforma. 2020, 69, e96. [Google Scholar] [CrossRef]
  25. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016; Available online: https://ggplot2.tidyverse.org (accessed on 1 April 2023).
  26. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.R-project.org/ (accessed on 1 April 2023).
  27. RStudio Team. RStudio: Integrated Development Environment for R; RStudio, PBC: Boston, MA, USA, 2020; Available online: http://www.rstudio.com/ (accessed on 1 April 2023).
  28. Caraballo, D.A.; Montani, M.E.; Martínez, L.M.; Antoniazzi, L.R.; Sambrana, T.C.; Fernández, C.; Cisterna, D.M.; Beltrán, F.J.; Colombo, V.C. Heterogeneous taxonomic resolution of cytochrome b gene identification of bats from Argentina: Implications for field studies. PLoS ONE 2020, 15, e0244750. [Google Scholar] [CrossRef]
  29. Parker, J.; Rambaut, A.; Pybus, O.G. Correlating viral phenotypes with phylogeny: Accounting for phylogenetic uncertainty. Infect. Genet. Evol. 2008, 8, 239–246. [Google Scholar] [CrossRef]
  30. Martin, D.P.; Murrell, B.; Golden, M.; Khoosal, A.; Muhire, B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015, 1, vev003. [Google Scholar] [CrossRef]
  31. Kosakovsky Pond, S.L.; Posada, D.; Gravenor, M.B.; Woelk, C.H.; Frost, S.D.W. GARD: A genetic algorithm for recombination detection. Bioinformatics 2006, 22, 3096–3098. [Google Scholar] [CrossRef] [PubMed]
  32. Weaver, S.; Shank, S.D.; Spielman, S.J.; Li, M.; Muse, S.V.; Kosakovsky Pond, S.L. Datamonkey 2.0: A Modern Web Application for Characterizing Selective and Other Evolutionary Processes. Mol. Biol. Evol. 2018, 35, 773–777. [Google Scholar] [CrossRef]
  33. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef] [PubMed]
  34. Caraballo, D.A.; Lombardo, M.A.; Becker, P.; Sabio, M.S.; Lema, C.; Martínez, L.M.; Beltrán, F.J.; Li, Y.; Cisterna, D.M. Evaluation of Two Real-Time, TaqMan Reverse Transcription-PCR Assays for Detection of Rabies Virus in Circulating Variants from Argentina: Influence of Sequence Variation. Viruses 2020, 13, 23. [Google Scholar] [CrossRef] [PubMed]
  35. Cisterna, D.; Bonaventura, R.; Caillou, S.; Pozo, O.; Andreau, M.L.; Dalla Fontana, L.; Echegoyen, C.; De Mattos, C.; De Mattos, C.; Russo, S.; et al. Antigenic and molecular characterization of rabies virus in Argentina. Virus Res. 2005, 109, 139–147. [Google Scholar] [CrossRef] [PubMed]
  36. Pinero, C.; Gury Dohmen, F.; Beltran, F.; Martinez, L.; Novaro, L.; Russo, S.; Palacios, G.; Cisterna, D.M. High Diversity of Rabies Viruses Associated with Insectivorous Bats in Argentina: Presence of Several Independent Enzootics. PLoS Negl. Trop. Dis. 2012, 6, e1635. [Google Scholar] [CrossRef]
  37. Jacquot, M.; Wallace, M.A.; Streicker, D.G.; Biek, R. Geographic Range Overlap Rather than Phylogenetic Distance Explains Rabies Virus Transmission among Closely Related Bat Species. Viruses 2022, 14, 2399. [Google Scholar] [CrossRef]
  38. Streicker, D.G.; Altizer, S.M.; Velasco-Villa, A.; Rupprecht, C.E. Variable evolutionary routes to host establishment across repeated rabies virus host shifts among bats. Proc. Natl. Acad. Sci. USA 2012, 109, 19715–19720. [Google Scholar] [CrossRef]
  39. Nadin-Davis, S.; Alnabelseya, N.; Knowles, M.K. The phylogeography of Myotis bat-associated rabies viruses across Canada. PLoS Negl. Trop. Dis. 2017, 11, e0005541. [Google Scholar] [CrossRef]
  40. Streicker, D.G.; Fallas González, S.L.; Luconi, G.; Barrientos, R.G.; Leon, B. Phylodynamics reveals extinction–recolonization dynamics underpin apparently endemic vampire bat rabies in Costa Rica. Proc. R. Soc. B Biol. Sci. 2019, 286, 20191527. [Google Scholar] [CrossRef]
  41. Gunnell, G.F.; Smith, R.; Smith, T. 33 million year old Myotis (Chiroptera, Vespertilionidae) and the rapid global radiation of modern bats. PLoS ONE 2017, 12, e0172621. [Google Scholar] [CrossRef]
  42. Armero, A.; Li, R.; Bienes, K.M.; Chen, X.; Li, J.; Xu, S.; Chen, Y.; Hughes, A.C.; Berthet, N.; Wong, G. Myotis fimbriatus Virome, a Window to Virus Diversity and Evolution in the Genus Myotis. Viruses 2022, 14, 1899. [Google Scholar] [CrossRef] [PubMed]
  43. Caraballo, D.A.; Sabio, M.S.; Colombo, V.C.; Piccirilli, M.G.; Vico, L.; Hirmas Riade, S.M.; Campos, J.; Martínez, G.; Beltrán, F.; Baumeister, E.; et al. The Role of Molossidae and Vespertilionidae in Shaping the Diversity of Alphacoronaviruses in the Americas. Microbiol. Spectr. 2022, 10, e03143-22. [Google Scholar] [CrossRef] [PubMed]
  44. Latinne, A.; Hu, B.; Olival, K.J.; Zhu, G.; Zhang, L.; Li, H.; Chmura, A.A.; Field, H.E.; Zambrana-Torrelio, C.; Epstein, J.H.; et al. Origin and cross-species transmission of bat coronaviruses in China. Nat. Commun. 2020, 11, 4235. [Google Scholar] [CrossRef] [PubMed]
  45. Novaes, R.L.M.; Cláudio, V.C.; Wilson, D.E.; Weber, M.M.; Weksler, M.; Moratelli, R. Integrative taxonomy and evolutionary history reveal cryptic diversity in Myotis riparius (Chiroptera: Vespertilionidae). Biol. J. Linn. Soc. 2023, 142, 91–112. [Google Scholar] [CrossRef]
  46. Ruedi, M.; Stadelmann, B.; Gager, Y.; Douzery, E.J.P.; Francis, C.M.; Lin, L.-K.; Guillén-Servent, A.; Cibois, A. Molecular phylogenetic reconstructions identify East Asia as the cradle for the evolution of the cosmopolitan genus Myotis (Mammalia, Chiroptera). Mol. Phylog. Evol. 2013, 69, 437–449. [Google Scholar] [CrossRef]
  47. Russell, A.L.; Medellín, R.A.; Mccracken, G.F. Genetic variation and migration in the Mexican free-tailed bat (Tadarida brasiliensis mexicana). Mol. Ecol. 2005, 14, 2207–2222. [Google Scholar] [CrossRef]
  48. Longdon, B.; Brockhurst, M.A.; Russell, C.A.; Welch, J.J.; Jiggins, F.M. The Evolution and Genetics of Virus Host Shifts. PLoS Pathog. 2014, 10, e1004395. [Google Scholar] [CrossRef]
  49. Oliveira, R.D.N.; De Souza, S.P.; Lobo, R.S.V.; Castilho, J.G.; Macedo, C.I.; Carnieli, P.; Fahl, W.O.; Achkar, S.M.; Scheffer, K.C.; Kotait, I.; et al. Rabies virus in insectivorous bats: Implications of the diversity of the nucleoprotein and glycoprotein genes for molecular epidemiology. Virology 2010, 405, 352–360. [Google Scholar] [CrossRef]
  50. de Souza, D.N.; Oliveira, R.N.; Asprino, P.F.; Bettoni, F.; Macedo, C.I.; Achkar, S.M.; Fahl, W.D.O.; Brandão, P.E.; Castilho, J.G. Evolution and divergence of the genetic lineage Desmodus rotundus/Artibeus lituratus of rabies virus in São Paulo State. Arch. Virol. 2023, 168, 266. [Google Scholar] [CrossRef]
  51. Fahl, W.O.; Carnieli Jr, P.; Castilho, J.G.; Carrieri, M.L.; Kotait, I.; Iamamoto, K.; Oliveira, R.N.; Brandão, P.E. Desmodus rotundus and Artibeus spp. bats might present distinct rabies virus lineages. Braz. J. Infect. Dis. 2012, 16, 545–551. [Google Scholar] [CrossRef]
  52. Favoretto, S.R.; De Mattos, C.C.; De Mattos, C.A.; Campos, A.C.A.; Sacramento, D.R.V.; Durigon, E.L. The emergence of wildlife species as a source of human rabies infection in Brazil. Epidemiol. Infect. 2013, 141, 1552–1561. [Google Scholar] [CrossRef] [PubMed]
  53. Kuzmina, N.A.; Lemey, P.; Kuzmin, I.V.; Mayes, B.C.; Ellison, J.A.; Orciari, L.A.; Hightower, D.; Taylor, S.T.; Rupprecht, C.E. The Phylogeography and Spatiotemporal Spread of South-Central Skunk Rabies Virus. PLoS ONE 2013, 8, e82348. [Google Scholar] [CrossRef]
  54. Puebla-Rodríguez, P.; Almazán-Marín, C.; Garcés-Ayala, F.; Rendón-Franco, E.; Chávez-López, S.; Gómez-Sierra, M.; Sandoval-Borja, A.; Martínez-Solís, D.; Escamilla-Ríos, B.; Sauri-González, I.; et al. Rabies virus in white-nosed coatis (Nasua narica) in Mexico: What do we know so far? Front. Vet. Sci. 2023, 10, 1090222. [Google Scholar] [CrossRef] [PubMed]
  55. Kotait, I.; Oliveira, R.D.N.; Carrieri, M.L.; Castilho, J.G.; Macedo, C.I.; Pereira, P.M.C.; Boere, V.; Montebello, L.; Rupprecht, C.E. Non-human primates as a reservoir for rabies virus in Brazil. Zoonoses Public. Health 2019, 66, 47–59. [Google Scholar] [CrossRef]
  56. Worsley-Tonks, K.E.; Escobar, L.E.; Biek, R.; Castaneda-Guzman, M.; Craft, M.E.; Streicker, D.G.; White, L.A.; Fountain-Jones, N.M. Using host traits to predict reservoir host species of rabies virus. PLoS Negl. Trop. Dis. 2020, 14, e0008940. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Sampling localities and bat host species of the RABV N sequences generated in this study.
Figure 1. Sampling localities and bat host species of the RABV N sequences generated in this study.
Viruses 16 01302 g001
Figure 2. Maximum clade credibility tree of the Bayesian phylogenetic analysis of bat RABV representing the whole diversity of the virus. Node sizes are proportional to Bayesian posterior probability (relevant nodes are shown with numeric values). Different colors were assigned to each bat genus. Pie charts represent uncertainty in the ancestral host reconstruction. The timescale is in years before the present.
Figure 2. Maximum clade credibility tree of the Bayesian phylogenetic analysis of bat RABV representing the whole diversity of the virus. Node sizes are proportional to Bayesian posterior probability (relevant nodes are shown with numeric values). Different colors were assigned to each bat genus. Pie charts represent uncertainty in the ancestral host reconstruction. The timescale is in years before the present.
Viruses 16 01302 g002
Figure 3. Graph representing host transition rates among bat genera. Line width indicates the switch significance level according to the Bayes factor test. Genus color indicates host family.
Figure 3. Graph representing host transition rates among bat genera. Line width indicates the switch significance level according to the Bayes factor test. Genus color indicates host family.
Viruses 16 01302 g003
Figure 4. Uncorrected pairwise genetic p-distances at the lineage level. Distances were calculated by considering all sequences belonging to a specific lineage, defined as members of either a monophyletic or paraphyletic group within a given host genus. Some genera comprised multiple lineages which were treated separately.
Figure 4. Uncorrected pairwise genetic p-distances at the lineage level. Distances were calculated by considering all sequences belonging to a specific lineage, defined as members of either a monophyletic or paraphyletic group within a given host genus. Some genera comprised multiple lineages which were treated separately.
Viruses 16 01302 g004
Figure 5. Uncorrected pairwise genetic p-distances of multilineage genera (Eptesicus, Myotis, and Tadarida). Violin plots and box plots show the density of data points, as well as the median (horizontal line). Dark-red dots indicate mean values.
Figure 5. Uncorrected pairwise genetic p-distances of multilineage genera (Eptesicus, Myotis, and Tadarida). Violin plots and box plots show the density of data points, as well as the median (horizontal line). Dark-red dots indicate mean values.
Viruses 16 01302 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Caraballo, D.A.; Vico, M.L.; Piccirilli, M.G.; Hirmas Riade, S.M.; Russo, S.; Martínez, G.; Beltrán, F.J.; Cisterna, D.M. Bat Rabies in the Americas: Is Myotis the Main Ancestral Spreader? Viruses 2024, 16, 1302. https://doi.org/10.3390/v16081302

AMA Style

Caraballo DA, Vico ML, Piccirilli MG, Hirmas Riade SM, Russo S, Martínez G, Beltrán FJ, Cisterna DM. Bat Rabies in the Americas: Is Myotis the Main Ancestral Spreader? Viruses. 2024; 16(8):1302. https://doi.org/10.3390/v16081302

Chicago/Turabian Style

Caraballo, Diego A., María Lorena Vico, María Guadalupe Piccirilli, Stella Maris Hirmas Riade, Susana Russo, Gustavo Martínez, Fernando J. Beltrán, and Daniel M. Cisterna. 2024. "Bat Rabies in the Americas: Is Myotis the Main Ancestral Spreader?" Viruses 16, no. 8: 1302. https://doi.org/10.3390/v16081302

APA Style

Caraballo, D. A., Vico, M. L., Piccirilli, M. G., Hirmas Riade, S. M., Russo, S., Martínez, G., Beltrán, F. J., & Cisterna, D. M. (2024). Bat Rabies in the Americas: Is Myotis the Main Ancestral Spreader? Viruses, 16(8), 1302. https://doi.org/10.3390/v16081302

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