Repeated evolution and the impact of evolutionary history on adaptation

Ancient biases in phenotype production drove the functional evolution of a protein family

Biological systems may be biased in the phenotypes they can access by mutation1-7, but the extent of these biases and their causal role in the evolution of extant phenotypic diversity remains unclear. There are three major challenges: it is difficult to isolate the effect of bias in the genotype-phenotype (GP) map from that of natural selection in producing natural diversity6,8-11, the universe of possible genotypes and phenotypes is so vast and complex that a direct characterization has been impossible, and most extant phenotypes evolved long ago in species whose GP maps cannot be recovered. Here we develop exhaustive multi-phenotype deep mutational scanning to experimentally characterize the complete GP maps of two reconstructed ancestral steroid receptor proteins, which existed during an ancient phylogenetic interval when a new phenotype-specific binding of a new DNA response element-evolved12. We measured all possible DNA specificity phenotypes encoded by all possible amino acid combinations at sites in the protein's DNA binding interface. We found that the ancestral GP maps are structured by strong global bias-unequal propensity to encode the various phenotypes-and extreme heterogeneity in the phenotypes accessible around each genotype, which strongly affect evolution on both long and short timescales. Distinct biases in the two ancestral maps steered evolution toward the lineage-specific functional phenotypes that evolved during history. Our findings establish that ancient biases in the GP relationship were causal factors in the evolutionary process that produced the present-day patterns of phenotypic conservation and diversity in this protein family.

Divergence time and environmental similarity predict the strength of morphological convergence in stick and leaf insects

Independent evolution of similar traits in lineages inhabiting similar environments (convergent or repeated evolution) is often taken as evidence for adaptation by natural selection, and used to illustrate the predictability of evolution. Yet convergence is rarely perfect for two reasons. First, environments may not be as similar as they appear. Second, responses to selection are contingent upon available genetic variation and independent lineages may differ in the alleles, genetic backgrounds, and even the developmental mechanisms responsible for the phenotypes in question. Both impediments to convergence are predicted to increase as the length of time separating two lineages increases, making it difficult to discern their relative importance. We quantified environmental similarity and the extent of convergence to show how habitat and divergence time each contribute to observed patterns of morphological evolution in 212 species of stick and leaf insects (order Phasmatodea). Dozens of phasmid lineages independently colonized similar habitats, repeatedly evolving in parallel directions on a 23-trait morphospace, though the magnitude and direction of these shifts varied. Lineages converging toward more similar environments ended up closer on the morphospace, as did closely related lineages, and closely related lineages followed more parallel evolutionary trajectories to arrive there than more distantly related ones. Remarkably, after accounting for habitat similarity, we show that divergence time reduced the extent of convergence at a constant rate across more than 100 My of separation, suggesting even the magnitude of contingency can be predictable, given sufficient spans of time.

Widespread adaptive evolution in angiosperm photosystems provides insight into the evolution of photosystem II repair

Oxygenic photosynthesis generates the initial energy source that fuels nearly all life on Earth. At the heart of the process are the photosystems, which are pigment binding multiprotein complexes that catalyze the first step of photochemical conversion of light energy into chemical energy. Here, we investigate the molecular evolution of the plastid-encoded photosystem subunits at single-residue resolution across 773 angiosperm species. We show that despite an extremely high level of conservation, 7% of residues in the photosystems, spanning all photosystem subunits, exhibit hallmarks of adaptive evolution. Through in silico modeling of these adaptive substitutions, we uncover the impact of these changes on the predicted properties of the photosystems, focusing on their effects on cofactor binding and intersubunit interface formation. By analyzing these cohorts of changes, we reveal that evolution has repeatedly altered the interaction between Photosystem II and its D1 subunit in a manner that is predicted to reduce the energetic barrier for D1 turnover and photosystem repair. Together, these results provide insight into the trajectory of photosystem adaptation during angiosperm evolution.

Learning genotype-phenotype associations from gaps in multi-species sequence alignments

Understanding the genetic basis of phenotypic variation is fundamental to biology. Here we introduce GAP, a novel machine learning framework for predicting binary phenotypes from gaps in multi-species sequence alignments. GAP employs a neural network to predict the presence or absence of phenotypes solely from alignment gaps, contrasting with existing tools that require additional and often inaccessible input data. GAP can be applied to three distinct problems: predicting phenotypes in species from known associated genomic regions, pinpointing positions within such regions that are important for predicting phenotypes, and extracting sets of candidate regions associated with phenotypes. We showcase the utility of GAP by exploiting the well-known association between the L-gulonolactone oxidase (Gulo) gene and vitamin C synthesis, demonstrating its perfect prediction accuracy in 34 vertebrates. This exceptional performance also applies more generally, with GAP achieving high accuracy and power on a large simulated dataset. Moreover, predictions of vitamin C synthesis in species with unknown status mirror their phylogenetic relationships, and positions with high predictive importance are consistent with those identified by previous studies. Last, a genome-wide application of GAP identifies many additional genes that may be associated with vitamin C synthesis, and analysis of these candidates uncovers functional enrichment for immunity, a widely recognized role of vitamin C. Hence, GAP represents a simple yet useful tool for predicting genotype-phenotype associations and addressing diverse evolutionary questions from data available in a broad range of study systems.

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.

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.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

Convergence and contingency in the evolution of a specialized mode of life: multiple origins and high disparity of rock-boring bivalves

