Comparative Genomics Reveals Accelerated Evolution in Conserved Pathways during the Diversification of Anole Lizards

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

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

Brief Communication: The Predictable Network Topology of Evolutionary Genomic Constraint

Large-scale comparative genomics studies offer valuable resources for understanding both functional and evolutionary rate constraints. It is suggested that constraint aligns with the topology of genomic networks, increasing toward the center, with intermediate nodes combining relaxed constraint with higher contributions to the phenotype due to pleiotropy. However, this pattern has yet to be demonstrated in vertebrates. This study shows that constraint intensifies toward the network's center in placental mammals. Genes with rate changes associated with emergence of hibernation cluster mostly toward intermediate positions, with higher constraint in faster-evolving genes, which is indicative of a "sweet spot" for adaptation. If this trend holds universally, network node metrics could predict high-constraint regions even in clades lacking empirical constraint data.

A high-quality genome for the slender anole (Anolis apletophallus): an emerging model for field studies of tropical ecology and evolution

The slender anole, Anolis apletophallus, is a small arboreal lizard of the rainforest understory of central and eastern Panama. This species has been the subject of numerous ecological and evolutionary studies over the past 60 years as a result of attributes that make it especially amenable to field and laboratory science. Slender anoles are highly abundant, short-lived (nearly 100% annual turnover), easy to manipulate in both the lab and field, and are ubiquitous in the forests surrounding the Smithsonian Tropical Research Institute in Panama, where researchers have access to high-quality laboratory facilities. Here, we present a high-quality genome for the slender anole, which is an important new resource for studying this model species. We assembled and annotated the slender anole genome by combining 3 technologies: Oxford Nanopore, 10× Genomics Linked-Reads, and Dovetail Omni-C. We compared this genome with the recently published brown anole (Anolis sagrei) and the canonical green anole (Anolis carolinensis) genomes. Our genome is the first assembled for an Anolis lizard from mainland Central or South America, the regions that host the majority of diversity in the genus. This new reference genome is one of the most complete genomes of any anole assembled to date and should facilitate deeper studies of slender anole evolution, as well as broader scale comparative genomic studies of both mainland and island species. In turn, such studies will further our understanding of the well-known adaptive radiation of Anolis lizards.

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.

Insights into the Genomics of Clownfish Adaptive Radiation: The Genomic Substrate of the Diversification

Clownfishes are an iconic group of coral reef fishes that evolved a mutualistic interaction with sea anemones, which triggered the rapid diversification of the group. Following the emergence of this mutualism, clownfishes diversified into different ecological niches and developed convergent phenotypes associated with their host use. The genetic basis of the initial acquisition of the mutualism with host anemones has been described, but the genomic architecture underlying clownfish diversification once the mutualism was established and the extent to which clownfish phenotypic convergence originated through shared genetic mechanisms are still unknown. Here, we investigated these questions by performing comparative genomic analyses on the available genomic data of five pairs of closely related but ecologically divergent clownfish species. We found that clownfish diversification was characterized by bursts of transposable elements, an overall accelerated coding evolution, incomplete lineage sorting, and ancestral hybridization events. Additionally, we detected a signature of positive selection in 5.4% of the clownfish genes. Among them, five presented functions associated with social behavior and ecology, and they represent candidate genes involved in the evolution of the size-based hierarchical social structure so particular to clownfishes. Finally, we found genes with patterns of either relaxation or intensification of purifying selection and signals of positive selection linked with clownfish ecological divergence, suggesting some level of parallel evolution during the diversification of the group. Altogether, this work provides the first insights into the genomic substrate of clownfish adaptive radiation and integrates the growing collection of studies investigating the genomic mechanisms governing species diversification.

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.

Draft genome of six Cuban Anolis lizards and insights into genetic changes during their diversification

Background

Detecting genomic variants and their accumulation processes during species diversification and adaptive radiation is important for understanding the molecular and genetic basis of evolution. Anolis lizards in the West Indies are good models for studying evolutionary mechanisms because of the repeated evolution of their morphology and the ecology. We performed de novo genome assembly of six Cuban Anolis lizards with different ecomorphs and thermal habitats (Anolis isolepis, Anolis allisoni, Anolis porcatus, Anolis allogus, Anolis homolechis, and Anolis sagrei). We carried out a comparative analysis of these genome assemblies to investigate the genetic changes that occurred during their diversification.

