Parallel sites implicate functional convergence of the hearing gene prestin among echolocating mammals

Recent advancement of sonogenetics: A promising noninvasive cellular manipulation by ultrasound

Recent advancements in biomedical research have underscored the importance of noninvasive cellular manipulation techniques. Sonogenetics, a method that uses genetic engineering to produce ultrasound-sensitive proteins in target cells, is gaining prominence along with optogenetics, electrogenetics, and magnetogenetics. Upon stimulation with ultrasound, these proteins trigger a cascade of cellular activities and functions. Unlike traditional ultrasound modalities, sonogenetics offers enhanced spatial selectivity, improving precision and safety in disease treatment. This technology broadens the scope of non-surgical interventions across a wide range of clinical research and therapeutic applications, including neuromodulation, oncologic treatments, stem cell therapy, and beyond. Although current literature predominantly emphasizes ultrasonic neuromodulation, this review offers a comprehensive exploration of sonogenetics. We discuss ultrasound properties, the specific ultrasound-sensitive proteins employed in sonogenetics, and the technique's potential in managing conditions such as neurological disorders, cancer, and ophthalmic diseases, and in stem cell therapies. Our objective is to stimulate fresh perspectives for further research in this promising field.

WGCCRR: a web-based tool for genome-wide screening of convergent indels and substitutions of amino acids

Summary

Genome-wide analyses of proteincoding gene sequences are being employed to examine the genetic basis of adaptive evolution in many organismal groups. Previous studies have revealed that convergent/parallel adaptive evolution may be caused by convergent/parallel amino acid changes. Similarly, detailed analysis of lineage-specific amino acid changes has shown correlations with certain lineage-specific traits. However, experimental validation remains the ultimate measure of causality. With the increasing availability of genomic data, a streamlined tool for such analyses would facilitate and expedite the screening of genetic loci that hold potential for adaptive evolution, while alleviating the bioinformatic burden for experimental biologists. In this study, we present a user-friendly web-based tool called WGCCRR (Whole Genome Comparative Coding Region Read) designed to screen both convergent/parallel and lineage-specific amino acid changes on a genome-wide scale. Our tool allows users to replicate previous analyses with just a few clicks, and the exported results are straightforward to interpret. In addition, we have also included amino acid indels that are usually neglected in previous work. Our website provides an efficient platform for screening candidate loci for downstream experimental tests.

Availability and Implementation

The tool is available at: https://fishevo.xmu.edu.cn/.

Convergence in hearing-related genes between echolocating birds and mammals

Echolocation, the detection of objects by means of sound waves, has evolved independently in diverse animals. Echolocators include not only mammals such as toothed whales and yangochiropteran and rhinolophoid bats but also Rousettus fruit bats, as well as two bird lineages, oilbirds and swiftlets. In whales and yangochiropteran and rhinolophoid bats, positive selection and molecular convergence has been documented in key hearing-related genes, such as prestin (SLC26A5), but few studies have examined these loci in other echolocators. Here, we examine patterns of selection and convergence in echolocation-related genes in echolocating birds and Rousettus bats. Fewer of these loci were under selection in Rousettus or birds compared with classically recognized echolocators, and elevated convergence (compared to outgroups) was not evident across this gene set. In certain genes, however, we detected convergent substitutions with potential functional relevance, including convergence between Rousettus and classic echolocators in prestin at a site known to affect hair cell electromotility. We also detected convergence between Yangochiroptera, Rhinolophidea, and oilbirds in TMC1, an important mechanosensory transduction channel in vertebrate hair cells, and observed an amino acid change at the same site within the pore domain. Our results suggest that although most proteins implicated in echolocation in specialized mammals may not have been recruited in birds or Rousettus fruit bats, certain hearing-related loci may have undergone convergent functional changes. Investigating adaptations in diverse echolocators will deepen our understanding of this unusual sensory modality.

PolyQ length-based molecular encoding of vocalization frequency in FOXP2

The transcription factor FOXP2, a regulator of vocalization- and speech/language-related phenotypes, contains two long polyQ repeats (Q1 and Q2) displaying marked, still enigmatic length variation across mammals. We found that the Q1/Q2 length ratio quantitatively encodes vocalization frequency ranges, from the infrasonic to the ultrasonic, displaying striking convergent evolution patterns. Thus, species emitting ultrasonic vocalizations converge with bats in having a low ratio, whereas species vocalizing in the low-frequency/infrasonic range converge with elephants and whales, which have higher ratios. Similar, taxon-specific patterns were observed for the FOXP2-related protein FOXP1. At the molecular level, we observed that the FOXP2 polyQ tracts form coiled coils, assembling into condensates and fibrils, and drive liquid-liquid phase separation (LLPS). By integrating evolutionary and molecular analyses, we found that polyQ length variation related to vocalization frequency impacts FOXP2 structure, LLPS, and transcriptional activity, thus defining a novel form of polyQ length-based molecular encoding of vocalization frequency.

Repeated evolution of similar phenotypes: Integrating comparative methods with developmental pathways

Repeated phenotypes, often referred to as 'homoplasies' in cladistic analyses, may evolve through changes in developmental processes. Genetic bases of recurrent evolution gained attention and have been studied in the past years using approaches that combine modern analytical phylogenetic tools with the stunning assemblage of new information on developmental mechanisms. In this review, we evaluated the topic under an integrated perspective, revisiting the classical definitions of convergence and parallelism and detailing comparative methods used to evaluate evolution of repeated phenotypes, which include phylogenetic inference, estimates of evolutionary rates and reconstruction of ancestral states. We provide examples to illustrate how a given methodological approach can be used to identify evolutionary patterns and evaluate developmental mechanisms associated with the intermittent expression of a given trait along the phylogeny. Finally, we address why repeated trait loss challenges strict definitions of convergence and parallelism, discussing how changes in developmental pathways might explain the high frequency of repeated trait loss in specific lineages.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

The evolutionary tuning of hearing