Evolutionary adaptation to novel, specialized modes of life is often associated with a close mapping of form to the new function, resulting in narrow morphological disparity. For bivalve molluscs, endolithy (rock-boring) has biomechanical requirements thought to diverge strongly from those of ancestral functions. However, endolithy in bivalves has originated at least eight times. Three-dimensional morphometric data representing 75 species from approximately 94% of extant endolithic genera and families, along with 310 non-endolithic species in those families, show that endolithy is evolutionarily accessible from many different morphological starting points. Although some endoliths appear to converge on certain shell morphologies, the range of endolith shell form is as broad as that belonging to any other bivalve substrate use. Nevertheless, endolithy is a taxon-poor function in Bivalvia today. This limited richness does not derive from origination within source clades having significantly low origination or high extinction rates, and today's endoliths are not confined to low-diversity biogeographic regions. Instead, endolithy may be limited by habitat availability. Both determinism (as reflected by convergence among distantly related taxa) and contingency (as reflected by the endoliths that remain close to the disparate morphologies of their source clades) underlie the occupation of endolith morphospace.

Multiple pathways to herbivory underpinned deep divergences in ornithischian evolution

The extent to which evolution is deterministic is a key question in biology,^{1,2,3,4,5,6,7,8,9} with intensive debate on how adaptation^{6,10,11,12,13} and constraints^14,15,16 might canalize solutions to ecological challenges.^4,5,6 Alternatively, unique adaptations^1,9,17 and phylogenetic contingency^1,3,18 may render evolution fundamentally unpredictable.³ Information from the fossil record is critical to this debate,^1,2,11 but performance data for extinct taxa are limited.⁷ This knowledge gap is significant, as general morphology may be a poor predictor of biomechanical performance.^17,19,20 High-fiber herbivory originated multiple times within ornithischian dinosaurs,²¹ making them an ideal clade for investigating evolutionary responses to similar ecological pressures.²² However, previous biomechanical modeling studies on ornithischian crania^17,23,24,25 have not compared early-diverging taxa spanning independent acquisitions of herbivory. Here, we perform finite-element analysis on the skull of five early-diverging members of the major ornithischian clades to characterize morphofunctional pathways to herbivory. Results reveal limited functional convergence among ornithischian clades, with each instead achieving comparable performance, in terms of reconstructed patterns and magnitudes of functionally induced stress, through different adaptations of the feeding apparatus. Thyreophorans compensated for plesiomorphic low performance through increased absolute size, heterodontosaurids expanded jaw adductor muscle volume, ornithopods increased jaw system efficiency, and ceratopsians combined these approaches. These distinct solutions to the challenges of herbivory within Ornithischia underpinned the success of this diverse clade. Furthermore, the resolution of multiple solutions to equivalent problems within a single clade through macroevolutionary time demonstrates that phenotypic evolution is not necessarily predictable, instead arising from the interplay of adaptation, innovation, contingency, and constraints.^{1,2,3,7,8,9,18}.

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.

Genomic evidence for contrasting patterns of host-associated genetic differentiation across shared host-plant species in leaf- and bud-galling sawflies

Resource specialization and ecological speciation arising through host-associated genetic differentiation (HAD) are frequently invoked as an explanation for the high diversity of plant-feeding insects and other organisms with a parasitic lifestyle. While genetic studies have demonstrated numerous examples of HAD in insect herbivores, the rarity of comparative studies means that we still lack an understanding of how deterministic HAD is, and whether patterns of host shifts can be predicted over evolutionary timescales. We applied genome-wide single nucleotide polymorphism and mitochondrial DNA sequence data obtained through genome resequencing to define species limits and to compare host-plant use in population samples of leaf- and bud-galling sawflies (Hymenoptera: Tenthredinidae: Nematinae) collected from seven shared willow (Salicaceae: Salix) host species. To infer the repeatability of long-term cophylogenetic patterns, we also contrasted the phylogenies of the two galler groups with each other as well as with the phylogeny of their Salix hosts estimated based on RADseq data. We found clear evidence for host specialization and HAD in both of the focal galler groups, but also that leaf gallers are more specialized to single host species compared with most bud gallers. In contrast to bud gallers, leaf gallers also exhibited statistically significant cophylogenetic signal with their Salix hosts. The observed discordant patterns of resource specialization and host shifts in two related galler groups that have radiated in parallel across a shared resource base indicate a lack of evolutionary repeatability in the focal system, and suggest that short- and long-term host use and ecological diversification in plant-feeding insects are dominated by stochasticity and/or lineage-specific effects.

Genes that are Used Together are More Likely to be Fused Together in Evolution by Mutational Mechanisms: A Bioinformatic Test of the Used-Fused Hypothesis

Cases of parallel or recurrent gene fusions in evolution as well as in genetic disease and cancer are difficult to explain, because unlike point mutations, they can require the repetition of a similar configuration of multiple breakpoints rather than the repetition of a single point mutation. The used-together-fused-together hypothesis holds that genes that are used together repeatedly and persistently in a specific context are more likely to undergo fusion mutation in the course of evolution for mechanistic reasons. This hypothesis offers to explain gene fusion in both evolution and disease under one umbrella. Using bioinformatic data, we tested this hypothesis against alternatives, including that all gene pairs can fuse by random mutation, but among pairs thus fused, those that had interacted previously are more likely to be favored by selection. Results show that across multiple measures of gene interaction, human genes whose orthologs are fused in one or more species are more likely to interact with each other than random pairs of genes of the same genomic distance between pair members; that an overlap exists between genes that fused in the course of evolution in non-human species and genes that undergo fusion in human cancers; and that across six primate species studied, fusions predominate over fissions and exhibit substantial evolutionary parallelism. Together, these results support the used-together-fused-together hypothesis over its alternatives. Multiple implications are discussed, including the relevance of mutational mechanisms to the evolution of genome organization, to the distribution of fitness effects of mutation, to evolutionary parallelism and more.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11692-022-09579-9.