Results

We reconstructed novel draft genomes with relatively long scaffolds and high gene completeness, with the scaffold N50 ranging from 5.56 to 39.79 Mb and vertebrate Benchmarking Universal Single-Copy Orthologs completeness ranging from 77.5% to 86.9%. Comparing the repeat element compositions and landscapes revealed differences in the accumulation process between Cuban trunk-crown and trunk-ground species and separate expansions of several families of LINE in each Cuban trunk-ground species. Duplicated gene analysis suggested that the proportional differences in duplicated gene numbers among Cuban Anolis lizards may be associated with differences in their habitat ranges. Additionally, Pairwise Sequentially Markovian Coalescent analysis suggested that the effective population sizes of each species may have been affected by Cuba’s geohistory.

Conclusions

We provide draft genomes of six Cuban Anolis lizards and detected species and lineage-specific transposon accumulation and gene copy number changes that may be involved in adaptive evolution. The change processes in the past effective population size was also estimated, and the factors involved were inferred. These results provide new insights into the genetic basis of Anolis lizard diversification and are expected to serve as a stepping stone for the further elucidation of their diversification mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02086-7.

Chromosome-scale genome assembly of the brown anole (Anolis sagrei), an emerging model species

Rapid technological improvements are democratizing access to high quality, chromosome-scale genome assemblies. No longer the domain of only the most highly studied model organisms, now non-traditional and emerging model species can be genome-enabled using a combination of sequencing technologies and assembly software. Consequently, old ideas built on sparse sampling across the tree of life have recently been amended in the face of genomic data drawn from a growing number of high-quality reference genomes. Arguably the most valuable are those long-studied species for which much is already known about their biology; what many term emerging model species. Here, we report a highly complete chromosome-scale genome assembly for the brown anole, Anolis sagrei - a lizard species widely studied across a variety of disciplines and for which a high-quality reference genome was long overdue. This assembly exceeds the vast majority of existing reptile and snake genomes in contiguity (N50 = 253.6 Mb) and annotation completeness. Through the analysis of this genome and population resequence data, we examine the history of repetitive element accumulation, identify the X chromosome, and propose a hypothesis for the evolutionary history of fusions between autosomes and the X that led to the sex chromosomes of A. sagrei.

Population Genomics of Variegated Toad-Headed Lizard Phrynocephalus versicolor and Its Adaptation to the Colorful Sand of the Gobi Desert

The variegated toad-headed agama, Phrynocephalus versicolor, lives in the arid landscape of the Chinese Gobi Desert. We analyzed populations from three different locations which vary in substrate color and altitude: Heishankou (HSK), Guazhou County (GZ), and Ejin Banner (EJN). The substrate color is either light-yellow (GZ-y), yellow (EJN-y), or black (HSK-b); the corresponding lizard population colors largely match their substrate in the degree of melanism. We assembled the P. versicolor genome and sequenced over 90 individuals from the three different populations. Genetic divergence between populations corresponds to their geographic distribution. We inferred the genetic relationships among these populations and used selection scans and differential expression to identify genes that show signatures of selection. Slc2a11 and akap12, among other genes, are highly differentiated and may be responsible for pigment adaptation to substrate color in P. versicolor.

Genomic adaptations for arboreal locomotion in Asian flying treefrogs

SignificanceTo adapt to arboreal lifestyles, treefrogs have evolved a suite of complex traits that support vertical movement and gliding, thus presenting a unique case for studying the genetic basis for traits causally linked to vertical niche expansion. Here, based on two de novo-assembled Asian treefrog genomes, we determined that genes involved in limb development and keratin cytoskeleton likely played a role in the evolution of their climbing systems. Behavioral and morphological evaluation and time-ordered gene coexpression network analysis revealed the developmental patterns and regulatory pathways of the webbed feet used for gliding in Rhacophorus kio.

A single locus regulates a female-limited color pattern polymorphism in a reptile

