Recurrent replacement of mtDNA and cryptic hybridization between two sibling bat species Myotis myotis and Myotis blythii

Distribution of bat species in Western Asia: Occurrence records from the Western Asia Bat Research Network (WAB-Net) project

Background

Western Asia represents a mixing pot of diverse bat species with distributions spanning across other geographic regions. Yet, relative to other regions, there is a significant gap in coordinated bat research in Western Asia, thereby resulting in a relatively limited number of curated occurrence records.

New information

The Western Asia Bat Research Network (WAB-Net) project was created to strengthen research capacity and knowledge of the diversity and distribution of bat species in a little-studied region, as well as to collect data to characterise the diversity and prevalence of coronaviruses in diverse bat species. Over a four-year period (2018-2022), we documented 4,278 individual records for 41 bat species using a cross-sectional survey approach at 50 sites in seven Western Asian countries, specifically Armenia, Azerbaijan, Georgia, Jordan, Oman, Pakistan and Turkiye. At each site, we captured, on average, 90 individual bats (range: 9-131) over multiple consecutive nights and used standardised methods to collect demographic and morphological data from captured individuals. We additionally completed a systematic evaluation of environmental characterisation and human-bat interactions at all 50 sites. Here, we report individual occurrence records and site conditions from this multi-country, multi-year sampling effort.

Cryptic hybridization between the ancient lineages of Natterer's bat (Myotis nattereri)

Studying hybrid zones that form between morphologically cryptic taxa offers valuable insights into the mechanisms of cryptic speciation and the evolution of reproductive barriers. Although hybrid zones have long been the focus of evolutionary studies, the awareness of cryptic hybrid zones increased recently due to rapidly growing evidence of biological diversity lacking obvious phenotypic differentiation. The characterization of cryptic hybrid zones with genome-wide analysis is in its early stages and offers new perspectives for studying population admixture and thus the impact of gene flow. In this study, we investigate the population genomics of the Myotis nattereri complex in one of its secondary contact zones, where a putative hybrid zone is formed between two of its cryptic lineages. By utilizing a whole-genome shotgun sequencing approach, we aim to characterize this cryptic hybrid zone in detail. Demographic analysis suggests that the cryptic lineages diverged during the Pliocene, c. 3.6 million years ago. Despite this ancient separation, the populations in the contact zone exhibit mitochondrial introgression and a considerable amount of mixing in nuclear genomes. The genomic structure of the populations corresponds to geographic locations and the genomic admixture changes along a geographic gradient. These findings suggest that there is no effective hybridization barrier between both lineages, nevertheless, their population structure is shaped by dispersal barriers. Our findings highlight how such deeply diverged cryptic lineages can still readily hybridize in secondary contact.

Rampant chloroplast capture in Sarracenia revealed by plastome phylogeny

Introgression can produce novel genetic variation in organisms that hybridize. Sympatric species pairs in the carnivorous plant genus Sarracenia L. frequently hybridize, and all known hybrids are fertile. Despite being a desirable system for studying the evolutionary consequences of hybridization, the extent to which introgression occurs in the genus is limited to a few species in only two field sites. Previous phylogenomic analysis of Sarracenia estimated a highly resolved species tree from 199 nuclear genes, but revealed a plastid genome that is highly discordant with the species tree. Such cytonuclear discordance could be caused by chloroplast introgression (i.e. chloroplast capture) or incomplete lineage sorting (ILS). To better understand the extent to which introgression is occurring in Sarracenia, the chloroplast capture and ILS hypotheses were formally evaluated. Plastomes were assembled de-novo from sequencing reads generated from 17 individuals in addition to reads obtained from the previous study. Assemblies of 14 whole plastomes were generated and annotated, and the remaining fragmented assemblies were scaffolded to these whole-plastome assemblies. Coding sequence from 79 homologous genes were aligned and concatenated for maximum-likelihood phylogeny estimation. The plastome tree is extremely discordant with the published species tree. Plastome trees were simulated under the coalescent and tree distance from the species tree was calculated to generate a null distribution of discordance that is expected under ILS alone. A t-test rejected the null hypothesis that ILS could cause the level of discordance seen in the plastome tree, suggesting that chloroplast capture must be invoked to explain the discordance. Due to the extreme level of discordance in the plastome tree, it is likely that chloroplast capture has been common in the evolutionary history of Sarracenia.

