Hybridization hotspots at bat swarming sites

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

Phylogenomic analyses unraveled the evolution of viral tolerance in bats

Bats host a range of viruses, exhibiting a coevolution process with many virus genera and a special capacity for viral tolerance. Foley et al.1 performed phylogenomic analyses for 60 bat species, finding that swarming behavior might facilitate cross-species introgression and the spread of anti-virus immunity gene loci across species.

Prehibernation swarming in temperate bats: a critical transition between summer activity and hibernation

In this contribution to Dr. Brock Fenton's Festschrift, we briefly reflect on Dr. Fenton's seminal works examining bat swarming behaviour in Ontario and use these reflections as a launch pad to conduct a global review on autumn swarming in bats, and underlying hypotheses to explain this behaviour. Our review frames the swarming period as a time of critical transitions, during which bats must balance multiple life history trade-offs, and we consider how various intrinsic and extrinsic factors may contribute to inter- and intraspecific differences in autumn behaviour. We discuss the transition away from summer residency, including maternity colony breakup, day roosting, and migration during autumn. We review key life history elements of swarming, including mating behaviours and associated reproductive condition, genetic exchange during swarming, and variation among sexes, ages, and species. Finally, we discuss the behaviours and physiological states of bats transitioning from the swarming period to hibernation. Throughout, we identify common patterns and also exceptions. Over 50 years of research has yielded many insights into autumn swarming, but knowledge gaps remain. Future research focus on a greater diversity of species will reveal general principles underlying the transition from summer active season, through the swarming period, and into winter hibernation.

Consequences of Hybridization in Mammals: A Systematic Review

Hybridization, defined as breeding between two distinct taxonomic units, can have an important effect on the evolutionary patterns in cross-breeding taxa. Although interspecific hybridization has frequently been considered as a maladaptive process, which threatens species genetic integrity and survival via genetic swamping and outbreeding depression, in some cases hybridization can introduce novel adaptive variation and increase fitness. Most studies to date focused on documenting hybridization events and analyzing their causes, while relatively little is known about the consequences of hybridization and its impact on the parental species. To address this knowledge gap, we conducted a systematic review of studies on hybridization in mammals published in 2010-2021, and identified 115 relevant studies. Of 13 categories of hybridization consequences described in these studies, the most common negative consequence (21% of studies) was genetic swamping and the most common positive consequence (8%) was the gain of novel adaptive variation. The total frequency of negative consequences (49%) was higher than positive (13%) and neutral (38%) consequences. These frequencies are biased by the detection possibilities of microsatellite loci, the most common genetic markers used in the papers assessed. As negative outcomes are typically easier to demonstrate than positive ones (e.g., extinction vs hybrid speciation), they may be over-represented in publications. Transition towards genomic studies involving both neutral and adaptive variation will provide a better insight into the real impacts of hybridization.

Genetic Introgression and Morphological Variation in Naked-Back Bats (Chiroptera: Mormoopidae: Pteronotus Species) along Their Contact Zone in Central America

Two sibling bare-backed bat species (Pteronotus fulvus and P. gymnonotus) have been traditionally differentiated by their size. However, intermediate specimens between the two species have been found in sympatric populations along southern Mexico and it has been suggested that they may be the outcome of a hybridization process between the two species. We used one mitochondrial (COI), three nuclear markers (PRKCL, STAT5A and RAG2) and 13 microsatellites to explore the evolutionary relationships between these two species and elucidate whether the intermediate morphotypes correspond to hybrid individuals. These markers have been analyzed in sympatric and allopatric populations of the two species plus the closely related species Pteronotus davyi. We confirmed the species-level differentiation of the three lineages (P. fulvus, P. davyi and P. gymnonotus), but the phylogenetic hypotheses suggested by the nuclear and mitochondrial markers were discordant. We confirm that the discordance between markers is due to genetic introgression through the mitochondrial capture of P. fulvus in P. gymnonotus populations. Such introgression was found in all P. gymnonotus specimens across its sympatric distribution range (Mexico to Costa Rica) and is related to expansion/retraction species distribution pulses associated with changes in forest distribution during the Quaternary climate cycles. Microsatellite analyses showed contemporary genetic contact between the two sympatric species and 3.0% of the samples studied were identified as hybrids. In conclusion, we found a historical and asymmetric genetic introgression (through mitochondrial capture) of P. fulvus into P. gymnonotus in Mexico and Central America and a limited contemporary gene exchange between the two species. However, no relationship was found between hybridization and the intermediate-sized specimens from southern Mexico, which might likely result from a clinal variation with latitude. These results confirm the need for caution when using forearm size to identify these species in the field and when differentiating them in the laboratory based on mitochondrial DNA alone.