Animal coloration is often expressed in periodic patterns that can arise from differential cell migration, yet how these processes are regulated remains elusive. We show that a female-limited polymorphism in dorsal patterning (diamond/chevron) in the brown anole is controlled by a single Mendelian locus. This locus contains the gene CCDC170 that is adjacent to, and coexpressed with, the Estrogen receptor-1 gene, explaining why the polymorphism is female limited. CCDC170 is an organizer of the Golgi-microtubule network underlying a cell's ability to migrate, and the two segregating alleles encode structurally different proteins. Our agent-based modeling of skin development demonstrates that, in principle, a change in cell migratory behaviors is sufficient to switch between the two morphs. These results suggest that CCDC170 might have been co-opted as a switch between color patterning morphs, likely by modulating cell migratory behaviors.

Genome-scale data reveal deep lineage divergence and a complex demographic history in the Texas horned lizard (Phrynosoma cornutum) throughout the southwestern and central US

The southwestern and central United States serve as an ideal region to test alternative hypotheses regarding biotic diversification. Genomic data can now be combined with sophisticated computational models to quantify the impacts of paleoclimate change, geographic features, and habitat heterogeneity on spatial patterns of genetic diversity. In this study, we combine thousands of genotyping-by-sequencing (GBS) loci with mtDNA sequences (ND1) from the Texas horned lizard (Phrynosoma cornutum) to quantify relative support for different catalysts of diversification. Phylogenetic and clustering analyses of the GBS data indicate support for at least three primary populations. The spatial distribution of populations appears concordant with habitat type, with desert populations in AZ and NM showing the largest genetic divergence from the remaining populations. The mtDNA data also support a divergent desert population, but other relationships differ and suggest mtDNA introgression. Genotype–environment association with bioclimatic variables supports divergence along precipitation gradients more than along temperature gradients. Demographic analyses support a complex history, with introgression and gene flow playing an important role during diversification. Bayesian multispecies coalescent analyses with introgression (MSci) analyses also suggest that gene flow occurred between populations. Paleo-species distribution models support two southern refugia that geographically correspond to contemporary lineages. We find that divergence times are underestimated and population sizes are overestimated when introgression occurred and is ignored in coalescent analyses, and furthermore, inference of ancient introgression events and demographic history is sensitive to inclusion of a single recently admixed sample. Our analyses cannot refute the riverine barrier or glacial refugia hypotheses. Results also suggest that populations are continuing to diverge along habitat gradients. Finally, the strong evidence of admixture, gene flow, and mtDNA introgression among populations suggests that P. cornutum should be considered a single widespread species under the General Lineage Species Concept.

Interspecific Gene Flow and Mitochondrial Genome Capture During the Radiation of Jamaican Anolis Lizards (Squamata; Iguanidae)

Gene flow and reticulation are increasingly recognized as important processes in the diversification of many taxonomic groups. With the increasing ease of collecting genomic data and the development of multispecies coalescent network approaches, such reticulations can be accounted for when inferring phylogeny and diversification. Caribbean Anolis lizards are a classic example of an adaptive radiation in which species have independently radiated on the islands of the Greater Antilles into the same ecomorph classes. Within the Jamaican radiation at least one species, A. opalinus, has been documented to be polyphyletic in its mitochondrial DNA, which could be the result of an ancient reticulation event or incomplete lineage sorting. Here we generate mtDNA and genotyping-by-sequencing (GBS) data and implement gene-tree, species-tree, and multispecies coalescent network methods to infer the diversification of this group. Our mtDNA gene-tree recovers the same relationships previously inferred for this group, which is strikingly different from the species-tree inferred from our GBS data. Posterior predictive simulations suggest that our genomic data violate commonly adopted assumptions of the multispecies coalescent model, so we use network approaches to infer phylogenetic relationships. The inferred network topology contains a reticulation event but does not explain the mtDNA polyphyly observed in this group, however coalescent simulations suggest that the observed mtDNA topology is likely the result of past introgression. How common a signature of gene flow and reticulation is across the radiation of Anolis is unknown; however, the reticulation events that we demonstrate here may have allowed for adaptive evolution, as has been suggested in other, more recent adaptive radiations.

Ocular elongation and retraction in foveated reptiles

BACKGROUND

Pronounced asymmetric changes in ocular globe size during eye development have been observed in a number of species ranging from humans to lizards. In contrast, largely symmetric changes in globe size have been described for other species like rodents. We propose that asymmetric changes in the three-dimensional structure of the developing eye correlate with the types of retinal remodeling needed to produce areas of high photoreceptor density. To test this idea, we systematically examined three-dimensional aspects of globe size as a function of eye development in the bifoveated brown anole, Anolis sagrei.