Convergence is Only Skin Deep: Craniofacial Evolution in Electric Fishes from South America and Africa (Apteronotidae and Mormyridae)

Apteronotidae and Mormyridae are species-rich clades of weakly electric fishes from Neotropical and Afrotropical freshwaters, respectively, known for their high morphological disparity and often regarded as a classic example of convergent evolution. Here, we use CT-imaging and 3D geometric morphometrics to quantify disparity in craniofacial morphologies, and to test the hypothesis of convergent skull-shape evolution in a phylogenetic context. For this study, we examined 391 specimens representing 78 species of Apteronotidae and Mormyridae including 30 of 37 (81%) of all valid genera with the goal to sample most of the craniofacial disparity known in these clades. We found no overlap between Apteronotidae and Mormyridae in the skull-shape morphospace using PCA and a common landmark scheme, and therefore no instances of complete phenotypic convergence. Instead, we found multiple potential instances of incomplete convergence, and at least one parallel shift among electric fish clades. The greatest components of shape variance in both families are the same as observed for most vertebrate clades: heterocephaly (i.e., opposite changes in relative sizes of the snout and braincase regions of the skull), and heterorhynchy (i.e., dorsoventral changes in relative snout flexion and mouth position). Mormyrid species examined here exhibit less craniofacial disparity than do apteronotids, potentially due to constraints associated with a larger brain size, ecological constraints related to food-type availability. Patterns of craniofacial evolution in these two clades depict a complex story of phenotypic divergence and convergence in which certain superficial similarities of external morphology obscure deeper osteological and presumably developmental differences of skull form and function. Among apteronotid and mormyrid electric fishes, craniofacial convergence is only skin deep.

Parallel evolutionary paths of rove beetle myrmecophiles: replaying a deep-time tape of life

The rise of ants over the past ~100 million years reshaped the biosphere, presenting ecological challenges for many organisms, but also opportunities. No insect group has been so adept at exploiting niches inside ant colonies as the rove beetles (Staphylinidae) - a global clade of>64,000 predominantly free-living predators from which numerous socially parasitic 'myrmecophile' lineages have emerged. Myrmecophilous staphylinids are specialized for colony life through changes in behavior, chemistry, anatomy, and life history that are often strikingly convergent, and hence potentially adaptive for this symbiotic way of life. Here, we examine how the interplay between ecological pressures and molecular, cellular, and neurobiological mechanisms shape the evolutionary trajectories of symbiotic lineages in this ancient, convergent system.

The repeated evolution of stripe patterns is correlated with body morphology in the adaptive radiations of East African cichlid fishes

Color patterns are often linked to the behavioral and morphological characteristics of an animal, contributing to the effectiveness of such patterns as antipredatory strategies. Species-rich adaptive radiations, such as the freshwater fish family Cichlidae, provide an exciting opportunity to study trait correlations at a macroevolutionary scale. Cichlids are also well known for their diversity and repeated evolution of color patterns and body morphology. To study the evolutionary dynamics between color patterns and body morphology, we used an extensive dataset of 461 species. A phylogenetic supertree of these species shows that stripe patterns evolved ~70 times independently and were lost again ~30 times. Moreover, stripe patterns show strong signs of correlated evolution with body elongation, suggesting that the stripes' effectiveness as antipredatory strategy might differ depending on the body shape. Using pedigree-based analyses, we show that stripes and body elongation segregate independently, indicating that the two traits are not genetically linked. Their correlation in nature is therefore likely maintained by correlational selection. Lastly, by performing a mate preference assay using a striped CRISPR-Cas9 mutant of a nonstriped species, we show that females do not differentiate between striped CRISPR mutant males and nonstriped wild-type males, suggesting that these patterns might be less important for species recognition and mate choice. In summary, our study suggests that the massive rates of repeated evolution of stripe patterns are shaped by correlational selection with body elongation, but not by sexual selection.

Phenotypic and genotypic parallel evolution in parapatric ecotypes of Senecio

The independent and repeated adaptation of populations to similar environments often results in the evolution of similar forms. This phenomenon creates a strong correlation between phenotype and environment and is referred to as parallel evolution. However, we are still largely unaware of the dynamics of parallel evolution, as well as the interplay between phenotype and genotype within natural systems. Here, we examined phenotypic and genotypic parallel evolution in multiple parapatric Dune‐Headland coastal ecotypes of an Australian wildflower, Senecio lautus. We observed a clear trait‐environment association in the system, with all replicate populations having evolved along the same phenotypic evolutionary trajectory. Similar phenotypes have arisen via mutational changes occurring in different genes, although many share the same biological functions. Our results shed light on how replicated adaptation manifests at the phenotypic and genotypic levels within populations, and highlight S. lautus as one of the most striking cases of phenotypic parallel evolution in nature.

Changes in the distribution of fitness effects and adaptive mutational spectra following a single first step towards adaptation

