Geographic variation in the acoustic traits of greater horseshoe bats: testing the importance of drift and ecological selection in evolutionary processes

Adaptive Evolution of the Greater Horseshoe Bat AANAT: Insights into the Link between AANAT and Hibernation Rhythms

Arylalkylamine N-acetyltransferase (AANAT) is a crucial rate-limiting enzyme in the synthesis of melatonin. AANAT has been confirmed to be independently duplicated and inactivated in different animal taxa in order to adapt to the environment. However, the evolutionary forces associated with having a single copy of AANAT remain unclear. The greater horseshoe bat has a single copy of AANAT but exhibits different hibernation rhythms in various populations. We analyzed the adaptive evolution at the gene and protein levels of AANAT from three distinct genetic lineages in China: northeast (NE), central east (CE), and southwest (SW). The results revealed greater genetic diversity in the AANAT loci of the NE and CE lineage populations that have longer hibernation times, and there were two positive selection loci. The catalytic capacity of AANAT in the Liaoning population that underwent positive selection was significantly higher than that of the Yunnan population (p < 0.05). This difference may be related to the lower proportion of α helix and the variation in two interface residues. The adaptive evolution of AANAT was significantly correlated with climate and environment (p < 0.05). After controlling for geographical factors (latitude and altitude), the evolution of AANAT by the negative temperature factor was represented by the monthly mean temperature (r = -0.6, p < 0.05). The results identified the gene level variation, functional adaptation, and evolutionary driving factors of AANAT, provide an important foundation for further understanding the adaptive evolution of the single copy of AANAT in pteropods, and may offer evidence for adaptive hibernation rhythms in bats.

Sex Differences in Cochlear Transcriptomes in Horseshoe Bats

Sexual dimorphism of calls is common in animals, whereas studies on the molecular basis underlying this phenotypic variation are still scarce. In this study, we used comparative transcriptomics of cochlea to investigate the sex-related difference in gene expression and alternative splicing in four Rhinolophus taxa. Based on 31 cochlear transcriptomes, we performed differential gene expression (DGE) and alternative splicing (AS) analyses between the sexes in each taxon. Consistent with the degree of difference in the echolocation pulse frequency between the sexes across the four taxa, we identified the largest number of differentially expressed genes (DEGs) and alternatively spliced genes (ASGs) in R. sinicus. However, we also detected multiple DEGs and ASGs in taxa without sexual differences in echolocation pulse frequency, suggesting that these genes might be related to other parameters of echolocation pulse rather than the frequency component. Some DEGs and ASGs are related to hearing loss or deafness genes in human or mice and they can be considered to be candidates associated with the sexual differences of echolocation pulse in bats. We also detected more than the expected overlap of DEGs and ASGs in two taxa. Overall, our current study supports the important roles of both DGE and AS in generating or maintaining sexual differences in animals.

Unraveling the macroevolution of horseshoe bats (Chiroptera: Rhinolophidae: Rhinolophus)

Unraveling the diversification mechanisms of organisms is a fundamental and important macroevolutionary question regarding the diversity, ecological niche, and morphological divergence of life. However, many studies have only explored diversification mechanisms via isolated factors. Here, based on comparative phylogenetic analysis, we performed a macroevolutionary examination of horseshoe bats (Chiroptera: Rhinolophidae: Rhinolophus), to reveal the inter-relationships among diversification, intrinsic/extrinsic factors, and climatic ecological niche characteristics. Results showed a general slowing trajectory during diversification, with two dispersal events from Asia into Southeast Asia and Africa playing key roles in shaping regional heterogeneous diversity. Morphospace expansions of the investigated traits (e.g., body size, echolocation, and climate niche) revealed a decoupled pattern between diversification trajectory and trait divergence, suggesting that other factors (e.g., biotic interactions) potentially played a key role in recent diversification. Based on ancestral traits and pathway analyses, most Rhinolophus lineages belonging to the same region overlapped with each other geographically and were positively associated with the diversification rate, implying a competitive prelude to speciation. Overall, our study showed that multiple approaches need to be integrated to address diversification history. Rather than a single factor, the joint effects of multiple factors (biogeography, environmental drivers, and competition) are responsible for the current diversity patterns in horseshoe bats, and a corresponding multifaceted strategy is recommended to study these patterns in the future.

