The roles of different gene expression regulators in acoustic variation in the intermediate horseshoe bat revealed by long-read and short-read RNA sequencing data

Gene expression changes contribute greatly to phenotypic variations in nature. Studying patterns of regulators of gene expression is important to fully understand the molecular mechanism underlying phenotypic variations. In horseshoe bats, the cochleae are finely tuned to echoes of call frequency. Here, using 2 recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis hainanus and R. a. himalayanus) with great acoustic variations as the system, we aim to explore relative roles of different regulators of gene expression (differential gene expression, alternative splicing (AS) and long non-coding RNAs (lncRNAs)) in phenotypic variation with a combination of Illumina short-read and Nanopore long-read RNA-seq data from the cochlea. Compared to R. a. hainanus, R. a. himalayanus exhibited much more upregulated differentially expressed genes (DEGs) and multiple of them may play important roles in the maintenance and damage repair of auditory hair cells. We identified 411 differentially expressed lncRNAs and their target DEGs upregulated in R. a. himalayanus were also mainly involved in a protective mechanism for auditory hair cells. Using 3 different methods of AS analysis, we identified several candidate alternatively spliced genes (ASGs) that expressed different isoforms which may be associated with acoustic divergence of the 2 subspecies. We observed significantly less overlap than expected between DEGs and ASGs, supporting complementary roles of differential gene expression and AS in generating phenotypic variations. Overall, our study highlights the importance of a combination of short-read and long-read RNA-seq data in examining the regulation of gene expression changes responsible for phenotypic variations.

Chromosome-level genome and population genomics of the intermediate horseshoe bat ( Rhinolophus affinis) reveal the molecular basis of virus tolerance in Rhinolophus and echolocation call frequency variation

Horseshoe bats (genus Rhinolophus, family Rhinolophidae) represent an important group within chiropteran phylogeny due to their distinctive traits, including constant high-frequency echolocation, rapid karyotype evolution, and unique immune system. Advances in evolutionary biology, supported by high-quality reference genomes and comprehensive whole-genome data, have significantly enhanced our understanding of species origins, speciation mechanisms, adaptive evolutionary processes, and phenotypic diversity. However, genomic research and understanding of the evolutionary patterns of Rhinolophus are severely constrained by limited data, with only a single published genome of R. ferrumequinum currently available. In this study, we constructed a high-quality chromosome-level reference genome for the intermediate horseshoe bat ( R. affinis). Comparative genomic analyses revealed potential genetic characteristics associated with virus tolerance in Rhinolophidae. Notably, we observed expansions in several immune-related gene families and identified various genes functionally associated with the SARS-CoV-2 signaling pathway, DNA repair, and apoptosis, which displayed signs of rapid evolution. In addition, we observed an expansion of the major histocompatibility complex class II (MHC-II) region and a higher copy number of the HLA- DQB2 gene in horseshoe bats compared to other chiropteran species. Based on whole-genome resequencing and population genomic analyses, we identified multiple candidate loci (e.g., GLI3) associated with variations in echolocation call frequency across R. affinis subspecies. This research not only expands our understanding of the genetic characteristics of the Rhinolophus genus but also establishes a valuable foundation for future research.

Introgression between highly divergent fungal sister species

To understand how species evolve and adapt to changing environments, it is important to study gene flow and introgression due to their influence on speciation and radiation events. Here, we apply a novel experimental system for investigating these mechanisms using natural populations. The system is based on two fungal sister species with morphological and ecological similarities occurring in overlapping habitats. We examined introgression between these species by conducting whole genome sequencing of individuals from populations in North America and Europe. We assessed genome-wide nucleotide divergence and performed crossing experiments to study reproductive barriers. We further used ABBA-BABA statistics together with a network analysis to investigate introgression, and conducted demographic modelling to gain insight into divergence times and introgression events. The results revealed that the species are highly divergent and incompatible in vitro. Despite this, small regions of introgression were scattered throughout the genomes and one introgression event likely involves a ghost population (extant or extinct). This study demonstrates that introgression can be found among divergent species and that population histories can be studied without collections of all the populations involved. Moreover, the experimental system is shown to be a useful tool for research on reproductive isolation in natural populations.