After the transition to life on land, tympanic middle ears emerged separately in different groups of tetrapods, facilitating the efficient detection of airborne sounds and paving the way for high frequency sensitivity. The processes that brought about high-frequency hearing in mammals are tightly linked to the accumulation of coding sequence changes in inner ear genes; many of which were selected during evolution. These include proteins involved in hair bundle morphology, mechanotransduction and high endolymphatic potential, somatic electromotility for sound amplification, ribbon synapses for high-fidelity transmission of sound stimuli, and efferent synapses for the modulation of sound amplification. Here, we review the molecular evolutionary processes behind auditory functional innovation. Overall, the evidence to date supports the hypothesis that changes in inner ear proteins were central to the fine tuning of mammalian hearing.

What Has Genomics Taught An Evolutionary Biologist?

Genomics, an interdisciplinary field of biology on the structure, function, and evolution of genomes, has revolutionized many subdisciplines of life sciences, including my field of evolutionary biology, by supplying huge data, bringing high-throughput technologies, and offering a new approach to biology. In this review, I describe what I have learned from genomics and highlight the fundamental knowledge and mechanistic insights gained. I focus on three broad topics that are central to evolutionary biology and beyond-variation, interaction, and selection-and use primarily my own research and study subjects as examples. In the next decade or two, I expect that the most important contributions of genomics to evolutionary biology will be to provide genome sequences of nearly all known species on Earth, facilitate high-throughput phenotyping of natural variants and systematically constructed mutants for mapping genotype-phenotype-fitness landscapes, and assist the determination of causality in evolutionary processes using experimental evolution.

Detecting macroevolutionary genotype-phenotype associations using error-corrected rates of protein convergence

On macroevolutionary timescales, extensive mutations and phylogenetic uncertainty mask the signals of genotype-phenotype associations underlying convergent evolution. To overcome this problem, we extended the widely used framework of non-synonymous to synonymous substitution rate ratios and developed the novel metric ω_C, which measures the error-corrected convergence rate of protein evolution. While ω_C distinguishes natural selection from genetic noise and phylogenetic errors in simulation and real examples, its accuracy allows an exploratory genome-wide search of adaptive molecular convergence without phenotypic hypothesis or candidate genes. Using gene expression data, we explored over 20 million branch combinations in vertebrate genes and identified the joint convergence of expression patterns and protein sequences with amino acid substitutions in functionally important sites, providing hypotheses on undiscovered phenotypes. We further extended our method with a heuristic algorithm to detect highly repetitive convergence among computationally non-trivial higher-order phylogenetic combinations. Our approach allows bidirectional searches for genotype-phenotype associations, even in lineages that diverged for hundreds of millions of years.

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.

Sensory biology: Tree mice use echolocation

A new study demonstrates that soft-furred tree mice orientate by using echolocation, emitting ultrasonic broadband chirps. Remarkable convergent evolution with distantly related bats and dolphins in ear bone morphology and sensory genes is evident.

Sonogenetic-Based Neuromodulation for the Amelioration of Parkinson's Disease

Sonogenetics is a promising strategy allowing the noninvasive and selective activation of targeted neurons in deep brain regions; nevertheless, its therapeutic outcome for neurodegeneration diseases that need long-term treatment remains to be verified. We previously enhanced the ultrasound (US) sensitivity of targeted cells by genetic modification with an engineered auditory-sensing protein, mPrestin (N7T, N308S). In this study, we expressed mPrestin in the dopaminergic neurons of the substantia nigra in Parkinson's disease (PD) mice and used 0.5 MHz US for repeated and localized brain stimulation. The mPrestin expression in dopaminergic neurons persisted for at least 56 days after a single shot of adeno-associated virus, suggesting that the period of expression was long enough for US treatment in mice. Compared to untreated mice, US stimulation ameliorated the dopaminergic neurodegeneration 10-fold and mitigated the PD symptoms of the mice 4-fold, suggesting that this sonogenetic strategy has the clinical potential to treat neurodegenerative diseases.

Cochlear hair cells of echolocating bats are immune to intense noise

Exposure to intense noise can damage cochlear hair cells, leading to hearing loss in mammals. To avoid this constraint, most mammals have evolved in relatively quiet environments. Echolocating bats, however, are naturally exposed to continuous intense sounds from their own and neighboring sonar emissions for maintaining sonar directionality and range. Here, we propose the presence of intense noise resistance in cochlear hair cells of echolocating bats against noise-induced hearing loss (NIHL). To test this hypothesis, we performed noise exposure experiments for laboratory mice, one nonecholocating bat species, and five echolocating bat species. Contrary to nonecholocating fruit bats and mice, the hearing and the cochlear hair cells of echolocating bats remained unimpaired after continuous intense noise exposure. The comparative analyses of cochleae transcriptomic data showed that several genes protecting cochlear hair cells from intense sounds were overexpressed in echolocating bats. Particularly, the experimental examinations revealed that ISL1 overexpression significantly improved the survival of cochlear hair cells. Our findings support the existence of protective effects in cochlear hair cells of echolocating bats against intense noises, which provides new insight into understanding the relationship between cochlear hair cells and intense noises, and preventing or ameliorating NIHL in mammals.

Echolocation in soft-furred tree mice

Echolocation is the use of reflected sound to sense features of the environment. Here, we show that soft-furred tree mice (Typhlomys) echolocate based on multiple independent lines of evidence. Behavioral experiments show that these mice can locate and avoid obstacles in darkness using hearing and ultrasonic pulses. The proximal portion of their stylohyal bone fuses with the tympanic bone, a form previously only seen in laryngeally echolocating bats. Further, we found convergence of hearing-related genes across the genome and of the echolocation-related gene prestin between soft-furred tree mice and echolocating mammals. Together, our findings suggest that soft-furred tree mice are capable of echolocation, and thus are a new lineage of echolocating mammals.

Functional Parameters of Prestin Are Not Correlated With the Best Hearing Frequency

