- Research article
- Open Access
Comparative inner ear transcriptome analysis between the Rickett’s big-footed bats (Myotis ricketti) and the greater short-nosed fruit bats (Cynopterus sphinx)
© Dong et al.; licensee BioMed Central Ltd. 2013
- Received: 8 June 2013
- Accepted: 18 December 2013
- Published: 23 December 2013
Bats have aroused great interests of researchers for the sake of their advanced echolocation system. However, this highly specialized trait is not characteristic of Old World fruit bats.
To comprehensively explore the underlying molecular basis between echolocating and non-echolocating bats, we employed a sequence-based approach to compare the inner ear expression difference between the Rickett’s big-footed bat (Myotis ricketti, echolocating bat) and the Greater short-nosed fruit bat (Cynopterus sphinx, non-echolocating bat). De novo sequence assemblies were developed for both species. The results showed that the biological implications of up-regulated genes in M. ricketti were significantly over-represented in biological process categories such as ‘cochlea morphogenesis’, ‘inner ear morphogenesis’ and ‘sensory perception of sound’, which are consistent with the inner ear morphological and physiological differentiation between the two bat species. Moreover, the expression of TMC1 gene confirmed its important function in echolocating bats.
Our work presents the first transcriptome comparison between echolocating and non-echolocating bats, and provides information about the genetic basis of their distinct hearing traits.
Bats belong to the order Chiroptera, which is one of the largest monophyletic clades in mammals. They constitute ~20% of living mammalian species with about 1,200 bat species worldwide [1–3]. Bats have long been regarded as special animals for the sake of being mysterious flyers of the night, and one of the few mammals that can use echolocation to navigate in the dark . Although many mammalian species are capable of detecting ultrasonic sounds (>20 kHz), only toothed whales and laryngeal echolocating bats have the most advanced auditory systems for detecting ultrasonic sound . Laryngeal echolocating bats can emit broad ultrasound frequency sounds and listen to their echoes, allowing them to orient in the dark and exploit food sources of the night sky. However, not all bats have laryngeal echolocation ability. Traditionally, bats can be divided into two suborders, the ‘Megachiroptera’ (megabats) and the ‘Microchiroptera’ (microbats), identified mainly on the basis of their morphology and echolocation characteristics . Megachiroptera, consisting of a single family, is a smaller group of bats that mainly live in the Old World tropical places. They have good dim light vision and do not have laryngeal echolocation ability, whereas Microchiropteran bats can use sophisticated laryngeal echolocation for orientation [6, 7].
With special morphological and physiological adaptations, echolocation allows bats to listen for their high frequency echoes to locate their prey in the dark. The cochleae of these bats are well-developed and show adaptations for detecting their own ultrasonic sounds . The cochlea is a key auditory system of the inner ear and is specialized for the use of high frequency sounds in the echolocating bats. It has been documented that the variation of cochlear size is associated with echolocation strategies [8, 9], and the cochlea diameters of laryngeal echolocating bats are higher than those of many non-echolocating megabats [8, 10, 11]. These phenotypic adaptations are thought to play important roles in the reception and funneling of high frequency echoes. Laryngeal echolocation is a highly technical and physiological adaption, however the molecular basis responsible for this phenotype is poorly characterized for the differences between echolocating and non-echolocating bats.
Divergence of gene expression is an important component of species evolution and essential means to generate biological diversity . Recent next-generation sequencing technologies provided us a large-scale platform to address evolutionary questions involving non-model organisms for which there are still limited genomic resources [13, 14]. In this work, we explored the utility of next-generation sequencing technologies for the comparative inner ear transcriptome analyses between the Rickett’s big-footed bat (Myotis ricketti) and the Greater short-nosed fruit bat (Cynopterus sphinx). The Rickett’s big-footed bats are species of microbats in the family Vespertilionidae. They have advanced laryngeal echolocation ability, which belongs to frequency modulated type; the Greater short-nosed fruit bats are species of megabats and come from the family Pteropodidae without laryngeal echolocation. Based on more recent molecular phylogenetic studies, the divergence time between these two bat species is around 60 million years ago (MYA) [3, 15]. The genome sequences of these two bat species are still lacking, and we employed a de novo assembly approach to gain insight into the genome-wide expression divergence patterns between the Rickett’s big-footed bat and Greater short-nosed fruit bat. The main purpose of this work is to elucidate the nature of transcriptomes of bat species and determine to what extent the two bat species with different hearing traits differ in gene expression. This study provided an initial step to comprehensively understand the inner ear transcriptome involved in the bat echolocation.
De novoassembly and functional annotation
Sequencing of the mRNA in the inner ear of Rickett’s big-footed bat (M. ricketti, echolocating bat) and the Greater short-nosed fruit bat (C. sphinx, non-echolocating bat) based on Illumina Genome Analyzer II generated a total of 6.8 Gbp of sequence from approximately 90 million paired-end 75 bp reads (3.3 GB and 3.5 GB for M. ricketti and C. sphinx, respectively).
Summary of the inner ear transcriptome sequencing assemblies in M. ricketti and C. sphinx
Max contig length
Mean contig length
Number of contigs
Number of contigs (> = 1 k bp)
Transcripts differentially expressed between M. ricketti and C. sphinx
Significantly enriched GO terms of differentially expressed genes between M. ricketti and C. sphinx
No. of genes
Up-regulated genes in the M. ricketti
Muscle filament sliding
Striated muscle contraction
Sensory perception of sound
Inner ear morphogenesis
Structural constituent of muscle
Enhancer sequence-specific DNA binding
Up-regulated genes in the C. sphinx
M phase of mitotic cell cycle
DNA strand elongation involved in DNA replication
Mitotic cell cycle
Zinc ion binding
Based on the alignments of 11,717 orthologous genes, d N /d S values were estimated between two bat species. d N /d S > 1 indicates that strong positive selection has acted to change the sequence, while lower d N /d S value means a more conservative trademark. The average d N /d S rate for all orthologous genes is 0.18 ± 0.013 (mean ± SD). When considering hearing related genes, the average d N /d S rate genes is 0.12 ± 0.041(mean ± SD), suggesting that hearing related genes are under strong functional constraint. We next examined the relationships between expression divergence and sequence evolution, and found the expression divergence was largely independent from the sequence evolution (r = 0.00017, P = 0.9). However, the expression divergences of hearing related genes are negatively correlated with sequence evolution (r = -0.16, P < 1e-23).
Quantitative RT-PCR analysis of TMC1gene in five mammalian species
The significantly up-regulated hearing related genes in M. ricketti
Transcript abundance ( M.ricketti)
Transcript abundance ( C.sphinx)
Transmembrane channel-like 1
Lipoxygenase homology domains 1
Collagen, type XI, alpha 2
Usher syndrome 1C (autosomal recessive, severe)
Carcinoembryonic antigen-related cell adhesion molecule 16
Caspase 3, apoptosis-related cysteine peptidase
Gap junction protein, beta 2, 26 kDa
Gap junction protein, beta 6, 30 kDa
Cholinergic receptor, nicotinic, alpha 9 (neuronal)
POU class 3 homeobox 4
NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 9, 22 kDa
SIX homeobox 1
SRY (sex determining region Y)-box 2
Hexosaminidase A (alpha polypeptide)
We have performed the novel inner ear transcriptome of two bat species (M. ricketti and C. sphinx) using Illumina sequencing technology. As was first proposed for over several decades ago, alterations (or innovations) in gene expressions were regarded as essential means to generate biological diversity [12, 25]. The characterization of transcriptome is essential to distinguish the functional implications of different species and to obtain a better understanding of their biological complexity. So, the analysis of differentially expressed genes can elucidate the molecular mechanisms underlying the morphological diversity and provide a better understanding of the relationship between gene expression patterns and the resultant morphologies. To our knowledge, this work represents the first effort to comparatively analyze the inner ear transcriptome of the bat species resulting in a reference transcriptome of more than 16,000 annotated genes from the M. ricketti and C. sphinx which represent the echolocating bat and non-echolocating bat, respectively.
Our transcriptome work based on next-generation sequencing technology attempted to reveal the underlying molecular mechanisms of the auditory system between the echolocating bat and non-echolocating bat at transcription level. Compared with other next-generation sequencing technology, such as Roche 454 technology, the Illumina solexa platform offers a higher sequencing depth with considerably less cost, which ensures more complete coverage of the transcriptome. In this work, 104,987 and 171,394 contigs were assembled for M. ricketti and C. sphinx, respectively. We evaluated the quality and quantity of short-read assemblies and proved that short-reads transcriptome assemblies are large in quantity and high in quality for further analysis, which can provide considerable utilities for non-model organisms.
Bats are amongst the few mammal species that use sophisticated echolocation. The generated high-frequency calls bounces off surrounding objects and the returning echoes allows them for the detection, localization and classification of these objects with extraordinary acuity. This astonishing sensory ability enables bats to navigate in the dark. Although bats are not the only mammals to have evolved ultrasonic echolocation, they are certainly the most developed species of the laryngeal echolocation. Fossil evidence suggests that bats might echolocate even at an early stage of their evolution . However, none of the Old World fruit bats have laryngeal echolocation ability, and they use their sense of vision to locate food . As we know, the sound waves enter the ear and travel until they reach the middle ear. The auditory ossicles will vibrate as a response to the sound waves, and the stirrup bone transmits vibrations into the inner ear. Then, the pressure waves flatten the hair of the inner ear and made them perceived as a sound. The inner ear consists of a cochlea which is a spiral-shaped cavity that functions as sound reception and processing apparatus for hearing. Cochlea functions as sound reception and processing apparatus in the inner ear, and its size is known to be correlated with echolocation behavior, and the cochlea structure significantly contributed to the diversification of bat species . A recent work reconstructed three-dimensional inner ear volumes of both echolocating and non-echolocating bats, and found that hearing in bats correlated with both measures of cochlea morphology . In this work, we evaluated and compared the expression difference between echolocating and non-echolocating bats, and found that hearing-related genes show significantly differential expression. Organs in the body always exhibit specialized forms that are essential for their functions. Although bat echolocation calls show a great diversity in duration and shape, all echolocating bats have specialized inner ears which allow them to hear sounds in the ultrasonic range. The expressions of inner ear/cochlea morphogenesis genes are responsible for the senses of hearing and balance, which might associate with high-frequency hearing ability. Based on morphological and anatomical analyses, the cochlea is enlarged relative to other skull structures in the echolocating bats, which makes them able to detect, and discriminate between high-frequency calls, whereas non-echolocating bats tend to have smaller cochleae than bats that use laryngeal echolocation . As we know, the generation of animal inner ear requires coordination between morphogenesis and cell fate specification . The inner ear/cochlea morphogenesis genes are highly expressed in echolocating bats, which are consistent with these findings that the cochlea of echolocating bats is enlarged relative to other skull structures.
Most previous works focused on morphological comparison between echolocating bats and non-echolocating bats. It has been documented that gene expression differences might be important contributors to echolocation-specific features [29, 30]. We compared the gene expression divergence between M. ricketti and C. sphinx, and aimed to put these expression differences in context by comparing them to the high-frequency hearing traits. Of the two bats we examined, the inner ear of M. ricketti was found to have higher expression level of hearing related genes and inner ear morphogenesis genes, such as TMC1 gene, which is probably an adaptation to their developed auditory tuning to their high-frequency calls. TMC1 gene is involved in hair cell structure and function, and show evidences of positive selection exerted on the some echolocating bat species [23, 24]. Thus, the study of the development of bat inner ear is important to understand the molecular mechanisms underlying the generation of ear with high-frequency hearing ability. We previously demonstrated the importance of Prestin gene in high frequency hearing mammals [31–33]. However, Prestin was shown a low expression value in both species, and no significant expression difference was found between M. ricketti and C. sphinx in this work (13.4 vs. 10.7, P-value = 0.21). Although the sequence evolution of Prestin contributes to the evolution of bat hearing, the result shows that echolocation of bats isn’t associated with the transcription level of Prestin gene.
It has been realized that stabilizing selection is likely to be the dominant signature of expression evolution . Although the expression evolution between two species involves different habitat use or ecological niches, the mode of expression divergence is largely explained by neutral evolution and not of direct adaptive significance. Moreover, the neutral view of gene expression divergence asserts that the functionally irrelevant component of gene expression evolves neutrally. Based on this scenario, the vast majority of gene expression divergences between two species should reflect neutral (as opposed to adaptive) variations. Based on this principle, most of expression changes between two species are likely to be of no direct significance. However, high-throughput expression analyses of the inner ear have shown that abnormal events of gene expression of deafness genes in the inner ear is associated with a varying degree of hearing loss . Furthermore, the finding of the enlarged cochlea relative to other skull structures in the echolocating bat can be well explained by the elevated expression of cochlea/inner ear morphogenesis genes. Although we still cannot accurately describe a more precise expression evolution trajectory history at the current form, our work provides a starting point for experimental follow-up.
The study of inner ear gene expression divergences between M. ricketti and C. sphinx using next generation sequencing technology revealed the extent of inner ear transcriptome evolution between the echolocating and non-echolocating bat, and identified a number of novel candidate genes associated with the echolocation ability. Revealing the difference of auditory systems between echolocating bats and non-echolocating bats not only provided better insight into understanding the causes of high-frequency sound hearing ability of echolocating bats, but has also opened an opportunity into exploring how the echolocation evolved in bats.
Sampling, RNA extraction and sequencing
All procedures were in accordance with the guidelines of Regulations for the Administration of Laboratory Animals (Decree No. 2 of the State Science and Technology Commission of the People’s Republic of China on November 14, 1988) approved by the Animal Ethics Committee of East China Normal University (ID no: 20090219). We captured the Rickett’s big-footed bats from a cave (39°42′N, 115°43′E) in Beijing in Oct. 2009, and the greater short-nosed fruit bats were captured during the period of Oct. 2009 from Yuexiu park (23°08′N, 113°20′E) located at Guangdong province, China. The inner ear of the Rickett’s big-footed bats and the Greater short-nosed fruit bats were collected. All tissues were flash frozen in liquid nitrogen and placed in a -80°C freezer until processed for total RNA isolation. Total RNA was isolated using TRIzol (Life Technologies Corp., Carlsbad, CA, USA) according to the manufacturer’s protocols and cleaned up using the RNeasy mini kit (Qiagen, Valencia, CA, USA). RNA samples were quantified by the 2100 Bioanalyzer (Agilent Technologies). We purified mRNA using RNA-Seq sample preparation Kit (Illumina, San Diego, CA). Four paired-end cDNA libraries of each tissue were generated using mRNA-Seq assay for transcriptome sequencing on Illumina Genome Analyzer II platform. Short sequence reads of 75 bp were generated. All these data have been deposited into the NCBI Sequence Read Archive database (SRA run accession numbers: M. ricketti: # SRR837386 and C. sphinx: #SRR837385).
At first, we removed the low quality reads prior to analyzing the data. Two criteria were used in this filtering step: removing reads with adaptors; removing reads with unknown ‘N’ bases. All subsequent analyses were based on these filtered reads. Next, de novo sequence assembly was carried out using Trinity software  designed for short read sequences assembly with default parameters. Only contigs with length greater than 200 bp were used for further analysis. To lower the redundancy in the dataset, low-coverage artifacts or redundancies were removed by using CD-HIT  with an identity threshold of 95%. The detailed work flow is described in Figure 1. To assess the quality of our assemblies, we downloaded all cDNA sequences (820 in total) of M. ricketti and C. sphinx from NCBI Genbank database, we searched the online Genbank database (non-redundant nt database) using species names and downloaded them in bulk (downloaded on April 13, 2013).
Functional annotation and identification of orthologous
Assembled contigs were annotated by using the best hits of BLASTX search against the non-redundant UniProt database with an E-value cutoff of 1E-6 for the annotation of these protein coding contigs that were conserved with other species. Next, we performed a pairwise multiple alignment was performed and the contigs with their the percentage of identity lower than 50% were discarded. The results of the best blast hits were extracted, and the open reading frames were subsequently determined. Next, a reciprocal best BLASTP search was conducted. The Gene Ontology (GO) categories for the non-redundant UniProt proteins (Release 2013 05) was used to assign the GO terms to these transcripts. In our study, only experiment based GO annotations were applied, which can provide a better quality data.
After the differentially expressed genes were classified into different GO categories, an in-house fisher’s exact test program was used to map differentially expressed genes to GO terms. The calculated P values were corrected through Bonferroni correction, and the corrected P values (false discovery rate, F.D.R) were taken as the thresholds of significance.
Measurement of gene expression
To estimate gene expression level, we measured the reads number derived from each contig or isoform using RSEM software package . Next, EBSeq method nested in RSEM package was performed to detect the differentially expressed genes between these two bat species. EBSeq is a Bayesian hierarchical model for the inference on the differential gene expression based on RNA-Seq data . We adjusted the significant level by a correction for false discovery rate (FDR) at < 0.01 with Benjamini-Hochberg correction method .
Estimation of substitution rates
Pair-wise and multiple alignments were generated for two bat species using Mafft software  based on protein sequences and back –translated to DNA sequences. We estimated the overall substitution rates (non-synonymous substitution to synonymous substitution, d N /d S ) using a maximum likelihood method implemented in the CODEML program nested in PAML package version 4.1 . To minimize statistical artefacts from short sequences and saturation effects in d S value, we excluded all alignments that were shorter than 100 bp or that d S value larger than 2 from the analysis.
To validate the expression patterns observed from our RNA-seq analyses, three replicates were analyzed per species using quantitative polymerase chain reaction (qPCR) using the SYBR Prime-Script RT-PCR Kit (TaKaRa) on Applied Biosystems 7300 Real-Time PCR System (Applied Biosystems). We sampled the inner ears from three bat species (M. ricketti, M. davidii, C. sphinx) and two rodents (M. musculus and R. norvegicus). Total RNA was extracted using TRIzol Reagent (Invitrogen) and treated with DNase I (Roche). PCR products were purified and sequenced in both directions with an ABI 3730 DNA sequencer (Applied Biosystems). We designed gene specific primers of TMC1 gene for qPCR. Student’s t-test was performed to examine the difference of gene expression between two different samples. Sequences of forward and reverse primers for normal PCR and qPCR are shown in Additional file 1: Table S1.
We thank Dr. Bin Shen and Dr. Yi-Hsuan Pan for her helpful comments and discussion. We also thank Jinyao Zeng, Junpeng Zhang, Libiao Zhang and Jinshuo Zhang for the help of field work. This work is supported by Key Construction Program of the National ‘985′ project of East China Normal University to Dong Dong (Grant No. is 79633006) and National Natural Science Foundation of China to Dong Dong (Grant No. is 44106480).
- 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 (9): 1869-1886. 10.1093/molbev/msi180.View ArticlePubMedGoogle Scholar
- Jones KE, Purvis A, MacLarnon A, Bininda-Emonds ORP, Simmons NB: A phylogenetic supertree of the bats (Mammalia : Chiroptera). Biol Rev. 2002, 77 (2): 223-259.View ArticlePubMedGoogle Scholar
- Murphy WJ, Eizirik E, O’Brien SJ, Madsen O, Scally M, Douady CJ, Teeling E, Ryder OA, Stanhope MJ, de Jong WW, et al: Resolution of the early placental mammal radiation using Bayesian phylogenetics. Science. 2001, 294 (5550): 2348-2351. 10.1126/science.1067179.View ArticlePubMedGoogle Scholar
- Speakman JR: The evolution of flight and echolocation in bats: another leap in the dark. Mammal Rev. 2001, 31 (2): 111-130. 10.1046/j.1365-2907.2001.00082.x.View ArticleGoogle Scholar
- Vater M, Kossl M: Comparative aspects of cochlear functional organization in mammals. Hearing Res. 2011, 273 (1–2): 89-99.View ArticleGoogle Scholar
- Jones G, Teeling EC: The evolution of echolocation in bats. Trends Ecol Evol. 2006, 21 (3): 149-156. 10.1016/j.tree.2006.01.001.View ArticlePubMedGoogle Scholar
- Teeling EC: Hear, hear: the convergent evolution of echolocation in bats?. Trends Ecol Evol. 2009, 24 (7): 351-354. 10.1016/j.tree.2009.02.012.View ArticlePubMedGoogle Scholar
- Altringham JD: Bats: Biology and Behaviour. 1996, Oxford: Oxford University PressGoogle Scholar
- Davies KT, Maryanto I, Rossiter SJ: Evolutionary origins of ultrasonic hearing and laryngeal echolocation in bats inferred from morphological analyses of the inner ear. Front Zool. 2013, 10 (1): 2-10.1186/1742-9994-10-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Pye A: The structure of the cochlea in chiroptera. I. Microchiroptera. Emballonuroidea and Rhinolophoidea. J Morphol. 1966, 118 (4): 495-510. 10.1002/jmor.1051180404.View ArticlePubMedGoogle Scholar
- Kossl M, Mayer F, Frank G, Faulstich M, Russell IJ: Evolutionary adaptations of cochlear function in Jamaican mormoopid bats. J Comp Physiol A. 1999, 185 (3): 217-228. 10.1007/s003590050381.View ArticlePubMedGoogle Scholar
- Carroll SB: Evolution at two levels: on genes and form. Plos Biol. 2005, 3 (7): 1159-1166.View ArticleGoogle Scholar
- Ekblom R, Galindo J: Applications of next generation sequencing in molecular ecology of non-model organisms. Heredity (Edinb). 2011, 107 (1): 1-15. 10.1038/hdy.2010.152.View ArticleGoogle Scholar
- Hudson ME: Sequencing breakthroughs for genomic ecology and evolutionary biology. Mol Ecol Resour. 2008, 8 (1): 3-17. 10.1111/j.1471-8286.2007.02019.x.View ArticlePubMedGoogle Scholar
- 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 (5709): 580-584. 10.1126/science.1105113.View ArticlePubMedGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng QD, et al: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29 (7): 644-10.1038/nbt.1883. U130PubMed CentralView ArticlePubMedGoogle Scholar
- Katoh K, Misawa K, Kuma K, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30 (14): 3059-3066. 10.1093/nar/gkf436.PubMed CentralView ArticlePubMedGoogle Scholar
- Li B, Dewey CN: RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011, 12: 323-10.1186/1471-2105-12-323.PubMed CentralView ArticlePubMedGoogle Scholar
- Leng N, Dawson JA, Thomson JA, Ruotti V, Rissman AI, Smits BM, Haag JD, Gould MN, Stewart RM, Kendziorski C: EBSeq: an empirical Bayes hierarchical model for inference in RNA-seq experiments. Bioinformatics. 2013, 29 (8): 1035-1043. 10.1093/bioinformatics/btt087.PubMed CentralView ArticlePubMedGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMed CentralView ArticlePubMedGoogle Scholar
- Novacek MJ: Evidence for echolocation in the oldest known bats. Nature. 1985, 315 (6015): 140-141. 10.1038/315140a0.View ArticlePubMedGoogle Scholar
- Marcotti W, Erven A, Johnson SL, Steel KP, Kros CJ: Tmc1 is necessary for normal functional maturation and survival of inner and outer hair cells in the mouse cochlea. J Physiol Lond. 2006, 574 (3): 677-698. 10.1113/jphysiol.2005.095661.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (5): 480-489. 10.1038/hdy.2011.119.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang G, Cowled C, Shi Z, Huang Z, Bishop-Lilly KA, Fang X, Wynne JW, Xiong Z, Baker ML, Zhao W, et al: Comparative analysis of bat genomes provides insight into the evolution of flight and immunity. Science. 2013, 339 (6118): 456-460. 10.1126/science.1230835.View ArticlePubMedGoogle Scholar
- Ohno S: Evolution by Gene Duplication. 1970, New York: Springer-VelagView ArticleGoogle Scholar
- 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 (11): 6241-6246. 10.1073/pnas.111551998.PubMed CentralView ArticlePubMedGoogle Scholar
- Simmons NB, Seymour KL, Habersetzer J, Gunnell GF: Primitive Early Eocene bat from Wyoming and the evolution of flight and echolocation. Nature. 2008, 451 (7180): 818-821. 10.1038/nature06549.View ArticlePubMedGoogle Scholar
- Chatterjee S, Kraus P, Lufkin T: A symphony of inner ear developmental control genes. BMC Genet. 2010, 11: 68-10.1186/1471-2156-11-68.PubMed CentralView ArticlePubMedGoogle Scholar
- Schwartz CP, Smotherman MS: Mapping vocalization-related immediate early gene expression in echolocating bats. Behav Brain Res. 2011, 224 (2): 358-368. 10.1016/j.bbr.2011.06.023.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (6): e1002788-10.1371/journal.pgen.1002788.PubMed CentralView ArticlePubMedGoogle Scholar
- Li G, Wang JH, Rossiter SJ, Jones G, Cottont JA, Zhang SY: The hearing gene Prestin reunites echolocating bats. Proc Natl Acad Sci USA. 2008, 105 (37): 13959-13964. 10.1073/pnas.0802097105.PubMed CentralView ArticlePubMedGoogle Scholar
- Li Y, Liu Z, Shi P, Zhang JZ: The hearing gene Prestin unites echolocating bats and whales. Curr Biol. 2010, 20 (2): R55-R56. 10.1016/j.cub.2009.11.042.View ArticlePubMedGoogle Scholar
- Liu Y, Cotton JA, Shen B, Han XQ, Rossiter SJ, Zhang SY: Convergent sequence evolution between echolocating bats and dolphins. Curr Biol. 2010, 20 (2): R53-R54. 10.1016/j.cub.2009.11.058.View ArticlePubMedGoogle Scholar
- Gilad Y, Oshlack A, Rifkin SA: Natural selection on gene expression. Trends Genet. 2006, 22 (8): 456-461. 10.1016/j.tig.2006.06.002.View ArticlePubMedGoogle Scholar
- Hertzano R, Elkon R: High throughput gene expression analysis of the inner ear. Hear Res. 2012, 288 (1–2): 77-88.View ArticlePubMedGoogle Scholar
- Li WZ, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22 (13): 1658-1659. 10.1093/bioinformatics/btl158.View ArticlePubMedGoogle Scholar
- Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I: Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001, 125 (1–2): 279-284.View ArticlePubMedGoogle Scholar
- Yang ZH: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.View ArticlePubMedGoogle Scholar
This article is published under license to 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 cited.