Diversifying selection and climatic effects on major histocompatibility complex class II gene diversity in the greater horseshoe bat

Heterogeneous pathogenic stress can shape major histocompatibility complex (MHC) diversity by influencing the functional plasticity of the immune response. Therefore, MHC diversity could reflect environmental stress, demonstrating its importance in uncovering the mechanisms of adaptive genetic variation. In this study, we combined neutral microsatellite loci, an immune-related MHC II-DRB locus, and climatic factors to unravel the mechanisms affecting the diversity and genetic differentiation of MHC genes in the greater horseshoe bat (Rhinolophus ferrumequinum), a species with a wide geographical distribution that has three distinct genetic lineages in China. First, increased genetic differentiation at the MHC locus among populations compared using microsatellites indicated diversifying selection. Second, the genetic differentiation of MHC and microsatellites were significantly correlated, suggesting that demographic processes exist. However, MHC genetic differentiation was significantly correlated with geographical distance among populations, even after controlling for the neutral markers, suggesting a major effect of selection. Third, although the MHC genetic differentiation was larger than that for microsatellites, there was no significant difference in the genetic differentiation between the two markers among genetic lineages, indicating the effect of balancing selection. Fourth, combined with climatic factors, MHC diversity and supertypes showed significant correlations with temperature and precipitation, but not with the phylogeographic structure of R. ferrumequinum, suggesting an effect of local adaptation driven by climate on MHC diversity. Moreover, the number of MHC supertypes varied between populations and lineages, suggesting regional characteristics and support for local adaptation. Taken together, the results of our study provide insights into the adaptive evolutionary driving forces at different geographic scales in R. ferrumequinum. In addition, climate factors may have played a vital role in driving adaptive evolution in this species.

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.

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

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

Correlation between acoustic divergence and phylogenetic distance in soniferous European gobiids (Gobiidae; Gobius lineage)

In fish, species identity can be encoded by sounds, which have been thoroughly investigated in European gobiids (Gobiidae, Gobius lineage). Recent evolutionary studies suggest that deterministic and/or stochastic forces could generate acoustic differences among related animal species, though this has not been investigated in any teleost group to date. In the present comparative study, we analysed the sounds from nine soniferous gobiids and quantitatively assessed their acoustic variability. Our interspecific acoustic study, incorporating for the first time the representative acoustic signals from the majority of soniferous gobiids, suggested that their sounds are truly species-specific (92% of sounds correctly classified into exact species) and each taxon possesses a unique set of spectro-temporal variables. In addition, we reconstructed phylogenetic relationships from a concatenated molecular dataset consisting of multiple molecular markers to track the evolution of acoustic signals in soniferous gobiids. The results of this study indicated that the genus Padogobius is polyphyletic, since P. nigricans was nested within the Ponto-Caspian clade, while the congeneric P. bonelli turned out to be a sister taxon to the remaining investigated soniferous species. Lastly, by extracting the acoustic and genetic distance matrices, sound variability and genetic distance were correlated for the first time to assess whether sound evolution follows a similar phylogenetic pattern. The positive correlation between the sound variability and genetic distance obtained here emphasizes that certain acoustic features from representative sounds could carry the phylogenetic signal in soniferous gobiids. Our study was the first attempt to evaluate the mutual relationship between acoustic variation and genetic divergence in any teleost fish.

Geographic variation in the skulls of the horseshoe bats, Rhinolophus simulator and R. cf. simulator: Determining the relative contributions of adaptation and drift using geometric morphometrics

