The evolution of sensory divergence in the context of limited gene flow in the bumblebee bat

Comparative Genome Microsynteny Illuminates the Fast Evolution of Nuclear Mitochondrial Segments (NUMTs) in Mammals

The escape of DNA from mitochondria into the nuclear genome (nuclear mitochondrial DNA, NUMT) is an ongoing process. Although pervasively observed in eukaryotic genomes, their evolutionary trajectories in a mammal-wide context are poorly understood. The main challenge lies in the orthology assignment of NUMTs across species due to their fast evolution and chromosomal rearrangements over the past 200 million years. To address this issue, we systematically investigated the characteristics of NUMT insertions in 45 mammalian genomes and established a novel, synteny-based method to accurately predict orthologous NUMTs and ascertain their evolution across mammals. With a series of comparative analyses across taxa, we revealed that NUMTs may originate from nonrandom regions in mtDNA, are likely found in transposon-rich and intergenic regions, and unlikely code for functional proteins. Using our synteny-based approach, we leveraged 630 pairwise comparisons of genome-wide microsynteny and predicted the NUMT orthology relationships across 36 mammals. With the phylogenetic patterns of NUMT presence-and-absence across taxa, we constructed the ancestral state of NUMTs given the mammal tree using a coalescent method. We found support on the ancestral node of Fereuungulata within Laurasiatheria, whose subordinal relationships are still controversial. This study broadens our knowledge on NUMT insertion and evolution in mammalian genomes and highlights the merit of NUMTs as alternative genetic markers in phylogenetic inference.

A new species of the Miniopterus australis species complex (Chiroptera: Miniopteridae) from the Western Ghats, India

The genus Miniopterus is highly diverse in cryptic species. Based on integrated approaches of morphometrics, echolocation call analysis, and molecular phylogenetics, we present evidence of a hitherto undescribed species, Srini's Bent-winged bat Miniopterus srinii sp. nov.. The new species, found in the Western Ghats of southern Karnataka, India, closely resembles the Small Bent-winged bat Miniopterus pusillus, found elsewhere in Karnataka, Kerala, and Tamil Nadu, India. Although it shows external similarity with Miniopterus pusillus sensu stricto from the Nicobar Islands, the new species is relatively larger and distinct in craniodental measurements and echolocation call parameters, despite slight overlaps. Importantly, Miniopterus srinii exhibits morphological convergence with Miniopterus pusillus as despite similarities in morphologies, there is a significant genetic distance of 10.84 ± 0.22%. The new cryptic species shares distribution with Miniopterus phillipsi and Miniopterus pusillus in the Western Ghats, further highlighting the need to study both the genus' cryptic diversity, and the region's conservation importance.

Genetic Characterization and Insular Habitat Enveloping of Endangered Leaf-Nosed Bat, Hipposideros nicobarulae (Mammalia: Chiroptera) in India: Phylogenetic Inference and Conservation Implication

The Nicobar leaf-nosed Bat (Hipposideros nicobarulae) was described in the early 20th century; however, its systematic classification has been debated for over 100 years. This endangered and endemic species has achieved species status through morphological data in the last 10 years. However, the genetic information and phylogenetic relationships of H. nicobarulae remain neglected. The generated mitochondrial cytochrome b gene (mtCytb) sequences (438 bp) of H. nicobarulae contains 53.42-53.65% AT composition and 1.82% variable sites. The studied species, H. nicobarulae maintains an 8.1% to 22.6% genetic distance from other Hipposideros species. The genetic divergence estimated in this study is congruent with the concept of gene speciation in bats. The Bayesian and Maximum-Likelihood phylogenies clearly discriminated all Hipposideros species and showed a sister relationship between H. nicobarulae and H. cf. antricola. Current mtCytb-based investigations of H. nicobarulae have confirmed the species status at the molecular level. Further, the MaxEnt-based species distribution modelling illustrates the most suitable habitat of H. nicobarulae (294 km²), of which the majority (171 km²) is located on Great Nicobar Island. The present study suggests rigorous sampling across the range, taxonomic coverage, the generation of multiple molecular markers (mitochondrial and nuclear), as well as more ecological information, which will help in understanding population genetic structure, habitat suitability, and the implementation of appropriate conservation action plans for H. nicobarulae and other Hipposideros species.

Evidence for rapid divergence of sensory systems between Texas populations of the Mexican tetra (Astyanax mexicanus)

