Instances of hybridization between mammalian taxa in the wild are rarely documented. To test for introgression between sibling species of horseshoe bat (Rhinolophus yunanensis and R. pearsoni) and two subspecies of the latter (R. p. pearsoni and R. p. chinensis), we sequenced two mtDNA and two ncDNA markers in individuals sampled from multiple localities within their overlapping ranges. The interspecific mtDNA gene tree corresponded to the expected taxonomic divisions, and coalescent-based analyses suggested divergence occurred around 4 MYA. However, these relationships strongly conflicted with those recovered from two independent nuclear gene trees, in which R. yunanensis clustered with R. p. pearsoni to the exclusion of R. p. chinensis. This geographically widespread discordance is best explained by large-scale historical introgression of ncDNA from R. yunanensis to R. pearsoni by male-mediated exchange in mixed species colonies during Pleistocene glacial periods, when ranges may have contracted and overlapped more than at present. Further species tree-gene tree conflicts were detected between R. p. pearsoni and R. p. chinensis, also indicating past and/or current introgression in their overlapping regions. However, here the patterns point to asymmetric mtDNA introgression without ncDNA introgression. Analyses of coalescence times indicate this exchange has occurred subsequent to the divergence of these subspecies from their common ancestor. Our work highlights the importance of using multiple data sets for reconstructing phylogeographic histories and resolving taxonomic relationships.

Rhinolophus (Rhinolophidae) is the second most speciose genus in Chiroptera and has extensively diversified diploid chromosome numbers (from 2n = 28 to 62). In spite of many attempts to explore the karyotypic evolution of this genus, most studies have been based on conventional Giemsa staining rather than G-banding. Here we have made a whole set of chromosome-specific painting probes from flow-sorted chromosomes of Aselliscus stoliczkanus (Hipposideridae). These probes have been utilized to establish the first genome-wide homology maps among six Rhinolophus species with four different diploid chromosome numbers (2n = 36, 44, 58, and 62) and three species from other families: Rousettus leschenaulti (2n = 36, Pteropodidae), Hipposideros larvatus (2n = 32, Hipposideridae), and Myotis altarium (2n = 44, Vespertilionidae) by fluorescence in situ hybridization. To facilitate integration with published maps, human paints were also hybridized to A. stoliczkanus chromosomes. Our painting results substantiate the wide occurrence of whole-chromosome arm conservation in Rhinolophus bats and suggest that Robertsonian translocations of different combinations account for their karyotype differences. Parsimony analysis using chromosomal characters has provided some new insights into the Rhinolophus ancestral karyotype and phylogenetic relationships among these Rhinolophus species so far studied. In addition to Robertsonian translocations, our results suggest that whole-arm (reciprocal) translocations involving multiple non-homologous chromosomes as well could have been involved in the karyotypic evolution within Rhinolophus, in particular those bats with low and medium diploid numbers.

The Vespertilionidae is the largest family in the order Chiroptera and has a worldwide distribution in the temperate and tropical regions. In order to further clarify the karyotype relationships at the lower taxonomic level in Vespertilionidae, genome-wide comparative maps have been constructed between Myotis myotis (MMY, 2n = 44) and six vesper bats from China: Myotis altarium (MAL, 2n = 44), Hypsugo pulveratus (HPU, 2n = 44), Nyctalus velutinus (NVE, 2n = 36), Tylonycteris robustula (TRO, 2n = 32), Tylonycteris sp. (TSP, 2n = 30)and Miniopterus fuliginosus (MFU, 2n = 46) by cross-species chromosome painting with a set of painting probes derived from flow-sorted chromosomes of Myotis myotis. Each Myotis myotis autosomal probe detected a single homologous chromosomal segment in the genomes of these six vesper bats except for MMY chromosome 3/4 paint which hybridized onto two chromosomes in the genome of M. fuliginosus. Our results show that Robertsonian translocation is the main mode of karyotype evolution in Vespertilionidae and that the addition of heterochromatic material also plays an important role in the karyotypic evolution of the genera Tylonycteris and Nyctalus. Two conserved syntenic associations (MMY9 + 23 and 18 + 19) could be the synapomorphic features for the genus Tylonycteris. The integration of our maps with the published maps has enabled us to deduce chromosomal homologies between human and these six vesper bats and provided new insight into the karyotype evolution of the family Vespertilionidae.

Although the monophyly of Chiroptera is well supported by many independent studies, higher-level systematics, e.g. the monophyly of microbats, remains disputed by morphological and molecular studies. Chromosomal rearrangements, as one type of rare genomic changes, have become increasingly popular in phylogenetic studies as alternatives to molecular and other morphological characters. Here, the representatives of families Megadermatidae and Emballonuridae are studied by comparative chromosome painting for the first time. The results have been integrated into published comparative maps, providing an opportunity to assess genome-wide chromosomal homologies between the representatives of eight bat families. Our results further substantiate the wide occurrence of Robertsonian translocations in bats, with the possible involvement of whole-arm reciprocal translocations (WARTs). In order to search for valid cytogenetic signature(s) for each family and superfamily, evolutionary chromosomal rearrangements identified by chromosomal painting and/or banding comparison are subjected to two independent analyses: (1) a cladistic analysis using parsimony and (2) the mapping of these chromosomal changes onto the molecularly defined phylogenetic tree available from the literature. Both analyses clearly indicate the prevalence of homoplasic events that reduce the reliability of chromosomal characters for resolving interfamily relationships in bats.

