Convergent and lineage-specific genomic differences in limb regulatory elements in limbless reptile lineages

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.

Blooming balloons: Searching for mechanisms of the inflated calyx

Studying morphological novelties offers special insights into developmental biology and evolution. The inflated calyx syndrome (ICS) is a largely unrecognized but fascinating feature of flower development, where sepals form balloon-like husks that encapsulate fruits. Despite its independent emergence in many lineages of flowering plants, the genetic and molecular mechanisms of ICS remain unknown. Early studies in the Solanaceae genus Physalis put forth key roles of MADS-box genes in ICS. However, recent work suggests these classical floral identity transcription factors were false leads. With newfound capabilities that allow rapid development of genetic systems through genomics and genome editing, Physalis has re-emerged as the most tractable model species for dissecting ICS. This review revisits current understanding of ICS and highlights how recent advancements enable a reset in the search for genetic and molecular mechanisms using unbiased, systematic approaches.

Conserved Noncoding Cis-Elements Associated with Hibernation Modulate Metabolic and Behavioral Adaptations in Mice

Our study elucidates functional roles for conserved cis-elements associated with the evolution of mammalian hibernation. Genomic analyses found topologically associated domains (TADs) that disproportionately accumulated convergent genomic changes in hibernators, including the TAD for the Fat Mass & Obesity (Fto) locus. Some hibernation-linked cis-elements in this TAD form regulatory contacts with multiple neighboring genes. Knockout mice for these cis-elements exhibit Fto, Irx3, and Irx5 gene expression changes, impacting hundreds of genes downstream. Profiles of pre-torpor, torpor, and post-torpor phenotypes found distinct roles for each cis-element in metabolic control, while a high caloric diet uncovered different obesogenic effects. One cis-element promoting a lean phenotype influences foraging behaviors throughout life, affecting specific behavioral sequences. Thus, convergent evolution in hibernators pinpoints functional genetic mechanisms of mammalian metabolic control.

Gene loss and relaxed selection of plaat1 in vertebrates adapted to low-light environments

Gene loss is an important mechanism for evolution in low-light or cave environments where visual adaptations often involve a reduction or loss of eyesight. The plaat gene family encodes phospholipases essential for the degradation of organelles in the lens of the eye. These phospholipases translocate to damaged organelle membranes, inducing them to rupture. This rupture is required for lens transparency and is essential for developing a functioning eye. Plaat3 is thought to be responsible for this role in mammals, while plaat1 is thought to be responsible in other vertebrates. We used a macroevolutionary approach and comparative genomics to examine the origin, loss, synteny and selection of plaat1 across bony fishes and tetrapods. We showed that plaat1 (probably ancestral to all bony fish + tetrapods) has been lost in squamates and is significantly degraded in lineages of low-visual-acuity and blind mammals and fishes. Our findings suggest that plaat1 is important for visual acuity across bony vertebrates, and that its loss through relaxed selection and pseudogenization may have played a role in the repeated evolution of visual systems in low-light environments. Our study sheds light on the importance of gene-loss in trait evolution and provides insights into the mechanisms underlying visual acuity in low-light environments.

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

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Differential Conservation and Loss of CR1 Retrotransposons in Squamates Reveals Lineage-Specific Genome Dynamics across Reptiles

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniotic vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; ~11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific dynamics have evolved over the course of squamate evolution to constrain genome size across the order. Thus, squamates may represent a prime model for investigations into TE diversity and evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the CR1 retrotransposon, a TE family found in most tetrapod genomes. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds, and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

Gene loss and relaxed selection of plaat1 in vertebrates adapted to low-light environments

Gene loss is an important mechanism for evolution in low-light or cave environments where visual adaptations often involve a reduction or loss of eyesight. The plaat gene family are phospholipases essential for the degradation of organelles in the lens of the eye. They translocate to damaged organelle membranes, inducing them to rupture. This rupture is required for lens transparency and is essential for developing a functioning eye. Plaat3 is thought to be responsible for this role in mammals, while plaat1 is thought to be responsible in other vertebrates. We used a macroevolutionary approach and comparative genomics to examine the origin, loss, synteny, and selection of plaat1 across bony fishes and tetrapods. We show that plaat1 (likely ancestral to all bony fish + tetrapods) has been lost in squamates and is significantly degraded in lineages of low-visual acuity and blind mammals and fish. Our findings suggest that plaat1 is important for visual acuity across bony vertebrates, and that its loss through relaxed selection and pseudogenization may have played a role in the repeated evolution of visual systems in low-light-environments. Our study sheds light on the importance of gene-loss in trait evolution and provides insights into the mechanisms underlying visual acuity in low-light environments.

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.