The relative contributions of adaptation and genetic drift to morphological diversification of the skulls of echolocating mammals were investigated using two horseshoe bat species, Rhinolophus simulator and R. cf. simulator, as test cases. We used 3D geometric morphometrics to compare the shapes of skulls of the two lineages collected at various localities in southern Africa. Size and shape variation was predominantly attributed to selective forces; the between-population variance (B) was not proportional to the within-population variance (W). Modularity was evident in the crania of R. simulator but absent in the crania of R. cf. simulator and the mandibles of both species. The skulls of the two lineages thus appeared to be under different selection pressures, despite the overlap in their distributions. Difference in the crania of R. cf. simulator was centered largely on the nasal dome region of R. cf. simulator but on the cranium and mandibles of R. simulator. It is likely that the size and shape of the nasal dome, which acts as a frequency-dependent acoustic horn, is more crucial in R. cf. simulator than in R. simulator because of the higher echolocation frequencies used by R. cf. simulator. A larger nasal dome in R. cf. simulator would allow the emission of higher intensity pulses, resulting in comparable detection distances to that of R. simulator. In contrast, selection pressure is probably more pronounced on the mandibles and cranium of R. simulator to compensate for the loss in bite force because of its elongated rostrum. The predominance of selection probably reflects the stringent association between environment and the optimal functioning of phenotypic characters associated with echolocation and feeding in bats.

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

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

Species specificity and intraspecific variation in the chemical profiles of Heliconius butterflies across a large geographic range

In many animals, mate choice is important for the maintenance of reproductive isolation between species. Traits important for mate choice and behavioral isolation are predicted to be under strong stabilizing selection within species; however, such traits can also exhibit variation at the population level driven by neutral and adaptive evolutionary processes. Here, we describe patterns of divergence among androconial and genital chemical profiles at inter- and intraspecific levels in mimetic Heliconius butterflies. Most variation in chemical bouquets was found between species, but there were also quantitative differences at the population level. We found a strong correlation between interspecific chemical and genetic divergence, but this correlation varied in intraspecific comparisons. We identified "indicator" compounds characteristic of particular species that included compounds already known to elicit a behavioral response, suggesting an approach for identification of candidate compounds for future behavioral studies in novel systems. Overall, the strong signal of species identity suggests a role for these compounds in species recognition, but with additional potentially neutral variation at the population level.

Evolutionary significance of the variation in acoustic communication of a cryptic nocturnal primate radiation (Microcebus spp.)

Acoustic phenotypic variation is of major importance for speciation and the evolution of species diversity. Whereas selective and stochastic forces shaping the acoustic divergence of signaling systems are well studied in insects, frogs, and birds, knowledge on the processes driving acoustic phenotypic evolution in mammals is limited. We quantified the acoustic variation of a call type exchanged during agonistic encounters across eight distinct species of the smallest-bodied nocturnal primate radiation, the Malagasy mouse lemurs. The species live in two different habitats (dry forest vs. humid forest), differ in geographic distance to each other, and belong to four distinct phylogenetic clades within the genus. Genetically defined species were discriminated reliably on the phenotypic level based on their acoustic distinctiveness in a discriminant function analysis. Acoustic variation was explained by genetic distance, whereas differences in morphology, forest type, or geographic distance had no effect. The strong impact of genetics was supported by a correlation between acoustic and genetic distance and the high agreement in branching pattern between the acoustic and molecular phylogenetic trees. In sum, stochastic factors such as genetic drift best explained acoustic diversification in a social communication call of mouse lemurs.

Whistle variation in Mediterranean common bottlenose dolphin: The role of geographical, anthropogenic, social, and behavioral factors

The studies on the variation of acoustic communication in different species have provided insight that genetics, geographic isolation, and adaptation to ecological and social conditions play important roles in the variability of acoustic signals. The dolphin whistles are communication signals that can vary significantly among and within populations. Although it is known that they are influenced by different environmental and social variables, the factors influencing the variation between populations have received scant attention. In the present study, we investigated the factors associated with the acoustic variability in the whistles of common bottlenose dolphin (Tursiops truncatus), inhabiting two Mediterranean areas (Sardinia and Croatia). We explored which factors, among (a) geographical isolation of populations, (b) different environments in terms of noise and boat presence, and (c) social factors (including group size, behavior, and presence of calves), were associated with whistle characteristics. We first applied a principal component analysis to reduce the number of collinear whistle frequency and temporal characteristics and then generalized linear mixed models on the first two principal components. The study revealed that both geographic distance/isolation and local environment are associated with whistle variations between localities. The prominent differences in the acoustic environments between the two areas, which contributed to the acoustic variability in the first principal component (PC1), were found. The calf's presence and foraging and social behavior were also found to be associated with dolphin whistle variation. The second principal component (PC2) was associated only with locality and group size, showing that longer and more complex tonal sound may facilitate individual recognition and cohesion in social groups. Thus, both social and behavioral context influenced significantly the structure of whistles, and they should be considered when investigating acoustic variability among distant dolphin populations to avoid confounding factors.