Historical contingency and diminishing returns epistasis have been typically studied for relatively divergent genotypes and/or over long evolutionary timescales. Here, we use Saccharomyces cerevisiae to study the extent of diminishing returns and the changes in the adaptive mutational spectra following a single first adaptive mutational step. We further evolve three clones that arose under identical conditions from a common ancestor. We follow their evolutionary dynamics by lineage tracking and determine adaptive outcomes using fitness assays and whole genome sequencing. We find that diminishing returns manifests as smaller fitness gains during the 2^nd step of adaptation compared to the 1^st step, mainly due to a compressed distribution of fitness effects. We also find that the beneficial mutational spectra for the 2^nd adaptive step are contingent on the 1^st step, as we see both shared and diverging adaptive strategies. Finally, we find that adaptive loss-of-function mutations, such as nonsense and frameshift mutations, are less common in the second step of adaptation than in the first step.

Analyses of both natural and experimental evolution suggest that adaptation depends on the evolutionary past and adaptive potential decreases over time. Here, by tracking yeast adaptation with DNA barcoding, the authors show that such evolutionary phenomena can be observed even after a single adaptive step.

Deep-time convergent evolution in animal communication presented by shared adaptations for coping with noise in lizards and other animals

Convergence in communication appears rare compared with other forms of adaptation. This is puzzling, given communication is acutely dependent on the environment and expected to converge in form when animals communicate in similar habitats. We uncover deep-time convergence in territorial communication between two groups of tropical lizards separated by over 140 million years of evolution: the Southeast Asian Draco and Caribbean Anolis. These groups have repeatedly converged in multiple aspects of display along common environmental gradients. Robot playbacks to free-ranging lizards confirmed that the most prominent convergence in display is adaptive, as it improves signal detection. We then provide evidence from a sample of the literature to further show that convergent adaptation among highly divergent animal groups is almost certainly widespread in nature. Signal evolution is therefore curbed towards the same set of adaptive solutions, especially when animals are challenged with the problem of communicating effectively in noisy environments.

Genomic basis of parallel adaptation varies with divergence in Arabidopsis and its relatives

Parallel adaptation provides valuable insight into the predictability of evolutionary change through replicated natural experiments. A steadily increasing number of studies have demonstrated genomic parallelism, yet the magnitude of this parallelism varies depending on whether populations, species, or genera are compared. This led us to hypothesize that the magnitude of genomic parallelism scales with genetic divergence between lineages, but whether this is the case and the underlying evolutionary processes remain unknown. Here, we resequenced seven parallel lineages of two Arabidopsis species, which repeatedly adapted to challenging alpine environments. By combining genome-wide divergence scans with model-based approaches, we detected a suite of 151 genes that show parallel signatures of positive selection associated with alpine colonization, involved in response to cold, high radiation, short season, herbivores, and pathogens. We complemented these parallel candidates with published gene lists from five additional alpine Brassicaceae and tested our hypothesis on a broad scale spanning ∼0.02 to 18 My of divergence. Indeed, we found quantitatively variable genomic parallelism whose extent significantly decreased with increasing divergence between the compared lineages. We further modeled parallel evolution over the Arabidopsis candidate genes and showed that a decreasing probability of repeated selection on the same standing or introgressed alleles drives the observed pattern of divergence-dependent parallelism. We therefore conclude that genetic divergence between populations, species, and genera, affecting the pool of shared variants, is an important factor in the predictability of genome evolution.

Population structure limits parallel evolution in sticklebacks

Population genetic theory predicts that small effective population sizes (N_e) and restricted gene flow limit the potential for local adaptation. In particular, the probability of evolving similar phenotypes based on shared genetic mechanisms (i.e., parallel evolution), is expected to be reduced. We tested these predictions in a comparative genomic study of two ecologically similar and geographically codistributed stickleback species (viz. Gasterosteus aculeatus and Pungitius pungitius). We found that P. pungitius harbors less genetic diversity and exhibits higher levels of genetic differentiation and isolation-by-distance than G. aculeatus. Conversely, G. aculeatus exhibits a stronger degree of genetic parallelism across freshwater populations than P. pungitius: 2,996 versus 379 single nucleotide polymorphisms located within 26 versus 9 genomic regions show evidence of selection in multiple freshwater populations of G. aculeatus and P. pungitius, respectively. Most regions involved in parallel evolution in G. aculeatus showed increased levels of divergence, suggestive of selection on ancient haplotypes. In contrast, haplotypes involved in freshwater adaptation in P. pungitius were younger. In accordance with theory, the results suggest that connectivity and genetic drift play crucial roles in determining the levels and geographic distribution of standing genetic variation, providing evidence that population subdivision limits local adaptation and therefore also the likelihood of parallel evolution.

Using Integrative Biology to Infer Adaptation from Comparisons of Two (or a Few) Species

AbstractPhylogenetic comparative methods represent a major advance in integrative and comparative biology and have allowed researchers to rigorously test for adaptation in a macroevolutionary framework. However, phylogenetic comparative methods require trait data for many species, which is impractical for certain taxonomic groups and trait types. We propose that the philosophical principle of severity can be implemented in an integrative framework to generate strong inference of adaptation in studies that compare only a few populations or species. This approach requires (1) ensuring that the study system contains species that are relatively closely related; (2) formulating a specific, clear, overarching hypothesis that can be subjected to integrative testing across levels of biological organization (e.g., ecology, behavior, morphology, physiology, and genetics); (3) collecting data that avoid statistical underdetermination and thus allow severe tests of hypotheses; and (4) systematically refining and refuting alternative hypotheses. Although difficult to collect for more than a few species, detailed, integrative data can be used to differentiate among several potential agents of selection. In this way, integrative studies of small numbers of closely related species can complement and even improve on broadscale phylogenetic comparative studies by revealing the specific drivers of adaptation.