Developmental regulation of conserved non-coding element evolution provides insights into limb loss in squamates

Limb loss shows recurrent phenotypic evolution across squamate lineages. Here, based on three de novo-assembled genomes of limbless lizards from different lineages, we showed that divergence of conserved non-coding elements (CNEs) played an important role in limb development. These CNEs were associated with genes required for limb initiation and outgrowth, and with regulatory signals in the early stage of limb development. Importantly, we identified the extensive existence of insertions and deletions (InDels) in the CNEs, with the numbers ranging from 111 to 756. Most of these CNEs with InDels were lineage-specific in the limbless squamates. Nearby genes of these InDel CNEs were important to early limb formation, such as Tbx4, Fgf10, and Gli3. Based on functional experiments, we found that nucleotide mutations and InDels both affected the regulatory function of the CNEs. Our study provides molecular evidence underlying limb loss in squamate reptiles from a developmental perspective and sheds light on the importance of regulatory element InDels in phenotypic evolution.

Whole-genome Comparisons Identify Repeated Regulatory Changes Underlying Convergent Appendage Evolution in Diverse Fish Lineages

Fins are major functional appendages of fish that have been repeatedly modified in different lineages. To search for genomic changes underlying natural fin diversity, we compared the genomes of 36 percomorph fish species that span over 100 million years of evolution and either have complete or reduced pelvic and caudal fins. We identify 1,614 genomic regions that are well-conserved in fin-complete species but missing from multiple fin-reduced lineages. Recurrent deletions of conserved sequences in wild fin-reduced species are enriched for functions related to appendage development, suggesting that convergent fin reduction at the organismal level is associated with repeated genomic deletions near fin-appendage development genes. We used sequencing and functional enhancer assays to confirm that PelA, a Pitx1 enhancer previously linked to recurrent pelvic loss in sticklebacks, has also been independently deleted and may have contributed to the fin morphology in distantly related pelvic-reduced species. We also identify a novel enhancer that is conserved in the majority of percomorphs, drives caudal fin expression in transgenic stickleback, is missing in tetraodontiform, syngnathid, and synbranchid species with caudal fin reduction, and alters caudal fin development when targeted by genome editing. Our study illustrates a broadly applicable strategy for mapping phenotypes to genotypes across a tree of vertebrate species and highlights notable new examples of regulatory genomic hotspots that have been used to evolve recurrent phenotypes across 100 million years of fish evolution.

The genetic architecture and evolution of the human skeletal form

The human skeletal form underlies bipedalism, but the genetic basis of skeletal proportions (SPs) is not well characterized. We applied deep-learning models to 31,221 x-rays from the UK Biobank to extract a comprehensive set of SPs, which were associated with 145 independent loci genome-wide. Structural equation modeling suggested that limb proportions exhibited strong genetic sharing but were independent of width and torso proportions. Polygenic score analysis identified specific associations between osteoarthritis and hip and knee SPs. In contrast to other traits, SP loci were enriched in human accelerated regions and in regulatory elements of genes that are differentially expressed between humans and great apes. Combined, our work identifies specific genetic variants that affect the skeletal form and ties a major evolutionary facet of human anatomical change to pathogenesis.

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.

The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order

Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.

Integrating gene annotation with orthology inference at scale

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

Whole-genome comparisons identify repeated regulatory changes underlying convergent appendage evolution in diverse fish lineages

Fins are major functional appendages of fish that have been repeatedly modified in different lineages. To search for genomic changes underlying natural fin diversity, we compared the genomes of 36 wild fish species that either have complete or reduced pelvic and caudal fins. We identify 1,614 genomic regions that are well-conserved in fin-complete species but missing from multiple fin-reduced lineages. Recurrent deletions of conserved sequences (CONDELs) in wild fin-reduced species are enriched for functions related to appendage development, suggesting that convergent fin reduction at the organismal level is associated with repeated genomic deletions near fin-appendage development genes. We used sequencing and functional enhancer assays to confirm that PelA , a Pitx1 enhancer previously linked to recurrent pelvic loss in sticklebacks, has also been independently deleted and may have contributed to the fin morphology in distantly related pelvic-reduced species. We also identify a novel enhancer that is conserved in the majority of percomorphs, drives caudal fin expression in transgenic stickleback, is missing in tetraodontiform, s yngnathid, and synbranchid species with caudal fin reduction, and which alters caudal fin development when targeted by genome editing. Our study illustrates a general strategy for mapping phenotypes to genotypes across a tree of vertebrate species, and highlights notable new examples of regulatory genomic hotspots that have been used to evolve recurrent phenotypes during 100 million years of fish evolution.

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.

