The ecological radiation of mammals to inhabit a variety of light environments is largely attributed to adaptive changes in their visual systems. Visual capabilities are conferred by anatomical features of the eyes as well as the combination and properties of their constituent light sensitive pigments. To test whether evolutionary switches to different niches characterized by dim-light conditions coincided with molecular adaptation of the rod pigment rhodopsin, we sequenced the rhodopsin gene in twenty-two mammals including several bats and subterranean mole-rats. We compared these to thirty-seven published mammal rhodopsin sequences, from species with divergent visual ecologies, including nocturnal, diurnal and aquatic groups. All taxa possessed an intact functional rhodopsin; however, phylogenetic tree reconstruction recovered a gene tree in which rodents were not monophyletic, and also in which echolocating bats formed a monophyletic group. These conflicts with the species tree appear to stem from accelerated evolution in these groups, both of which inhabit low light environments. Selection tests confirmed divergent selection pressures in the clades of subterranean rodents and bats, as well as in marine mammals that live in turbid conditions. We also found evidence of divergent selection pressures among groups of bats with different sensory modalities based on vision and echolocation. Sliding window analyses suggest most changes occur in transmembrane domains, particularly obvious within the pinnipeds; however, we found no obvious pattern between photopic niche and predicted spectral sensitivity based on known critical amino acids. This study indicates that the independent evolution of rhodopsin vision in ecologically specialised groups of mammals has involved molecular evolution at the sequence level, though such changes might not mediate spectral sensitivity directly.

The phylogenetic and taxonomic relationships among the Old World leaf-nosed bats (Hipposideridae) and the closely related horseshoe bats (Rhinolophidae) remain unresolved. In this study, we generated a novel approximately 10-kb molecular data set of 19 nuclear exon and intron gene fragments for 40 bat species to elucidate the phylogenetic relationships within the families Rhinolophidae and Hipposideridae. We estimated divergence times and explored potential reasons for any incongruent phylogenetic signal. We demonstrated the effects of outlier taxa and genes on phylogenetic reconstructions and compared the relative performance of intron and exon data to resolve phylogenetic relationships. Phylogenetic analyses produced a well-resolved phylogeny, supporting the familial status of Hipposideridae and demonstrated the paraphyly of the largest genus, Hipposideros. A fossil-calibrated timetree and biogeographical analyses estimated that Rhinolophidae and Hipposideridae diverged in Africa during the Eocene approximately 42 Ma. The phylogram, the timetree, and a unique retrotransposon insertion supported the elevation of the subtribe Rhinonycterina to family level and which is diagnosed herein. Comparative analysis of diversification rates showed that the speciose genera Rhinolophus and Hipposideros underwent diversification during the Mid-Miocene Climatic Optimum. The intron versus exon analyses demonstrated the improved nodal support provided by introns for our optimal tree, an important finding for large-scale phylogenomic studies, which typically rely on exon data alone. With the recent outbreak of Middle East respiratory syndrome, caused by a novel coronavirus, the study of these species is urgent as they are considered the natural reservoir for emergent severe acute respiratory syndrome (SARS)-like coronaviruses. It has been shown that host phylogeny is the primary factor that determines a virus's persistence, replicative ability, and can act as a predictor of new emerging disease. Therefore, this newly resolved phylogeny can be used to direct future assessments of viral diversity and to elucidate the origin and development of SARS-like coronaviruses in mammals.

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.

The Piwi-interacting RNA (piRNA) pathway is responsible for germline specification, gametogenesis, transposon silencing, and genome integrity. Transposable elements can disrupt genome and its functions. However, piRNA pathway evolution and its adaptation to transposon diversity in the teleost fish remain unknown. This article unveils evolutionary scene of piRNA pathway and its association with diverse transposons by systematically comparative analysis on diverse teleost fish genomes. Selective pressure analysis on piRNA pathway and miRNA/siRNA (microRNA/small interfering RNA) pathway genes between teleosts and mammals showed an accelerated evolution of piRNA pathway genes in the teleost lineages, and positive selection on functional PAZ (Piwi/Ago/Zwille) and Tudor domains involved in the Piwi-piRNA/Tudor interaction, suggesting that the amino acid substitutions are adaptive to their functions in piRNA pathway in the teleost fish species. Notably five piRNA pathway genes evolved faster in the swamp eel, a kind of protogynous hermaphrodite fish, than the other teleosts, indicating a differential evolution of piRNA pathway between the swamp eel and other gonochoristic fishes. In addition, genome-wide analysis showed higher diversity of transposons in the teleost fish species compared with mammals. Our results suggest that rapidly evolved piRNA pathway in the teleost fish is likely to be involved in the adaption to transposon diversity.

