Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Highly Accessed Research article

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

Lizelle J Odendaal*, David S Jacobs and Jacqueline M Bishop

Author Affiliations

Department of Biological Sciences, University of Cape Town, 7701 Cape Town, South Africa

For all author emails, please log on.

BMC Evolutionary Biology 2014, 14:60  doi:10.1186/1471-2148-14-60

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/14/60


Received:24 October 2013
Accepted:20 March 2014
Published:27 March 2014

© 2014 Odendaal et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Abstract

Background

Across heterogeneous environments selection and gene flow interact to influence the rate and extent of adaptive trait evolution. This complex relationship is further influenced by the rarely considered role of phenotypic plasticity in the evolution of adaptive population variation. Plasticity can be adaptive if it promotes colonization and survival in novel environments and in doing so may increase the potential for future population differentiation via selection. Gene flow between selectively divergent environments may favour the evolution of phenotypic plasticity or conversely, plasticity itself may promote gene flow, leading to a pattern of trait differentiation in the presence of gene flow. Variation in sensory traits is particularly informative in testing the role of environment in trait and population differentiation. Here we test the hypothesis of ‘adaptive differentiation with minimal gene flow’ in resting echolocation frequencies (RF) of Cape horseshoe bats (Rhinolophus capensis) across a gradient of increasingly cluttered habitats.

Results

Our analysis reveals a geographically structured pattern of increasing RF from open to highly cluttered habitats in R. capensis; however genetic drift appears to be a minor player in the processes influencing this pattern. Although Bayesian analysis of population structure uncovered a number of spatially defined mitochondrial groups and coalescent methods revealed regional-scale gene flow, phylogenetic analysis of mitochondrial sequences did not correlate with RF differentiation. Instead, habitat discontinuities between biomes, and not genetic and geographic distances, best explained echolocation variation in this species. We argue that both selection for increased detection distance in relatively less cluttered habitats and adaptive phenotypic plasticity may have influenced the evolution of matched echolocation frequencies and habitats across different populations.

Conclusions

Our study reveals significant sensory trait differentiation in the presence of historical gene flow and suggests roles for both selection and plasticity in the evolution of echolocation variation in R. capensis. These results highlight the importance of population level analyses to i) illuminate the subtle interplay between selection, plasticity and gene flow in the evolution of adaptive traits and ii) demonstrate that evolutionary processes may act simultaneously and that their relative influence may vary across different environments.

Keywords:
Resting echolocation frequency; Neutral evolution; Mitochondrial DNA; Gene flow; Sensory ecology; Rhinolophidae

Background

In populations inhabiting heterogeneous environments adaptive trait divergence in the absence of gene flow is both predicted by theoretical models [1-3] and supported by a large number of studies demonstrating a generally negative relationship between levels of population connectivity and the degree of adaptive divergence in the traits under study [4-8]. This phenomenon is centred on the prediction that gene flow limits adaptive trait variation by homogenising gene pools that would otherwise diverge in response to selection under different ecological regimes [9-11]. The relationship is, however, clearly more subtle than this because empirical evidence also supports extensive adaptive variation in the presence of gene flow e.g. [12-15] highlighting the role of local selection gradients across different environments. If the selection gradient is steep enough to reduce immigrant fitness the homogenising effects of gene flow on trait variation are likely to be minimised [16]. Alternatively, if selection is weak and there are minimal costs to immigration, even nominal gene flow could constrain the evolution of trait variation and lineage divergence [17].

Undoubtedly there is a complex and dynamic relationship between selection, fitness and gene flow, but the interaction among these factors together with the role of trait plasticity is seldom considered in studies of adaptive population variation [18,19]. Phenotypic plasticity, environmentally mediated phenotypic changes without changes in allele frequency [20,21], has traditionally been viewed as a minor evolutionary process. This is so largely because environmentally induced phenotypic change does not influence the genes that an individual transfers to its offspring [19]; but see [22] for comment on the role of heritable epigenetic variation on phenotype evolution. Divergent selection cannot therefore act on genetic variants in a population and consequently phenotypic plasticity is thought to minimise rather than enhance selection [18]. Recently, however, there has been renewed interest in elucidating the role of plasticity in promoting diversification within and between populations and species [20,23-27]. This renaissance is due, in part, to the recognition that selection for adaptive plasticity can occur if it promotes successful dispersal and survival in different environments, increasing the potential for adaptive diversification under divergent selection regimes [19,28]. High gene flow between heterogeneous environments may then favour the evolution of increased phenotypic plasticity, over adaptive genetic divergence, because it would promote adaptation to new conditions within one or two generations [18,29]. Plasticity may in turn promote gene flow if dispersers are not selected against in their new environments [30]. Plasticity would allow a population or individuals to persist long enough in a new environment for selection to bring about the evolution of adaptation to the new environment by acting on the standing genetic variation among individuals [18,23].

Geographic variation in sensory traits i.e. traits used by organisms to perceive and respond to information about their environment [31] is well documented [32] and directly impacts on individual fitness e.g. via their role in resource acquisition and species and mate recognition [33]. Acoustic signalling for general communication as well as mate and food acquisition [34] is used by a range of organisms including fish [35], insects [36], anurans [37], birds [38] and mammals [39]. Extensive geographic variation in acoustic signals is common, and the role of both phenotypic plasticity [40,41] and sexually mediated selection have been implicated in promoting the evolution of population divergence [32,34]. In bats (Chiroptera) ultrasonic acoustic signalling in the form of echolocation is used for spatial orientation, prey detection and communication [42-44]. Bats echolocate over a wide range of frequencies [45] and the trait can be highly flexible both within [46,47] and between species [48]. This flexibility has been attributed to numerous factors including differences in foraging habitat [49,50], prey size [51], variation in body size [52], age [53] and sexual dimorphism [54]. In the horseshoe (Rhinolophidae) and leaf-nosed (Hipposideridae) bats a unique echolocation system has evolved that functions in both resource acquisition and intraspecific communication [55]. The signal design of horseshoe bats is characterised by long, high duty cycle calls which have a prominent constant frequency (CF) component preceded and followed by a brief frequency modulated (FM) component (45). During flight individuals compensate for Doppler shifts induced by their own flight speed by lowering the frequency of their emitted pulse. This ensures that the returning echo falls within the narrow frequency range of their acoustic fovea–a region of the cochlea with sharply tuned neurons sensitive to a unique frequency called the reference frequency [56]. Horseshoe bats are able to couple the frequency of the returning echo to their reference frequency independently of the size of Doppler shifts with extreme accuracy [57]. The reference frequency is always 150-200 Hz higher than the frequency these bats emit when stationary, often referred to as the ‘resting frequency’ (RF). Horseshoe bats typically forage for insect prey in narrow-space (i.e. highly cluttered) environments either on the wing (aerial hawking), or from a perch (flycatcher style) [58,59]. In both cases, the long CF component of their echolocation calls generate echoes from the flapping wings of flying insects that are characterised by amplitude and frequency modulations resulting in acoustic glints against a background of unmodulated echoes from background vegetation. Thus, the long CF calls emitted at a high duty cycle, together with Doppler shift compensation and the auditory fovea allow horseshoe bats to detect fluttering insects in cluttered space (reviewed in [60]). In horseshoe bats the RF is largely genetically determined [61]. However the peak frequency of juveniles is also partially influenced by that of their mothers [62,63] and the fine tuning of the frequency of the acoustic fovea of young horseshoe bats may have a learnt component via mother-to-offspring transmission [63,64]. Several recent studies investigating sensory variation in rhinolophid and hipposiderid bats [65-69] attribute observed variation predominantly to differences in diet [70] and environmental humidity (Humidity Hypothesis [71,72]). Such variation may be correlated with morphological features directly involved in echolocation production [47,73] and reception [74]. Sensory divergence may however also be an indirect result of adaptive changes in body size to prevailing environmental conditions (e.g. [71]). Among sympatric species call divergence may have evolved under regimes of social selection where character displacement maintains private bandwidths for intraspecific communication (Acoustic Communication Hypothesis [48,67,75]). Within a particular community, the echolocation frequencies used by any one species may therefore also be a consequence of the frequencies used by other species in that community [48,67]. Lastly, the physical structure of habitats (vegetation cover) can impose significant constraints on the range of frequencies bats use to successfully navigate within their environment and to detect prey, purely because of the physics of sound and sound transmission. Low frequency calls with longer wavelengths enable long-range detection of prey or other targets and are generally associated with open habitats. Conversely, higher frequency calls are more directional and provide greater resolution of targets, at least for bats using low duty cycle echolocation, at short detection ranges typical of highly cluttered habitats [45]. The Foraging Habitat Hypothesis [76] thus proposes that habitat clutter influences echolocation frequency and predicts that a gradient of increasing vegetation clutter selects for a gradient of increasing echolocation frequency. Most studies generally assume that observed variation in RF results from micro-evolutionary processes selecting for calls that are sharply tuned to the narrowband range of the acoustic fovea (± 1 kHz RF) [56,77,78]. Recent evidence, however, suggests that horseshoe bats and other high duty cycle echolocating bats exhibit some degree of temporal flexibility in RFs and associated changes in cochlear tuning [63,79] in response to both conspecifics [46,80] and local ambient noise [81].

More recently neutral evolutionary processes have become a focus of research centred on understanding echolocation divergence in bats [66,82]. Using the theoretical framework of population genetics, phylogeography and phylogenetics, the effects of gene flow and genetic drift on divergence in peak frequencies between divergent lineages has revealed relationships between echolocation divergence and e.g. female philopatry [66] geographic distance [7,69], and palaeoclimatic lineage divergence [82,83]. Here we use multiple data sets to investigate variation in echolocation frequency across the geographic range of the South African endemic Cape horseshoe bat, Rhinolophus capensis. We assess echolocation variation within a neutral evolutionary framework to (i) quantify the extent to which population genetic structure and gene flow contribute to variation in echolocation frequencies, and (ii) using a model-based approach explore the roles of habitat structure (Foraging Habitat Hypothesis [76]) and humidity (Humidity Hypothesis) in echolocation divergence. Our analyses reveal significant sensory diversification among R. capensis populations in the presence of extensive historic gene flow, and we discuss how selection, plasticity and gene flow may interact to characterize patterns of trait variation in this system.

Methods

Regional sampling of echolocation frequency and genetic variation

We measured echolocation parameters and collected tissue biopsies from populations across the distribution (Figure 1) of Rhinolophus capensis. Eleven populations were sampled across several major biomes in the region (Additional file 1: Figure S1, adapted from [84], no permission to reproduce maps was required) including Desert at the extreme northern edge, through Succulent Karoo and Fynbos in the centre of the distribution, and finally to Forest and areas of transition between Fynbos and neighbouring biomes (Albany Thicket, Nama-Karoo and Grassland) at the eastern limit of its range (Additional file 1: Figure S1). The geographic range of R. capensis is characterised by two distinct environmental gradients that influence regional rainfall patterns. The first is a latitudinal gradient of increasing aridity northwards. The second is a longitudinal seasonality shift from a predominantly winter to an aseasonal rainfall regime from west to east, and another shift to a summer rainfall regime at the extreme eastern edge of the distribution of the species [85] (Additional file 1: Figure S1). The resulting rainfall gradients translate into a clear habitat gradient from relatively open and sparse habitats in the north and west (Desert and Karoo Biomes), to more cluttered habitats in the east (Fynbos, Albany Thicket and Forest). Bats were captured at their roosts during the day with hand nets, or as they emerged from roosts at dusk using mist nets and/or harp traps. The age (adult or juvenile) and sex of each bat was recorded; juveniles were distinguished from adults by the presence of cartilaginous epiphyseal plates in their finger bones [86] and excluded from subsequent analyses. Individual body mass was measured. Seasonal and diurnal variation in body mass was controlled for by excluding pregnant females, sampling only in the southern hemisphere spring and summer and measuring bats only after their gut was emptied by keeping them overnight in soft cotton bags.