Characterization of microRNA and gene expression in the cochlea of an echolocating bat (Rhinolophus affinis)

MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression and play key roles in many biological processes, such as development and response to multiple stresses. However, little is known about their roles in generating novel phenotypes and phenotypic variation during the course of animal evolution. Here, we, for the first time, characterized the miRNAs of the cochlea in an echolocating bat (Rhinolophus affinis). We sampled eight individuals from two R. affinis subspecies with significant echolocation call frequency differences. We identified 365 miRNAs and 121 of them were novel. By searching sequences of these miRNAs precursors in multiple high-quality mammal genomes, we found one specific miRNA shared by all echolocating bats but not present in all other nonecholocating mammals. The targeted genes of this miRNA included several known hearing genes (e.g., KCNQ4 and GJB6). Together with the matched mRNA-seq data, we identified 1766 differentially expressed genes (DEGs) between the two subspecies and 555 of them were negatively regulated by differentially expressed miRNAs (DEMs). We found that almost half of known hearing genes in the list of all DEGs were regulated negatively by DEMs, suggesting an important role of miRNAs in call frequency variation of the two subspecies. These targeted DEGs included several important hearing genes (e.g., Piezo1, Piezo2, and CDH23) that have been shown to be important in ultrasonic hearing of echolocating mammals.

Next-generation phylogeography of the banded newts (Ommatotriton): A phylogenetic hypothesis for three ancient species with geographically restricted interspecific gene flow and deep intraspecific genetic structure

Technological developments now make it possible to employ many markers for many individuals in a phylogeographic setting, even for taxa with large and complex genomes such as salamanders. The banded newt (genus Ommatotriton) from the Near East has been proposed to contain three species (O. nesterovi, O. ophryticus and O. vittatus) with unclear phylogenetic relationships, apparently limited interspecific gene flow and deep intraspecific geographic mtDNA structure. We use parallel tagged amplicon sequencing to obtain 177 nuclear DNA markers for 35 banded newts sampled throughout the range. We determine population structure (with Bayesian clustering and principal component analysis), interspecific gene flow (by determining the distribution of species-diagnostic alleles) and phylogenetic relationships (by maximum likelihood inference of concatenated sequence data and based on a summary-coalescent approach). We confirm that the three proposed species are genetically distinct. A sister relationship between O. nesterovi and O. ophryticus is suggested. We find evidence for introgression between O. nesterovi and O. ophryticus, but this is geographically limited. Intraspecific structuring is extensive, with the only recognized banded newt subspecies, O. vittatus cilicensis, representing the most distinct lineage below the species level. While mtDNA mostly mirrors the pattern observed in nuclear DNA, all banded newt species show mito-nuclear discordance as well.

Mitonuclear mismatch alters nuclear gene expression in naturally introgressed Rhinolophus bats

Background

Mitochondrial function involves the interplay between mitochondrial and nuclear genomes. Such mitonuclear interactions can be disrupted by the introgression of mitochondrial DNA between taxa or divergent populations. Previous studies of several model systems (e.g. Drosophila) indicate that the disruption of mitonuclear interactions, termed mitonuclear mismatch, can alter nuclear gene expression, yet few studies have focused on natural populations.

Results

Here we study a naturally introgressed population in the secondary contact zone of two subspecies of the intermediate horseshoe bat (Rhinolophus affinis), in which individuals possess either mitonuclear matched or mismatched genotypes. We generated transcriptome data for six tissue types from five mitonuclear matched and five mismatched individuals. Our results revealed strong tissue-specific effects of mitonuclear mismatch on nuclear gene expression with the largest effect seen in pectoral muscle. Moreover, consistent with the hypothesis that genes associated with the response to oxidative stress may be upregulated in mitonuclear mismatched individuals, we identified several such gene candidates, including DNASE1L3, GPx3 and HSPB6 in muscle, and ISG15 and IFI6 in heart.

Conclusion