RESULTS

During embryonic development, the anole eye undergoes dynamic changes in ocular shape. Initially spherical, the eye elongates in the presumptive foveal regions of the retina and then proceeds through a period of retraction that returns the eye to its spherical shape. During this period of retraction, pit formation and photoreceptor cell packing are observed. We found a similar pattern of elongation and retraction associated with the single fovea of the veiled chameleon, Chamaeleo calyptratus.

CONCLUSIONS

These results, together with those reported for other foveated species, support the idea that areas of high photoreceptor packing occur in regions where the ocular globe asymmetrically elongates and retracts during development.

The genomic basis of evolutionary differentiation among honey bees

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

A highly conserved ontogenetic limb allometry and its evolutionary significance in the adaptive radiation of Anolis lizards

Diversifications often proceed along highly conserved, evolutionary trajectories. These patterns of covariation arise in ontogeny, which raises the possibility that adaptive morphologies are biased towards trait covariations that resemble growth trajectories. Here, we test this prediction in the diverse clade of Anolis lizards by investigating the covariation of embryonic growth of 13 fore- and hindlimb bones in 15 species, and compare these to the evolutionary covariation of these limb bones across 267 Anolis species. Our results demonstrate that species differences in relative limb length are established already at hatching, and are resulting from both differential growth and differential sizes of cartilaginous anlagen. Multivariate analysis revealed that Antillean Anolis share a common ontogenetic allometry that is characterized by positive allometric growth of the long bones relative to metapodial and phalangeal bones. This major axis of ontogenetic allometry in limb bones deviated from the major axis of evolutionary allometry of the Antillean Anolis and the two clades of mainland Anolis lizards. These results demonstrate that the remarkable diversification of locomotor specialists in Anolis lizards are accessible through changes that are largely independent from ontogenetic growth trajectories, and therefore likely to be the result of modifications that manifest at the earliest stages of limb development.

Monarch Butterfly Migration as an Integrative Model of Complex Trait Evolution

AbstractUnderstanding the genetic architecture of complex trait adaptation in natural populations requires the continued development of tractable models that explicitly confront organismal and environmental complexity. A decade of high-throughput sequencing-based investigations into the genomic basis of migration points to an integrative framework that incorporates quantitative genetics, evolutionary developmental biology, phenotypic plasticity, and epigenetics to explain migration evolution. In this perspective, I argue that the transcontinental migration of the monarch butterfly (Danaus plexippus) can serve as a compelling system to study the mechanism of evolutionary lability of a complex trait. Monarchs show significant phenotypic and genotypic diversity across their global range, with phenotypic switching that allows for explicit study of evolutionary lability. A developmental approach for elucidating how migratory traits are generated and functionally integrated will be important for understanding the evolution of monarch migration traits. I propose a plasticity threshold model to describe migration lability, and I describe novel functional techniques that will help resolve open questions and model assumptions. I conclude by considering the relationships between adaptive genetic architecture, anthropogenic climate change, and conservation management practice and the timeliness of the monarch migration model to illuminate these connections given the rapid decline of the North American migration.

Molecular Evolution of Ecological Specialisation: Genomic Insights from the Diversification of Murine Rodents

Adaptive radiations are characterized by the diversification and ecological differentiation of species, and replicated cases of this process provide natural experiments for understanding the repeatability and pace of molecular evolution. During adaptive radiation, genes related to ecological specialization may be subject to recurrent positive directional selection. However, it is not clear to what extent patterns of lineage-specific ecological specialization (including phenotypic convergence) are correlated with shared signatures of molecular evolution. To test this, we sequenced whole exomes from a phylogenetically dispersed sample of 38 murine rodent species, a group characterized by multiple, nested adaptive radiations comprising extensive ecological and phenotypic diversity. We found that genes associated with immunity, reproduction, diet, digestion, and taste have been subject to pervasive positive selection during the diversification of murine rodents. We also found a significant correlation between genome-wide positive selection and dietary specialization, with a higher proportion of positively selected codon sites in derived dietary forms (i.e., carnivores and herbivores) than in ancestral forms (i.e., omnivores). Despite striking convergent evolution of skull morphology and dentition in two distantly related worm-eating specialists, we did not detect more genes with shared signatures of positive or relaxed selection than in a nonconvergent species comparison. Although a small number of the genes we detected can be incidentally linked to craniofacial morphology or diet, protein-coding regions are unlikely to be the primary genetic basis of this complex convergent phenotype. Our results suggest a link between positive selection and derived ecological phenotypes, and highlight specific genes and general functional categories that may have played an integral role in the extensive and rapid diversification of murine rodents.

