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

Detection distances in desert dwelling, high duty cycle echolocators: A test of the foraging habitat hypothesis

High Duty Cycle (HDC) echolocating bats use high frequency echolocation pulses that are clutter resistant, but their high frequencies give them limited range. Despite their unique ability to reject background clutter while simultaneously detecting fluttering prey, the frequency of their echolocation pulses has a strong correlation with level of environmental clutter, lower frequency pulses of HDC bats being associated with more open environments. The Foraging Habitat Hypothesis (FHH) proposes that the ecological significance of these lower frequency pulses in HDC bats in open environments is that they allow longer prey detection distances. To test the FHH, we compared the frequencies, Source Levels (SLs) and detection distances of Rhinolophus capensis, a HDC bat that has been shown to vary its call frequency in relation to habitat structure. As a further test of the FHH we investigated the SLs and detection distances of Rhinolophus damarensis (a heterospecific species that occurs in the same open desert environment as R. capensis but echolocates at a higher dominant pulse frequency). In the open desert, R. capensis emitted both lower frequency and higher SL pulses giving them longer detection distances than R. capensis in the cluttered fynbos. SL contributed more to differences in detection distances in both R. capensis and R. damarensis than frequency. In a few instances, R. damarensis achieved similar detection distances to desert–inhabiting R. capensis by emitting much higher SLs despite their average SLs being lower. These results suggest that lower frequency echolocation pulses are not a prerequisite for open desert living but may increase detection distance while avoiding energetic costs required for high SLs.

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.

The role of common ancestry and gene flow in the evolution of human-directed play behaviour in dogs

Among-population variance of phenotypic traits is of high relevance for understanding evolutionary mechanisms that operate in relatively short timescales, but various sources of nonindependence, such as common ancestry and gene flow, can hamper the interpretations. In this comparative analysis of 138 dog breeds, we demonstrate how such confounders can independently shape the evolution of a behavioural trait (human-directed play behaviour from the Dog Mentality Assessment project). We combined information on genetic relatedness and haplotype sharing to reflect common ancestry and gene flow, respectively, and entered these into a phylogenetic mixed model to partition the among-breed variance of human-directed play behaviour while also accounting for within-breed variance. We found that 75% of the among-breed variance was explained by overall genetic relatedness among breeds, whereas 15% could be attributed to haplotype sharing that arises from gene flow. Therefore, most of the differences in human-directed play behaviour among breeds have likely been caused by constraints of common ancestry as a likely consequence of past selection regimes. On the other hand, gene flow caused by crosses among breeds has played a minor, but not negligible role. Our study serves as an example of an analytical approach that can be applied to comparative situations where the effects of shared origin and gene flow require quantification and appropriate statistical control in a within-species/among-population framework. Altogether, our results suggest that the evolutionary history of dog breeds has left remarkable signatures on the among-breed variation of a behavioural phenotype.

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.

It's not all about the Soprano: Rhinolophid bats use multiple acoustic components in echolocation pulses to discriminate between conspecifics and heterospecifics

Acoustic communication plays a pivotal role in conspecific recognition in numerous animal taxa. Vocalizations must therefore have discrete acoustic signatures to facilitate intra-specific communication and to avoid misidentification. Here we investigate the potential role of echolocation in communication in horseshoe bats. Although it has been demonstrated that echolocation can be used to discriminate among con- and hetero-specifics, the specific acoustic cues used in discrimination are still relatively unknown. Furthermore, the Acoustic Communication Hypothesis proposes that in multispecies assemblages, in which echolocation frequencies are likely to overlap, bats partition acoustic space along several dimensions so that each species occupies a discrete communication domain. Thus, multiple echolocation variables may be used in discrimination. The objective of this study was to investigate the potential of various echolocation variables to function as discriminatory cues in echolocation-based species discrimination. Using habituation-dishabituation playback experiments, we firstly tested the ability of Rhinolophus clivosus to discriminate between echolocation pulses of heterospecifics with either discrete or overlapping frequencies. Secondly, to determine whether R. clivosus could use echolocation variables other than frequency, we investigated its ability to discriminate among echolocation pulses differing in only one manipulated parameter. These test variables were identified by their contribution to the dissimilarity among pulses. Our results suggest that R. clivosus could discriminate readily between species using echolocation pulses with discrete frequencies. When frequencies overlapped, the ability of bats to discriminate was dependant on additional acoustic variables that defined the acoustic space occupied by the test signal. These additional acoustic variables included, but may not be restricted to, sweep rate of the FM and duty cycle. Thus, when echolocation pulses share a similar acoustic domain, bats use several cues to reliably discriminate among heterospecifics.