Population divergence is often quantified using phenotypic variation. However, because sensory abilities are more difficult to discern, we have little information on the plasticity and rate of sensory change between different environments. The Mexican tetra (Astyanax mexicanus) is a fish distributed throughout Southern Texas and Northern Mexico and has evolved troglomorphic phenotypes, such as vestigial eyes and reduced pigmentation, when surface ancestors invaded caves in the past several hundred thousand years. In the early 1900s, surface A. mexicanus were introduced to the karstic Edwards-Trinity Aquifer in Texas. Subsequent cave colonization of subterranean environments resulted in fish with phenotypic and behavioral divergence from their surface counterparts, allowing examination of how new environments lead to sensory changes. We hypothesized that recently introduced cave populations would be more sensitive to light and sound when compared to their surface counterparts. We quantified divergence using auditory evoked potentials (AEPs) and particle acceleration levels (PALs) to measure differences in sound sensitivity, and electroretinography (ERGs) to measure light sensitivity. We also compared these results to measurements taken from native populations and lab-born individuals of the introduced populations. Honey Creek Cave fish were significantly more sensitive than proximate Honey Creek surface fish to sound pressure levels between 0.6 and 0.8 kHz and particle acceleration levels between 0.4 and 0.8 kHz. Pairwise differences were found between San Antonio Zoo surface and the facultative subterranean San Pedro Springs and Blue Hole populations, which exhibited more sensitivity to particle acceleration levels between 0.5 and 0.7 kHz. Electroretinography results indicate no significant differences between populations, although Honey Creek Cave fish may be trending toward reduced visual sensitivity. Auditory thresholds between wild-caught and lab-raised populations of recently invaded fish show significant differences in sensitivity, suggesting that these traits are plastic. Collectively, while these results may point to the rapid divergence of A. mexicanus in cave habitats, it also highlights the responsive plasticity of A. mexicanus auditory system to disparate environments.

The EDGE2 protocol: Advancing the prioritisation of Evolutionarily Distinct and Globally Endangered species for practical conservation action

The conservation of evolutionary history has been linked to increased benefits for humanity and can be captured by phylogenetic diversity (PD). The Evolutionarily Distinct and Globally Endangered (EDGE) metric has, since 2007, been used to prioritise threatened species for practical conservation that embody large amounts of evolutionary history. While there have been important research advances since 2007, they have not been adopted in practice because of a lack of consensus in the conservation community. Here, building from an interdisciplinary workshop to update the existing EDGE approach, we present an "EDGE2" protocol that draws on a decade of research and innovation to develop an improved, consistent methodology for prioritising species conservation efforts. Key advances include methods for dealing with uncertainty and accounting for the extinction risk of closely related species. We describe EDGE2 in terms of distinct components to facilitate future revisions to its constituent parts without needing to reconsider the whole. We illustrate EDGE2 by applying it to the world's mammals. As we approach a crossroads for global biodiversity policy, this Consensus View shows how collaboration between academic and applied conservation biologists can guide effective and practical priority-setting to conserve biodiversity.

Divergence in interspecific and intersubspecific gene expression between two closely related horseshoe bats (Rhinolophus)

Closely related species have been used as representative systems to investigate the genetic mechanisms involved in the early stages of species differentiation. Previous studies have indicated that variation in gene expression might be a sensitive indicator of initial species divergence, although the role of expression divergence, and especially that associated with phenotypic variation remained relatively undefined. For three organs (cochlea, brain, and liver) from two closely related bat species (Rhinolophus siamensis and R. episcopus), the interspecific and intersubspecific gene expression profiles were compared using transcriptomics in this study. Striking organ specificity of expression was observed, and expression profiles exhibited similarities between cochlea and brain tissues. Numerous differentially expressed genes (DEGs) were identified for each organ in the interspecific comparison (cochlea/brain/liver: 1,069/647/692) and intersubspecific comparison (608/528/368). Functional enrichment analysis indicated vital variation in expression related to the immune system, ion activities, neuronal function, and multisensory system regulation in both comparisons. DEGs relevant to the variation in echolocation calls (RF) were found, and some of them were involved in the pivotal patterns of expression variation. The regulation of immune, ion channel, neural activity, and sophisticated sensory functions at the expression level might be key mechanisms in the early species divergence of bats, and the expression variation related to acoustical signal could have played a crucial part. This study expands our knowledge of gene expression and patterns of variation for three key organs to echolocation at both the interspecific and intersubspecific levels. Further, the framework described here provides insight into the genetic basis of phenotypic variation during the incipient stage of species differentiation.

Extensive Adaptive Variation in Gene Expression within and between Closely Related Horseshoe Bats (Chiroptera, Rhinolophus) Revealed by Three Organs