Inferring genetic structure when there is little: population genetics versus genomics of the threatened bat Miniopterus schreibersii across Europe

Despite their paramount importance in molecular ecology and conservation, genetic diversity and structure remain challenging to quantify with traditional genotyping methods. Next-generation sequencing holds great promises, but this has not been properly tested in highly mobile species. In this article, we compared microsatellite and RAD-sequencing (RAD-seq) analyses to investigate population structure in the declining bent-winged bat (Miniopterus schreibersii) across Europe. Both markers retrieved general patterns of weak range-wide differentiation, little sex-biased dispersal, and strong isolation by distance that associated with significant genetic structure between the three Mediterranean Peninsulas, which could have acted as glacial refugia. Microsatellites proved uninformative in individual-based analyses, but the resolution offered by genomic SNPs illuminated on regional substructures within several countries, with colonies sharing migrators of distinct ancestry without admixture. This finding is consistent with a marked philopatry and spatial partitioning between mating and rearing grounds in the species, which was suspected from marked-recaptured data. Our study advocates that genomic data are necessary to properly unveil the genetic footprints left by biogeographic processes and social organization in long-distant flyers, which are otherwise rapidly blurred by their high levels of gene flow.

Whole-genome Analysis Reveals Contrasting Relationships Among Nuclear and Mitochondrial Genomes Between Three Sympatric Bat Species

Understanding mechanisms involved in speciation can be challenging, especially when hybridization or introgression blurs species boundaries. In bats, resolving relationships of some closely related groups has proved difficult due subtle interspecific variation both in morphometrics and molecular data sets. The endemic South American Histiotus bats, currently considered a subgenus of Eptesicus, harbor unresolved phylogenetic relationships and of those is a trio consisting of two closely related species: Eptesicus (Histiotus) macrotus and Eptesicus (Histiotus) montanus, and their relationship with a third, Eptesicus (Histiotus) magellanicus. The three sympatric species bear marked resemblance to each other, but can be differentiated morphologically. Furthermore, previous studies have been unable to differentiate the species from each other at a molecular level. In order to disentangle the phylogenetic relationships of these species, we examined the differentiation patterns and evolutionary history of the three Eptesicus (H.) species at the whole-genome level. The nuclear DNA statistics between the species suggest strong gene flow and recent hybridization between E. (H.) montanus and E. (H.) macrotus, whereas E. (H.) magellanicus shows a higher degree of isolation. In contrast, mitochondrial DNA shows a closer relationship between E. (H.) magellanicus and E. (H.) montanus. Opposing patterns in mtDNA and nuclear markers are often due to differences in dispersal, and here it could be both as a result of isolation in refugia during the last glacial maximum and female philopatry and male-biased dispersal. In conclusion, this study shows the importance of both the nuclear and mitochondrial DNA in resolving phylogenetic relationships and species histories.

Species delineation and genetic structure of two Chaerephon species (C. pusillus and C. leucogaster) on Madagascar and the Comoro archipelago

Cryptic species diversity is known to be common in bats but remains challenging to study in these mammals, whose natural history traits render their sampling and monitoring challenging. For these animals, indirect genetic approaches provide a powerful tool to gain insight into the evolutionary history and ecology of cryptic bat species. The speciation history of the polyphyletic Chaerephon pumilus species group (Molossidae) is poorly understood, including those found on western Indian Ocean islands. Two species in this complex have been identified in the Comoros: C. pusillus and C. leucogaster. Here, we aim to genetically characterize these two species and investigate their spatial population genetic structure. Analyzing five nuclear microsatellite markers from 200 individuals and one mitochondrial DNA gene (Cyt-b) from 161 (out of the 200) individuals sampled on Madagascar and the Comoros, our findings indicated that these species are genetically differentiated. We observed mitonuclear discordance in numerous individuals (33% of the 161 mtDNA-sequenced individuals). Based on ABC analyses, we found that this pattern could potentially be the result of asymmetric introgressive hybridization from C. leucogaster to C. pusillus and calls for further studies on the demographic history of these species. Moreover, at the intra-specific level, analyses of the microsatellite loci suggested the evidence of a more pronounced, although weak, geographically based genetic structure in C. pusillus than in C. leucogaster. Altogether, our findings provide preliminary insights into the eco-evolutionary aspects of this species complex and warrant further research to understand hybridization dynamics and mechanisms responsible for mitonuclear discordance.