The mammalian sense of smell is governed by the largest gene family, which encodes the olfactory receptors (ORs). The gain and loss of OR genes is typically correlated with adaptations to various ecological niches. Modern humans have 853 OR genes but 55% of these have lost their function. Here we show evidence of additional OR loss of function in the Neanderthal and Denisovan hominin genomes using comparative genomic methodologies. Ten Neanderthal and 8 Denisovan ORs show evidence of loss of function that differ from the reference modern human OR genome. Some of these losses are also present in a subset of modern humans, while some are unique to each lineage. Morphological changes in the cranium of Neanderthals suggest different sensory arrangements to that of modern humans. We identify differences in functional olfactory receptor genes among modern humans, Neanderthals and Denisovans, suggesting varied loss of function across all three taxa and we highlight the utility of using genomic information to elucidate the sensory niches of extinct species.

Retinoic acid-related orphan receptor alpha gene (RORa) and the microRNA MIR137 have both recently been identified as novel candidate genes for neuropsychiatric disorders. RORa encodes a ligand-dependent orphan nuclear receptor that acts as a transcriptional regulator and miR-137 is a brain enriched small non-coding RNA that interacts with gene transcripts to control protein levels. Given the mounting evidence for RORa in autism spectrum disorders (ASD) and MIR137 in schizophrenia and ASD, we investigated if there was a functional biological relationship between these two genes. Herein, we demonstrate that miR-137 targets the 3'UTR of RORa in a site specific manner. We also provide further support for MIR137 as an autism candidate by showing that a large number of previously implicated autism genes are also putatively targeted by miR-137. This work supports the role of MIR137 as an ASD candidate and demonstrates a direct biological link between these previously unrelated autism candidate genes.

The vestibular system maintains the body's sense of balance and, therefore, was probably subject to strong selection during evolutionary transitions in locomotion. Among mammals, bats possess unique traits that place unusual demands on their vestibular systems. First, bats are capable of powered flight, which in birds is associated with enlarged semicircular canals. Second, many bats have enlarged cochleae associated with echolocation, and both cochleae and semicircular canals share a space within the petrosal bone. To determine how bat vestibular systems have evolved in the face of these pressures, we used micro-CT scans to compare canal morphology across species with contrasting flight and echolocation capabilities. We found no increase in canal radius in bats associated with the acquisition of powered flight, but canal radius did correlate with body mass in bat species from the suborder Yangochiroptera, and also in non-echolocating Old World fruit bats from the suborder Yinpterochiroptera. No such trend was seen in members of the Yinpterochiroptera that use laryngeal echolocation, although canal radius was associated with wing-tip roundedness in this group. We also found that the vestibular system scaled with cochlea size, although the relationship differed in species that use constant frequency echolocation. Across all bats, the shape of the anterior and lateral canals was associated with large cochlea size and small body size respectively, suggesting differential spatial constraints on each canal depending on its orientation within the skull. Thus in many echolocating bats, it seems that the combination of small body size and enlarged cochlea together act as a principal force on the vestibular system. The two main groups of echolocating bats displayed different canal morphologies, in terms of size and shape in relation to body mass and cochlear size, thus suggesting independent evolutionary pathways and offering tentative support for multiple acquisitions of echolocation.

Heterozygous mutations in CNTNAP2 have been identified in patients with a range of complex phenotypes including intellectual disability, autism and schizophrenia. However heterozygous CNTNAP2 mutations are also found in the normal population. Conversely, homozygous mutations are rare in patient populations and have not been found in any unaffected individuals.