Behaviour, biology and evolution of vocal learning in bats

The comparative approach can provide insight into the evolution of human speech, language and social communication by studying relevant traits in animal systems. Bats are emerging as a model system with great potential to shed light on these processes given their learned vocalizations, close social interactions, and mammalian brains and physiology. A recent framework outlined the multiple levels of investigation needed to understand vocal learning across a broad range of non-human species, including cetaceans, pinnipeds, elephants, birds and bats. Here, we apply this framework to the current state-of-the-art in bat research. This encompasses our understanding of the abilities bats have displayed for vocal learning, what is known about the timing and social structure needed for such learning, and current knowledge about the prevalence of the trait across the order. It also addresses the biology (vocal tract morphology, neurobiology and genetics) and evolution of this trait. We conclude by highlighting some key questions that should be answered to advance our understanding of the biological encoding and evolution of speech and spoken communication. This article is part of the theme issue 'What can animal communication teach us about human language?'

Characterization of the acoustic community of vocal fishes in the Azores

Sounds produced by teleost fishes are an important component of marine soundscapes, making passive acoustic monitoring (PAM) an effective way to map the presence of vocal fishes with a minimal impact on ecosystems. Based on a literature review, we list the known soniferous fish species occurring in Azorean waters and compile their sounds. We also describe new fish sounds recorded in Azores seamounts. From the literature, we identified 20 vocal fish species present in Azores. We analysed long-term acoustic recordings carried out since 2008 in Condor and Princesa Alice seamounts and describe 20 new putative fish sound sequences. Although we propose candidates as the source of some vocalizations, this study puts into evidence the myriad of fish sounds lacking species identification. In addition to identifying new sound sequences, we provide the first marine fish sound library for Azores. Our acoustic library will allow to monitor soniferous fish species for conservation and management purposes.

Gene expression vs. sequence divergence: comparative transcriptome sequencing among natural Rhinolophus ferrumequinum populations with different acoustic phenotypes

BACKGROUND

Although the sensory drive hypothesis can explain the geographic variation in echolocation frequencies of some bat species, the molecular mechanisms underlying this phenomenon are still unclear. The three lineages of greater horseshoe bat (Rhinolophus ferrumequinum) in China (northeast, central-east, and southwest) have significant geographic variation in resting frequencies (RF) of echolocation calls. Because their cochleae have an acoustic fovea that is highly sensitive to a narrow range of frequencies, we reported the transcriptomes of cochleae collected from three genetic lineages of R. ferrumequinum, which is an ideal organism for studying geographic variation in echolocation signals, and tried to understand the mechanisms behind this bat phenomenon by analyzing gene expression and sequence variation.

RESULTS

A total of 8190 differentially expressed genes (DEGs) were identified. We identified five modules from all DEGs that were significantly related to RF or forearm length (FL). DEGs in the RF-related modules were significantly enriched in the gene categories involved in neural activity, learning, and response to sound. DEGs in the FL-related modules were significantly enriched in the pathways related to muscle and actin functions. Using 21,945 single nucleotide polymorphisms, we identified 18 candidate unigenes associated with hearing, five of which were differentially expressed among the three populations. Additionally, the gene ERBB4, which regulates diverse cellular processes in the inner ear such as cell proliferation and differentiation, was in the largest module. We also found 49 unigenes that were under positive selection from 4105 one-to-one orthologous gene pairs between the three R. ferrumequinum lineages and three other Chiroptera species.

CONCLUSIONS

The variability of gene expression and sequence divergence at the molecular level might provide evidence that can help elucidate the genetic basis of geographic variation in echolocation signals of greater horseshoe bats.