Evolution of silk anchor structure as the joint effect of spinning behavior and spinneret morphology

Spider web anchors are attachment structures composed of the bi-phasic glue-fiber secretion from the piriform silk glands. The mechanical performance of the anchors strongly correlates with the structural assembly of the silk lines, which makes spider silk anchors an ideal system to study the biomechanical function of extended phenotypes and its evolution. It was proposed that silk anchor function guided the evolution of spider web architectures, but its fine-structural variation and whether its evolution was rather determined by changes of the shape of the spinneret tip or in the innate spinning choreography remained unresolved. Here, we comparatively studied the micro-structure of silk anchors across the spider tree of life, and set it in relation to spinneret morphology, spinning behavior and the ecology of the spider. We identified a number of apomorphies in the structure of silk anchors that may positively affect anchor function: 1. bundled dragline, 2. dragline envelope, and 3. dragline suspension ('bridge'). All these characters were apomorphic and evolved repeatedly in multiple lineages, supporting the notion that they are adaptive. The occurrence of these structural features can be explained with changes in the shape and mobility of the spinneret tip, the spinning behavior or both. Spinneret shapes generally varied less than their fine-tuned movements, indicating that changes in construction behavior play a more important role in the evolution of silk anchor assembly. However, the morphology of the spinning apparatus is also a major constraint to the evolution of the spinning choreography. These results highlight changes in behavior as the proximate and in morphology as the ultimate causes of extended phenotype evolution. Further, this research provides a roadmap for future bioprospecting research to design high-performance instant line anchors.

Multigenerational exposure to elevated temperatures leads to a reduction in standard metabolic rate in the wild

In light of global climate change, there is a pressing need to understand and predict the capacity of populations to respond to rising temperatures. Metabolic rate is a key trait that is likely to influence the ability to cope with climate change. Yet, empirical and theoretical work on metabolic rate responses to temperature changes has so far produced mixed results and conflicting predictions.Our study addresses this issue using a novel approach of comparing fish populations in geothermally warmed lakes and adjacent ambient-temperature lakes in Iceland. This unique 'natural experiment' provides repeated and independent examples of populations experiencing contrasting thermal environments for many generations over a small geographic scale, thereby avoiding the confounding factors associated with latitudinal or elevational comparisons. Using Icelandic sticklebacks from three warm and three cold habitats, we measured individual metabolic rates across a range of acclimation temperatures to obtain reaction norms for each population.We found a general pattern for a lower standard metabolic rate (SMR) in sticklebacks from warm habitats when measured at a common temperature, as predicted by Krogh's rule. Metabolic rate differences between warm- and cold-habitat sticklebacks were more pronounced at more extreme acclimation temperatures, suggesting the release of cryptic genetic variation upon exposure to novel conditions, which can reveal hidden evolutionary potential. We also found a stronger divergence in metabolic rate between thermal habitats in allopatry than sympatry, indicating that gene flow may constrain physiological adaptation when dispersal between warm and cold habitats is possible.In sum, our study suggests that fish may diverge toward a lower SMR in a warming world, but this might depend on connectivity and gene flow between different thermal habitats. A free Plain Language Summary can be found within the Supporting Information of this article.

A Similar Genetic Architecture Underlies the Convergent Evolution of the Selfing Syndrome in Capsella

Whether, and to what extent, phenotypic evolution follows predictable genetic paths remains an important question in evolutionary biology. Convergent evolution of similar characters provides a unique opportunity to address this question. The transition to selfing and the associated changes in flower morphology are among the most prominent examples of repeated evolution in plants. In this study, we take advantage of the independent transitions to self-fertilization in the genus Capsella to compare the similarities between parallel modifications of floral traits and test for genetic and developmental constraints imposed on flower evolution in the context of the selfing syndrome. Capsella rubella and Capsella orientalis emerged independently but evolved almost identical flower characters. Not only is the evolutionary outcome identical but the same developmental strategies underlie the convergent reduction of flower size. This has been associated with convergent evolution of gene expression changes. The transcriptomic changes common to both selfing lineages are enriched in genes with low network connectivity and with organ-specific expression patterns. Comparative genetic mapping also suggests that, at least in the case of petal size evolution, these similarities have a similar genetic basis. Based on these results, we hypothesize that the limited availability of low-pleiotropy paths predetermines closely related species to similar evolutionary outcomes.

Phenotypic and Genomic Local Adaptation across Latitude and Altitude in Populus trichocarpa