The sense of smell, or olfaction, is fundamental in the life of animals. However, penguins (Aves: Sphenisciformes) possess relatively small olfactory bulbs compared with most other waterbirds such as Procellariiformes and Gaviiformes. To test whether penguins have a reduced reliance on olfaction, we analyzed the draft genome sequences of the two penguins, which diverged at the origin of the order Sphenisciformes; we also examined six closely related species with available genomes, and identified 29 one-to-one orthologous olfactory receptor genes (i.e. ORs) that are putatively functionally conserved and important across the eight birds. To survey the 29 one-to-one orthologous ORs in penguins and their relatives, we newly generated 34 sequences that are missing from the draft genomes. Through the analysis of totaling 378 OR sequences, we found that, of these functionally important ORs common to other waterbirds, penguins have a significantly greater percentage of OR pseudogenes than other waterbirds, suggesting a reduction of olfactory capability. The penguin-specific reduction of olfactory capability arose in the common ancestor of penguins between 23 and 60 Ma, which may have resulted from the aquatic specializations for underwater vision. Our study provides genetic evidence for a possible reduction of reliance on olfaction in penguins.

A fundamental condition for any work with free-ranging animals is correct species identification. However, in case of bats, information on local species assemblies is frequently limited especially in regions with high biodiversity such as the Neotropics. The bat genus Molossus is a typical example of this, with morphologically similar species often occurring in sympatry. We used a multi-method approach based on molecular, morphometric and acoustic information collected from 962 individuals of Molossus bondae, M. coibensis, and M. molossus captured in Panama. We distinguished M. bondae based on size and pelage coloration. We identified two robust species clusters composed of M. molossus and M. coibensis based on 18 microsatellite markers but also on a more stringently determined set of four markers. Phylogenetic reconstructions using the mitochondrial gene co1 (DNA barcode) were used to diagnose these microsatellite clusters as M. molossus and M. coibensis. To differentiate species, morphological information was only reliable when forearm length and body mass were combined in a linear discriminant function (95.9% correctly identified individuals). When looking in more detail at M. molossus and M. coibensis, only four out of 13 wing parameters were informative for species differentiation, with M. coibensis showing lower values for hand wing area and hand wing length and higher values for wing loading. Acoustic recordings after release required categorization of calls into types, yielding only two informative subsets: approach calls and two-toned search calls. Our data emphasizes the importance of combining morphological traits and independent genetic data to inform the best choice and combination of discriminatory information used in the field. Because parameters can vary geographically, the multi-method approach may need to be adjusted to local species assemblies and populations to be entirely informative.

We previously discovered that mutations of the human FOXP2 gene cause a monogenic communication disorder, primarily characterized by difficulties in learning to make coordinated sequences of articulatory gestures that underlie speech. Affected people have deficits in expressive and receptive linguistic processing and display structural and/or functional abnormalities in cortical and subcortical brain regions. FOXP2 provides a unique window into neural processes involved in speech and language. In particular, its role as a transcription factor gene offers powerful functional genomic routes for dissecting critical neurogenetic mechanisms. Here, we employ chromatin immunoprecipitation coupled with promoter microarrays (ChIP-chip) to successfully identify genomic sites that are directly bound by FOXP2 protein in native chromatin of human neuron-like cells. We focus on a subset of downstream targets identified by this approach, showing that altered FOXP2 levels yield significant changes in expression in our cell-based models and that FOXP2 binds in a specific manner to consensus sites within the relevant promoters. Moreover, we demonstrate significant quantitative differences in target expression in embryonic brains of mutant mice, mediated by specific in vivo Foxp2-chromatin interactions. This work represents the first identification and in vivo verification of neural targets regulated by FOXP2. Our data indicate that FOXP2 has dual functionality, acting to either repress or activate gene expression at occupied promoters. The identified targets suggest roles in modulating synaptic plasticity, neurodevelopment, neurotransmission, and axon guidance and represent novel entry points into in vivo pathways that may be disturbed in speech and language disorders.

The application of high-throughput sequencing (HTS) for metabarcoding of mixed samples offers new opportunities in conservation biology. Recently, the successful detection of prey DNA from the guts of leeches has raised the possibility that these, and other blood-feeding invertebrates, might serve as useful samplers of mammals. Yet little is known about whether sympatric leech species differ in their feeding preferences, and whether this has a bearing on their relative suitability for monitoring local mammalian diversity. To address these questions, we collected spatially matched samples of two congeneric leech species Haemadipsa picta and Haemadipsa sumatrana from lowland rainforest in Borneo. For each species, we pooled ~500 leeches into batches of 10 individuals, performed PCR to target a section of the mammalian 16S rRNA locus and undertook sequencing of amplicon libraries using an Illumina MiSeq. In total, we identified sequences from 14 mammalian genera, spanning nine families and five orders. We found greater numbers of detections, and higher diversity of OTUs, in H. picta compared with H. sumatrana, with rodents only present in the former leech species. However, comparison of samples from across the landscape revealed no significant difference in mammal community composition between the leech species. We therefore suggest that H. picta is the more suitable iDNA sampler in this degraded Bornean forest. We conclude that the choice of invertebrate sampler can influence the detectability of different mammal groups and that this should be accounted for when designing iDNA studies.