Understanding sources of variance and correlation among features of Bornean gibbon (Hylobates muelleri) female calls

Acoustic signals serve important functions in mate choice, resource defense, and species recognition. Quantifying patterns and sources of variation in acoustic signals can advance understanding of the evolutionary processes that shape behavioral diversity more broadly. Animal vocalization datasets are inherently multivariate and hierarchical, wherein multiple features are estimated from calls of many individuals across different recording locations. Patterns of variation within different hierarchical levels-notwithstanding the challenges they present for modeling and inference-can provide insight into processes shaping vocal variation. The current work presents a multivariate, variance components model to investigate three levels of variance (within-female, between-female, and between-site) in Bornean gibbon calls. For six of the eight features estimated from call spectrograms, between-female variance was the most important contributor to total variance. For one feature, trill rate, there were site-level differences, which may be related to geographic isolation of certain gibbon populations. There was also a negative relationship between trill rate and duration of the introduction, suggesting trade-offs in the production of gibbon calls. Given substantial inter-individual variation in gibbon calls, it seems likely that there has been selection to confer information regarding caller identity, but mechanisms leading to site-level variation in trill rate remain to be determined.

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.

Testing the Sensory Drive Hypothesis: Geographic variation in echolocation frequencies of Geoffroy's horseshoe bat (Rhinolophidae: Rhinolophus clivosus)

Geographic variation in sensory traits is usually influenced by adaptive processes because these traits are involved in crucial life-history aspects including orientation, communication, lineage recognition and mate choice. Studying this variation can therefore provide insights into lineage diversification. According to the Sensory Drive Hypothesis, lineage diversification may be driven by adaptation of sensory systems to local environments. It predicts that acoustic signals vary in association with local climatic conditions so that atmospheric attenuation is minimized and transmission of the signals maximized. To test this prediction, we investigated the influence of climatic factors (specifically relative humidity and temperature) on geographic variation in the resting frequencies of the echolocation pulses of Geoffroy's horseshoe bat, Rhinolophus clivosus. If the evolution of phenotypic variation in this lineage tracks climate variation, human induced climate change may lead to decreases in detection volumes and a reduction in foraging efficiency. A complex non-linear interaction between relative humidity and temperature affects atmospheric attenuation of sound and principal components composed of these correlated variables were, therefore, used in a linear mixed effects model to assess their contribution to observed variation in resting frequencies. A principal component composed predominantly of mean annual temperature (factor loading of -0.8455) significantly explained a proportion of the variation in resting frequency across sites (P < 0.05). Specifically, at higher relative humidity (around 60%) prevalent across the distribution of R. clivosus, increasing temperature had a strong negative effect on resting frequency. Climatic factors thus strongly influence acoustic signal divergence in this lineage, supporting the prediction of the Sensory Drive Hypothesis. The predicted future increase in temperature due to climate change is likely to decrease the detection volume in echolocating bats and adversely impact their foraging efficiency.

Spatial and temporal trends in fin whale vocalizations recorded in the NE Pacific Ocean between 2003-2013

In order to study the long-term stability of fin whale (Balaenoptera physalus) singing behavior, the frequency and inter-pulse interval of fin whale 20 Hz vocalizations were observed over 10 years from 2003–2013 from bottom mounted hydrophones and seismometers in the northeast Pacific Ocean. The instrument locations extended from 40°N to 48°N and 130°W to 125°W with water depths ranging from 1500–4000 m. The inter-pulse interval (IPI) of fin whale song sequences was observed to increase at a rate of 0.54 seconds/year over the decade of observation. During the same time period, peak frequency decreased at a rate of 0.17 Hz/year. Two primary call patterns were observed. During the earlier years, the more commonly observed pattern had a single frequency and single IPI. In later years, a doublet pattern emerged, with two dominant frequencies and IPIs. Many call sequences in the intervening years appeared to represent a transitional state between the two patterns. The overall trend was consistent across the entire geographical span, although some regional differences exist. Understanding changes in acoustic behavior over long time periods is needed to help establish whether acoustic characteristics can be used to help determine population identity in a widely distributed, difficult to study species such as the fin whale.