Local adaptation to climate allows plants to cope with temporally and spatially heterogeneous environments, and parallel phenotypic clines provide a natural experiment to uncover the genomic architecture of adaptation. Though extensive effort has been made to investigate the genomic basis of local adaptation to climate across the latitudinal range of tree species, less is known for altitudinal clines. We used exome capture to genotype 451 Populus trichocarpa genotypes across altitudinal and latitudinal gradients spanning the natural species range, and phenotyped these trees for a variety of adaptive traits in two common gardens. We observed clinal variation in phenotypic traits across the two transects, which indicates climate-driven selection, and coupled gene-based genotype–phenotype and genotype–environment association scans to identify imprints of climatic adaptation on the genome. Although many of the phenotype- and climate-associated genes were unique to one transect, we found evidence of parallelism between latitude and altitude, as well as significant convergence when we compared our outlier genes with those putatively involved in climatic adaptation in two gymnosperm species. These results suggest that not only genomic constraint during adaptation to similar environmental gradients in poplar but also different environmental contexts, spatial scale, and perhaps redundant function among potentially adaptive genes and polymorphisms lead to divergent adaptive architectures.

What can be learned by scanning the genome for molecular convergence in wild populations?

Convergent evolution, where independent lineages evolve similar phenotypes in response to similar challenges, can provide valuable insight into how selection operates and the limitations it encounters. However, it has only recently become possible to explore how convergent evolution is reflected at the genomic level. The overlapping outlier approach (OOA), where genome scans of multiple independent lineages are used to find outliers that overlap and therefore identify convergently evolving loci, is becoming popular. Here, we present a quantitative analysis of 34 studies that used this approach across many sampling designs, taxa, and sampling intensities. We found that OOA studies with increased biological sampling power within replicates have increased likelihood of finding overlapping, "convergent" signals of adaptation between them. When identifying convergent loci as overlapping outliers, it is tempting to assume that any false-positive outliers derived from individual scans will fail to overlap across replicates, but this cannot be guaranteed. We highlight how population demographics and genomic context can contribute toward both true convergence and false positives in OOA studies. We finish with an exploration of emerging methods that couple genome scans with phenotype and environmental measures, leveraging added information from genome data to more directly test hypotheses of the likelihood of convergent evolution.

The shadow of the past: Convergence of young and old South American desert lizards as measured by head shape traits

Convergence is a pervasive phenomenon in the Tree of Life, and evolution of similar phenotypes sharing the same environmental conditions is expected in phylogenetically closely related species. In contrast, contingent factors are probably more influential in shaping phenotypic diversity for distantly related taxa. Here, we test putative convergent evolution of lizard head morphologies among relatively closely related desert dwelling Liolaemus species, and the very distantly related Ctenoblepharys adspersa. We estimated a multilocus time-calibrated phylogeny of 57 species of South American liolaemus lizards, based on seven molecular markers. We collected head shape data for 468 specimens, and used three phylogenetic comparative methods (SURFACE, CONVEVOL, and WHEATSHEAF index) to test for and estimate the strength of convergence. We found strong evidence for convergence among Pacific desert lizard C. adspersa, Liolaemus audivetulatus, Liolaemus insolitus, Liolaemus poconchilensis, Liolaemus stolzmanni, and a candidate species (Liolaemus "Moquegua"). Our results suggest that, despite the long divergence and phylogenetic distance of C. adspersa with respect to convergent Liolaemus species, natural selection was probably more important than historical contingency in shaping phenotypic evolution in these desert lizards.

Contingency and determinism in evolution: Replaying life’s tape

Historical processes display some degree of “contingency,” meaning their outcomes are sensitive to seemingly inconsequential events that can fundamentally change the future. Contingency is what makes historical outcomes unpredictable. Unlike many other natural phenomena, evolution is a historical process. Evolutionary change is often driven by the deterministic force of natural selection, but natural selection works upon variation that arises unpredictably through time by random mutation, and even beneficial mutations can be lost by chance through genetic drift. Moreover, evolution has taken place within a planetary environment with a particular history of its own. This tension between determinism and contingency makes evolutionary biology a kind of hybrid between science and history. While philosophers of science examine the nuances of contingency, biologists have performed many empirical studies of evolutionary repeatability and contingency. Here, we review the experimental and comparative evidence from these studies. Replicate populations in evolutionary “replay” experiments often show parallel changes, especially in overall performance, although idiosyncratic outcomes show that the particulars of a lineage’s history can affect which of several evolutionary paths is taken. Comparative biologists have found many notable examples of convergent adaptation to similar conditions, but quantification of how frequently such convergence occurs is difficult. On balance, the evidence indicates that evolution tends to be surprisingly repeatable among closely related lineages, but disparate outcomes become more likely as the footprint of history grows deeper. Ongoing research on the structure of adaptive landscapes is providing additional insight into the interplay of fate and chance in the evolutionary process.

(Non)Parallel Evolution

Parallel evolution across replicate populations has provided evolutionary biologists with iconic examples of adaptation. When multiple populations colonize seemingly similar habitats, they may evolve similar genes, traits, or functions. Yet, replicated evolution in nature or in the laboratory often yields inconsistent outcomes: Some replicate populations evolve along highly similar trajectories, whereas other replicate populations evolve to different extents or in distinct directions. To understand these heterogeneous outcomes, biologists are increasingly treating parallel evolution not as a binary phenomenon but rather as a quantitative continuum ranging from parallel to nonparallel. By measuring replicate populations’ positions along this (non)parallel continuum, we can test hypotheses about evolutionary and ecological factors that influence the extent of repeatable evolution. We review evidence regarding the manifestation of (non)parallel evolution in the laboratory, in natural populations, and in applied contexts such as cancer. We enumerate the many genetic, ecological, and evolutionary processes that contribute to variation in the extent of parallel evolution.

Genomic divergence between Spanish Littorina saxatilis ecotypes unravels limited admixture and extensive parallelism associated with population history