Our study reveals how mitonuclear mismatch arising from introgression in natural populations is likely to have fitness consequences. Underlying the processes that maintain mitonuclear discordance is a step forward to understand the role of mitonuclear interactions in population divergence and speciation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00424-x.

Extensive alternative splicing triggered by mitonuclear mismatch in naturally introgressed Rhinolophus bats

Mitochondrial function needs strong interactions of mitochondrial and nuclear (mitonuclear) genomes, which can be disrupted by mitonuclear mismatch due to mitochondrial DNA (mtDNA) introgression between two formerly isolated populations or taxa. This mitonuclear disruption may cause severe cellular stress in mismatched individuals. Gene expression changes and alternative splicing (AS) are two important transcriptional regulations to respond to environmental or cellular stresses. We previously identified a naturally introgressed population in the intermediate horseshoe bat (Rhinolophus affinis). Individuals from this population belong to R. a. himalayanus and share almost identical nuclear genetic background; however, some of them had mtDNA from another subspecies (R. a. macrurus). With this unique natural system, we examined gene expression changes in six tissues between five mitonuclear mismatched and five matched individuals. A small number of differentially expressed genes (DEGs) were identified, and functional enrichment analysis revealed that most DEGs were related to immune response although some may be involved in response to oxidative stress. In contrast, we identified extensive AS events and alternatively spliced genes (ASGs) between mismatched and matched individuals. Functional enrichment analysis revealed that multiple ASGs were directly or indirectly associated with energy production in mitochondria which is vital for survival of organism. To our knowledge, this is the first study to examine the role of AS in responding to cellular stress caused by mitonuclear mismatch in natural populations. Our results suggest that AS may play a more important role than gene expression regulation in responding to severe environmental or cellular stresses.

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.

The complete mitochondrial genome of Rhinolophus affinis himalayanus

Here we generated the complete mitochondrial genome of one subspecies of R. affinis (R. affinis himalayanus) using next generation sequencing and Sanger sequencing. The length of the complete mitochondrial genome was 16,886 bp, containing 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and a non-coding control region. A maximum-likelihood tree based on the 13 concatenated mitochondrial protein-coding genes of 16 Rhinolophus taxon and one outgroup Hipposideros armiger indicates that R. affinis shows a closer relationship with R. sinicus complex than with other taxa.

Genetic Divergence Across Glacial Refugia Despite Interglacial Gene Flow in a Crested Newt

MtDNA-based phylogeography has illuminated the impact of the Pleistocene Ice Age on species distribution dynamics and the build-up of genetic divergence. The well-known shortcomings of mtDNA in biogeographical inference can be compensated by integrating multilocus data and species distribution modelling into phylogeography. We re-visit the phylogeography of the Italian crested newt (Triturus carnifex), a species distributed in two of Europe’s main glacial refugia, the Balkan and Italian Peninsulas. While a new 51 nuclear DNA marker dataset supports the existence of three lineages previously suggested by mtDNA (Balkan, northern Italy and southern Italy), the nuclear DNA dataset also provides improved resolution where these lineages have obtained secondary contact. We observe geographically restricted admixture at the contact between the Balkan and northern Italy gene pools and identify a potential mtDNA ghost lineage here. At the contact between the northern and southern Italy gene pools we find admixture over a broader area, as well as asymmetric mtDNA introgression. Our species distribution model is in agreement with a distribution restricted to distinct refugia during Pleistocene glacial cycles and postglacial expansion with secondary contact. Our study supports: (1) the relevance of the north-western Balkan Peninsula as a discrete glacial refugium; (2) the importance of north-eastern Italy and the northern Apennine as suture zones; and (3) the applicability of a refugia-within-refugia scenario within the Italian Peninsula.

Transcriptome sequencing of cochleae from constant-frequency and frequency-modulated echolocating bats