Ecological speciation by sympatric host shifts in a clade of herbivorous sea slugs, with introgression and localized mitochondrial capture between species

Host shifting in insect-plant systems was historically important to the development of ecological speciation theory, yet surprisingly few studies have examined whether host shifting drives diversification of marine herbivores. When small-bodied consumers feed and also mate on a preferred host, disruptive selection can split a population into host races despite gene flow. Support for host shifts is notably lacking for invertebrates associated with macroalgae, where the scale of dispersal by planktonic larvae often far exceeds the grain of host patchiness, and adults are typically less specialized than terrestrial herbivores. Here, we present a candidate example of ecological speciation in a clade of sea slugs that primarily consume green algae in the genus Caulerpa, including highly invasive species. Ancestral character state reconstructions supported 'sea grapes' (C. racemosa, C. lentillifera) as the ancestral host for a tropical radiation of 12 Elysia spp., with one shift onto alternative Caulerpa spp. in the Indo-Pacific. A Caribbean radiation of three species included symaptric host shifts to Rhipocephalus brevicaulis in the ancestor of E. pratensis Ortea & Espinosa, 1996, and to C. prolifera in E. hamanni Krug, Vendetti & Valdes 2016, plus a niche expansion to a range of Caulerpa spp. in E. subornata Verrill, 1901. All three species are broadly sympatric across the Caribbean but are host-partitioned at a fine grain, and distinct by morphology and at nuclear loci. However, non-recombining mtDNA revealed a history of gene flow between E. pratensis and E. subornata: COI haplotypes from E. subornata were 10.4% divergent from E. pratensis haplotypes from four sites, but closely related to all E. pratensis haplotypes sampled from six Bahamian islands, indicating historical introgression and localized "mitochondrial capture." Disruptive selective likely fueled divergence and adaptation to distinct host environments, indicating ecological speciation may be an under-appreciated driver of diversification for marine herbivores as well as epibionts and other resource specialists.

Geological and climatic influences on population differentiation of the Phrynocephalus vlangalii species complex (Sauria: Agamidae) in the northern Qinghai-Tibet Plateau

Extremely heterogeneous topography and complex paleoclimate history of the Qinghai-Tibet Plateau (QTP) have a key role in promoting genetic divergence and lineage/species formation. Here, we sequenced one nuclear and three mitochondrial markers of 532 individuals from the entire range of the Phrynocephalus vlangalii species complex including two species, P. putjatai and P. vlangalii, endemic to the northern QTP. We integrated multilocus phylogeny, demographic analysis and geographic barrier detection to evaluate the population structure and dynamics. We found a new mitochondrial clade (PV-I) in the Gonghe County population of P. vlangalii, partial mitochondrial DNA replacement within P. vlangalii and complete mitochondrial DNA replacement between P. putjatai and P. vlangalii. Neutrality test, mismatch distribution analysis and Extended Bayesian Skyline Plot (EBSP) analysis all supported a significant expansion of the Qaidam Basin population of P. vlangalii (PV-II-2) from 0.091 to 0.026 Ma after Penultimate Glaciation. The uplift of the Arjin and Anyemanqen Mountains during the Kunhuang Movement (∼1.2 Ma) split populations of P. vlangalii in Akesai, Qaidam Basin and source of the Yellow River. The uplift of the Elashan Mountains during the second phase of the Qingzang Movement (∼2.5 Ma) contributed to the divergence of the Gonghe County population of P. vlangalii from other conspecific populations. The third phase of the Qingzang Movement (∼1.7 Ma) contributed to the divergence of the Xinghai population of P. vlangalii from P. putjatai and to the divergence of the northern populations of P. putjatai from the southern conspecific populations. Our data support the idea that the geological and climatic changes following the orogeny of the QTP may have promoted population differentiation and shaped the current population patterns of the P. vlangalii species complex in the northeastern QTP.

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.

A reassessment of the diversity of green toads (Bufotes) in the circum-Aegean region