It is commonly, although not universally, accepted that most intra and interspecific genome sequence variations are more or less neutral, whereas a large fraction of organism-level phenotypic variations are adaptive. Gene expression levels are molecular phenotypes that bridge the gap between genotypes and corresponding organism-level phenotypes. Yet, it is unknown whether natural variations in gene expression levels are mostly neutral or adaptive. Here we address this fundamental question by genome-wide profiling and comparison of gene expression levels in nine yeast strains belonging to three closely related Saccharomyces species and originating from five different ecological environments. We find that the transcriptome-based clustering of the nine strains approximates the genome sequence-based phylogeny irrespective of their ecological environments. Remarkably, only ∼0.5% of genes exhibit similar expression levels among strains from a common ecological environment, no greater than that among strains with comparable phylogenetic relationships but different environments. These and other observations strongly suggest that most intra and interspecific variations in yeast gene expression levels result from the accumulation of random mutations rather than environmental adaptations. This finding has profound implications for understanding the driving force of gene expression evolution, genetic basis of phenotypic adaptation, and general role of stochasticity in evolution.

Divalent and trivalent salts exhibit a complex taste profile. They are perceived as being astringent/drying, sour, bitter, and metallic. We hypothesized that human bitter-taste receptors may mediate some taste attributes of these salts. Using a cell-based functional assay, we found that TAS2R7 responds to a broad range of divalent and trivalent salts, including zinc, calcium, magnesium, copper, manganese, and aluminum, but not to potassium, suggesting TAS2R7 may act as a metal cation receptor mediating bitterness of divalent and trivalent salts. Molecular modeling and mutagenesis analysis identified 2 residues, H943.37 and E2647.32, in TAS2R7 that appear to be responsible for the interaction of TAS2R7 with metallic ions. Taste receptors are found in both oral and extraoral tissues. The responsiveness of TAS2R7 to various mineral salts suggests it may act as a broad sensor, similar to the calcium-sensing receptor, for biologically relevant metal cations in both oral and extraoral tissues.

Understanding aging is a grand challenge in biology. Exceptionally long-lived animals have mechanisms that underpin extreme longevity. Telomeres are protective nucleotide repeats on chromosome tips that shorten with cell division, potentially limiting life span. Bats are the longest-lived mammals for their size, but it is unknown whether their telomeres shorten. Using >60 years of cumulative mark-recapture field data, we show that telomeres shorten with age in Rhinolophus ferrumequinum and Miniopterus schreibersii, but not in the bat genus with greatest longevity, Myotis. As in humans, telomerase is not expressed in Myotis myotis blood or fibroblasts. Selection tests on telomere maintenance genes show that ATM and SETX, which repair and prevent DNA damage, potentially mediate telomere dynamics in Myotis bats. Twenty-one telomere maintenance genes are differentially expressed in Myotis, of which 14 are enriched for DNA repair, and 5 for alternative telomere-lengthening mechanisms. We demonstrate how telomeres, telomerase, and DNA repair genes have contributed to the evolution of exceptional longevity in Myotis bats, advancing our understanding of healthy aging.

The evolution of taste perception is usually associated with the ecology and dietary changes of organisms. However, the association between feeding ecology and taste receptor evolution is unclear in some lineages of vertebrate animals. One example is the sweet taste receptor gene Tas1r2 Previous analysis of partial sequences has revealed that Tas1r2 has undergone equally strong purifying selection between insectivorous and frugivorous bats. To test whether the sweet taste function is also important in bats with contrasting diets, we examined the complete coding sequences of both sweet taste receptor genes (Tas1r2 and Tas1r3) in 34 representative bat species. Although these two genes are highly conserved between frugivorous and insectivorous bats at the sequence level, our behavioral experiments revealed that an insectivorous bat (Myotis ricketti) showed no preference for natural sugars, whereas the frugivorous species (Rousettus leschenaultii) showed strong preferences for sucrose and fructose. Furthermore, while both sweet taste receptor genes are expressed in the taste tissue of insectivorous and frugivorous bats, our cell-based assays revealed striking functional divergence: the sweet taste receptors of frugivorous bats are able to respond to natural sugars whereas those of insectivorous bats are not, which is consistent with the behavioral preference tests, suggesting that functional evolution of sweet taste receptors is closely related to diet. This comprehensive study suggests that using sequence conservation alone could be misleading in inferring protein and physiological function and highlights the power of combining behavioral experiments, expression analysis, and functional assays in molecular evolutionary studies.

