RERconverge: an R package for associating evolutionary rates with convergent traits

ERCnet: Phylogenomic Prediction of Interaction Networks in the Presence of Gene Duplication

Assigning gene function from genome sequences is a rate-limiting step in molecular biology research. A protein's position within an interaction network can potentially provide insights into its molecular mechanisms. Phylogenetic analysis of evolutionary rate covariation (ERC) in protein sequence has been shown to be effective for large-scale prediction of functional relationships and interactions. However, gene duplication, gene loss, and other sources of phylogenetic incongruence are barriers for analyzing ERC on a genome-wide basis. Here, we developed ERCnet, a bioinformatic program designed to overcome these challenges, facilitating efficient all-versus-all ERC analyses for large protein sequence datasets. We simulated proteome datasets and found that ERCnet achieves combined false positive and negative error rates well below 10% and that our novel "branch-by-branch" length measurements outperforms "root-to-tip" approaches in most cases, offering a valuable new strategy for performing ERC. We also compiled a sample set of 35 angiosperm genomes to test the performance of ERCnet on empirical data, including its sensitivity to user-defined analysis parameters such as input dataset size and branch-length measurement strategy. We investigated the overlap between ERCnet runs with different species samples to understand how species number and composition affect predicted interactions and to identify the protein sets that consistently exhibit ERC across angiosperms. Our systematic exploration of the performance of ERCnet provides a roadmap for design of future ERC analyses to predict functional interactions in a wide array of genomic datasets. ERCnet code is freely available at https://github.com/EvanForsythe/ERCnet.

The Genomic Landscape, Causes, and Consequences of Extensive Phylogenomic Discordance in Murine Rodents

A species tree is a central concept in evolutionary biology whereby a single branching phylogeny reflects relationships among species. However, the phylogenies of different genomic regions often differ from the species tree. Although tree discordance is widespread in phylogenomic studies, we still lack a clear understanding of how variation in phylogenetic patterns is shaped by genome biology or the extent to which discordance may compromise comparative studies. We characterized patterns of phylogenomic discordance across the murine rodents-a large and ecologically diverse group that gave rise to the laboratory mouse and rat model systems. Combining recently published linked-read genome assemblies for seven murine species with other available rodent genomes, we first used ultraconserved elements (UCEs) to infer a robust time-calibrated species tree. We then used whole genomes to examine finer-scale patterns of discordance across ∼12 million years of divergence. We found that proximate chromosomal regions tended to have more similar phylogenetic histories. There was no clear relationship between local tree similarity and recombination rates in house mice, but we did observe a correlation between recombination rates and average similarity to the species tree. We also detected a strong influence of linked selection whereby purifying selection at UCEs led to appreciably less discordance. Finally, we show that assuming a single species tree can result in substantial deviation from the results with gene trees when testing for positive selection under different models. Collectively, our results highlight the complex relationship between phylogenetic inference and genome biology and underscore how failure to account for this complexity can mislead comparative genomic studies.

A phylogenetic approach to comparative genomics

Comparative genomics, whereby the genomes of different species are compared, has the potential to address broad and fundamental questions at the intersection of genetics and evolution. However, species, genomes and genes cannot be considered as independent data points within statistical tests. Closely related species tend to be similar because they share genes by common descent, which must be accounted for in analyses. This problem of non-independence may be exacerbated when examining genomes or genes but can be addressed by applying phylogeny-based methods to comparative genomic analyses. Here, we review how controlling for phylogeny can change the conclusions of comparative genomics studies. We address common questions on how to apply these methods and illustrate how they can be used to test causal hypotheses. The combination of rapidly expanding genomic datasets and phylogenetic comparative methods is set to revolutionize the biological insights possible from comparative genomic studies.

Convergent relaxation of molecular constraint in herbivores reveals the changing role of liver and kidney functions across mammalian diets

Mammalia comprises a great diversity of diet types and associated adaptations. An understanding of the genomic mechanisms underlying these adaptations may offer insights for improving human health. Comparative genomic studies of diet that employ taxonomically restricted analyses or simplified diet classifications may suffer reduced power to detect molecular convergence associated with diet evolution. Here, we use a quantitative carnivory score-indicative of the amount of animal protein in the diet-for 80 mammalian species to detect significant correlations between the relative evolutionary rates of genes and changes in diet. We have identified six genes-ACADSB, CLDN16, CPB1, PNLIP, SLC13A2, and SLC14A2-that experienced significant changes in evolutionary constraint alongside changes in carnivory score, becoming less constrained in lineages evolving more herbivorous diets. We further consider the biological functions associated with diet evolution and observe that pathways related to amino acid and lipid metabolism, biological oxidation, and small molecule transport experienced reduced purifying selection as lineages became more herbivorous. Liver and kidney functions show similar patterns of constraint with dietary change. Our results indicate that these functions are important for the consumption of animal matter and become less important with the evolution of increasing herbivory. So, genes expressed in these tissues experience a relaxation of evolutionary constraint in more herbivorous lineages.

Comparative genomics of the sexually transmitted parasite Trichomonas vaginalis reveals relaxed and convergent evolution and genes involved in spillover from birds to humans

Trichomonas vaginalis is the causative agent of the venereal disease trichomoniasis which infects men and women globally and is associated with serious outcomes during pregnancy and cancers of the human reproductive tract. Trichomonads parasitize a range of hosts in addition to humans including birds, livestock, and domesticated animals. Recent genetic analysis of trichomonads recovered from columbid birds has provided evidence that these parasite species undergo frequent host-switching, and that a current epoch spillover event from columbids likely gave rise to T. vaginalis in humans. We undertook a comparative evolutionary genomics study of seven trichomonad species, generating chromosome-scale reference genomes for T. vaginalis and its avian sister species Trichomonas stableri, and assemblies of five other species that infect birds and mammals. Human-infecting trichomonad lineages have undergone recent and convergent genome size expansions compared to their avian sister species, and the major contributor to their increased genome size is increased repeat expansions, especially multicopy gene families and transposable elements, with genetic drift likely a driver due to relaxed selection. Trichomonads have independently host-switched twice from birds to humans, and genes implicated in the transition to the human host include those associated with host tissue adherence and phagocytosis, extracellular vesicles, and CAZyme virulence factors.

