Comparative inner ear transcriptome analysis between the Rickett's big-footed bats (Myotis ricketti) and the greater short-nosed fruit bats (Cynopterus sphinx)

Sex Differences in Cochlear Transcriptomes in Horseshoe Bats

Sexual dimorphism of calls is common in animals, whereas studies on the molecular basis underlying this phenotypic variation are still scarce. In this study, we used comparative transcriptomics of cochlea to investigate the sex-related difference in gene expression and alternative splicing in four Rhinolophus taxa. Based on 31 cochlear transcriptomes, we performed differential gene expression (DGE) and alternative splicing (AS) analyses between the sexes in each taxon. Consistent with the degree of difference in the echolocation pulse frequency between the sexes across the four taxa, we identified the largest number of differentially expressed genes (DEGs) and alternatively spliced genes (ASGs) in R. sinicus. However, we also detected multiple DEGs and ASGs in taxa without sexual differences in echolocation pulse frequency, suggesting that these genes might be related to other parameters of echolocation pulse rather than the frequency component. Some DEGs and ASGs are related to hearing loss or deafness genes in human or mice and they can be considered to be candidates associated with the sexual differences of echolocation pulse in bats. We also detected more than the expected overlap of DEGs and ASGs in two taxa. Overall, our current study supports the important roles of both DGE and AS in generating or maintaining sexual differences in animals.

Effects of Climate Change on the Distribution of Threatened Fishing Bat Myotis pilosus in China

Climate change and biodiversity loss are two severe challenges that the world is facing. Studying the distribution shifts of species in response to climate change could provide insights into long-term conservation and biodiversity maintenance. Myotis pilosus is the only known fishing bat in East Asia, whereas its population has been decreasing in recent years and it is listed as a "Vulnerable" species. To assess the impact of climate change on the distribution of M. pilosus, we obtained 33 M. pilosus occurrence records within China where they are mainly distributed, and extracted 30 environmental variables. MaxEnt was applied to assess the habitat suitability, recognize the important environmental variables, predict future distribution changes, and identify the potential future climate refugia. The prediction result based on eleven dominant environmental variables was excellent. The Jackknife test showed that the "minimum temperature of coldest month", "precipitation of wettest quarter", "percent tree cover", and "precipitation of driest month" were the main factors affecting the distribution of M. pilosus. The current suitable areas were predicted to be mainly located in southwest and southeast China with a total area of about 160.54 × 104 km2, accounting for 16.72% of China's land area. Based on the CCSM4, it was predicted that the future (2050 and 2070) suitable areas of M. pilosus will expand and shift to high latitudes and altitudes with global warming, but the area of moderately and highly suitable habitats will be small. Considering the dispersal capacity of M. pilosus, the area of colonized suitable habitats in 2050 and 2070 was predicted to be only ca. 94 × 104 km2 and 155 × 104 km2, respectively. The central and southern parts of Hainan, southern Guangdong, central Guizhou, and southern Beijing were identified as potential climate refugia and could be considered as priority conservation areas for M. pilosus. Thus, we suggest long-term monitoring of the priority conservation areas, especially the areas at high latitudes and altitudes. These results contribute to our knowledge of the possible spatial distribution pattern of M. pilosus under current and future climate scenarios, which is important for the population protection and habitat management of this special piscivorous bat species.

Untargeted metabolomics of the cochleae from two laryngeally echolocating bats

High-frequency hearing is regarded as one of the most functionally important traits in laryngeally echolocating bats. Abundant candidate hearing-related genes have been identified to be the important genetic bases underlying high-frequency hearing for laryngeally echolocating bats, however, extensive metabolites presented in the cochleae have not been studied. In this study, we identified 4,717 annotated metabolites in the cochleae of two typical laryngeally echolocating bats using the liquid chromatography-mass spectroscopy technology, metabolites classified as amino acids, peptides, and fatty acid esters were identified as the most abundant in the cochleae of these two echolocating bat species, Rhinolophus sinicus and Vespertilio sinensis. Furthermore, 357 metabolites were identified as significant differentially accumulated (adjusted p-value <0.05) in the cochleae of these two bat species with distinct echolocating dominant frequencies. Downstream KEGG enrichment analyses indicated that multiple biological processes, including signaling pathways, nervous system, and metabolic process, were putatively different in the cochleae of R. sinicus and V. sinensis. For the first time, this study investigated the extensive metabolites and associated biological pathways in the cochleae of two laryngeal echolocating bats and expanded our knowledge of the metabolic molecular bases underlying high-frequency hearing in the cochleae of echolocating bats.