Accelerated Evolution of Tissue-Specific Genes Mediates Divergence Amidst Gene Flow in European Green Lizards

The European green lizards of the Lacerta viridis complex consist of two closely related species, L. viridis and Lacerta bilineata that split less than 7 million years ago in the presence of gene flow. Recently, a third lineage, referred to as the “Adriatic” was described within the L. viridis complex distributed from Slovenia to Greece. However, whether gene flow between the Adriatic lineage and L. viridis or L. bilineata has occurred and the evolutionary processes involved in their diversification are currently unknown. We hypothesized that divergence occurred in the presence of gene flow between multiple lineages and involved tissue-specific gene evolution. In this study, we sequenced the whole genome of an individual of the Adriatic lineage and tested for the presence of gene flow amongst L. viridis, L. bilineata, and Adriatic. Additionally, we sequenced transcriptomes from multiple tissues to understand tissue-specific effects. The species tree supports that the Adriatic lineage is a sister taxon to L. bilineata. We detected gene flow between the Adriatic lineage and L. viridis suggesting that the evolutionary history of the L. viridis complex is likely shaped by gene flow. Interestingly, we observed topological differences between the autosomal and Z-chromosome phylogenies with a few fast evolving genes on the Z-chromosome. Genes highly expressed in the ovaries and strongly co-expressed in the brain experienced accelerated evolution presumably contributing to establishing reproductive isolation in the L. viridis complex.

Adaptive Radiation Genomics of Two Ecologically Divergent Hawai'ian Honeycreepers: The 'akiapōlā'au and the Hawai'i 'amakihi

The Hawai'ian honeycreepers (drepanids) are a classic example of adaptive radiation: they adapted to a variety of novel dietary niches, evolving a wide range of bill morphologies. Here we investigated genomic diversity, demographic history, and genes involved in bill morphology phenotypes in 2 honeycreepers: the 'akiapōlā'au (Hemignathus wilsoni) and the Hawai'i 'amakihi (Chlorodrepanis virens). The 'akiapōlā'au is an endangered island endemic, filling the "woodpecker" niche by using a unique bill morphology, while the Hawai'i 'amakihi is a dietary generalist common on the islands of Hawai'i and Maui. We de novo sequenced the 'akiapōlā'au genome and compared it to the previously sequenced 'amakihi genome. The 'akiapōlā'au is far less heterozygous and has a smaller effective population size than the 'amakihi, which matches expectations due to its smaller census population and restricted ecological niche. Our investigation revealed genomic islands of divergence, which may be involved in the honeycreeper radiation. Within these islands of divergence, we identified candidate genes (including DLK1, FOXB1, KIF6, MAML3, PHF20, RBP1, and TIMM17A) that may play a role in honeycreeper adaptations. The gene DLK1, previously shown to influence Darwin's finch bill size, may be related to honeycreeper bill morphology evolution, while the functions of the other candidates remain unknown.

Phenotypic Convergence Is Not Mirrored at the Protein Level in a Lizard Adaptive Radiation

There are many compelling examples of molecular convergence at individual genes. However, the prevalence and the relative importance of adaptive genome-wide convergence remain largely unknown. Many recent works have reported striking examples of excess genome-wide convergence, but some of these studies have been called into question because of the use of inappropriate null models. Here, we sequenced and compared the genomes of 12 species of anole lizards that have independently converged on suites of adaptive behavioral and morphological traits. Despite extensive searches for a genome-wide signature of molecular convergence, we found no evidence supporting molecular convergence at specific amino acids either at individual genes or at genome-wide comparisons; we also uncovered no evidence supporting an excess of adaptive convergence in the rates of amino acid substitutions within genes. Our findings indicate that comprehensive phenotypic convergence is not mirrored at genome-wide protein-coding levels in anoles, and therefore, that adaptive phenotypic convergence is likely not constrained by the evolution of many specific protein sequences or structures.