The dynamic biogeography of glacial refugia may cause complex patterns of genetic admixture between parapatric taxa, which in turn can mislead their systematics, diversity, and distributions. We investigated this issue for green toads (Bufotes) inhabiting the circum-Aegean region, a biodiversity hotspot of the Eastern Mediterranean. A previous phylogeographic study based on mitochondrial and microsatellite loci identified the hybrid zone between the European (viridis) and Anatolian (sitibundus) lineages of B. viridis all over the Balkan Peninsula, but subsequent range-wide genomic analyses (>1000 SNPs) located this transition in Turkey, a thousand kilometers eastwards. In order to clarify the diversity and taxonomy of the circum-Aegean populations, we reconciled these conflicting findings by integrating previous data with pure sitibundus individuals. Our results confirmed that the viridis/sitibundus hybrid zone extends in Western Anatolia, but that southeastern European populations feature cytonuclear discordances and a high and structured microsatellite diversity. This remarkable situation may stem from a massive geographic displacement of the hybrid zone during the last glacial fluctuations, an underappreciated yet seemingly common feature among the herpetofauna of the region. Our study thus contributes to the rising view that mitochondrial DNA can be a poor predictor of current phylogeographic structure, hence the need for genomic data, especially for narrowly distributed taxa. Finally, the analyses unambiguously support the distinction of a micro-endemic clade of green toads unique to some Cyclades islands, for which we provide a formal taxonomic description.

Genetic variation and relationships among Afrotropical species of Myotis (Chiroptera: Vespertilionidae)

The genus Myotis is nearly cosmopolitan and the second-most speciose genus of mammals, but its Afrotropical members are few and poorly known. We analyzed phylogenetic and phylogeographic relationships of six of the eight known Afrotropical species using Cytb and sequences from four nuclear introns. Using Bayesian and maximum-likelihood approaches to generate single-locus, concatenated, and species trees, we confirmed prior evidence that the clade containing Afrotropical Myotis also contains both Palearctic and Indomalayan members. Additionally, we demonstrate that M. bocagii is sister to the Indian Ocean species M. anjouanensis, that this group is sister to M. tricolor and the Palearctic M. emarginatus, and find evidence suggesting that M. welwitschii is the earliest-diverging Afrotropical species and sister to the remainder. Although M. tricolor and M. welwitschii are both currently regarded as monotypic, both mitochondrial and nuclear data sets document significant, largely concordant geographic structure in each. Evidence for the distinction of two lineages within M. tricolor is particularly strong. On the other hand, geographic structure is lacking in M. bocagii, despite the current recognition of two subspecies in that species. Additional geographic sampling (especially at or near type localities), finer-scale sampling (especially in zones of sympatry), and integrative taxonomic assessments will be needed to better document this radiation and refine its nomenclature.

Porous barriers? Assessment of gene flow within and among sympatric long-eared bat species

Species are the basic units for measuring biodiversity and for comprehending biological interactions. Yet, their delineation is often contentious, especially in groups that are both diverse and phenotypically conservative. Three cryptic species of long-eared bats, Plecotus auritus, P. austriacus, and P. macrobullaris, co-occur over extensive areas of Western Europe. The latter is a fairly recent discovery, questioning the overall diversity of the entire Plecotus complex. Yet, high morphological and acoustic similarities compromise the reliable identification of long-eared bats in the field. We postulate that such extensive phenotypic overlap, along with the recurrent observation of morphologically intermediate individuals, may hide rampant interspecific hybridization. Based on a geographic sampling centered on areas of sympatry in the Alps and Corsica, we assessed the level of reproductive isolation of these three Plecotus species with mitochondrial and nuclear markers, looking at both inter- and intraspecific genetic population structuring. No sign of hybridization was detected between these three species that appear well separated biologically. Genetic structuring of populations, however, reflected different species-specific responses to environmental connectivity, that is, to the presence of orographic or sea barriers. While the Alpine range and the Ligurian Sea coincided with sharp genetic discontinuities in P. macrobullaris and P. austriacus, the more ubiquitous P. auritus showed no significant population structuration. There were clear phylogeographic discrepancies between microsatellite and mitochondrial markers at the intraspecific level, however, which challenges the reliance on simple barcoding approaches for the delineation of sound conservation units.

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.

