- Open Access
Evolutionary origins of ultrasonic hearing and laryngeal echolocation in bats inferred from morphological analyses of the inner ear
Frontiers in Zoology volume 10, Article number: 2 (2013)
Many mammals have evolved highly adapted hearing associated with ecological specialisation. Of these, bats possess the widest frequency range of vocalisations and associated hearing sensitivities, with frequencies of above 200 kHz in some lineages that use laryngeal echolocation. High frequency hearing in bats appears to have evolved via structural modifications of the inner ear, however, studying these minute features presents considerable challenges and hitherto few such attempts have been made. To understand these adaptations more fully, as well as gain insights into the evolutionary origins of ultrasonic hearing and echolocation in bats, we undertook micro-computed tomography (μCT) scans of the cochleae of representative bat species from 16 families, encompassing their broad range of ecological diversity. To characterise cochlear gross morphology, we measured the relative basilar membrane length and number of turns, and compared these values between echolocating and non-echolocating bats, as well as other mammals.
We found that hearing and echolocation call frequencies in bats correlated with both measures of cochlear morphology. In particular, relative basilar membrane length was typically longer in echolocating species, and also correlated positively with the number of cochlear turns. Ancestral reconstructions of these parameters suggested that the common ancestor of all extant bats was probably capable of ultrasonic hearing; however, we also found evidence of a significant decrease in the rate of morphological evolution of the basilar membrane in multiple ancestral branches within the Yangochiroptera suborder. Within the echolocating Yinpterochiroptera, there was some evidence of an increase in the rate of basilar membrane evolution in some tips of the tree, possibly associated with reported shifts in call frequency associated with recent speciation events.
The two main groups of echolocating bat were found to display highly variable inner ear morphologies. Ancestral reconstructions and rate shift analyses of ear morphology point to a complex evolutionary history, with the former supporting ultrasonic hearing in the common bat ancestor but the latter suggesting that morphological changes associated with echolocation might have occurred later. These findings are consistent with theories that sophisticated laryngeal echolocation, as seen in modern lineages, evolved following the divergence of the two main suborders.
The evolutionary success of mammals can in part be attributed to remarkable sensory diversification e.g.. While some lineages show evidence of changes across multiple sensory modalities e.g. others have evolved one highly specialised system e.g.. Mammalian auditory systems are particularly well-developed compared to those of many other vertebrate groups, furthermore, mammalian hearing can be characterised by high sensitivity and selectivity and in particular by broad frequency ranges with high upper frequency limits[4, 5]. Multidisciplinary evidence suggests that these auditory adaptations in mammals e.g.[6, 7] can be linked to three principal adaptations: the evolution of three ossicles in the middle ears (malleus, incus and stapes), elongation of the basilar membrane in the cochlea which provides a supportive base for the sensory hair cells and the evolutionary innovation of outer hair cells (OHC)[8–10]. However, genetic studies suggest additional molecular changes have occurred in the motor protein of the OHC, known as Prestin, that are probably related to the acquisition of high frequency hearing seen in extant therian mammals. Therefore, it has recently been summarized that true high frequency hearing (i.e. >20 kHz) likely evolved approximately 125 million years (Ma) ago within the therian lineage and required additional structural modifications of the inner ear (as reviewed in).
The mammalian cochlea is a coiled cavity in which high and low frequency sounds are perceived by the basal and upper turns respectively, and this tonotopic organisation is partly achieved by a decrease in basilar membrane stiffness from base to apex[14, 15]. Mammals exhibit considerable variation in hearing capabilities and cochlear morphology, although some consistent relationships link these two traits. Basilar membrane length is positively correlated with body mass, and absolute basilar membrane length is negatively correlated with both high and low frequency hearing limits e.g.[18, 19]. It has been hypothesised that cochlear coiling evolved in response to selection pressures relating to the accommodation of elongated auditory sensory membranes of the inner ear, however, a previous study found no significant relationship between the number of cochlear turns and basilar membrane length. Recent evidence suggests that coiled cochleae may play a mechanical role in low frequency hearing limit. Across auditory ‘generalists’ high frequency hearing limits correlate with inter-aural distance, whereas auditory specialists such as subterranean mole rats and echolocating cetaceans deviate from this relationship. Other factors determining the morphology of auditory systems include physical and mechanical constraints, as well as phylogeny e.g.[16, 22, 23].
Although several mammalian taxa, including some rodents, carnivores and primates are capable of either detecting or producing ultrasonic sounds (>20 kHz)[4, 24] the most highly developed auditory systems for perceiving ultrasonic sound are seen in toothed whales and laryngeal echolocating bats[16, 25]. Bats possess some of the widest frequency ranges of vocalisations and, therefore, assumed associated hearing sensitivities of any mammal group, with recorded vocalisations ranging from below 20 kHz to over 200 kHz across the order (as reviewed in[26, 27]). Of 19 currently recognised bat families, all but one (the Old World fruit bats) use laryngeal echolocation for orientation, obstacle avoidance and, in most taxa, prey detection. The inner ears of laryngeal echolocating bats show several structural adaptations for detecting ultrasonic echoes; in particular, their cochleae are often enlarged and contain 2.5 to 3.5 turns compared to only an average of 1.75 in non-echolocating fruit bats[29–31]. Furthermore, different forms of echolocation appear to have led to different and sometimes convergent inner ear adaptations[32, 33]. For example, both Old World horseshoe bats and the New World moustached bat possess a greatly enlarged cochlear basal turn, which allows exquisite tuning of the inner ear to the echoes of the specialised constant frequency (CF) calls produced by these taxa[34, 35]. Previous studies suggest specific adaptations of the anchoring system and the width of the basilar membrane in echolocating bats.
Well-supported phylogenies of bats show that laryngeal echolocation is distributed across two highly divergent suborders of bats, termed Yinpterochiroptera and Yangochiroptera, the former of which also contains the non-echolocating Old World fruit bats. To account for this pattern, two main evolutionary scenarios have been proposed; first, that echolocation evolved once in the common ancestor of all modern bats with subsequent loss in Old World fruit bats, and two, that laryngeal echolocation evolved multiple times across the order[38, 39]. While fossils bats from the early Eocene have been taxonomically classified as falling outside of the modern bats, and thus might not inform this issue[40, 41], recent reports that echolocating members of the two suborders have undergone convergent amino acid replacements in several ‘hearing genes’ would appear to support the multiple origin hypothesis see[42–45].
Evolutionary modifications of the inner ear are likely to have arisen from selection acting on numerous loci, and thus a comparative analysis of morphology might offer a powerful means of reconstructing the origins of ultrasonic hearing and thus laryngeal echolocation in bats. Here we reconstructed three-dimensional bat inner ear volumes of a range of bat species from 16 families, encompassing the broad diversity of echolocation call types and ecological traits seen in the order. We compared patterns of cochlear morphological variation - as defined by relative basilar membrane length and number of turns - among bats, and also between bats and other mammals for which data were available. Correlations were investigated between high and low hearing frequency limits and also echolocation call parameters and morphological characters of the cochlea. We predicted that echolocating bat species would show specific adaptations in aspects of cochlear gross morphology compared to both non-echolocating bats and other mammal species, due to the particular demands associated with receiving the high frequency sounds produced during laryngeal echolocation.
If significant inner ear adaptations for ultrasonic hearing do occur, then these patterns might help us to distinguish between the two most parsimonious alternate scenarios of the evolution of echolocation. For example, it may be possible to determine whether Old World fruit bats show evidence of a loss of echolocation. Alternatively there may be evidence of functional adaptation in the inner ears of the two groups of echolocating bats. Specifically, if Old World fruit bats have lost echolocation, then we might expect to find signatures of this in their cochleae such as intermediate forms between those of echolocating bats and non-echolocating mammals or increased morphological variability, consistent with morphological relaxation.
The two alternate evolutionary scenarios that have been proposed to account for the distribution of laryngeal echolocation across divergent clades of bats might also be expected to have left different traces of inner ear morphological evolution across the bat phylogenetic tree. Firstly, given a single origin in the bat common ancestor we might expect to see increased rates of morphological change in the ancestral bat branch, possibly coupled with a subsequent rate shift in the non-echolocating Old World fruit bats corresponding to a loss of structures associated with sophisticated echolocation capability. In contrast, multiple origins might be expected to leave a signal of morphological rate shifts on specific branches after the point that the two main bat suborders diverged.
Basilar membrane length
A plot of log basilar membrane length (mm) versus log body mass0.33 revealed that a substantial proportion of laryngeal echolocating bats fell above the upper 95% prediction intervals (PI) of the regression line based on non-echolocating placental mammals (excluding baleen whales and Old World fruit bats) (see Figure 1). Most of these outliers were bat species that use constant frequency (CF) echolocation, with the exceptions of Taphozous peli and Cheiromeles torquatus. Just one laryngeal echolocating bat - Macroderma gigas - fell below the lower 95% PI. In contrast to the laryngeal echolocating bats, all non-echolocating Old World fruit bats fell close to the placental mammal regression line and within the PIs. Of the cetacean species plotted, all toothed whales fell within the PI, however, several baleen whales fell outside of the lower 95% PI, and only one placental mammal - the Californian sea lion, Zalophus californianus - fell above the upper 95% PI. Of the two marsupial species, one, Didelphis virginiana, fell within the placental mammal distribution and the other, Monodelphis domestica, just below. Values for the platypus, Ornithorhynchus anatinus, and the echidna, Tachyglossus aculeatus, clearly fell below the placental mammal prediction interval consistent with them having shorter relative basilar membrane lengths. However, the small sample size of both monotremes and marsupials (n = 2 in each case) should be noted.
Bayesian phylogenetic mixed model comparisons revealed no consistent difference in the basilar membrane length of either laryngeal echolocating bats or CF echolocating bats compared to those of other mammals. Across all taxa sampled laryngeal echolocation was not a significant factor (P MCMC = 0.059) and did not result in improved model fit (ΔDIC = −0.852). Although CF echolocation was a significant factor (P MCMC = 0.016), it too did not result in an improved model fit (ΔDIC = −0.602). Similar results were obtained in an analysis with a reduced dataset, in which only non-whale placental mammals were included (see supplementary results Additional file1: Table S5). Baleen and toothed cetaceans were omitted due to their specialised low- and high-frequency hearing respectively, and also due to their much larger body mass.
Number of spiral turns
To test the hypothesis that increased cochlear coiling evolved in response to selection pressures relating to the accommodation of elongated auditory sensory membranes, correlations were investigated between numbers of turns and relative basilar membrane lengths (Additional file1: Table S2 and Figures 1 and2). Since cochlear coiling might be predicted to be more important in smaller bodied taxa, the relative membrane length was calculated by dividing membrane length by body mass0.33. Across all taxa sampled, a significant positive correlation was found between log number of turns and log relative membrane length (log turns = 0.123 log relative membrane + 0.342, RSE = 0.122 (98 d.f.), R = 0.45, adjusted R2 = 0.19, F = 24.9 (1, 98 d.f., P = 2.68 × 10-6), which remained significant after accounting for phylogeny (log turns = 0.099 log relative membrane + 0.246, DIC: -321.79, P MCMC (Intercept) < 5 × 10-4, P MCMC (log relative membrane) = 0.024).
Echolocation call analysis
Step-wise multiple regression analyses were used to assess whether echolocation call parameters (i.e. minimum frequency, maximum frequency and peak frequency) relate to both cochlear morphology (basilar membrane length and number of cochlear turns) and body mass across 62 laryngeal echolocating bats and one species of Old World fruit bat that uses tongue-clicking as a rudimentary form of echolocation (Rousettus aegyptiacus). After accounting for phylogeny, basilar membrane length and number of cochlear turns both showed significant relationships with all three echolocation call parameters, while body mass did not (Additional file2: Figure S2 for plots and multiple regression statistics). Furthermore, the two cochlear variables were both significant when fitted together, improving fit (ΔDIC) over univariate models. For minimum echolocation frequency, ΔDIC = 9.48 (Min. frequency = −0.89 log membrane + 2.25 log number of turns + 1.66, DIC: -65.98, P MCMC (log membrane) <1 × 10-4, P MCMC (log turns) <1 × 10-4. P MCMC (intercept) <1 × 10-4). For peak echolocation frequency, ΔDIC = 4.83, (Peak frequency = −0.88 log membrane + 2.00 log number of turns + 1.88, DIC: -53.44, P MCMC (log membrane) <1 × 10-4, P MCMC (log turns) < 5 × 10-4, P MCMC (intercept) <1 × 10-4). For maximum echolocation frequency, ΔDIC = 4.89, (Max. frequency = −0.78 log membrane + 1.44 log number of turns + 2.09, DIC: -50.31, P MCMC (log membrane) = 8.90 × 10-4, P MCMC (log turns) = 0.01, PMCM C(intercept) <1 × 10-4). Full details are given in Additional file1: Table S6.
Ancestral reconstructions of inner ear morphology, hearing and echolocation
Phylogenetic independent contrast and maximum likelihood ancestral reconstructions of both log relative membrane length and log number of cochlear turns are shown in Figure 3 and Additional file3: Figure S3, respectively. Both methods indicated that the log relative basilar membrane length of the hypothetical common ancestor of modern bats was 0.64 mm/g0.33, which corresponds to a basilar membrane length of 11.75 mm based on a mass of 19.22 g. Similarly, both methods gave an approximately equal number of cochlear turns: 2.45 (anti-logged). These values suggest a high frequency hearing limit of ~100 kHz at 60 dB and ~60 kHz at 30 dB, with low frequency hearing limit of ~1 kHz at either 30 or 60 dB (see Figure 4).
Based on our ancestral reconstructions of basilar membrane length and cochlear turns, it is not possible to predict whether or not the common ancestral bat was capable of sophisticated laryngeal echolocation. However, under the proposed scenario of a single origin of echolocation, then the observed relationship between inner ear characters and echolocation call frequencies in modern echolocating bats would suggest that the hypothetical call frequencies of the ancestor would have been around ~40 kHz for the minimum call frequency, ~65 kHz for the maximum call frequency, and ~50 kHz for the frequency of peak energy (Additional file2: Figure S2 for comparison with the echolocating bats included in the study).
Rate shifts in inner ear morphology
An analysis of the rate of morphological change in relative basilar membrane length revealed no change in the rate on the branch leading to all modern bats, as might be expected under the hypothesis of a single gain of echolocation with a subsequent loss in Old World fruit bats. At the same time, most Yangochiroptera branches were shown to have undergone a decrease in rate; this shift was detected with a significant posterior probability (> 0.95) at the basal ancestral Yangochiroptera node (Figure 5A). The only exception to this slowing down was seen in the terminal branches, such as those leading to two Molossus species. A much clearer signature of accelerated basilar membrane evolution was seen in the sub-clade of Rhinolophus philippinensis morphs, albeit with low probability (see Figure 5A). This latter result is of particular interest as these bats are thought to have undergone recent divergence, putatively linked to a harmonic shift in their echolocation call frequency. Our analysis of echolocation calls of the three R. philippinensis morphs suggested that the large morph’s frequency is lower than expected given its forearm size, and the small morph’s is lower than expected given its basilar membrane length (Additional file4: Figure S4). Apart from within the bats, a sub-clade of baleen whales and California sea lion, Zalophus californianus, were also shown to have undergone an increase in the rate of morphological change along their branches.
In contrast to membrane length, rates of change in the number of cochlear turns did not reveal any consistent significant shifts across any bat sub-clade, and instead several increases in rate were detected in scattered clusters of mainly terminal branches (Figure 5B). All Cetartiodactyla branches displayed signatures of accelerated evolution in the number of cochlear turns.
Tests of phylogenetic signal for the three morphological traits showed that the calculated λ value of 0.899 for log basilar membrane length was significantly different from both 0 and 1 (P = 1.87 × 10-28 and 1.30 × 10-24, respectively), as was the λ of 0.963 for log number of turns (P = 6.55 × 10-29 and 3.71 × 10-05, respectively). The λ value of 0.995 for log body mass was significantly different from 0 but not from 1 (P = 4.09 × 10-53 and 0.30, respectively), thus only this trait showed significant phylogenetic signal consistent with Brownian motion (Additional file1: Table S7 for full results).
We undertook three-dimensional reconstructions of the cochleae of 56 bat species, and compared relative basilar membrane length and number of turns among bats and non-echolocating mammals. By relating these structures to echolocation abilities and call parameters, we assessed whether the acquisition of high-frequency hearing and echolocation in bats was associated with morphological change of the inner ear.
Cochlear gross morphology
Across the bat order, cochleae were found to have highly variable morphologies, with evidence of both a phylogenetic signal as well as one relating to echolocation ability. The extent of this morphological variation differed greatly across families, with some (e.g. Rhinolophidae and Pteropodidae) displaying low inter-specific differences, and others (e.g. Hipposideridae and Mormoopidae) showing higher variation. This pattern might in part reflect variation in the echolocation call parameters of our focal species (see Figure 2); for example the hipposiderid Cloeotis percivali has one of the highest known echolocation call frequencies at ~212 kHz, compared to values of 64–157 kHz for other members of this family. The Cloeotis cochlea contains a highly modified basal turn compared to the other species studied, and these modifications are consistent with the tonotopic organisation of the cochlea in which the basal area corresponds to the highest frequencies. This gradation is thought to be at least partly achieved by a decrease in stiffness of the basilar membrane from base to apex. Similar to Cloeotis, Pteronotus parnellii – which along with the Old World horseshoe bats has independently evolved CF echolocation – also shows considerable expansion of the cochlear basal turn (see Figure 2), which probably relates to the well-developed auditory fovea in this taxon (as reviewed by).
Relative basilar membrane length, number of cochlear turns and echolocation
Basilar membrane lengths also varied widely across the bat species studied. All Old World fruit bats were found to have relative basilar membrane lengths similar to those of non-echolocating mammals. Although echolocating bats typically had elongated basilar membranes compared to non-echolocating mammals, we found no consistent significant difference between these groups after accounting for phylogenetic relatedness. Some echolocating bats showed surprisingly short basilar membranes; for example, Macroderma gigas had a similarly proportioned basilar membrane to a marsupial, with values for both species falling below the placental mammal distribution. This is noteworthy given the documented development and elongation of therian inner ear features[5, 7], and suggests that more in-depth sampling is necessary to quantify inter-specific variation accurately within groups. Previous studies have also found that some members of the Megadermatidae have small cochleae for their body mass, and suggest not all echolocating bats have enlarged cochleae. The number of cochlear turns was found to be typically, but not absolutely, higher in the echolocating bats examined (2 - 3.75 turns) compared to the Old World fruit bats (1.75 – 2.25), and published values for other placental mammals (1.5 – 4.25) (Additional file1: Table S2;[49, 50]).
Basilar membrane length displayed an overall negative allometric relationship with body mass, indicating that small mammals have proportionally longer basilar membranes than large mammals. Previously no relationship was found between absolute basilar membrane length and number of turns, however, we found a positive significant correlation after accounting for body mass and phylogeny using Bayesian phylogenetic mixed models. This result suggests that cochlear coiling in order to accommodate the basilar membrane may be more important in small-bodied species. For example, in Old World horseshoe bats, low body mass and long basilar membranes might explain their characteristically high number of cochlear turns.
We found that the number of cochlear turns and basilar membrane length were both correlated with echolocation call parameters after taking phylogeny into account. All measures of echolocation call frequency displayed a significant negative relationship with basilar membrane length, and significant positive association with the number of cochlear turns. Body mass was not found to significantly improve model fit which was somewhat surprising given the previously established relationship between body mass and echolocation call frequency[51, 52]. It is as yet unclear whether the inferred hearing frequency is influenced by the number of turns per se – perhaps due to some mechanical properties – or whether this result is simply an artefact of the relationship between turn number and membrane length. Previously it was shown that cochlear coiling may be an important factor in determining the lower frequency limit of hearing, therefore, potential mechanical implications of cochlear shape on hearing should not be overlooked. In fact, cochlear width, basilar membrane length and number of turns are all likely to be interrelated. The tonotopic organisation of the cochlea means it is not necessarily essential for an animal with good high frequency hearing to have either many turns and/or long basilar membranes; for example, echolocating cetaceans cochleae are not characterised by high numbers of turns but instead have a greatly expanded basal turn. Therefore, given that the optimal cochlear form is likely to vary between different mammalian groups, as well as across species comparative studies, within-group comparisons are also critical for understanding inner ear adaptations.
Of the bats we examined, the horseshoe bats (Rhinolophidae) were found to have the longest relative basilar membrane lengths, which is probably an adaptation to their fine auditory tuning to their CF calls. While call frequencies in horseshoe bats have previously been shown to correlate with cochlear width and body mass e.g.[51, 52], we found that echolocation call frequencies of three size morphs of Rhinolophus philippinensis were lower than expected given their forearm, body mass and basilar membrane length, based on the relationship calculated from other congeners. This supports genetic data suggesting these bats have undergone very recent divergence via call frequency shifts, which could lead to a decoupling between call parameters and morphology also see.
Our study sought to build on previous methods that attempted to relate inner ear measurements with auditory thresholds in vertebrates e.g.[18, 20, 56]. These studies have previously focused on taxa considered to have ‘generalist’ hearing, therefore, it is of particular interest to see how these previously documented relationships may apply to species considered to have ‘specialist’ hearing. However, as discussed above, from correlation studies alone it is impossible to make functional inferences, for example, regarding cochlear coiling. Furthermore although a precedent has been set to correlate low frequency hearing with the product of membrane length and number of turns, the functional rationale behind this has not been explained[18, 20]. It is therefore vital that as three dimensional inner ear datasets continue to be collected, increasingly sophisticated models are developed that will more accurately reflect the functional aspects of vertebrate auditory systems.
Ancestral reconstructions and origins of echolocation
Inferred hearing characteristics based on ancestral reconstructions of cochlear morphology suggested that the ancestor of modern bats was likely capable of perceiving higher frequencies than many other mammals. However, the reconstructed upper hearing limit cannot be taken as proof that this species could echolocate, nor even hear well at very high frequencies; in fact, the audiograms of many non-echolocating mammals, including Old World fruit bats, show hearing limits of >40 KHz, whereas some species of laryngeal echolocating bats (e.g. some members of Molossidae) utilise low frequency echolocation calls, with either the entire FM sweep or just part of the call, within the frequency range audible to humans. Previous morphological examination of the earliest known fossil bat, Onychonycteris finneyi, suggested an inability to echolocate; nevertheless, this taxon was almost certainly not the last common ancestor of the bat ‘crown’ group considered in our study and thus direct comparisons with our results may not be meaningful. Others have attempted to predict the echolocation capabilities of ‘primitive’ fossil bats from their inner ear dimensions[31, 59], however, until the exact placement of the early fossil bat taxa with respect to modern bats is resolved, then their ability to inform our knowledge of the origin of echolocation remains limited.
A persistent problem with attempts at exploring the correlations between bat inner ears and their auditory capabilities is the lack of published audiograms. As echolocation call frequencies are expected to correlate with hearing capabilities (as previously demonstrated) it should be possible to obtain functionally meaningful correlations between inner ear morphology and echolocation parameters. Yet such associations will be less straightforward in bats that use frequency modulated (FM) echolocation, where single calls might cover a range of frequencies and so make parameterisation difficult. Furthermore, the harmonic structure of echolocation calls must be taken into account, for example, the Old World CF bats, utilise calls with the most energy in the second harmonic.
Rate shifts in inner ear morphology
The only significant rate shift (a decrease) in the rate of morphological evolution was found in the relative membrane length of the ancestral Yangochiroptera. Typically, members of the Yangochiroptera were also shown to have longer relative basilar membrane lengths compared to non-echolocating Old World fruit bats. Therefore, this apparent slowing down could suggest that the increase in basilar membrane length was an adaptive trait from an early point in Yangochiroptera evolution. Given the critical role of the basilar membrane in supporting the organ of Corti (which contains the inner and outer hair cells) it might be expected that these structures will show adaptations for processing ultrasonic echoes. Indeed, three genes (Tmc1, Kcnq4 and Pcdh15,) involved in hair cell structure and function show evidence of positive and/or divergent selection acting on the ancestral Yangochiroptera branch[43–45]. The basilar membranes themselves are also known to be long in echolocating members of the Yinpterochiroptera, although this clade did not display evidence of a consistent decrease in the rate of morphological evolution. Instead we found some evidence of an increase in basilar membrane evolution rate that may be related to recent speciation events in horseshoe bats, especially since this group contains among the youngest taxa in our study. Surprisingly, given their highly modified cochleae, no significant rate shifts were found across all branches of the Old World CF bats, which could reflect poor taxonomic sampling, or might mean that neither of the two morphological characters examined in this study fully describes their modified inner ears. For example, in CF bats the presence of auditory foveae, and more generally basilar membrane width, will both be important considerations for modelling bat cochlear mechanics[29, 34, 56].
Neither relative basilar membrane length nor the number of turns showed evidence of a positive shift in morphological change on the ancestral branch of all modern bats, as might be expected if the acquisition of sophisticated laryngeal echolocation had occurred rapidly at this stage of bat evolution. Instead echolocation must have evolved either gradually in this branch without showing a detectable elevated rate of morphological change compared to the background rate, or otherwise it could have evolved later in bat evolution following the divergence of the two main suborders (~64 Ma[39, 62]). The absence of any rate change in the branch of the non-echolocating Old World fruit bats, as might be expected if there had been a loss, would appear to support the latter scenario also see[39, 42–45]. The finding that O. finneyi did not possess an enlarged cochlea whereas other Eocene fossil bats did, provides additional evidence that echolocation was not present in all early bat lineages, and might have arisen over a short evolutionary timeframe.
Of the three morphological characters examined (basilar membrane length, number of turns and body mass) for which phylogenetic signal was estimated, only body mass showed consistent variation across the tree (i.e. following Brownian motion, BM). In comparison, observed species values of both inner ear characters deviated from the pattern expected given the phylogeny (branch lengths and topology). Consequently variation shown across the species is not consistent with BM, indicating that certain types of traditional phylogenetic corrections, such as independent contrasts, may not be suitable using the untransformed tree. Given the sophisticated high-frequency hearing possessed by some laryngeal echolocating bats it is perhaps not surprising that the morphological variation of the two inner ear features studied here were found not to fit the pattern expected under BM. Previous studies that have documented morphological variation in the inner ears of key taxa have focused principally on primates, rodents and cetaceans[22, 49, 63]. In contrast to the results found by this study, it has previously been concluded that primate inner ears have evolved under BM and have been used as informative phylogenetic characters[63, 64], so corroborating assertions that primates have unspecialised ‘generalist’ hearing.
Our study focused on two key parameters of the inner ear, both of which have putatively interacted to play a crucial role in the development of the high-frequency sensitivity that is characteristic of mammalian hearing. While the inner ears of Old World non-echolocating fruit bats did not deviate significantly from other non-echolocating mammals, laryngeal echolocating bats were shown to display highly variable cochleae, and correlations with echolocation call parameters suggest that inner ear morphology is mechanistically linked to call structure in echolocating species. At the same time, patterns in trait variation associated with echolocation were not universal across all echolocating species, and, furthermore, some were not robust after taking phylogeny into account. Ancestral reconstructions suggest that the common ancestor of extant bats had well-developed high-frequency hearing; however, shifts in the rate of morphological evolution suggest that significant changes in inner ear morphology occurred after the two main bat suborders diverged, consistent with multiple origins of echolocation in bats. Finally, our study examined morphological variation of the inner ears from ~5% of extant bat diversity based on recent species estimates. Further fine-scale studies are therefore necessary to fully understand the remarkable morphological diversity of the bat order.
Materials and methods
We studied 56 bat species (n = 68 individuals) from 16 families, with broad taxonomic, geographic, ecological and echolocation call type coverage (Additional file1: Table S1 for species list). Our dataset included three documented size-morphs of Rhinolophus philippinensis, which appear to be incipient species. Specimens were scanned in the frontal plane using the Metris X-Tek HMX ST 225 CT System at the Department of Mineralogy, EMMA Division, NHM, London. Volumes were reconstructed using CT PRO (Metris X-Tek, UK), and following reconstruction volumes were visualized using VG Studio Max 2.0 (Volume Graphics, Heidelberg, Germany). Internal voids of bony labyrinth were digitally dissected to produce digital endocasts, and converted into shells describing the surface geometry with MeshLab v.1.2.2 (MeshLab Visual Computing Lab - ISTI - CNR).
1) Basilar membrane length
Using the ‘single point’ feature in Landmark v3.6, a series of 86 approximately equidistantly placed landmark points were placed along the length of the depression between the scala media and the scala tympani. These points approximated the position of the outer edge of the basilar membrane, beginning at the lowest point of the base of the basilar membrane (where the depression between the two scala is first visible), and ending at the apex of the cochlea (Additional file5: Figure S1A). This number of landmarks should adequately describe the path of the membrane (Additional file5: Figure S1B). The 3D coordinates were exported into Microsoft Excel, where the total Euclidean distance was calculated by summing the distance between each set of consecutive points. Where the distance (x) between points (p1, p2, p3) and (q1, q2, q3) is calculated using the formula:
The measurements collected by this study were combined with those from previous studies (Additional file1: Table S2 for values and sources). In order to compare the relationship between the linear measurements of basilar membrane length and body mass, the cube root of the latter was calculated. Values were then log10 transformed so that the linear regression between the variables could be studied. To test whether data from bats fell into the same distribution as published data from non-bats, we calculated the 95% prediction intervals (PI) of the latter.
2. Number of cochlear spiral turns
The number of cochlear turns was measured in each bat species following West, in which the cochlea was viewed apically and a line drawn from the point of the round window (where the cochlear duct initially begins to curl) to the apex. The number of times the line was crossed by the path of the duct was then recorded. Measurements were taken to the nearest one quarter of a complete turn. We supplemented our measurements with those from literature sources to include additional bat and non-bat species (Additional file1: Table S2).
Constructing the phylogeny and estimating branch lengths
Published cytochrome b sequences, for 131 ingroup species and 3 outgroup species (Additional file1: Table S2) were aligned using ClustalW2 and checked by eye. Branch lengths were estimated using a constrained tree topology based on published phylogenies[69–82]. Cytochrome b sequences were not available for most members of Megadermatidae; therefore we used the phylogeny proposed by Griffiths et al.. We were able to estimate diversification times for the division between Macroderma gigas with Megaderma lyra, and this value and those from Jones et al. were used to date the remaining nodes. Nine fossil calibration points were used (Additional file1: Table S3), each following a normal prior distribution with mean and standard deviation set so that the 5th and 95th quantiles correspond to the published lower and upper suggested node ages respectively. Analyses were run in BEAST v.1.5.4 using an uncorrelated log-normal relaxed molecular clock, a Yule speciation prior and a GTR+I+Γ model, for 10,000,000 generations, with every 1000 parameters logged. Tracer v.1.5 was used to check for appropriate burn-in length and run convergence. The maximum clade credibility tree was produced using TreeAnnotator v.1.5.4, with a sample burn-in of 500 and node heights set to mean-heights.
To test whether observed inter-species variation in cochlear parameters remained after accounting for phylogenetic relatedness, we implemented Bayesian phylogenetic mixed models (BPMMs) in ‘MCMCglmm’ in R v.2.11.1 (http://www.r-project.org) (see Additional file1: Supplementary Methods). Model fit was assessed based the Deviance Information Criterion (DIC) in which ΔDIC values of ≥2 was used to denote significant statistical improvement. For significance tests of fixed effects, we report the PMCMC value, which is twice the posterior probability that a model parameter is greater or less than zero (whichever is lower), as estimated by the Markov chain, and is one possible Bayesian analogue to a two-tailed frequentist p-value.
Inferring ancestral auditory and echolocation capabilities
Ancestral reconstructions of relative basilar membrane length and number of cochlear turns were estimated using the maximum likelihood approach, ‘ace’ function, and phylogenetic independent contrasts, both undertaken in the R package ‘PICANTE’. Low and high frequency auditory thresholds, at 30 and 60 decibels (dB), were extracted from published audiograms for 14 bat and 24 non-bat mammal species (Additional file1: Table S4). The estimated body mass of the ancestral bat, 19.22 g, was taken from Safi et al.. The upper frequency hearing limits of the hypothetical ancestral bat at 30 and 60 dB were estimated using the relationship between log relative membrane length and hearing limits in extant taxa. Lower frequency hearing limits at 30 and 60 dB were estimated using the relationship between the product of log basilar membrane length and the number of cochlear turns with low frequency hearing limits in extant taxa, following published methods[18, 20].
Additionally, to test for an association between bat inner ear morphology and echolocation call frequency, we collected values for three echolocation call parameters: mean frequency at maximum intensity, minimum frequency and maximum frequency. Due to species specific differences in call intensities e.g.[89, 90], echolocation calls frequencies are not recorded at consistent sound levels, and instead by convention the absolute minimum, maximum and peak energy frequencies are typically measured. The audiograms that are available for echolocating bat species suggest that the bats are capable of hearing the majority, if not the entire, range of frequencies covered by the call e.g.. As many values as possible were collected for each echolocation call parameter, and in cases where more than one value was available, we took the mean.
Rates of morphological change
To assess the relative support for the two scenarios proposed to explain the absence of echolocating in Old World fruit bats, we characterised the rate of cochlear evolution (based on coils and length) across the Laurasiatheria. For this we performed a Bayesian analysis of rate shifts in the morphological traits using the package ‘AUTEUR’. Analyses were run twice, with 4,000,000 generations sampled every 4,000 generations. Convergence was assessed using Tracer v1.5. Shift plots were drawn with a burn in of 25%. In this analysis, each branch of the phylogeny is coloured according to the calculated model-averaged rate estimate for that particular branch and the posterior probability of a rate shift occurring at each node is also indicated. For our study we defined significant rate shifts as those with a posterior probability greater than 0.95.
As a further test of character evolution, the phylogenetic signal (Pagel’s λ) of basilar membrane length, body mass and the number of cochlear turns was estimated using the ‘fitContinuous’ argument of the ‘Geiger’ package. The estimated values were then tested to see if they were significantly different to either 0 or 1, where a λ estimate of 1 signifies an exact fit between the phylogeny and a given characters under Brownian motion, and a λ of 0 signifies no phylogenetic signal and thus all species approximate independent points.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its additional files.
Barton RA, Purvis A, Harvey PH: Evolutionary radiations of visual and olfactory brain systems in primates, bats and insectivores. Phil Trans R Soc B. 1995, 348: 381-392. 10.1098/rstb.1995.0076.
Nevo E, In Subterranean Rodents: Mosaic evolution of subterranean mammals: Tinkering, regression, progression, and global convergence. 2007, Springer Berlin Heidelberg: News from Underground. Edited by Begall S, Burda H, Schleich CE, 375-388.
Hendrickson A, Djajadi HR, Nakamura L, Possin DE, Sajuthi D: Nocturnal tarsier retina has both short and long/medium-wavelength cones in an unusual topography. J Comp Neurol. 2000, 424: 719-730.
Heffner HE, Heffner RS: High-frequency hearing. Handbook of the senses: Audition. Edited by: Dallos P, Oertel D, Hoy R. 2008, New York: Elsevier, 55-60.
Manley GA: Cochlear mechanisms from a phylogenetic viewpoint. Proc Natl Acad Sci USA. 2000, 97: 11736-11743. 10.1073/pnas.97.22.11736.
Cantos R, Cole LK, Acampora D, Simeone A, Wu DK: Patterning of the mammalian cochlea. Proc Natl Acad Sci USA. 2000, 97: 11707-11713. 10.1073/pnas.97.22.11707.
Luo ZX, Ruf I, Schultz JA, Martin T: Fossil evidence on evolution of inner ear cochlea in Jurassic mammals. Proc R Soc B. 2011, 278: 28-34. 10.1098/rspb.2010.1148.
Meng J, Fox RC: Osseous inner ear structures and hearing in early marsupials and placentals. Zool J Linn Soc. 1995, 115: 47-71.
Ruggero MA, Temchin AN: The roles of the external, middle, and inner ears in determining the bandwidth of hearing. Proc Natl Acad Sci USA. 2002, 99: 13206-13210. 10.1073/pnas.202492699.
Manley GA: An evolutionary perspective on middle ears. Hear Res. 2010, 263: 3-8. 10.1016/j.heares.2009.09.004.
Zheng J, Shen WX, He DZZ, Kevin BL, Madison LD, Dallos P: Prestin is the motor protein of cochlear outer hair cells. Nature. 2000, 405: 149-155. 10.1038/35012009.
Franchini LF, Elgoyhen AB: Adaptive evolution in mammalian proteins involved in cochlear outer hair cell electromotility. Mol Phylogenet Evol. 2006, 41: 622-635. 10.1016/j.ympev.2006.05.042.
Manley GA: Evolutionary paths to mammalian cochleae. JARO J Assoc Res Otolaryngol. 2012, 13: 733-743. 10.1007/s10162-012-0349-9.
Robles L, Ruggero MA: Mechanics of the mammalian cochlea. Physiol Rev. 2001, 81: 1305-1352.
Amma LL, Goodyear R, Faris JS, Jones I, Ng L, Richardson G, Forrest D: An emilin family extracellular matrix protein identified in the cochlear basilar membrane. Mol Cell Neurosci. 2003, 23: 460-472. 10.1016/S1044-7431(03)00075-7.
Vater M, Kossl M: Comparative aspects of cochlear functional organization in mammals. Hear Res. 2011, 273: 89-99. 10.1016/j.heares.2010.05.018.
Kirk EC, Gosselin-Ildari AD: Cochlear labyrinth volume and hearing abilities in primates. Anat Rec. 2009, 292: 765-776. 10.1002/ar.20907.
West CD: The relationship of the spiral turns of the cochlea and the length of the basilar membrane to the range of audible frequencies in ground dwelling mammals. J Acoust Soc Am. 1985, 77: 1091-1101. 10.1121/1.392227.
Rosowski JJ, Graybeal A: What did Morganucodon hear?. Zool J Linn Soc. 1991, 101: 131-168. 10.1111/j.1096-3642.1991.tb00890.x.
Manoussaki D, Chadwick RS, Ketten DR, Arruda J, Dimitriadis EK, O'Malley JT: The influence of cochlear shape on low-frequency hearing. Proc Natl Acad Sci USA. 2008, 105: 6162-6166. 10.1073/pnas.0710037105.
Heffner R, Heffner H: Hearing in the elephant (Elephas maximus). Science. 1980, 208: 518-520. 10.1126/science.7367876.
Webster DB, Webster M: Morphological adaptations of the ear in the rodent family Heteromyidae. Am Zool. 1980, 20: 247-254.
Nummela S: Scaling of the mammalian middle-ear. Hear Res. 1995, 85: 18-30. 10.1016/0378-5955(95)00030-8.
Ramsier MA, Cunningham AJ, Moritz GL, Finneran JJ, Williams CV, Ong PS, Gursky-Doyen SL, Dominy NJ: Primate communication in the pure ultrasound. Biol Lett. 2012, 8: 508-511. 10.1098/rsbl.2011.1149.
Pye JD: Echolocation and Countermeasures. Bioacoustics: A Comparative Approach. Edited by: Lewis B. 1983, New York: Academic Press, 407-429.
Jones G, Holderied MW: Bat echolocation calls: adaptation and convergent evolution. Proc R Soc B. 2007, 274: 905-912. 10.1098/rspb.2006.0200.
Schmieder DA, Kingston T, Hashim R, Siemers BM: Breaking the trade-off: rainforest bats maximize bandwidth and repetition rate of echolocation calls as they approach prey. Biol Lett. 2010, 6: 604-609. 10.1098/rsbl.2010.0114.
Schnitzler HU, Kalko EKV: Echolocation by insect-eating bats. Bioscience. 2001, 51: 557-569. 10.1641/0006-3568(2001)051[0557:EBIEB]2.0.CO;2.
Pye A: The structure of the cochlea in chiroptera I. Microchiroptera: Emballonuroidea and Rhinolophoidea. J Morphol. 1966, 118: 495-510. 10.1002/jmor.1051180404.
Pye A: The Megachiroptera and Vespertilionoidea of the Microchiroptera. J Morphol. 1966, 119: 101-119. 10.1002/jmor.1051190202.
Habersetzer J, Storch G: Cochlea size in extant Chiroptera and Middle Eocene microchiropterans from Messel. Naturwissenschaften. 1992, 79: 462-466. 10.1007/BF01139198.
Kossl M, Mayer F, Frank G, Faulstich M, Russell IJ: Evolutionary adaptations of cochlear function in Jamaican mormoopid bats. J Comp Physiol Sens Neural Behav Physiol. 1999, 185: 217-228. 10.1007/s003590050381.
Kossl M, Foeller E, Faulstich M: Otoacoustic emissions and cochlear mechanisms in echolocating bats. Echolocation in Bats and Dolphins. Edited by: Thomas JA, Moss CF, Vater M. 2004, Chicago and London: The University of Chicago Press, 104-108.
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. 10.1007/BF00617731.
Russell IJ, Drexl M, Foeller E, Vater M, Kossl M: The development of a single frequency place in the mammalian cochlea: The cochlear resonance in the mustached bat Pteronotus parnellii. J Neurosci. 2003, 23: 10971-10981.
Kucuk B, Abe K: Microstructures of the osseous spiral laminae in the bat cochlea - a scanning electron-microscopic study. Arch Histol Cytol. 1992, 55: 315-319. 10.1679/aohc.55.315.
Teeling EC, Madsen O, Van den Bussche RA, de Jong WW, Stanhope MJ, Springer MS: Microbat paraphyly and the convergent evolution of a key innovation in Old World rhinolophoid microbats. Proc Natl Acad Sci USA. 2002, 99: 1431-1436. 10.1073/pnas.022477199.
Springer MS, Teeling EC, Madsen O, Stanhope MJ, de Jong WW: Integrated fossil and molecular data reconstruct bat echolocation. Proc Natl Acad Sci USA. 2001, 98: 6241-6246. 10.1073/pnas.111551998.
Eick GN, Jacobs DS, Matthee CA: A nuclear DNA phylogenetic perspective on the evolution of echolocation and historical biogeography of extant bats (Chiroptera). Mol Biol Evol. 2005, 22: 1869-1886. 10.1093/molbev/msi180.
Simmons NB, Seymour KL, Habersetzer J, Gunnell GF: Primitive Early Eocene bat from Wyoming and the evolution of flight and echolocation. Nature. 2008, 451: 818-821. 10.1038/nature06549.
Habersetzer J, Richter G, Storch G: Paleoecology of early middle Eocene bats from Messel, FRG. aspects of flight, feeding and echolocation. Hist Biol. 1994, 8: 235-260.
Li G, Wang JH, Rossiter SJ, Jones G, Cotton JA, Zhang SY: The hearing gene Prestin reunites echolocating bats. Proc Natl Acad Sci USA. 2008, 105: 13959-13964. 10.1073/pnas.0802097105.
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. 2012, 108: 480-489. 10.1038/hdy.2011.119.
Liu Y, Han N, Franchini LF, Xu H, Pisciottano F, Elgoyhen AB, Rajan KE, Zhang S: The voltage-gated potassium channel subfamily KQT Member 4 (KCNQ4) displays parallel evolution in echolocating bats. Mol Biol Evol. 2012, 29: 1441-1450. 10.1093/molbev/msr310.
Shen YY, Liang L, Li GS, Murphy RW, Zhang YP: Parallel evolution of auditory genes for echolocation in bats and toothed whales. PLoS Genet. 2012, 8: e1002788-10.1371/journal.pgen.1002788.
Fenton MB, Bell GP: Recognition of species of insectivorous bats by their echolocation calls. J Mammal. 1981, 62: 233-243. 10.2307/1380701.
Kossl M, Foeller E, Drexl M, Vater M, Mora E, Coro F, Russell IJ: Postnatal development of cochlear function in the mustached bat Pteronotus parnellii. J Neurophysiol. 2003, 90: 2261-2273. 10.1152/jn.00100.2003.
Neuweiler G: Evolutionary aspects of bat echolocation. J Comp Physiol Sens Neural Behav Physiol. 2003, 189: 245-256.
Ketten DR: Cetacean Ears. Hearing by Whales and Dolphins. Edited by: Au WWL, Popper AN, Fay RR. 2000, New York: Springer, 43-108. [Fay RR, popper AN (Series Editor): Springer Handbook of Auditory Research]
Begall S, Burda H: Acoustic communication and burrow acoustics are reflected in the ear morphology of the coruro (Spalacopus cyanus, Octodontidae), a social fossorial rodent. J Morphol. 2006, 267: 382-390. 10.1002/jmor.10411.
Jones G: Scaling of echolocation call parameters in bats. J Exp Biol. 1999, 202: 3359-3367.
Francis CM, Habersetzer J: Interspecific and intraspecific variation in echolocation call frequency and morphology of horseshoe bats, Rhinolophus and Hipposideros. Bat Biology and Conservation. Edited by: Kunz TH, Racey P. 1998, Washington and London: Smithsonian Institution Press
Ketten DR: Structure and function in whale ears. Bioacoustics. 1997, 8: 103-135. 10.1080/09524622.1997.9753356.
Kingston T, Rossiter SJ: Harmonic-hopping in Wallacea's bats. Nature. 2004, 429: 654-657. 10.1038/nature02487.
Odendaal LJ, Jacobs DS: Morphological correlates of echolocation frequency in the endemic Cape horseshoe bat, Rhinolophus capensis (Chiroptera: Rhinolophidae). J Comp Physiol A Neuroethol Sens Neural Behav Physiol. 2011, 197: 435-446. 10.1007/s00359-010-0601-x.
Manley GA: Some aspects of the evolution of hearing in vertebrates. Nature. 1971, 230: 506-509. 10.1038/230506a0.
Heffner RS, Koay G, Heffner HE: Hearing in large (Eidolon helvum) and small (Cynopterus brachyotis) non-echolocating fruit bats. Hear Res. 2006, 221: 17-25. 10.1016/j.heares.2006.06.008.
Fenton MB, Jacobs DS, Richardson EJ, Taylor PJ, White W: Individual signatures in the frequency-modulated sweep calls of African large-eared, free-tailed bats Otomops martiensseni (Chiroptera: Molossidae). J Zool. 2004, 262: 11-19. 10.1017/S095283690300431X.
Habersetzer J, Schlosser-Sturm E, Storch G, Sige B: Shoulder joint and inner ear of Tachypteron franzeni, an emballonurid bat from the middle Eocene of Messel. Evolutionary history of bats fossils, molecules and morphology. Edited by: Gunnell G, Simmons N. 2012, Cambridge: Cambridge Studies in Morphology and Molecules: New Paradigms in Evolutionary Bio]
Bohn KM, Moss CF, Wilkinson GS: Correlated evolution between hearing sensitivity and social calls in bats. Biol Lett. 2006, 2: 561-564. 10.1098/rsbl.2006.0501.
Jones G, Teeling EC: The evolution of echolocation in bats. Trends Ecol Evol. 2006, 21: 149-156. 10.1016/j.tree.2006.01.001.
Teeling EC, Springer MS, Madsen O, Bates P, O'Brien SJ, Murphy WJ: A molecular phylogeny for bats illuminates biogeography and the fossil record. Science. 2005, 307: 580-584. 10.1126/science.1105113.
Lebrun R, de Leon MP, Tafforeau P, Zollikofer C: Deep evolutionary roots of strepsirrhine primate labyrinthine morphology. J Anat. 2010, 216: 368-380. 10.1111/j.1469-7580.2009.01177.x.
Gunz P, Ramsier MA, Kuhrig M, Hublin JJ, Spoor F: The mammalian bony labyrinth reconsidered, introducing a comprehensive geometric morphometric approach. J Anat. 2012, 220: 529-543. 10.1111/j.1469-7580.2012.01493.x.
Heffner RS: Primate hearing from a mammalian perspective. The Anatomical Record Part A. 2004, 281A: 1111-1122. 10.1002/ar.a.20117.
Reeder DM, Helgen KM, Wilson DE: Global trends and biases in new mammal species discoveries. Occasional Papers, Museum of Texas Tech University. 2007, 269: 1-36.
Wiley DF: Landmark. 2007, Institute for Data Analysis and Visualization, University of California, Davishttp://www.idav.ucdavis.edu/research/projects/EvoMorph,
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.
Csorba G, Ujhelyi P, Thomas N: Horseshoe bats of the world (Chiroptera: Rhinolophidae). 2003, Shropshire, United Kingdom: Alana Books
Guillén-Servent A, Francis CM: A new species of bat of the Hipposideros bicolor group (Chiroptera: Hipposideridae) from Central Laos, with evidence of convergent evolution with Sundaic taxa. Acta Chiropt. 2006, 8: 39-61. 10.3161/1733-5329(2006)8[39:ANSOBO]2.0.CO;2.
Davalos LM: The geography of diversification in the mormoopids (Chiroptera: Mormoopidae). Biol J Linn Soc. 2006, 88: 101-118. 10.1111/j.1095-8312.2006.00605.x.
Lewis-Oritt N, Porter CA, Baker RJ: Molecular systematics of the family Mormoopidae (Chiroptera) based on cytochrome b and recombination activating gene 2 sequences. Mol Phylogenet Evol. 2001, 20: 426-436. 10.1006/mpev.2001.0978.
Lim BK: Divergence times and origin of neotropical sheath-tailed bats (Tribe Diclidurini) in South America. Mol Phylogenet Evol. 2007, 45: 777-791. 10.1016/j.ympev.2007.09.003.
Ammerman LK, Lee DN, Tipps TM: First molecular phylogenetic insights into the evolution of free-tailed bats in the subfamily Molossinae (Molossidae, Chiroptera). J Mammal. 2012, 93: 12-28. 10.1644/11-MAMM-A-103.1.
Rojas D, Vale A, Ferrero V, Navarro L: When did plants become important to leafnosed bats? Diversification of feeding habits in the family Phyllostomidae. Mol Ecol. 2011, 20: 2217-2228. 10.1111/j.1365-294X.2011.05082.x.
Thabah A, Li G, Wang YN, Liang B, Hu KL, Zhang SY, Jones G: Diet, echolocation calls, and phylogenetic affinities of the great evening bat (Ia io; Vespertilionidae): Another carnivorous bat. J Mammal. 2007, 88: 728-735. 10.1644/06-MAMM-A-167R1.1.
McGowen MR, Spaulding M, Gatesy J: Divergence date estimation and a comprehensive molecular tree of extant cetaceans. Mol Phylogenet Evol. 2009, 53: 891-906. 10.1016/j.ympev.2009.08.018.
Roehrs ZP, Lack JB, Van Den Bussche RA: Tribal phylogenetic relationships within Vespertilioninae (Chiroptera: Vespertilionidae) based on mitochondrial and nuclear sequence data. J Mammal. 2010, 91: 1073-1092. 10.1644/09-MAMM-A-325.1.
Kawai K, Nikaido M, Harada M, Matsumura S, Lin L-K, Wu Y, Hasegawa M, Okada N: Intra- and interfamily relationships of Vespertilionidae inferred by various molecular markers including SINE insertion data. J Mol Evol. 2002, 55: 284-301. 10.1007/s00239-002-2326-0.
Blanga-Kanfi S, Miranda H, Penn O, Pupko T, DeBry RW, Huchon D: Rodent phylogeny revised: analysis of six nuclear genes from all major rodent clades. BMC Evol Biol. 2009, 9: 71-10.1186/1471-2148-9-71.
Zhou X, Xu S, Xu J, Chen B, Zhou K, Yang G: Phylogenomic analysis resolves the interordinal relationships and rapid diversification of the Laurasiatherian mammals. Syst Biol. 2012, 61: 150-164. 10.1093/sysbio/syr089.
Meredith RW, Janecka JE, Gatesy J, Ryder OA, Fisher CA, Teeling EC, Goodbla A, Eizirik E, Simao TLL, Stadler T: Impacts of the Cretaceous Terrestrial Revolution and KPg extinction on mammal diversification. Science. 2011, 334: 521-524. 10.1126/science.1211028.
Griffiths TA, Truckenbrod A, Sponholtz PJ: Systematics of megadermatid bats (Chiroptera, Megadermatidae): based on hyoid morphology. Am Mus Novit. 1992, 3041: 1-21.
Jones KE, Bininda-Emonds ORP, Gittleman JL: Bats, clocks, and rocks: Diversification patterns in chiroptera. Evolution. 2005, 59: 2243-2255.
Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: 699-710.
Hadfield JD: MCMC methods for multi-response generalized linear mixed models: The MCMCglmm R Package. J Stat Softw. 2010, 33: 1-22.
Kembel SW, Cowan PD, Helmus MR, Cornwell WK, Morlon H, Ackerly DD, Blomberg SP, Webb CO: Picante: R tools for integrating phylogenies and ecology. Bioinformatics. 2010, 26: 1463-1464. 10.1093/bioinformatics/btq166.
Safi K, Seid MA, Dechmann DKN: Bigger is not always better - when brains get smaller. Biol Lett. 2005, 1: 283-286. 10.1098/rsbl.2005.0333.
Waters DA, Jones G: Echolocation call structure and intensity in 5 species of insectivorous bats. J Exp Biol. 1995, 198: 475-489.
Surlykke A, Kalko EKV: Echolocating bats cry out loud to detect their prey. PLoS One. 2008, 3: e2036-2010.1371/journal.pone.0002036.
Neuweiler G, Singh S, Sripathi K: Audiograms of a South Indian bat community. J Comp Physiol Sens Neural Behav Physiol. 1984, 154: 133-142. 10.1007/BF00605398.
Eastman JM, Alfaro ME, Joyce P, Hipp AL, Harmon LJ: A novel comparative method for identifying shifts in the rate of character evolution on trees. Evol Dev. 2011, 65: 3578-3589. 10.1111/j.1558-5646.2011.01401.x.
Rambaut A, Drummond AJ: Tracer v1.5. 2009, Available fromhttp://beast.bio.ed.ac.uk/Tracer,
Harmon LJ, Weir JT, Brock CD, Glor RE, Challenger W: GEIGER: investigating evolutionary radiations. Bioinformatics. 2008, 24: 129-131. 10.1093/bioinformatics/btm538.
From the Natural History Museum (London) we are grateful to Norman MacLeod for facilitating data collection, Richard Abel and Stig Walsh for training in CT techniques, and to curators for help with accessing specimens. We also thank Steven Le Comber and Matthew Struebig (QMUL), Paul Bates from the Harrison Zoological Museum (HZM), and the Bogor Zoological Museum (BZM) for access to and loaning specimens. This work was funded by a NERC studentship awarded to K.T.J.D. with additional support from the Centre for Ecology and Evolution and the University of London Central Research Fund.
The authors declare that they have no competing interests.
KTJD and SJR conceived and designed the study. KTJD undertook data collection and analysed the data with input from SJR. IM provided material. KTJD and SJR wrote the paper with input from IM. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Figure S2: Multiple regression plots of echolocation call parameters, basilar membrane length, number of cochlear turns and body mass. Stepwise multiple regressions suggest that equations with only inner ear parameters were the best fitting models: log maximum frequency = -0.96 log basilar membrane + 1.58 log turns + 2.22, multiple R2= 0.26, F = 9.01 (2, 51 d.f.), P = 4 x 10-4; log peak energy frequency = -1.01 log basilar membrane + 2.26 log turns + 1.92, multiple R2= 0.36, F = 14.05 (2, 51 d.f.), P = 1.38 x 10-5; log minimum frequency = -1.02 log basilar membrane + 2.76 log turns + 1.60, multiple R2= 0.40, F = 17.09 (2, 51 d.f.), P = 2.09 x 10-6. (PDF 168 KB)
Additional file 3: Figure S3: Maximum likelihood ancestral reconstructions of bat inner ears - (A) relative basilar membrane length and (B) number of cochlear turns. Phylogenies and character values are depicted as ‘Traitgrams’, whereby the position along the y-axis corresponds to node age in millions of years and position along the x-axis corresponds to the reconstructed character value. Coloured bars indicate key subdivisions within bats: Old World fruit bats (orange); echolocating Yinpterochiroptera (blue); Yangochiroptera (green). Keys nodes: bat common ancestor (a); Yangochiroptera common ancestor (b); Yinpterochiroptera common ancestor (c); Old World fruit bat common ancestor (d); echolocating Yinpterochiroptera common ancestor (e). (PDF 280 KB)
Additional file 4: Figure S4: Morphological parameters versus echolocation call frequency in Rhinolophus species. Published values for species taken from literature (black points), the three Rhinolophus philippinensis size morphs measured by this study: small (red), medium (orange), large (yellow) and values published for one R. philippinensis values taken from  (blue). (A) Average forearm length, body mass and echolocation call frequency for Rhinolophus spp. from values obtained from literature sources. A significant negative relationship was found (log CF = -1.54 log forearm + 4.44; R2 = 0.42, F = 36.57, P = 1.98 x 10-7 and log CF = -0.38 log body mass + 2.24; R2 = 0.25, F = 12.34, P = 0.001). (B) The relationship between basilar membrane length and constant frequency echolocation call (log CF = -0.98 log basilar membrane + 3.09; R2 = 0.63, F = 21.74, P < 0.001). (PDF 310 KB)
Additional file 5: Figure S1: Measuring basilar membrane length from reconstructed inner ear endocasts. (A) Left: Medial view of the right cochlear endocast of Craseonycteris thonglongyai (specimen number HZM.1.34982, ref. Table S1). A representation of the path of the basilar membrane measured by this study is shown by the dotted line. Right upper: Apical view of the right cochlear endocast of C. thonglongyai (HZM.1.34982). Black arrows correspond to the end point of the representation of the path of the basilar membrane (dotted line) measured by this study. Right lower: Medial view of the right cochlear endocast of C. thonglongyai (HZM.1.34982). Black arrows correspond to the start point of the representation of the path of the basilar membrane (dotted line) measured by this study. (B) Two-dimensional plots showing the representative paths of the basilar membrane for the right cochlea of Pipistrellus pipistrellus, using either 86 or 43 landmark points, connected with straight connecting lines. The basilar membrane path as depicted by a smoothed curvilinear path is also superimposed over these points. The estimated length calculated from the subset of 43 coordinates was only 8.304 mm, compared to 8.473 mm from 86 coordinates. This corresponds to a negative difference of 0.170 mm or a 2% underestimate of membrane length. Furthermore, the path traced by the straight lines connecting the 86 points much more faithfully follows that of the curved path. Therefore, 86 landmark points were deemed to be a suitable compromise between efficiency and accuracy and was used to collect all basilar membrane estimates. The 86 landmarks used in this study (circles); curved path between points (black line); straight lines between points used to estimate basilar membrane length (dotted line); 43 points (white circles), and the dashed line the straight line distance between white circles (dashed line). (PDF 427 KB)
About this article
Cite this article
Davies, K.T., Maryanto, I. & Rossiter, S.J. Evolutionary origins of ultrasonic hearing and laryngeal echolocation in bats inferred from morphological analyses of the inner ear. Front Zool 10, 2 (2013) doi:10.1186/1742-9994-10-2
- Ancestral reconstruction