In the process of species differentiation and adaption, the relative influence of natural selection on gene expression variation often remains unclear (especially its impact on phenotypic divergence). In this study, we used differentially expressed genes from brain, cochlea, and liver samples collected from two species of bats to determine the gene expression variation forced by natural selection when comparing at the interspecific (Rhinolophus siamensis and R. episcopus episcopus) and the intraspecific (R. e. episcopus and R. episcopus spp.) levels. In both cases, gene expression variation was extensively adaptive (>66.0%) and mainly governed by directional selection, followed by stabilizing selection, and finally balancing selection. The expression variation related to acoustic signals (resting frequency, RF) and body size (forearm length, FA) was also widely governed by natural selection (>69.1%). Different functional patterns of RF- or FA-related adaptive expression variation were found between the two comparisons, which manifested as abundant immune-related regulations between subspecies (indicating a relationship between immune response and phenotypic adaption). Our study verifies the extensive adaptive expression variation between both species and subspecies and provides insight into the effects of natural selection on species differentiation and adaptation as well as phenotypic divergence at the expression level.

Geographic Variation in Social Vocalizations of the Great Himalayan Leaf-Nosed Bat, Hipposideros armiger: Acoustic Overflow Across Population Boundaries

Bat populations employ rich vocal repertoires for social communication in addition to emitting sound pulses for echolocation. Acoustic parameters of echolocation pulses can vary with the context in which they are emitted, and also with the individual and across populations as a whole. The acoustic parameters of social vocalizations, or “calls”, also vary with the individual and context, but not much is known about their variation across populations at different geographic locations. Here, we leveraged the detailed acoustic classification of social vocalizations available for the Great Himalayan leaf-nosed bat, Hipposideros armiger, to examine geographic variation in five commonly emitted simple syllable types. We hypothesized that individuals within geographically dispersed populations communicate using spectrographically similar constructs or “syllable types”. We also examined whether call syllables vary discordantly with the correlation pattern observed for echolocation pulses across those same geographic regions. Furthermore, we postulated that the acoustic boundaries of a syllable type are not uniquely constrained to its variation within a particular population of the same subspecies. To test our hypotheses, we obtained recordings of social calls of H. a. armiger from nine locations within the oriental region. These locations were consolidated into five geographic regions based on previously established region-specific differences in the peak frequency of echolocation pulses. A multivariate cluster analysis established that unlike echolocation pulses, syllable types exhibit a relatively large variance. Analysis of this variance showed significant differences in Least Squares Means estimates, establishing significant population-level differences in the multiparametric means of individual syllable types across geographic regions. Multivariate discriminant analysis confirmed the presence of region-specific centroids for different syllable constructs, but also showed a large overlap of their multiparametric boundaries across geographic regions. We propose that despite differences in the population-specific core construct of a syllable type, bats maximize acoustic variation across individuals within a population irrespective of its overflow and overlap with other populations.

Micro-endemic species of snails and amphipods show population genetic structure across very small geographic ranges

Understanding variation in population genetic structure, even across small distances and for species with extremely limited ranges, is critical for conservation planning and the development of effective management strategies for imperiled species. Organisms that occupy the same geographic extent can maintain different population structures, ranging from highly diverged to panmictic. Such differences can result from differences in biological characteristics such as dispersal ability or demographic history. We used microsatellite loci to evaluate population genetic structure and variation of four desert spring invertebrates having high to low dispersal ability: the lung snail Physa acuta, two species of gilled snails (Juturnia kosteri and Pyrgulopsis roswellensis; family Hydrobiidae) and the amphipod Gammarus desperatus. The study location represents entire species ranges for the micro-endemic hydrobiids and G. desperatus, while P. acuta is ubiquitous throughout much of North America. We found little evidence of significant population genetic structure for P. acuta and J. kosteri, but much more for P. roswellensis and G. desperatus. Our results demonstrate differences in habitat preference and/or dispersal ability between the species. This information provides insight into how gene flow shapes varying population genetic structure between species across small spatial scales (<100 km2). Most importantly, our results suggest that conservation agencies should not consider these micro-endemic species to be composed of single populations, but rather, that management plans for such species should account for population genetic variation across the species' ranges.

Morphological and molecular confirmation of the common pipistrelle bat, Pipistrellus pipistrellus Schreber, 1774 (Vespertilionidae: Chiroptera), in Xinjiang, China

Species identification is pivotal in taxonomy, systematics, evolutionary biology and conservation biology. We collected bats that died of natural causes in Shihezi city, Xinjiang, China, and carried out morphological and genetic identification. Morphologically, all individuals were adults/subadults or juveniles of Pipistrellus pipistrellus. We found one haplotype for the mitochondrial gene ND1 and five for the mitochondrial gene cytochrome b (Cytb) among six specimens. Phylogenetically, all the Cytb sequences grouped with P. pipistrellus. We confirm this species’ occurrence in Xinjiang, China.

Using historical genome-wide DNA to unravel the confused taxonomy in a songbird lineage that is extinct in the wild