Among the vertebrate lineages with different hearing frequency ranges, humans lie between the low-frequency hearing (1 kHz) of fish and amphibians and the high-frequency hearing (100 kHz) of bats and dolphins. Little is known about the mechanism underlying such a striking difference in the limits of hearing frequency. Prestin, responsible for cochlear amplification and frequency selectivity in mammals, seems to be the only candidate to date. Mammalian prestin is densely expressed in the lateral plasma membrane of the outer hair cells (OHCs) and functions as a voltage-dependent motor protein. To explore the molecular basis for the contribution of prestin in hearing frequency detection, we collected audiogram data from humans, dogs, gerbils, bats, and dolphins because their average hearing frequency rises in steps. We generated stable cell lines transfected with human, dog, gerbil, bat, and dolphin prestins (hPres, dPres, gPres, bPres, and nPres, respectively). The non-linear capacitance (NLC) of different prestins was measured using a whole-cell patch clamp. We found that the Q_max/C_lin of bPres and nPres was significantly higher than that of humans. The V_1/2 of hPres was more hyperpolarized than that of nPres. The z values of hPres and bPres were higher than that of nPres. We further analyzed the relationship between the high-frequency hearing limit (F_max) and the functional parameters (V_1/2, z, and Q_max/C_lin) of NLC among five prestins. Interestingly, no significant correlation was found between the functional parameters and F_max. Additionally, a comparative study showed that the amino acid sequences and tertiary structures of five prestins were quite similar. There might be a common fundamental mechanism driving the function of prestins. These findings call for a reconsideration of the leading role of prestin in hearing frequency perception. Other intriguing kinetics underlying the hearing frequency response of auditory organs might exist.

Molecular convergent and parallel evolution among four high-elevation anuran species from the Tibetan region

BACKGROUND

To date, evidence for the relative prevalence or rarity of molecular convergent and parallel evolution is conflicting, and understanding of how these processes contribute to adaptation is limited. We compared four high-elevation anuran species (Bufo tibetanus, Nanorana parkeri, Rana kukunoris and Scutiger boulengeri) from the Tibetan region, and examined convergent and parallel amino acid substitutions between them and how they may have contributed to high-elevation adaptation.

RESULTS

Genomic data of the four high-elevation species and eight of their low-elevation close relatives were gathered. A total of 1098 orthologs shared by all species were identified. We first conducted pairwise comparisons using Zhang and Kumar's test. Then, the R_conv index was calculated and convergence/divergence correlation plotting was conducted. Furthermore, genes under positive selection and with elevated evolutionary rate were examined. We detected a large number of amino acid sites with convergent or parallel substitutions. Several pairs of high-elevation species, in particular, R. kukunoris vs N. parkeri and B. tibetanus vs S. boulengeri, had excessive amounts of convergent substitutions compared to neutral expectation. Nevertheless, these sites were mostly concentrated in a small number of genes (3-32), and no genome-wide convergence was detected. Furthermore, the majority of these convergent genes were neither under detectable positive selection nor had elevated evolutionary rates, although functional prediction analysis suggested some of the convergent genes could potentially contribute to high-elevation adaptation.

CONCLUSIONS

There is a substantial amount of convergent evolution at the amino-acid level among high-elevation amphibians, although these sites are concentrated in a few genes, not widespread across the genomes. This may attribute to the fact that all the target species are from the same environment. The relative prevalence of convergent substitutions among high-elevation amphibians provides an excellent opportunity for further study of molecular convergent evolution.

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.

Functional evolution of vertebrate sensory receptors

Sensory receptors enable animals to perceive their external world, and functional properties of receptors evolve to detect the specific cues relevant for an organism's survival. Changes in sensory receptor function or tuning can directly impact an organism's behavior. Functional tests of receptors from multiple species and the generation of chimeric receptors between orthologs with different properties allow for the dissection of the molecular basis of receptor function and identification of the key residues that impart functional changes in different species. Knowledge of these functionally important sites facilitates investigation into questions regarding the role of epistasis and the extent of convergence, as well as the timing of sensory shifts relative to other phenotypic changes. However, as receptors can also play roles in non-sensory tissues, and receptor responses can be modulated by numerous other factors including varying expression levels, alternative splicing, and morphological features of the sensory cell, behavioral validation can be instrumental in confirming that responses observed in heterologous systems play a sensory role. Expression profiling of sensory cells and comparative genomics approaches can shed light on cell-type specific modifications and identify other proteins that may affect receptor function and can provide insight into the correlated evolution of complex suites of traits. Here we review the evolutionary history and diversity of functional responses of the major classes of sensory receptors in vertebrates, including opsins, chemosensory receptors, and ion channels involved in temperature-sensing, mechanosensation and electroreception.

Echolocating Whales and Bats Express the Motor Protein Prestin in the Inner Ear: A Potential Marker for Hearing Loss

Prestin is an integral membrane motor protein located in outer hair cells of the mammalian cochlea. It is responsible for electromotility and required for cochlear amplification. Although prestin works in a cycle-by-cycle mode up to frequencies of at least 79 kHz, it is not known whether or not prestin is required for the extreme high frequencies used by echolocating species. Cetaceans are known to possess a prestin coding gene. However, the expression and distribution pattern of the protein in the cetacean cochlea has not been determined, and the contribution of prestin to echolocation has not yet been resolved. Here we report the expression of the protein prestin in five species of echolocating whales and two species of echolocating bats. Positive labeling in the basolateral membrane of outer hair cells, using three anti-prestin antibodies, was found all along the cochlear spiral in echolocating species. These findings provide morphological evidence that prestin can have a role in cochlear amplification in the basolateral membrane up to 120–180 kHz. In addition, labeling of the cochlea with a combination of anti-prestin, anti-neurofilament, anti-myosin VI and/or phalloidin and DAPI will be useful for detecting potential recent cases of noise-induced hearing loss in stranded cetaceans. This study improves our understanding of the mechanisms involved in sound transduction in echolocating mammals, as well as describing an optimized methodology for detecting cases of hearing loss in stranded marine mammals.

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.

Complex nonlinear capacitance in outer hair cell macro-patches: effects of membrane tension

Outer hair cell (OHC) nonlinear capacitance (NLC) represents voltage sensor charge movements of prestin (SLC26a5), the protein responsible for OHC electromotility. Previous measures of NLC frequency response have employed methods which did not assess the influence of dielectric loss (sensor charge movements out of phase with voltage) that may occur, and such loss conceivably may influence prestin's frequency dependent activity. Here we evaluate prestin's complex capacitance out to 30 kHz and find that prestin's frequency response determined using this approach coincides with all previous estimates. We also show that membrane tension has no effect on prestin's frequency response, despite substantial shifts in its voltage operating range, indicating that prestin transition rate alterations do not account for the shifts. The magnitude roll-off of prestin activity across frequency surpasses the reductions of NLC caused by salicylate treatments that are known to abolish cochlear amplification. Such roll-off likely limits the effectiveness of prestin in contributing to cochlear amplification at the very high acoustic frequencies processed by some mammals.