Sperm competition intensity shapes divergence in both sperm morphology and reproductive genes across murine rodents

It remains unclear how variation in the intensity of sperm competition shapes phenotypic and molecular evolution across clades. Mice and rats in the subfamily Murinae are a rapid radiation exhibiting incredible diversity in sperm morphology and production. We combined phenotypic and genomic data to perform phylogenetic comparisons of male reproductive traits and genes across 78 murine species. We identified several shifts towards smaller relative testes mass (RTM), presumably reflecting reduced sperm competition. Several sperm traits were associated with RTM, suggesting that mating system evolution selects for convergent suites of traits related to sperm competitive ability. We predicted that sperm competition would also drive more rapid molecular divergence in species with large testes. Contrary to this, we found that many spermatogenesis genes evolved more rapidly in species with smaller RTM due to relaxed purifying selection. While some reproductive genes evolved rapidly under recurrent positive selection, relaxed selection played a greater role in underlying rapid evolution in small testes species. Our work demonstrates that postcopulatory sexual selection can impose strong purifying selection shaping the evolution of male reproduction and that broad patterns of molecular evolution may help identify genes that contribute to male fertility.

Convergent Evolution Associated with the Loss of Developmental Diapause May Promote Extended Lifespan in Bees

Diapause has long been proposed to play a significant role in the evolution of eusociality in Hymenoptera. Recent studies have shown that shifts in the diapause stage precede social evolution in wasps and bees; however, the genomic basis remains unknown. Given the overlap in molecular pathways that regulate diapause and lifespan, we hypothesized that the evolutionary loss of developmental diapause may lead to extended lifespan among adults, which is a prerequisite for the evolution of eusociality. To test whether the loss of prepupal diapause is followed by genomic changes associated with lifespan extension, we compared 27 bee genomes with or without prepupal diapause. Our results point to several potential mechanisms for lifespan extension in species lacking prepupal diapause, including the loss of the growth hormone PTTH and its receptor TORSO, along with convergent selection in genes known to regulate lifespan in animals. Specifically, we observed purifying selection of prolongevity genes and relaxed selection of antilongevity genes within the IIS/TOR pathway in species that have lost prepupal diapause. Changes in selection pressures on this pathway may lead to the evolution of new phenotypes, such as lifespan extension and altered responses to nutritional signals that are crucial for social evolution.

Reticulate evolution: Detection and utility in the phylogenomics era

Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.

Extensive longevity and DNA virus-driven adaptation in nearctic Myotis bats

The genus Myotis is one of the largest clades of bats, and exhibits some of the most extreme variation in lifespans among mammals alongside unique adaptations to viral tolerance and immune defense. To study the evolution of longevity-associated traits and infectious disease, we generated near-complete genome assemblies and cell lines for 8 closely related species of Myotis. Using genome-wide screens of positive selection, analyses of structural variation, and functional experiments in primary cell lines, we identify new patterns of adaptation contributing to longevity, cancer resistance, and viral interactions in bats. We find that Myotis bats have some of the most significant variation in cancer risk across mammals and demonstrate a unique DNA damage response in primary cells of the long-lived M. lucifugus. We also find evidence of abundant adaptation in response to DNA viruses - but not RNA viruses - in Myotis and other bats in sharp contrast with other mammals, potentially contributing to the role of bats as reservoirs of zoonoses. Together, our results demonstrate how genomics and primary cells derived from diverse taxa uncover the molecular bases of extreme adaptations in non-model organisms.

A versatile pipeline to identify convergently lost ancestral conserved fragments associated with convergent evolution of vocal learning

Molecular convergence in convergently evolved lineages provides valuable insights into the shared genetic basis of converged phenotypes. However, most methods are limited to coding regions, overlooking the potential contribution of regulatory regions. We focused on the independently evolved vocal learning ability in multiple avian lineages, and developed a whole-genome-alignment-free approach to identify genome-wide Convergently Lost Ancestral Conserved fragments (CLACs) in these lineages, encompassing noncoding regions. We discovered 2711 CLACs that are overrepresented in noncoding regions. Proximal genes of these CLACs exhibit significant enrichment in neurological pathways, including glutamate receptor signaling pathway and axon guidance pathway. Moreover, their expression is highly enriched in brain tissues associated with speech formation. Notably, several have known functions in speech and language learning, including ROBO family, SLIT2, GRIN1, and GRIN2B. Additionally, we found significantly enriched motifs in noncoding CLACs, which match binding motifs of transcriptional factors involved in neurogenesis and gene expression regulation in brain. Furthermore, we discovered 19 candidate genes that harbor CLACs in both human and multiple avian vocal learning lineages, suggesting their potential contribution to the independent evolution of vocal learning in both birds and humans.

Practical Guidance and Workflows for Identifying Fast Evolving Non-Coding Genomic Elements Using PhyloAcc

Comparative genomics provides ample ways to study genome evolution and its relationship to phenotypic traits. By developing and testing alternate models of evolution throughout a phylogeny, one can estimate rates of molecular evolution along different lineages in a phylogeny and link these rates with observations in extant species, such as convergent phenotypes. Pipelines for such work can help identify when and where genomic changes may be associated with, or possibly influence, phenotypic traits. We recently developed a set of models called PhyloAcc, using a Bayesian framework to estimate rates of nucleotide substitution on different branches of a phylogenetic tree and evaluate their association with pre-defined or estimated phenotypic traits. PhyloAcc-ST and PhyloAcc-GT both allow users to define a priori a set of target lineages and then compare different models to identify loci accelerating in one or more target lineages. Whereas ST considers only one species tree across all input loci, GT considers alternate topologies for every locus. PhyloAcc-C simultaneously models molecular rates and rates of continuous trait evolution, allowing the user to ask whether the two are associated. Here, we describe these models and provide tips and workflows on how to prepare the input data and run PhyloAcc.

RERconverge Expansion: Using Relative Evolutionary Rates to Study Complex Categorical Trait Evolution

Comparative genomics approaches seek to associate molecular evolution with the evolution of phenotypes across a phylogeny. Many of these methods lack the ability to analyze non-ordinal categorical traits with more than two categories. To address this limitation, we introduce an expansion to RERconverge that associates shifts in evolutionary rates with the convergent evolution of categorical traits. The categorical RERconverge expansion includes methods for performing categorical ancestral state reconstruction, statistical tests for associating relative evolutionary rates with categorical variables, and a new method for performing phylogeny-aware permutations, "permulations", on categorical traits. We demonstrate our new method on a three-category diet phenotype, and we compare its performance to binary RERconverge analyses and two existing methods for comparative genomic analyses of categorical traits: phylogenetic simulations and a phylogenetic signal based method. We present an analysis of how the categorical permulations scale with the number of species and the number of categories included in the analysis. Our results show that our new categorical method outperforms phylogenetic simulations at identifying genes and enriched pathways significantly associated with the diet phenotypes and that the categorical ancestral state reconstruction drives an improvement in our ability to capture diet-related enriched pathways compared to binary RERconverge when implemented without user input on phenotype evolution. The categorical expansion to RERconverge will provide a strong foundation for applying the comparative method to categorical traits on larger data sets with more species and more complex trait evolution than have previously been analyzed.

Sperm competition intensity shapes divergence in both sperm morphology and reproductive genes across murine rodents

It remains unclear how variation in the intensity of sperm competition shapes phenotypic and molecular evolution across clades. Mice and rats in the subfamily Murinae are a rapid radiation exhibiting incredible diversity in sperm morphology and production. We combined phenotypic and genomic data to perform phylogenetic comparisons of male reproductive traits and genes across 78 murine species. We identified several shifts towards smaller relative testes mass, presumably reflecting reduced sperm competition. Several sperm traits were associated with relative testes mass, suggesting that mating system evolution selects for convergent suites of traits related to sperm competitive ability. We predicted that sperm competition would also drive more rapid molecular divergence in species with large testes. Contrary to this, we found that many spermatogenesis genes evolved more rapidly in species with smaller relative testes mass due to relaxed purifying selection. While some reproductive genes evolved rapidly under recurrent positive selection, relaxed selection played a greater role in underlying rapid evolution in small testes species. Our work demonstrates that postcopulatory sexual selection can impose strong purifying selection shaping the evolution of male reproduction, and that broad patterns of molecular evolution may help identify genes that contribute to male fertility.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these 2 fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we lay out a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., genome-wide association studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur analytically and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study, we re-examine an analysis testing for coevolution of expression levels between genes across a fungal phylogeny and show that including eigenvectors of the covariance matrix as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

Distinct Genes with Similar Functions Underlie Convergent Evolution in Myotis Bat Ecomorphs

Convergence offers an opportunity to explore to what extent evolution can be predictable when genomic composition and environmental triggers are similar. Here, we present an emergent model system to study convergent evolution in nature in a mammalian group, the bat genus Myotis. Three foraging strategies-gleaning, trawling, and aerial hawking, each characterized by different sets of phenotypic features-have evolved independently multiple times in different biogeographic regions in isolation for millions of years. To investigate the genomic basis of convergence and explore the functional genomic changes linked to ecomorphological convergence, we sequenced and annotated 17 new genomes and screened 16,426 genes for positive selection and associations between relative evolutionary rates and foraging strategies across 30 bat species representing all Myotis ecomorphs across geographic regions as well as among sister groups. We identify genomic changes that describe both phylogenetic and ecomorphological trends. We infer that colonization of new environments may have first required changes in genes linked to hearing sensory perception, followed by changes linked to fecundity and development, metabolism of carbohydrates, and heme degradation. These changes may be linked to prey acquisition and digestion and match phylogenetic trends. Our findings also suggest that the repeated evolution of ecomorphs does not always involve changes in the same genes but rather in genes with the same molecular functions such as developmental and cellular processes.

Comparative genomics sheds new light on the convergent evolution of infrared vision in snakes

Infrared vision is a highly specialized sensory system that evolved independently in three clades of snakes. Apparently, convergent evolution occurred in the transient receptor potential ankyrin 1 (TRPA1) proteins of infrared-sensing snakes. However, this gene can only explain how infrared signals are received, and not the transduction and processing of those signals. We sequenced the genome of Xenopeltis unicolor, a key outgroup species of pythons, and performed a genome-wide analysis of convergence between two clades of infrared-sensing snakes. Our results revealed pervasive molecular adaptation in pathways associated with neural development and other functions, with parallel selection on loci associated with trigeminal nerve structural organization. In addition, we found evidence of convergent amino acid substitutions in a set of genes, including TRPA1 and TRPM2. The analysis also identified convergent accelerated evolution in non-coding elements near 12 genes involved in facial nerve structural organization and optic nerve development. Thus, convergent evolution occurred across multiple dimensions of infrared vision in vipers and pythons, as well as amino acid substitutions, non-coding elements, genes and functions. These changes enabled independent groups of snakes to develop and use infrared vision.

Comparative genomics reveals convergent signals associated with the high metabolism and longevity in birds and bats

Birds and bats have long lifespans relative to their body size compared with non-flying animals. However, the genomic basis associated with longer lifespan of flying species despite their higher metabolism was unclear. In this study, we hypothesized that genes involved in the regulation of metabolism and lifespan changed with the acquisition of flight and searched for genes that show specific evolutionary patterns in flying species. As a result, we identified several genes that show different evolutionary rates in bird and bat lineages. Genes in pathways involved in lifespan regulation were conserved in birds, while they evolved at an accelerated rate in bats. We also searched for genes in which convergent amino acid substitutions occurred in birds and bats and found such substitutions in genes involved in cancer, reactive oxygen species control and immunity. Our study revealed genomic changes associated with the acquisition of flight in birds and bats and suggested that multiple genes involved in the regulation of lifespan and metabolism support both high metabolism and longevity in flying species.

The genome of African manatee Trichechus senegalensis reveals secondary adaptation to the aquatic environment

Sirenians exhibit unique aquatic adaptations, showcasing both convergent adaptive features shared with cetaceans and unique characteristics such as cold sensitivity and dense bones. Here, we report a chromosome-level genome of the African manatee (Trichechus senegalensis) with high continuity, completeness, and accuracy. We found that genes associated with osteopetrosis have undergone positive selection (CSF1R and LRRK1) or pseudogenized (FAM111A and IGSF23) in the African manatee, potentially contributing to the dense bone formation. The loss of KCNK18 may have increased their sensitivity to cold water temperatures. Moreover, we identified convergent evolutionary signatures in 392 genes among fully aquatic mammals, primarily enriched in skin or skeletal system development and circadian rhythm, which contributed to the transition from terrestrial to fully aquatic lifestyles. The African manatee currently possesses a small effective population size and low genome-wide heterozygosity. Overall, our study provides genetic resources for understanding the evolutionary characteristics and conservation efforts of this species.

Gene Protein Sequence Evolution Can Predict the Rapid Divergence of Ovariole Numbers in the Drosophila melanogaster Subgroup

Ovaries play key roles in fitness and evolution: they are essential female reproductive structures that develop and house the eggs in sexually reproducing animals. In Drosophila, the mature ovary contains multiple tubular egg-producing structures known as ovarioles. Ovarioles arise from somatic cellular structures in the larval ovary called terminal filaments (TFs), formed by TF cells and subsequently enclosed by sheath (SH) cells. As in many other insects, ovariole number per female varies extensively in Drosophila. At present, however, there is a striking gap of information on genetic mechanisms and evolutionary forces that shape the well-documented rapid interspecies divergence of ovariole numbers. To address this gap, here we studied genes associated with Drosophila melanogaster ovariole number or functions based on recent experimental and transcriptional datasets from larval ovaries, including TFs and SH cells, and assessed their rates and patterns of molecular evolution in five closely related species of the melanogaster subgroup that exhibit species-specific differences in ovariole numbers. From comprehensive analyses of protein sequence evolution (dN/dS), branch-site positive selection, expression specificity (tau), and phylogenetic regressions (phylogenetic generalized least squares), we report evidence of 42 genes that showed signs of playing roles in the genetic basis of interspecies evolutionary change of Drosophila ovariole number. These included the signaling genes upd2 and Ilp5 and extracellular matrix genes vkg and Col4a1, whose dN/dS predicted ovariole numbers among species. Together, we propose a model whereby a set of ovariole-involved gene proteins have an enhanced evolvability, including adaptive evolution, facilitating rapid shifts in ovariole number among Drosophila species.

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/.

Evolutionary Rate Shifts in Coding and Regulatory Regions Underpin Repeated Adaptation to Sulfidic Streams in Poeciliid Fishes

Adaptation to extreme environments often involves the evolution of dramatic physiological changes. To better understand how organisms evolve these complex phenotypic changes, the repeatability and predictability of evolution, and possible constraints on adapting to an extreme environment, it is important to understand how adaptive variation has evolved. Poeciliid fishes represent a particularly fruitful study system for investigations of adaptation to extreme environments due to their repeated colonization of toxic hydrogen sulfide-rich springs across multiple species within the clade. Previous investigations have highlighted changes in the physiology and gene expression in specific species that are thought to facilitate adaptation to hydrogen sulfide-rich springs. However, the presence of adaptive nucleotide variation in coding and regulatory regions and the degree to which convergent evolution has shaped the genomic regions underpinning sulfide tolerance across taxa are unknown. By sampling across seven independent lineages in which nonsulfidic lineages have colonized and adapted to sulfide springs, we reveal signatures of shared evolutionary rate shifts across the genome. We found evidence of genes, promoters, and putative enhancer regions associated with both increased and decreased convergent evolutionary rate shifts in hydrogen sulfide-adapted lineages. Our analysis highlights convergent evolutionary rate shifts in sulfidic lineages associated with the modulation of endogenous hydrogen sulfide production and hydrogen sulfide detoxification. We also found that regions with shifted evolutionary rates in sulfide spring fishes more often exhibited convergent shifts in either the coding region or the regulatory sequence of a given gene, rather than both.

Vocal learning-associated convergent evolution in mammalian proteins and regulatory elements

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

Bringing CAM photosynthesis to the table: Paving the way for resilient and productive agricultural systems in a changing climate

Modern agricultural systems are directly threatened by global climate change and the resulting freshwater crisis. A considerable challenge in the coming years will be to develop crops that can cope with the consequences of declining freshwater resources and changing temperatures. One approach to meeting this challenge may lie in our understanding of plant photosynthetic adaptations and water use efficiency. Plants from various taxa have evolved crassulacean acid metabolism (CAM), a water-conserving adaptation of photosynthetic carbon dioxide fixation that enables plants to thrive under semi-arid or seasonally drought-prone conditions. Although past research on CAM has led to a better understanding of the inner workings of plant resilience and adaptation to stress, successful introduction of this pathway into C3 or C4 plants has not been reported. The recent revolution in molecular, systems, and synthetic biology, as well as innovations in high-throughput data generation and mining, creates new opportunities to uncover the minimum genetic tool kit required to introduce CAM traits into drought-sensitive crops. Here, we propose four complementary research avenues to uncover this tool kit. First, genomes and computational methods should be used to improve understanding of the nature of variations that drive CAM evolution. Second, single-cell 'omics technologies offer the possibility for in-depth characterization of the mechanisms that trigger environmentally controlled CAM induction. Third, the rapid increase in new 'omics data enables a comprehensive, multimodal exploration of CAM. Finally, the expansion of functional genomics methods is paving the way for integration of CAM into farming systems.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these two fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we derive a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., Genome-Wide Association Studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur using analytical theory and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate this by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study of this, we re-examine an analysis testing for co-evolution of expression levels between genes across a fungal phylogeny, and show that including covariance matrix eigenvectors as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

Detection of evolutionary conserved and accelerated genomic regions related to adaptation to thermal niches in Anolis lizards

Understanding the genetic basis for adapting to thermal environments is important due to serious effects of global warming on ectothermic species. Various genes associated with thermal adaptation in lizards have been identified mainly focusing on changes in gene expression or the detection of positively selected genes using coding regions. Only a few comprehensive genome-wide analyses have included noncoding regions. This study aimed to identify evolutionarily conserved and accelerated genomic regions using whole genomes of eight Anolis lizard species that have repeatedly adapted to similar thermal environments in multiple lineages. Evolutionarily conserved genomic regions were extracted as regions with overall sequence conservation (regions with fewer base substitutions) across all lineages compared with the neutral model. Genomic regions that underwent accelerated evolution in the lineage of interest were identified as those with more base substitutions in the target branch than in the entire background branch. Conserved elements across all branches were relatively abundant in "intergenic" genomic regions among noncoding regions. Accelerated regions (ARs) of each lineage contained a significantly greater proportion of noncoding RNA genes than the entire multiple alignment. Common genes containing ARs within 5 kb of their vicinity in lineages with similar thermal habitats were identified. Many genes associated with circadian rhythms and behavior were found in hot-open and cool-shaded habitat lineages. These genes might play a role in contributing to thermal adaptation and assist future studies examining the function of genes involved in thermal adaptation via genome editing.

Evolutionary rate covariation is a reliable predictor of co-functional interactions but not necessarily physical interactions

Co-functional proteins tend to have rates of evolution that covary over time. This correlation between evolutionary rates can be measured over the branches of a phylogenetic tree through methods such as evolutionary rate covariation (ERC), and then used to construct gene networks by the identification of proteins with functional interactions. The cause of this correlation has been hypothesized to result from both compensatory coevolution at physical interfaces and nonphysical forces such as shared changes in selective pressure. This study explores whether coevolution due to compensatory mutations has a measurable effect on the ERC signal. We examined the difference in ERC signal between physically interacting protein domains within complexes compared to domains of the same proteins that do not physically interact. We found no generalizable relationship between physical interaction and high ERC, although a few complexes ranked physical interactions higher than nonphysical interactions. Therefore, we conclude that coevolution due to physical interaction is weak, but present in the signal captured by ERC, and we hypothesize that the stronger signal instead comes from selective pressures on the protein as a whole and maintenance of the general function.

Convergent rates of protein evolution identify novel targets of sexual selection in primates

Sexual selection is the differential reproductive success of individuals, resulting from competition for mates, mate choice, or success in fertilization. In primates, this selective pressure often leads to the development of exaggerated traits which play a role in sexual competition and successful reproduction. In order to gain insight into the mechanisms driving the development of sexually selected traits, we used an unbiased genome-wide approach across 21 primate species to correlate individual rates of protein evolution to relative testes size and sexual dimorphism in body size, 2 anatomical hallmarks of sexual selection in mammals. Among species with presumed high levels of sperm competition, we detected strong conservation of testes-specific proteins responsible for spermatogenesis and ciliary form and function. In contrast, we identified accelerated evolution of female reproductive proteins expressed in the vagina, cervix, and fallopian tubes in these same species. Additionally, we found accelerated protein evolution in lymphoid tissue, indicating that adaptive immune functions may also be influenced by sexual selection. This study demonstrates the distinct complexity of sexual selection in primates revealing contrasting patterns of protein evolution between male and female reproductive tissues.

RERconverge Expansion: Using Relative Evolutionary Rates to Study Complex Categorical Trait Evolution

Comparative genomics approaches seek to associate evolutionary genetic changes with the evolution of phenotypes across a phylogeny. Many of these methods, including our evolutionary rates based method, RERconverge, lack the capability of analyzing non-ordinal, multicategorical traits. To address this limitation, we introduce an expansion to RERconverge that associates shifts in evolutionary rates with the convergent evolution of multi-categorical traits. The categorical RERconverge expansion includes methods for performing categorical ancestral state reconstruction, statistical tests for associating relative evolutionary rates with categorical variables, and a new method for performing phylogenetic permulations on multi-categorical traits. In addition to demonstrating our new method on a three-category diet phenotype, we compare its performance to naive pairwise binary RERconverge analyses and two existing methods for comparative genomic analyses of categorical traits: phylogenetic simulations and a phylogenetic signal based method. We also present a diagnostic analysis of the new permulations approach demonstrating how the method scales with the number of species and the number of categories included in the analysis. Our results show that our new categorical method outperforms phylogenetic simulations at identifying genes and enriched pathways significantly associated with the diet phenotype and that the new ancestral reconstruction drives an improvement in our ability to capture diet-related enriched pathways. Our categorical permulations were able to account for non-uniform null distributions and correct for non-independence in gene rank during pathway enrichment analysis. The categorical expansion to RERconverge will provide a strong foundation for applying the comparative method to categorical traits on larger data sets with more species and more complex trait evolution.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Roles of the crucial mitochondrial DNA in hypertrophic cardiomyopathy prognosis and diagnosis: A review

Mitochondrial DNA is implicated in hypertrophic cardiomyopathy (HCM) development. We aimed to identify valuable mtDNAs that contribute to the development of HCM. Differentially expressed mitochondrial DNAs (DEMGs) between HCM and controls were screened. GO and KEGG functional enrichment analyses were performed, and the optimum genes were explored using the LASSO regression mode and SVM-RFE model. A diagnostic scoring model was constructed and verified using ROC curves. Mitochondria-based subtypes were identified. Immune performance among the subtypes including immune cells, immune checkpoint genes, and HLA family genes was analyzed. Finally, an mRNA-transcription factor (TF)-miRNA network was constructed using Cytoscape software. Twelve DEMGs in HCM were selected. Among them, 6 DEMGs, including PDK4, MGST1, TOMM40, LYPLAL1, GATM, and CPT1B were demonstrated as DEMGs at the point of intersection of Lasso regression and SVM-RFE. The ROC of the model for the training and validation datasets was 0.999 and 0.958, respectively. Two clusters were divided, and 4 immune cell types were significantly different between the 2 clusters, including resting mast cells, macrophages M2, and plasma cells. Nine upregulated KEGG pathways were enriched in cluster 1 vs. cluster 2 including O-glycan biosynthesis, the ErbB signaling pathway, and the GnRH signaling pathway. Meanwhile, 49 down-regulated pathways were enriched such as the toll-like signaling pathway and natural killer cell-mediated cytotoxicity pathway. The 6 gene-based mRNA-TF-miRNA networks included other 133 TFs and 18 miRNAs. Six DEMGs in HCM, including PDK4, MGST1, TOMM40, LYPLAL1, GATM, and CPT1B, can be indicative of HCM prognosis or disease progression.

Molecular Footprints on Osmoregulation-Related Genes Associated with Freshwater Colonization by Cetaceans and Sirenians

The genetic basis underlying adaptive physiological mechanisms has been extensively explored in mammals after colonizing the seas. However, independent lineages of aquatic mammals exhibit complex patterns of secondary colonization in freshwater environments. This change in habitat represents new osmotic challenges, and additional changes in key systems, such as the osmoregulatory system, are expected. Here, we studied the selective regime on coding and regulatory regions of 20 genes related to the osmoregulation system in strict aquatic mammals from independent evolutionary lineages, cetaceans, and sirenians, with representatives in marine and freshwater aquatic environments. We identified positive selection signals in genes encoding the protein vasopressin (AVP) in mammalian lineages with secondary colonization in the fluvial environment and in aquaporins for lineages inhabiting the marine and fluvial environments. A greater number of sites with positive selection signals were found for the dolphin species compared to the Amazonian manatee. Only the AQP5 and AVP genes showed selection signals in more than one independent lineage of these mammals. Furthermore, the vasopressin gene tree indicates greater similarity in river dolphin sequences despite the independence of their lineages based on the species tree. Patterns of distribution and enrichment of Transcription Factors in the promoter regions of target genes were analyzed and appear to be phylogenetically conserved among sister species. We found accelerated evolution signs in genes ACE, AQP1, AQP5, AQP7, AVP, NPP4, and NPR1 for the fluvial mammals. Together, these results allow a greater understanding of the molecular bases of the evolution of genes responsible for osmotic control in aquatic mammals.

A multifunctional Wnt regulator underlies the evolution of rodent stripe patterns

Animal pigment patterns are excellent models to elucidate mechanisms of biological organization. Although theoretical simulations, such as Turing reaction-diffusion systems, recapitulate many animal patterns, they are insufficient to account for those showing a high degree of spatial organization and reproducibility. Here, we study the coat of the African striped mouse (Rhabdomys pumilio) to uncover how periodic stripes form. Combining transcriptomics, mathematical modelling and mouse transgenics, we show that the Wnt modulator Sfrp2 regulates the distribution of hair follicles and establishes an embryonic prepattern that foreshadows pigment stripes. Moreover, by developing in vivo gene editing in striped mice, we find that Sfrp2 knockout is sufficient to alter the stripe pattern. Strikingly, mutants exhibited changes in pigmentation, revealing that Sfrp2 also regulates hair colour. Lastly, through evolutionary analyses, we find that striped mice have evolved lineage-specific changes in regulatory elements surrounding Sfrp2, many of which may be implicated in modulating the expression of this gene. Altogether, our results show that a single factor controls coat pattern formation by acting both as an orienting signalling mechanism and a modulator of pigmentation. More broadly, our work provides insights into how spatial patterns are established in developing embryos and the mechanisms by which phenotypic novelty originates.

The evolution of aging and lifespan

Aging is a nearly inescapable trait among organisms yet lifespan varies tremendously across different species and spans several orders of magnitude in vertebrates alone. This vast phenotypic diversity is driven by distinct evolutionary trajectories and tradeoffs that are reflected in patterns of diversification and constraint in organismal genomes. Age-specific impacts of selection also shape allele frequencies in populations, thus impacting disease susceptibility and environment-specific mortality risk. Further, the mutational processes that spawn this genetic diversity in both germline and somatic cells are strongly influenced by age and life history. We discuss recent advances in our understanding of the evolution of aging and lifespan at organismal, population, and cellular scales, and highlight outstanding questions that remain unanswered.

Genomic signatures of convergent shifts to plunge-diving behavior in birds

Understanding the genetic basis of convergence at broad phylogenetic scales remains a key challenge in biology. Kingfishers (Aves: Alcedinidae) are a cosmopolitan avian radiation with diverse colors, diets, and feeding behaviors-including the archetypal plunge-dive into water. Given the sensory and locomotor challenges associated with air-water transitions, kingfishers offer a powerful opportunity to explore the effects of convergent behaviors on the evolution of genomes and phenotypes, as well as direct comparisons between continental and island lineages. Here, we use whole-genome sequencing of 30 diverse kingfisher species to identify the genomic signatures associated with convergent feeding behaviors. We show that species with smaller ranges (i.e., on islands) have experienced stronger demographic fluctuations than those on continents, and that these differences have influenced the dynamics of molecular evolution. Comparative genomic analyses reveal positive selection and genomic convergence in brain and dietary genes in plunge-divers. These findings enhance our understanding of the connections between genotype and phenotype in a diverse avian radiation.

Live-bearing cockroach genome reveals convergent evolutionary mechanisms linked to viviparity in insects and beyond

Live birth (viviparity) has arisen repeatedly and independently among animals. We sequenced the genome and transcriptome of the viviparous Pacific beetle-mimic cockroach and performed comparative analyses with two other viviparous insect lineages, tsetse flies and aphids, to unravel the basis underlying the transition to viviparity in insects. We identified pathways undergoing adaptive evolution for insects, involved in urogenital remodeling, tracheal system, heart development, and nutrient metabolism. Transcriptomic analysis of cockroach and tsetse flies revealed that uterine remodeling and nutrient production are increased and the immune response is altered during pregnancy, facilitating structural and physiological changes to accommodate and nourish the progeny. These patterns of convergent evolution of viviparity among insects, together with similar adaptive mechanisms identified among vertebrates, highlight that the transition to viviparity requires changes in urogenital remodeling, enhanced tracheal and heart development (corresponding to angiogenesis in vertebrates), altered nutrient metabolism, and shifted immunity in animal systems.

Genomic Signatures Associated with Transitions to Viviparity in Cyprinodontiformes

The transition from oviparity to viviparity has occurred independently over 150 times across vertebrates, presenting one of the most compelling cases of phenotypic convergence. However, whether the repeated, independent evolution of viviparity is driven by redeployment of similar genetic mechanisms and whether these leave a common signature in genomic divergence remains largely unknown. Although recent investigations into the evolution of viviparity have demonstrated striking similarity among the genes and molecular pathways involved across disparate vertebrate groups, quantitative tests for genome-wide convergent have provided ambivalent answers. Here, we investigate the potential role of molecular convergence during independent transitions to viviparity across an order of ray-finned freshwater fish (Cyprinodontiformes). We assembled de novo genomes and utilized publicly available genomes of viviparous and oviparous species to test for molecular convergence across both coding and noncoding regions. We found no evidence for an excess of molecular convergence in amino acid substitutions and in rates of sequence divergence, implying independent genetic changes are associated with these transitions. However, both statistical power and biological confounds could constrain our ability to detect significant correlated evolution. We therefore identified candidate genes with potential signatures of molecular convergence in viviparous Cyprinodontiformes lineages. Motif enrichment and gene ontology analyses suggest transcriptional changes associated with early morphogenesis, brain development, and immunity occurred alongside the evolution of viviparity. Overall, however, our findings indicate that independent transitions to viviparity in these fish are not strongly associated with an excess of molecular convergence, but a few genes show convincing evidence of convergent evolution.

CAAStools: a toolbox to identify and test Convergent Amino Acid Substitutions

MOTIVATION

Coincidence of Convergent Amino Acid Substitutions (CAAS) with phenotypic convergences allow pinpointing genes and even individual mutations that are likely to be associated with trait variation within their phylogenetic context. Such findings can provide useful insights into the genetic architecture of complex phenotypes.

RESULTS

Here we introduce CAAStools, a set of bioinformatics tools to identify and validate CAAS in orthologous protein alignments for predefined groups of species representing the phenotypic values targeted by the user.

AVAILABILITY AND IMPLEMENTATION

CAAStools source code is available at http://github.com/linudz/caastools, along with documentation and examples.

Diverse signatures of convergent evolution in cacti-associated yeasts

Many distantly related organisms have convergently evolved traits and lifestyles that enable them to live in similar ecological environments. However, the extent of phenotypic convergence evolving through the same or distinct genetic trajectories remains an open question. Here, we leverage a comprehensive dataset of genomic and phenotypic data from 1,049 yeast species in the subphylum Saccharomycotina (Kingdom Fungi, Phylum Ascomycota) to explore signatures of convergent evolution in cactophilic yeasts, ecological specialists associated with cacti. We inferred that the ecological association of yeasts with cacti arose independently ~17 times. Using machine-learning, we further found that cactophily can be predicted with 76% accuracy from functional genomic and phenotypic data. The most informative feature for predicting cactophily was thermotolerance, which is likely associated with duplication and altered evolutionary rates of genes impacting the cell envelope in several cactophilic lineages. We also identified horizontal gene transfer and duplication events of plant cell wall-degrading enzymes in distantly related cactophilic clades, suggesting that putatively adaptive traits evolved through disparate molecular mechanisms. Remarkably, multiple cactophilic lineages and their close relatives are emerging human opportunistic pathogens, suggesting that the cactophilic lifestyle-and perhaps more generally lifestyles favoring thermotolerance-may preadapt yeasts to cause human disease. This work underscores the potential of a multifaceted approach involving high throughput genomic and phenotypic data to shed light onto ecological adaptation and highlights how convergent evolution to wild environments could facilitate the transition to human pathogenicity.

PhyloAcc-GT: A Bayesian Method for Inferring Patterns of Substitution Rate Shifts on Targeted Lineages Accounting for Gene Tree Discordance

An important goal of evolutionary genomics is to identify genomic regions whose substitution rates differ among lineages. For example, genomic regions experiencing accelerated molecular evolution in some lineages may provide insight into links between genotype and phenotype. Several comparative genomics methods have been developed to identify genomic accelerations between species, including a Bayesian method called PhyloAcc, which models shifts in substitution rate in multiple target lineages on a phylogeny. However, few methods consider the possibility of discordance between the trees of individual loci and the species tree due to incomplete lineage sorting, which might cause false positives. Here, we present PhyloAcc-GT, which extends PhyloAcc by modeling gene tree heterogeneity. Given a species tree, we adopt the multispecies coalescent model as the prior distribution of gene trees, use Markov chain Monte Carlo (MCMC) for inference, and design novel MCMC moves to sample gene trees efficiently. Through extensive simulations, we show that PhyloAcc-GT outperforms PhyloAcc and other methods in identifying target lineage-specific accelerations and detecting complex patterns of rate shifts, and is robust to specification of population size parameters. PhyloAcc-GT is usually more conservative than PhyloAcc in calling convergent rate shifts because it identifies more accelerations on ancestral than on terminal branches. We apply PhyloAcc-GT to two examples of convergent evolution: flightlessness in ratites and marine mammal adaptations, and show that PhyloAcc-GT is a robust tool to identify shifts in substitution rate associated with specific target lineages while accounting for incomplete lineage sorting.

The genomic landscape, causes, and consequences of extensive phylogenomic discordance in Old World mice and rats

A species tree is a central concept in evolutionary biology whereby a single branching phylogeny reflects relationships among species. However, the phylogenies of different genomic regions often differ from the species tree. Although tree discordance is often widespread in phylogenomic studies, we still lack a clear understanding of how variation in phylogenetic patterns is shaped by genome biology or the extent to which discordance may compromise comparative studies. We characterized patterns of phylogenomic discordance across the murine rodents (Old World mice and rats) - a large and ecologically diverse group that gave rise to the mouse and rat model systems. Combining new linked-read genome assemblies for seven murine species with eleven published rodent genomes, we first used ultra-conserved elements (UCEs) to infer a robust species tree. We then used whole genomes to examine finer-scale patterns of discordance and found that phylogenies built from proximate chromosomal regions had similar phylogenies. However, there was no relationship between tree similarity and local recombination rates in house mice, suggesting that genetic linkage influences phylogenetic patterns over deeper timescales. This signal may be independent of contemporary recombination landscapes. We also detected a strong influence of linked selection whereby purifying selection at UCEs led to less discordance, while genes experiencing positive selection showed more discordant and variable phylogenetic signals. Finally, we show that assuming a single species tree can result in high error rates when testing for positive selection under different models. Collectively, our results highlight the complex relationship between phylogenetic inference and genome biology and underscore how failure to account for this complexity can mislead comparative genomic studies.

The genomic basis of temporal niche evolution in a diurnal rodent

