Abstract The ability of bats and toothed whales to echolocate is a remarkable case of convergent evolution. Previous genetic studies have documented parallel evolution of nucleotide sequences in Prestin and KCNQ4, both of which are associated with voltage motility during the cochlear amplification of signals. Echolocation involves complex mechanisms. The most important factors include cochlear amplification, nerve transmission, and signal re-coding. Herein, we screen three genes that play different roles in this auditory system. Cadherin 23 (Cdh23) and its ligand, protocadherin 15 (Pcdh15), are essential for bundling motility in the sensory hair. Otoferlin (Otof) responds to nerve signal transmission in the auditory inner hair cell. Signals of parallel evolution occur in all three genes in the three groups of echolocators—two groups of bats (Yangochiroptera and Rhinolophoidea) plus the dolphin. Significant signals of positive selection also occur in Cdh23 in the Rhinolophoidea and dolphin, and Pcdh15 in Yangochiroptera. In addition, adult echolocating bats have higher levels of Otof expression in the auditory cortex than do their embryos and non-echolocation bats. Cdh23 and Pcdh15 encode the upper and lower parts of tip-links, and both genes show signals of convergent evolution and positive selection in echolocators, implying that they may co-evolve to optimize cochlear amplification. Convergent evolution and expression patterns of Otof suggest the potential role of nerve and brain in echolocation. Our synthesis of gene sequence and gene expression analyses reveals that positive selection, parallel evolution, and perhaps co-evolution and gene expression affect multiple hearing genes that play different roles in audition, including voltage and bundle motility in cochlear amplification, nerve transmission, and brain function.

Author Summary The convergent origins of laryngeal echolocation in two groups of bats (Yangochiroptera and Rhinolophoidea) and toothed whales have long been a focus of interest for biologists. We screened three candidate genes—Cdh23, Pcdh15, and Otof—involved in different steps in the echolocation system. Signals of parallel evolution occurred in all three genes in the three groups of echolocators. Cdh23 and Pcdh15 constitute part of the mechanical link within the hair bundle of the ear. Both genes showed signals of both convergent evolution and positive selection, which implied they may have co-evolved to optimize cochlear amplification. Further, three lines of evidence suggest that Otof plays an important role in the transmission of signals in the brain during echolocation. First, the gene is more highly expressed in the auditory cortex of the brain in echolocating adult bats than in other cerebral cortexes. Second, this expression is higher in adult echolocating female bats than in their embryos, which do not use echolocation. Third, echolocators also have a higher level of expression in their cerebral cortexes than do non-echolocating bats. Taken with other evidence, the independent origins of echolocation involve the same genes that have evolved in precisely identical ways.

Citation: Shen Y-Y, Liang L, Li G-S, Murphy RW, Zhang Y-P (2012) Parallel Evolution of Auditory Genes for Echolocation in Bats and Toothed Whales. PLoS Genet 8(6): e1002788. https://doi.org/10.1371/journal.pgen.1002788 Editor: Jianzhi Zhang, University of Michigan, United States of America Received: February 16, 2012; Accepted: May 12, 2012; Published: June 28, 2012 Copyright: © 2012 Shen et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: This work was supported by the National Natural Science Foundation of China (30621092 and 31172080) and the Bureau of Science and Technology of Yunnan Province. RWM was supported by a Visiting Professorship for Senior International Scientists from the Chinese Academy of Sciences, and manuscript preparation was supported by a Natural Sciences and Engineering Research Council of Canada Discovery Grant A3148. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist.