The genetic architecture of the human skeletal form

The human skeletal form underlies our ability to walk on two legs, but unlike standing height, the genetic basis of limb lengths and skeletal proportions is less well understood. Here we applied a deep learning model to 31,221 whole body dual-energy X-ray absorptiometry (DXA) images from the UK Biobank (UKB) to extract 23 different image-derived phenotypes (IDPs) that include all long bone lengths as well as hip and shoulder width, which we analyzed while controlling for height. All skeletal proportions are highly heritable (∼40-50%), and genome-wide association studies (GWAS) of these traits identified 179 independent loci, of which 102 loci were not associated with height. These loci are enriched in genes regulating skeletal development as well as associated with rare human skeletal diseases and abnormal mouse skeletal phenotypes. Genetic correlation and genomic structural equation modeling indicated that limb proportions exhibited strong genetic sharing but were genetically independent of width and torso proportions. Phenotypic and polygenic risk score analyses identified specific associations between osteoarthritis (OA) of the hip and knee, the leading causes of adult disability in the United States, and skeletal proportions of the corresponding regions. We also found genomic evidence of evolutionary change in arm-to-leg and hip-width proportions in humans consistent with striking anatomical changes in these skeletal proportions in the hominin fossil record. In contrast to cardiovascular, auto-immune, metabolic, and other categories of traits, loci associated with these skeletal proportions are significantly enriched in human accelerated regions (HARs), and regulatory elements of genes differentially expressed through development between humans and the great apes. Taken together, our work validates the use of deep learning models on DXA images to identify novel and specific genetic variants affecting the human skeletal form and ties a major evolutionary facet of human anatomical change to pathogenesis.

Chromosome-level genome assembly of the Muscovy duck provides insight into fatty liver susceptibility

The Muscovy duck (Cairina moschata) is an economically important poultry species, which is susceptible to fatty liver. Thus, the Muscovy duck may serve as an excellent candidate animal model of non-alcoholic fatty liver disease. However, the mechanisms underlying fatty liver development in this species are poorly understood. In this study, we report a chromosome-level genome assembly of the Muscovy duck, with a contig N50 of 11.8 Mb and scaffold N50 of 83.16 Mb. The susceptibility of Muscovy duck to fatty liver was mainly attributed to weak lipid catabolism capabilities (fatty acid β-oxidation and lipolysis). Furthermore, conserved noncoding elements (CNEs) showing accelerated evolution contributed to fatty liver formation by down-regulating the expression of genes involved in hepatic lipid catabolism. We propose that the susceptibility of Muscovy duck to fatty liver is an evolutionary by-product. In conclusion, this study revealed the potential mechanisms underlying the susceptibility of Muscovy duck to fatty liver.

Attempting genetic inference from directional asymmetry during convergent hindlimb reduction in squamates

Loss and reduction in paired appendages are common in vertebrate evolution. How often does such convergent evolution depend on similar developmental and genetic pathways? For example, many populations of the threespine stickleback and ninespine stickleback (Gasterosteidae) have independently evolved pelvic reduction, usually based on independent mutations that caused reduced Pitx1 expression. Reduced Pitx1 expression has also been implicated in pelvic reduction in manatees. Thus, hindlimb reduction stemming from reduced Pitx1 expression has arisen independently in groups that diverged tens to hundreds of millions of years ago, suggesting a potential for repeated use of Pitx1 across vertebrates. Notably, hindlimb reduction based on the reduction in Pitx1 expression produces left-larger directional asymmetry in the vestiges. We used this phenotypic signature as a genetic proxy, testing for hindlimb directional asymmetry in six genera of squamate reptiles that independently evolved hindlimb reduction and for which genetic and developmental tools are not yet developed: Agamodon anguliceps, Bachia intermedia, Chalcides sepsoides, Indotyphlops braminus, Ophisaurus attenuatuas and O. ventralis, and Teius teyou. Significant asymmetry occurred in one taxon, Chalcides sepsoides, whose left-side pelvis and femur vestiges were 18% and 64% larger than right-side vestiges, respectively, suggesting modification in Pitx1 expression in that species. However, there was either right-larger asymmetry or no directional asymmetry in the other five taxa, suggesting multiple developmental genetic pathways to hindlimb reduction in squamates and the vertebrates more generally.