Urgent conservation action for terminally endangered species is sometimes hampered by taxonomic uncertainty, especially in illegally traded animals that are often cross-bred in captivity. To overcome these problems, we used a genomic approach to analyze historical DNA from museum samples across the Asian Pied Starling (Gracupica contra) complex in tropical Asia, a popular victim of the ongoing songbird crisis whose distinct Javan population ("Javan Pied Starling") is extinct in the wild and subject to admixture in captivity. Comparing genomic profiles across the entire distribution, we detected three deeply diverged lineages at the species level characterized by a lack of genomic intermediacy near areas of contact. Our study demonstrates that the use of historical DNA can be instrumental in delimiting species in situations of taxonomic uncertainty, especially when modern admixture may obfuscate species boundaries. Results of our research will enable conservationists to commence a dedicated ex situ breeding program for the Javan Pied Starling, and serve as a blueprint for similar conservation problems involving terminally endangered species subject to allelic infiltration from close congeners.

Extensive Genetic Diversity of Polyomaviruses in Sympatric Bat Communities: Host Switching versus Coevolution

Polyomaviruses (PyVs) are small DNA viruses carried by diverse vertebrates. The evolutionary relationships of viruses and hosts remain largely unclear due to very limited surveillance in sympatric communities. In order to investigate whether PyVs can transmit among different mammalian species and to identify host-switching events in the field, we conducted a systematic study of a large collection of bats (n = 1,083) from 29 sympatric communities across China which contained multiple species with frequent contact. PyVs were detected in 21 bat communities, with 192 PyVs identified in 186 bats from 15 species within 6 families representing at least 28 newly described PyVs. Surveillance results and phylogenetic analyses surprisingly revealed three interfamily PyV host-switching events in these sympatric bat communities: two distinct PyVs were identified in two bat species in restricted geographical locations, while another PyV clustered phylogenetically with PyVs carried by bats from a different host family. Virus-host relationships of all discovered PyVs were also evaluated, and no additional host-switching events were found. PyVs were identified in different horseshoe bat species in sympatric communities without observation of host-switching events, showed high genomic identities, and clustered with each other. This suggested that even for PyVs with high genomic identities in closely related host species, the potential for host switching is low. In summary, our findings revealed that PyV host switching in sympatric bat communities can occur but is limited and that host switching of bat-borne PyVs is relatively rare on the predominantly evolutionary background of codivergence with their hosts.IMPORTANCE Since the discovery of murine polyomavirus in the 1950s, polyomaviruses (PyVs) have been considered highly host restricted in mammals. Sympatric bat communities commonly contain several different bat species in an ecological niche facilitating viral transmission, and they therefore represent a model to identify host-switching events of PyVs. In this study, we screened PyVs in a large number of bats in sympatric communities from diverse habitats across China. We provide evidence that cross-species bat-borne PyV transmission exists, though is limited, and that host-switching events appear relatively rare during the evolutionary history of these viruses. PyVs with close genomic identities were also identified in different bat species without host-switching events. Based on these findings, we propose an evolutionary scheme for bat-borne PyVs in which limited host-switching events occur on the background of codivergence and lineage duplication, generating the viral genetic diversity in bats.

Integrating population genetics to define conservation units from the core to the edge of Rhinolophus ferrumequinum western range

The greater horseshoe bat (Rhinolophus ferrumequinum) is among the most widespread bat species in Europe but it has experienced severe declines, especially in Northern Europe. This species is listed Near Threatened in the European IUCN Red List of Threatened Animals, and it is considered to be highly sensitive to human activities and particularly to habitat fragmentation. Therefore, understanding the population boundaries and demographic history of populations of this species is of primary importance to assess relevant conservation strategies. In this study, we used 17 microsatellite markers to assess the genetic diversity, the genetic structure, and the demographic history of R. ferrumequinum colonies in the western part of its distribution. We identified one large population showing high levels of genetic diversity and large population size. Lower estimates were found in England and northern France. Analyses of clustering and isolation by distance suggested that the Channel and the Mediterranean seas could impede R. ferrumequinum gene flow. These results provide important information to improve the delineation of R. ferrumequinum management units. We suggest that a large management unit corresponding to the population ranging from Spanish Basque Country to northern France must be considered. Particular attention should be given to mating territories as they seem to play a key role in maintaining high levels of genetic mixing between colonies. Smaller management units corresponding to English and northern France colonies must also be implemented. These insular or peripheral colonies could be at higher risk of extinction in the near future.

A continent-scale study of the social structure and phylogeography of the bent-wing bat, Miniopterus schreibersii (Mammalia: Chiroptera), using new microsatellite data