thumbnailFigure 1. Sampling localities and geographic distribution of sympatric horseshoe bats. Geographic ranges of R. capensis, R. damarensis and R. clivosus. The ranges of RF of R. damarensis[87] and R. clivosus[77] in areas of range overlap with R. capensis are also indicated.

Additional file 1: Figure S1. Map of the biomes of South Africa together with the geographic locations of the 11 populations of Rhinolophus capensis sampled in this study Biomes are from Rutherford et al. [84] and lines indicate the approximate positions of the different rainfall zones of South Africa: Solid line indicates the winter rainfall zone; dashed line indicates all year rainfall zone; and the rest of South Africa receives summer rainfall. The key to population acronyms are found in Table 1.

Format: PDF Size: 1.1MB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Echolocation calls were recorded from hand held bats positioned 10 cm in front of an Avisoft Ultrasound Gate 416 (Avisoft Bioacoustics, Germany) microphone connected to a laptop running Avisoft SasLab Pro software (sampling rate 500 kHz). Resting peak frequency (RF, where the CF component is stable and inter-pulse variation is low [47]) was recorded from hand-held individuals because it eliminates differences in peak frequency that may be due to Doppler shift compensation during flight [45]. RF is a reliable indicator of the reference frequency because the difference between the reference and resting frequency is stable in horseshoe bats (within 150-200 Hz) [60] Furthermore, when hunting from a perch rhinolophids emit their resting frequencies [88]. Recordings were slowed down by ten times and analysed using BatSound Pro software (Pettersson Elektronik AB, Sweden) with a sampling rate of 500 kHz. We measured the peak frequency (kHz) (frequency of maximum intensity) of the dominant second harmonic of the CF component determined from the fast Fourier transformation power spectrum (FFT = 1024, frequency resolution 684 Hz) with a Hanning window. To identify the typical echolocation parameters for each bat we calculated means of call duration, RF, lowest frequency of the FM component, bandwidth and inter-pulse interval from 5-10 high quality calls (amplitude of the signal at least three times higher than that of the background noise as displayed on the oscillogram). We chose the RF from an original call sequence that was most similar to the calculated mean parameters for all subsequent analyses. Thus the ‘true’ RF instead of a constructed statistical value was used in our analyses [73].

We also collected a tissue sample for DNA extraction and sequencing from most individuals for which we had echolocation recordings. Biopsy punches (3mm) were taken from the wing or tail membrane following [89]. Membranes were illuminated to ensure that no blood vessels were ruptured during sampling and tissues were stored in molecular grade (99%) ethanol at room temperature in the field, and at 4°C until extraction. Our capture, handling and tissue collection methods were approved by the Science Faculty Animal Ethics Committee (approval number: 2008/V18/LO) of the University of Cape Town. Total DNA was extracted using a DNeasy Blood and Tissue Kit (Qiagen) and stored at 4°C. We amplified a 520 base pair (bp) fragment of the mitochondrial d-loop using the primers N777 (5′TACACTGGTCTTGTAAAACC 3′) and E (3′ CCTGAAGTAGGAACCAGATG 5′) from [90] and [91] respectively. PCR conditions consisted of an initial cycle of 94°C for 5 minutes, followed by 35 cycles of 94°C, 50-55°C and 72°C each for 30 seconds and a final step of 72°C for 7 minutes. PCR products were checked on 1% agrose and gel purified using a Wizard SV Gel and PCR Clean-up System (Promega). Samples were sequenced in both directions using BigDye 3.1 chemistry on an ABI 3730 XL DNA Analyzer (Applied Biosystems). Chromatograms were edited and aligned using BioEdit v7.1.3.0 [92]. For each population the number of unique haplotypes, haplotype diversity (Hd) and nucleotide diversity (π) [93] were calculated using DnaSP v5.10.01 [94]. All unique haplotype sequences obtained in this study were deposited in GenBank (accession numbers: KF 175232-175270).

Influence of non-genetic factors on echolocation divergence

We conducted a Kolmogorov-Smirnov test for normality, Levene’s test for homogeneity of variances and regressed means against standard deviations for all grouping levels to ensure the data met the assumptions of subsequent analyses.

To quantify variation in RF within and between populations we calculated the mean ± standard deviation (SD) and the coefficient of variation (CV = SD/mean × 100) of each population.

Various intrinsic (e.g. body size and sex) and extrinsic (e.g. habitat characteristics) factors can influence the evolution of geographic variation in acoustic signals across a species range. To tease apart the effects of these factors on RF divergence we first conducted a general linear model (GLM) with RF as the dependent variable and population, sex and their interaction as categorical predictors. We then used SPSS (version 21, IBM, USA) to conduct linear mixed models (LMM’s) performed with type III sums of squares to determine the ecological factors which best predict RF in R. capensis. Mixed models are robust to unbalanced designs and thus data from single-sex populations could be included [95]. LMM’s estimate the effects of both fixed and random factors to models of data that are normally distributed [95]. Our GLM results revealed significant differences between populations and sexes, but no significant interaction between them (see Results). Population and sex were therefore included as random factors to control for spatial clustering of samples and intrinsic sexual differences in RF [95]. Fixed factors included biome category and body mass. Even though the presence or absence of congeneric species may also influence RF variation across different populations of horseshoe bats, we were unable to explicitly test the Acoustic Communication Hypothesis for several reasons. First, the range of Rhinolophus damarensis (mean RF of 85.1 kHz in the region of overlap [87]), only overlaps with that of R. capensis at LS (Figure 1) and therefore its presence cannot be replicated across populations. Second, the presence of R. damarensis and the Desert Biome category are collinear and thus cannot be included in the same model [96]. Lastly, with the exception of the LS site, R. capensis co-occurs with a larger congeneric species, Rhinolophus clivosus, over much of its range. Five geographical lineages of R. clivosus corresponding to previously described sub-species based on genetic, acoustic and morphological differences have been identified, and of these, R. capensis overlaps with two lineages which echolocate between 92.5 and 92.2 kHz [83] (Figure 1). At these frequencies the RF of R. clivosus is unlikely to influence the RF of R. capensis and R. clivosus was therefore excluded from our analyses.

We employed a model selection approach based on Akaike’s information criterion (AIC) to determine which model, out of a range of candidate models, best explained RF divergence in R. capensis[97,98]. The model with the lowest AIC value was considered the most parsimonious and the difference in AIC scores (Δi) were calculated to determine the likelihood that a given model was the best model relative to other candidate models [98]. A Δi value of zero indicates the best fit model; values up to two indicate models with substantial empirical support; values between four and seven indicate less support and models with values > 10 have essentially no support [98]. We also used Akaike weights (wi) to calculate the probability that a given model is the best among a candidate set of models. Thus, the best model has the lowest Δi and highest wi[98]. Once the best fit model was identified, we estimated the effects of each factor level on RF variation in R. capensis using the restricted maximum likelihood method. Then, to determine the importance of each variable included in the best model, we calculated the summed Akaike weight (w+) for all models containing that particular variable. The variable with the largest w+ is likely the most important variable in the model [99].

Our model selection results led us to explore the influence of habitat structure (typically vegetation cover) and body size on RF variation in greater detail. The Normalized Difference Vegetation Index (NDVI) was used as a measure of vegetation cover. NDVI is a suitable measure because it provides an estimate of above ground primary productivity [100] and it has been shown to be associated with a wide range of vegetation properties including photosynthetic activity [101] vegetation cover [102,103] and vegetation biomass [104]. Thus NDVI is commonly used to link vegetation dynamics to various aspects of animal ecology (reviewed in [105,106]). NDVI is a measure of the density of chlorophyll contained in vegetation and it is calculated as (NIR–RED)/(NIR + RED), where NIR is the near-infrared light, and RED is the visible-red light, reflected by vegetation and captured by the satellite. The values of NDVI range from -1 to 1, where negative values correspond to an absence of vegetation. Green and/or dense vegetation has high RED absorption together with high NIR, leading to high, positive NDVI values. In contrast, sparse vegetation absorbs substantially more NIR, leading to lower NDVI values [101]. Bare soils, snow and cloud have NDVI values close to zero [102,107]. We used the Expedited Moderate Resolution Imaging Spectroradiometer (eMODIS) Vegetation Indices dataset from NASA, which provides the maximum value for NDVI images for composites over a 10-day period at a resolution of 250 m from the year 2000 to present. Average NDVI values from a 20 km radius around each sampling site were extracted and compiled in ArcGIS 9.3.1 (ESRI®). Because echolocation frequency scales negatively with body size in horseshoe bats [52], variation in body size could cause concomitant changes in RF. Thus, to test whether RF variation is associated with differences in NDVI while controlling for the potential influence of body size, we used SPSS to conduct a hierarchal multiple regression analysis (HMRA). The first step of a HMRA is to add the independent variable that you wish to control for (in our case, body size). The second step examines the relationship between an independent (NDVI) and dependent variable (RF) while controlling for the effect of the first independent variable. We limited this analysis to males (n = 153) and evaluated the correlations between RF, body mass and NDVI as well as the change in the correlation coefficients between the model including body size and the model including NDVI but controlling for body size.

We also investigated whether there was a negative relationship between relative humidity and mean RF as predicted by the Humidity Hypothesis [108]. Relative humidity (%) data for each sampling site were obtained from the literature or from the nearest weather stations and provided by the South African Weather Service.

Like most other rhinolophids, R. capensis forage in or near highly cluttered habitats [67,109]. While the long CF portion of their echolocation calls allow the detection of fluttering prey against structurally complex backgrounds like dense vegetation [88], differences in vegetation cover between sampling sites could still influence the sound transmission properties of the echolocation call frequency used in each habitat. To better understand the effects of vegetation cover on variation in echolocation frequencies we calculated the mean detection distances for prey and vertical background targets (e.g. leafy vegetation edge) for each population according to the method developed in Stilz and Schnitzler [110]. This method depends on the dynamic range and frequency of the sonar system, local atmospheric conditions and target type. The dynamic range was calculated as the difference between peak intensity (dB SPL) at 1m (79.1 dB SPL for R. capensis at De Hoop; Jacobs and Parsons, unpublished data) and the auditory threshold of the bat (assumed to be 0 dB SPL for horseshoe bats [70,111]). We assumed the different populations had similar dynamic ranges. The different prey size categories tested in [110] were derived from Houston et al. [51] and included small, medium and large categories–all within the size range of prey consumed by R. capensis (2 mm-19 mm [67]). Mean minimum temperature (°C) data for each population were obtained from the nearest weather stations and provided by the South African Weather Service.

Spatial population structure and patterns of gene flow