Introduction The ability of echolocation using ultrahigh frequency sounds occurs in two groups of bats (Yangochiroptera and Rhinolophoidea) and in toothed whales including dolphins [1]–[3]. These mammals use this complex bio-sonar system to assist with orientation and feeding [4], [5]. Echolocation by bats and dolphins provides an iconic example of either parallel or convergent evolution via natural selection. Previous molecular studies on echolocation have mainly focused on the Organ of Corti. In this organ, the motor protein prestin plays a key role in voltage motility [6]–[8]. It appears to have undergone sequence convergence between bats and dolphins [1], [2], as well as within laryngeal echolocating bats [9]. Further, the voltage-gated potassium channel gene KCNQ4 underwent parallel evolution in echolocating bats [3], [10]. Mammalian audition requires not only voltage motility, but also hair bundle motility, which is executed by outer hair cells in the cochlea [11]. Proteins encoded by the genes Cdh23 and Pcdh15 are essential to hair bundle motility [12]–[14], and their malfunctions in humans cause deafness in newborns and progressive retinitis pigmentosa (Usher syndrome type I) [15]. Homodimers of Cdh23 and Pcdh15 directly link to each other via their amino termini; they form the upper and lower part of tip-links, respectively (Figure S1), which lie between the stereocilia within the hair bundle [14], [16], [17]. The auditory system involves the perception and enhancement of sound signals, as well as transformation of the mechanical signals to ion fluxes in inner hair cells. Management of the electric signals to the brain involves a series of nerve channel openings [18]. Genetic mutations in the gene encoding otoferlin (Otof) cause a clinical, autosomal recessive nonsyndromic form of prelingual and sensorineural deafness [19]–[21]. This protein that may act as the major Ca2+ sensor that triggers membrane fusion at the ribbon synapse of the auditory inner hair cell [22]. Although the above functions are involved in the conversion of sound signals into electrical impulses in the inner ear, the expression of Otof also occurs in neurons and nerve fibers in the brain [23]. The molecular mechanism of voltage motility in echolocation is widely studied. Echolocation is a complex system that includes signal reception by hair cells in the Organ of Corti, nerve transmission, and signal processing in the brain [24]. Therefore, herein we investigate the gene sequence evolution of Cdh23, Pcdh15, and Otof. These proteins function in different steps during echolocation. Because the brain modulates sensory information from peripheral sensory organs [25], and because Otof is involved in transferring sound signal by electrical impulses, we also examine the expression patterns of Otof in the cerebral cortexes of different species. We synthesize evidence from sequences and expressions to study the convergent evolution of echolocation in bats and dolphins.