Full-Length Transcriptome of Myotis pilosus as a Reference Resource and Mining of Auditory and Immune Related Genes

Rickett's big-footed bat, Myotis pilosus, which belongs to the family Vespertilionida, is the only known piscivorous bat in East Asia. Accurate whole genome and transcriptome annotations are essential for the study of bat biological evolution. The lack of a whole genome for M. pilosus has limited our understanding of the molecular mechanisms underlying the species' evolution, echolocation, and immune response. In the present work, we sequenced the entire transcriptome using error-corrected PacBio single-molecule real-time (SMRT) data. Then, a total of 40 GB of subreads were generated, including 29,991 full-length non-chimeric (FLNC) sequences. After correction by Illumina short reads and de-redundancy, we obtained 26,717 error-corrected isoforms with an average length of 3018.91 bp and an N50 length of 3447 bp. A total of 1528 alternative splicing (AS) events were detected by transcriptome structural analysis. Furthermore, 1032 putative transcription factors (TFs) were identified, with additional identification of several long non-coding RNAs (lncRNAs) with high confidence. Moreover, several key genes, including PRL-2, DPP4, Glul, and ND1 were also identified as being associated with metabolism, immunity, nervous system processes, and auditory perception. A multitude of pattern recognition receptors was identified, including NLR, RLR, SRCR, the antiviral molecule IRF3, and the IFN receptor subunit IFNAR1. High-quality reference genomes at the transcriptome level may be used to quantify gene or transcript expression, evaluate alternative splicing levels, identify novel transcripts, and enhance genome annotation in bats.

Transcriptome Analyses Provide Insights into the Auditory Function in Trachemys scripta elegans

An auditory ability is essential for communication in vertebrates, and considerable attention has been paid to auditory sensitivity in mammals, birds, and frogs. Turtles were thought to be deaf for a long time; however, recent studies have confirmed the presence of an auditory ability in Trachemys scripta elegans as well as sex-related differences in hearing sensitivity. Earlier studies mainly focused on the morphological and physiological functions of the hearing organ in turtles; thus, the gene expression patterns remain unclear. In this study, 36 transcriptomes from six tissues (inner ear, tympanic membrane, brain, eye, lung, and muscle) were sequenced to explore the gene expression patterns of the hearing system in T. scripta elegans. A weighted gene co-expression network analysis revealed that hub genes related to the inner ear and tympanic membrane are involved in development and signal transduction. Moreover, we identified six differently expressed genes (GABRA1, GABRG2, GABBR2, GNAO1, SLC38A1, and SLC12A5) related to the GABAergic synapse pathway as candidate genes to explain the differences in sexually dimorphic hearing sensitivity. Collectively, this study provides a critical foundation for genetic research on auditory functions in turtles.

Comparative transcriptome analysis provides insights into the molecular mechanisms of high-frequency hearing differences between the sexes of Odorrana tormota

Background

Acoustic communication is important for the survival and reproduction of anurans and masking background noise is a critical factor for their effective acoustic communication. Males of the concave-eared frog (Odorrana tormota) have evolved an ultrasonic communication capacity to avoid masking by the widespread background noise of local fast-flowing streams, whereas females exhibit no ultrasonic sensitivity. However, the molecular mechanisms underlying the high-frequency hearing differences between the sexes of O. tormota are still poorly understood.

Results

In this study, we sequenced the brain transcriptomes of male and female O. tormota, and compared their differential gene expression. A total of 4,605 differentially expressed genes (DEGs) between the sexes of O. tormota were identified and eleven of them were related to auditory based on the annotation and enrichment analysis. Most of these DEGs in males showed a higher expression trend than females in both quantity and expression quantity. The highly expressed genes in males were relatively concentrated in neurogenesis, signal transduction, ion transport and energy metabolism, whereas the up-expressed genes in females were mainly related to the growth and development regulation of specific auditory cells.