Echolocating bats are fascinating for their ability to ‘see’ the world in the darkness. Ultrahigh frequency hearing is essential for echolocation. In this study we collected cochlear tissues from constant-frequency (CF) bats (two subspecies of Rhinolophus affinis, Rhinolophidae) and frequency-modulated (FM) bats (Myotis ricketti, Vespertilionidae) and applied PacBio single-molecule real-time isoform sequencing (Iso-seq) technology to generate the full-length (FL) transcriptomes for the three taxa. In total of 10103, 9676 and 10504 non-redundant FL transcripts for R. a. hainanus, R. a. himalayanus and Myotis ricketti were obtained respectively. These data present a comprehensive list of transcripts involved in ultrahigh frequency hearing of echolocating bats including 26342 FL transcripts, 24833 of which are annotated by public databases. No further comparative analyses were performed on the current data in this study. This data can be reused to quantify gene or transcript expression, assess the level of alternative splicing, identify novel transcripts and improve genome annotation of bat species.

Measurement(s)	cochlea • transcriptome • sequence feature annotation
Technology Type(s)	isoform sequencing • sequence annotation
Factor Type(s)	species
Sample Characteristic - Organism	Rhinolophus affinis • Myotis ricketti

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12999614

Echolocation call frequency variation in horseshoe bats: molecular basis revealed by comparative transcriptomics

Recently diverged taxa with contrasting phenotypes offer opportunities for unravelling the genetic basis of phenotypic variation in nature. Horseshoe bats are a speciose group that exhibit a derived form of high-duty cycle echolocation in which the inner ear is finely tuned to echoes of the narrowband call frequency. Here, by focusing on three recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis) that display divergent echolocation call frequencies, we aim to identify candidate loci putatively involved in hearing frequency variation. We used de novo transcriptome sequencing of two mainland taxa (himalayanus and macrurus) and one island taxon (hainanus) to compare expression profiles of thousands of genes. By comparing taxa with divergent call frequencies (around 15 kHz difference), we identified 252 differentially expressed genes, of which six have been shown to be involved in hearing or deafness in human/mouse. To obtain further validation of these results, we applied quantitative reverse transcription-PCR to the candidate gene FBXL15 and found a broad association between the level of expression and call frequency across taxa. The genes identified here represent strong candidate loci associated with hearing frequency variation in bats.

Genome-wide data reveal discordant mitonuclear introgression in the intermediate horseshoe bat (Rhinolophus affinis)

Closely related taxa often exhibit mitonuclear discordance attributed to introgression of mitochondrial DNA (mtDNA), yet few studies have considered the underlying causes of mtDNA introgression. Here we test for demographic versus adaptive processes as explanations for mtDNA introgression in three subspecies of the intermediate horseshoe bat (Rhinolophus affinis). We generated sequences of 1692 nuclear genes and 13 mitochondrial protein-coding genes for 48 individuals. Phylogenetic reconstructions based on 320 exon sequences and 2217 single nucleotide polymorphisms (SNPs) both revealed conflicts between the species tree and mtDNA tree. These results, together with geographic patterns of mitonuclear discordance, and shared identical or near-identical mtDNA sequences, suggest extensive introgression of mtDNA between the two parapatric mainland subspecies. Under demographic hypotheses, we would also expect to uncover traces of ncDNA introgression, however, population structure and gene flow analyses revealed little nuclear admixture. Furthermore, we found inconsistent estimates of the timing of population expansion and that of the most recent common ancestor for the clade containing introgressed haplotypes. Without a clear demographic explanation, we also examined whether introgression likely arises from adaptation. We found that five mtDNA genes contained fixed amino acid differences between introgressed and non-introgressed individuals, including putative positive selection found in one codon, although this did not show introgression. While our evidence for rejecting demographic hypotheses is arguably stronger than that for rejecting adaptation, we find no definitive support for either explanation. Future efforts will focus on larger-scale resequencing to decipher the underlying causes of discordant mitonuclear introgression in this system.

The complex evolutionary history of big-eared horseshoe bats (Rhinolophus macrotis complex): insights from genetic, morphological and acoustic data