The rough periwinkle, Littorina saxatilis, is a model system for studying parallel ecological speciation in microparapatry. Phenotypically parallel wave-adapted and crab-adapted ecotypes that hybridize within the middle shore are replicated along the northwestern coast of Spain and have likely arisen from two separate glacial refugia. We tested whether greater geographic separation corresponding to reduced opportunity for contemporary or historical gene flow between parallel ecotypes resulted in less parallel genomic divergence. We sequenced double-digested restriction-associated DNA (ddRAD) libraries from individual snails from upper, mid, and low intertidal levels of three separate sites colonized from two separate refugia. Outlier analysis of 4256 SNP markers identified 34.4% sharing of divergent loci between two geographically close sites; however, these sites each shared only 9.9%-15.1% of their divergent loci with a third more-distant site. STRUCTURE analysis revealed that genotypes from only three of 166 phenotypically intermediate mid-shore individuals appeared to result from recent hybridization, suggesting that hybrids cannot be reliably identified using shell traits. Hierarchical AMOVA indicated that the primary source of genomic differentiation was geographic separation, but also revealed greater similarity of the same ecotype across the two geographically close sites than previously estimated with dominant markers. These results from a model system for ecological speciation suggest that genomic parallelism is affected by the opportunity for historical or contemporary gene flow between populations.

Genomic comparison of two independent seagrass lineages reveals habitat-driven convergent evolution

Seagrasses are marine angiosperms that live fully submerged in the sea. They evolved from land plant ancestors, with multiple species representing at least three independent return-to-the-sea events. This raises the question of whether these marine angiosperms followed the same adaptation pathway to allow them to live and reproduce under the hostile marine conditions. To compare the basis of marine adaptation between seagrass lineages, we generated genomic data for Halophila ovalis and compared this with recently published genomes for two members of Zosteraceae, as well as genomes of five non-marine plant species (Arabidopsis, Oryza sativa, Phoenix dactylifera, Musa acuminata, and Spirodela polyrhiza). Halophila and Zosteraceae represent two independent seagrass lineages separated by around 30 million years. Genes that were lost or conserved in both lineages were identified. All three species lost genes associated with ethylene and terpenoid biosynthesis, and retained genes related to salinity adaptation, such as those for osmoregulation. In contrast, the loss of the NADH dehydrogenase-like complex is unique to H. ovalis. Through comparison of two independent return-to-the-sea events, this study further describes marine adaptation characteristics common to seagrass families, identifies species-specific gene loss, and provides molecular evidence for convergent evolution in seagrass lineages.

Divergent evolution peaks under intermediate population bottlenecks during bacterial experimental evolution

There is growing evidence that parallel molecular evolution is common, but its causes remain poorly understood. Demographic parameters such as population bottlenecks are predicted to be major determinants of parallelism. Here, we test the hypothesis that bottleneck intensity shapes parallel evolution by elucidating the genomic basis of adaptation to antibiotic-supplemented media in hundreds of populations of the bacterium Pseudomonas fluorescens Pf0-1. As expected, bottlenecking decreased the rate of phenotypic and molecular adaptation. Surprisingly, bottlenecking had no impact on the likelihood of parallel adaptive molecular evolution at a genome-wide scale. However, bottlenecking had a profound impact on the genes involved in antibiotic resistance. Specifically, under either intense or weak bottlenecking, resistance predominantly evolved by strongly beneficial mutations which provide high levels of antibiotic resistance. In contrast with intermediate bottlenecking regimes, resistance evolved by a greater diversity of genetic mechanisms, significantly reducing the observed levels of parallel genetic evolution. Our results demonstrate that population bottlenecking can be a major predictor of parallel evolution, but precisely how may be more complex than many simple theoretical predictions.

Approaches to Macroevolution: 1. General Concepts and Origin of Variation

Approaches to macroevolution require integration of its two fundamental components, i.e. the origin and the sorting of variation, in a hierarchical framework. Macroevolution occurs in multiple currencies that are only loosely correlated, notably taxonomic diversity, morphological disparity, and functional variety. The origin of variation within this conceptual framework is increasingly understood in developmental terms, with the semi-hierarchical structure of gene regulatory networks (GRNs, used here in a broad sense incorporating not just the genetic circuitry per se but the factors controlling the timing and location of gene expression and repression), the non-linear relation between magnitude of genetic change and the phenotypic results, the evolutionary potential of co-opting existing GRNs, and developmental responsiveness to nongenetic signals (i.e. epigenetics and plasticity), all requiring modification of standard microevolutionary models, and rendering difficult any simple definition of evolutionary novelty. The developmental factors underlying macroevolution create anisotropic probabilities-i.e., an uneven density distribution-of evolutionary change around any given phenotypic starting point, and the potential for coordinated changes among traits that can accommodate change via epigenetic mechanisms. From this standpoint, "punctuated equilibrium" and "phyletic gradualism" simply represent two cells in a matrix of evolutionary models of phenotypic change, and the origin of trends and evolutionary novelty are not simply functions of ecological opportunity. Over long timescales, contingency becomes especially important, and can be viewed in terms of macroevolutionary lags (the temporal separation between the origin of a trait or clade and subsequent diversification); such lags can arise by several mechanisms: as geological or phylogenetic artifacts, or when diversifications require synergistic interactions among traits, or between traits and external events. The temporal and spatial patterns of the origins of evolutionary novelties are a challenge to macroevolutionary theory; individual events can be described retrospectively, but a general model relating development, genetics, and ecology is needed. An accompanying paper (Jablonski in Evol Biol 2017) reviews diversity dynamics and the sorting of variation, with some general conclusions.