Bats are still not birds in the digital era: echolocation call variation and why it matters for bat species identification

The recording and analysis of echolocation calls are fundamental methods used to study bat distribution, ecology, and behavior. However, the goal of identifying bats in flight from their echolocation calls is not always possible. Unlike bird songs, bat calls show large variation that often makes identification challenging. The problem has not been fully overcome by modern digital-based hardware and software for bat call recording and analysis. Besides providing fundamental insights into bat physiology, ecology, and behavior, a better understanding of call variation is therefore crucial to best recognize limits and perspectives of call classification. We provide a comprehensive overview of sources of interspecific and intraspecific echolocation call variations, illustrating its adaptive significance and highlighting gaps in knowledge. We remark that further research is needed to better comprehend call variation and control for it more effectively in sound analysis. Despite the state-of-art technology in this field, combining acoustic surveys with capture and roost search, as well as limiting identification to species with distinctive calls, still represent the safest way of conducting bat surveys.

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.

Environmental correlates of geographic divergence in a phenotypic trait: A case study using bat echolocation

Divergence in phenotypic traits may arise from the interaction of different evolutionary forces, including different kinds of selection (e.g., ecological), genetic drift, and phenotypic plasticity. Sensory systems play an important role in survival and reproduction, and divergent selection on such systems may result in lineage diversification. Such diversification could be largely influenced by selection in different environments as a result of isolation by environment (IbE). We investigated this process using geographic variation in the resting echolocation frequency of the horseshoe bat species, Rhinolophus damarensis, as a test case. Bats were sampled along a latitudinal gradient ranging from 16°S to 32°S in the arid western half of southern Africa. We measured body size and peak resting frequencies (RF) from handheld individual bats. Three hypotheses for the divergence in RF were tested: (1) James' Rule, (2) IbE, and (3) genetic drift through isolation by distance (IbD) to isolate the effects of body size, local climatic conditions, and geographic distance, respectively, on the resting frequency of R. damarensis. Our results did not support genetic drift because there was no correlation between RF variation and geographic distance. Our results also did not support James' Rule because there was no significant relationship between (1) geographic distances and RF, (2) body size and RF, or (3) body size and climatic variables. Instead, we found support for IbE in the form of a correlation between RF and both region and annual mean temperature, suggesting that RF variation may be the result of environmental discontinuities. The environmental discontinuities coincided with previously reported genetic divergence. Climatic gradients in conjunction with environmental discontinuities could lead to local adaptation in sensory signals and directed dispersal such that gene flow is restricted, allowing lineages to diverge. However, our study cannot exclude the role of processes like phenotypic plasticity in phenotypic variation.

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.

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.

Sensory Drive Mediated by Climatic Gradients Partially Explains Divergence in Acoustic Signals in Two Horseshoe Bat Species, Rhinolophus swinnyi and Rhinolophus simulator