Palaeoclimatic oscillations and different landscapes frequently result in complex population-level structure or the evolution of cryptic species. Elucidating the potential mechanisms is vital to understanding speciation events. However, such complex evolutionary patterns have rarely been reported in bats. In China, the Rhinolophus macrotis complex contains a large form and a small form, suggesting the existence of a cryptic bat species. Our field surveys found these two sibling species have a continuous and widespread distribution with partial sympatry. However, their evolutionary history has received little attention. Here, we used extensive sampling, morphological and acoustic data, as well as different genetic markers to investigate their evolutionary history. Genetic analyses revealed discordance between the mitochondrial and nuclear data. Mitochondrial data identified three reciprocally monophyletic lineages: one representing all small forms from Southwest China, and the other two containing all large forms from Central and Southeast China, respectively. The large form showed paraphyly with respect to the small form. However, clustering analyses of microsatellite and Chd1 gene sequences support two divergent clusters separating the large form and the small form. Moreover, morphological and acoustic analyses were consistent with nuclear data. This unusual pattern in the R. macrotis complex might be accounted for by palaeoclimatic oscillations, shared ancestral polymorphism and/or interspecific hybridization.

Historical and current introgression in a Mesoamerican hummingbird species complex: a biogeographic perspective

The influence of geologic and Pleistocene glacial cycles might result in morphological and genetic complex scenarios in the biota of the Mesoamerican region. We tested whether berylline, blue-tailed and steely-blue hummingbirds, Amazilia beryllina, Amazilia cyanura and Amazilia saucerottei, show evidence of historical or current introgression as their plumage colour variation might suggest. We also analysed the role of past and present climatic events in promoting genetic introgression and species diversification. We collected mitochondrial DNA (mtDNA) sequence data and microsatellite loci scores for populations throughout the range of the three Amazilia species, as well as morphological and ecological data. Haplotype network, Bayesian phylogenetic and divergence time inference, historical demography, palaeodistribution modelling, and niche divergence tests were used to reconstruct the evolutionary history of this Amazilia species complex. An isolation-with-migration coalescent model and Bayesian assignment analysis were assessed to determine historical introgression and current genetic admixture. mtDNA haplotypes were geographically unstructured, with haplotypes from disparate areas interdispersed on a shallow tree and an unresolved haplotype network. Assignment analysis of the nuclear genome (nuDNA) supported three genetic groups with signs of genetic admixture, corresponding to: (1) A. beryllina populations located west of the Isthmus of Tehuantepec; (2) A. cyanura populations between the Isthmus of Tehuantepec and the Nicaraguan Depression (Nuclear Central America); and (3) A. saucerottei populations southeast of the Nicaraguan Depression. Gene flow and divergence time estimates, and demographic and palaeodistribution patterns suggest an evolutionary history of introgression mediated by Quaternary climatic fluctuations. High levels of gene flow were indicated by mtDNA and asymmetrical isolation-with-migration, whereas the microsatellite analyses found evidence for three genetic clusters with distributions corresponding to isolation by the Isthmus of Tehuantepec and the Nicaraguan Depression and signs of admixture. Historical levels of migration between genetically distinct groups estimated using microsatellites were higher than contemporary levels of migration. These results support the scenario of secondary contact and range contact during the glacial periods of the Pleistocene and strongly imply that the high levels of structure currently observed are a consequence of the limited dispersal of these hummingbirds across the isthmus and depression barriers.

Differential introgression among loci across a hybrid zone of the intermediate horseshoe bat (Rhinolophus affinis)

Background

Hybrid zones formed by the secondary contact of divergent lineages represent natural laboratories for studying the genetic basis of speciation. Here we tested for patterns of differential introgression among three X-linked and 11 autosomal regions to identify candidate loci related to either reproductive isolation or adaptive introgression across a hybrid zone between two Chinese mainland subspecies of the intermediate horseshoe bat Rhinolophus affinis: R. a. himalayanus and R. a. macrurus.

Results