Discussion Morphological development is very complex and involves a suite of genes. Some mammals have independently developed similar features, such as the ability to echolocate objects. Echolocation by bats and dolphins provides an extreme example of parallel or convergent evolution. Although the morphologies of sending and receiving sonic signals differ greatly, the hearing of ultrasonic sounds and the mechanisms of decoding signals are shared [29], [30]. Thus, the genes coding for the auditory system become ideal candidates for adaptive evolution at the molecular level during echolocation. The genes Cdh23 and Pcdh15 are essential to hair bundle motility [12]–[14], [31], [32]. Otof encodes a protein that may act as the major Ca2+ sensor to trigger membrane fusion at the auditory inner hair cell ribbon synapse [22]. The three genes are involved in different steps of the hearing system. Their nucleotide gene trees are largely congruent with the species tree and the nodes enjoy high support. In contrast, the amino acid trees conflict with the species tree; all unite the echolocators. Further, trees based on nonsynonymous sites have topologies similar to those based on corresponding amino acids, and yet the trees based on synonymous mutations do not show this pattern. Clearly, the difference in branching order is the result of amino acid changes (nonsynonymous mutations). Three independent lineages of mammals can echolocate and it is important to identify whether this involves convergent or parallel evolution. If echolocation involves convergent evolution, then similar traits or functions independently emerge in two or more lineages from different ancestral states. Parallel evolution differs in that similar ancestral traits descend into similar extant states in different lineages [33]. Reconstructions of the ancestral sequences of the internal nodes can detect which amino acid changes cause the incongruence between the amino acid and species trees. A suite of sites show parallel changes in these echolocators, and these changes have a statistically significant signal of parallel evolution. Twenty-one amino acid sites converge on the same residue in branch g and branch b and with a highly significant probability (P<0.001). One convergent site appears in branches g and d (P<0.001), and four parallel sites occur in branches b and d (P<0.001) in Cdh23. In Pcdh15, 22 and three convergent sites occur between branch g and branches b and d, respectively (P<0.001 in both cases) and eight parallel sites have been identified between branches b and d (P<0.001). Four parallel sites in Otof occur on branches g and b (P<0.001). Convergent evolution is not accommodated by current phylogenetic methods and it can strongly mislead phylogenetic inference [34]. Upon excluding the parallel-evolved amino acid sites, the reconstructed amino acid trees did not incorrectly unite all echolocators (Figures S3, S5, S7). Thus, the discovery of multiple parallel-evolving amino acid sites explains the unnatural uniting of echolocators in the amino acid trees. In addition to the parallel evolution of Prestin [1], [2], [9] and KCNQ4 [3], our analyses document a high level of complexity in a large-scale, multigene adaptation. Functional assays have proven that parallel and convergent amino acid changes are responsible for parallel and convergent functional changes [35], [36]. Thus, parallel evolution of gene sequences may have driven phenotypic and functional convergence in echolocating bats and dolphins. Further functional assays are needed to affirm this association. Our analyses of selection pressure detect signals of positive selection in Cdh23 on the branches leading to the Rhinolophoidea and dolphin. The same occurs in Pcdh15 along the branch leading to Yinpterochiroptera (see Table 1 and Table S1). Positive selection appears to have acted on these genes to fit the requirements for echolocation. Whereas outer hair cells amplify sound by somatic [37], [38] or ciliary [39], [40] mechanisms, inner hair cells are passive detectors of the amplified vibratory signal. The signals are converted into electrical impulses by activating fibers of the cochlear (auditory) nerve [41], which then sends the signals to higher auditory processing centers in the brain. Prestin mediates the voltage somatic motility unique to mammals [42]–[44]. Part of the cochlear amplification, Cdh23 and Pcdh15 participate in hair bundle motility [39], [45], [46]. Otof participates by releasing neurotransmitter to nerves [22]. Convergent evolution and positive selection on these genes reflect the pathway from receipt of signal to signal amplification, and then to neural transduction. All parts of the pathway are involved in the auditory system, and thus may play important roles in echolocation. Homodimers of Cdh23 and Pcdh15 directly link to each other via their amino termini, and they constitute the upper and lower part of tip-links, respectively, that lie between the stereocilia within the hair bundle [14], [16], [17]. Sequence alignments suggest that Cdh23 and Pcdh15 have 27 and 11 extracellular cadherin repeats, respectively [47]. Functional evidence encompassing a classical genetic approach shows that mutations at these two cadherin proteins can interact to cause hearing loss in digenic heterozygotes of both mice and humans [48]. Cdh23 and Pcdh15 appear to have undergone both convergent functional evolution and positive selection in echolocators. This finding suggests a strong interaction between the proteins and co-evolution is likely to have optimized their function in cochlear amplification for the development of echolocation. Dolphins and bats employ different tools to echolocate. Whereas the larynx generates sound in echolocating bats [29], the monkey lip/dorsal bursa complex does the same in dolphins [30]. In bats, the sound involves a constant frequency (CF) and frequency modulation (FM), but dolphins use FM and amplitude modulation (AM) [49]. Convergence occurs via both auditory systems being adapted to receiving and processing ultrahigh frequency sounds. Our study discovers evidence of parallel sequence evolution in three genes involved in hearing and this may indicate a genetic basis for echolocation. Sequence evolution only tells a part of the story. Expression pattern is usually a strong indicator of protein-demand and function [50], [51]. The central auditory system plays the crucial role of receiving impulses from auditory nerves and sending messages back to the cochlea [52]. Because Otof plays a role in neural signal transmission, we have evaluated its patterns of expression in different cerebral cortexes. Otof is most highly expressed in the auditory cortex of echolocating adult female Common Bent-Wing Bats. Adult females, which frequently use ultrasonic sounds to explore their environments, have higher levels of expression than their embryos, which do not use echolocation. Further, a comparison of expression of Otof between the brains of an adult echolocating bat (Common Bent-Wing Bat, Miniopterus schreibersii) and a non-echolocating bat (Old World Fruit Bat, Rousettus leschenaultii) reveals a higher level of expression in the former. Indeed, the gene exhibits great differences in levels of expression in both different cerebral cortexes and in species with or without the ability to echolocate. Echolocation signals begin at the hair cells in the Organ of Corti, continue along the auditory nerve, and terminate in the auditory cortex of the brain [24]. Combined with sequence data and expression data, we conclude that multiple instances of parallel sequence evolution are involved in genes in different parts of auditory system between the three groups of echolocators. This occurs not only in well-studied voltage motility, but also bundle motility, and not only in cochlear amplification, but also in neural transduction. Further, co-evolution optimizes the function of homodimers of Cdh23 and Pcdh15 in cochlear amplification. The expression pattern of Otof in different cerebral cortexes implies that the evolution of gene expression might be required for echolocation. In conclusion, we synthesize gene sequence and gene expression analyses and conclude that positive selection, convergent evolution, and perhaps co-evolution and gene expression evolution play roles in audition (voltage motility and bundle motility in cochlear amplification, nerve transmission, and brain) during the independent origins of echolocation in bats and dolphins.