Miniopterus schreibersii is a cave-dwelling bat species with a wide distribution in the western Palearctic spanning southern and central Europe, North Africa, Anatolia, the Middle East, and the Caucasus. We investigated the social structure and its effects on the genetic makeup of this species, using 10 nuclear microsatellite markers and a partial fragment of the mitochondrial cytochrome b gene. Samples were examined from the species' entire circum-Mediterranean range. Local structuring that was previously detected among populations of M. schreibersii using mitochondrial markers was not observed for microsatellite markers, indicating male-biased dispersal for the species. Some support was found for postglacial expansions in Europe, with Anatolia potentially acting as the primary refugium during the Last Glacial Maximum (LGM). However, support for this hypothesis is not as strong as that previously detected using mitochondrial DNA markers. This is likely due to the diminishing effect of male-mediated dispersal, replenishing the nuclear diversity faster than the mitochondrial diversity in regions that are relatively far from the glacial refugia.

Scientists' Warning on the Conservation of Subterranean Ecosystems

In light of recent alarming trends in human population growth, climate change, and other environmental modifications, a “Warning to humanity” manifesto was published in BioScience in 2017. This call reiterated most of the ideas originally expressed by the Union of Concerned Scientists in 1992, including the fear that we are “pushing Earth's ecosystems beyond their capacities to support the web of life.” As subterranean biologists, we take this opportunity to emphasize the global importance and the conservation challenges associated with subterranean ecosystems. They likely represent the most widespread nonmarine environments on Earth, but specialized subterranean organisms remain among the least documented and studied. Largely overlooked in conservation policies, subterranean habitats play a critical role in the function of the web of life and provide important ecosystem services. We highlight the main threats to subterranean ecosystems and propose a set of effective actions to protect this globally important natural heritage.

A coalescent-based estimator of genetic drift, and acoustic divergence in the Pteronotus parnellii species complex

Determining the processes responsible for phenotypic variation is one of the central tasks of evolutionary biology. While the importance of acoustic traits for foraging and communication in echolocating mammals suggests adaptation, the seldom-tested null hypothesis to explain trait divergence is genetic drift. Here we derive FST values from multi-locus coalescent isolation-with-migration models, and couple them with estimates of quantitative trait divergence, or PST, to test drift as the evolutionary process responsible for phenotypic divergence in island populations of the Pteronotus parnellii species complex. Compared to traditional comparisons of PST to FST, the migration-based estimates of FST are unidirectional instead of bidirectional, simultaneously integrate variation among loci and individuals, and posterior densities of PST and FST can be compared directly. We found the evolution of higher call frequencies is inconsistent with genetic drift for the Hispaniolan population, despite many generations of isolation from its Puerto Rican counterpart. While the Hispaniolan population displays dimorphism in call frequencies, the higher frequency of the females is incompatible with sexual selection. Instead, cultural drift toward higher frequencies among Hispaniolan females might explain the divergence. By integrating Bayesian coalescent and trait analyses, this study demonstrates a powerful approach to testing genetic drift as the default evolutionary mechanism of trait differentiation between populations.

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.

High Duty Cycle Echolocation May Constrain the Evolution of Diversity within Horseshoe Bats (Family: Rhinolophidae)

The phenotype of organisms is the net result of various evolutionary forces acting upon their lineages over time. When an innovative trait arises that confers a substantial advantage in terms of survival and reproduction, the evolution of adaptive complexes between such an innovation and other traits may constrain diversification of that lineage. The specialized echolocation system of the Rhinolophidae may represent such an innovation which affects other parts of the phenotype. We investigated a potential constraint on the diversity of phenotypes of several species of horseshoe bats within a phylogenetic framework. If phenotypic convergence stems from stasis as a result of the specialized echolocation system, phenotypes should converge not only among members of the same species and between sexes but also among species. We analyzed the phenotypic diversity of >800 individuals of 13–16 species. The phenotypes in the horseshoe bats did indeed converge. There was no sexual size dimorphism in mass, forearm length and wingspan within species and there was marked interspecific similarity in both wing and echolocation variables but marked variability in body mass. Furthermore, correlations of wing and echolocation variables with mass suggest that variability within horseshoe bats was largely the result of selection on body size with allometric responses in wing and echolocation parameters, a potential consequence of constraints imposed by their specialized echolocation.

Using Approximate Bayesian Computation to infer sex ratios from acoustic data