Deep-Time Convergence in Rove Beetle Symbionts of Army Ants

Recent adaptive radiations provide striking examples of convergence [1-4], but the predictability of evolution over much deeper timescales is controversial, with a scarcity of ancient clades exhibiting repetitive patterns of phenotypic evolution [5, 6]. Army ants are ecologically dominant arthropod predators of the world's tropics, with large nomadic colonies housing diverse communities of socially parasitic myrmecophiles [7]. Remarkable among these are many species of rove beetle (Staphylinidae) that exhibit ant-mimicking "myrmecoid" body forms and are behaviorally accepted into their aggressive hosts' societies: emigrating with colonies and inhabiting temporary nest bivouacs, grooming and feeding with workers, but also consuming the brood [8-11]. Here, we demonstrate that myrmecoid rove beetles are strongly polyphyletic, with this adaptive morphological and behavioral syndrome having evolved at least 12 times during the evolution of a single staphylinid subfamily, Aleocharinae. Each independent myrmecoid clade is restricted to one zoogeographic region and highly host specific on a single army ant genus. Dating estimates reveal that myrmecoid clades are separated by substantial phylogenetic distances-as much as 105 million years. All such groups arose in parallel during the Cenozoic, when army ants diversified into modern genera [12] and rose to ecological dominance [13, 14]. This work uncovers a rare example of an ancient system of complex morphological and behavioral convergence, with replicate beetle lineages following a predictable phenotypic trajectory during their parasitic adaptation to host colonies.

Getting a head in hard soils: Convergent skull evolution and divergent allometric patterns explain shape variation in a highly diverse genus of pocket gophers (Thomomys)

Background

High morphological diversity can occur in closely related animals when selection favors morphologies that are subject to intrinsic biological constraints. A good example is subterranean rodents of the genus Thomomys, one of the most taxonomically and morphologically diverse mammalian genera. Highly procumbent, tooth-digging rodent skull shapes are often geometric consequences of increased body size. Indeed, larger-bodied Thomomys species tend to inhabit harder soils. We used geometric morphometric analyses to investigate the interplay between soil hardness (the main extrinsic selection pressure on fossorial mammals) and allometry (i.e. shape change due to size change; generally considered the main intrinsic factor) on crania and humeri in this fast-evolving mammalian clade.

Results

Larger Thomomys species/subspecies tend to have more procumbent cranial shapes with some exceptions, including a small-bodied species inhabiting hard soils. Counter to earlier suggestions, cranial shape within Thomomys does not follow a genus-wide allometric pattern as even regional subpopulations differ in allometric slopes. In contrast, humeral shape varies less with body size and with soil hardness. Soft-soil taxa have larger humeral muscle attachment sites but retain an orthodont (non-procumbent) cranial morphology. In intermediate soils, two pairs of sister taxa diverge through differential modifications on either the humerus or the cranium. In the hardest soils, both humeral and cranial morphology are derived through large muscle attachment sites and a high degree of procumbency.

Conclusions

Our results show that conflict between morphological function and intrinsic allometric patterning can quickly and differentially alter the rodent skeleton, especially the skull. In addition, we found a new case of convergent evolution of incisor procumbency among large-, medium-, and small-sized species inhabiting hard soils. This occurs through different combinations of allometric and non-allometric changes, contributing to shape diversity within the genus. The strong influence of allometry on cranial shape appears to confirm suggestions that developmental change underlies mammalian cranial shape divergences, but this requires confirmation from ontogenetic studies. Our findings illustrate how a variety of intrinsic processes, resulting in species-level convergence, could sustain a genus-level range across a variety of extrinsic environments. This might represent a mechanism for observations of genus-level niche conservation despite species extinctions in mammals.

Electronic supplementary material

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

Evolutionary assembly patterns of prokaryotic genomes

Evolutionary innovation must occur in the context of some genomic background, which limits available evolutionary paths. For example, protein evolution by sequence substitution is constrained by epistasis between residues. In prokaryotes, evolutionary innovation frequently happens by macrogenomic events such as horizontal gene transfer (HGT). Previous work has suggested that HGT can be influenced by ancestral genomic content, yet the extent of such gene-level constraints has not yet been systematically characterized. Here, we evaluated the evolutionary impact of such constraints in prokaryotes, using probabilistic ancestral reconstructions from 634 extant prokaryotic genomes and a novel framework for detecting evolutionary constraints on HGT events. We identified 8228 directional dependencies between genes and demonstrated that many such dependencies reflect known functional relationships, including for example, evolutionary dependencies of the photosynthetic enzyme RuBisCO. Modeling all dependencies as a network, we adapted an approach from graph theory to establish chronological precedence in the acquisition of different genomic functions. Specifically, we demonstrated that specific functions tend to be gained sequentially, suggesting that evolution in prokaryotes is governed by functional assembly patterns. Finally, we showed that these dependencies are universal rather than clade-specific and are often sufficient for predicting whether or not a given ancestral genome will acquire specific genes. Combined, our results indicate that evolutionary innovation via HGT is profoundly constrained by epistasis and historical contingency, similar to the evolution of proteins and phenotypic characters, and suggest that the emergence of specific metabolic and pathological phenotypes in prokaryotes can be predictable from current genomes.