Conclusions

The transcriptome of male and female O. tormota has been sequenced and de novo assembled, which will provide gene reference for further genomic studies. In addition, this is the first research to reveal the molecular mechanisms of sex differences in ultrasonic hearing between the sexes of O. tormota and will provide new insights into the genetic basis of the auditory adaptation in amphibians during their transition from water to land.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08536-2.

Novel de Novo Genome of Cynopterus brachyotis Reveals Evolutionarily Abrupt Shifts in Gene Family Composition across Fruit Bats

Major novel physiological or phenotypic adaptations often require accompanying modifications at the genic level. Conversely, the detection of considerable contractions and/or expansions of gene families can be an indicator of fundamental but unrecognized physiological change. We sequenced a novel fruit bat genome (Cynopterus brachyotis) and adopted a comparative approach to reconstruct the evolution of fruit bats, mapping contractions and expansions of gene families along their evolutionary history. Despite a radical change in life history as compared with other bats (e.g., loss of echolocation, large size, and frugivory), fruit bats have undergone surprisingly limited change in their genic composition, perhaps apart from a potentially novel gene family expansion relating to telomere protection and longevity. In sharp contrast, within fruit bats, the new Cynopterus genome bears the signal of unusual gene loss and gene family contraction, despite its similar morphology and lifestyle to two other major fruit bat lineages. Most missing genes are regulatory, immune-related, and olfactory in nature, illustrating the diversity of genomic strategies employed by bats to contend with responses to viral infection and olfactory requirements. Our results underscore that significant fluctuations in gene family composition are not always associated with obvious examples of novel physiological and phenotypic adaptations but may often relate to less-obvious shifts in immune strategies.

Transcriptome sequencing of cochleae from constant-frequency and frequency-modulated echolocating bats

Echolocating bats are fascinating for their ability to ‘see’ the world in the darkness. Ultrahigh frequency hearing is essential for echolocation. In this study we collected cochlear tissues from constant-frequency (CF) bats (two subspecies of Rhinolophus affinis, Rhinolophidae) and frequency-modulated (FM) bats (Myotis ricketti, Vespertilionidae) and applied PacBio single-molecule real-time isoform sequencing (Iso-seq) technology to generate the full-length (FL) transcriptomes for the three taxa. In total of 10103, 9676 and 10504 non-redundant FL transcripts for R. a. hainanus, R. a. himalayanus and Myotis ricketti were obtained respectively. These data present a comprehensive list of transcripts involved in ultrahigh frequency hearing of echolocating bats including 26342 FL transcripts, 24833 of which are annotated by public databases. No further comparative analyses were performed on the current data in this study. This data can be reused to quantify gene or transcript expression, assess the level of alternative splicing, identify novel transcripts and improve genome annotation of bat species.

Measurement(s)	cochlea • transcriptome • sequence feature annotation
Technology Type(s)	isoform sequencing • sequence annotation
Factor Type(s)	species
Sample Characteristic - Organism	Rhinolophus affinis • Myotis ricketti

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12999614

Echolocation call frequency variation in horseshoe bats: molecular basis revealed by comparative transcriptomics

Recently diverged taxa with contrasting phenotypes offer opportunities for unravelling the genetic basis of phenotypic variation in nature. Horseshoe bats are a speciose group that exhibit a derived form of high-duty cycle echolocation in which the inner ear is finely tuned to echoes of the narrowband call frequency. Here, by focusing on three recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis) that display divergent echolocation call frequencies, we aim to identify candidate loci putatively involved in hearing frequency variation. We used de novo transcriptome sequencing of two mainland taxa (himalayanus and macrurus) and one island taxon (hainanus) to compare expression profiles of thousands of genes. By comparing taxa with divergent call frequencies (around 15 kHz difference), we identified 252 differentially expressed genes, of which six have been shown to be involved in hearing or deafness in human/mouse. To obtain further validation of these results, we applied quantitative reverse transcription-PCR to the candidate gene FBXL15 and found a broad association between the level of expression and call frequency across taxa. The genes identified here represent strong candidate loci associated with hearing frequency variation in bats.