Sequencing platform shifts provide opportunities but pose challenges for combining genomic data sets

Technological advances in DNA sequencing over the last decade now permit the production and curation of large genomic data sets in an increasing number of nonmodel species. Additionally, these new data provide the opportunity for combining data sets, resulting in larger studies with a broader taxonomic range. Whilst the development of new sequencing platforms has been beneficial, resulting in a higher throughput of data at a lower per-base cost, shifts in sequencing technology can also pose challenges for those wishing to combine new sequencing data with data sequenced on older platforms. Here, we outline the types of studies where the use of curated data might be beneficial, and highlight potential biases that might be introduced by combining data from different sequencing platforms. As an example of the challenges associated with combining data across sequencing platforms, we focus on the impact of the shift in Illumina's base calling technology from a four-channel system to a two-channel system. We caution that when data are combined from these two systems, erroneous guanine base calls that result from the two-channel chemistry can make their way through a bioinformatic pipeline, eventually leading to inaccurate and potentially misleading conclusions. We also suggest solutions for dealing with such potential artefacts, which make samples sequenced on different sequencing platforms appear more differentiated from one another than they really are. Finally, we stress the importance of archiving tissue samples and the associated sequences for the continued reproducibility and reusability of sequencing data in the face of ever-changing sequencing platform technology.

Deep Convergence, Shared Ancestry, and Evolutionary Novelty in the Genetic Architecture of Heliconius Mimicry

Convergent evolution can occur through different genetic mechanisms in different species. It is now clear that convergence at the genetic level is also widespread, and can be caused by either (i) parallel genetic evolution, where independently evolved convergent mutations arise in different populations or species, or (ii) collateral evolution in which shared ancestry results from either ancestral polymorphism or introgression among taxa. The adaptive radiation of Heliconius butterflies shows color pattern variation within species, as well as mimetic convergence between species. Using comparisons from across multiple hybrid zones, we use signals of shared ancestry to identify and refine multiple putative regulatory elements in Heliconius melpomene and its comimics, Heliconius elevatus and Heliconius besckei, around three known major color patterning genes: optix, WntA, and cortex While we find that convergence between H. melpomene and H. elevatus is caused by a complex history of collateral evolution via introgression in the Amazon, convergence between these species in the Guianas appears to have evolved independently. Thus, we find adaptive convergent genetic evolution to be a key driver of regulatory changes that lead to rapid phenotypic changes. Furthermore, we uncover evidence of parallel genetic evolution at some loci around optix and WntA in H. melpomene and its distant comimic Heliconius erato Ultimately, we show that all three of convergence, conservation, and novelty underlie the modular architecture of Heliconius color pattern mimicry.

Disentangling the determinants of transposable elements dynamics in vertebrate genomes using empirical evidences and simulations

The interactions between transposable elements (TEs) and their hosts constitute one of the most profound co-evolutionary processes found in nature. The population dynamics of TEs depends on factors specific to each TE families, such as the rate of transposition and insertional preference, the demographic history of the host and the genomic landscape. How these factors interact has yet to be investigated holistically. Here we are addressing this question in the green anole (Anolis carolinensis) whose genome contains an extraordinary diversity of TEs (including non-LTR retrotransposons, SINEs, LTR-retrotransposons and DNA transposons). We observed a positive correlation between recombination rate and frequency of TEs and densities for LINEs, SINEs and DNA transposons. For these elements, there was a clear impact of demography on TE frequency and abundance, with a loss of polymorphic elements and skewed frequency spectra in recently expanded populations. On the other hand, some LTR-retrotransposons displayed patterns consistent with a very recent phase of intense amplification. To determine how demography, genomic features and intrinsic properties of TEs interact we ran simulations using SLiM3. We determined that i) short TE insertions are not strongly counter-selected, but long ones are, ii) neutral demographic processes, linked selection and preferential insertion may explain positive correlations between average TE frequency and recombination, iii) TE insertions are unlikely to have been massively recruited in recent adaptation. We demonstrate that deterministic and stochastic processes have different effects on categories of TEs and that a combination of empirical analyses and simulations can disentangle these mechanisms.