Population sex ratios are of high ecological relevance, but are challenging to determine in species lacking conspicuous external cues indicating their sex. Acoustic sexing is an option if vocalizations differ between sexes, but is precluded by overlapping distributions of the values of male and female vocalizations in many species. A method allowing the inference of sex ratios despite such an overlap will therefore greatly increase the information extractable from acoustic data. To meet this demand, we developed a novel approach using Approximate Bayesian Computation (ABC) to infer the sex ratio of populations from acoustic data. Additionally, parameters characterizing the male and female distribution of acoustic values (mean and standard deviation) are inferred. This information is then used to probabilistically assign a sex to a single acoustic signal. We furthermore develop a simpler means of sex ratio estimation based on the exclusion of calls from the overlap zone. Applying our methods to simulated data demonstrates that sex ratio and acoustic parameter characteristics of males and females are reliably inferred by the ABC approach. Applying both the ABC and the exclusion method to empirical datasets (echolocation calls recorded in colonies of lesser horseshoe bats, Rhinolophus hipposideros) provides similar sex ratios as molecular sexing. Our methods aim to facilitate evidence-based conservation, and to benefit scientists investigating ecological or conservation questions related to sex- or group specific behaviour across a wide range of organisms emitting acoustic signals. The developed methodology is non-invasive, low-cost and time-efficient, thus allowing the study of many sites and individuals. We provide an R-script for the easy application of the method and discuss potential future extensions and fields of applications. The script can be easily adapted to account for numerous biological systems by adjusting the type and number of groups to be distinguished (e.g. age, social rank, cryptic species) and the acoustic parameters investigated.

25 Years of sensory drive: the evidence and its watery bias

It has been 25 years since the formalization of the Sensory Drive hypothesis was published in the American Naturalist (1992). Since then, there has been an explosion of research identifying its utility in contributing to our understanding of inter- and intra-specific variation in sensory systems and signaling properties. The main tenet of Sensory Drive is that environmental characteristics will influence the evolutionary trajectory of both sensory (detecting capabilities) and signaling (detectable features and behaviors) traits in predictable directions. We review the accumulating evidence in 154 studies addressing these questions and categorized their approach in terms of testing for environmental influence on sensory tuning, signal characteristics, or both. For the subset of studies that examined sensory tuning, there was greater support for Sensory Drive processes shaping visual than auditory tuning, and it was more prevalent in aquatic than terrestrial habitats. Terrestrial habitats and visual traits were the prevalent habitat and sensory modality in the 104 studies showing support for environmental influence on signaling properties. An additional 19 studies that found no supporting evidence for environmental influence on signaling traits were all based in terrestrial ecosystems and almost exclusively involved auditory signals. Only 29 studies examined the complete coevolutionary process between sensory and signaling traits and were dominated by fish visual communication. We discuss biophysical factors that may contribute to the visual and aquatic bias for Sensory Drive evidence, as well as biotic factors that may contribute to the lack of Sensory Drive processes in terrestrial acoustic signaling systems.

Male-biased dispersal and the potential impact of human-induced habitat modifications on the Neotropical bat Trachops cirrhosus

Gene flow, maintained through natal dispersal and subsequent mating events, is one of the most important processes in both ecology and population genetics. Among mammalian populations, gene flow is strongly affected by a variety of factors, including the species' ability to disperse, and the composition of the environment which can limit dispersal. Information on dispersal patterns is thus crucial both for conservation management and for understanding the social system of a species. We used 16 polymorphic nuclear microsatellite loci in addition to mitochondrial DNA sequences (1.61 kbp) to analyse the population structure and the sex-specific pattern of natal dispersal in the frog-eating fringe-lipped bat, Trachops cirrhosus, in Central Panama. Our study revealed that-unlike most of the few other investigated Neotropical bats-gene flow in this species is mostly male-mediated. Nevertheless, distinct genetic clusters occur in both sexes. In particular, the presence of genetic differentiation in the dataset only consisting of the dispersing sex (males) indicates that gene flow is impeded within our study area. Our data are in line with the Panama Canal in connection with the widening of the Río Chagres during the canal construction acting as a recent barrier to gene flow. The sensitivity of T. cirrhosus to human-induced habitat modifications is further indicated by an extremely low capture success in highly fragmented areas. Taken together, our genetic and capture data provide evidence for this species to be classified as less mobile and thus vulnerable to habitat change, information that is important for conservation management.

Geographical variation in the echolocation calls of bent-winged bats, Miniopterus fuliginosus