Positive Selection and Inactivation in the Vision and Hearing Genes of Cetaceans

The transition to an aquatic lifestyle in cetaceans (whales and dolphins) resulted in a radical transformation in their sensory systems. Toothed whales acquired specialized high-frequency hearing tied to the evolution of echolocation, whereas baleen whales evolved low-frequency hearing. More generally, all cetaceans show adaptations for hearing and seeing underwater. To determine the extent to which these phenotypic changes have been driven by molecular adaptation, we performed large-scale targeted sequence capture of 179 sensory genes across the Cetacea, incorporating up to 54 cetacean species from all major clades as well as their closest relatives, the hippopotamuses. We screened for positive selection in 167 loci related to vision and hearing and found that the diversification of cetaceans has been accompanied by pervasive molecular adaptations in both sets of genes, including several loci implicated in nonsyndromic hearing loss. Despite these findings, however, we found no direct evidence of positive selection at the base of odontocetes coinciding with the origin of echolocation, as found in studies examining fewer taxa. By using contingency tables incorporating taxon- and gene-based controls, we show that, although numbers of positively selected hearing and nonsyndromic hearing loss genes are disproportionately high in cetaceans, counts of vision genes do not differ significantly from expected values. Alongside these adaptive changes, we find increased evidence of pseudogenization of genes involved in cone-mediated vision in mysticetes and deep-diving odontocetes.

Gene expression vs. sequence divergence: comparative transcriptome sequencing among natural Rhinolophus ferrumequinum populations with different acoustic phenotypes

BACKGROUND

Although the sensory drive hypothesis can explain the geographic variation in echolocation frequencies of some bat species, the molecular mechanisms underlying this phenomenon are still unclear. The three lineages of greater horseshoe bat (Rhinolophus ferrumequinum) in China (northeast, central-east, and southwest) have significant geographic variation in resting frequencies (RF) of echolocation calls. Because their cochleae have an acoustic fovea that is highly sensitive to a narrow range of frequencies, we reported the transcriptomes of cochleae collected from three genetic lineages of R. ferrumequinum, which is an ideal organism for studying geographic variation in echolocation signals, and tried to understand the mechanisms behind this bat phenomenon by analyzing gene expression and sequence variation.

RESULTS

A total of 8190 differentially expressed genes (DEGs) were identified. We identified five modules from all DEGs that were significantly related to RF or forearm length (FL). DEGs in the RF-related modules were significantly enriched in the gene categories involved in neural activity, learning, and response to sound. DEGs in the FL-related modules were significantly enriched in the pathways related to muscle and actin functions. Using 21,945 single nucleotide polymorphisms, we identified 18 candidate unigenes associated with hearing, five of which were differentially expressed among the three populations. Additionally, the gene ERBB4, which regulates diverse cellular processes in the inner ear such as cell proliferation and differentiation, was in the largest module. We also found 49 unigenes that were under positive selection from 4105 one-to-one orthologous gene pairs between the three R. ferrumequinum lineages and three other Chiroptera species.

CONCLUSIONS

The variability of gene expression and sequence divergence at the molecular level might provide evidence that can help elucidate the genetic basis of geographic variation in echolocation signals of greater horseshoe bats.

A metaanalysis of bat phylogenetics and positive selection based on genomes and transcriptomes from 18 species

This work represents a large, order-wide evolutionary analysis of the order Chiroptera (bats). Our pipeline for assembling sequence data and curating orthologous multiple sequence alignments includes methods for improving results when combining genomic and transcriptomic data sources. The resulting phylogenetic tree divides the order Chiroptera into Yinpterochiroptera and Yangochiroptera, in disagreement with the previous division into Megachiroptera and Microchiroptera and in agreement with some other recent molecular studies, and also provides evidence for other contested branch placements. We also performed a genome-wide analysis of positive selection and found 181 genes with signatures of positive selection. Enrichment analysis shows these positively selected genes to be primarily related to immune responses but also, surprisingly, collagen formation.