Genetic diversity, population structure, and effective population size in two yellow bat species in south Texas

There are increasing concerns regarding bat mortality at wind energy facilities, especially as installed capacity continues to grow. In North America, wind energy development has recently expanded into the Lower Rio Grande Valley in south Texas where bat species had not previously been exposed to wind turbines. Our study sought to characterize genetic diversity, population structure, and effective population size in Dasypterus ega and D. intermedius, two tree-roosting yellow bats native to this region and for which little is known about their population biology and seasonal movements. There was no evidence of population substructure in either species. Genetic diversity at mitochondrial and microsatellite loci was lower in these yellow bat taxa than in previously studied migratory tree bat species in North America, which may be due to the non-migratory nature of these species at our study site, the fact that our study site is located at a geographic range end for both taxa, and possibly weak ascertainment bias at microsatellite loci. Historical effective population size (NEF) was large for both species, while current estimates of Ne had upper 95% confidence limits that encompassed infinity. We found evidence of strong mitochondrial differentiation between the two putative subspecies of D. intermedius (D. i. floridanus and D. i. intermedius) which are sympatric in this region of Texas, yet little differentiation using microsatellite loci. We suggest this pattern is due to secondary contact and hybridization and possibly incomplete lineage sorting at microsatellite loci. We also found evidence of some hybridization between D. ega and D. intermedius in this region of Texas. We recommend that our data serve as a starting point for the long-term genetic monitoring of these species in order to better understand the impacts of wind-related mortality on these populations over time.

Hidden in our pockets: building of a DNA barcode library unveils the first record of Myotis alcathoe for Portugal

Background

The advent and boom of DNA barcoding technologies have provided a powerful tool for the fields of ecology and systematics. Here, we present the InBIO Barcoding Initiative Database: Portuguese Bats (Chiroptera) dataset containing DNA sequences of 63 specimens representing the 25 bat species currently known for continental Portugal. For that, we sequenced tissues samples obtained in a vast array of projects spanning the last two decades.

New information

We added four new Barcoding Index Numbers (BINs) to existing Chiroptera barcodes on BOLD, two belonging to Myotisescalerai, one to Plecotusauritus and the other to Rhinolophushipposideros. Surprisingly, one of the samples initially identified in the field as Myotismystacinus turned out to be Myotisalcathoe, which represents the first record of this species for Portugal. The presence of Nyctalusnoctula in Portugal was also genetically confirmed for the first time. This case study shows the power and value of DNA barcoding initiatives to unravel new data that may be hidden on biological collections.

Species assignment in forensics and the challenge of hybrids

Forensic identification of species is in growing demand, particularly from law enforcement authorities in the areas of wildlife, fisheries and hunting as well as food authentication. Within the non-human forensic genetics expanding applications' field, the major current difficulties result from the lack of standards and genetic databases as well as the poor or absent taxonomic definition of several groups. Here we focus on a forensically important and overlooked problem in species identification: the exclusive use of uniparental markers, a common practice in current genetic barcoding methodologies, may lead to incorrect or impossible assignment whenever hybrids can occur (frequently, not only in domesticates, but also in the wild). For example, if one of these cases involves a mammal, and mitochondrial DNA alone is used (which in instances may be the only type of DNA sequence available in databases), the sample will be wrongfully assigned to the female parental species, completely missing the detection of a possible hybrid animal. The importance of this issue in the forensic contributions to food authentication, wildlife and conservation genetics is analyzed. We present a cautionary guidance on the forensic reporting of results avoiding this error.

Species delimitation and evolutionary reconstruction within an integrative taxonomic framework: A case study on Rhinolophus macrotis complex (Chiroptera: Rhinolophidae)