A brief review on the natural history, venomics and the medical importance of bushmaster (Lachesis) pit viper snakes

Snakes of the genus Lachesis, commonly known as bushmasters, are the largest venomous snakes in the Americas. Because these snakes have their habitats in areas of remote forests they are difficult to find, and consequently there are few studies of Lachesis taxa in their natural ecosystems. Bushmasters are distributed in tropical forest areas of South and Central America. In Brazil they can be found in the Amazon Rainforest and the Atlantic Forest. Despite the low incidence of cases, laquetic envenoming causes severe permanent sequelae due to the high amount of inoculated venom. These accidents are characterized by local pain, hemorrhage and myonecrosis that can be confused with bothropic envenomings. However, victims of Lachesis bites develop symptoms characteristic of Lachesis envenoming, known as vagal syndrome. An important message of this bibliographic synthesis exercise is that, despite having the proteomic profiles of all the taxa of the genus available, very few structure-function correlation studies have been carried out. Therefore the motivation for this review was to fill a gap in the literature on the genus Lachesis, about which there is no recent review. Here we discuss data scattered in a number of original articles published in specialized journals, spanning the evolutionary history and extant phylogeographic distribution of the bushmasters, their venom composition and diet, as well as the pathophysiology of their bites to humans and the biological activities and possible biotechnological applicability of their venom toxins.

Phylogenomics and molecular species delimitation reveals great cryptic diversity of leaf-toed geckos (Phyllodactylidae: Phyllodactylus), ancient origins, and diversification in Mexico

We utilize the efficient GBS technique to obtain thousands of nuclear loci and SNPs to reconstruct the evolutionary history of Mexican leaf-toed geckos (Phyllodactylus). Through the incorporation of unprecedented sampling for this group of geckos, in combination with genomic data analysis, we generate mostly consistent phylogenetic hypotheses using two approaches: supermatrix and coalescent-based inference. All topologies depict three, mutually exclusive major clades. Clade I comprises P. bordai and all species closer to P. bordai than to any other Phyllodactylus. Clade II comprises P. nocticolus and all species closer to P. nocticolus than to any other Phyllodactylus. Clade III comprises P. tuberculosus and all species closer to P. tuberculosus than to any other Phyllodactylus. Analyses estimate the age for the most recent common ancestor of Phyllodactylus in the Eocene (~43 mya), and the ancestors of each major clade date to the Eocene-Oligocene transition (32-36 mya). This group includes one late-Eocene lineage (P. bordai), Oligocene lineages (P. paucituberculatus, P. delcampi), but also topological patterns that indicate a recent radiation occurred during the Pleistocene on islands in the Gulf of California. The wide spatial and temporal scale indicates a complex and unique biogeographic history for each major clade. The 33 species delimited by BPP and stepping-stone BFD*coalescent based genomic approaches reflect this history. This diversity delimited for Mexican leaf-toed geckos demonstrates a vast underestimation in the number of species based on morphological data alone.

Evolutionary lability in Hox cluster structure and gene expression in Anolis lizards

Hox genes orchestrate development by patterning the embryonic axis. Vertebrate Hox genes are arranged in four compact clusters, and the spacing between genes is assumed to be crucial for their function. The genomes of squamate reptiles are unusually rich and variable in transposable elements (TEs), and it has been suggested that TE invasion is responsible for the Hox cluster expansion seen in snakes and lizards. Using de novo TE prediction on 17 genomes of lizards and snakes, I show that TE content of Hox clusters are generally 50% lower than genome‐wide TE levels. However, two distantly related lizards of the species‐rich genus Anolis have Hox clusters with a TE content that approaches genomic levels. The age distribution of TEs in Anolis lizards revealed that peaks of TE activity broadly coincide with speciation events. In accordance with theoretical models of Hox cluster regulation, I find that Anolis species with many TEs in their Hox clusters show aberrant Hox gene expression patterns, suggesting a functional link between TE accumulation and embryonic development. These results are consistent with the hypothesis that TEs play a role in developmental processes as well as in evolutionary diversifications.

Widespread cis-regulatory convergence between the extinct Tasmanian tiger and gray wolf