Targeted delivery of engineered auditory sensing protein for ultrasound neuromodulation in the brain

Sonogenetics is a promising approach for in vivo neuromodulation using ultrasound (US) to non-invasively stimulate cells in deep tissue. However, sonogenetics requires accurate transduction of US-responsive proteins into target cells. Here, we introduce a non-invasive and non-viral approach for intracerebral gene delivery. This approach utilizes temporary ultrasonic disruption of the blood-brain barrier (BBB) to transfect neurons at specific sites in the brain via DNA that encodes engineered US-responsive protein (murine Prestin (N7T, N308S))-loaded microbubbles (pPrestin-MBs). Prestin is a transmembrane protein that exists in the mammalian auditory system and functions as an electromechanical transducer. We further improved the US sensitivity of Prestin by introducing specific amino acid substitutions that frequently occur in sonar species into the mouse Prestin protein. We demonstrated this concept in mice using US with pPrestin-MBs to non-invasively modify and activate neurons within the brain for spatiotemporal neuromodulation.

Method: MBs composed of cationic phospholipid and C₃F₈ loaded with mouse Prestin plasmid (pPrestin) via electrostatic interactions. The mean concentration and size of the pPrestin-MBs were (16.0 ± 0.2) × 10⁹ MBs/mL and 1.1 ± 0.2 μm, respectively. SH-SY5Y neuron-like cells and C57BL mice were used in this study. We evaluated the gene transfection efficiency and BBB-opening region resulting from pPrestin-MBs with 1-MHz US (pressure = 0.1-0.5 MPa, cycle = 50-10000, pulse repetition frequency (PRF): 0.5-5 Hz, sonication time = 60 s) using green fluorescence protein (Venus) and Evans blue staining.

Results: The maximum pPrestin expression with the highest cell viability occurred at a pressure of 0.5 MPa, cycle number of 5000, and PRF of 1 Hz. The cellular transfection rate with pPrestin-MBs and US was 20.2 ± 2.5%, which was 1.5-fold higher than that of commercial transfection agents (LT-1). In vivo data suggested that the most profound expression of pPrestin occurred at 2 days after performing pPrestin-MBs with US (0.5 MPa, 240 s sonication time). In addition, no server erythrocyte extravasations and apoptosis cells were observed at US-sonicated region. We further found that with 0.5-MHz US stimulation, cells with Prestin expression were 6-fold more likely to exhibit c-Fos staining than cells without Prestin expression.

Conclusion: Successful activation of Prestin-expressing neurons suggests that this technology provides non-invasive and spatially precise selective modulation of one or multiple specific brain regions.

Sonogenetic Modulation of Cellular Activities Using an Engineered Auditory-Sensing Protein

Biomolecules that respond to different external stimuli enable the remote control of genetically modified cells. We report herein a sonogenetic approach that can manipulate target cell activities by focused ultrasound stimulation. This system requires an ultrasound-responsive protein derived from an engineered auditory-sensing protein prestin. Heterologous expression of mouse prestin containing two parallel amino acid substitutions, N7T and N308S, that frequently exist in prestins from echolocating species endowed transfected mammalian cells with the ability to sense ultrasound. An ultrasound pulse of low frequency and low pressure efficiently evoked cellular calcium responses after transfecting with prestin(N7T, N308S). Moreover, pulsed ultrasound can also noninvasively stimulate target neurons expressing prestin(N7T, N308S) in deep regions of mouse brains. Our study delineates how an engineered auditory-sensing protein can cause mammalian cells to sense ultrasound stimulation. Moreover, our sonogenetic tools will serve as new strategies for noninvasive therapy in deep tissues.

Convergent evolution in toothed whale cochleae

BACKGROUND

Odontocetes (toothed whales) are the most species-rich marine mammal lineage. The catalyst for their evolutionary success is echolocation - a form of biological sonar that uses high-frequency sound, produced in the forehead and ultimately detected by the cochlea. The ubiquity of echolocation in odontocetes across a wide range of physical and acoustic environments suggests that convergent evolution of cochlear shape is likely to have occurred. To test this, we used SURFACE; a method that fits Ornstein-Uhlenbeck (OU) models with stepwise AIC (Akaike Information Criterion) to identify convergent regimes on the odontocete phylogeny, and then tested whether convergence in these regimes was significantly greater than expected by chance.

RESULTS

We identified three convergent regimes: (1) True's (Mesoplodon mirus) and Cuvier's (Ziphius cavirostris) beaked whales; (2) sperm whales (Physeter macrocephalus) and all other beaked whales sampled; and (3) pygmy (Kogia breviceps) and dwarf (Kogia sima) sperm whales and Dall's porpoise (Phocoenoides dalli). Interestingly the 'river dolphins', a group notorious for their convergent morphologies and riverine ecologies, do not have convergent cochlear shapes. The first two regimes were significantly convergent, with habitat type and dive type significantly correlated with membership of the sperm whale + beaked whale regime.

CONCLUSIONS

The extreme acoustic environment of the deep ocean likely constrains cochlear shape, causing the cochlear morphology of sperm and beaked whales to converge. This study adds support for cochlear morphology being used to predict the ecology of extinct cetaceans.

A functional enrichment test for molecular convergent evolution finds a clear protein-coding signal in echolocating bats and whales

Distantly related species entering similar biological niches often adapt by evolving similar morphological and physiological characters. How much genomic molecular convergence (particularly of highly constrained coding sequence) contributes to convergent phenotypic evolution, such as echolocation in bats and whales, is a long-standing fundamental question. Like others, we find that convergent amino acid substitutions are not more abundant in echolocating mammals compared to their outgroups. However, we also ask a more informative question about the genomic distribution of convergent substitutions by devising a test to determine which, if any, of more than 4,000 tissue-affecting gene sets is most statistically enriched with convergent substitutions. We find that the gene set most overrepresented (q-value = 2.2e-3) with convergent substitutions in echolocators, affecting 18 genes, regulates development of the cochlear ganglion, a structure with empirically supported relevance to echolocation. Conversely, when comparing to nonecholocating outgroups, no significant gene set enrichment exists. For aquatic and high-altitude mammals, our analysis highlights 15 and 16 genes from the gene sets most affected by molecular convergence which regulate skin and lung physiology, respectively. Importantly, our test requires that the most convergence-enriched set cannot also be enriched for divergent substitutions, such as in the pattern produced by inactivated vision genes in subterranean mammals. Showing a clear role for adaptive protein-coding molecular convergence, we discover nearly 2,600 convergent positions, highlight 77 of them in 3 organs, and provide code to investigate other clades across the tree of life.