Geographic variation in phenotypic traits may be a consequence of neutral evolutionary processes, particularly when dispersal distances result in a pattern of predominantly nearest-neighbour gene exchange [112,113]. To better understand the role of random genetic drift in the evolution of RF variation we used a number of statistical approaches to (i) determine the degree to which genetic variation is spatially structured in R. capensis and (ii) quantify levels of historic gene flow among populations. We first explored the evolutionary relationships among mtDNA haplotypes to determine whether observed relationships reflected either the geographic sampling of populations or specific biome discontinuities across the range of the species by constructing a Neighbour-net network in SplitsTree v4.12.6 [114] using uncorrected ‘p’ distances. Network analysis allows reticulations among branches and is highly suitable for analysing evolutionary relationships among populations [115].

The spatial clustering of samples and distribution of genetic variation were also investigated using BAPS v6.0 [116] and Analysis of Molecular Variance (AMOVA) in GenAlEx v6.5 [117]. The AMOVA quantified partitioning of genetic variation within and among populations and biomes. The significance of variance components and ΦST (a measure of population genetic differentiation analogous to the fixation index FST[118,119]) were estimated with 1000 random permutations. To test for the presence of discrete evolutionary lineages we estimated the most probable number of genetic clusters (K) among populations using a spatially explicit Bayesian clustering mixture model in BAPS. We performed spatial clustering for groups with K ranging from 4-11. The analysis was repeated ten times for each maximum K and the log marginal likelihood value for each genetic partition was evaluated.

To obtain estimates of historic gene flow among populations we used a maximum likelihood method based on the coalescent implemented in MIGRATE-N v3.3.2 [120]; MIGRATE uses an equilibrium model that estimates migration rates averaged across the coalescent history and simultaneously estimates Θ, the effective population size scaled by mutation rate where Θ = Neμ, together with pairwise migration rates summarised as M = m/μ, where m is the effective immigration rate per generation between populations. Banding data from European horseshoe bats reveal generally small home ranges where maximum dispersal distances rarely exceed 100km over the course of an individual’s life time [121]. If similar, dispersal in R. capensis likely occurs over relatively short distances. Thus we estimated Θ and M using a custom designed migration matrix where migration was only allowed between neighbouring populations. Populations separated by large geographic distances were not directly connected except in cases where unique haplotypes were shared. Starting values for all parameter estimates were initially obtained using FST[122] and the following search parameters were used: 10 short chains with 500 000 gene trees sampled and 5000 trees recorded; three long chains with 50 million sampled trees of which 50 000 were recorded. The first 10 000 trees were discarded as burn-in and a static heating scheme with six temperatures and a swapping interval of one was used. The results were averaged over five replicate runs.

Comparing echolocation divergence and genetic distance

Genetic divergence among populations is dependent on both the degree of physical isolation and levels of connectivity between them [123]. To test whether genetic and RF differentiation is associated with geographic isolation (isolation by distance) we used Mantel [124] and partial Mantel tests to explore the relationships between genetic structure, RF divergence and geographic distance [125]. Matrix correlations were calculated in GenAlex v6.5 [117] and XLStat (v2013, Addinsoft) with 1000 random permutations. The three matrices analysed were pairwise geographic distances (km) calculated as straight line distances from geographic coordinates using the program Geographic Distance Matrix Generator v1.2.1 (Ersts, Internet); genetic distance using Slatkin’s linearized ΦST[126] and RF differences (kHz) among populations. Log-transformed geographic distances were regressed against genetic distance and peak frequency difference. Our partial Mantel test determined whether RF difference was associated with genetic distance while controlling for the effect of geographic distance (isolation by adaption).

Results

Effects of morphology, sex and ecology on RF divergence

We measured the RF of 248 individuals across multiple biomes (Additional file 1: Figure S1) and observed a clinal increase in mean RF across the distribution of R. capensis ranging from 75.7 kHz (LS) in the west, to 86.5 kHz (BAV) in the east (Figure 2, Table 1). Resting frequencies differed significantly among populations (GLM: F18, 486 = 120.5, P < 0.001; Tukey HSD tests: P < 0.005) with the exception of KNY which used similar frequencies to HDH and DHC (Tukey HSD tests: P > 0.05; Table 1). Sex significantly influenced echolocation variation within populations with females emitting higher frequencies than males (GLM: F3, 172 = 33.5, P < 0.001; Table 1). However, we detected no significant interaction between sex and population (GLM, F18, 486 = 1.1, P > 0.05), suggesting that the degree of sexual dimorphism in RF was identical across populations. The CV of call frequencies for each population was low and at most 1% (Table 1).

thumbnailFigure 2. Box and whisker plots of median resting frequencies across sampled populations of Rhinolophus capensis. Colours represent biome category and a key to acronyms is given in Table 1.

Table 1. Biome category and mean (± SD) body mass and RF for each sampled population

The most parsimonious model explaining variation in RF included the factors biome category (F5, 12 = 40.3, P < 0.005) and body mass (F1, 232 = 10.4, P < 0.001) of which the former was the most important variable in the model (w+ biome = 0.99, w+ body mass = 0.95) (Table 2). Generally, bats from the Desert Biome use significantly lower frequencies than Succulent Karoo bats (point of reference selected by SPSS). Also, bats inhabiting the Forest Biome or regions comprised of multiple biomes use significantly higher frequencies (Table 3) because long–range detection is not an advantage in such habitats. In the first stage of our HMRA, body size significantly influenced RF (R2 = 0.121, F1, 150 = 20.6, P < 0.0001) but it only explained 12% of the variation in RF (Figure 3). Bats from LS use relatively lower RFs given their body size compared to other sampled populations. The inclusion of NDVI in the second stage of the regression model significantly increased the proportion of variance explained in the model (Δ R2 = 0.68, F1, 149 = 528.7.9, P < 0.0001), with NDVI accounting for 80% of the variation in RF (R2 = 0.80, F2, 149 = 310.9, P < 0.0001) after controlling for the effect of body size (Figure 4).

Table 2. Model selection results for three candidate models explaining variation in resting frequencies

Table 3. Restricted maximum likelihood estimates and confidence intervals for the best-fit linear mixed effects model

thumbnailFigure 3. The regression of mass (g) on resting frequency (kHz) for male R. capensis (n = 153) from 11 populations. Colours represent the biome category of each population and shapes represent different populations. The key to the abbreviations of population names are given in Table 1.

thumbnailFigure 4. The relationship between RF and NDVI (as a proxy for habitat clutter). The regression of RF (kHz) and NDVI across populations of R. capensis. Habitat photographs show the vegetation cover and structure for each population, A: Lekkersing (Desert), B: Steenkampskraal (Succulent Karoo), C: De Hoop (Fynbos), D: Knysna (Forest). A key to acronyms used are given in Table 1.

We also found no relationship between relative humidity and RF (R2 = 0.04, P > 0.5) but we did discover significant differences in detection distances for large prey and vegetation edge (GLM: F16, 306 = 275, P < 0.01, Tukey HSD tests: P’s < 0.05; Additional file 2: Table S2), but not for small or medium sized prey (Tukey HSD tests: P’s > 0.05; Additional file 2: Table S2), between populations. There was, however, no difference in detection ranges associated with vegetation or large prey between ZPK, SKK and DHL and between TF and BKL (Tukey HSD tests: P’s > 0.05; Additional file 2: Table S2). Bats inhabiting more open habitats used lower frequencies and had longer detection distances than bats in more cluttered habitats, with LS bats (Desert Biome) having the longest detection distances of both large insects and background vegetation (Additional file 2: Table S2).