Species delimitation and evolutionary reconstruction remain challenging for non-model species that have experienced reticulate evolution and exhibit conflicting patterns of differentiation among multiple lines of evidence, such as mitochondrial and nuclear data and phenotypes. Here, we applied an integrative taxonomic approach to a case study of Rhinolophus macrotis complex, whose taxonomic status remains controversial, to provide insight into the systematics and evolutionary history of these species. By integrating traditional genetic markers with different modes of inheritance, genome-wide SNPs as well as phenotypic characteristics, we clarified the presence of three closely related species, R. episcopus, R. siamensis, and R. osgoodi, within this complex, and proposed a new taxonomic treatment for R. osgoodi. Our results suggested that hybridization and introgression are the main causes of low mtDNA divergence in these species. Combined with the demographic inference, we deduced that glacial-interglacial cycles drove geographic isolation and secondary contacts of these species, then promoted hybridization and lineage fusion among them, finally resulting in a reticulate evolutionary pattern. Overall, our study highlights the importance of combining multiple types of data to delimit species, especially those with conserved morphology, and to reveal the sophisticated processes of speciation.

Coexistence of two sympatric cryptic bat species in French Guiana: insights from genetic, acoustic and ecological data

Background

The distinction between lineages of neotropical bats from the Pteronotus parnellii species complex has been previously made according to mitochondrial DNA, and especially morphology and acoustics, in order to separate them into two species. In these studies, either sample sizes were too low when genetic and acoustic or morphological data were gathered on the same individuals, or genetic and other data were collected on different individuals. In this study, we intensively sampled bats in 4 caves and combined all approaches in order to analyse genetic, morphologic, and acoustic divergence between these lineages that live in the same caves in French Guiana.

Results

A multiplex of 20 polymorphic microsatellite markers was developed using the 454-pyrosequencing technique to investigate for the first time the extent of reproductive isolation between the two lineages and the population genetic structure within lineages. We genotyped 748 individuals sampled between 2010 and 2015 at the 20 nuclear microsatellite loci and sequenced a portion of the cytochrome c oxydase I gene in a subset of these. Two distinct, non-overlapping haplogroups corresponding to cryptic species P. alitonus and P. rubiginosus were revealed, in accordance with previous findings. No spatial genetic structure between caves was detected for both species. Hybridization appeared to be quite limited (0.1–4%) using microsatellite markers whereas introgression was more common (7.5%) and asymmetric for mitochondrial DNA (mtDNA).

Conclusions

The extremely low rate of hybridization could be explained by differences in life cycle phenology between species as well as morphological and acoustical distinction between sexes in one or the other species. Taken together, these results add to our growing understanding of the nature of species boundaries in Pteronotus parnelli, but deserve more in-depth studies to understand the evolutionary processes underlying asymmetric mtDNA introgression in this group of cryptic species.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1289-8) contains supplementary material, which is available to authorized users.

Conflicting Evolutionary Histories of the Mitochondrial and Nuclear Genomes in New World Myotis Bats

The rapid diversification of Myotis bats into more than 100 species is one of the most extensive mammalian radiations available for study. Efforts to understand relationships within Myotis have primarily utilized mitochondrial markers and trees inferred from nuclear markers lacked resolution. Our current understanding of relationships within Myotis is therefore biased towards a set of phylogenetic markers that may not reflect the history of the nuclear genome. To resolve this, we sequenced the full mitochondrial genomes of 37 representative Myotis, primarily from the New World, in conjunction with targeted sequencing of 3648 ultraconserved elements (UCEs). We inferred the phylogeny and explored the effects of concatenation and summary phylogenetic methods, as well as combinations of markers based on informativeness or levels of missing data, on our results. Of the 294 phylogenies generated from the nuclear UCE data, all are significantly different from phylogenies inferred using mitochondrial genomes. Even within the nuclear data, quartet frequencies indicate that around half of all UCE loci conflict with the estimated species tree. Several factors can drive such conflict, including incomplete lineage sorting, introgressive hybridization, or even phylogenetic error. Despite the degree of discordance between nuclear UCE loci and the mitochondrial genome and among UCE loci themselves, the most common nuclear topology is recovered in one quarter of all analyses with strong nodal support. Based on these results, we re-examine the evolutionary history of Myotis to better understand the phenomena driving their unique nuclear, mitochondrial, and biogeographic histories.

Babesia vesperuginis, a neglected piroplasmid: new host and geographical records, and phylogenetic relations

BACKGROUND