Bat Biology, Genomes, and the Bat1K Project: To Generate Chromosome-Level Genomes for All Living Bat Species

Bats are unique among mammals, possessing some of the rarest mammalian adaptations, including true self-powered flight, laryngeal echolocation, exceptional longevity, unique immunity, contracted genomes, and vocal learning. They provide key ecosystem services, pollinating tropical plants, dispersing seeds, and controlling insect pest populations, thus driving healthy ecosystems. They account for more than 20% of all living mammalian diversity, and their crown-group evolutionary history dates back to the Eocene. Despite their great numbers and diversity, many species are threatened and endangered. Here we announce Bat1K, an initiative to sequence the genomes of all living bat species (n∼1,300) to chromosome-level assembly. The Bat1K genome consortium unites bat biologists (>148 members as of writing), computational scientists, conservation organizations, genome technologists, and any interested individuals committed to a better understanding of the genetic and evolutionary mechanisms that underlie the unique adaptations of bats. Our aim is to catalog the unique genetic diversity present in all living bats to better understand the molecular basis of their unique adaptations; uncover their evolutionary history; link genotype with phenotype; and ultimately better understand, promote, and conserve bats. Here we review the unique adaptations of bats and highlight how chromosome-level genome assemblies can uncover the molecular basis of these traits. We present a novel sequencing and assembly strategy and review the striking societal and scientific benefits that will result from the Bat1K initiative.

Pleistocene range shifts, refugia and the origin of widespread species in western Palaearctic water beetles

Quaternary glacial cycles drove major shifts in both the extent and location of the geographical ranges of many organisms. During glacial maxima, large areas of central and northern Europe were inhospitable to temperate species, and these areas are generally assumed to have been recolonized during interglacials by range expansions from Mediterranean refugia. An alternative is that this recolonization was from non-Mediterranean refugia, in central Europe or western Asia, but data on the origin of widespread central and north European species remain fragmentary, especially for insects. We studied three widely distributed lineages of freshwater beetles (the Platambus maculatus complex, the Hydraena gracilis complex, and the genus Oreodytes), all restricted to running waters and including both narrowly distributed southern endemics and widespread European species, some with distributions spanning the Palearctic. Our main goal was to determine the role of the Pleistocene glaciations in shaping the diversification and current distribution of these lineages. We sequenced four mitochondrial and two nuclear genes in populations drawn from across the ranges of these taxa, and used Bayesian probabilities and Maximum Likelihood to reconstruct their phylogenetic relationships, age and geographical origin. Our results suggest that all extant species in these groups are of Pleistocene origin. In the H. gracilis complex, the widespread European H. gracilis has experienced a rapid, recent range expansion from northern Anatolia, to occupy almost the whole of Europe. However, in the other two groups widespread central and northern European taxa appear to originate from central Asia, rather than the Mediterranean. These widespread species of eastern origin typically have peripherally isolated forms in the southern Mediterranean peninsulas, which may be remnants of earlier expansion-diversification cycles or result from incipient isolation of populations during the most recent Holocene expansion. The accumulation of narrow endemics of such lineages in the Mediterranean may result from successive cycles of range expansion, with subsequent speciation (and local extinction in glaciated areas) through multiple Pleistocene climatic cycles.

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.

Introgression and selection shaped the evolutionary history of sympatric sister-species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister-species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site-associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister-species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174-0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.

Complex within a Complex: Integrative Taxonomy Reveals Hidden Diversity in Cicadetta brevipennis (Hemiptera: Cicadidae) and Unexpected Relationships with a Song Divergent Relative