Materials and Methods Ethics statement All research involving animals used in this study followed the guidelines and bylaws on animal experimentation. Bats were anesthetized using an intraperitoneal injection of sodium pentobarbital (C 11 H 17 N 2 NaO 3 ) at a dosage of 100 mg per kg body weight. Following anesthesia, bats were euthanized and then their brains were sampled. Protocols were approved by the Ethics and Experimental Animal Committee of the Kunming Institute of Zoology, Chinese Academy of Sciences. Source of data and primary treatments Sixteen species of bats were used in our analysis (Table S3). Total RNA was isolated from the brain using a RNAiso Plus Kit (Takara, China), and RT-PCR was performed on 2 ug of RNA using the PrimeScript RT-PCR Kit (Takara, China) to obtain cDNA. Subsequently, genes were amplified from cDNA using gene-specific primers (Table S4). PCR products were purified on a 1% agarose gel and a Watson Gel Purification Kit (Watson BioTechnology, Shanghai), and finally transformed into the pMD18-T vector (Takara, China). Each strand was sequenced in both directions with an ABI 3730 sequencer. RNA samples of brain were not available for the dolphin (Tursiops truncatus), so genomic sequences were amplified from total genomic DNA, which was extracted from muscle tissue using a standard 3-step phenol/chloroform extraction method [53]. Raw nucleotide sequences were edited using Lasergene SeqMan software (DNASTAR Inc., Madison, WI, USA). Newly determined sequences were deposited in GenBank (Accession numbers JF808081–JF808094, JQ284400–JQ284430). The sequences of background species came from the Ensembl database (Release 66) and those of high quality were used (Table S3). All the sequences were aligned using ClustalX 1.81 [54] and then visually checked for accuracy (the aligned sequences are available by request). Phylogenetic and molecular evolutionary analyses The best-fit models were selected by jModeltest v0.1.1 [55], [56] for nucleotide sequences and ProtTest 3.0 beta [57] for amino acid sequences. Maximum likelihood (ML) trees were reconstructed by PAUP* [58] with 1,000 replications, and Bayesian inference (BI) trees were reconstructed by MrBayes 3.1.2 with 1,000,000 replications [59], [60]. Neighbor-joining (NJ) phenograms were based on Kimura 2-parameter corrected distances for nucleotide sequences and uncorrected P-distances for amino acid sequences, each with 1,000 bootstrap replications. We implemented the Li-Wu-Luo method [61] to reconstruct NJ trees based on both synonymous and nonsynonymous sites. The sequences of the internal nodes were reconstructed using distance-based Bayesian methods, which included the branch lengths estimated by the least squares method and the ancestral amino acids inferred by the Bayesian approach. These data were used to obtain an unbiased estimate of the true probability [62]. Convergent and parallel amino acid substitutions along each lineage were detected. The statistical significance of these amino acid changes was tested with the method developed by Zhang and Kumar [63]. The CODEML program in PAML 4 [64] was used to detect selective pressure. The species tree [26]–[28] was used as guide tree for analysis. Four models of evolution were used: one-ratio model, free-ratio model, branch models, and branch-site models. For the last two models, three groups of echolocators were set as the foreground branch to detect whether or not they had undergone positive selection. Expressional data analysis We sampled auditory cortex, visual cortex, motor and sensory cortex, olfactory bulb, and cerebellum from euthanized adult and embryonic Miniopterus schreibersii, and adult Rousettus leschenaultii, which represented the adult echolocating bats, embryonic echolocating bats, and adult non-echolocating bats, respectively. Tissue was selected according the human brain atlas [65]. These tissues were stored in liquid nitrogen. Total RNA was isolated using a RNAiso Plus kit (Takara, China). Next, DNA-free RNA samples were condensed using a RNAqueous-4PCR Kit (Applied Biosystems, US). We synthesized cDNA using a PrimeScript RT-PCR Kit (Takara, China), and then used it in Real-Time PCR. TaqMan Gene Expression Assays were custom designed by Applied Biosystems based on our sequencing data of Otof and Actb in bats. The assay details were listed in Table S5. Sequences of Otof and Actb were amplified and detected using an ABI PRISM 7000 Sequence Detection System with a PCR profile as follows: 50°C for 2 min, 95°C for 10 min, followed by 40 cycles at 95°C for 15 s, and 60°C for 1 min. The products were purified and sequenced in both directions with an ABI 3730 sequencer to insure the assays' specificity. Real-Time PCR was performed on 96-well reaction plates in a 20 µl reaction volume containing 100 ng of cDNA per reaction with TaqMan Universal Master Mix II (Applied Biosystems, US). For each group, three individuals were used and each sample was performed three times. Expression data for the target gene Otof were normalized relative to the housekeeping gene Actb. Raw data were obtained and analyzed using the 7000 SDS 1.1 software (Applied Biosystems, US). The comparative C T method (ΔΔCt) was chosen to further calculate the relative expressions between different groups.

Acknowledgments We thank Jing-Kuan Wei from the Laboratory of Sensory Motor Integration in the Kunming Institute of Zoology for his help in sampling different cerebral cortexes from bats. We thank the Ethics and Experimental Animal Committee of the Kunming Institute of Zoology for approving our experiments. We also thank the three anonymous reviewers, Christopher Blair, and Ross MacCulloch for their help in improving the manuscript.

Author Contributions Conceived and designed the experiments: Y-YS Y-PZ. Performed the experiments: LL G-SL. Analyzed the data: Y-YS LL. Contributed reagents/materials/analysis tools: Y-PZ. Wrote the paper: Y-YS LL RWM Y-PZ.