Comparative Transcriptome Analyses Reveal the Role of Conserved Function in Electric Organ Convergence Across Electric Fishes

The independent origins of multiple electric organs (EOs) of fish are fascinating examples of convergent evolution. However, comparative transcriptomics of different electric fish lineages are scarce. In this study, we found that the gene expression of EOs and skeletal muscles from three lineages (Mormyroidea, Siluriformes, and Gymnotiformes) tended to cluster together based on the species of origin, irrespective of the organ from which they are derived. A pairwise comparison of differentially expressed genes (DEGs) revealed that no less than half of shared DEGs exhibited parallel expression differentiation, indicating conserved directionality of differential expression either in or between lineages, but only a few shared DEGs were identified across all focal species. Nevertheless, the functional enrichment analysis of DEGs indicated that there were more parallel gene expression changes at the level of pathways and biological functions. Therefore, we may conclude that there is no parallel evolution of the entire transcriptomes of EOs among different lineages. Further, our results support the hypothesis that it is not different genes but conserved biological functions that play a crucial role in the convergence of complex phenotypes. This study provides insight into the genetic basis underlying the EO convergent evolution; however, more studies in different cases will be needed to demonstrate whether this pattern can be extended to other cases to derive a general rule for convergent evolution.

Functional Shifts in Bat Dim-Light Visual Pigment Are Associated with Differing Echolocation Abilities and Reveal Molecular Adaptation to Photic-Limited Environments

Bats are excellent models for studying the molecular basis of sensory adaptation. In Chiroptera, a sensory trade-off has been proposed between the visual and auditory systems, though the extent of this association has yet to be fully examined. To investigate whether variation in visual performance is associated with echolocation, we experimentally assayed the dim-light visual pigment rhodopsin from bat species with differing echolocation abilities. While spectral tuning properties were similar among bats, we found that the rate of decay of their light-activated state was significantly slower in a nonecholocating bat relative to species that use distinct echolocation strategies, consistent with a sensory trade-off hypothesis. We also found that these rates of decay were remarkably slower compared with those of other mammals, likely indicating an adaptation to dim light. To examine whether functional changes in rhodopsin are associated with shifts in selection intensity upon bat Rh1 sequences, we implemented selection analyses using codon-based likelihood clade models. While no shifts in selection were identified in response to diverse echolocation abilities of bats, we detected a significant increase in the intensity of evolutionary constraint accompanying the diversification of Chiroptera. Taken together, this suggests that substitutions that modulate the stability of the light-activated rhodopsin state were likely maintained through intensified constraint after bats diversified, being finely tuned in response to novel sensory specializations. Our study demonstrates the power of combining experimental and computational approaches for investigating functional mechanisms underlying the evolution of complex sensory adaptations.

The role of ecological factors in shaping bat cone opsin evolution

Bats represent one of the largest and most striking nocturnal mammalian radiations, exhibiting many visual system specializations for performance in light-limited environments. Despite representing the greatest ecological diversity and species richness in Chiroptera, Neotropical lineages have been undersampled in molecular studies, limiting the potential for identifying signatures of selection on visual genes associated with differences in bat ecology. Here, we investigated how diverse ecological pressures mediate long-term shifts in selection upon long-wavelength (Lws) and short-wavelength (Sws1) opsins, photosensitive cone pigments that form the basis of colour vision in most mammals, including bats. We used codon-based likelihood clade models to test whether ecological variables associated with reliance on visual information (e.g. echolocation ability and diet) or exposure to varying light environments (e.g. roosting behaviour and foraging habitat) mediated shifts in evolutionary rates in bat cone opsin genes. Using additional cone opsin sequences from newly sequenced eye transcriptomes of six Neotropical bat species, we found significant evidence for different ecological pressures influencing the evolution of the cone opsins. While Lws is evolving under significantly lower constraint in highly specialized high-duty cycle echolocating lineages, which have enhanced sonar ability to detect and track targets, variation in Sws1 constraint was significantly associated with foraging habitat, exhibiting elevated rates of evolution in species that forage among vegetation. This suggests that increased reliance on echolocation as well as the spectral environment experienced by foraging bats may differentially influence the evolution of different cone opsins. Our study demonstrates that different ecological variables may underlie contrasting evolutionary patterns in bat visual opsins, and highlights the suitability of clade models for testing ecological hypotheses of visual evolution.

Genomic Insights into the Adaptive Convergent Evolution

Adaptive convergent evolution, which refers to the same or similar phenotypes produced by species from independent lineages under similar selective pressures, has been widely examined for a long time. Accumulating studies on the adaptive convergent evolution have been reported from many different perspectives (cellular, anatomical, morphological, physiological, biochemical, and behavioral). Recent advances in the genomic technologies have demonstrated that adaptive convergence can arise from specific genetic mechanisms in different hierarchies, ranging from the same nucleotide or amino acid substitutions to the biological functions or pathways. Among these genetic mechanisms, the same amino acid changes in protein-coding genes play an important role in adaptive phenotypic convergence. Methods for detecting adaptive convergence at the protein sequence level have been constantly debated and developed. Here, we review recent progress on using genomic approaches to evaluate the genetic mechanisms of adaptive convergent evolution, summarize the research methods for identifying adaptive amino acid convergence, and discuss the future perspectives for researching adaptive convergent evolu-tion.

Genomic and functional evidence reveals molecular insights into the origin of echolocation in whales