The extinct marsupial Tasmanian tiger, or thylacine, and the eutherian gray wolf are among the most widely recognized examples of convergent evolution in mammals. Despite being distantly related, these large predators independently evolved extremely similar craniofacial morphologies, and evidence suggests that they filled similar ecological niches. Previous analyses revealed little evidence of adaptive convergence between their protein-coding genes. Thus, the genetic basis of their convergence is still unclear. Here, we identified candidate craniofacial cis-regulatory elements across vertebrates and compared their evolutionary rates in the thylacine and wolf, revealing abundant signatures of convergent positive selection. Craniofacial thylacine-wolf accelerated regions were enriched near genes involved in TGF beta (TGFB) and BMP signaling, both of which are key morphological signaling pathways with critical roles in establishing the identities and boundaries between craniofacial tissues. Similarly, enhancers of genes involved in craniofacial nerve development showed convergent selection and involvement in these pathways. Taken together, these results suggest that adaptation in cis-regulators of TGF beta and BMP signaling may provide a mechanism to explain the coevolution of developmentally and functionally integrated craniofacial structures in these species. We also found that despite major structural differences in marsupial and eutherian brains, accelerated regions in both species were common near genes with roles in brain development. Our findings support the hypothesis that, relative to protein-coding genes, positive selection on cis-regulatory elements is likely to be an essential driver of adaptive convergent evolution and may underpin thylacine-wolf phenotypic similarities.

Selection at behavioural, developmental and metabolic genes is associated with the northward expansion of a successful tropical colonizer

What makes a species able to colonize novel environments? This question is key to understand the dynamics of adaptive radiations and ecological niche shifts, but the mechanisms that underlie expansion into novel habitats remain poorly understood at a genomic scale. Lizards from the genus Anolis are typically tropical, and the green anole (Anolis carolinensis) constitutes an exception since it expanded into temperate North America from subtropical Florida. Thus, we used the green anole as a model to investigate signatures of selection associated with colonization of a new environment, namely temperate North America. To this end, we analysed 29 whole-genome sequences, covering the entire native range of the species. We used a combination of recent methods to quantify both positive and balancing selection in northern populations, including F_ST outlier methods, machine learning and ancestral recombination graphs. We naively scanned for genes of interest and assessed the overlap between multiple tests. Strikingly, we identified many genes involved in behaviour, suggesting that the recent successful colonization of northern environments may have been linked to behavioural shifts as well as physiological adaptation. Using a candidate genes strategy, we determined that genes involved in response to cold or behaviour displayed more frequently signals of selection, while controlling for local recombination rate, gene clustering and gene length. In addition, we found signatures of balancing selection at immune genes in all investigated genetic groups, but also at genes involved in neuronal and anatomical development.

Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation

BACKGROUND

Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes.

FINDINGS

Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated.

CONCLUSION

The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.

Molecular Adaptations for Sensing and Securing Prey and Insight into Amniote Genome Diversity from the Garter Snake Genome

Colubridae represents the most phenotypically diverse and speciose family of snakes, yet no well-assembled and annotated genome exists for this lineage. Here, we report and analyze the genome of the garter snake, Thamnophis sirtalis, a colubrid snake that is an important model species for research in evolutionary biology, physiology, genomics, behavior, and the evolution of toxin resistance. Using the garter snake genome, we show how snakes have evolved numerous adaptations for sensing and securing prey, and identify features of snake genome structure that provide insight into the evolution of amniote genomes. Analyses of the garter snake and other squamate reptile genomes highlight shifts in repeat element abundance and expansion within snakes, uncover evidence of genes under positive selection, and provide revised neutral substitution rate estimates for squamates. Our identification of Z and W sex chromosome-specific scaffolds provides evidence for multiple origins of sex chromosome systems in snakes and demonstrates the value of this genome for studying sex chromosome evolution. Analysis of gene duplication and loss in visual and olfactory gene families supports a dim-light ancestral condition in snakes and indicates that olfactory receptor repertoires underwent an expansion early in snake evolution. Additionally, we provide some of the first links between secreted venom proteins, the genes that encode them, and their evolutionary origins in a rear-fanged colubrid snake, together with new genomic insight into the coevolutionary arms race between garter snakes and highly toxic newt prey that led to toxin resistance in garter snakes.