Historically, the evolution of bats has been analyzed using a small number of genetic loci for many species or many genetic loci for a few species. Here we present a phylogeny of 18 bat species, each of which is represented in 1,107 orthologous gene alignments used to build the tree. We generated a transcriptome sequence of Hypsignathus monstrosus, the African hammer-headed bat, and additional transcriptome sequence for Rousettus aegyptiacus, the Egyptian fruit bat. We then combined these data with existing genomic and transcriptomic data from 16 other bat species. In the analysis of such datasets, there is no clear consensus on the most reliable computational methods for the curation of quality multiple sequence alignments since these public datasets represent multiple investigators and methods, including different source materials (chromosomal DNA or expressed RNA). Here we lay out a systematic analysis of parameters and produce an advanced pipeline for curating orthologous gene alignments from combined transcriptomic and genomic data, including a software package: the Mismatching Isoform eXon Remover (MIXR). Using this method, we created alignments of 11,677 bat genes, 1,107 of which contain orthologs from all 18 species. Using the orthologous gene alignments created, we assessed bat phylogeny and also performed a holistic analysis of positive selection acting in bat genomes. We found that 181 genes have been subject to positive natural selection. This list is dominated by genes involved in immune responses and genes involved in the production of collagens.

De Novo Transcriptome Assembly and Functional Annotation in Five Species of Bats

High-throughput RNA sequencing is a powerful tool that allows us to perform gene prediction and analyze tissue-specific overexpression of genes, but also at species level comparisons can be performed, although in a more restricted manner. In the present study complete liver transcriptomes of five tropical bat species were De novo assembled and annotated. Highly expressed genes in the five species were involved in glycolysis and lipid metabolism pathways. Cross-species differential expression analysis was conducted using single copy orthologues shared across the five species. Between 22 and 29 orthologs were upregulated for each species. We detected upregulated expression in Artibeus jamaicensis genes related to fructose metabolism pathway. Such findings can be correlated with A. jamaicensis dietary habits, as it was the unique frugivorous species included. This is the first report of transcriptome assembly by RNA-seq in these species, except for A. jamaicensis and as far as our knowledge is the first cross-species comparisons of transcriptomes and gene expression in tropical bats.

Comparative cochlear transcriptomics of echolocating bats provides new insights into different nervous activities of CF bat species

The molecular mechanisms used by echolocating bats to deal with different ultrasonic signals remain to be revealed. Here, we utilised RNA-Seq data to conduct comparative cochlear transcriptomics to assess the variation of gene expression among bats with three types of echolocation: constant-frequency (CF) bats, frequency-modulated (FM) bats and click bats. Our results suggest larger differences in gene expression between CF and click bats than between CF and FM bats and small differences between FM and click bats. We identified 426 and 1,504 differentially expressed genes (DEGs) by the different methods as functionally important for CF bats, in that they showed consistent upregulation in the cochlea of two CF bats, relative to the levels in click and FM bats. Subsequently, downstream GO and KEGG enrichment analyses indicated that both the 426 and 1,504 gene sets were associated with changes in nervous activities in the cochleae of CF bats. In addition, another set of 1,764 DEGs were identified to have crucial hearing related physiological functions for laryngeally echolocating bats. Our study provides a comprehensive overview of the genetic basis of differences among echolocating bats, revealing different nervous system activities during auditory perception in the cochlea particularly in CF bats.

Age-dependent gene expression in the inner ear of big brown bats (Eptesicus fuscus)

For echolocating bats, hearing is essential for survival. Specializations for detecting and processing high frequency sounds are apparent throughout their auditory systems. Recent studies on echolocating mammals have reported evidence of parallel evolution in some hearing-related genes in which distantly related groups of echolocating animals (bats and toothed whales), cluster together in gene trees due to apparent amino acid convergence. However, molecular adaptations can occur not only in coding sequences, but also in the regulation of gene expression. The aim of this study was to examine the expression of hearing-related genes in the inner ear of developing big brown bats, Eptesicus fuscus, during the period in which echolocation vocalizations increase dramatically in frequency. We found that seven genes were significantly upregulated in juveniles relative to adults, and that the expression of four genes through development correlated with estimated age. Compared to available data for mice, it appears that expression of some hearing genes is extended in juvenile bats. These results are consistent with a prolonged growth period required to develop larger cochlea relative to body size, a later maturation of high frequency hearing, and a greater dependence on high frequency hearing in echolocating bats.