Factors ranging from ecological opportunity to genome composition might explain why only some lineages form adaptive radiations. While being rare, particular systems can provide natural experiments within an identical ecological setting where species numbers and phenotypic divergence in two closely related lineages are notably different. We investigated one such natural experiment using two de novo assembled and 40 resequenced genomes and asked why two closely related Neotropical cichlid fish lineages, the Amphilophus citrinellus species complex (Midas cichlids; radiating) and Archocentrus centrarchus (Flyer cichlid; nonradiating), have resulted in such disparate evolutionary outcomes. Although both lineages inhabit many of the same Nicaraguan lakes, whole-genome inferred demography suggests that priority effects are not likely to be the cause of the dissimilarities. Also, genome-wide levels of selection, transposable element dynamics, gene family expansion, major chromosomal rearrangements and the number of genes under positive selection were not markedly different between the two lineages. To more finely investigate particular subsets of the genome that have undergone adaptive divergence in Midas cichlids, we also examined if there was evidence for 'molecular pre-adaptation' in regions identified by QTL mapping of repeatedly diverging adaptive traits. Although most of our analyses failed to pinpoint substantial genomic differences, we did identify functional categories containing many genes under positive selection that provide candidates for future studies on the propensity of Midas cichlids to radiate. Our results point to a disproportionate role of local, rather than genome-wide factors underlying the propensity for these cichlid fishes to adaptively radiate.

Morphological variation between individuals can increase niche segregation and decrease intraspecific competition when heterogeneous individuals explore their environment in different ways. Among bat species, wing shape correlates with flight maneuverability and habitat use, with species that possess broader wings typically foraging in more cluttered habitats. However, few studies have investigated the role of morphological variation in bats for niche partitioning at the individual level. To determine the relationship between wing shape and diet, we studied a population of the insectivorous bat species Pteronotus mesoamericanus in the dry forest of Costa Rica. Individual diet was resolved using DNA metabarcoding, and bat wing shape was assessed using geometric morphometric analysis. Inter-individual variation in wing shape showed a significant relationship with both dietary dissimilarity based on Bray-Curtis estimates, and nestedness derived from an ecological network. Individual bats with broader and more rounded wings were found to feed on a greater diversity of arthropods (less nested) in comparison to individuals with triangular and pointed wings (more nested). We conclude that individual variation in bat wing morphology can impact foraging efficiency leading to the observed overall patterns of diet specialization and differentiation within the population.

Adapted to different ecological niches, moth species belonging to the Hyles genus exhibit a spectacular diversity of larval color patterns. These species diverged ∼7.5 million years ago, making this rather young genus an interesting system to study a wide range of questions including the process of speciation, ecological adaptation, and adaptive radiation.

Observations that rates of introgression between taxa can vary across loci are increasingly common. Here, we test for differential locus-wise introgression in 2 parapatric subspecies of Pearson's horseshoe bat (Rhinolophus pearsoni chinensis and R. p. pearsoni). To efficiently identify putative speciation genes and/or beneficial genes in our current system, we used a candidate gene approach by including loci from X chromosome that are suggested to be more likely involved in reproductive isolation in other organisms and loci underlying hearing that have been suggested to spread across the hybrid zone in another congeneric species. Phylogenetic and coalescent analyses were performed at 2 X-linked, 4 hearing genes, as well as 2 other autosomal loci individually. Likelihood ratio tests could not reject the model of zero gene flow at 2 X-linked and 2 autosomal genes. In contrast, gene flow was supported at 3 of 4 hearing genes. While this introgression could be adaptive, we cannot rule out stochastic processes. Our results highlight the utility of the candidate gene approach in searching for speciation genes and/or beneficial genes across the species boundary in natural populations.