Additional file 2: Table S2. Climatic variables for each population and the mean detection distances for prey and background vertical targets (leafy vegetation edge). Detection ranges calculated from the method of Stilz and Schnitzler [110] (http://134.2.91.93/~peter/calculator/range.php webcite). The size range of prey taken by R. capensis at De Hoop is 2-19 mm [67] and this covers the range of small, medium and large prey in [110]. Climatic data were obtained from the nearest weather stations and provided by the South African Weather Service.

Format: DOCX Size: 16KB Download fileOpen Data

Spatial population genetic structure and patterns of gene flow

A total of 39 unique mtDNA haplotypes were identified from 203 individuals (Additional file 3: Figure S2). Haplotype diversity ranged from 0.57 (SPH) to 0.88 (BKL) and the highest number of unique haplotypes (n = 12) was found in the Baviaanskloof area (BAV) (Additional file 4: Table S2). Most populations shared haplotypes with their nearest neighbours and a few haplotypes were shared between more distant populations (e.g. between ZPK and BAV) (Additional file 3: Figure S2). We identified three genetically isolated populations (i.e. LS, TF, and KNY) consisting largely of unique mitochondrial lineages (Additional file 3: Figure S2).

Additional file 3: Figure S2. The distribution of 39 unique haplotypes across the 11 populations of Rhinolophus capensis sampled in this study. Key to acronyms given in Table 1.

Format: PDF Size: 610KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 4: Table S2. Genetic variability in 11 populations of Rhinolophus capensis based on 519 bp of the mitochondrial control region. Haplotype diversity (Hd), number of haplotypes, nucleotide diversity (π) and number of polymorphic sites are shown.

Format: DOCX Size: 15KB Download fileOpen Data

Complex network relationships characterise the evolutionary history of R. capensis populations sampled in our study. The neighbour-net network revealed numerous reticulations between haplotypes, suggesting several alternative evolutionary pathways among them ([127], Additional file 5: Figure S3). While the network recovered a number of clear clades these were not generally structured by biome or geographic proximity. Only haplotypes from the desert biome and populations occurring in regions where multiple biomes are connected formed discrete clusters (Additional file 5: Figure S3).

Additional file 5: Figure S3. Neighbour-net network based on p-corrected distances of the 39 unique haplotypes isolated in this study. Each circle represents a unique haplotype coloured according to the biome/s in which it occurs. Pie graphs indicate where haplotypes are shared across several biomes.

Format: PDF Size: 43KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Our investigation of hierarchal population genetic structure revealed significant partitioning of genetic variation at all three levels but most variation occurred within populations (ΦST = 0.54) rather than between populations (ΦST = 0.33) or biomes (ΦST = 0.13) (P’s < 0.005). Bayesian clustering identified four spatially explicit genetic clusters in our data (Figure 5). Only individuals situated at opposite ends of the species distribution (populations LS and TF), were assigned to unique clusters while the nine populations between them were assigned to two spatially structured clusters (Figure 5). Furthermore, estimates of historical migration rates (M) revealed generally asymmetrical patterns of gene flow between populations situated in different biomes (Figure 5, Additional file 6: Table S3). The main source populations were ZPK, BKL and BAV and gene flow generally occurred between neighbouring populations. There was also evidence for long distance gene flow from BAV to BKL (430 km straight line distance) and from DHC to BAV (370 km) but this occurred at relatively low levels (Figure 5). Together, these results support a pattern of relatively high regional connectivity and hence minimal fine-scale population structure in the R. capensis.

thumbnailFigure 5. Estimated migration patterns among populations of Cape horseshoe bats. The thickness of the arrows indicates the relative migration rates; values indicate M, the number of immigrants per generation scaled by mutation rate. The 95% confidence intervals are provided in Additional file 6: Table S3. Inset: Voronoi tessellation showing the assignment of populations to four spatially defined genetic clusters. A key to acronyms used are given in Table 1.

Additional file 6: Table S3. Lower and upper profile likelihood percentiles of M, the number of immigrants per generation scaled by mutation rate, calculated in Migrate-N (Beerli [120]).

Format: DOCX Size: 14KB Download fileOpen Data

Influence of genetic distance and geography

Pairwise differences in both RF and genetic distance were characterised by significantly positive relationships with geographic distance (Figure 6A, B). RF and genetic distance co-varied with geographic distance and followed a general model of isolation by distance (IBD). Geographic distance, however, explained only a relatively small proportion of the variation in RF and genetic distance (44% and 30%, respectively; Figure 6A, B). We also found a significant correlation between RF difference and genetic distance among populations (Figure 6C) and this remained significant when we controlled for the effect of geographic distance (Partial Mantel Test: R2 = 0.075, P < 0.001). Despite the significance of the correlation, genetic distance only explained a small proportion of the variation in echolocation frequencies (7.5%).

thumbnailFigure 6. Results from Mantel and Partial Mantel tests. Pairwise geographic distance versus A: RF difference, and B: genetic distance; and C: genetic distance versus RF difference.

Discussion

Our results reveal significant geographic variation in the resting echolocation frequencies (RF) of R. capensis despite substantial historical gene flow across the distribution of this species. Body size, genetic distance and geographic distance play minor roles in the evolution of geographic variation in RF in R. capensis. Further, despite a lack of strong geographic structure in mitochondrial lineages, biome was identified as the best predictor of RF. In support of this we found a significant relationship between RF and increasing vegetation clutter from west to east across the range of this species. Our results suggest that the evolution of geographic variation in RF in the face of homogenizing gene flow in R. capensis was probably influenced by selection for lower echolocation frequencies in less cluttered habitats. However, the relationship we found between RF and habitat type is solely correlative, and therefore a thorough comparison of foraging behaviours of bats between different habitats is required. Thus, alternative evolutionary processes e.g. selection for discrete frequency bands or phenotypic plasticity cannot be excluded at this stage.

Disparities between RF and population genetic structure: trait diversification in the presence of gene flow

Studies investigating the various evolutionary forces shaping phenotypic and genetic divergence between populations often describe trait variation in the context of substantial population structure and limited gene flow [7]. While phylogeographic patterns of maternally (mtDNA) and bi-parentally (e.g. microsatellites) inherited genetic markers usually concur, discordance between phenotypic and genetic structure is also reported. This is usually attributed to demographic processes that characterize species, such as sex-biased dispersal [128] or secondary contact following historical isolation [47] (reviewed in [129]). Similar discordance has been reported in several European and Asian horseshoe bats either as a result of male-biased dispersal and female philopatry (e.g. R. pumilus: [66]) or historical introgression of mtDNA (R. pearsoni: [130]; R. sinicus: [131]; or nuclear genomes (R. yunanensis to R. pearsoni: [130]) between sister lineages. Despite recent advances in testing the evolutionary processes that shape contemporary population genetic structure in horseshoe bats (e.g. [132]), few studies have specifically evaluated sensory variation within a phylogeographic or population genetic framework in this genus [66,69]. In our study, we expected RF divergence in the R. capensis to reflect mtDNA structure given that (i) the fine tuning of the echolocation frequency of young horseshoe bats are partly learned from their mothers[63], (ii) female philopatry and male dispersal characterize other horseshoe bats studied to date [66,133], and (iii) the degree of RF divergence observed among populations (range: 1-11 kHz) is similar to that reported for other high duty cycle bats which have corresponding significant genetic structure among maternal lineages [77]. Our results instead reveal minimal genetic structure amongst R. capensis populations. A Bayesian clustering analysis identified four dominant genetic lineages which broadly reflect estimated patterns of regional gene flow that are not limited by biomes; notably only LS and TF (populations at the opposite extremes of the distribution of the species) are identified as unique genetic clusters (Figure 5). A number of mitochondrial haplotypes were also found to be shared between geographically distant populations; these may reflect long-distance dispersal events or the retention of common ancestral polymorphisms. Mitochondrial data clearly reveals a recent evolutionary history of complex reticulations in R. capensis, suggesting that gene flow and not incomplete lineage sorting is responsible for the observed genetic structure or lack thereof.

One caveat of our analytical approach is the statistical non-independence of our populations due to the substantial gene flow we detected between them [134]. Populations that are more closely related to one another or exchange higher numbers of migrants are likely to have similar phenotypic trait values irrespective of the local selection pressures they experience [134]. Ideally, the inclusion of our calculated migration matrices between populations in our ecological models would control for the effect of gene flow, but the development of effective computational methods to include the complex reticulate relationships among population’s remains challenging [134,135]. Nonetheless, the spatial distribution of the four dominant genetic lineages does not reflect the broad pattern of RF variation among populations of R. capensis. For example, we detected considerable gene flow between populations at SKK, ZPK and DHL, but no direct gene flow between these populations and the population at DHC, even though they were all assigned to the same genetic cluster (Figure 5). Bats from these genetically similar populations use RF’s ranging from 80.6-84.6 kHz. Thus the discordance between patterns of genetic and RF variation among R. capensis lineages suggest that gene flow does not significantly contribute to RF variation across the distribution of R. capensis, although this remains to be explicitly tested. Furthermore, although RF is also characterised by IBD, geographic distance only accounted for a minor proportion of the variation in RF (Figure 6A and B). Interestingly, RF divergence is also positively correlated to some degree with genetic distance, even after controlling for geographic distance (Figure 6C), suggesting that local environmental conditions may influence sensory differentiation in areas with some degree of restricted gene flow [136].

However, estimates of maternal gene flow support significant regional connectivity in the recent past and is at odds with the pattern of structuring in RF we observe across populations. This is in contrast to previous studies of high duty cycle bats where population genetic structure generally reflects the variation of sonar frequencies across often widespread species (e.g. Pteronotus parnellii: [77], R. clivosus: [73], R. hildebrandtii: [82], R. rouxii: [137]). As we discuss below, this anomaly may be a consequence of the complex interactions between gene flow, diversifying selection and environmentally mediated selection for phenotypic plasticity. Phenotypic plasticity can be advantageous if it results in the expression of different phenotypes that increase an individual’s fitness in diverse environments [19,138]. In this way plasticity can minimise the costs incurred from dispersal into environments with different selection regimes [18,29]. If plasticity influences RF diversification, it may lead to a situation in which selection does not significantly constrain gene flow among populations, leading to adaptive phenotypic variation in the presence of gene flow [18].

Body size and habitat structure predict RF: influence of diversifying selection and adaptive phenotypic plasticity

The frequency of acoustic signals scales negatively with body size in a range of organisms [139-141] including bats [52]; larger bat species produce echolocation calls of lower frequencies than smaller bat species. This relationship also exists within R. capensis with larger bats generally using lower RF’s. However, body size only explained a minor proportion of the variation in RF because concomitant changes in body size and RF were not evident across all populations. For example, although BAV bats are not the smallest, they use the highest echolocation frequencies. In contrast, individuals in populations at both ends of the distribution of this species (LS in the west and TF in the east) were similar in size but echolocated at the lowest and second highest frequencies respectively (Table 1, Figure 3). This may explain why previous research on echolocation variation in the species did not support a relationship between body size and RF, and instead suggested the decoupling of echolocation divergence from the evolution of body size [73]. This result was likely an artefact of under sampling trait variation in R. capensis and highlights the difficulties of inferring evolutionary processes from data sets which inadequately sample the true distribution of a trait. Our results here suggest that the allometric relationship between body size and RF collapses in populations situated towards the edge of the range of this species. This is perhaps not unexpected given that range edges and ecotones provide novel environments to which species can become locally adapted, leading to significant phenotypic divergence of edge populations from those at the centre of the distribution of a species [142]. This may explain why populations at either end of the species’ range (LS in the west, and TF in the east) were identified as unique genetic clusters, while the nine sampled populations between them were assigned to only two genetic clusters (Figure 5).

Paleoenvironmental change from the Miocene to Pleistocene profoundly impacted the evolution of population divergence and speciation in a wide range of southern African taxa; several studies report a strong link between divergent genetic lineages and the biomes or ecogeographical regions of southern Africa for various organisms including reptiles [143], invertebrates [144] and small mammals [145,146] including bats [82,87]. In our study we did not find any clear population genetic divergence with topographical features. Neither did we find any support for the humidity hypothesis which proposes that divergence in RF is the result of selection against higher frequencies in humid environments [108]. Instead we found that bats in the Desert Biome used significantly lower frequencies than those occupying Forest and areas of transition between multiple biomes (Figure 4, Table 3). This clinal increase in RF across populations of R. capensis from west to east may be the result of the gradient of increasing vegetation cover and density from west to east (Figure 4). At LS in the west the vegetation was very sparse and consisted of low shrubs (< 1m in height). In the east the vegetation ranged from dense Fynbos to Forest (Figure 4). Thus, there is a steep habitat gradient between LS (genetic cluster 1) and its nearest neighbours (SKK and ZPK: genetic cluster 2), and selection may therefore dominate in this region. In contrast, plasticity may be favoured as an explanation for RF variation amongst the other populations because the habitat gradients are not as steep. A notable exception is that of the population at TF. At TF (genetic cluster 4) bats use similar RF’s to other populations situated in ecotones (genetic cluster 3), and yet appears to be relatively isolated genetically (Figure 5). TF is situated at the interface between aseasonal and summer rainfall zones [147] which, together with the intrinsic habitat heterogeneity of ecotones, may serve as a significant barrier to gene flow. This pattern has been observed for various other species in this region (e.g. four-striped mouse, Rhabdomys pumilio[145]; forest shrew, Myosorex varius[146]).

The positive correlation we found between RF and NDVI suggests that variation in the degree of habitat clutter might explain variation in RF. Across our sampling sites, the mean detection distance for large prey and background vegetation edge was significantly different between populations; bats occupying more open habitats have lower RF’s and thus longer detection distances than those in more cluttered habitats, allowing them to detect larger prey or background targets at greater distances. While this result was statistically significant, LS bats only had a 10 cm and 50 cm greater detection distance than their nearest neighbours for large prey and vegetation edge, respectively (Additional file 2: Table S1). Although all rhinolophids supposedly fly close to vegetation and may experience even relatively open habitats as cluttered, the 50 cm greater detection distances of the vegetation edge may be advantageous during orientation and commuting flight in the sparse vegetation of LS where the distance between clumps of vegetation are greater than in the Fynbos or the Forest (Figure 4). Further, prey density is also likely to be lower at LS and selection is likely to favour the evolution of lower RF’s that would allow the detection of larger prey at greater distances. There is sufficient empirical evidence to suggest that even rhinolophids, constrained by their echolocation to hunt in narrow space [58], display some degree of flexibility in the foraging habitat they exploit [148-151] or the foraging style they adopt (aerial hawking vs. perch hunting) [59,152] as a result of resource partitioning [150,153], habitat structure [148] or seasonal changes in prey resources [150]. At least one species of rhinolophid also appears to vary its echolocation frequency in response to different degrees of clutter [148]. Within the same nature reserve, greater horseshoe bats (R. ferrumequinum) use a variety of habitats with differing degrees of clutter and use significantly lower echolocation frequencies in relatively open habitats than in cluttered habitats [148]. It is likely therefore that R. capensis uses both aerial hawking and perch hunting styles to different degrees in the different habitats perhaps also altering its call frequency to deal with different degrees of clutter. The distinctly lower call frequency at LS could possibly be explained by the pronounced clutter gradient between LS and the other habitats occupied by R. capensis. However, our detection distance calculations must be interpreted with caution because different populations may have different echolocation call intensities, which could greatly influence the calculation of maximum detection distance [154].

Although differences in habitat structure provide a compelling explanation for different RF’s in R. capensis, an interesting anomaly needs to be accommodated; R. damarensis (forearm length 49.5 ± 1.7 mm; n = 20) is sympatric with R. capensis in the extremely arid area of LS but uses frequencies (mean ± SD = 85.4 ± 1.4 kHz; range = 82–89 kHz; [87]) as high as those used by R. capensis in highly cluttered habitats e.g. TF. It seems reasonable to assume that selection would have favoured similar call frequencies, and therefore detection distances, in R. damarensis and R. capensis where they are sympatric in the arid region of LS. However, it is possible that R. damarensis exploits a different foraging niche to R. capensis in this area of sympatry. Alternatively, this species may compensate by using higher intensity calls to achieve more or less the same detection distances as R. capensis[154]. Of course it is possible that the presence of R. damarensis at LS has selected for the low frequencies we observe in R. capensis. Observed frequencies are much lower in the LS population than in the other arid populations (Figure 4) despite similarities in body size to the other arid populations (Table 1). On the basis of its body size, R. capensis at LS should echolocate at approximately 82 kHz (Figure 3). Instead we observe a mean RF of 75.7 kHz. Bats at LS may have shifted their frequency below 82 kHz to avoid acoustic overlap with R. damarensis (82-89 kHz), and maintain effective intraspecific communication. However, the closest known roosts of these two species are 40 km apart and it is not known if the foraging areas of these two species overlap, i.e. if they are syntopic, as required by the Acoustic Communication Hypothesis.

Nevertheless, character displacement, mediated by some form of resource partitioning or local adaptation, in sonar parameters can contribute to the initial stages of lineage divergence by causing populations in sympatry with heterospecifics to diverge from their respective conspecific populations in allopatry. The result would be reduced gene flow between populations found in different selective regimes or heterospecific assemblages [155]. For example, Lemmon [37] showed that acoustic traits important for female preference and mate choice in populations of the chorus frog, Pseudacris feriarum, diverged to maximize differences from the heterospecific assemblage present, ultimately promoting reproductive isolation between conspecific populations via sexual selection. Ecologically adaptive traits can also promote divergence if divergence has a pleiotropic effect on reproductive isolation via assortative mating; so called ‘magic traits’ [34]. In R. philippinensis assortative mating has evolved between size morphs as a by-product of selection for different frequencies used to exploit different prey sizes [70]. The significant and consistent sexual dimorphism we observe in the RFs of R. capensis (female’s echolocate at higher frequencies than males: Table 1), may serve a role in sex-specific communication in the species. Recent experimental evidence reveals that horseshoe bats (R. euryale and R. mehelyi) and emballonurid bats (Saccopteryx bilineata) are indeed able to recognise the sex of conspecifics based on their echolocation calls [44,156]. Thus, the limited gene flow between LS and other populations may be a consequence of LS bats not effectively recognizing other R. capensis as potential mates. Divergence in RF may have under-appreciated consequences for the evolution of reproductive isolation via female preference for male RFs in different populations of horseshoe bats. Evaluating female preference in LS bats for local versus allopatric RFs may provide intriguing insights into the causes and consequences of sexual selection in horseshoe bats. Finally, because R. damarensis only occurs in sympatry with R. capensis at a single site (LS), we were not able to test whether habitat or the presence of R. damarensis offered a better explanation for RF differentiation in R. capensis. An experimental approach is required to test whether R. capensis has shifted its RF to avoid acoustic overlap with R. damarensis. Such an approach would evaluate whether R. capensis from LS responds differentially to the echolocation calls of heterospecifics and acoustically divergent conspecifics using playback experiments.

At a ‘local to regional’ scale geographic distance is not a significant barrier to gene flow in R. capensis, and the evolution of sensory divergence in the presence of this gene flow may also reflect a degree of adaptive phenotypic plasticity in RF. Despite the tight coupling between RF and the acoustic fovea in high duty cycle bats [45], empirical studies have shown that species are able to shift their RF’s in response to both neighbouring conspecifics (maximum shift 3.9 kHz: [46]) and different ambient noise conditions (maximum shift <0.5 kHz: [81]). Such small shifts in frequency may explain the range of RF variation in the southern and eastern populations of our study (approximately 3 kHz) where plasticity in response to slightly varying degrees of vegetation clutter towards the east might occur. However, it is unlikely to explain the 9 kHz shift in the Lekkersing bats. Nonetheless, it appears that southern and eastern populations of R. capensis use RF’s within the best hearing range of the acoustic fovea of their nearest neighbours, possibly facilitating gene flow and promoting relatively flexible RF’s in these populations. While small shifts in the acoustic fovea and its corresponding reference frequency are possible in high duty cycle bats [46,81], we do not know the precise limits of the flexibility of the acoustic fovea. Long-term experimental studies evaluating the change in RFs in response to the RFs of bats from acoustically divergent populations may shed light on the degree of plasticity in this system. Alternatively, the relative influence of plasticity versus selection can be evaluated indirectly. For example, selection may better explain our observations if variation in functional genes involved in hearing co-varies with RF variation across populations. Recent studies reveal a wide range of candidate hearing genes which show strong signals of ancestral positive selection in the evolution of echolocation in bats and cetaceans [157,158]. Selection may also be favoured as an explanation for variation in RF if there is a strong correlation, independent of variation in body size, between RF and morphological features directly involved in echolocation production and emission (such as dorsal nasal chambers) [47]. Previous research investigating morphological correlates of RF in R. capensis revealed that RF is best predicted by nasal chamber length [73], but a thorough evaluation of skull morphology variation across the distribution of the species is required. The social life of bats may also influence the relative roles of diversifying selection versus plasticity across the distribution of a species. If bats are able to recognise conspecific calls from a range of acoustically divergent populations this might suggest that selection for some degree of plasticity in the trait is also favoured. Classic playback experiments can be used to assess the sensitivity of individuals to the range of frequencies exhibited by a species. These two hypotheses are clearly not mutually exclusive; and their relative influences across the highly heterogeneous environments of R. capensis certainly merit further attention.

Conclusions

This study reveals significant sensory trait variation in the Cape horseshoe bat despite substantial historical gene flow. While genetic and geographic distances do influence sensory variation to some extent, our results suggest that differences in habitat complexity across the range of the Cape horseshoe bat may be the dominant driver of sensory differentiation in this system. Classical divergent selection together with some degree of phenotypic plasticity may be responsible for RF variation in the presence of gene flow. However, an investigation of the variation in foraging behaviour within and between populations of the Cape horseshoe bat is required to support our results. Nonetheless, our findings highlight the importance of population level analyses to elucidate the complex interactions among selection, plasticity and gene flow in the evolution of adaptive trait variation and reveal a number of interesting avenues for future research into the evolution of this remarkable sensory system.

Competing interests

All authors have declared that no competing interests exist.

Authors’ contributions

LJO, DSJ and JMB conceived of and designed the study and carried out the fieldwork. LJO generated the molecular data, and conducted all genetic and ecological analyses with input from JMB and DSJ. All authors contributed to the writing of the manuscript and read and approved the final manuscript.

Acknowledgements

We thank B. D. Kruger, H. Babiker, C. Hartnick, T. Haupt, H. Craven, S. Stoffberg, E. Kelly, C. Sirami, D. Peinke, O. Smit, M. Cunnama, S. Rakers, L. Lange, A. Mead, A. Thomas, S. Jordaan, F. P. D. Cotterill, M. Brigham, A. Brigham and A. Doty, for assistance in the field. Thanks also to G. Mutumi for assistance with ArcGIS, and H. Craven for permission to publish her photograph of Rhinolophus capensis taken while assisting with fieldwork. This research was funded by grants to LJO (National Research Foundation, Harry Crossley Foundation Fellowship and University of Cape Town Research Associateship) and DSJ from the South African Research Chair Initiative of the Department of Science and Technology and administered by the National Research Foundation of South Africa. We also thank two anonymous reviewers for their constructive comments on previous versions of this manuscript.

References

  1. Haldane JBS: The theory of a cline.

    J Genet 1948, 48:277-284. PubMed Abstract | Publisher Full Text OpenURL

  2. Slatkin M: Gene flow and selection in a cline.

    Genetics 1973, 75:733-756. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Hendry AP, Day T, Taylor EB: Population mixing and the adaptive divergence of quantitative traits in discrete populations: a theoretical framework for empirical tests.

    Evolution 2001, 55:459-466. PubMed Abstract | Publisher Full Text OpenURL

  4. Tobias JA, Aben J, Brumfield RT, Derryberry EP, Halfwerk W, Slabbekoorn H, Seddon N: Song divergence by sensory drive in Amazonian birds.

    Evolution 2010, 64:2820-2839. PubMed Abstract OpenURL

  5. Wang IJ, Summers K: Genetic structure is correlated with phenotypic divergence rather than geographic isolation in the highly polymorphic strawberry poison‒dart frog.

    Mol Ecol 2010, 19:447-458. PubMed Abstract | Publisher Full Text OpenURL

  6. González C, Ornelas J, Gutiérrez-Rodríguez C: Selection and geographic isolation influence hummingbird speciation: genetic, acoustic and morphological divergence in the wedge-tailed sabrewing (Campylopterus curvipennis).

    BMC Evol Biol 2011, 11:38. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  7. Puechmaille SJ, Gouilh MA, Piyapan P, Yokubol M, Mie KM, Bates PJ, Teeling EC: The evolution of sensory divergence in the context of limited gene flow in the bumblebee bat.

    Nat Commun 2011, 2:573. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Edelaar P, Alonso D, Lagerveld S, Senar JC, Björklund M: Population differentiation and restricted gene flow in Spanish crossbills: not isolation‒by‒distance but isolation‒by‒ecology.

    J Evol Biol 2012, 25:417-430. PubMed Abstract | Publisher Full Text OpenURL

  9. Antonovics J: Evolution in closely adjacent plant populations: VI: manifold effects of gene flow.

    Heredity 1968, 23:507-524. OpenURL

  10. Slatkin M: Gene flow and the geographic structure of natural populations.

    Science 1987, 236:787-792. PubMed Abstract | Publisher Full Text OpenURL

  11. Räsänen K, Hendry AP: Disentangling interactions between adaptive divergence and gene flow when ecology drives diversification.

    Ecol Lett 2008, 11:624-636. PubMed Abstract | Publisher Full Text OpenURL

  12. Niemiller ML, Fitzpatrick BM, Miller BT: Recent divergence with gene flow in Tennessee cave salamanders (Plethodontidae: Gyrinophilus) inferred from gene genealogies.

    Mol Ecol 2008, 17:2258-2275. PubMed Abstract | Publisher Full Text OpenURL

  13. Milá B, Wayne RK, Fitze P, Smith TB: Divergence with gene flow and fine‒scale phylogeographical structure in the wedge‒billed woodcreeper, Glyphorynchus spirurus, a Neotropical rainforest bird.

    Mol Ecol 2009, 18:2979-2995. PubMed Abstract | Publisher Full Text OpenURL

  14. Ribeiro ÂM, Lopes RJ, Bowie RC: Historical demographic dynamics underlying local adaptation in the presence of gene flow.

    Ecol Evol 2012, 2:2710-2721. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Muñoz MM, Crawford NG, McGreevy TJ, Messana NJ, Tarvin RD, Revell LJ, Zandvliet RM, Hopwood JM, Mock E, Schneider AL, Schneider CJ: Divergence in coloration and ecological speciation in the Anolis marmoratus species complex.

    Mol Ecol 2013, 22:2668-2682. PubMed Abstract | Publisher Full Text OpenURL

  16. Nosil P: The role of selection and gene flow in the evolution of sexual isolation in Timema walking sticks and other Orthopteroids.

    J Orthop Res 2005, 14:247-253. Publisher Full Text OpenURL

  17. Nosil P, Crespi BJ: Does gene flow constrain adaptive divergence or vice versa? A test using ecomorphology and sexual isolation in Timema cristinae walking‒sticks.

    Evolution 2004, 58:102-112. PubMed Abstract OpenURL

  18. Crispo E: Modifying effects of phenotypic plasticity on interactions among natural selection, adaptation and gene flow.

    J Evol Biol 2008, 21:1460-1469. PubMed Abstract | Publisher Full Text OpenURL

  19. Pfennig DW, Wund MA, Snell-Rood EC, Cruickshank T, Schlichting CD, Moczek AP: Phenotypic plasticity’s impacts on diversification and speciation.

    Trends Ecol Evol 2010, 25:459-467. PubMed Abstract | Publisher Full Text OpenURL

  20. Via S, Gomulkiewicz R, De Jong G, Scheiner SM, Schlichting CD, Van Tienderen PH: Adaptive phenotypic plasticity: consensus and controversy.

    Trends Ecol Evol 1995, 10:212-217. PubMed Abstract | Publisher Full Text OpenURL

  21. West-Eberhard MJ: Developmental plasticity and evolution. New York: Oxford University Press; 2003. OpenURL

  22. Richards CL, Bossdorf O, Pigliucci M: What role does heritable epigenetic variation play in phenotypic evolution?

    Bioscience 2010, 60:232-237. Publisher Full Text OpenURL

  23. Ghalambor CK, McKay JK, Carroll SP, Reznick DN: Adaptive versus non‒adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments.

    Funct Ecol 2007, 21:394-407. Publisher Full Text OpenURL

  24. Ord TJ, Stamps JA, Losos JB: Adaptation and plasticity of animal communication in fluctuating environments.

    Evolution 2010, 64:3134-3148. PubMed Abstract | Publisher Full Text OpenURL

  25. Richter‒Boix A, Teplitsky C, Rogell B, Laurila A: Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow.

    Mol Ecol 2010, 19:716-731. PubMed Abstract | Publisher Full Text OpenURL

  26. Thibert‒Plante X, Hendry AP: The consequences of phenotypic plasticity for ecological speciation.

    J Evol Biol 2011, 24:326-342. PubMed Abstract | Publisher Full Text OpenURL

  27. Lundsgaard‒Hansen B, Matthews B, Vonlanthen P, Taverna A, Seehausen O: Adaptive plasticity and genetic divergence in feeding efficiency during parallel adaptive radiation of whitefish (Coregonus spp.).

    J Evol Biol 2013, 26:483-498. PubMed Abstract | Publisher Full Text OpenURL

  28. Pelletier F, Reale D, Garant D, Coltman DW, Festa‒Bianchet M: Selection on heritable seasonal phenotypic plasticity of body mass.

    Evolution 2007, 61:1969-1979. PubMed Abstract | Publisher Full Text OpenURL

  29. Jourdan‒Pineau H, David P, Crochet PA: Phenotypic plasticity allows the Mediterranean parsley frog Pelodytes punctatus to exploit two temporal niches under continuous gene flow.

    Mol Ecol 2012, 21:876-886. PubMed Abstract | Publisher Full Text OpenURL

  30. Crispo E, Chapman LJ: Population genetic structure across dissolved oxygen regimes in an African cichlid fish.

    Mol Ecol 2008, 17:2134-2148. PubMed Abstract | Publisher Full Text OpenURL

  31. Ali MA: Sensory Ecology. New York: Plenum Press; 1978. OpenURL

  32. Boughman JW: How sensory drive can promote speciation.

    Trends Ecol Evol 2002, 17:571-577. Publisher Full Text OpenURL

  33. Dangles O, Irschick D, Chittka L, Casas J: Variability in sensory ecology: expanding the bridge between physiology and evolutionary biology.

    Q Rev Biol 2009, 84:51-74. PubMed Abstract | Publisher Full Text OpenURL

  34. Wilkins MR, Seddon N, Safran RJ: Evolutionary divergence in acoustic signals: causes and consequences.

    Trends Ecol Evol 2013, 28:156-166. PubMed Abstract | Publisher Full Text OpenURL

  35. Phillips CT, Johnston CE: Geographical divergence of acoustic signals in Cyprinella galactura the whitetail shiner (Cyprinidae).

    Anim Behav 2008, 75:617-626. Publisher Full Text OpenURL

  36. Pinto-Juma G, Simões PC, Seabra SG, Quartau JA: Calling song structure and geographic variation in Cicada orni Linnaeus (Hemiptera: Cicadidae).

    Zool Stud 2005, 44:81-94. OpenURL

  37. Lemmon EM: Diversification of conspecific signals in sympatry: geographic overlap drives multidimensional reproductive character displacement in frogs.

    Evolution 2009, 63:1155-1170. PubMed Abstract | Publisher Full Text OpenURL

  38. Podos J: Acoustic discrimination of sympatric morphs in Darwin’s finches: a behavioural mechanism for assortative mating?

    Phil Trans R Soc B 2010, 365:1031-1039. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Thinh VN, Hallam C, Roos C, Hammerschmidt K: Concordance between vocal and genetic diversity in crested gibbons.

    BMC Evol Biol 2011, 11:36. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  40. Beckers OM, Schul J: Developmental plasticity of mating calls enables acoustic communication in diverse environments.

    Proc R Soc B Biol Sci 2008, 275:1243-1248. Publisher Full Text OpenURL

  41. Ziegler L, Arim M, Narins PM: Linking amphibian call structure to the environment: the interplay between phenotypic flexibility and individual attributes.

    Behav Ecol 2011, 22:520-526. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Griffin DR: Echolocation in blind men, bats, and radar.

    Science 1944, 100:589-590. PubMed Abstract | Publisher Full Text OpenURL

  43. Griffin DR: Bat sounds under natural conditions, with evidence for echolocation of insect prey.

    J Exp Zool 1953, 123:435-465. Publisher Full Text OpenURL

  44. Knörnschild M, Jung K, Nagy M, Metz M, Kalko E: Bat echolocation calls facilitate social communication.

    Proc R Soc B 2012, 279:4827-4835. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Neuweiler G: Foraging, echolocation and audition in bats.

    Naturwissenschaften 1984, 71:446-455. Publisher Full Text OpenURL

  46. Hiryu S, Katsura K, Nagato T, Yamazaki H, Lin LK, Watanabe Y, Riquimaroux H: Intra-individual variation in the vocalized frequency of the Taiwanese leaf-nosed bat, Hipposideros terasensis, influenced by conspecific colony members.

    J Comp Physiol A 2006, 192:807-815. Publisher Full Text OpenURL

  47. Armstrong KN, Coles RB: Echolocation call frequency differences between geographic isolates of Rhinonicteris aurantia (Chiroptera: Hipposideridae): implications of nasal chamber size.

    J Mammal 2007, 88:94-104. Publisher Full Text OpenURL

  48. Heller K-G, Von Helversen O: Resource partitioning of sonar frequency bands in rhinolophoid bats.

    Oecologia 1989, 80:178-186. OpenURL

  49. Barclay RMR, Fullard JH, Jacobs DS: Variation in the echolocation calls of the hoary bat (Lasiurus cinereus): influence of body size, habitat structure, and geographic location.

    Can J Zool 1999, 77:530-534. Publisher Full Text OpenURL

  50. Wund MA: Variation in the echolocation calls of little brown bats (Myotis lucifugus) in response to different habitats.

    Am Midl Nat 2006, 156:99-108. Publisher Full Text OpenURL

  51. Houston RD, Boonman AM, Jones G: Do echolocation signal parameters restrict bats’ choice of prey? In Echolocation in Bats and Dolphins. Edited by Thomas JA, Moss CF, Vater M. Chicago: University of Chicago Press; 2004:339-345. OpenURL

  52. Jones G: Does echolocation constrain the evolution of body size in bats?

    Symp Zool Soc Lond 1996, 69:111-128. OpenURL

  53. Russo D, Jones G, Mucedda M: Influence of age, sex and body size on the echolocation calls of Mediterranean and Mehely’s horseshoe bats, Rhinolophus euryale and R. mehelyi (Chiroptera: Rhinolophidae).

    Mammalia 2001, 65:429-436. OpenURL

  54. Siemers BM, Beedholm K, Dietz C, Dietz I, Ivanova T: Is species identity, sex, age or individual quality conveyed by echolocation call frequency in European horseshoe bats?

    Acta Chiropterol 2005, 7:259-274. Publisher Full Text OpenURL

  55. Jones G, Teeling EC: The evolution of echolocation in bats.

    Trends Ecol Evol 2006, 21:149-156. PubMed Abstract | Publisher Full Text OpenURL

  56. Schuller G, Pollak G: Disproportionate frequency representation in the inferior colliculus of doppler-compensating greater horseshoe bats: evidence for an acoustic fovea.

    J Comp Physiol 1979, 132:47-54. Publisher Full Text OpenURL

  57. Schuller G, Beuter K, Schnitzler HU: Response to frequency shifted artificial echoes in the bat Rhinolophus ferrumequinum.

    J Comp Physiol A 1974, 89:275-286. Publisher Full Text OpenURL

  58. Schnitzler HU, Kalko EKV: Echolocation by insect-eating bats.

    Bioscience 2001, 51:557-569. Publisher Full Text OpenURL

  59. Jones G, Rayner JMV: Foraging behaviour and echolocation of wild horseshoe bats Rhinolophus ferrumequinum and R. hipposideros (Chiroptera, Rhinolophidae).

    Behav Ecol Sociobiol 1989, 25:183-191. Publisher Full Text OpenURL

  60. Schnitzler HU, Denzinger A: Auditory fovea and Doppler shift compensation: adaptations for flutter detection in echolocating bats using CF-FM signals.

    J Comp Physiol A 2011, 197:541-559. Publisher Full Text OpenURL

  61. Rübsamen R: Ontogenesis of the echolocation system in the Rufous horseshoe bat, Rhinolophus rouxi (audition and vocalisation in early postnatal development).

    J Comp Physiol A 1987, 161:899-913. PubMed Abstract | Publisher Full Text OpenURL

  62. Matsumura S: Mother-infant communication in the horseshoe bat (Rhinolophus ferrumequinum nippon): development of vocalization.

    J Mammal 1979, 60:76-84. Publisher Full Text OpenURL

  63. Jones G, Ransome RD: Echolocation calls of bats are influenced by maternal effects and change over a lifetime.

    Proc R Soc Lond B 1993, 252:125-128. Publisher Full Text OpenURL

  64. Boughman JW, Moss CF: Social sounds: vocal learning and development of mammal and bird calls. In Acoustic Communication. Edited by Simmons AM, Popper AN, Fay RR. New York: Springer-Verlag; 2003:138-224. OpenURL

  65. Kingston T, Lara MC, Jones G, Akbar Z, Kunz TH, Schneider CJ: Acoustic divergence in two cryptic Hipposideros species: a role for social selection?

    Proc R Soc Lond B 2001, 268:1381-1386. Publisher Full Text OpenURL

  66. Yoshino H, Armstrong KN, Izawa M, Yokoyama J, Kawata M: Genetic and acoustic population structuring in the Okinawa least horseshoe bat: are intercolony acoustic differences maintained by vertical maternal transmission?

    Mol Ecol 2008, 17:4978-4991. PubMed Abstract | Publisher Full Text OpenURL

  67. Jacobs DS, Barclay RMR, Walker MH: The allometry of echolocation call frequencies of insectivorous bats: why do some species deviate from the pattern?

    Oecologia 2007, 152:583-594. PubMed Abstract | Publisher Full Text OpenURL

  68. Russo D, Mucedda M, Bello M, Biscardi S, Pidinchedda E, Jones G: Divergent echolocation call frequencies in insular rhinolophids (Chiroptera): a case of character displacement?

    J Biogeogr 2007, 34:2129-2138. Publisher Full Text OpenURL

  69. Chen SF, Jones G, Rossiter SJ: Determinants of echolocation call frequency variation in the Formosan lesser horseshoe bat (Rhinolophus monoceros).

    Proc R Soc B 2009, 276:3901-3909. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  70. Kingston T, Rossiter SJ: Harmonic-hopping in Wallacea’s bats.

    Nature 2004, 429:654-657. PubMed Abstract | Publisher Full Text OpenURL

  71. Guillén A, Juste BJ, Ibáñez C: Variation in the frequency of echolocation calls of Hipposideros ruber in the Gulf of Guinea: an exploration of the adaptive meaning of the constant frequency value in Rhinolophoid CF bats.

    J Evol Biol 2000, 13:70-80. Publisher Full Text OpenURL

  72. Jiang T, Metzner W, You Y, Liu S, Lu G, Li S, Feng J: Variation in the resting frequency of Rhinolophus pusillus in Mainland China: effect of climate and implications for conservation.

    J Acoust Soc Am 2010, 128:2204-2211. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  73. Odendaal LJ, Jacobs DS: Morphological correlates of echolocation frequency in the endemic Cape horseshoe bat, Rhinolophus capensis (Chiroptera: Rhinolophidae).

    J Comp Physiol A 2011, 197:435-446. Publisher Full Text OpenURL

  74. Francis CM, Habersetzer J: Interspecific and intraspecific variation in echolocation call frequency and morphology of horseshoe bats, Rhinolophus and Hipposideros. In Bat Biology and Conservation. Edited by Kunz TH, Racey PA. Washington DC: Smithsonian Institute Press; 1998:169-179. OpenURL

  75. Duellman WE, Pyles RA: Acoustic resource partitioning in anuran communities.

    Copeia 1983, 1983:639-649. Publisher Full Text OpenURL

  76. Jones G, Barlow KE: Cryptic species of echolocating bats. In Echolocation in Bats and Dolphins. Edited by Thomas JA, Moss CF, Vater M. Chicago: University of Chicago Press; 2004:345-349. OpenURL

  77. Clare EL, Adams AM, Maya-Simões AZ, Eger JL, Hebert PD, Fenton MB: Diversification and reproductive isolation: cryptic species in the only New World high duty cycle bat: Pteronotus parnellii.

    BMC Evol Biol 2013, 13:26. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  78. Thabah A, Rossiter SJ, Kingston T, Zhang S, Parsons S, Mya KM, Jones G: Genetic divergence and echolocation call frequency in cryptic species of Hipposideros larvatus sl (Chiroptera: Hipposideridae) from the Indo‒Malayan region.

    Biol J Linn Soc 2006, 88:119-130. Publisher Full Text OpenURL

  79. Liu Y, Feng J, Metzner W: Different auditory feedback control for echolocation and communication in horseshoe bats.

    PloS One 2013, 8:e62710. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  80. Furusawa Y, Hiryu S, Kobayasi KI, Riquimaroux H: Convergence of reference frequencies by multiple CF–FM bats (Rhinolophus ferrumequinum nippon) during paired flights evaluated with onboard microphones.

    J Comp Physiol A 2012, 198:683-693. Publisher Full Text OpenURL

  81. Hage SR, Jiang T, Berquist SW, Feng J, Metzner W: Ambient noise induces independent shifts in call frequency and amplitude within the Lombard effect in echolocating bats.

    PNAS 2013, 110:4063-4068. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  82. Taylor PJ, Stoffberg S, Monadjem A, Schoeman MC, Bayliss J, Cotterill FP: Four new bat species (Rhinolophus hildebrandtii complex) reflect Plio-Pleistocene divergence of dwarfs and giants across an Afromontane archipelago.

    PloS One 2012, 7:e41744. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  83. Stoffberg S, Schoeman MC, Matthee CA: Correlated genetic and ecological diversification in a widespread southern African horseshoe bat.

    PloS One 2012, 7:e31946. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  84. Rutherford MC, Mucina L, Powrie LW: Biomes and bioregions of southern Africa. In The Vegetation of South Africa, Lesotho and Swaziland: Strelitzia 19. Edited by Mucina L, Rutherford MC. Pretoria: South African National Biodiversity Institute; 2006:30-51. OpenURL

  85. Linder HP, Johnson SD, Kuhlmann M, Matthee CA, Nyffeler R, Swartz E: Biotic diversity in the Southern African winter-rainfall region.

    Curr Opin Environ Sustain 2010, 2:109-116. Publisher Full Text OpenURL

  86. Anthony EL: Age determination in bats. In Ecological and Behavioural Methods for the Study of Bats. Edited by Kunz TH. Washington DC: Smithsonian Institute Press; 1988:47-58. OpenURL

  87. Jacobs DS, Babiker H, Bastian A, Kearney T, van Eeden R, Bishop JM: Phenotypic convergence in genetically distinct lineages of a Rhinolophus species complex (Mammalia, Chiroptera).

    PloS One 2013, 8:e82614. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  88. Neuweiler G: Foraging ecology and audition in echolocating bats.

    Trends Ecol Evol 1989, 4:160-166. PubMed Abstract | Publisher Full Text OpenURL

  89. Worthington WJ, Barratt E: A non-lethal method of tissue sampling for genetic studies of chiropterans.

    Bat Res News 1996, 37:1-4. OpenURL

  90. Hoelzel AR, Hancock JM, Dover GA: Evolution of the cetacean mitochondrial d-loop region.

    Mol Biol Evol 1991, 8:475-493. PubMed Abstract | Publisher Full Text OpenURL

  91. Wilkinson GS, Chapman AM: Length and sequence variation in evening bat D-loop mtDNA.

    Genetics 1991, 128:607-617. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  92. Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT.

    Nucleic Acids Symp Ser 1999, 41:95-98. OpenURL

  93. Nei M: Molecular evolutionary genetics. New York: Columbia University Press; 1987. OpenURL

  94. Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data.

    Bioinformatics 2009, 25:1451-1452. PubMed Abstract | Publisher Full Text OpenURL

  95. McCulloch CE, Searle SR: Generalized, linear and mixed models. New York: Willey; 2000. OpenURL

  96. Freckleton RP: Dealing with collinearity in behavioural and ecological data: model averaging and the problems of measurement error.

    Behav Ecol Sociobiol 2011, 65:91-101. Publisher Full Text OpenURL

  97. Akaike H: Information theory and an extension of the maximum likelihood principle. In Second International Symposium on Information Theory. Edited by Petrov BN, Czàke F. Budapest: Akademiai Kiadó; 1973:267-281. OpenURL

  98. Burnham KP, Anderson DR: Model selection and multi-model inference: a practical information-theoretic approach. New York: Springer-Verlag; 2002. OpenURL

  99. Burnham KP, Anderson DR: Multimodel inference: understanding AIC and BIC in model selection.

    Sociol Methods Res 2004, 33:261-304. Publisher Full Text OpenURL

  100. Turner W, Spector S, Gardiner N, Fladeland M, Sterling E, Steininger M: Remote sensing for biodiversity science and conservation.

    Trends Ecol Evol 2003, 18:306-314. Publisher Full Text OpenURL

  101. Myneni R, Hall F, Sellers P, Marshak A: The interpretation of spectral vegetation indexes.

    IEEE Trans Geosci Rem Sens 1995, 33:481-486. Publisher Full Text OpenURL

  102. Carlson TN, Ripley DA: On the relation between NDVI, fractional vegetation cover, and leaf area index.

    Remote Sens Environ 1997, 62:241-252. Publisher Full Text OpenURL

  103. Purevdorj TS, Tateishi R, Ishiyama T, Honda Y: Relationships between percent vegetation cover and vegetation indices.

    Int J Remote Sens 1998, 19:3519-3535. Publisher Full Text OpenURL

  104. Borowik T, Pettorelli N, Sönnichsen L, Jędrzejewska B: Normalized difference vegetation index (NDVI) as a predictor of forage availability for ungulates in forest and field habitats.

    Eur J Wildl Res 2013, 59:675-682. Publisher Full Text OpenURL

  105. Pettorelli N, Vik JO, Mysterud A, Gaillard JM, Tucker CJ, Stenseth NC: Using the satellite-derived normalized difference vegetation index (NDVI) to assess ecological effects of environmental change.

    Trends Ecol Evol 2005, 20:503-510. PubMed Abstract | Publisher Full Text OpenURL

  106. Pettorelli N, Ryan S, Mueller T, Bunnefeld N, Jedrzejewska B, Lima M, Kausrud K: The normalized difference vegetation index (NDVI): unforeseen successes in animal ecology.

    Clim Res 2011, 46:5-27. OpenURL

  107. Neigh CSR, Tucker CJ, Townshend JRG: North American vegetation dynamics observed with multi-resolution satellite data.

    Remote Sens Environ 2008, 112:1749-1772. Publisher Full Text OpenURL

  108. Armstrong KN, Kerry LJ: Modelling the prey detection performance of Rhinonicteris aurantia (Chiroptera: Hipposideridae) in different atmospheric conditions discounts the notional role of relative humidity in adaptive evolution.

    J Theor Biol 2011, 278:44-54. PubMed Abstract | Publisher Full Text OpenURL

  109. Csorba G, Ujhelyi P, Thomas N: Horseshoe bats of the world: (Chiroptera: Rhinolophidae). Shropshire: Alana Books; 2003. OpenURL

  110. Stilz WP, Schnitzler HU: Estimation of the acoustic range of bat echolocation for extended targets.

    J Acoust Soc Am 2012, 132:1765. PubMed Abstract | Publisher Full Text OpenURL

  111. Long GR, Schnitzler HU: Behavioural audiograms from the bat, Rhinolophus ferrumequinum.

    J Comp Physiol A 1975, 100:211-219. Publisher Full Text OpenURL

  112. Meyer CF, Kalko EK, Kerth G: Small‒scale fragmentation effects on local genetic diversity in two Phyllostomid bats with different dispersal abilities in Panama.

    Biotropica 2009, 41:95-102. Publisher Full Text OpenURL

  113. Petren K, Grant PR, Grant BR, Keller LF: Comparative landscape genetics and the adaptive radiation of Darwin’s finches: the role of peripheral isolation.

    Mol Ecol 2005, 14:2943-2957. PubMed Abstract | Publisher Full Text OpenURL

  114. Huson DH, Bryant D: Application of phylogenetic networks in evolutionary studies.

    Mol Biol Evol 2006, 23:254-267. PubMed Abstract | Publisher Full Text OpenURL

  115. Cassens I, Van Waerebeek K, Best PB, Crespo EA, Reyes J, Milinkovitch MC: The phylogeography of dusky dolphins (Lagenorhynchus obscurus): a critical examination of network methods and rooting procedures.

    Mol Ecol 2003, 12:1781-1792. PubMed Abstract | Publisher Full Text OpenURL

  116. Cheng L, Connor TR, Sirén J, Aanensen DM, Corander J: Hierarchical and spatially explicit clustering of DNA sequences with BAPS software.

    Mol Biol Evol 2013, 30:1224-1228. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  117. Peakall R, Smouse PE: GenAlEx 6.5: genetic analysis in excel: population genetic software for teaching and research—an update.

    Bioinformatics 2012, 28:2537-2539. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  118. Wright S: The interpretation of population structure by F-statistics with special regard to systems of mating.

    Evolution 1965, 19:395-420. Publisher Full Text OpenURL

  119. Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data.

    Genetics 1992, 131:479-491. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  120. Beerli P: How to use MIGRATE or why are Markov Chain Monte Carlo programs difficult to use? In Population Genetics for Animal Conservation. Edited by Bertorelle G, Bruford MW, Hauffe HC, Rizzoli A, Vernesi C. Cambridge: Cambridge University Press; 2009:42-79. OpenURL

  121. Hutterer R, Ivanova T, Meyer-Cords C, Rodrigues L: Bat migration in Europe: a review of banding data and literature. Bonn: Federal Agency for Nature Conservation; 2005. OpenURL

  122. Beerli P, Felsenstein J: Maximum-likelihood estimation of migration rates and effective population numbers in two populations using a coalescent approach.

    Genetics 1999, 152:763-773. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  123. Rudh A, Rogell B, Höglund J: Non‒gradual variation in colour morphs of the strawberry poison frog Dendrobates pumilio: genetic and geographical isolation suggest a role for selection in maintaining polymorphism.

    Mol Ecol 2007, 16:4284-4294. PubMed Abstract | Publisher Full Text OpenURL

  124. Mantel N: The detection of disease clustering and a generalized regression approach.

    Cancer Res 1967, 27:209-220. PubMed Abstract | Publisher Full Text OpenURL

  125. Smouse PE, Long JC, Sokal RR: Multiple regression and correlation extensions of the Mantel test of matrix correspondence.

    Syst Zool 1986, 35:627-632. Publisher Full Text OpenURL

  126. Slatkin M: Isolation by distance in equilibrium and non-equilibrium populations.

    Evolution 1993, 47:264-279. Publisher Full Text OpenURL

  127. Bryant D, Moulton V: Neighbor-net: an agglomerative method for the construction of phylogenetic networks.

    Mol Biol Evol 2004, 21:255-265. PubMed Abstract | Publisher Full Text OpenURL

  128. Turmelle AS, Kunz TH, Sorenson MD: A tale of two genomes: contrasting patterns of phylogeographic structure in a widely distributed bat.

    Mol Ecol 2011, 20:357-375. PubMed Abstract | Publisher Full Text OpenURL

  129. Toews DP, Brelsford A: The biogeography of mitochondrial and nuclear discordance in animals.

    Mol Ecol 2012, 21:3907-3930. PubMed Abstract | Publisher Full Text OpenURL

  130. Mao XG, Zhu GJ, Zhang S, Rossiter SJ: Pleistocene climatic cycling drives intra‒specific diversification in the intermediate horseshoe bat (Rhinolophus affinis) in Southern China.

    Mol Ecol 2010, 19:2754-2769. PubMed Abstract | Publisher Full Text OpenURL

  131. Mao X, He G, Hua P, Jones G, Zhang S, Rossiter SJ: Historical introgression and the persistence of ghost alleles in the intermediate horseshoe bat (Rhinolophus affinis).

    Mol Ecol 2013, 22:1035-1050. PubMed Abstract | Publisher Full Text OpenURL

  132. Dool SE, Puechmaille SJ, Dietz C, Juste J, Ibáñez C, Hulva P, Roué SG, Petit EJ, Jones G, Russo D, Toffoli R, Viglino A, Martinoli A, Rossiter SJ, Teeling EC: Phylogeography and postglacial recolonisation of Europe by Rhinolophus hipposideros: evidence from multiple genetic markers.

    Mol Ecol 2013, 22:4055-4070. PubMed Abstract | Publisher Full Text OpenURL

  133. Chen SF, Jones G, Rossiter SJ: Sex‒biased gene flow and colonization in the Formosan lesser horseshoe bat: inference from nuclear and mitochondrial markers.

    J Zool 2008, 274:207-215. Publisher Full Text OpenURL

  134. Stone GN, Nee S, Felsenstein J: Controlling for non-independence in comparative analysis of patterns across populations within species.

    Phil Trans R Soc B 2011, 366:1410-1424. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  135. Felsenstein J: Contrasts for a within-species comparative method. In Modern developments in theoretical population genetics. Edited by Slatkin M, Veuille M. Oxford: Oxford University Press; 2002:118-129. OpenURL

  136. Nosil P: Speciation with gene flow could be common.

    Mol Ecol 2008, 17:2103-2106. PubMed Abstract | Publisher Full Text OpenURL

  137. Chattopadhyay B, Garg KM, Ramakrishnan U, Kandula S: Sibling species in South Indian populations of the rufous horse-shoe bat Rhinolophus rouxii.

    Conserv Genet 2012, 13:1435-1445. Publisher Full Text OpenURL

  138. Fitzpatrick BM: Underappreciated consequences of phenotypic plasticity for ecological speciation.

    Int J Ecol 2012, 2012:256017. OpenURL

  139. Genevois F, Bretagnolle V: Male blue petrels reveal body mass when calling.

    Ethol Ecol Evol 1994, 6:377-383. Publisher Full Text OpenURL

  140. Castellano S, Cuatt B, Rinella R, Rosso A, Giacoma C: The advertisement call of European treefrogs (Hyla arborea): a multi-level study of variation.

    Ethology 2002, 108:75-89. Publisher Full Text OpenURL

  141. Brown WD, Wideman J, Andrade MCB, Mason AC, Gwynne DT: Female choice for an indicator of male size in the song of the black-horned tree cricket, Oecanthus nigricornis (Orthoptera: Gryllidae: Oecanthinae).

    Evolution 1996, 50:2400-2411. Publisher Full Text OpenURL

  142. Smith TB, Wayne RK, Girman DJ, Bruford MW: A role for ecotones in generating rainforest biodiversity.

    Science 1997, 276:1855-1857. Publisher Full Text OpenURL

  143. Tolley KA, Chase BM, Forest F: Speciation and radiations track climate transitions since the miocene climatic optimum: a case study of southern African chameleons.

    J Biogeogr 2008, 35:1402-1414. Publisher Full Text OpenURL

  144. Daniels SR, Gouws G, Crandall KA: Phylogeographic patterning in a freshwater crab species (Decapoda: Potamonautidae: Potamonautes) reveals the signature of historical climatic oscillations.

    J Biogeogr 2006, 33:1538-1549. Publisher Full Text OpenURL

  145. du Toit ND, Van Vuuren BJ, Matthee S, Matthee CA: Biome specificity of distinct genetic lineages within the four-striped mouse Rhabdomys pumilio (Rodentia: Muridae) from southern Africa with implications for taxonomy.

    Mol Phylogenet Evol 2012, 65:75-86. PubMed Abstract | Publisher Full Text OpenURL

  146. Willows‒Munro S, Matthee CA: Linking lineage diversification to climate and habitat heterogeneity: phylogeography of the southern African shrew Myosorex varius.

    J Biogeogr 2011, 38:1976-1991. Publisher Full Text OpenURL

  147. Chase BM, Meadows ME: Late Quaternary dynamics of southern Africa’s winter rainfall zone.

    Earth-Sci Rev 2007, 84:103-138. Publisher Full Text OpenURL

  148. Xu Z, Jing W, Keping S, Tinglei J, Yunlei J, Jiang F: Echolocation calls of Rhinolophus ferrumequinum in relation to habitat type and environmental factors.

    Acta Ecol Sin 2008, 28:5248-5258. Publisher Full Text OpenURL

  149. Goiti U, Aihartza J, Garin I, Zabala J: Influence of habitat on the foraging behaviour of the Mediterranean horseshoe bat, Rhinolophus euryale.

    Acta Chiropterol 2003, 5:75-84. Publisher Full Text OpenURL

  150. Goiti U, Garin I, Almenar D, Salsamendi E, Aihartza J: Foraging by Mediterranean horseshoe bats (Rhinolophus euryale) in relation to prey distribution and edge habitat.

    J Mammal 2008, 89:493-502. Publisher Full Text OpenURL

  151. Salsamendi E, Arostegui I, Aihartza J, Almenar D, Goiti U, Garin I: Foraging ecology in Mehely’s horseshoe bats: influence of habitat structure and water availability.

    Acta Chiropterol 2012, 14:121-132. Publisher Full Text OpenURL

  152. Lee YF, Kuo YM, Chu WC, Lin YH, Chang HY, Chen WM: Ecomorphology, differentiated habitat use, and nocturnal activities of Rhinolophus and Hipposideros species in East Asian tropical forests.

    Zoology 2012, 115:22-29. PubMed Abstract | Publisher Full Text OpenURL

  153. Russo D, Almenar D, Aihartza J, Goiti U, Salsamendi E, Garin I: Habitat selection in sympatric Rhinolophus mehelyi and R. euryale (Mammalia: Chiroptera).

    J Zool (Lond) 2005, 266:327-332. Publisher Full Text OpenURL

  154. Surlykke A, Kalko EK: Echolocating bats cry out loud to detect their prey.

    PLoS One 2008, 3:e2036. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  155. Rice AM, Pfennig DW: Does character displacement initiate speciation? Evidence of reduced gene flow between populations experiencing divergent selection.

    J Evol Biol 2010, 23:854-865. PubMed Abstract | Publisher Full Text OpenURL

  156. Schuchmann M, Puechmaille SJ, Siemers BM: Horseshoe bats recognise the sex of conspecifics from their echolocation calls.

    Acta Chiropterol 2012, 14:161-166. Publisher Full Text OpenURL

  157. Li Y, Liu Z, Shi P, Zhang J: The hearing gene Prestin unites echolocating bats and whales.

    Curr Biol 2010, 20:R55-R56. PubMed Abstract | Publisher Full Text OpenURL

  158. Davies KTJ, Cotton JA, Kirwan JD, Teeling EC, Rossiter SJ: Parallel signatures of sequence evolution among hearing genes in echolocating mammals: an emerging model of genetic convergence.

    Heredity 2011, 108:480-489. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL