Evolutionary repeatability of chromatin accessibility in the hippocampus between echolocating bats and soft-furred tree mice

Echolocation is a compelling example of behavioral convergence in mammals and has evolved independently in at least five mammalian lineages. Despite increasing evidence of convergent evolution in protein-coding regions among echolocating mammals, the evolution of gene regulatory activity underlying mammalian echolocation remains undetermined. Here, we systematically investigate chromatin accessibility and transcriptomic data from the hippocampi of echolocating bats and soft-furred tree mice, as the hippocampus plays a crucial role in processing echolocation signals. Our results showed significant evolutionary repeatability in accessible chromatin regions among echolocating mammals compared with their non-echolocating counterparts. An analysis that paired gene expression with chromatin accessibility revealed that echolocating mammals have more complex gene regulatory networks in the hippocampus than non-echolocating mammals. This complexity is primarily attributed to the observed evolutionary repeatability in accessible chromatin regions among echolocating mammals. The shared gene regulations among echolocating mammals are particularly enriched in functional terms related to synaptic function, which is supported by the higher synaptic and mitochondrial densities in soft-furred tree mice than in laboratory mice. These findings demonstrate significant evolutionary repeatability of gene regulatory activity in the hippocampus among echolocating mammalian species, suggesting crucial contributions of gene regulatory activities to convergent phenotypes in mammals.

Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES

Current genome sequencing initiatives across a wide range of life forms offer significant potential to enhance our understanding of evolutionary relationships and support transformative biological and medical applications. Species trees play a central role in many of these applications; however, despite the widespread availability of genome assemblies, accurate inference of species trees remains challenging due to the limited automation, substantial domain expertise, and computational resources required by conventional methods. To address this limitation, we present ROADIES, a fully automated pipeline to infer species trees starting from raw genome assemblies. In contrast to the prominent approach, ROADIES incorporates a unique strategy of randomly sampling segments of the input genomes to generate gene trees. This eliminates the need for predefining a set of loci, limiting the analyses to a fixed number of genes, and performing the cumbersome gene annotation and/or whole genome alignment steps. ROADIES also eliminates the need to infer orthology by leveraging existing discordance-aware methods that allow multicopy genes. Using the genomic datasets from large-scale sequencing efforts across four diverse life forms (placental mammals, pomace flies, birds, and budding yeasts), we show that ROADIES infers species trees that are comparable in quality to the state-of-the-art studies but in a fraction of the time and effort, including on challenging datasets with rampant gene tree discordance and complex polyploidy. With its speed, accuracy, and automation, ROADIES has the potential to vastly simplify species tree inference, making it accessible to a broader range of scientists and applications.

Unveiling cell-type-specific mode of evolution in comparative single-cell expression data

While methodology for determining the mode of evolution in coding sequences has been well established, evaluation of adaptation events in emerging types of phenotype data needs further development. Here we propose an analysis framework (expression variance decomposition, EVaDe) for comparative single-cell expression data based on phenotypic evolution theory. After decomposing the gene expression variance into separate components, we use two strategies to identify genes exhibiting large between-taxon expression divergence and small within-cell-type expression noise in certain cell types, attributing this pattern to putative adaptive evolution. In a dataset of primate prefrontal cortex, we find that such human-specific key genes enrich with neurodevelopment-related functions, while most other genes exhibit neutral evolution patterns. Specific neuron types are found to harbor more of these key genes than other cell types, thus likely to have experienced more extensive adaptation. Reassuringly, at molecular sequence level, the key genes are significantly associated with the rapidly evolving conserved non-coding elements. An additional case analysis comparing the naked mole-rat (NMR) with the mouse suggests that innate-immunity-related genes and cell types have undergone putative expression adaptation in NMR. Overall, the EVaDe framework may effectively probe adaptive evolution mode in single-cell expression data.

Multiple transitions to high l-DOPA 4,5-dioxygenase activity reveal molecular pathways to convergent betalain pigmentation in Caryophyllales

Many specialized metabolic pathways have evolved convergently in plants, but distinguishing multiple origins from alternative evolutionary scenarios can be difficult. Here, we explore the evolution of l-3,4-dihydroxyphenylalanine (l-DOPA) 4,5-dioxygenase (DODA) enzymes to better resolve the convergent evolution of the betalain biosynthetic pathway within the flowering plant order Caryophyllales. We use yeast-based heterologous assays to quantify enzymatic activity of extant proteins and then employ ancestral sequence reconstruction to resurrect and assay ancestral DODA enzymes. We use a combination of ancestral sequence reconstruction, model-based methods, and structural modelling to describe patterns of molecular convergence. We confirm that high l-DOPA 4,5-dioxygenase activity is polyphyletic and show that high activity DODAs evolved at least three times from ancestral proteins with low activity. We show that molecular convergence is concentrated proximally to the binding pockets but also appears distally to active sites. Moreover, our analysis also suggests that many unique and divergent substitutions contribute to the evolution of DODA. Given the key role of DODA in betalain biosynthesis, our analysis further supports the convergent origins of betalains and illustrates how the iterative evolution of betalain biosynthesis has drawn on a complex mixture of convergent, divergent, and unique variation.

Drivers of avian genomic change revealed by evolutionary rate decomposition

Modern birds have diversified into a striking array of forms, behaviours and ecological roles. Analyses of molecular evolutionary rates can reveal the links between genomic and phenotypic change1-4, but disentangling the drivers of rate variation at the whole-genome scale has been difficult. Using comprehensive estimates of traits and evolutionary rates across a family-level phylogeny of birds5,6, we find that genome-wide mutation rates across lineages are predominantly explained by clutch size and generation length, whereas rate variation across genes is driven by the content of guanine and cytosine. Here, to find the subsets of genes and lineages that dominate evolutionary rate variation in birds, we estimated the influence of individual lineages on decomposed axes of gene-specific evolutionary rates. We find that most of the rate variation occurs along recent branches of the tree, associated with present-day families of birds. Additional tests on axes of rate variation show rapid changes in microchromosomes immediately after the Cretaceous-Palaeogene transition. These apparent pulses of evolution are consistent with major changes in the genetic machineries for meiosis, heart performance, and RNA splicing, surveillance and translation, and correlate with the ecological diversity reflected in increased tarsus length. Collectively, our analyses paint a nuanced picture of avian evolution, revealing that the ancestors of the most diverse lineages of birds underwent major genomic changes related to mutation, gene usage and niche expansion in the early Palaeogene period.

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Low-dimensional genotype-fitness mapping across divergent environments suggests a limiting functions model of fitness

A central goal in evolutionary biology is to be able to predict the effect of a genetic mutation on fitness. This is a major challenge because fitness depends both on phenotypic changes due to the mutation, and how these phenotypes map onto fitness in a particular environment. Genotype, phenotype, and environment spaces are all extremely complex, rendering bottom-up prediction unlikely. Here we show, using a large collection of adaptive yeast mutants, that fitness across a set of lab environments can be well-captured by top-down, low-dimensional linear models that generate abstract genotype-phenotype-fitness maps. We find that these maps are low-dimensional not only in the environment where the adaptive mutants evolved, but also in more divergent environments. We further find that the genotype-phenotype-fitness spaces implied by these maps overlap only partially across environments. We argue that these patterns are consistent with a "limiting functions" model of fitness, whereby only a small number of limiting functions can be modified to affect fitness in any given environment. The pleiotropic side-effects on non-limiting functions are effectively hidden from natural selection locally, but can be revealed globally. These results combine to emphasize the importance of environmental context in genotype-phenotype-fitness mapping, and have implications for the predictability and trajectory of evolution in complex environments.

Multi-level genomic convergence of secondary aquatic adaptation in marine mammals

Marine mammals provide a valuable model for studying the molecular basis of convergent evolution during secondary aquatic adaptation. Using multi-omics data and functional experiments, including CRISPR-Cas9 mouse models and luciferase reporter assays, this study explored the molecular mechanisms driving this transition across coding regions, regulatory elements, and genomic architecture. Convergent amino acid substitutions in APPL1 P378L and NEIL1 E71G were found to promote lipid accumulation and suppress cancer cell proliferation, likely contributing to the evolution of extensive blubber layers and cancer resistance. Convergently evolved conserved non-exonic elements (CNEs) and lineage-specific regulatory variations were shown to influence the activity of nearby genes (e.g., NKX3-2, SOX9, HAND2), shaping cetacean limb phenotypes. Additionally, convergent shifts in topologically associating domains (TADs) across cetaceans and pinnipeds were implicated in the regulation of ASXL3 and FAM43B expression, playing a role in the formation of thickened blubber layers and mitigating cancer susceptibility. Structural variations within conserved TADs were associated with the expression of neuronal genes, including NUP153 and ID4, potentially driving cognitive and social adaptations. These findings provide novel insights into the molecular foundations of the convergent evolution of secondary aquatic adaptations in mammals.

CASTER: Direct species tree inference from whole-genome alignments

Genomes contain mosaics of discordant evolutionary histories, challenging the accurate inference of the tree of life. Although genome-wide data are routinely used for discordance-aware phylogenomic analyses, because of modeling and scalability limitations, the current practice leaves out large chunks of genomes. As more high-quality genomes become available, we urgently need discordance-aware methods to infer the tree directly from a multiple genome alignment. In this study, we introduce Coalescence-Aware Alignment-Based Species Tree Estimator (CASTER), a theoretically justified site-based method that eliminates the need to predefine recombination-free loci. CASTER is scalable to hundreds of mammalian whole genomes. We demonstrate the accuracy and scalability of CASTER in simulations that include recombination and apply CASTER to several biological datasets, showing that its per-site scores can reveal both biological and artifactual patterns of discordance across the genome.

Comparative genomics uncovers evolutionary drivers of locust migratory adaptation

BACKGROUND

Locust migration is one of the main causes of locust plagues. While existing research has highlighted the adaptive migratory capabilities of locusts, the evolutionary patterns of their migration remain elusive. This study aims to explore these evolutionary patterns of locust migratory behavior at the genomic level. To achieve this, we conducted comparative genomics analysis using genomic data from 10 locust species with diverse migratory tendencies.

RESULTS

We identified 1064 genes showing signatures of positive selection in five migratory locust species using a dN/dS model. The BUSTED-PH model revealed 116 genes associated with migratory phenotypes. Gene ontology enrichment analysis indicated that these genes were predominantly related to metabolism and mitochondria-related pathways through both methods. Additionally, the evolutionary rate (RER) analysis between migratory and non-migratory locusts revealed significant divergence in energy metabolism pathways. Notably, of the genes analyzed, the SETX gene consistently showed evidence of positive selection across all five migratory species.

CONCLUSIONS

The findings suggest that the evolution of migratory behavior is associated with increased selective pressure on metabolism and mitochondria-related pathways. Hundreds of genes undergo selective changes during repetitive transitions to migratory behavior. These findings enhance our understanding of the genetic and phenotypic relationships underlying different locust migratory behaviors, providing important data for understanding the biological mechanisms behind locust outbreaks.

Palaeognath birds innovate to solve a novel foraging problem

The ability to innovate implies flexible cognition, and is used as a broad metric of intelligence. Innovation in birds has been intensively studied in the larger and more taxonomically diverse Neognathae clade (particularly crows and parrots) and overlooked in the smaller and more ancestral Palaeognathae clade. The current study provides the first known evidence of technical innovation in palaeognath birds. We tested the ability of nine individuals of three species to move a hole towards a chamber to access a food reward. This problem was different to traditional innovation puzzle-boxes where an obstacle is moved away from a food chamber. Three emus and one rhea produced a wheel-turning innovation, moving the hole in the most efficient direction (closer to the nearest food item) in 90% of cases. One rhea dismantled the task twice by removing the central bolt, which we suggest is a second type of innovation, and it did not persist once they innovated the wheel turning solution. Ostriches did not innovate. We classify innovation in palaeognaths as low level/simplistic, relying on general exploration and asocial trial and error learning. Our research suggests that technical innovation may have evolved far earlier in birds than previously thought, and palaeognath birds are a compelling taxonomic group for further cognitive research.

Molecular Innovations Shaping Beak Morphology in Birds

The beak, a pivotal evolutionary trait characterized by high morphological diversity and plasticity, has enabled birds to survive mass extinction events and subsequently radiate into diverse ecological niches worldwide. This remarkable ecological adaptability underscores the importance of uncovering the molecular mechanisms shaping avian beak morphology, particularly benefiting from the rapidly advancing archives of genomics and epigenomics. We review the latest advancements in understanding how genetic and epigenetic innovations control or regulate beak development and drive beak morphological adaptation and diversification over the past two decades. We conclude with several recommendations for future endeavors, expanding to more bird lineages, with a focus on beak shape and the lower beak, and conducting functional experiments. By directing research efforts toward these aspects and integrating advanced omics techniques, the complex molecular mechanisms involved in avian beak evolution and morphogenesis will be deeply interpreted.

Avian Island Radiations Shed Light on the Dynamics of Adaptive and Nonadaptive Radiation

Understanding the mechanisms underlying species formation and differentiation is a central goal of evolutionary biology and a formidable challenge. This understanding can provide valuable insights into the origins of the astonishing diversity of organisms living on our planet. Avian evolutionary radiations on islands have long fascinated biologists as they provide the ideal variation to study the ecological and evolutionary forces operating on the continuum between incipient lineages to complete speciation. In this review, we summarize the key insights gained from decades of research on adaptive and nonadaptive radiations of both extant and extinct insular bird species. We present a new comprehensive global list of potential avian radiations on oceanic islands, based on published island species checklists, taxonomic studies, and phylogenetic analyses. We demonstrate that our understanding of evolutionary processes is being greatly enhanced through the use of genomic tools. However, to advance the field, it is critical to complement this information with a solid understanding of the ecological and behavioral traits of both extinct and extant avian island species.

The macroevolutionary dynamics of pharyngognathy in fishes fail to support the key innovation hypothesis

Key innovations, traits that provide species access to novel niches, are thought to be a major generator of biodiversity. One commonly cited example of key innovation is pharyngognathy, a set of modifications to the pharyngeal jaws found in some highly species-rich fish clades such as cichlids and wrasses. Here, using comparative phylogenomics and phylogenetic comparative methods, we investigate the genomic basis of pharyngognathy and the impact of this innovation on diversification. Whole genomes resolve the relationships of fish clades with this innovation and their close relatives, but high levels of topological discordance suggest the innovation may have evolved fewer times than previously thought. Closer examination of the topology of noncoding elements accelerated in clades with the pharyngognathy innovation reveals hidden patterns of shared ancestry across putatively independent transitions to pharyngognathy. When our updated phylogenomic relationships are used alongside large-scale phylogenetic and ecological datasets, we find no evidence pharyngognathy consistently modifies the macroevolutionary landscape of trophic ecology nor does it increase diversification. Our results highlight the necessity of incorporating genomic information in studies of key innovation.

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.

Effects of Testosterone on Gene Expression Are Concordant between Sexes but Divergent across Species of Sceloporus Lizards

AbstractHormones mediate sexual dimorphism by regulating sex-specific patterns of gene expression, but it is unclear how much of this regulation involves sex-specific hormone levels versus sex-specific transcriptomic responses to the same hormonal signal. Moreover, transcriptomic responses to hormones can evolve, but the extent to which hormonal pleiotropy in gene regulation is conserved across closely related species is not well understood. We addressed these issues by elevating testosterone levels in juvenile females and males of three Sceloporus lizard species before sexual divergence in circulating testosterone and then characterizing transcriptomic responses in the liver. In each species, more genes were responsive to testosterone in males than in females, suggesting that early developmental processes prime sex-specific transcriptomic responses to testosterone later in life. However, overall transcriptomic responses to testosterone were concordant between sexes, with no genes exhibiting sex-by-treatment interactions. By contrast, hundreds of genes exhibited species-by-treatment interactions, particularly when comparing distantly related species with different patterns of sexual dimorphism, suggesting evolutionary lability in gene regulation by testosterone. Collectively, our results indicate that early organizational effects may lead to sex-specific differences in the magnitude, but not the direction, of transcriptomic responses to testosterone and that the hormone-genome interface accrues regulatory changes over evolutionary time.

Insights into avian molecular cytogenetics-with reptilian comparisons

In last 100 years or so, much information has been accumulated on avian karyology, genetics, physiology, biochemistry and evolution. The chicken genome project generated genomic resources used in comparative studies, elucidating fundamental evolutionary processes, much of it funded by the economic importance of domestic fowl (which are also excellent model species in many areas). Studying karyotypes and whole genome sequences revealed population processes, evolutionary biology, and genome function, uncovering the role of repetitive sequences, transposable elements and gene family expansion. Knowledge of the function of many genes and non-expressed or identified regulatory components is however still lacking. Birds (Aves) are diverse, have striking adaptations for flight, migration and survival and inhabit all continents most islands. They also have a unique karyotype with ~ 10 macrochromosomes and ~ 30 microchromosomes that are smaller than other reptiles. Classified into Palaeognathae and Neognathae they are evolutionarily close, and a subset of reptiles. Here we overview avian molecular cytogenetics with reptilian comparisons, shedding light on their karyotypes and genome structure features. We consider avian evolution, then avian (followed by reptilian) karyotypes and genomic features. We consider synteny disruptions, centromere repositioning, and repetitive elements before turning to comparative avian and reptilian genomics. In this context, we review comparative cytogenetics and genome mapping in birds as well as Z- and W-chromosomes and sex determination. Finally, we give examples of pivotal research areas in avian and reptilian cytogenomics, particularly physical mapping and map integration of sex chromosomal genes, comparative genomics of chicken, turkey and zebra finch, California condor cytogenomics as well as some peculiar cytogenetic and evolutionary examples. We conclude that comparative molecular studies and improving resources continually contribute to new approaches in population biology, developmental biology, physiology, disease ecology, systematics, evolution and phylogenetic systematics orientation. This also produces genetic mapping information for chromosomes active in rearrangements during the course of evolution. Further insights into mutation, selection and adaptation of vertebrate genomes will benefit from these studies including physical and online resources for the further elaboration of comparative genomics approaches for many fundamental biological questions.

The case against simplistic genetic explanations of evolution

Humans are curious to understand the causes of traits that distinguish us from other animals and that distinguish vastly different species from one another. We also have a proclivity for simple stories and sometimes tend toward seeking and accepting simple genetic explanations for large evolutionary shifts, even to a single gene. Here, I reveal how a biased expectation of mechanistic simplicity threads through the long history of evolutionary and developmental genetics. I argue, however, that expecting a simple mechanism threatens a deeper understanding of evolution, and I define the limitations for interpreting experimental evidence in evolutionary developmental genetics.

Ancient developmental genes underlie evolutionary novelties in walking fish

A critical question in biology is how new traits evolve, but studying this in wild animals remains challenging. Here, we probe the genetic basis of trait gain in sea robin fish, which have evolved specialized leg-like appendages for locomotion and digging along the ocean floor. We use genome sequencing, transcriptional profiling, and interspecific hybrid analysis to explore the molecular and developmental basis of leg formation. We identified the ancient, conserved transcription factor tbx3a as a major determinant of sensory leg development. Genome editing confirms that tbx3a is required for normal leg formation in sea robins, and for formation of enlarged central nervous system lobes, sensory papillae, and adult digging behavior. Our study establishes sea robins as a model organism for studying the evolution of major trait gain and illustrates how ancient developmental control genes can underlie novel organ formation.

Evolution: Sea robins get a leg up

Sea robins, fish with legs, walk on the ocean bottom. They have evolved taste receptors on their legs that direct digging to access prey. Examining these structures and behaviors advances our understanding of the origin of novel phenotypes.

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.

Immobilization secondary to cell death of muscle precursors with a dual transcriptional signature contributes to the emu wing skeletal pattern

Limb reduction has occurred multiple times in tetrapod history. Among ratites, wing reductions range from mild vestigialization to complete loss, with emus (Dromaius novaehollandiae) serving as a model for studying the genetic mechanisms behind limb reduction. Here, we explore the developmental mechanisms underlying wing reduction in emu. Our analyses reveal that immobilization resulting from the absence of distal muscles contributes to skeletal shortening, fusion and left-right intraindividual variation. Expression analysis and single cell-RNA sequencing identify muscle progenitors displaying a dual lateral plate mesodermal and myogenic signature. These cells aggregate at the proximal region of wing buds and undergo cell death. We propose that this cell death, linked to the lack of distal muscle masses, underlines the morphological features and variability in skeletal elements due to reduced mechanical loading. Our results demonstrate that differential mobility during embryonic development may drive morphological diversification in vestigial structures.

Evolutionary changes of noncoding elements associated with transition of sexual mode in Caenorhabditis nematodes

The transition of the sexual mode occurs widely in animal evolution. In Caenorhabditis nematodes, androdioecy, a sexual polymorphism composed of males and hermaphrodites having the ability to self-fertilize, has evolved independently multiple times. While the modification of noncoding regulatory elements likely contributed to the evolution of hermaphroditism, little is known about these changes. Here, we conducted a genome-wide analysis of conserved noncoding elements (CNEs) focusing on the evolution of hermaphroditism in Caenorhabditis nematodes. We found that, in androdioecious nematodes, mutations rapidly accumulated in CNEs' neighboring genes associated with sexual traits. Expression analysis indicate that the identified CNEs are involved in spermatogenesis in hermaphrodites and associated with the transition of gene expression from dioecious to androdioecious nematodes. Last, genome editing of a CNE neighboring laf-1 resulted in a change in its expression in the gonadal region undergoing spermatogenesis. Our bioinformatic and experimental analyses highlight the importance of CNEs in gene regulation associated with the development of hermaphrodites.

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.

Alternative splicing in parallel evolution and the evolutionary potential in sticklebacks

Repeatability of adaptation to similar environments provides opportunity to evaluate the predictability of natural selection. While many studies have investigated gene expression differences between populations adapted to contrasting environments, the role of post-transcriptional processes such as alternative splicing has rarely been evaluated in the context of parallel adaptation. To address the aforementioned knowledge gap, we reanalysed transcriptomic data from three pairs of threespine stickleback (Gasterosteus aculeatus) ecotypes adapted to marine or freshwater environment. First, we identified genes with repeated expression or splicing divergence across ecotype pairs, and compared the genetic architecture and biological processes between parallelly expressed and parallelly spliced loci. Second, we analysed the extent to which parallel adaptation was reflected at gene expression and alternative splicing levels. Finally, we tested how the two axes of transcriptional variation differed in their potential for evolutionary change. Although both repeated differential splicing and differential expression across ecotype pairs showed tendency for parallel divergence, the degree of parallelism was lower for splicing than expression. Furthermore, parallel divergences in splicing and expression were likely to be associated with distinct cis-regulatory genetic variants and functionally unique set of genes. Finally, we found that parallelly spliced genes showed higher nucleotide diversity than parallelly expressed genes, indicating splicing is less susceptible to genetic variation erosion during parallel adaptation. Our results provide novel insight into the role of splicing in parallel adaptation, and underscore the contribution of splicing to the evolutionary potential of wild populations under environmental change.

Parallel Evolution at the Regulatory Base-Pair Level Contributes to Mammalian Interspecific Differences in Polygenic Traits

Parallel evolution occurs when distinct lineages with similar ancestral states converge on a new phenotype. Parallel evolution has been well documented at the organ, gene pathway, and amino acid sequence level but in theory, it can also occur at individual nucleotides within noncoding regions. To examine the role of parallel evolution in shaping the biology of mammalian complex traits, we used data on single-nucleotide polymorphisms (SNPs) influencing human intraspecific variation to predict trait values in other species for 11 complex traits. We found that the alleles at SNP positions associated with human intraspecific height and red blood cell (RBC) count variation are associated with interspecific variation in the corresponding traits across mammals. These associations hold for deeper branches of mammalian evolution as well as between strains of collaborative cross mice. While variation in RBC count between primates uses both ancient and more recently evolved genomic regions, we found that only primate-specific elements were correlated with primate body size. We show that the SNP positions driving these signals are flanked by conserved sequences, maintain synteny with target genes, and overlap transcription factor binding sites. This work highlights the potential of conserved but tunable regulatory elements to be reused in parallel to facilitate evolutionary adaptation in mammals.

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.

Genome-Wide Gene Birth-Death Dynamics Are Associated with Diet Breadth Variation in Lepidoptera

Comparative analyses of gene birth-death dynamics have the potential to reveal gene families that played an important role in the evolution of morphological, behavioral, or physiological variation. Here, we used whole genomes of 30 species of butterflies and moths to identify gene birth-death dynamics among the Lepidoptera that are associated with specialist or generalist feeding strategies. Our work advances this field using a uniform set of annotated proteins for all genomes, investigating associations while correcting for phylogeny, and assessing all gene families rather than a priori subsets. We discovered that the sizes of several important gene families (e.g. those associated with pesticide resistance, xenobiotic detoxification, and/or protein digestion) are significantly correlated with diet breadth. We also found 22 gene families showing significant shifts in gene birth-death dynamics at the butterfly (Papilionoidea) crown node, the most notable of which was a family of pheromone receptors that underwent a contraction potentially linked with a shift to visual-based mate recognition. Our findings highlight the importance of uniform annotations, phylogenetic corrections, and unbiased gene family analyses in generating a list of candidate genes that warrant further exploration.

Paleoneurology of stem palaeognaths clarifies the plesiomorphic condition of the crown bird central nervous system

Lithornithidae, an assemblage of volant Palaeogene fossil birds, provide our clearest insights into the early evolutionary history of Palaeognathae, the clade that today includes the flightless ratites and volant tinamous. The neotype specimen of Lithornis vulturinus, from the early Eocene (approximately 53 million years ago) of Europe, includes a partial neurocranium that has never been thoroughly investigated. Here, we describe these cranial remains including the nearly complete digital endocasts of the brain and bony labyrinth. The telencephalon of Lithornis is expanded and its optic lobes are ventrally shifted, as is typical for crown birds. The foramen magnum is positioned caudally, rather than flexed ventrally as in some crown birds, with the optic lobes, cerebellum, and foramen magnum shifted further ventrally. The overall brain shape is similar to that of tinamous, the only extant clade of flying palaeognaths, suggesting that several aspects of tinamou neuroanatomy may have been evolutionarily conserved since at least the early Cenozoic. The estimated ratio of the optic lobe's surface area relative to the total brain suggests a diurnal ecology. Lithornis may provide the clearest insights to date into the neuroanatomy of the ancestral crown bird, combining an ancestrally unflexed brain with a caudally oriented connection with the spinal cord, a moderately enlarged telencephalon, and ventrally shifted, enlarged optic lobes.

Current status of community resources and priorities for weed genomics research

Weeds are attractive models for basic and applied research due to their impacts on agricultural systems and capacity to swiftly adapt in response to anthropogenic selection pressures. Currently, a lack of genomic information precludes research to elucidate the genetic basis of rapid adaptation for important traits like herbicide resistance and stress tolerance and the effect of evolutionary mechanisms on wild populations. The International Weed Genomics Consortium is a collaborative group of scientists focused on developing genomic resources to impact research into sustainable, effective weed control methods and to provide insights about stress tolerance and adaptation to assist crop breeding.

A nuclear genome assembly of an extinct flightless bird, the little bush moa

We present a draft genome of the little bush moa (Anomalopteryx didiformis)-one of approximately nine species of extinct flightless birds from Aotearoa, New Zealand-using ancient DNA recovered from a fossil bone from the South Island. We recover a complete mitochondrial genome at 249.9× depth of coverage and almost 900 megabases of a male moa nuclear genome at ~4 to 5× coverage, with sequence contiguity sufficient to identify more than 85% of avian universal single-copy orthologs. We describe a diverse landscape of transposable elements and satellite repeats, estimate a long-term effective population size of ~240,000, identify a diverse suite of olfactory receptor genes and an opsin repertoire with sensitivity in the ultraviolet range, show that the wingless moa phenotype is likely not attributable to gene loss or pseudogenization, and identify potential function-altering coding sequence variants in moa that could be synthesized for future functional assays. This genomic resource should support further studies of avian evolution and morphological divergence.

Convergent accelerated evolution of mammal-specific conserved non-coding elements in hibernators

Mammals maintain their body temperature, yet hibernators can temporarily lower their metabolic rate as an energy-saving strategy. It has been proposed that hibernators evolved independently from homeotherms, and it is possible that the convergent evolution of hibernation involved common genomic changes among hibernator-lineages. Since hibernation is a seasonal trait, the evolution of gene regulatory regions in response to changes in season may have been important for the acquisition of hibernation traits. High-frequency accumulation of mutations in conserved non-coding elements (CNEs) could, in principle, alter the expression of neighboring genes and thereby contribute to the acquisition of new traits. To address this possibility, we performed a comparative genomic analysis of mammals to identify accelerated CNEs commonly associated with hibernation. We found that accelerated CNEs are common to hibernator-lineages and could be involved with hibernation. We also found that common factors of genes that located near accelerated CNEs and are differentially expressed between normal and hibernation periods related to gene regulation and cell-fate determination. It suggests that the molecular mechanisms controlling hibernation have undergone convergent evolution. These results help broaden our understanding of the genetic adaptations that facilitated hibernation in mammals and may offer insights pertaining to stress responses and energy conservation.

Conserved regulatory switches for the transition from natal down to juvenile feather in birds

The transition from natal downs for heat conservation to juvenile feathers for simple flight is a remarkable environmental adaptation process in avian evolution. However, the underlying epigenetic mechanism for this primary feather transition is mostly unknown. Here we conducted time-ordered gene co-expression network construction, epigenetic analysis, and functional perturbations in developing feather follicles to elucidate four downy-juvenile feather transition events. We report that extracellular matrix reorganization leads to peripheral pulp formation, which mediates epithelial-mesenchymal interactions for branching morphogenesis. α-SMA (ACTA2) compartmentalizes dermal papilla stem cells for feather renewal cycling. LEF1 works as a key hub of Wnt signaling to build rachis and converts radial downy to bilateral symmetry. Novel usage of scale keratins strengthens feather sheath with SOX14 as the epigenetic regulator. We show that this primary feather transition is largely conserved in chicken (precocial) and zebra finch (altricial) and discuss the possibility that this evolutionary adaptation process started in feathered dinosaurs.

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.

Parthenogenetic Stick Insects Exhibit Signatures of Preservation in the Molecular Architecture of Male Reproduction

After the loss of a trait, theory predicts that the molecular machinery underlying its phenotypic expression should decay. Yet, empirical evidence is contrasting. Here, we test the hypotheses that (i) the molecular ground plan of a lost trait could persist due to pleiotropic effects on other traits and (ii) that gene co-expression network architecture could constrain individual gene expression. Our testing ground has been the Bacillus stick insect species complex, which contains close relatives that are either bisexual or parthenogenetic. After the identification of genes expressed in male reproductive tissues in a bisexual species, we investigated their gene co-expression network structure in two parthenogenetic species. We found that gene co-expression within the male gonads was partially preserved in parthenogens. Furthermore, parthenogens did not show relaxed selection on genes upregulated in male gonads in the bisexual species. As these genes were mostly expressed in female gonads, this preservation could be driven by pleiotropic interactions and an ongoing role in female reproduction. Connectivity within the network also played a key role, with highly connected-and more pleiotropic-genes within male gonad also having a gonad-biased expression in parthenogens. Our findings provide novel insight into the mechanisms which could underlie the production of rare males in parthenogenetic lineages; more generally, they provide an example of the cryptic persistence of a lost trait molecular architecture, driven by gene pleiotropy on other traits and within their co-expression network.

Emx2 underlies the development and evolution of marsupial gliding membranes

Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium-an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

Evolution of bird sex chromosomes: a cytogenomic approach in Palaeognathae species

BACKGROUND

Different patterns of sex chromosome differentiation are seen in Palaeognathae birds, a lineage that includes the ratites (Struthioniformes, Rheiformes, Apterygiformes, Casuariiformes, and the sister group Tinamiformes). While some Tinamiform species have well-differentiated W chromosomes, both Z and W of all the flightless ratites are still morphologically undifferentiated. Here, we conducted a comprehensive analysis of the ZW differentiation in birds using a combination of cytogenetic, genomic, and bioinformatic approaches. The whole set of satDNAs from the emu (Dromaius novaehollandiae) was described and characterized. Furthermore, we examined the in situ locations of these satDNAs alongside several microsatellite repeats and carried out Comparative Genomic Hybridizations in two related species: the greater rhea (Rhea americana) and the tataupa tinamou (Crypturellus tataupa).

RESULTS

From the 24 satDNA families identified (which represent the greatest diversity of satDNAs ever uncovered in any bird species), only three of them were found to accumulate on the emu's sex chromosomes, with no discernible accumulation observed on the W chromosome. The W chromosomes of both the greater rhea and the emu did not exhibit a significant buildup of either C-positive heterochromatin or repetitive DNAs, indicating their large undifferentiation both at morphological and molecular levels. In contrast, the tataupa tinamou has a highly differentiated W chromosome that accumulates several DNA repeats.

CONCLUSION

The findings provide new information on the architecture of the avian genome and an inside look at the starting points of sex chromosome differentiation in birds.

The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

Ostrich (Struthio camelus) syrinx morphology and vocal repertoire across postnatal ontogeny and sex: Implications for understanding vocal evolution in birds

Vocal production in birds has been the target of considerable research that mostly has focused on phylogenetically well-nested songbirds. Anatomical descriptions and recordings of many non-songbirds have often only focused on a single ontogenetic stage or sex. While basic morphology of the vocal organ (syrinx) of ostrich (Palaeognathae, Struthio camelus) has been known since the 1800s, descriptions of its vocal repertoire and syrinx anatomy since then have been incomplete or inconsistent. New toolkits now enable detailed qualitative description of internal anatomy and meristic data and allow it to be compared to vocal production. Here we describe the anatomy of the syrinx in Struthio camelus for three post-hatching ontogenetic stages and both an adult male and female utilizing dissection and contrast enhanced X-ray computed tomography (diceCT). We find changes in ring geometry and spacing through ontogeny as well as lateral labia thickness. We document a small unpaired, midline, cartilaginous structure, a "pessuliform process" at the tracheobronchial juncture present throughout ontogeny and in both males and females. Investigation of the vocal repertoire of ostriches across ontogeny using a new dataset of 77 recordings led to identification of four vocalizations not previously reported in the literature, including the simultaneous production of a hiss and tonal. We find syrinx morphology largely consistent across ontogeny and in male and female adults. Both are capable of producing long duration tonal calls, but these may be more frequent in male birds. Closed-mouth boom calls remain unique to males. A detailed understanding of diversity in parts of early diverging clades is pivotal in attempting to estimate features of the ancestral syrinx in birds and how avian vocalization evolved.

Divergence time shapes gene reuse during repeated adaptation

When diverse lineages repeatedly adapt to similar environmental challenges, the extent to which the same genes are involved (gene reuse) varies across systems. We propose that divergence time among lineages is a key factor driving this variability: as lineages diverge, the extent of gene reuse should decrease due to reductions in allele sharing, functional differentiation among genes, and restructuring of genome architecture. Indeed, we show that many genomic studies of repeated adaptation find that more recently diverged lineages exhibit higher gene reuse during repeated adaptation, but the relationship becomes less clear at older divergence time scales. Thus, future research should explore the factors shaping gene reuse and their interplay across broad divergence time scales for a deeper understanding of evolutionary repeatability.

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.

Convergent genomic signatures associated with vertebrate viviparity

BACKGROUND

Viviparity-live birth-is a complex and innovative mode of reproduction that has evolved repeatedly across the vertebrate Tree of Life. Viviparous species exhibit remarkable levels of reproductive diversity, both in the amount of care provided by the parent during gestation, and the ways in which that care is delivered. The genetic basis of viviparity has garnered increasing interest over recent years; however, such studies are often undertaken on small evolutionary timelines, and thus are not able to address changes occurring on a broader scale. Using whole genome data, we investigated the molecular basis of this innovation across the diversity of vertebrates to answer a long held question in evolutionary biology: is the evolution of convergent traits driven by convergent genomic changes?

RESULTS

We reveal convergent changes in protein family sizes, protein-coding regions, introns, and untranslated regions (UTRs) in a number of distantly related viviparous lineages. Specifically, we identify 15 protein families showing evidence of contraction or expansion associated with viviparity. We additionally identify elevated substitution rates in both coding and noncoding sequences in several viviparous lineages. However, we did not find any convergent changes-be it at the nucleotide or protein level-common to all viviparous lineages.

CONCLUSIONS

Our results highlight the value of macroevolutionary comparative genomics in determining the genomic basis of complex evolutionary transitions. While we identify a number of convergent genomic changes that may be associated with the evolution of viviparity in vertebrates, there does not appear to be a convergent molecular signature shared by all viviparous vertebrates. Ultimately, our findings indicate that a complex trait such as viviparity likely evolves with changes occurring in multiple different pathways.

Comparative phylogenomics and phylotranscriptomics provide insights into the genetic complexity of nitrogen-fixing root-nodule symbiosis

Plant root-nodule symbiosis (RNS) with mutualistic nitrogen-fixing bacteria is restricted to a single clade of angiosperms, the Nitrogen-Fixing Nodulation Clade (NFNC), and is best understood in the legume family. Nodulating species share many commonalities, explained either by divergence from a common ancestor over 100 million years ago or by convergence following independent origins over that same time period. Regardless, comparative analyses of diverse nodulation syndromes can provide insights into constraints on nodulation-what must be acquired or cannot be lost for a functional symbiosis-and the latitude for variation in the symbiosis. However, much remains to be learned about nodulation, especially outside of legumes. Here, we employed a large-scale phylogenomic analysis across 88 species, complemented by 151 RNA-seq libraries, to elucidate the evolution of RNS. Our phylogenomic analyses further emphasize the uniqueness of the transcription factor NIN as a master regulator of nodulation and identify key mutations that affect its function across the NFNC. Comparative transcriptomic assessment revealed nodule-specific upregulated genes across diverse nodulating plants, while also identifying nodule-specific and nitrogen-response genes. Approximately 70% of symbiosis-related genes are highly conserved in the four representative species, whereas defense-related and host-range restriction genes tend to be lineage specific. Our study also identified over 900 000 conserved non-coding elements (CNEs), over 300 000 of which are unique to sampled NFNC species. NFNC-specific CNEs are enriched with the active H3K9ac mark and are correlated with accessible chromatin regions, thus representing a pool of candidate regulatory elements for genes involved in RNS. Collectively, our results provide novel insights into the evolution of nodulation and lay a foundation for engineering of RNS traits in agriculturally important crops.

Description of natal down of the ostrich (Struthio camelus) and comparison with common quail (Coturnix coturnix): Developmental and evolutionary implications

Natal down is a feather stage that differs in both form and function from the definitive feathers of adult birds. It has a simpler structure that has been speculated to be similar to the body coverings of non-avian dinosaurs. However, inference of the evolution of natal down has been limited by our understanding of its structural variation in extant birds. Most descriptive work has focused on neognathous birds, limiting our knowledge of the full diversity of feathers in extant taxa. Here, we describe the natal down of a post-hatch ostrich (Struthio camelus) and compare it to that of a post-hatch quail (Coturnix coturnix). We confirm the presence of featherless spaces (apteria) in S. camelus and the lack of barbules on the tips of natal down in both species. We also find differences between dorsal and ventral natal down structures, such as barbule density in S. camelus and the extent of the bare portion of the barb in both species. Surprisingly, we do not find that the neoptiles of either species follow the ideal morphologies for increasing insulation. Finally, we hypothesize that the different barb types present in S. camelus natal down result from a large addition of new barb ridges during development, which is not known except in feathers with a rachis. These results have implications for our understanding of how structure informs function and development in understudied feather types, such as those shared by non-avian dinosaurs.

Species differences in hormonally mediated gene expression underlie the evolutionary loss of sexually dimorphic coloration in Sceloporus lizards

Phenotypic sexual dimorphism often involves the hormonal regulation of sex-biased expression for underlying genes. However, it is generally unknown whether the evolution of hormonally mediated sexual dimorphism occurs through upstream changes in tissue sensitivity to hormone signals, downstream changes in responsiveness of target genes, or both. Here, we use comparative transcriptomics to explore these possibilities in 2 species of Sceloporus lizards exhibiting different patterns of sexual dichromatism. Sexually dimorphic S. undulatus develops blue and black ventral coloration in response to testosterone, while sexually monomorphic S. virgatus does not, despite exhibiting similar sex differences in circulating testosterone levels. We administered testosterone implants to juveniles of each species and used RNAseq to quantify gene expression in ventral skin. Transcriptome-wide responses to testosterone were stronger in S. undulatus than in S. virgatus, suggesting species differences in tissue sensitivity to this hormone signal. Species differences in the expression of genes for androgen metabolism and sex hormone-binding globulin were consistent with this idea, but expression of the androgen receptor gene was higher in S. virgatus, complicating this interpretation. Downstream of androgen signaling, we found clear species differences in hormonal responsiveness of genes related to melanin synthesis, which were upregulated by testosterone in S. undulatus, but not in S. virgatus. Collectively, our results indicate that hormonal regulation of melanin synthesis pathways contributes to the development of sexual dimorphism in S. undulatus, and that changes in the hormonal responsiveness of these genes in S. virgatus contribute to the evolutionary loss of ventral coloration.

3D atlas of tinamou (Neornithes: Tinamidae) pectoral morphology: Implications for reconstructing the ancestral neornithine flight apparatus

Palaeognathae, the extant avian clade comprising the flightless ratites and flight-capable tinamous (Tinamidae), is the sister group to all other living birds, and recent phylogenetic studies illustrate that tinamous are phylogenetically nested within a paraphyletic assemblage of ratites. As the only extant palaeognaths that have retained the ability to fly, tinamous may provide key information on the nature of the flight apparatus of ancestral crown palaeognaths-and, in turn, crown birds-as well as insight into convergent modifications to the wing apparatus among extant ratite lineages. To reveal new information about the musculoskeletal anatomy of tinamous and facilitate development of computational biomechanical models of tinamou wing function, we generated a three-dimensional musculoskeletal model of the flight apparatus of the extant Andean tinamou (Nothoprocta pentlandii) using diffusible iodine-based contrast-enhanced computed tomography (diceCT). Origins and insertions of the pectoral flight musculature of N. pentlandii are generally consistent with those of other extant volant birds specialized for burst flight, and the entire suite of presumed ancestral neornithine flight muscles are present in N. pentlandii with the exception of the biceps slip. The pectoralis and supracoracoideus muscles are robust, similar to the condition in other extant burst-flying birds such as many extant Galliformes. Contrary to the condition in most extant Neognathae (the sister clade to Palaeognathae), the insertion of the pronator superficialis has a greater distal extent than the pronator profundus, although most other anatomical observations are broadly consistent with the conditions observed in extant neognaths. This work will help form a basis for future comparative studies of the avian musculoskeletal system, with implications for reconstructing the flight apparatus of ancestral crown birds and clarifying musculoskeletal modifications underlying the convergent origins of ratite flightlessness.

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.

The genetic basis of novel trait gain in walking fish

A major goal in biology is to understand how organisms evolve novel traits. Multiple studies have identified genes contributing to regressive evolution, the loss of structures that existed in a recent ancestor. However, fewer examples exist for genes underlying constructive evolution, the gain of novel structures and capabilities in lineages that previously lacked them. Sea robins are fish that have evolved enlarged pectoral fins, six mobile locomotory fin rays (legs) and six novel macroscopic lobes in the central nervous system (CNS) that innervate the corresponding legs. Here, we establish successful husbandry and use a combination of transcriptomics, CRISPR-Cas9 editing, and behavioral assays to identify key transcription factors that are required for leg formation and function in sea robins. We also generate hybrids between two sea robin species with distinct leg morphologies and use allele-specific expression analysis and gene editing to explore the genetic basis of species-specific trait diversity, including a novel sensory gain of function. Collectively, our study establishes sea robins as a new model for studying the genetic basis of novel organ formation, and demonstrates a crucial role for the conserved limb gene tbx3a in the evolution of chemosensory legs in walking fish.

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.