Our results support the previous suggestion that macrurus formed when a third subspecies (R. a. hainanus) recolonized the mainland from Hainan Island, and that himalayanus is the ancestral taxon. However, this overall evolutionary history was not reflected in all loci examined, with considerable locus-wise heterogeneity seen in gene tree topologies, levels of polymorphism, genetic differentiation and rates of introgression. Coalescent simulations suggested levels of lineage mixing seen at some nuclear loci might result from incomplete lineage sorting. Isolation with migration models supported evidence of gene flow across the hybrid zone at one intronic marker of the hearing gene Prestin.

Conclusions

We suggest that phylogenetic discordance with respect to the species tree seen here is likely to arise via a combination of incomplete lineage sorting and a low incidence of introgression although we cannot rule out other explanations such as selection and recombination. Two X-linked loci and one autosomal locus were identified as candidate regions related to reproductive isolation across the hybrid zone. Our work highlights the importance of including multiple genomic regions in characterizing patterns of divergence and gene flow across a hybrid zone.

Sensory trait variation in an echolocating bat suggests roles for both selection and plasticity

Background

Across heterogeneous environments selection and gene flow interact to influence the rate and extent of adaptive trait evolution. This complex relationship is further influenced by the rarely considered role of phenotypic plasticity in the evolution of adaptive population variation. Plasticity can be adaptive if it promotes colonization and survival in novel environments and in doing so may increase the potential for future population differentiation via selection. Gene flow between selectively divergent environments may favour the evolution of phenotypic plasticity or conversely, plasticity itself may promote gene flow, leading to a pattern of trait differentiation in the presence of gene flow. Variation in sensory traits is particularly informative in testing the role of environment in trait and population differentiation. Here we test the hypothesis of ‘adaptive differentiation with minimal gene flow’ in resting echolocation frequencies (RF) of Cape horseshoe bats (Rhinolophus capensis) across a gradient of increasingly cluttered habitats.

Results

Our analysis reveals a geographically structured pattern of increasing RF from open to highly cluttered habitats in R. capensis; however genetic drift appears to be a minor player in the processes influencing this pattern. Although Bayesian analysis of population structure uncovered a number of spatially defined mitochondrial groups and coalescent methods revealed regional-scale gene flow, phylogenetic analysis of mitochondrial sequences did not correlate with RF differentiation. Instead, habitat discontinuities between biomes, and not genetic and geographic distances, best explained echolocation variation in this species. We argue that both selection for increased detection distance in relatively less cluttered habitats and adaptive phenotypic plasticity may have influenced the evolution of matched echolocation frequencies and habitats across different populations.

Conclusions

Our study reveals significant sensory trait differentiation in the presence of historical gene flow and suggests roles for both selection and plasticity in the evolution of echolocation variation in R. capensis. These results highlight the importance of population level analyses to i) illuminate the subtle interplay between selection, plasticity and gene flow in the evolution of adaptive traits and ii) demonstrate that evolutionary processes may act simultaneously and that their relative influence may vary across different environments.

Phylogeography and postglacial recolonization of Europe by Rhinolophus hipposideros: evidence from multiple genetic markers

The demographic history of Rhinolophus hipposideros (lesser horseshoe bat) was reconstructed across its European, North African and Middle-Eastern distribution prior to, during and following the most recent glaciations by generating and analysing a multimarker data set. This data set consisted of an X-linked nuclear intron (Bgn; 543 bp), mitochondrial DNA (cytb-tRNA-control region; 1630 bp) and eight variable microsatellite loci for up to 373 individuals from 86 localities. Using this data set of diverse markers, it was possible to determine the species' demography at three temporal stages. Nuclear intron data revealed early colonization into Europe from the east, which pre-dates the Quaternary glaciations. The mtDNA data supported multiple glacial refugia across the Mediterranean, the largest of which were found in the Ibero-Maghreb region and an eastern location (Anatolia/Middle East)-that were used by R. hipposideros during the most recent glacial cycles. Finally, microsatellites provided the most recent information on these species' movements since the Last Glacial Maximum and suggested that lineages that had diverged into glacial refugia, such as in the Ibero-Maghreb region, have remained isolated. These findings should be used to inform future conservation management strategies for R. hipposideros and show the power of using a multimarker data set for phylogeographic studies.