Babesia spp. are hemoparasites which infect the red blood cells of a large variety of mammals. In bats, the only known species of the genus is Babesia vesperuginis. However, except a few old reports, the host range and geographical distribution of this bat parasite have been poorly studied. This study aimed to investigate the presence of piroplasms in tissues of bats collected in four different countries from eastern and central Europe: Austria, Czech Republic, Hungary and Romania.

METHODS

A total of 461 bat carcasses (24 species) were collected between 2001 and 2016 from caves, mines and buildings. PCR was performed using specific primers targeting a portion of the 18S rDNA nuclear gene and cytochrome c oxidase subunit 1 mitochondrial gene, followed by sequencing.

RESULTS

The results of this study show for the first time the presence of B. vesperuginis in bats in central and eastern Europe. The phylogenetic analysis of the 18S rDNA nuclear gene revealed no variability between the sequences and the phylogenetic analysis of the cox1 mitochondrial gene proved that B. vesperuginis could be divided into two subclades.

CONCLUSION

Our study showed a broad geographical distribution of B. vesperuginis in European bats, reporting its presence in five new host species (M. cf. alcathoe, M. bechsteinii, M. myotis, Pi. nathusii and V. murinus) and three new countries.

Investigating Hybridization between the Two Sibling Bat Species Myotis myotis and M. blythii from Guano in a Natural Mixed Maternity Colony

Because they can form seasonal mixed-species groups during mating and maternal care, bats are exciting models for studying interspecific hybridization. Myotis myotis and M. blythii are genetically close and morphologically almost identical, but they differ in some aspects of their ecology and life-history traits. When they occur in sympatry, they often form large mixed maternity colonies, in which their relative abundance can vary across time due to a shift in the timing of parturition. For the first time, we used non-invasive genetic methods to assess the hybridization rate and colony composition in a maternity colony of M. myotis and M. blythii located in the French Alps. Bat guano was collected on five sampling dates spread across the roost occupancy period and was analysed for individual genotype. We investigated whether the presence of hybrids followed the pattern of one of the parental species or if it was intermediate. We identified 140 M. myotis, 12 M. blythii and 13 hybrids among 250 samples. Parental species appeared as genetically well-differentiated clusters, with an asymmetrical introgression towards M. blythii. By studying colony parameters (effective size, sex ratio and proportion of the three bat types) across the sampling dates, we found that the abundances of hybrid and M. blythii individuals were positively correlated. Our study provides a promising non-invasive method to study hybridization in bats and raises questions about the taxonomic status of the two Myotis species. We discuss the contribution of this study to the knowledge of hybrid ecology, and we make recommendations for possible future research to better understand the ecology and behaviour of hybrid individuals.

Autumn flight activity of the greater horseshoe bat at hibernacula

The activity patterns of vespertilionid bats at the entrances of hibernacula in temperate zones have been investigated on many occasions. However, little is known about nocturnal flight activity in autumn and its function in horseshoe bats so far. Here, we have recorded the autumn flight activity in and around the entrances of hibernacula of greater horseshoe bats,Rhinolophus ferrumequinumby behavioural observations, mistnetting, and acoustic surveys for two consecutive years. Greater horseshoe bats not only chased each other around cave entrances but also entered and left the entrance frequently. However, no mating events were directly observed in this study, normally a prominent part of swarming behaviour in vespertilionids. Therefore, it is difficult to deduce the function of the autumn flight activity of the greater horseshoe bats at hibernacula from the present data. Additionally, the recorded flight activity was positively correlated with ambient temperature and humidity. To our knowledge, our study provides the first detailed record of the pattern of autumn flight activity in horseshoe bats, and suggests that weather conditions may impact the bats’ nocturnal flight activity.

Speciation processes in putative island endemic sister bat species: false impressions from mitochondrial DNA and microsatellite data