Territorial choruses of giant otter groups (Pteronura brasiliensis) encode information on group identity

Group living animals often engage in corporate territorial defence. Territorial group vocalizations can provide information about group identity, size and composition. Neighbouring groups may use this information to avoid unfavourable direct conflicts. Giant otters are highly social and territorial animals with an elaborate vocal repertoire. They produce long-range screams when they are alert or excited, i.e. in an alarm, isolation or begging context. Long-range screams are not only produced by one individual at a time ('single screams') but also by multiple group members simultaneously, resulting in a highly conspicuous 'group chorus'. Wild giant otters regularly produce group choruses during interactions with predators, when they detect intruders in their territory or before group reunions after separation. Since single screams and especially group choruses probably contribute to the groups' corporate territorial defence, we hypothesized that group identity is encoded in single screams and group choruses. We analysed vocalizations from five wild and three captive giant otter groups and found statistical evidence for a group signature in group choruses. Results for single screams were less conclusive, which might have been caused by the comparatively lower sample size. We suggest that giant otters may gain information on group identity by listening to group choruses. Group identity likely constitutes important social information for giant otters since territory boundaries of neighbouring groups can overlap and direct inter-group conflicts are severe. Therefore, group chorusing may contribute to the mutual avoidance of members from different groups.

The relative contribution of drift and selection to phenotypic divergence: A test case using the horseshoe bats Rhinolophus simulator and Rhinolophus swinnyi

Natural selection and drift can act on populations individually, simultaneously or in tandem and our understanding of phenotypic divergence depends on our ability to recognize the contribution of each. According to the quantitative theory of evolution, if an organism has diversified through neutral evolutionary processes (mutation and drift), variation of phenotypic characteristics between different geographic localities (B) should be directly proportional to the variation within localities (W), that is, B ∝ W. Significant deviations from this null model imply that non-neutral forces such as natural selection are acting on a phenotype. We investigated the relative contributions of drift and selection to intraspecific diversity using southern African horseshoe bats as a test case. We characterized phenotypic diversity across the distributional range of Rhinolophus simulator (n = 101) and Rhinolophus swinnyi (n = 125) using several traits associated with flight and echolocation. Our results suggest that geographic variation in both species was predominantly caused by disruptive natural selection (B was not directly proportional to W). Evidence for correlated selection (co-selection) among traits further confirmed that our results were not compatible with drift. Selection rather than drift is likely the predominant evolutionary process shaping intraspecific variation in traits that strongly impact fitness.

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.

Phylogenetic relationships and evolutionary history of the greater horseshoe bat, Rhinolophus ferrumequinum, in Northeast Asia

The greater horseshoe bat, Rhinolophus ferrumequinum, is an important model organism for studies on chiropteran phylogeographic patterns. Previous studies revealed the population history of R. ferrumequinum from Europe and most Asian regions, yet there continue to be arguments about their evolutionary process in Northeast Asia. In this study, we obtained mitochondrial DNA cyt b and D-loop data of R. ferrumequinum from Northeast China, South Korea and Japan to clarify their phylogenetic relationships and evolutionary process. Our results indicate a highly supported monophyletic group of Northeast Asian greater horseshoe bats, in which Japanese populations formed a single clade and clustered into the mixed branches of Northeast Chinese and South Korean populations. We infer that R. ferrumequinum in Northeast Asia originated in Northeast China and South Korea during a cold glacial period, while some ancestors likely arrived in Japan by flying or land bridge and subsequently adapted to the local environment. Consequently, during the warm Eemian interglaciation, the Korea Strait, between Japan and South Korea, became a geographical barrier to Japanese and inland populations, while the Changbai Mountains, between China and North Korea, did not play a significant role as a barrier between Northeast China and South Korea populations.

Geographic variation in advertisement calls of a Microhylid frog - testing the role of drift and ecology