Evolutionary biologists had a long-standing interest in the evolutionary forces underlying geographical variation in the acoustic signals of animals. However, the evolutionary forces driving acoustic variation are still unclear. In this study, we quantified the geographical variation in the peak frequencies of echolocation calls in eight Miniopterus fuliginosus bat colonies, and assessed the forces that drive acoustic divergence. Our results demonstrated that seven of the colonies had very similar peak frequencies, while only one colony was significantly higher than the others. This similarity in echolocation call frequency among the seven colonies was likely due to frequent dispersal and migration, leading to male-mediated infiltration of nuclear genes. This infiltration enhances gene flow and weakens ecological selection, and also increases interactions in the presence of conspecifics. Significant correlations were not observed between acoustic distances and morphological distances, climatic differences, geographic distances or mtDNA genetic distances. However, variation in acoustic distances was significantly positive correlated with nDNA genetic distance, even after controlling for geographic distance. Interestingly, the relationship between call divergence and genetic distance was no longer significant after excluding the colony with the highest call frequency, which may be due to the minimal genetic distance among the other seven colonies. The highest frequencies of echolocation calls observed in the one colony may be shaped by selection pressure due to loud background noise in the area. Taken together, these results suggest that geographic divergence of echolocation calls may not be subject to genetic drift, but rather, that the strong selective pressure induced by background noise may lead to acoustic and genetic differentiation between JXT and the other colonies.

Growing old, yet staying young: The role of telomeres in bats' exceptional longevity

Understanding aging is a grand challenge in biology. Exceptionally long-lived animals have mechanisms that underpin extreme longevity. Telomeres are protective nucleotide repeats on chromosome tips that shorten with cell division, potentially limiting life span. Bats are the longest-lived mammals for their size, but it is unknown whether their telomeres shorten. Using >60 years of cumulative mark-recapture field data, we show that telomeres shorten with age in Rhinolophus ferrumequinum and Miniopterus schreibersii, but not in the bat genus with greatest longevity, Myotis. As in humans, telomerase is not expressed in Myotis myotis blood or fibroblasts. Selection tests on telomere maintenance genes show that ATM and SETX, which repair and prevent DNA damage, potentially mediate telomere dynamics in Myotis bats. Twenty-one telomere maintenance genes are differentially expressed in Myotis, of which 14 are enriched for DNA repair, and 5 for alternative telomere-lengthening mechanisms. We demonstrate how telomeres, telomerase, and DNA repair genes have contributed to the evolution of exceptional longevity in Myotis bats, advancing our understanding of healthy aging.

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.

Stabilization of a bat-pitcher plant mutualism

Despite the long persistence of many mutualisms, it is largely unknown which mechanisms stabilize these interactions. This is especially true if only one mutualism partner can choose alternative partners while the other cannot, resulting in a power asymmetry. According to biological market theory the choosing partner should prefer the more dependent partner if the latter offers commodities of higher quality than its competitors. We tested this prediction using Bornean carnivorous pitcher plants (Nepenthes hemsleyana) that strongly rely on faecal nitrogen of bats (Kerivoula hardwickii) which roost inside the pitchers. The bats also roost in furled leaves of various plants. Surprisingly, during field observations the bats did not always choose N. hemsleyana pitchers despite their superior quality but were generally faithful either to pitchers or to furled leaves. In behavioural experiments 21% of the leaf-roosting bats switched to pitchers, while the majority of these bats and all pitcher-roosting individuals were faithful to the roost type in which we had found them. Genetic differentiation cannot explain this faithfulness, which likely results from different roosting traditions. Such traditions could have stabilizing or destabilizing effects on various mutualisms and should be investigated in more detail.

Wing Membrane Biopsies for Bat Cytogenetics: Finding of 2n = 54 in Irish Rhinolophushipposideros (Rhinolophidae, Chiroptera, Mammalia) Supports Two Geographically Separated Chromosomal Variants in Europe

In Europe, 2 different diploid chromosome numbers, 2n = 54 and 2n = 56, have been described in the lesser horseshoe bat (Rhinolophushipposideros). The eastern form with 2n = 56 extends from the Czech Republic to Greece. To date, specimens with 54 chromosomes have been reported only from Spain and Germany. This study expands the distributional area of the western variant to Ireland. Strikingly, this distribution of European chromosomal variants is in contrast to the available molecular data that indicate little genetic differentiation of R. hipposideros populations spanning Northwestern to Central Europe. Further, we have developed an optimized protocol for establishing fibroblast cell cultures, suitable for karyotype analyses, from 3-mm wing membrane biopsies. This is a useful technique for cytogenetic studies of endangered bat species, as this non-lethal sampling method imposes only minimum stress to the animal without lasting adverse effects and is routinely used to sample tissue probes for molecular genetic studies in bats.

Patterns and causes of geographic variation in bat echolocation pulses