Transcriptomic Signatures of Tacaribe Virus-Infected Jamaican Fruit Bats

As reservoir hosts of viruses associated with human disease, little is known about the interactions between bats and viruses. Using Jamaican fruit bats infected with Tacaribe virus (TCRV) as a model, we characterized the gene expression responses to infection in different tissues and identified pathways involved with the response to infection. This report is the most detailed gene discovery work in the species to date and the first to describe immune gene expression responses in bats during a pathogenic viral infection.

Tacaribe virus (TCRV) is a mammalian arenavirus that was first isolated from artibeus bats in the 1950s. Subsequent experimental infection of Jamaican fruit bats (Artibeus jamaicensis) caused a disease similar to that of naturally infected bats. Although substantial attention has focused on bats as reservoir hosts of viruses that cause human disease, little is known about the interactions between bats and their pathogens. We performed a transcriptome-wide study to illuminate the response of Jamaican fruit bats experimentally infected with TCRV. Differential gene expression analysis of multiple tissues revealed global and organ-specific responses associated with innate antiviral responses, including interferon alpha/beta and Toll-like receptor signaling, activation of complement cascades, and cytokine signaling, among others. Genes encoding proteins involved in adaptive immune responses, such as gamma interferon signaling and costimulation of T cells by the CD28 family, were also altered in response to TCRV infection. Immunoglobulin gene expression was also elevated in the spleens of infected bats, including IgG, IgA, and IgE isotypes. These results indicate an active innate and adaptive immune response to TCRV infection occurred but did not prevent fatal disease. This de novo assembly provides a high-throughput data set of the Jamaican fruit bat and its host response to TCRV infection, which remains a valuable tool to understand the molecular signatures involved in antiviral responses in bats.

IMPORTANCE As reservoir hosts of viruses associated with human disease, little is known about the interactions between bats and viruses. Using Jamaican fruit bats infected with Tacaribe virus (TCRV) as a model, we characterized the gene expression responses to infection in different tissues and identified pathways involved with the response to infection. This report is the most detailed gene discovery work in the species to date and the first to describe immune gene expression responses in bats during a pathogenic viral infection.

Blood miRNomes and transcriptomes reveal novel longevity mechanisms in the long-lived bat, Myotis myotis

BACKGROUND

Chiroptera, the bats, are the only order of mammals capable of true self-powered flight. Bats exhibit a number of other exceptional traits such as echolocation, viral tolerance and, perhaps most puzzlingly, extreme longevity given their body size. Little is known about the molecular mechanisms driving their extended longevity particularly at the levels of gene expression and post-transcriptional regulation. To elucidate the molecular mechanisms that may underlie their unusual longevity, we have deep sequenced 246.5 million small RNA reads from whole blood of the long-lived greater mouse-eared bats, Myotis myotis, and conducted a series of genome-wide comparative analyses between bat and non-bat mammals (human, pig and cow) in both blood miRNomes and transcriptomes, for the first time.

RESULTS

We identified 539 miRNA gene candidates from bats, of which 468 unique mature miRNA were obtained. More than half of these miRNA (65.1 %) were regarded as bat-specific, regulating genes involved in the immune, ageing and tumorigenesis pathways. We have also developed a stringent pipeline for genome-wide miRNome comparisons across species, and identified 37 orthologous miRNA groups shared with bat, human, pig and cow, 6 of which were differentially expressed. For bats, 3 out of 4 up-regulated miRNA (miR-101-3p, miR-16-5p, miR-143-3p) likely function as tumor suppressors against various kinds of cancers, while one down-regulated miRNA (miR-221-5p) acts as a tumorigenesis promoter in human breast and pancreatic cancers. Additionally, a genome-wide comparison of mRNA transcriptomes across species also revealed specific gene expression patterns in bats. 127 up-regulated genes were enriched mainly in mitotic cell cycle and DNA repair mechanisms, while 364 down-regulated genes were involved primarily in mitochondrial activity.