Multiple sources of data in combination are essential for species delimitation and classification of difficult taxonomic groups. Here we investigate a cicada taxon with unusual cryptic diversity and we attempt to resolve seemingly contradictory data sets. Cicada songs act as species-specific premating barriers and have been used extensively to reveal hidden taxonomic diversity in morphologically similar species. The Palaearctic Cicadetta montana species complex is an excellent example where distinct song patterns have disclosed multiple recently described species. Indeed, two taxa turned out to be especially diverse in that they form a “complex within the complex”: the Cicadetta cerdaniensis song group (four species studied previously) and Cicadetta brevipennis (examined in details here). Based on acoustic, morphological, molecular, ecological and spatial data sampled throughout their broad European distribution, we find that Cicadetta brevipennis s. l. comprises five lineages. The most distinct lineage is identified as Cicadetta petryi Schumacher, 1924, which we re-assign to the species level. Cicadetta brevipennis litoralis Puissant & Hertach ssp. n. and Cicadetta brevipennis hippolaidica Hertach ssp. n. are new to science. The latter hybridizes with Cicadetta brevipennis brevipennis Fieber, 1876 at a zone inferred from intermediate song patterns. The fifth lineage requires additional investigation. The C. cerdaniensis and the C. brevipennis song groups exhibit characteristic, clearly distinct basic song patterns that act as reproductive barriers. However, they remain completely intermixed in the Bayesian and maximum likelihood COI and COII mitochondrial DNA phylogenies. The closest relative of each of the four cerdaniensis group species is a brevipennis group taxon. In our favoured scenario the phylogenetic pairs originated in common Pleistocene glacial refuges where the taxa speciated and experienced sporadic inter-group hybridization leading to extensive introgression and mitochondrial capture.

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.

The Value of Molecular vs. Morphometric and Acoustic Information for Species Identification Using Sympatric Molossid Bats

A fundamental condition for any work with free-ranging animals is correct species identification. However, in case of bats, information on local species assemblies is frequently limited especially in regions with high biodiversity such as the Neotropics. The bat genus Molossus is a typical example of this, with morphologically similar species often occurring in sympatry. We used a multi-method approach based on molecular, morphometric and acoustic information collected from 962 individuals of Molossus bondae, M. coibensis, and M. molossus captured in Panama. We distinguished M. bondae based on size and pelage coloration. We identified two robust species clusters composed of M. molossus and M. coibensis based on 18 microsatellite markers but also on a more stringently determined set of four markers. Phylogenetic reconstructions using the mitochondrial gene co1 (DNA barcode) were used to diagnose these microsatellite clusters as M. molossus and M. coibensis. To differentiate species, morphological information was only reliable when forearm length and body mass were combined in a linear discriminant function (95.9% correctly identified individuals). When looking in more detail at M. molossus and M. coibensis, only four out of 13 wing parameters were informative for species differentiation, with M. coibensis showing lower values for hand wing area and hand wing length and higher values for wing loading. Acoustic recordings after release required categorization of calls into types, yielding only two informative subsets: approach calls and two-toned search calls. Our data emphasizes the importance of combining morphological traits and independent genetic data to inform the best choice and combination of discriminatory information used in the field. Because parameters can vary geographically, the multi-method approach may need to be adjusted to local species assemblies and populations to be entirely informative.

Genome-wide data reveal cryptic diversity and genetic introgression in an Oriental cynopterine fruit bat radiation

Background

The Oriental fruit bat genus Cynopterus, with several geographically overlapping species, presents an interesting case study to evaluate the evolutionary significance of coexistence versus isolation. We examined the morphological and genetic variability of congeneric fruit bats Cynopterus sphinx and C. brachyotis using 405 samples from two natural contact zones and 17 allopatric locations in the Indian subcontinent; and investigated the population differentiation patterns, evolutionary history, and the possibility of cryptic diversity in this species pair.

Results

Analysis of microsatellites, cytochrome b gene sequences, and restriction digestion based genome-wide data revealed that C. sphinx and C. brachyotis do not hybridize in contact zones. However, cytochrome b gene sequences and genome-wide SNP data helped uncover a cryptic, hitherto unrecognized cynopterine lineage in northeastern India coexisting with C. sphinx. Further analyses of shared variation of SNPs using Patterson’s D statistics suggest introgression between this lineage and C. sphinx. Multivariate analyses of morphology using genetically classified grouping confirmed substantial morphological overlap between C. sphinx and C. brachyotis, specifically in the high elevation contact zones in southern India.

Conclusion

Our results uncover novel diversity and detect a pattern of genetic introgression in a cryptic radiation of bats, demonstrating the complicated nature of lineage diversification in this poorly understood taxonomic group. Our results highlight the importance of genome-wide data to study evolutionary processes of morphologically similar species pairs. Our approach represents a significant step forward in evolutionary research on young radiations of non-model species that may retain the ability of interspecific gene flow.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0599-y) contains supplementary material, which is available to authorized users.

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.