Geographic variation can be an indicator of still poorly understood evolutionary processes such as adaptation and drift. Sensory systems used in communication play a key role in mate choice and species recognition. Habitat-mediated (i.e. adaptive) differences in communication signals may therefore lead to diversification. We investigated geographic variation in echolocation calls of African horseshoe bats, Rhinolophus simulator and R. swinnyi in the context of two adaptive hypotheses: 1) James’ Rule and 2) the Sensory Drive Hypothesis. According to James’ Rule body-size should vary in response to relative humidity and temperature so that divergence in call frequency may therefore be the result of climate-mediated variation in body size because of the correlation between body size and call frequency. The Sensory Drive Hypothesis proposes that call frequency is a response to climate-induced differences in atmospheric attenuation and predicts that increases in atmospheric attenuation selects for calls of lower frequency. We measured the morphology and resting call frequency (RF) of 111 R. simulator and 126 R. swinnyi individuals across their distributional range to test the above hypotheses. Contrary to the prediction of James’ Rule, divergence in body size could not explain the variation in RF. Instead, acoustic divergence in RF was best predicted by latitude, geography and climate-induced differences in atmospheric attenuation, as predicted by the Sensory Drive Hypothesis. Although variation in RF was strongly influenced by temperature and humidity, other climatic variables (associated with latitude and altitude) as well as drift (as suggested by a positive correlation between call variation and geographic distance, especially in R. simulator) may also play an important role.

Development and characterization of 10 microsatellite markers in the Cape horseshoe bat, Rhinolophus capensis (Chiroptera, Rhinolophidae) and cross-amplification in southern African Rhinolophus species

Background

The Cape horseshoe bat, Rhinolophus capensis, is endemic to the Cape region of South Africa. Coalescent analysis of mitochondrial DNA sequence data suggests extensive historical gene flow between populations despite strong geographic variation of their echolocation call phenotype. Nevertheless the fine-scale genetic structure and evolutionary ecology of R. capensis remains poorly understood. Here we describe the development of 10 novel polymorphic microsatellite loci to investigate of the dispersal ecology of R. capensis and to facilitate taxonomic studies of Rhinolophus species in southern Africa.

Findings

We report 10 microsatellite primer pairs that consistently amplify scorable and polymorphic loci across 12 African rhinolophid species. Initial analysis of two populations of R. capensis from South Africa revealed moderate to high levels of allelic variation with 4–14 alleles per locus and observed heterozygosities of 0.450–0.900. No evidence of linkage disequilibrium was observed and eight of the loci showed no departure from Hardy–Weinberg equilibrium. Cross-species utility of these markers revealed consistently amplifiable polymorphic loci in eleven additional rhinolophid species.

Conclusions

The cross-amplification success of the microsatellites developed here provides a cost-effective set of population genetic marker for the study of rhinolophid evolutionary ecology and conservation in southern Africa.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1465-5) contains supplementary material, which is available to authorized users.

Effect of echolocation behavior-related constant frequency-frequency modulation sound on the frequency tuning of inferior collicular neurons in Hipposideros armiger

In constant frequency-frequency modulation (CF-FM) bats, the CF-FM echolocation signals include both CF and FM components, yet the role of such complex acoustic signals in frequency resolution by bats remains unknown. Using CF and CF-FM echolocation signals as acoustic stimuli, the responses of inferior collicular (IC) neurons of Hipposideros armiger were obtained by extracellular recordings. We tested the effect of preceding CF or CF-FM sounds on the shape of the frequency tuning curves (FTCs) of IC neurons. Results showed that both CF-FM and CF sounds reduced the number of FTCs with tailed lower-frequency-side of IC neurons. However, more IC neurons experienced such conversion after adding CF-FM sound compared with CF sound. We also found that the Q 20 value of the FTC of IC neurons experienced the largest increase with the addition of CF-FM sound. Moreover, only CF-FM sound could cause an increase in the slope of the neurons' FTCs, and such increase occurred mainly in the lower-frequency edge. These results suggested that CF-FM sound could increase the accuracy of frequency analysis of echo and cut-off low-frequency elements from the habitat of bats more than CF sound.