Echolocation allows toothed whales to adapt to underwater habitats where vision is ineffective. Because echolocation requires the ability to detect exceptional high-frequency sounds, fossils related to the auditory system can help to pinpoint the origin of echolocation in whales. However, because of conflicting interpretations of archaeocete fossils, when and how whales evolved the high-frequency hearing correlated with echolocation remain unclear. We address these questions at the molecular level by systematically investigating the convergent evolution of 7206 orthologs across 16 mammals and find that convergent genes between the last common ancestor of all whales (LCAW) and echolocating bats are not significantly enriched in functional categories related to hearing, and that convergence in hearing-related proteins between them is not stronger than that between nonecholocating mammalian lineages and echolocating bats. However, these results contrast with those of parallel analyses between the LCA of toothed whales (LCATW) and echolocating bats. Furthermore, we reconstruct the ancestral genes for the hearing protein prestin for the LCAW and LCATW; we show that the LCAW prestin exhibits the same function as that of nonecholocating mammals, but the LCATW prestin shows functional convergence with that of extant echolocating mammals. Mutagenesis shows that functional convergence of prestin is driven by convergent changes in the prestins S392A and L497M in the LCATW and echolocating bats. Our results provide genomic and functional evidence supporting the origin of high-frequency hearing in the LCAW, not the LCATW, and reveal molecular insights into the origin and evolutionary trajectories of echolocation in whales.

Limited Evidence for Parallel Molecular Adaptations Associated with the Subterranean Niche in Mammals: A Comparative Study of Three Superorders

Among mammals, several lineages have independently adapted to a subterranean niche and possess similar phenotypic traits for burrowing (e.g., cylindrical bodies, short limbs, and absent pinnae). Previous research on mole-rats has revealed molecular adaptations for coping with reduced oxygen, elevated carbon dioxide, and the absence of light. In contrast, almost nothing is known regarding molecular adaptations in other subterranean lineages (e.g., true moles and golden moles). Therefore, the extent to which the recurrent phenotypic adaptations of divergent subterranean taxa have arisen via parallel routes of molecular evolution remains untested. To address these issues, we analyzed ∼8,000 loci in 15 representative subterranean taxa of four independent transitions to an underground niche for signatures of positive selection and convergent amino acid substitutions. Complementary analyses were performed in nonsubterranean "control" taxa to assess the biological significance of results. We found comparable numbers of positively selected genes in each of the four subterranean groups; however, correspondence in terms of gene identity between gene sets was low. Furthermore, we did not detect evidence of more convergent amino acids among subterranean species pairs compared with levels found between nonsubterranean controls. Comparisons with nonsubterranean taxa also revealed loci either under positive selection or with convergent substitutions, with similar functional enrichment (e.g., cell adhesion, immune response, and coagulation). Given the limited indication that positive selection and convergence occurred in the same loci, we conclude that selection may have acted on different loci across subterranean mammal lineages to produce similar phenotypes.

Molecular parallelism in fast-twitch muscle proteins in echolocating mammals

Detecting associations between genomic changes and phenotypic differences is fundamental to understanding how phenotypes evolved. By systematically screening for parallel amino acid substitutions, we detected known as well as novel cases (Strc, Tecta, and Cabp2) of parallelism between echolocating bats and toothed whales in proteins that could contribute to high-frequency hearing adaptations. Our screen also showed that echolocating mammals exhibit an unusually high number of parallel substitutions in fast-twitch muscle fiber proteins. Both echolocating bats and toothed whales produce an extremely rapid call rate when homing in on their prey, which was shown in bats to be powered by specialized superfast muscles. We show that these genes with parallel substitutions (Casq1, Atp2a1, Myh2, and Myl1) are expressed in the superfast sound-producing muscle of bats. Furthermore, we found that the calcium storage protein calsequestrin 1 of the little brown bat and the bottlenose dolphin functionally converged in its ability to form calcium-sequestering polymers at lower calcium concentrations, which may contribute to rapid calcium transients required for superfast muscle physiology. The proteins that our genomic screen detected could be involved in the convergent evolution of vocalization in echolocating mammals by potentially contributing to both rapid Ca²⁺ transients and increased shortening velocities in superfast muscles.

Epistatic interactions influence terrestrial-marine functional shifts in cetacean rhodopsin

Like many aquatic vertebrates, whales have blue-shifting spectral tuning substitutions in the dim-light visual pigment, rhodopsin, that are thought to increase photosensitivity in underwater environments. We have discovered that known spectral tuning substitutions also have surprising epistatic effects on another function of rhodopsin, the kinetic rates associated with light-activated intermediates. By using absorbance spectroscopy and fluorescence-based retinal release assays on heterologously expressed rhodopsin, we assessed both spectral and kinetic differences between cetaceans (killer whale) and terrestrial outgroups (hippo, bovine). Mutation experiments revealed that killer whale rhodopsin is unusually resilient to pleiotropic effects on retinal release from key blue-shifting substitutions (D83N and A292S), largely due to a surprisingly specific epistatic interaction between D83N and the background residue, S299. Ancestral sequence reconstruction indicated that S299 is an ancestral residue that predates the evolution of blue-shifting substitutions at the origins of Cetacea. Based on these results, we hypothesize that intramolecular epistasis helped to conserve rhodopsin's kinetic properties while enabling blue-shifting spectral tuning substitutions as cetaceans adapted to aquatic environments. Trade-offs between different aspects of molecular function are rarely considered in protein evolution, but in cetacean and other vertebrate rhodopsins, may underlie multiple evolutionary scenarios for the selection of specific amino acid substitutions.

Adaptive sequence convergence of the tumor suppressor ADAMTS9 between small-bodied mammals displaying exceptional longevity

Maximum lifespan varies by two orders of magnitude across mammals. How such divergent lifespans have evolved remains an open question, with ramifications that may potentially lead to therapies for age-related diseases in humans. Several species of microbats as well as the naked mole-rat live much longer than expected given their small sizes, show reduced susceptibility to neoplasia, and largely remain healthy and reproductively capable throughout the majority of their extended lifespans. The convergent evolution of extreme longevity in these two groups allows for the opportunity to identify potentially important aging related genes that have undergone adaptive sequence convergence in these long-lived, yet small-bodied species. Here, we have tested 4,628 genes for evidence of convergence between the microbats and naked mole-rat. We find a strong signal of adaptive sequence convergence in the gene A disintegrin-like and metalloprotease with thrombospondin type 1 motifs 9 (ADAMTS9). We also provide evidence that the shared substitutions were driven by selection. Intriguingly, ADAMTS9 is a known inhibitor of the mTor pathway and has been implicated in several aging related processes.