Limpets of the genus Cellana (Patellogastropoda) from Pakistan, North Arabian Sea: species identification based on DNA sequencing

The true limpets are found in the intertidal zone of the rocky shores of Pakistan, North Arabian Sea. Partial sequence of the mitochondrial cytochrome oxidase I was used to estimate the degree of genetic differentiation among the morphological forms of Cellana, which were considered as three separate species earlier in Pakistan. The study revealed that the three morphs of Cellana on COI sequence generated a single haplotype and matched with the COI sequence of Cellana karachiensis. This point out the phenotypic plasticity between the proposed species.

Parasite diversity of European Myotis species with special emphasis on Myotis myotis (Microchiroptera, Vespertilionidae) from a typical nursery roost

Background

Bats belong to one of the most species-rich orders within the Mammalia. They show a worldwide distribution, a high degree of ecological diversification as well as a high diversity of associated parasites and pathogens. Despite their prominent and unique role, the knowledge of their parasite-host-relationships as well as the mechanisms of co-evolutionary processes are, partly due to strict conservation regulations, scarce.

Methods

Juvenile specimens of the greater mouse-eared bat (Myotis myotis) from a roosting colony in Gladenbach (Hesse, Germany) were examined for their metazoan endo-and ectoparasite infections and pathogens. Morphometric data were recorded and the individuals were checked for Lyssavirus-specific antigen using a direct immunofluorescence test. For unambiguous species identification, the bats were analysed by cyt-b sequence comparison.

Results

Myotis myotis were parasitized by the six insect and arachnid ectoparasite species, i.e. Ixodes ricinus, Ischnopsyllus octactenus, Ichoronyssus scutatus, Steatonyssus periblepharus, Spinturnix myoti and Cimex dissimilis. Additionally, the nematode Molinostrongylus alatus and the cestode Vampirolepis balsaci were recorded. Each bat was parasitized by at least four species. The parasites showed partially extreme rates of infection, never recorded before, with more than 1,440 parasites per single host. Ichoronyssus scutatus, Steatonyssus periblepharus, Vampirolepis balsaci and Molinostrongylus alatus are recorded for the first time in Germany. A checklist for Europe is presented containing records of 98 parasite species of 14 Myotis species.

Conclusions

The Myotis myotis from Gladenbach (Hesse, Germany) were parasitized by a diverse parasite fauna with high infestation rates. We assume that in juvenile Myotis the number of parasites is generally higher than in adults due to only later acquired immune competence and behavioural adaptations. Our results revealed new insights into parasite fauna of M. myotis and European bats in general. The finding of endoparasitic cyclophyllidean cestodes that have a two-host lifecycle is, considering the stationary behaviour of the juvenile bats, rather unusual and suggests a non-predatory transmission mechanism (e.g. via autoinfection).

A new insight gained from the collated literature was that the European wide composition of the Myotis parasite fauna is dominated by a few specific taxonomic groups in Europe.

The tail plays a major role in the differing manoeuvrability of two sibling species of mouse-eared bats (Myotis myotis and Myotis blythii)

Two sympatrically occurring bat species, the greater mouse-eared bat (Myotis myotis (Borkhausen, 1797)) and the lesser mouse-eared bat (Myotis blythii (Tomes, 1857)) (Chiroptera, Vespertillionidae), share numerous similarities in morphology, roosting behaviour, and echolocation and are often difficult to distinguish. However, despite these similarities, their foraging behaviour is noticeably different. Our aim was to examine the extent to which these different foraging strategies reflect morphological adaptation. We assessed whether the morphology of the wing, body, and tail differed between M. myotis and M. blythii. In addition, in a laboratory experiment involving an obstacle course, we compared differences in manoeuvrability by relating them to our morphological measurements. The two species differed in their overall size, wing-tip shape, and tail-to-body length ratio. The generally smaller sized M. blythii performed better in the obstacle course and was therefore considered to be more manoeuvrable. Although differences in wing-tip shape were observed, we found the most important characteristic affecting manoeuvrability in both species to be the tail-to-body length ratio. Additionally, when we compared two bats with injured wing membranes with unharmed bats of the same species, we found no difference in manoeuvrability, even when the wing shape was asymmetric. We therefore postulate that morphometric differences between the two species in their overall size and, more importantly, in their tail-to-body length ratio are the main physical characteristics providing proof of adaptation to different foraging and feeding strategies.