Acoustic signals for mating are important traits that could drive population differentiation and speciation. Ecology may play a role in acoustic divergence through direct selection (e.g., local adaptation to abiotic environment), constraint of correlated traits (e.g., acoustic traits linked to another trait under selection), and/or interspecific competition (e.g., character displacement). However, genetic drift alone can also drive acoustic divergence. It is not always easy to differentiate the role of ecology versus drift in acoustic divergence. In this study, we tested the role of ecology and drift in shaping geographic variation in the advertisement calls of Microhyla fissipes. We examined three predictions based on ecological processes: (1) the correlation between temperature and call properties across M. fissipes populations; (2) the correlation between call properties and body size across M. fissipes populations; and (3) reproductive character displacement (RCD) in call properties between M. fissipes populations that are sympatric with and allopatric to a congener M. heymonsi. To test genetic drift, we examined correlations among call divergence, geographic distance, and genetic distance across M. fissipes populations. We recorded the advertisement calls from 11 populations of M. fissipes in Taiwan, five of which are sympatrically distributed with M. heymonsi. We found geographic variation in both temporal and spectral properties of the advertisement calls of M. fissipes. However, the call properties were not correlated with local temperature or the callers' body size. Furthermore, we did not detect RCD. By contrast, call divergence, geographic distance, and genetic distance between M. fissipes populations were all positively correlated. The comparisons between phenotypic Q_st (P_st) and F_st values did not show significant differences, suggesting a role of drift. We concluded that genetic drift, rather than ecological processes, is the more likely driver for the geographic variation in the advertisement calls of M. fissipes.

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.

Gene expression and adaptive evolution of ZBED1 in the hibernating greater horseshoe bat (Rhinolophus ferrumequinum)

Mammalian hibernators experience physiological extremes, e.g. ischemia, muscle disuse and hypothermia, which are lethal to non-hibernators, implying the existence of underlying mechanisms that allow hibernators to withstand these physiological extremes. Increased cell proliferation is suggested to be such a strategy, but its molecular basis remains unknown. In this study, we characterized the expression pattern of ZBED1 (zinc finger, BED-type containing 1), a transcription factor that plays a crucial role in regulating cell proliferation, in five tissues of the greater horseshoe bat (Rhinolophus ferrumequinum) during pre-hibernation, deep hibernation and post-hibernation. Moreover, we investigated the ZBED1 genetic divergence from individuals with variable hibernation phenotypes that cover all three known mtDNA lineages of the species. Expression analyses showed that ZBED1 is overexpressed only in brain and skeletal muscle, not in the other three tissues, suggesting an increased cell proliferation in these two tissues during deep hibernation. Evolutionary analyses showed that ZBED1 sequences were clustered into two well-supported clades with each one dominated by hibernating and non-hibernating individuals, respectively. Positive selection analyses further showed some positively selected sites and a divergent selection pressure among hibernating and non-hibernating groups of R. ferrumequinum. Our results suggest that ZBED1 as a potential candidate gene that regulates cell proliferation for hibernators to face physiological extremes during hibernation.

Natural epigenetic variation in bats and its role in evolution

When facing the challenges of environmental change, such as habitat fragmentation, organisms have to adjust their phenotype to adapt to various environmental stresses. Recent studies show that epigenetic modifications could mediate environmentally induced phenotypic variation, and this epigenetic variance could be inherited by future generations, indicating that epigenetic processes have potential evolutionary effects. Bats living in diverse environments show geographic variations in phenotype, and the females usually have natal philopatry, presenting an opportunity to explore how environments shape epigenetic marks on the genome and the evolutionary potential of epigenetic variance in bat populations for adaptation. We have explored the natural epigenetic diversity and structure of female populations of the great roundleaf bat (Hipposideros armiger), the least horseshoe bat (Rhinolophus pusillus) and the eastern bent-winged bat (Miniopterus fuliginosus) using a methylation-sensitive amplified polymorphism technique. We have also estimated the effects of genetic variance and ecological variables on epigenetic diversification. All three bat species have a low level of genomic DNA methylation and extensive epigenetic diversity that exceeds the corresponding genetic variance. DNA sequence divergence, epigenetic drift and environmental variables contribute to the epigenetic diversities of each species. Environment-induced epigenetic variation may be inherited as a result of both mitosis and meiosis, and their potential roles in evolution for bat populations are also discussed in this review.