Patterns of diel activity-how animals allocate their activity throughout the 24-h daily cycle-play key roles in shaping the internal physiology of an animal and its relationship with the external environment.1,2,3,4,5 Although shifts in diel activity patterns have occurred numerous times over the course of vertebrate evolution,6 the genomic correlates of such transitions remain unknown. Here, we use the African striped mouse (Rhabdomys pumilio), a species that transitioned from the ancestrally nocturnal diel niche of its close relatives to a diurnal one,7,8,9,10,11 to define patterns of naturally occurring molecular variation in diel niche traits. First, to facilitate genomic analyses, we generate a chromosome-level genome assembly of the striped mouse. Next, using transcriptomics, we show that the switch to daytime activity in this species is associated with a realignment of daily rhythms in peripheral tissues with respect to the light:dark cycle and the central circadian clock. To uncover selection pressures associated with this temporal niche shift, we perform comparative genomic analyses with closely related rodent species and find evidence of relaxation of purifying selection on striped mouse genes in the rod phototransduction pathway. In agreement with this, electroretinogram measurements demonstrate that striped mice have functional differences in dim-light visual responses compared with nocturnal rodents. Taken together, our results show that striped mice have undergone a drastic change in circadian organization and provide evidence that the visual system has been a major target of selection as this species transitioned to a novel temporal niche.

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.

A switch in jaw form-function coupling during the evolution of mammals

The evolutionary shift from a single-element ear, multi-element jaw to a multi-element ear, single-element jaw during the transition to crown mammals marks one of the most dramatic structural transformations in vertebrates. Research on this transformation has focused on mammalian middle-ear evolution, but a mandible comprising only the dentary is equally emblematic of this evolutionary radiation. Here, we show that the remarkably diverse jaw shapes of crown mammals are coupled with surprisingly stereotyped jaw stiffness. This strength-based morphofunctional regime has a genetic basis and allowed mammalian jaws to effectively resist deformation as they radiated into highly disparate forms with markedly distinct diets. The main functional consequences for the mandible of decoupling hearing and mastication were a trade-off between higher jaw stiffness versus decreased mechanical efficiency and speed compared with non-mammals. This fundamental and consequential shift in jaw form-function underpins the ecological and taxonomic diversification of crown mammals. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Relating enhancer genetic variation across mammals to complex phenotypes using machine learning

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Genomic convergence underlying high-altitude adaptation in alpine plants

Evolutionary convergence is one of the most striking examples of adaptation driven by natural selection. However, genomic evidence for convergent adaptation to extreme environments remains scarce. Here, we assembled reference genomes of two alpine plants, Saussurea obvallata (Asteraceae) and Rheum alexandrae (Polygonaceae), with 37,938 and 61,463 annotated protein-coding genes. By integrating an additional five alpine genomes, we elucidated genomic convergence underlying high-altitude adaptation in alpine plants. Our results detected convergent contractions of disease-resistance genes in alpine genomes, which might be an energy-saving strategy for surviving in hostile environments with only a few pathogens present. We identified signatures of positive selection on a set of genes involved in reproduction and respiration (e.g., MMD1, NBS1, and HPR), and revealed signatures of molecular convergence on genes involved in self-incompatibility, cell wall modification, DNA repair and stress resistance, which may underlie adaptation to extreme cold, high ultraviolet radiation and hypoxia environments. Incorporating transcriptomic data, we further demonstrated that genes associated with cuticular wax and flavonoid biosynthetic pathways exhibit higher expression levels in leafy bracts, shedding light on the genetic mechanisms of the adaptive "greenhouse" morphology. Our integrative data provide novel insights into convergent evolution at a high-taxonomic level, aiding in a deep understanding of genetic adaptation to complex environments.

The GRN concept as a guide for evolutionary developmental biology

Organismal phenotypes result largely from inherited developmental programs, usually executed during embryonic and juvenile life stages. These programs are not blank slates onto which natural selection can draw arbitrary forms. Rather, the mechanisms of development play an integral role in shaping phenotypic diversity and help determine the evolutionary trajectories of species. Modern evolutionary biology must, therefore, account for these mechanisms in both theory and in practice. The gene regulatory network (GRN) concept represents a potent tool for achieving this goal whose utility has grown in tandem with advances in "omic" technologies and experimental techniques. However, while the GRN concept is widely utilized, it is often less clear what practical implications it has for conducting research in evolutionary developmental biology. In this Perspective, we attempt to provide clarity by discussing how experiments and projects can be designed in light of the GRN concept. We first map familiar biological notions onto the more abstract components of GRN models. We then review how diverse functional genomic approaches can be directed toward the goal of constructing such models and discuss current methods for functionally testing evolutionary hypotheses that arise from them. Finally, we show how the major steps of GRN model construction and experimental validation suggest generalizable workflows that can serve as a scaffold for project design. Taken together, the practical implications that we draw from the GRN concept provide a set of guideposts for studies aiming at unraveling the molecular basis of phenotypic diversity.

Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles

Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.

Convergent genomics of longevity in rockfishes highlights the genetics of human life span variation

Longevity is a defining, heritable trait that varies dramatically between species. To resolve the genetic regulation of this trait, we have mined genomic variation in rockfishes, which range in longevity from 11 to over 205 years. Multiple shifts in rockfish longevity have occurred independently and in a short evolutionary time frame, thus empowering convergence analyses. Our analyses reveal a common network of genes under convergent evolution, encompassing established aging regulators such as insulin signaling, yet also identify flavonoid (aryl-hydrocarbon) metabolism as a pathway modulating longevity. The selective pressures on these pathways indicate the ancestral state of rockfishes was long lived and that the changes in short-lived lineages are adaptive. These pathways were also used to explore genome-wide association studies of human longevity, identifying the aryl-hydrocarbon metabolism pathway to be significantly associated with human survival to the 99th percentile. This evolutionary intersection defines and cross-validates a previously unappreciated genetic architecture that associates with the evolution of longevity across vertebrates.

Genomic signatures of recent convergent transitions to social life in spiders

The transition from solitary to social life is a major phenotypic innovation, but its genetic underpinnings are largely unknown. To identify genomic changes associated with this transition, we compare the genomes of 22 spider species representing eight recent and independent origins of sociality. Hundreds of genes tend to experience shifts in selection during the repeated transition to social life. These genes are associated with several key functions, such as neurogenesis, behavior, and metabolism, and include genes that previously have been implicated in animal social behavior and human behavioral disorders. In addition, social species have elevated genome-wide rates of molecular evolution associated with relaxed selection caused by reduced effective population size. Altogether, our study provides unprecedented insights into the genomic signatures of social evolution and the specific genetic changes that repeatedly underpin the evolution of sociality. Our study also highlights the heretofore unappreciated potential of transcriptomics using ethanol-preserved specimens for comparative genomics and phylotranscriptomics.