Mutational Biases Influence Parallel Adaptation

While mutational biases strongly influence neutral molecular evolution, the role of mutational biases in shaping the course of adaptation is less clear. Here we consider the frequency of transitions relative to transversions among adaptive substitutions. Because mutation rates for transitions are higher than those for transversions, if mutational biases influence the dynamics of adaptation, then transitions should be overrepresented among documented adaptive substitutions. To test this hypothesis, we assembled two sets of data on putatively adaptive amino acid replacements that have occurred in parallel during evolution, either in nature or in the laboratory. We find that the frequency of transitions in these data sets is much higher than would be predicted under a null model where mutation has no effect. Our results are qualitatively similar even if we restrict ourself to changes that have occurred, not merely twice, but three or more times. These results suggest that the course of adaptation is biased by mutation.

Molecular Convergent Evolution of the MYBPC2 Gene Among Three High-Elevation Amphibian Species

We report a strong pattern of molecular-level convergent/parallel evolution of the MYBPC2 gene. Three high-elevation amphibian species, Bufo gargarizans minshanicus, Nanorana pleskei, Rana kukunoris, revealed remarkable numbers of convergent and parallel amino acid substitutions. On the MYBPC2 gene tree of eleven anurans, the three distantly related species formed a strongly supported clade that was away from their respective relatives. Furthermore, we generated both model-based and empirical data-based null distributions for neutral convergent evolution. All three pairwise comparisons among the three species showed significantly more convergent and parallel substitutions than the null distributions. This study adds to the very small roster of clear cases of non-neutral molecular convergent evolution (e.g. prestin, rhodopsin). Molecular convergent evolution has significant implications in biology and detailed case studies will likely provide more insight into its genetic mechanisms.

Morphological variation among the inner ears of extinct and extant baleen whales (Cetacea: Mysticeti)

Living mysticetes (baleen whales) and odontocetes (toothed whales) differ significantly in auditory function in that toothed whales are sensitive to high-frequency and ultrasonic sound vibrations and mysticetes to low-frequency and infrasonic noises. Our knowledge of the evolution and phylogeny of cetaceans, and mysticetes in particular, is at a point at which we can explore morphological and physiological changes within the baleen whale inner ear. Traditional comparative anatomy and landmark-based 3D-geometric morphometric analyses were performed to investigate the anatomical diversity of the inner ears of extinct and extant mysticetes in comparison with other cetaceans. Principal component analyses (PCAs) show that the cochlear morphospace of odontocetes is tangential to that of mysticetes, but odontocetes are completely separated from mysticetes when semicircular canal landmarks are combined with the cochlear data. The cochlea of the archaeocete Zygorhiza kochii and early diverging extinct mysticetes plot within the morphospace of crown mysticetes, suggesting that mysticetes possess ancestral cochlear morphology and physiology. The PCA results indicate variation among mysticete species, although no major patterns are recovered to suggest separate hearing or locomotor regimes. Phylogenetic signal was detected for several clades, including crown Cetacea and crown Mysticeti, with the most clades expressing phylogenetic signal in the semicircular canal dataset. Brownian motion could not be excluded as an explanation for the signal, except for analyses combining cochlea and semicircular canal datasets for Balaenopteridae. J. Morphol. 277:1599-1615, 2016. © 2016 Wiley Periodicals, Inc.

Controlling for Phylogenetic Relatedness and Evolutionary Rates Improves the Discovery of Associations Between Species' Phenotypic and Genomic Differences

The growing number of sequenced genomes allows us now to address a key question in genetics and evolutionary biology: which genomic changes underlie particular phenotypic changes between species? Previously, we developed a computational framework called Forward Genomics that associates phenotypic to genomic differences by focusing on phenotypes that are independently lost in different lineages. However, our previous implementation had three main limitations. Here, we present two new Forward Genomics methods that overcome these limitations by (1) directly controlling for phylogenetic relatedness, (2) controlling for differences in evolutionary rates, and (3) computing a statistical significance. We demonstrate on large-scale simulated data and on real data that both new methods substantially improve the sensitivity to detect associations between phenotypic and genomic differences. We applied these new methods to detect genomic differences involved in the loss of vision in the blind mole rat and the cape golden mole, two independent subterranean mammals. Forward Genomics identified several genes that are enriched in functions related to eye development and the perception of light, as well as genes involved in the circadian rhythm. These new Forward Genomics methods represent a significant advance in our ability to discover the genomic basis underlying phenotypic differences between species. Source code: https://github.com/hillerlab/ForwardGenomics/.

Selection on different genes with equivalent functions: the convergence story told by Hox genes along the evolution of aquatic mammalian lineages

Background

Convergent evolution has been a challenging topic for decades, being cetaceans, pinnipeds and sirenians textbook examples of three independent origins of equivalent phenotypes. These mammalian lineages acquired similar anatomical features correlated to an aquatic life, and remarkably differ from their terrestrial counterparts. Whether their molecular evolutionary history also involved similar genetic mechanisms underlying such morphological convergence nevertheless remained unknown. To test for the existence of convergent molecular signatures, we studied the molecular evolution of Hox genes in these three aquatic mammalian lineages, comparing their patterns to terrestrial mammals. Hox genes are transcription factors that play a pivotal role in specifying embryonic regional identity of nearly any bilateral animal, and are recognized major agents for diversification of body plans.

Results

We detected few signatures of positive selection on Hox genes across the three aquatic mammalian lineages and verified that purifying selection prevails in these sequences, as expected for pleiotropic genes. Genes found as being positively selected differ across the aquatic mammalian lineages, but we identified a substantial overlap of their developmental functions. Such pattern likely resides on the duplication history of Hox genes, which probably provided different possible evolutionary routes for achieving the same phenotypic solution.

Conclusions

Our results indicate that convergence occurred at a functional level of Hox genes along three independent origins of aquatic mammals. This conclusion reinforces the idea that different changes in developmental genes may lead to similar phenotypes, probably due to the redundancy provided by the participation of Hox paralogous genes in several developmental functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0682-4) contains supplementary material, which is available to authorized users.