Bats are a unique but enigmatic group of mammals and have a world-wide distribution. The phylogenetic relationships of extant bats are far from being resolved. Here, we investigated the karyotypic relationships of representative species from four families of the order Chiroptera by comparative chromosome painting and banding. A complete set of painting probes derived from flow-sorted chromosomes of Myotis myotis (family Vespertilionidae) were hybridized onto metaphases of Cynopterus sphinx (2n = 34, family Pteropodidae), Rhinolophus sinicus (2n=36, family Rhinolophidae) and Aselliscus stoliczkanus (2n=30, family Hipposideridae) and delimited 27, 30 and 25 conserved chromosomal segments in the three genomes, respectively. The results substantiate that Robertsonian translocation is the main mode of chromosome evolution in the order Chiroptera, with extensive conservation of whole chromosomal arms. The use of M. myotis (2n=44) probes has enabled the integration of C. sphinx, R. sinicus and A. stoliczkanus chromosomes into the previously established comparative maps between human and Eonycteris spelaea (2n=36), Rhinolophus mehelyi (2n=58), Hipposideros larvatus (2n=32), and M. myotis. Our results provide the first cytogenetic signature rearrangement that supports the grouping of Pteropodidae and Rhinolophoidea in a common clade (i.e. Pteropodiformes or Yinpterochiroptera) and thus improve our understanding on the karyotypic relationships and genome phylogeny of these bat species.

Closely related taxa living in sympatry provide good opportunities to investigate the origin of barriers to gene flow as well as the extent of reproductive isolation. The only two recognized subspecies of the Chinese rufous horseshoe bat Rhinolophus sinicus are characterized by unusual relative distributions in which R. s. septentrionalis is restricted to a small area within the much wider range of its sister taxon R. s. sinicus. To determine the history of lineage divergence and gene flow between these taxa, we applied phylogenetic, demographic and coalescent analyses to multi-locus datasets. MtDNA gene genealogies and microsatellite-based clustering together revealed three divergent lineages of sinicus, corresponding to Central China, East China and the offshore Hainan Island. However, the central lineage of sinicus showed a closer relationship with septentrionalis than with other lineages of R. s. sinicus, in contrary to morphological data. Paraphyly of sinicus could result from either past asymmetric mtDNA introgression between these two taxa, or could suggest septentrionalis evolved in situ from its more widespread sister subspecies. To test between these hypotheses, we applied coalescent-based phylogenetic reconstruction and Approximate Bayesian Computation (ABC). We found that septentrionalis is likely to be the ancestral taxon and therefore a recent origin of this subspecies can be ruled out. On the other hand, we found a clear signature of asymmetric mtDNA gene flow from septentrionalis into central populations of sinicus yet no nuclear gene flow, thus strongly pointing to historical mtDNA introgression. We suggest that the observed deeply divergent lineages within R. sinicus probably evolved in isolation in separate Pleistocene refugia, although their close phylogeographic correspondence with distinct eco-environmental zones suggests that divergent selection might also have promoted broad patterns of population genetic structure.

Phylogenetic conflicts between genetic markers can help to disentangle complex histories of phylogeography and introgression among taxa. We previously proposed that the Chinese mainland subspecies of the intermediate horseshoe bat Rhinolophus affinis himalayanus colonized Hainan Island to form the subspecies R. a. hainanus. Subsequent recolonization of the mainland formed a third taxon, R. a macrurus, and a secondary contact zone with the ancestral himalayanus. To test for historical and recurrent genetic exchange between these mainland subspecies, we sampled populations of each from two parapatric zones and undertook analyses using one mtDNA marker, three nuclear genes and 14 microsatellites. Nuclear DNA, echolocation call frequencies and morphological data all recovered two taxa; however, a mtDNA phylogeny revealed two himalayanus clades, of which one clustered with macrurus, as well as some shared or related mtDNA haplotypes in eastern populations. Isolation-with-migration (IM) models suggested some mtDNA gene flow from macrurus to himalayanus. However, strong population structure in himalayanus raises the possibility that macrurus captured mtDNA from a coastal population of himalayanus that has since become rare or extinct. To reconcile these two sets of results, we suggest that the IM estimates might reflect historical mtDNA gene flow among populations of himalayanus, before mtDNA was subsequently captured by macrurus. Finally, microsatellite-based ABC analyses supported the island origin of macrurus; however, mtDNA-based ABC analyses suggest this taxon might have evolved on the mainland. Our findings highlight the importance of understanding population history and structure for interpreting hybridization and introgression events.

Museums hold most of the world's most valuable biological specimens and tissues collected, including type material that is often decades or even centuries old. Unfortunately, traditional museum collection and storage methods were not designed to preserve the nucleic acids held within the material, often reducing its potential viability and value for many genetic applications. High‐throughput sequencing technologies and associated applications offer new opportunities for obtaining sequence data from museum samples. In particular, target sequence capture offers a promising approach for recovering large numbers of orthologous loci from relatively small amounts of starting material. In the present study, we test the utility of target sequence capture for obtaining data from museum‐held material from a speciose mammalian genus: the horseshoe bats (Rhinolophidae: Chiroptera). We designed a ‘bait’ for capturing > 3600 genes and applied this to 10 species of horseshoe bat that had been collected between 93 and 7 years ago and preserved using a range of methods. We found that the mean recovery rate per species was approximately 89% of target genes with partial sequence coverage, ranging from 3024 to 3186 genes recovered. On average, we recovered 1206 genes with ≥ 90% sequence coverage, per species. Our findings provide good support for the application of large‐scale bait capture across congeneric species spanning approximately 15 Myr of evolution. On the other hand, we observed no clear association between the success of capture and the phylogenetic distance from the bait model, although sample sizes precluded a formal test.