CONCLUSIONS

Our comprehensive and integrative analyses revealed bat-specific and differentially expressed miRNA and mRNA that function in key longevity pathways, producing a distinct bat gene expression pattern. For the first time, we show that bats may possess unique regulatory mechanisms for resisting tumorigenesis, repairing cellular damage and preventing oxidative stresses, all of which likely contribute to the extraordinary lifespan of Myotis myotis.

Bat Accelerated Regions Identify a Bat Forelimb Specific Enhancer in the HoxD Locus

The molecular events leading to the development of the bat wing remain largely unknown, and are thought to be caused, in part, by changes in gene expression during limb development. These expression changes could be instigated by variations in gene regulatory enhancers. Here, we used a comparative genomics approach to identify regions that evolved rapidly in the bat ancestor, but are highly conserved in other vertebrates. We discovered 166 bat accelerated regions (BARs) that overlap H3K27ac and p300 ChIP-seq peaks in developing mouse limbs. Using a mouse enhancer assay, we show that five Myotis lucifugus BARs drive gene expression in the developing mouse limb, with the majority showing differential enhancer activity compared to the mouse orthologous BAR sequences. These include BAR116, which is located telomeric to the HoxD cluster and had robust forelimb expression for the M. lucifugus sequence and no activity for the mouse sequence at embryonic day 12.5. Developing limb expression analysis of Hoxd10-Hoxd13 in Miniopterus natalensis bats showed a high-forelimb weak-hindlimb expression for Hoxd10-Hoxd11, similar to the expression trend observed for M. lucifugus BAR116 in mice, suggesting that it could be involved in the regulation of the bat HoxD complex. Combined, our results highlight novel regulatory regions that could be instrumental for the morphological differences leading to the development of the bat wing.

Divergent evolutionary rates in vertebrate and mammalian specific conserved non-coding elements (CNEs) in echolocating mammals

BACKGROUND

The majority of DNA contained within vertebrate genomes is non-coding, with a certain proportion of this thought to play regulatory roles during development. Conserved Non-coding Elements (CNEs) are an abundant group of putative regulatory sequences that are highly conserved across divergent groups and thus assumed to be under strong selective constraint. Many CNEs may contain regulatory factor binding sites, and their frequent spatial association with key developmental genes - such as those regulating sensory system development - suggests crucial roles in regulating gene expression and cellular patterning. Yet surprisingly little is known about the molecular evolution of CNEs across diverse mammalian taxa or their role in specific phenotypic adaptations. We examined 3,110 vertebrate-specific and ~82,000 mammalian-specific CNEs across 19 and 9 mammalian orders respectively, and tested for changes in the rate of evolution of CNEs located in the proximity of genes underlying the development or functioning of auditory systems. As we focused on CNEs putatively associated with genes underlying the development/functioning of auditory systems, we incorporated echolocating taxa in our dataset because of their highly specialised and derived auditory systems.

RESULTS

Phylogenetic reconstructions of concatenated CNEs broadly recovered accepted mammal relationships despite high levels of sequence conservation. We found that CNE substitution rates were highest in rodents and lowest in primates, consistent with previous findings. Comparisons of CNE substitution rates from several genomic regions containing genes linked to auditory system development and hearing revealed differences between echolocating and non-echolocating taxa. Wider taxonomic sampling of four CNEs associated with the homeobox genes Hmx2 and Hmx3 - which are required for inner ear development - revealed family-wise variation across diverse bat species. Specifically within one family of echolocating bats that utilise frequency-modulated echolocation calls varying widely in frequency and intensity high levels of sequence divergence were found.

CONCLUSIONS

Levels of selective constraint acting on CNEs differed both across genomic locations and taxa, with observed variation in substitution rates of CNEs among bat species. More work is needed to determine whether this variation can be linked to echolocation, and wider taxonomic sampling is necessary to fully document levels of conservation in CNEs across diverse taxa.