Cases of geographically restricted co-occurring sister taxa are rare and may point to potential divergence with gene flow. The two bat species Murina gracilis and Murina recondita are both endemic to Taiwan and are putative sister species. To test for nonallopatric divergence and gene flow in these taxa, we generated sequences using Sanger and next-generation sequencing, and combined these with microsatellite data for coalescent-based analyses. MtDNA phylogenies supported the reciprocally monophyletic sister relationship between M. gracilis and M. recondita; however, clustering of microsatellite genotypes revealed several cases of species admixture suggesting possible introgression. Sequencing of microsatellite flanking regions revealed that admixture signatures stemmed from microsatellite allele homoplasy rather than recent introgressive hybridization, and also uncovered an unexpected sister relationship between M. recondita and the continental species Murina eleryi, to the exclusion of M. gracilis. To dissect the basis of these conflicts between ncDNA and mtDNA, we analysed sequences from 10 anonymous ncDNA loci with *beast and isolation-with-migration and found two distinct clades of M. eleryi, one of which was sister to M. recondita. We conclude that Taiwan was colonized by the ancestor of M. gracilis first, followed by the ancestor of M. recondita after a period of allopatric divergence. After colonization, the mitochondrial genome of M. recondita was replaced by that of the resident M. gracilis. This study illustrates how apparent signatures of sympatric divergence can arise from complex histories of allopatric divergence, colonization and hybridization, thus highlighting the need for rigorous analyses to distinguish between such scenarios.

Temporal Changes in Body Mass and Body Condition of Cave-Hibernating Bats During Staging and Swarming

The rapid colonization of the Pseudogymnoascus destructans fungus across cave systems in eastern North America and the associated bat mortalities (white-nose syndrome; WNS), necessitates studies of cave-hibernating bats that remain unaffected by, or in close proximity to, the leading edge of the fungal distribution to provide baseline predisturbance data from which to assess changes due to fungal effects. Studies of the physiological ecology of cave-hibernating bats during the spring staging and autumn swarming seasons are few, and an understanding of patterns in body condition of bats associated with entry into and emergence from hibernation is incomplete. We sampled bats at the entrance to a cave in Mammoth Cave National Park, Kentucky, during swarming and staging, prior to (2011 and 2012), concurrent with (2013), and following (2014) the arrival of the WNS fungus. We evaluated seasonal and annual changes in body mass and body condition of bats entering and leaving the cave. We captured 1,232 bats of eight species. Sex ratios of all species were male-biased. Capture success was substantially reduced in 2014, following the second winter after arrival of the WNS fungus. Significant temporal variation in body mass and body mass index was observed for little brown bats Myotis lucifugus, northern long-eared bats M. septentrionalis, and tri-colored bats Perimyotis subflavus, but not Indiana bats M. sodalis. Little brown bats and northern long-eared bats demonstrated significant increases in mean body mass index in 2014; this pattern likely reflected a relatively better body condition in bats that survived exposure to the WNS fungus. Most species demonstrated highest body mass and body mass index values in late swarming compared with other sampling periods, with tri-colored bats showing the greatest percent increase in body mass (42.5%) and body mass index (42.9%) prior to entering hibernation. These data indicate significant intraspecific variation in body condition of cave-hibernating bat species, both among years and between the seasons of autumn swarming and spring staging. We suggest this variation is likely to have implications for the relative vulnerability of species to WNS infection across the distribution of the Pseudogymnoascus fungus.

Genetic connectivity among swarming sites in the wide ranging and recently declining little brown bat (Myotis lucifugus)

Characterizing movement dynamics and spatial aspects of gene flow within a species permits inference on population structuring. As patterns of structuring are products of historical and current demographics and gene flow, assessment of structure through time can yield an understanding of evolutionary dynamics acting on populations that are necessary to inform management. Recent dramatic population declines in hibernating bats in eastern North America from white-nose syndrome have prompted the need for information on movement dynamics for multiple bat species. We characterized population genetic structure of the little brown bat, Myotis lucifugus, at swarming sites in southeastern Canada using 9 nuclear microsatellites and a 292-bp region of the mitochondrial genome. Analyses of F ST, ΦST, and Bayesian clustering (STRUCTURE) found weak levels of genetic structure among swarming sites for the nuclear and mitochondrial genome (Global F ST = 0.001, P < 0.05, Global ΦST = 0.045, P < 0.01, STRUCTURE K = 1) suggesting high contemporary gene flow. Hierarchical AMOVA also suggests little structuring at a regional (provincial) level. Metrics of nuclear genetic structure were not found to differ between males and females suggesting weak asymmetries in gene flow between the sexes. However, a greater degree of mitochondrial structuring does support male-biased dispersal long term. Demographic analyses were consistent with past population growth and suggest a population expansion occurred from approximately 1250 to 12,500 BP, following Pleistocene deglaciation in the region. Our study suggests high gene flow and thus a high degree of connectivity among bats that visit swarming sites whereby mainland areas of the region may be best considered as one large gene pool for management and conservation.