Evolutionary biologists have a long-standing interest in how acoustic signals in animals vary geographically, because divergent ecology and sensory perception play an important role in speciation. Geographic comparisons are valuable in determining the factors that influence divergence of acoustic signals. Bats are social mammals and they depend mainly on echolocation pulses to locate prey, to navigate and to communicate. Mounting evidence shows that geographic variation of bat echolocation pulses is common, with a mean 5-10 kHz differences in peak frequency, and a high level of individual variation may be nested in this geographical variation. However, understanding the geographic variation of echolocation pulses in bats is very difficult, because of differences in sample and statistical analysis techniques as well as the variety of factors shaping the vocal geographic evolution. Geographic differences in echolocation pulses of bats generally lack latitudinal, longitudinal and elevational patterns, and little is known about vocal dialects. Evidence is accumulating to support the fact that geographic variation in echolocation pulses of bats may be caused by genetic drift, cultural drift, ecological selection, sexual selection and social selection. Future studies could relate geographic differences in echolocation pulses to social adaptation, vocal learning strategies and patterns of dispersal. In addition, new statistical techniques and acoustic playback experiments may help to illustrate the causes and consequences of the geographic evolution of echolocation pulse in bats.

Nuclear introns outperform mitochondrial DNA in inter-specific phylogenetic reconstruction: Lessons from horseshoe bats (Rhinolophidae: Chiroptera)

Despite many studies illustrating the perils of utilising mitochondrial DNA in phylogenetic studies, it remains one of the most widely used genetic markers for this purpose. Over the last decade, nuclear introns have been proposed as alternative markers for phylogenetic reconstruction. However, the resolution capabilities of mtDNA and nuclear introns have rarely been quantified and compared. In the current study we generated a novel ∼5kb dataset comprising six nuclear introns and a mtDNA fragment. We assessed the relative resolution capabilities of the six intronic fragments with respect to each other, when used in various combinations together, and when compared to the traditionally used mtDNA. We focused on a major clade in the horseshoe bat family (Afro-Palaearctic clade; Rhinolophidae) as our case study. This old, widely distributed and speciose group contains a high level of conserved morphology. This morphological stasis renders the reconstruction of the phylogeny of this group with traditional morphological characters complex. We sampled multiple individuals per species to represent their geographic distributions as best as possible (122 individuals, 24 species, 68 localities). We reconstructed the species phylogeny using several complementary methods (partitioned Maximum Likelihood and Bayesian and Bayesian multispecies-coalescent) and made inferences based on consensus across these methods. We computed pairwise comparisons based on Robinson-Foulds tree distance metric between all Bayesian topologies generated (27,000) for every gene(s) and visualised the tree space using multidimensional scaling (MDS) plots. Using our supported species phylogeny we estimated the ancestral state of key traits of interest within this group, e.g. echolocation peak frequency which has been implicated in speciation. Our results revealed many potential cryptic species within this group, even in taxa where this was not suspected a priori and also found evidence for mtDNA introgression. We demonstrated that by using just two introns one can recover a better supported species tree than when using the mtDNA alone, despite the shorter overall length of the combined introns. Additionally, when combining any single intron with mtDNA, we showed that the result is highly similar to the mtDNA gene tree and far from the true species tree and therefore this approach should be avoided. We caution against the indiscriminate use of mtDNA in phylogenetic studies and advocate for pilot studies to select nuclear introns. The selection of marker type and number is a crucial step that is best based on critical examination of preliminary or previously published data. Based on our findings and previous publications, we recommend the following markers to recover phylogenetic relationships between recently diverged taxa (<20 My) in bats and other mammals: ACOX2, COPS7A, BGN, ROGDI and STAT5A.

A nonlethal sampling method to obtain, generate and assemble whole blood transcriptomes from small, wild mammals

The acquisition of tissue samples from wild populations is a constant challenge in conservation biology, especially for endangered species and protected species where nonlethal sampling is the only option. Whole blood has been suggested as a nonlethal sample type that contains a high percentage of bodywide and genomewide transcripts and therefore can be used to assess the transcriptional status of an individual, and to infer a high percentage of the genome. However, only limited quantities of blood can be nonlethally sampled from small species and it is not known if enough genetic material is contained in only a few drops of blood, which represents the upper limit of sample collection for some small species. In this study, we developed a nonlethal sampling method, the laboratory protocols and a bioinformatic pipeline to sequence and assemble the whole blood transcriptome, using Illumina RNA-Seq, from wild greater mouse-eared bats (Myotis myotis). For optimal results, both ribosomal and globin RNAs must be removed before library construction. Treatment of DNase is recommended but not required enabling the use of smaller amounts of starting RNA. A large proportion of protein-coding genes (61%) in the genome were expressed in the blood transcriptome, comparable to brain (65%), kidney (63%) and liver (58%) transcriptomes, and up to 99% of the mitogenome (excluding D-loop) was recovered in the RNA-Seq data. In conclusion, this nonlethal blood sampling method provides an opportunity for a genomewide transcriptomic study of small, endangered or critically protected species, without sacrificing any individuals.