Mitochondrial data are not suitable for resolving placental mammal phylogeny

Mitochondrial data have traditionally been used in reconstructing a variety of species phylogenies. The low rates of recombination and thorough characterization of mitochondrial data across vertebrate species make it a particularly attractive phylogenetic marker. The relatively low number of fully sequenced mammal genomes and the lack of extensive sampling within Superorders have posed a serious problem for reaching agreement on the placement mammal species. The use of mitochondrial data sequences from large numbers of mammals could serve to circumvent the taxon-sampling deficit. Here we assess the suitability of mitochondrial data as a phylogenetic marker in mammal phylogenetics. MtDNA datasets of mammal origin have been filtered as follows: (i) we have sampled sparsely across the phylogenetic tree, (ii) we have constrained our sampling to genes with high taxon coverage, (iii) we have categorised rates across sites in a phylogeny independent manner and have removed fast evolving sites, and (iv), we have sampled from very shallow divergence times to reduce phylogenetic conflict. However, topologies obtained using these filters are not consistent with previous studies and are discordant across different genes. Individual mitochondrial genes, and indeed all mitochondrial genes analysed as a supermatrix, resulted in poor resolution of the species phylogeny. Overall, our study highlights the limitations of mitochondrial data, not only for resolving deep divergences and but also for shallow divergences in the mammal phylogeny.

A cryptic species of the Tylonycteris pachypus complex (Chiroptera: Vespertilionidae) and its population genetic structure in southern China and nearby regions

Three distinct bamboo bat species (Tylonycteris) are known to inhabit tropical and subtropical areas of Asia, i.e., T. pachypus, T. robustula, and T. pygmaeus. This study performed karyotypic examinations of 4 specimens from southern Chinese T. p. fulvidus populations and one specimen from Thai T. p. fulvidus population, which detected distinct karyotypes (2n=30) compared with previous karyotypic descriptions of T. p. pachypus (2n=46) and T. robustula (2n=32) from Malaysia. This finding suggested a cryptic Tylonycteris species within T. pachypus complex in China and Thailand. Morphometric studies indicated the difficulty in distinguishing the cryptic species and T. p. pachypus from Indonesia apart from the external measurements, which might be the reason for their historical misidentification. Based on 623 bp mtDNA COI segments, a phylogeographic examination including T. pachypus individuals from China and nearby regions, i.e., Vietnam, Laos, and Cambodia, was conducted to examine the population genetic structure. Genealogical and phylogeographical results indicated that at least two diverged lineages existed in these regions (average 3.4 % of Kimura 2-parameter distances) and their population structure did not match the geographic pattern. These results suggested that at least two historical colonizations have occurred by the cryptic species. Furthermore, through integration of traditional and geometric morphological results, morphological differences on zygomatic arches, toothrows and bullae were detected between two lineages in China. Given the similarity of vegetation and climate of Guangdong and Guangxi regions, we suggested that such differences might be derived from their historical adaptation or distinct evolutionary history rather than the differences of habitats they occurred currently.

Social learning within and across species: information transfer in mouse-eared bats

Social learning describes information transfer between individuals through observation or direct interaction. Bats can live and forage in large groups, sometimes comprising several species, and are thus well suited for investigations of both intraspecific and interspecific information transfer. Although social learning has been documented within several bat species, it has not been shown to occur between species. Furthermore, it is not fully understood what level of interaction between individuals is necessary for social learning in bats. We address these questions by comparing the efficiency of observation versus interaction in intraspecific social learning and by considering interspecific social learning in sympatric bat species. Observers learned from demonstrators to identify food sources using a light cue. We show that intraspecific social learning exists in the greater mouse-eared bat (Myotis myotis (Borkhausen, 1797)) and that direct interaction with a demonstrator more efficiently leads to information transfer than observational learning alone. We also found evidence for interspecific information transfer from M. myotis to the lesser mouse-eared bat (Myotis oxygnathus Monticelli, 1885). Additionally, we opportunistically retested one individual that we recaptured from the wild 1 year after initial learning and found long-term memory of the trained association. Our study adds to the understanding of learning, information transfer, and long-term memory in wild-living animals.