Evolution of mammalian sound localization circuits: A developmental perspective

Localization of sound sources is a central aspect of auditory processing. A unique feature of mammals is the smooth, tonotopically organized extension of the hearing range to high frequencies (HF) above 10kHz, which likely induced positive selection for novel mechanisms of sound localization. How this change in the auditory periphery is accompanied by changes in the central auditory system is unresolved. I will argue that the major VGlut2(+) excitatory projection neurons of sound localization circuits (dorsal cochlear nucleus (DCN), lateral and medial superior olive (LSO and MSO)) represent serial homologs with modifications, thus being paramorphs. This assumption is based on common embryonic origin from an Atoh1(+)/Wnt1(+) cell lineage in the rhombic lip of r5, same cell birth, a fusiform cell morphology, shared genetic components such as Lhx2 and Lhx9 transcription factors, and similar projection patterns. Such a parsimonious evolutionary mechanism likely accelerated the emergence of neurons for sound localization in all three dimensions. Genetic analyses indicate that auditory nuclei in fish, birds, and mammals receive contributions from the same progenitor lineages. Anatomical and physiological differences and the independent evolution of tympanic ears in vertebrate groups, however, argue for convergent evolution of sound localization circuits in tetrapods (amphibians, reptiles, birds, and mammals). These disparate findings are discussed in the context of the genetic architecture of the developing hindbrain, which facilitates convergent evolution. Yet, it will be critical to decipher the gene regulatory networks underlying development of auditory neurons across vertebrates to explore the possibility of homologous neuronal populations.

Convergent Evolution of Head Crests in Two Domesticated Columbids Is Associated with Different Missense Mutations in EphB2

Head crests are important display structures in wild bird species and are also common in domesticated lineages. Many breeds of domestic rock pigeon (Columba livia) have crests of reversed occipital feathers, and this recessive trait is associated with a nonsynonymous coding mutation in the intracellular kinase domain of EphB2 (Ephrin receptor B2). The domestic ringneck dove (Streptopelia risoria) also has a recessive crested morph with reversed occipital feathers, and interspecific crosses between crested doves and pigeons produce crested offspring, suggesting a similar genetic basis for this trait in both species. We therefore investigated EphB2 as a candidate for the head crest phenotype of ringneck doves and identified a nonsynonymous coding mutation in the intracellular kinase domain that is significantly associated with the crested morph. This mutation is over 100 amino acid positions away from the crest mutation found in rock pigeons, yet both mutations are predicted to negatively affect the function of ATP-binding pocket. Furthermore, bacterial toxicity assays suggest that "crest" mutations in both species severely impact kinase activity. We conclude that head crests are associated with different mutations in the same functional domain of the same gene in two different columbid species, thereby representing striking evolutionary convergence in morphology and molecules.

Are Convergent and Parallel Amino Acid Substitutions in Protein Evolution More Prevalent Than Neutral Expectations?

Convergent and parallel amino acid substitutions in protein evolution, collectively referred to as molecular convergence here, have small probabilities under neutral evolution. For this reason, molecular convergence is commonly viewed as evidence for similar adaptations of different species. The surge in the number of reports of molecular convergence in the last decade raises the intriguing question of whether molecular convergence occurs substantially more frequently than expected under neutral evolution. We here address this question using all one-to-one orthologous proteins encoded by the genomes of 12 fruit fly species and those encoded by 17 mammals. We found that the expected amount of molecular convergence varies greatly depending on the specific neutral substitution model assumed at each amino acid site and that the observed amount of molecular convergence is explainable by neutral models incorporating site-specific information of acceptable amino acids. Interestingly, the total number of convergent and parallel substitutions between two lineages, relative to the neutral expectation, decreases with the genetic distance between the two lineages, regardless of the model used in computing the neutral expectation. We hypothesize that this trend results from differences in the amino acids acceptable at a given site among different clades of a phylogeny, due to prevalent epistasis, and provide simulation as well as empirical evidence for this hypothesis. Together, our study finds no genomic evidence for higher-than-neutral levels of molecular convergence, but suggests the presence of abundant epistasis that decreases the likelihood of molecular convergence between distantly related lineages.

Determining the Null Model for Detecting Adaptive Convergence from Genomic Data: A Case Study using Echolocating Mammals

Convergent evolution occurs when the same trait arises independently in multiple lineages. In most cases of phenotypic convergence such transitions are adaptive, so finding the underlying molecular causes of convergence can provide insight into the process of adaptation. Convergent evolution at the genomic level also lends itself to study by comparative methods, although molecular convergence can also occur by chance, adding noise to this process. Parker et al. studied convergence across the genomes of several mammals, including echolocating bats and dolphins (Parker J, Tsagkogeorga G, Cotton JA, Liu Y, Provero P, Stupka E, Rossiter SJ. 2013. Genome-wide signatures of convergent evolution in echolocating mammals. Nature 502:228-231). On the basis of a null distribution of site-specific likelihood support (SSLS) generated using simulated topologies, they concluded that there was evidence for genome-wide adaptive convergence between echolocating taxa. Here, we demonstrate that methods based on SSLS do not adequately measure convergence, and reiterate the use of an empirical null model that directly compares convergent substitutions between all pairs of species. We find that when the proper comparisons are made there is no surprising excess of convergence between echolocating mammals, even in sensory genes.

No genome-wide protein sequence convergence for echolocation

Toothed whales and two groups of bats independently acquired echolocation, the ability to locate and identify objects by reflected sound. Echolocation requires physiologically complex and coordinated vocal, auditory, and neural functions, but the molecular basis of the capacity for echolocation is not well understood. A recent study suggested that convergent amino acid substitutions widespread in the proteins of echolocators underlay the convergent origins of mammalian echolocation. Here, we show that genomic signatures of molecular convergence between echolocating lineages are generally no stronger than those between echolocating and comparable nonecholocating lineages. The same is true for the group of 29 hearing-related proteins claimed to be enriched with molecular convergence. Reexamining the previous selection test reveals several flaws and invalidates the asserted evidence for adaptive convergence. Together, these findings indicate that the reported genomic signatures of convergence largely reflect the background level of sequence convergence unrelated to the origins of echolocation.