Parallel molecular evolution in an herbivore community

No detectable effect of urbanization on genetic drift or gene flow in specialist herbivorous insects of milkweed

Urbanization is hypothesized to isolate populations and restrict dispersal, leading to reduced genetic diversity and increased genetic differentiation. We tested this hypothesis in specialist herbivorous insects of milkweed, positing that higher dispersal ability would mitigate the negative effects of urbanization on genetic drift and gene flow, and that these effects would vary with city size. In this study, we collected 383 milkweed insects from urban and rural sites in Toronto, Canada, and five surrounding cities. Using ddRADseq, we generated 145,000 SPNs for monarchs, 10,000 SNPs for beetles, 6,000 SNPs for weevils to quantify genetic diversity, demographic history and population genetic structure. Contrary to our hypotheses, our results indicated no effect of urbanization or dispersal ability on diversity or genetic differentiation. Genetic diversity, measured as π, varied between 0.0013 and 0.0044 across species, with no urban vs. rural component, but with monarchs having >2 X higher diversity compared to beetles and weevils. Similarly, genetic differentiation was generally low, FST varying between 0.01 and 0.28, but there are no consistent trends among urban vs. rural samples for any of the three species. However, demographic analyses revealed a consistent decline in effective population size for all three sampled species, beginning around the last glacial maximum and intensifying over the past 1,000 years. Our findings suggest that both urbanization and dispersal ability have not been a major factor in reducing gene flow or increasing genetic drift among milkweed's herbivorous insect populations. Instead, historical events such as climatic change since the last glacial maximum, and large-scale anthropogenic disturbance in general, have had a more pronounced impact on demography. These results highlight the importance of considering the combined effects of natural and anthropogenic long-term historical processes when studying population genetics in the context of urbanization.

Entomopathogenic Nematode Species Vary in Their Behavior and Virulence in Response to Cardiac Glycosides Within and Around Insect Hosts

Plants produce defensive toxins to deter herbivores. In response, some specialized herbivores evolved resistance and even the capacity to sequester toxins, affecting interactions at higher trophic levels. Here, we test the hypothesis that potential natural enemies of specialized herbivores are differentially affected by plant toxins depending on their level of adaptation to the plant-herbivore system. We focus on toxic cardiac glycosides (CGs) from milkweeds (Asclepias spp.), which inhibit animal Na+/K+-ATPases, and two CG-resistant insects, the large milkweed bug Oncopeltus fasciatus and a CRISPR-edited Drosophila melanogaster. Both have CG-resistant Na+/K+-ATPases through a set of key amino acid substitutions, which facilitate CG sequestration. We conducted infection experiments with entomopathogenic nematodes (Steinernema carpocapsae, S. feltiae, and S. hermaphroditum) as natural enemies on host insects containing mixtures of milkweed-derived CGs or purified CGs (ouabain, digoxin, and digitoxin) that vary in toxicity. The nematode S. carpocapsae is known to occur in soil near milkweed plants and naturally has several of the same Na+/K+-ATPase substitutions as the milkweed bug O. fasciatus and our Drosophila mutant. This nematode not only exhibited higher fecundity in hosts that carried CGs relative to the other nematode species (which have sensitive Na+/K+-ATPases), but also showed attraction to mixtures of CGs in milkweed root extracts and to purified ouabain when tested on agar plates. A coiling phenotype, which is a symptom of neurotoxicity, was observed more frequently in S. feltiae and S. hermaphroditum upon exposure to milkweed root extracts than in S. carpocapsae. Nematode behavior was further tested in sand, and while attraction to CGs was found for S. carpocapsae, nematodes of the other species tended to migrate away from milkweed root chemicals. Thus, S. carpocapsae can tolerate CGs and may use these as chemical cues to locate insect hosts that live on or around milkweed plants.

Recurrent evolution of small body size and loss of the sword ornament in Northern swordtail fish

Across the tree of life, species have repeatedly evolved similar phenotypes. While well-studied for ecological traits, there is also evidence for recurrent evolution of sexually selected traits. Swordtail fish (Xiphophorus) is a classic model system for studying sexual selection, and female Xiphophorus exhibit strong mate preferences for large male body sizes and a range of sexually dimorphic ornaments. Interestingly, sexually selected traits have also been lost multiple times in the genus. However, there has been uncertainty over the number of losses of ornamentation and large body size because phylogenetic relationships between species in this group have historically been controversial, partially due to prevalent gene flow. Here, we use whole-genome sequencing approaches to reexamine phylogenetic relationships within a Xiphophorus clade that varies in the presence and absence of sexually selected traits. Using wild-caught individuals, we determine the phylogenetic placement of a small, unornamented species, X. continens, confirming an additional loss of ornamentation and large body size in the clade. With these revised phylogenetic relationships, we analyze evidence for coevolution between body size and other sexually selected traits using phylogenetic comparative methods. These results provide insights into the evolutionary pressures driving the recurrent loss of suites of sexually selected traits.

Substrate Specificity of ABCB Transporters Predicted by Docking Simulations Can Be Confirmed by Experimental Tests

ATP-binding cassette (ABC) transporters, particularly those of subfamily B, are involved in cell detoxification, multidrug resistance, drug treatment pharmacodynamics, and also ecological adaptation. In this regard, ABCB transporters may play a decisive role in the co-evolution between plants and herbivores. Cardenolides, toxic steroid glycosides, are secondary plant metabolites that defend plants against herbivores by targeting their sodium-potassium ATPase. Despite their toxicity, several herbivorous insects such as the large milkweed bug (Oncopeltus fasciatus) have evolved adaptations to tolerate cardenolides and sequester them for their own defense. We investigate the role of two ABCB transporters of O. fasciatus for the paracellular transport of cardenolides by docking simulations and ATPase assays. Cardenolide binding of OfABCB1 and OfABCB2 is predicted by docking simulations and calculated binding energies are compared with substrate specificities determined in ATPase assays. Both tested ABCB transporters showed activity upon exposure to cardenolides and Km values that agreed well with the predictions of our docking simulations. We conclude that docking simulations can help identify transporter binding regions and predict substrate specificity, as well as provide deeper insights into the structural basis of ABC transporter function.

Convergence and Divergence among Herbivorous Insects Specialized on Toxic Plants: Revealing Syndromes among the Cardenolide Feeders across the Insect Tree of Life

AbstractRepeatable macroevolutionary patterns provide hope for rules in biology, especially when we can decipher the underlying mechanisms. Here we synthesize natural history, genetic adaptations, and toxin sequestration in herbivorous insects that specialize on plants with cardiac glycoside defenses. Work on the monarch butterfly provided a model for evolution of the "sequestering specialist syndrome," where specific amino acid substitutions in the insect's Na+/K+-ATPase are associated with (1) high toxin resistance (target site insensitivity [TSI]), (2) sequestration of toxins, and (3) aposematic coloration. We evaluate convergence for these traits within and between Lepidoptera, Coleoptera, Diptera, Hemiptera, Hymenoptera, and Orthoptera, encompassing hundreds of toxin-adapted species. Using new and existing data on ∼28 origins of specialization, we show that the monarch model evolved independently in five taxonomic orders (but not Diptera). An additional syndrome occurs in five orders (all but Hymenoptera): aposematic sequesterers with modest to medium TSI. Indeed, all sequestering species were aposematic, and all but one had at least modest TSI. Additionally, several species were aposematic nonsequesterers (potential Batesian mimics), and this combination evolved in species with a range of TSI levels. Finally, we identified some biases among these strategies within taxonomic orders. Biodiversity in this microcosm of life evolved repeatedly with a high degree of similarity across six taxonomic orders, yet we identified alternative trait combinations as well as lineage-specific outcomes.

Escalation by duplication: Milkweed bug trumps Monarch butterfly

The iconic Monarch butterfly is probably the best-known example of chemical defence against predation, as pictures of vomiting naive blue jays in countless textbooks vividly illustrate. Larvae of the butterfly take up toxic cardiac glycosides from their milkweed hostplants and carry them over to the adult stage. These compounds (cardiotonic steroids, including cardenolides and bufadienolides) inhibit the animal transmembrane sodium-potassium ATPase (Na,K-ATPase), but the Monarch enzyme resists this inhibition thanks to amino acid substitutions in its catalytic alpha-subunit. Some birds also have substitutions and can feast on cardiac glycoside-sequestering insects with impunity. A flurry of recent work has shown how the alpha-subunit gene has been duplicated multiple times in separate insect lineages specializing in cardiac glycoside-producing plants. In this issue of Molecular Ecology, Herbertz et al. toss the beta-subunit into the mix, by expressing all nine combinations of three alpha- and three beta-subunits of the milkweed bug Na,K-ATPase and testing their response to a cardenolide from the hostplant. The findings suggest that the diversification and subfunctionalization of genes allow milkweed bugs to balance trade-offs between resistance towards sequestered host plant toxins that protect the bugs from predators, and physiological costs in terms of Na,K-ATPase activity.

Coevolutionary escalation led to differentially adapted paralogs of an insect's Na,K-ATPase optimizing resistance to host plant toxins

Cardiac glycosides are chemical defence toxins known to fatally inhibit the Na,K-ATPase (NKA) throughout the animal kingdom. Several animals, however, have evolved target-site insensitivity through substitutions in the otherwise highly conserved cardiac glycoside binding pocket of the NKA. The large milkweed bug, Oncopeltus fasciatus, shares a long evolutionary history with cardiac glycoside containing plants that led to intricate adaptations. Most strikingly, several duplications of the bugs' NKA1α gene provided the opportunity for differential resistance-conferring substitutions and subsequent sub-functionalization of the enzymes. Here, we analysed cardiac glycoside resistance and ion pumping activity of nine functional NKA α/β-combinations of O. fasciatus expressed in cell culture. We tested the enzymes with two structurally distinct cardiac glycosides, calotropin, a host plant compound, and ouabain, a standard cardiac glycoside. The identity and number of known resistance-conferring substitutions in the cardiac glycoside binding site significantly impacted activity and toxin resistance in the three α-subunits. The β-subunits also influenced the enzymes' characteristics, yet to a lesser extent. Enzymes containing the more ancient αC-subunit were inhibited by both compounds but much more strongly by the host plant toxin calotropin than by ouabain. The sensitivity to calotropin was diminished in enzymes containing the more derived αB and αA, which were only marginally inhibited by both cardiac glycosides. This trend culminated in αAβ1 having higher resistance against calotropin than against ouabain. These results support the coevolutionary escalation of plant defences and herbivore tolerance mechanisms. The possession of multiple paralogs additionally mitigates pleiotropic effects by compromising between ion pumping activity and resistance.

Increased Positive Selection in Highly Recombining Genes Does not Necessarily Reflect an Evolutionary Advantage of Recombination

It is commonly thought that the long-term advantage of meiotic recombination is to dissipate genetic linkage, allowing natural selection to act independently on different loci. It is thus theoretically expected that genes with higher recombination rates evolve under more effective selection. On the other hand, recombination is often associated with GC-biased gene conversion (gBGC), which theoretically interferes with selection by promoting the fixation of deleterious GC alleles. To test these predictions, several studies assessed whether selection was more effective in highly recombining genes (due to dissipation of genetic linkage) or less effective (due to gBGC), assuming a fixed distribution of fitness effects (DFE) for all genes. In this study, I directly derive the DFE from a gene's evolutionary history (shaped by mutation, selection, drift, and gBGC) under empirical fitness landscapes. I show that genes that have experienced high levels of gBGC are less fit and thus have more opportunities for beneficial mutations. Only a small decrease in the genome-wide intensity of gBGC leads to the fixation of these beneficial mutations, particularly in highly recombining genes. This results in increased positive selection in highly recombining genes that is not caused by more effective selection. Additionally, I show that the death of a recombination hotspot can lead to a higher dN/dS than its birth, but with substitution patterns biased towards AT, and only at selected positions. This shows that controlling for a substitution bias towards GC is therefore not sufficient to rule out the contribution of gBGC to signatures of accelerated evolution. Finally, although gBGC does not affect the fixation probability of GC-conservative mutations, I show that by altering the DFE, gBGC can also significantly affect nonsynonymous GC-conservative substitution patterns.

Firefly toxin lucibufagins evolved after the origin of bioluminescence

Fireflies were believed to originally evolve their novel bioluminescence as warning signals to advertise their toxicity to predators, which was later adopted in adult mating. Although the evolution of bioluminescence has been investigated extensively, the warning signal hypothesis of its origin has not been tested. In this study, we test this hypothesis by systematically determining the presence or absence of firefly toxin lucibufagins (LBGs) across firefly species and inferring the time of origin of LBGs. We confirm the presence of LBGs in the subfamily Lampyrinae, but more importantly, we reveal the absence of LBGs in other lineages, including the subfamilies of Luciolinae, Ototretinae, and Psilocladinae, two incertae sedis lineages, and the Rhagophthalmidae family. Ancestral state reconstructions for LBGs based on firefly phylogeny constructed using genomic data suggest that the presence of LBGs in the common ancestor of the Lampyrinae subfamily is highly supported but unsupported in more ancient nodes, including firefly common ancestors. Our results suggest that firefly LBGs probably evolved much later than the evolution of bioluminescence. We thus conclude that firefly bioluminescence did not originally evolve as direct warning signals for toxic LBGs and advise that future studies should focus on other hypotheses. Moreover, LBG toxins are known to directly target and inhibit the α subunit of Na+, K+-ATPase (ATPα). We further examine the effects of amino acid substitutions in firefly ATPα on its interactions with LBGs. We find that ATPα in LBG-containing fireflies is relatively insensitive to LBGs, which suggests that target-site insensitivity contributes to LBG-containing fireflies' ability to deal with their own toxins.

Comparative transcriptomics revealed parallel evolution and innovation of photosymbiosis molecular mechanisms in a marine bivalve

Photosymbioses between heterotrophic hosts and autotrophic symbionts are evolutionarily prevalent and ecologically significant. However, the molecular mechanisms behind such symbioses remain less elucidated, which hinders our understanding of their origin and adaptive evolution. This study compared gene expression patterns in a photosymbiotic bivalve (Fragum sueziense) and a closely related non-symbiotic species (Trigoniocardia granifera) under different light conditions to detect potential molecular pathways involved in mollusc photosymbiosis. We discovered that the presence of algal symbionts greatly impacted host gene expression in symbiont-containing tissues. We found that the host immune functions were suppressed under normal light compared with those in the dark. In addition, we found that cilia in the symbiont-containing tissues play important roles in symbiont regulation or photoreception. Interestingly, many potential photosymbiosis genes could not be annotated or do not exhibit orthologues in T. granifera transcriptomes, indicating unique molecular functions in photosymbiotic bivalves. Overall, we found both novel and known molecular mechanisms involved in animal-algal photosymbiosis within bivalves. Given that many of the molecular pathways are shared among distantly related host lineages, such as molluscs and cnidarians, it indicates that parallel and/or convergent evolution is instrumental in shaping host-symbiont interactions and responses in these organisms.

Emx2 underlies the development and evolution of marsupial gliding membranes

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

Accurate Detection of Convergent Mutations in Large Protein Alignments With ConDor

Evolutionary convergences are observed at all levels, from phenotype to DNA and protein sequences, and changes at these different levels tend to be correlated. Notably, convergent mutations can lead to convergent changes in phenotype, such as changes in metabolism, drug resistance, and other adaptations to changing environments. We propose a two-component approach to detect mutations subject to convergent evolution in protein alignments. The "Emergence" component selects mutations that emerge more often than expected, while the "Correlation" component selects mutations that correlate with the convergent phenotype under study. With regard to Emergence, a phylogeny deduced from the alignment is provided by the user and is used to simulate the evolution of each alignment position. These simulations allow us to estimate the expected number of mutations in a neutral model, which is compared to the observed number of mutations in the data studied. In Correlation, a comparative phylogenetic approach, is used to measure whether the presence of each of the observed mutations is correlated with the convergent phenotype. Each component can be used on its own, for example Emergence when no phenotype is available. Our method is implemented in a standalone workflow and a webserver, called ConDor. We evaluate the properties of ConDor using simulated data, and we apply it to three real datasets: sedge PEPC proteins, HIV reverse transcriptase, and fish rhodopsin. The results show that the two components of ConDor complement each other, with an overall accuracy that compares favorably to other available tools, especially on large datasets.

Fecal Deployment: An Alternative Way of Defensive Host Plant Cardenolide Use by Lilioceris merdigera Larvae

The brilliant red Lilioceris merdigera (Coleoptera, Chrysomelidae) can spend its entire life cycle on the cardenolide-containing plant Convallaria majalis (lily of the valley) and forms stable populations on this host. Yet, in contrast to many other insects on cardenolide-containing plants L. merdigera does not sequester these plant toxins in the body but rather both adult beetles and larvae eliminate ingested cardenolides with the feces. Tracer feeding experiments showed that this holds true for radioactively labeled ouabain and digoxin, a highly polar and a rather apolar cardenolide. Both compounds or their derivatives are incorporated in the fecal shields of the larvae. The apolar digoxin, but not the polar ouabain, showed a deterrent effect on the generalist predatory ant Myrmica rubra, which occurs in the habitat of L. merdigera. The deterrent effect was detected for digoxin both in choice and feeding time assays. In a predator choice assay, a fecal shield derived from a diet of cardenolide-containing C. majalis offered L. merdigera larvae better protection from M. rubra than one derived from non-cardenolide Allium schoenoprasum (chives) or no fecal shield at all. Thus, we here present data suggesting a new way how insects may gain protection by feeding on cardenolide-containing plants.

Multiple Origins of Bioluminescence in Beetles and Evolution of Luciferase Function

Bioluminescence in beetles has long fascinated biologists, with diverse applications in biotechnology. To date, however, our understanding of its evolutionary origin and functional variation mechanisms remains poor. To address these questions, we obtained high-quality reference genomes of luminous and nonluminous beetles in 6 Elateroidea families. We then reconstructed a robust phylogenetic relationship for all luminous families and related nonluminous families. Comparative genomic analyses and biochemical functional experiments suggested that gene evolution within Elateroidea played a crucial role in the origin of bioluminescence, with multiple parallel origins observed in the luminous beetle families. While most luciferase-like proteins exhibited a conserved nonluminous amino acid pattern (TLA346 to 348) in the luciferin-binding sites, luciferases in the different luminous beetle families showed divergent luminous patterns at these sites (TSA/CCA/CSA/LVA). Comparisons of the structural and enzymatic properties of ancestral, extant, and site-directed mutant luciferases further reinforced the important role of these sites in the trade-off between acyl-CoA synthetase and luciferase activities. Furthermore, the evolution of bioluminescent color demonstrated a tendency toward hypsochromic shifts and variations among the luminous families. Taken together, our results revealed multiple parallel origins of bioluminescence and functional divergence within the beetle bioluminescent system.

Predatory fireflies and their toxic firefly prey have evolved distinct toxin resistance strategies

Toxic cardiotonic steroids (CTSs) act as a defense mechanism in many firefly species (Lampyridae) by inhibiting a crucial enzyme called Na+,K+-ATPase (NKA). Although most fireflies produce these toxins internally, species of the genus Photuris acquire them from a surprising source: predation on other fireflies. The contrasting physiology of toxin exposure and sequestration between Photuris and other firefly genera suggests that distinct strategies may be required to prevent self-intoxication. Our study demonstrates that both Photuris and their firefly prey have evolved highly resistant NKAs. Using an evolutionary analysis of the specific target of CTS (ATPα) in fireflies and gene editing in Drosophila, we find that the initial steps toward resistance were shared among Photuris and other firefly lineages. However, the Photuris lineage subsequently underwent multiple rounds of gene duplication and neofunctionalization, resulting in the development of ATPα paralogs that are differentially expressed and exhibit increasing resistance to CTS. By contrast, other firefly species have maintained a single copy. Our results implicate gene duplication as a facilitator in the transition of Photuris to its distinct ecological role as a predator of toxic firefly prey.

Urbanization and a green corridor do not impact genetic divergence in common milkweed (Asclepias syriaca L.)

Urbanization is altering landscapes globally at an unprecedented rate. While ecological differences between urban and rural environments often promote phenotypic divergence among populations, it is unclear to what degree these trait differences arise from genetic divergence as opposed to phenotypic plasticity. Furthermore, little is known about how specific landscape elements, such as green corridors, impact genetic divergence in urban environments. We tested the hypotheses that: (1) urbanization, and (2) proximity to an urban green corridor influence genetic divergence in common milkweed (Asclepias syriaca) populations for phenotypic traits. Using seeds from 52 populations along three urban-to-rural subtransects in the Greater Toronto Area, Canada, one of which followed a green corridor, we grew ~ 1000 plants in a common garden setup and measured > 20 ecologically-important traits associated with plant defense/damage, reproduction, and growth over four years. We found significant heritable variation for nine traits within common milkweed populations and weak phenotypic divergence among populations. However, neither urbanization nor an urban green corridor influenced genetic divergence in individual traits or multivariate phenotype. These findings contrast with the expanding literature demonstrating that urbanization promotes rapid evolutionary change and offer preliminary insights into the eco-evolutionary role of green corridors in urban environments.

Predatory fireflies and their toxic firefly prey have evolved distinct toxin resistance strategies

Toxic cardiotonic steroids (CTS) act as a defense mechanism in many firefly species (Lampyridae) by inhibiting a crucial enzyme called Na+,K+-ATPase (NKA). While most fireflies produce these toxins internally, species of the genus Photuris acquire them from a surprising source: predation on other fireflies. The contrasting physiology of toxin exposure and sequestration between Photuris and other firefly genera suggests that distinct strategies may be required to prevent self-intoxication. Our study demonstrates that both Photuris and their firefly prey have evolved highly-resistant NKAs. Using an evolutionary analysis of the specific target of CTS (ATPα) in fireflies, and gene-editing in Drosophila, we find that the initial steps towards resistance were shared among Photuris and other firefly lineages. However, the Photuris lineage subsequently underwent multiple rounds of gene duplication and neofunctionalization, resulting in the development of ATPα paralogs that are differentially expressed and exhibit increasing resistance to CTS. In contrast, other firefly species have maintained a single copy. Our results implicate gene duplication as a facilitator in the transition of Photuris to its distinct ecological role as predator of toxic firefly prey.

Understanding natural selection and similarity: Convergent, parallel and repeated evolution

Parallel and convergent evolution offer some of the most compelling evidence for the significance of natural selection in evolution, as the emergence of similar adaptive solutions is unlikely to occur by random chance alone. However, these terms are often employed inconsistently, leading to misinterpretation and confusion, and recently proposed definitions have unintentionally diminished the emphasis on the evolution of similar adaptive solutions. Here, I examine various conceptual frameworks and definitions related to parallel and convergent evolution and propose a consolidated framework that enhances our comprehension of these evolutionary patterns. The primary aim of this framework is to harmonize the concepts of parallel and convergent evolution together with natural selection and the idea of similarity. Both concepts involve the evolution of similar adaptive solutions as a result of environmental challenges. The distinction lies in ancestral phenotypes. Parallel evolution takes place when the ancestral phenotypes (before selection) of the lineages are similar. Convergent evolution happens when the lineages have distinct ancestral phenotypes (before selection). Because an ancestral-based distinction will inevitably lead to cases where uncertainty in the distinction may arise, the framework includes a general term, repeated evolution, which can be used as a term applying to the evolution of similar phenotypes and genotypes as well as similar responses to environmental pressures. Based on the argument that genetic similarity may frequently arise without selection, the framework posits that the similarity of genetic sequences is not of great interest unless linked to the actions of natural selection or to the origins (mutation, standing genetic variation, gene flow) and locations of the similar sequences.

Recent parallel speciation in Antirrhinum involved complex haplotypes and multiple adaptive characters

A role of ecological adaptation in speciation can be obscured by stochastic processes and differences that species accumulate after genetic isolation. One way to identify adaptive characters and their underlying genes is to study cases of speciation involving parallel adaptations. Recently resolved phylogenies reveal that alpine morphology has evolved in parallel in the genus Antirrhinum (snapdragons): first in an early split of an alpine from a lowland lineage and, more recently, from within the lowland lineage to produce closely related sympatric species with contrasting alpine and lowland forms. Here, we find that two of these later diverged sympatric species are differentiated by only around 2% of nuclear loci. Though showing evidence of recent gene flow, the species remain distinct for a suite of morphological characters typical of earlier-diverged alpine or lowland lineages and their morphologies correlate with features of the local landscape, as expected of ecological adaptations. Morphological differences between the two species involve multiple, unlinked genes so that parental character combinations are readily broken up by recombination in hybrids. We detect little evidence for post-pollination barriers to gene flow or recombination, suggesting that genetic isolation related to ecological adaptation is important in maintaining character combinations and might have contributed to parallel speciation. We also find evidence that genes involved in the earlier alpine-lowland split were reused in parallel evolution of alpine species, consistent with introgressive hybridisation, and speculate that many non-ecological barriers to gene flow might have been purged during the process.

Evolution of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae)

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

Mutation bias and the predictability of evolution

Predicting evolutionary outcomes is an important research goal in a diversity of contexts. The focus of evolutionary forecasting is usually on adaptive processes, and efforts to improve prediction typically focus on selection. However, adaptive processes often rely on new mutations, which can be strongly influenced by predictable biases in mutation. Here, we provide an overview of existing theory and evidence for such mutation-biased adaptation and consider the implications of these results for the problem of prediction, in regard to topics such as the evolution of infectious diseases, resistance to biochemical agents, as well as cancer and other kinds of somatic evolution. We argue that empirical knowledge of mutational biases is likely to improve in the near future, and that this knowledge is readily applicable to the challenges of short-term prediction. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Knockdown of Na,K-ATPase β-subunits in Oncopeltus fasciatus induces molting problems and alterations in tracheal morphology

The ubiquitously expressed transmembrane enzyme Na,K-ATPase (NKA) is vital in maintaining functionality of cells. The association of α- and β-subunits is believed to be essential for forming a functional enzyme. In the large milkweed bug Oncopeltus fasciatus four α1-paralogs and four β-subunits exist that can associate into NKA complexes. This diversity raises the question of possible tissue-specific distribution and function. While the α1-subunits are known to modulate cardenolide-resistance and ion-transport efficiency, the functional importance of the β-subunits needed further investigation. We here characterize all four different β-subunits at the cellular, tissue, and whole organismal scales. A knockdown of different β-subunits heavily interferes with molting success resulting in strongly hampered phenotypes. The failure of ecdysis might be related to disrupted septate junction (SJ) formation, also reflected in β2-suppression-induced alteration in tracheal morphology. Our data further suggest the existence of isolated β-subunits forming homomeric or β-heteromeric complexes. This possible standalone and structure-specific distribution of the β-subunits predicts further, yet unknown pump-independent functions. The different effects caused by β knockdowns highlight the importance of the various β-subunits to fulfill tissue-specific requirements.

Polygenic routes lead to parallel altitudinal adaptation in Heliosperma pusillum (Caryophyllaceae)

Understanding how organisms adapt to the environment is a major goal of modern biology. Parallel evolution-the independent evolution of similar phenotypes in different populations-provides a powerful framework to investigate the evolutionary potential of populations, the constraints of evolution, its repeatability and therefore its predictability. Here, we quantified the degree of gene expression and functional parallelism across replicated ecotype formation in Heliosperma pusillum (Caryophyllaceae), and gained insights into the architecture of adaptive traits. Population structure analyses and demographic modelling support a previously formulated hypothesis of parallel polytopic divergence of montane and alpine ecotypes. We detect a large proportion of differentially expressed genes (DEGs) underlying divergence within each replicate ecotype pair, with a strikingly low number of shared DEGs across pairs. Functional enrichment of DEGs reveals that the traits affected by significant expression divergence are largely consistent across ecotype pairs, in strong contrast to the nonshared genetic basis. The remarkable redundancy of differential gene expression indicates a polygenic architecture for the diverged adaptive traits. We conclude that polygenic traits appear key to opening multiple routes for adaptation, widening the adaptive potential of organisms.

Genomic convergence underlying high-altitude adaptation in alpine plants

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

Identification of the NA+/K+-ATPase α-Isoforms in Six Species of Poison Dart Frogs and their Sensitivity to Cardiotonic Steroids

Cardiotonic steroids (CTS) are a group of compounds known to be toxic due to their ability to inhibit the Na⁺/K⁺-ATPase (NKA), which is essential to maintain the balance of ions in animal cells. An evolutionary strategy of molecular adaptation to avoid self-intoxication acquired by CTS defended organisms and their predators is the structural modification of their NKA where specific amino acid substitutions confer resistant phenotypes. Several lineages of poison dart frogs (Dendrobatidae) are well known to sequester a wide variety of lipophilic alkaloids from their arthropod diet, however there is no evidence of CTS-sequestration or dietary exposure. Interestingly this study identified the presence of α-NKA isoforms (α₁ and α₂) with amino acid substitutions indicative of CTS-resistant phenotypes in skeletal muscle transcriptomes obtained from six species of dendrobatids: Phyllobates aurotaenia, Oophaga anchicayensis, Epipedobates boulengeri, Andinobates bombetes, Andinobates minutus, and Leucostethus brachistriatus, collected in the Valle del Cauca (Colombia). P. aurotaenia, A. minutus, and E. boulengeri presented two variants for α₁-NKA, with one of them having these substitutions. In contrast, O. anchicayensis and A. bombetes have only one α₁-NKA isoform with an amino acid sequence indicative of CTS susceptibility and an α₂-NKA with one substitution that could confer a reduced affinity for CTS. The α₁ and α₂ isoforms of L. brachistriatus do not contain substitutions imparting CTS resistance. Our findings indicate that poison dart frogs express α-NKA isoforms with different affinities for CTS and the pattern of this expression might be influenced by factors related to evolutionary, physiological, ecological, and geographical burdens.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Detecting macroevolutionary genotype-phenotype associations using error-corrected rates of protein convergence

On macroevolutionary timescales, extensive mutations and phylogenetic uncertainty mask the signals of genotype-phenotype associations underlying convergent evolution. To overcome this problem, we extended the widely used framework of non-synonymous to synonymous substitution rate ratios and developed the novel metric ω_C, which measures the error-corrected convergence rate of protein evolution. While ω_C distinguishes natural selection from genetic noise and phylogenetic errors in simulation and real examples, its accuracy allows an exploratory genome-wide search of adaptive molecular convergence without phenotypic hypothesis or candidate genes. Using gene expression data, we explored over 20 million branch combinations in vertebrate genes and identified the joint convergence of expression patterns and protein sequences with amino acid substitutions in functionally important sites, providing hypotheses on undiscovered phenotypes. We further extended our method with a heuristic algorithm to detect highly repetitive convergence among computationally non-trivial higher-order phylogenetic combinations. Our approach allows bidirectional searches for genotype-phenotype associations, even in lineages that diverged for hundreds of millions of years.

Evaluation of Methods to Detect Shifts in Directional Selection at the Genome Scale

Identifying the footprints of selection in coding sequences can inform about the importance and function of individual sites. Analyses of the ratio of nonsynonymous to synonymous substitutions (dN/dS) have been widely used to pinpoint changes in the intensity of selection, but cannot distinguish them from changes in the direction of selection, that is, changes in the fitness of specific amino acids at a given position. A few methods that rely on amino-acid profiles to detect changes in directional selection have been designed, but their performances have not been well characterized. In this paper, we investigate the performance of six of these methods. We evaluate them on simulations along empirical phylogenies in which transition events have been annotated and compare their ability to detect sites that have undergone changes in the direction or intensity of selection to that of a widely used dN/dS approach, codeml's branch-site model A. We show that all methods have reduced performance in the presence of biased gene conversion but not CpG hypermutability. The best profile method, Pelican, a new implementation of Tamuri AU, Hay AJ, Goldstein RA. (2009. Identifying changes in selective constraints: host shifts in influenza. PLoS Comput Biol. 5(11):e1000564), performs as well as codeml in a range of conditions except for detecting relaxations of selection, and performs better when tree length increases, or in the presence of persistent positive selection. It is fast, enabling genome-scale searches for site-wise changes in the direction of selection associated with phenotypic changes.

Epistatic Effects Between Amino Acid Insertions and Substitutions Mediate Toxin resistance of Vertebrate Na+,K+-ATPases

The recurrent evolution of resistance to cardiotonic steroids (CTS) across diverse animals most frequently involves convergent amino acid substitutions in the H1-H2 extracellular loop of Na+,K+-ATPase (NKA). Previous work revealed that hystricognath rodents (e.g., chinchilla) and pterocliform birds (sandgrouse) have convergently evolved amino acid insertions in the H1-H2 loop, but their functional significance was not known. Using protein engineering, we show that these insertions have distinct effects on CTS resistance in homologs of each of the two species that strongly depend on intramolecular interactions with other residues. Removing the insertion in the chinchilla NKA unexpectedly increases CTS resistance and decreases NKA activity. In the sandgrouse NKA, the amino acid insertion and substitution Q111R both contribute to an augmented CTS resistance without compromising ATPase activity levels. Molecular docking simulations provide additional insight into the biophysical mechanisms responsible for the context-specific mutational effects on CTS insensitivity of the enzyme. Our results highlight the diversity of genetic substrates that underlie CTS insensitivity in vertebrate NKA and reveal how amino acid insertions can alter the phenotypic effects of point mutations at key sites in the same protein domain.

Genomic signatures of recent convergent transitions to social life in spiders

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

Shared Patterns of Gene Expression and Protein Evolution Associated with Adaptation to Desert Environments in Rodents

Desert specialization has arisen multiple times across rodents and is often associated with a suite of convergent phenotypes, including modification of the kidneys to mitigate water loss. However, the extent to which phenotypic convergence in desert rodents is mirrored at the molecular level is unknown. Here, we sequenced kidney mRNA and assembled transcriptomes for three pairs of rodent species to search for shared differences in gene expression and amino acid sequence associated with adaptation to deserts. We conducted phylogenetically independent comparisons between a desert specialist and a non-desert relative in three families representing ∼70 million years of evolution. Overall, patterns of gene expression faithfully recapitulated the phylogeny of these six taxa providing a strong evolutionary signal in levels of mRNA abundance. We also found that 8.6% of all genes showed shared patterns of expression divergence between desert and non-desert taxa, much of which likely reflects convergent evolution, and representing more than expected by chance under a model of independent gene evolution. In addition to these shared changes, we observed many species-pair-specific changes in gene expression indicating that instances of adaptation to deserts include a combination of unique and shared changes. Patterns of protein evolution revealed a small number of genes showing evidence of positive selection, the majority of which did not show shared changes in gene expression. Overall, our results suggest that convergent changes in gene regulation play an important role in the complex trait of desert adaptation in rodents.

Marker-free co-selection for successive rounds of prime editing in human cells

Prime editing enables the introduction of precise point mutations, small insertions, or short deletions without requiring donor DNA templates. However, efficiency remains a key challenge in a broad range of human cell types. In this work, we design a robust co-selection strategy through coediting of the ubiquitous and essential sodium/potassium pump (Na⁺/K⁺ ATPase). We readily engineer highly modified pools of cells and clones with homozygous modifications for functional studies with minimal pegRNA optimization. This process reveals that nicking the non-edited strand stimulates multiallelic editing but often generates tandem duplications and large deletions at the target site, an outcome dictated by the relative orientation of the protospacer adjacent motifs. Our approach streamlines the production of cell lines with multiple genetic modifications to create cellular models for biological research and lays the foundation for the development of cell-type specific co-selection strategies.

Prime editing enables the introduction of precise point mutations, small insertions, or short deletions without requiring donor DNA templates. Here the authors develop a co-selection strategy to facilitate prime editing in human cells and provide design principles to prevent the formation of undesired editing byproducts at the target site.

Defence mitigation by predators of chemically defended prey integrated over the predation sequence and across biological levels with a focus on cardiotonic steroids

Predator-prey interactions have long served as models for the investigation of adaptation and fitness in natural environments. Anti-predator defences such as mimicry and camouflage provide some of the best examples of evolution. Predators, in turn, have evolved sensory systems, cognitive abilities and physiological resistance to prey defences. In contrast to prey defences which have been reviewed extensively, the evolution of predator counter-strategies has received less attention. To gain a comprehensive view of how prey defences can influence the evolution of predator counter-strategies, it is essential to investigate how and when selection can operate. In this review we evaluate how predators overcome prey defences during (i) encounter, (ii) detection, (iii) identification, (iv) approach, (v) subjugation, and (vi) consumption. We focus on prey that are protected by cardiotonic steroids (CTS)-defensive compounds that are found in a wide range of taxa, and that have a specific physiological target. In this system, coevolution is well characterized between specialist insect herbivores and their host plants but evidence for coevolution between CTS-defended prey and their predators has received less attention. Using the predation sequence framework, we organize 574 studies reporting predators overcoming CTS defences, integrate these counter-strategies across biological levels of organization, and discuss the costs and benefits of attacking CTS-defended prey. We show that distinct lineages of predators have evolved dissecting behaviour, changes in perception of risk and of taste perception, and target-site insensitivity. We draw attention to biochemical, hormonal and microbiological strategies that have yet to be investigated as predator counter-adaptations to CTS defences. We show that the predation sequence framework will be useful for organizing future studies of chemically mediated systems and coevolution.

Defence mitigation by predators of chemically defended prey integrated over the predation sequence and across biological levels with a focus on cardiotonic steroids

Predator-prey interactions have long served as models for the investigation of adaptation and fitness in natural environments. Anti-predator defences such as mimicry and camouflage provide some of the best examples of evolution. Predators, in turn, have evolved sensory systems, cognitive abilities and physiological resistance to prey defences. In contrast to prey defences which have been reviewed extensively, the evolution of predator counter-strategies has received less attention. To gain a comprehensive view of how prey defences can influence the evolution of predator counter-strategies, it is essential to investigate how and when selection can operate. In this review we evaluate how predators overcome prey defences during (i) encounter, (ii) detection, (iii) identification, (iv) approach, (v) subjugation, and (vi) consumption. We focus on prey that are protected by cardiotonic steroids (CTS)-defensive compounds that are found in a wide range of taxa, and that have a specific physiological target. In this system, coevolution is well characterized between specialist insect herbivores and their host plants but evidence for coevolution between CTS-defended prey and their predators has received less attention. Using the predation sequence framework, we organize 574 studies reporting predators overcoming CTS defences, integrate these counter-strategies across biological levels of organization, and discuss the costs and benefits of attacking CTS-defended prey. We show that distinct lineages of predators have evolved dissecting behaviour, changes in perception of risk and of taste perception, and target-site insensitivity. We draw attention to biochemical, hormonal and microbiological strategies that have yet to be investigated as predator counter-adaptations to CTS defences. We show that the predation sequence framework will be useful for organizing future studies of chemically mediated systems and coevolution.

Constraints on the evolution of toxin-resistant Na,K-ATPases have limited dependence on sequence divergence

A growing body of theoretical and experimental evidence suggests that intramolecular epistasis is a major determinant of rates and patterns of protein evolution and imposes a substantial constraint on the evolution of novel protein functions. Here, we examine the role of intramolecular epistasis in the recurrent evolution of resistance to cardiotonic steroids (CTS) across tetrapods, which occurs via specific amino acid substitutions to the α-subunit family of Na,K-ATPases (ATP1A). After identifying a series of recurrent substitutions at two key sites of ATP1A that are predicted to confer CTS resistance in diverse tetrapods, we then performed protein engineering experiments to test the functional consequences of introducing these substitutions onto divergent species backgrounds. In line with previous results, we find that substitutions at these sites can have substantial background-dependent effects on CTS resistance. Globally, however, these substitutions also have pleiotropic effects that are consistent with additive rather than background-dependent effects. Moreover, the magnitude of a substitution's effect on activity does not depend on the overall extent of ATP1A sequence divergence between species. Our results suggest that epistatic constraints on the evolution of CTS-resistant forms of Na,K-ATPase likely depend on a small number of sites, with little dependence on overall levels of protein divergence. We propose that dependence on a limited number sites may account for the observation of convergent CTS resistance substitutions observed among taxa with highly divergent Na,K-ATPases (See S1 Text for Spanish translation).

Constraints on the evolution of toxin-resistant Na,K-ATPases have limited dependence on sequence divergence

A growing body of theoretical and experimental evidence suggests that intramolecular epistasis is a major determinant of rates and patterns of protein evolution and imposes a substantial constraint on the evolution of novel protein functions. Here, we examine the role of intramolecular epistasis in the recurrent evolution of resistance to cardiotonic steroids (CTS) across tetrapods, which occurs via specific amino acid substitutions to the α-subunit family of Na,K-ATPases (ATP1A). After identifying a series of recurrent substitutions at two key sites of ATP1A that are predicted to confer CTS resistance in diverse tetrapods, we then performed protein engineering experiments to test the functional consequences of introducing these substitutions onto divergent species backgrounds. In line with previous results, we find that substitutions at these sites can have substantial background-dependent effects on CTS resistance. Globally, however, these substitutions also have pleiotropic effects that are consistent with additive rather than background-dependent effects. Moreover, the magnitude of a substitution's effect on activity does not depend on the overall extent of ATP1A sequence divergence between species. Our results suggest that epistatic constraints on the evolution of CTS-resistant forms of Na,K-ATPase likely depend on a small number of sites, with little dependence on overall levels of protein divergence. We propose that dependence on a limited number sites may account for the observation of convergent CTS resistance substitutions observed among taxa with highly divergent Na,K-ATPases (See S1 Text for Spanish translation).

Convergent copy number increase of genes associated with freshwater colonization in fishes

Copy number variation (CNV) can cause phenotypic changes. However, in contrast to amino acid substitutions and cis-regulatory changes, little is known about the functional categories of genes in which CNV is important for adaptation to novel environments. It is also unclear whether the same genes repeatedly change the copy numbers for adapting to similar environments. Here, we investigate CNV associated with freshwater colonization in fishes, which was observed multiple times across different lineages. Using 48 ray-finned fishes across diverse orders, we identified 23 genes whose copy number increases were associated with freshwater colonization. These genes showed enrichment for peptide receptor activity, hexosyltransferase activity and unsaturated fatty acid metabolism. We further revealed that three of the genes showed copy number increases in freshwater populations compared to marine ancestral populations of the stickleback genus Gasterosteus. These results indicate that copy number increases of genes involved in fatty acid metabolism (FADS2), immune function (PSMB8a) and thyroid hormone metabolism (UGT2) may be important for freshwater colonization at both the inter-order macroevolutionary scale and at the intra-genus microevolutionary scale. Further analysis across diverse taxa will help to understand the role of CNV in the adaptation to novel environments. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

An enhancer of Agouti contributes to parallel evolution of cryptically colored beach mice

Identifying the genetic basis of repeatedly evolved traits provides a way to reconstruct their evolutionary history and ultimately investigate the predictability of evolution. Here, we focus on the oldfield mouse (Peromyscus polionotus), which occurs in the southeastern United States, where it exhibits considerable color variation. Dorsal coats range from dark brown in mainland mice to near white in mice inhabiting sandy beaches; this light pelage has evolved independently on Florida's Gulf and Atlantic coasts as camouflage from predators. To facilitate genomic analyses, we first generated a chromosome-level genome assembly of Peromyscus polionotus subgriseus. Next, in a uniquely variable mainland population (Peromyscus polionotus albifrons), we scored 23 pigment traits and performed targeted resequencing in 168 mice. We find that pigment variation is strongly associated with an ∼2-kb region ∼5 kb upstream of the Agouti signaling protein coding region. Using a reporter-gene assay, we demonstrate that this regulatory region contains an enhancer that drives expression in the dermis of mouse embryos during the establishment of pigment prepatterns. Moreover, extended tracts of homozygosity in this Agouti region indicate that the light allele experienced recent and strong positive selection. Notably, this same light allele appears fixed in both Gulf and Atlantic coast beach mice, despite these populations being separated by >1,000 km. Together, our results suggest that this identified Agouti enhancer allele has been maintained in mainland populations as standing genetic variation and from there, has spread to and been selected in two independent beach mouse lineages, thereby facilitating their rapid and parallel evolution.

Ehrlich and Raven escape and radiate coevolution hypothesis at different levels of organization: Past and future perspectives

The classic paper by Ehrlich and Raven on coevolution will soon be 60 years old. Although they were not the first to develop the idea of coevolution, their thought-provoking paper certainly popularized this idea and inspired several generations of scientists interested in coevolution. Here, we describe some of their main contributions, quantitatively measure the impact of their seminal paper on different fields of research, and discuss how ideas related to their original paper might push the study of coevolution forward. To guide our discussion, we explore their original hypothesis into three research fields that are associated with distinct scales/levels of organization: (1) the genetic mechanisms underlying coevolutionary interactions; (2) the potential association between coevolutionary diversification and the organization of ecological networks; and (3) the micro- and macroevolutionary mechanisms and expected patterns under their hypothesis. By doing so, we discuss potentially overlooked aspects and future directions for the study of coevolutionary dynamics and diversification.

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.

Convergent evolution of toxin resistance in animals

Convergence is the phenomenon whereby similar phenotypes evolve independently in different lineages. One example is resistance to toxins in animals. Toxins have evolved many times throughout the tree of life. They disrupt molecular and physiological pathways in target species, thereby incapacitating prey or deterring a predator. In response, molecular resistance has evolved in many species exposed to toxins to counteract their harmful effects. Here, we review current knowledge on the convergence of toxin resistance using examples from a wide range of toxin families. We explore the evolutionary processes and molecular adaptations driving toxin resistance. However, resistance adaptations may carry a fitness cost if they disrupt the normal physiology of the resistant animal. Therefore, there is a trade‐off between maintaining a functional molecular target and reducing toxin susceptibility. There are relatively few solutions that satisfy this trade‐off. As a result, we see a small set of molecular adaptations appearing repeatedly in diverse animal lineages, a phenomenon that is consistent with models of deterministic evolution. Convergence may also explain what has been called ‘autoresistance’. This is often thought to have evolved for self‐protection, but we argue instead that it may be a consequence of poisonous animals feeding on toxic prey. Toxin resistance provides a unique and compelling model system for studying the interplay between trophic interactions, selection pressures and the molecular mechanisms underlying evolutionary novelties.

Two Apriona Species Sharing a Host Niche Have Different Gut Microbiome Diversity

The adaptability of herbivorous insects to toxic plant defense compounds is partly related to the structure of the gut microbiome. To overcome plant resistance, the insect gut microbiome should respond to a wide range of allelochemicals derived from dietary niches. Nevertheless, for sibling herbivorous insect species, whether the gut microbiome contributes to success in food niche competition is unclear. Based on 16S rDNA high-throughput sequencing, the gut microbiomes of two Apriona species that share the same food niche were investigated in this study to determine whether the gut microbiome contributes to insect success in food-niche competition. Our observations indicated that the gut microbiome tended to play a part in host niche competition between the two Apriona species. The gut microbiome of Apriona swainsoni had many enriched pathways that can help degrade plant toxic secondary compounds, including xenobiotic biodegradation and metabolism, terpenoid and polyketide metabolism, and secondary metabolite biosynthesis. Meanwhile, A. swainsoni hosted a much greater variety of microorganisms and had more viable bacteria than A. germari. We conclude that gut microbes may influence the coevolution of herbivores and host plants. Gut bacteria may not only serve to boost nutritional relationships, but may also play an important role in insect food niche competition.

Evidence for tissue-specific defense-offense interactions between milkweed and its community of specialized herbivores

Coevolution between plants and herbivores often involves escalation of defense-offense strategies, but attack by multiple herbivores may obscure the match of plant defense to any one attacker. As herbivores often specialize on distinct plant parts, we hypothesized that defense-offense interactions in coevolved systems may become physiologically and evolutionarily compartmentalized between plant tissues. We report that roots, leaves, flower buds and seeds of the tropical milkweed (Asclepias curassavica) show increasing concentrations of cardenolide toxins acropetally, with latex showing the highest concentration. In vitro assays of the physiological target of cardenolides, the Na⁺ /K⁺ -ATPase (hereafter 'sodium pump'), of three specialized milkweed herbivores (root-feeding Tetraopes tetrophthalmus, leaf-feeding Danaus plexippus, and seed-feeding Oncopeltus fasciatus) show that they are proportionally tolerant to the cardenolide concentrations of the tissues they eat. Indeed, molecular substitutions in the insects' sodium pumps predicted their tolerance to toxins from their target tissues. Nonetheless, the relative inhibition of the sodium pumps of these specialists by the concentration vs. composition (inhibition controlled for concentration, what we term 'potency') of cardenolides from their target vs. non-target plant tissues revealed different degrees of insect adaptation to tissue-specific toxins. In addition, a trade-off between toxin concentration and potency emerged across plant tissues, potentially reflecting coevolutionary history or plant physiological constraints. Our findings suggest that tissue-specific coevolutionary dynamics may be proceeding between the plant and its specialized community of herbivores. This novel finding may be common in nature, contributing to ways in which coevolution proceeds in multi-species communities.

Na,K-ATPase α1 and β-subunits show distinct localizations in the nervous tissue of the large milkweed bug

The Na,K-ATPase (NKA) is an essential ion transporter and signaling molecule in all animal tissues and believed to consist at least one α and one ß-subunit to form a functional enzyme. In the large milkweed bug, Oncopeltus fasciatus, adaptation to dietary cardiac glycosides (CGs), which can fatally block the NKA, has resulted in gene duplications leading to four α1-subunits. These differ in sensitivity to CGs, but resistance trades off against ion pumping activity, thus influencing the α1-subunits’ suitability for specific tissues. Besides, O. fasciatus possesses four different ß-subunits that can alter the NKA's kinetics and should play an essential role in the formation of cellular junctions.

Proteomic analyses revealed the distribution and composition of α1/ß-complexes in the nervous tissue of O. fasciatus. The highly CG-resistant, but less active α1B and the highly active, but less resistant α1C predominated in the nervous tissue and co-occurred with ß2 and ß3, partly forming larger complexes than just heterodimers. Immunohistochemical analyses provided a fine scale resolution of the subunits’ distribution in different morphological structures of the nervous tissue. This may suggest that α1 as well as ß-subunits occur in isolation without the other subunit, which contradicts the present understanding that the two types of subunits have to associate to form functional complexes. An isolated occurrence was especially prominent for ß3 and βx, the enigmatic fourth and N-terminally largely truncated ß-subunit. We hypothesize that dimerization of these ß-subunits plays a role in cell–cell contacts.

Evolution in small steps and giant leaps

The first Editor of Evolution was Ernst Mayr. His foreword to the first issue of Evolution published in 1947 framed evolution as a "problem of interaction" that was just beginning to be studied in this broad context. First, I explore progress and prospects on understanding the subsidiary interactions identified by Mayr, including interactions between parts of organisms, between individuals and populations, between species, and between the organism and its abiotic environment. Mayr's overall "problem of interaction" framework is examined in the context of coevolution within and among levels of biological organization. This leads to a comparison in the relative roles of biotic versus abiotic agents of selection and fluctuating versus directional selection, followed by stabilizing selection in shaping the genomic architecture of adaptation. Oligogenic architectures may be typical for traits shaped more by fluctuating selection and biotic selection. Conversely, polygenic architectures may be typical for traits shaped more by directional followed by stabilizing selection and abiotic selection. The distribution of effect sizes and turnover dynamics of adaptive alleles in these scenarios deserves further study. Second, I review two case studies on the evolution of acquired toxicity in animals, one involving cardiac glycosides obtained from plants and one involving bacterial virulence factors horizontally transferred to animals. The approaches used in these studies and the results gained directly flow from Mayr's vision of an evolutionary biology that revolves around the "problem of interaction."

Sequestration of Plant Defense Compounds by Insects: From Mechanisms to Insect-Plant Coevolution

Plant defense compounds play a key role in the evolution of insect-plant associations by selecting for behavioral, morphological, and physiological insect adaptations. Sequestration, the ability of herbivorous insects to accumulate plant defense compounds to gain a fitness advantage, represents a complex syndrome of adaptations that has evolved in all major lineages of herbivorous insects and involves various classes of plant defense compounds. In this article, we review progress in understanding how insects selectively accumulate plant defense metabolites and how the evolution of specific resistance mechanisms to these defense compounds enables sequestration. These mechanistic considerations are further integrated into the concept of insect-plant coevolution. Comparative genome and transcriptome analyses, combined with approaches based on analytical chemistry that are centered in phylogenetic frameworks, will help to reveal adaptations underlying the sequestration syndrome, which is essential to understanding the influence of sequestration on insect-plant coevolution.

Parallel genomic responses to historical climate change and high elevation in East Asian songbirds

Parallel evolution can be expected among closely related taxa exposed to similar selective pressures. However, parallelism is typically stronger at the phenotypic level, while genetic solutions to achieve these phenotypic similarities may differ. For polygenic traits, the availability of standing genetic variation (i.e., heterozygosity) may influence such genetic nonparallelism. Here, we examine the extent to which high-elevation adaptation is parallel-and whether the level of parallelism is affected by heterozygosity-by analyzing genomes of 19 Paridae species distributed across East Asia with a dramatic east-west elevation gradient. We find that western highlands endemic parids have consistently lower levels of heterozygosity-likely the result of late-Pleistocene demographic contraction-than do parids found exclusively in eastern lowlands, which remained unglaciated during the late Pleistocene. Three widespread species (east to west) have high levels of heterozygosity similar to that observed in eastern species, although their western populations are less variable than eastern ones. Comparing genomic responses to extreme environments of the Qinghai-Tibet Plateau, we find that the most differentiated genomic regions between each high-elevation taxon and its low-elevation relative are significantly enriched for genes potentially related to the oxygen transport cascade and/or thermogenesis. Despite no parallelism at particular genes, high similarity in gene function is found among comparisons. Furthermore, parallelism is not higher in more heterozygous widespread parids than in highland endemics. Thus, in East Asian parids, parallel functional response to extreme elevation appears to rely on different genes, with differences in heterozygosity having no effect on the degree of genetic parallelism.

Dietary cardenolides enhance growth and change the direction of the fecundity-longevity trade-off in milkweed bugs (Heteroptera: Lygaeinae)

Sequestration, that is, the accumulation of plant toxins into body tissues for defense, was predicted to incur physiological costs and may require resistance traits different from those of non-sequestering insects. Alternatively, sequestering species could experience a cost in the absence of toxins due to selection on physiological homeostasis under permanent exposure of sequestered toxins in body tissues. Milkweed bugs (Heteroptera: Lygaeinae) sequester high amounts of plant-derived cardenolides. Although being potent inhibitors of the ubiquitous animal enzyme Na+/K+-ATPase, milkweed bugs can tolerate cardenolides by means of resistant Na+/K+-ATPases. Both adaptations, resistance and sequestration, are ancestral traits of the Lygaeinae. Using four milkweed bug species (Heteroptera: Lygaeidae: Lygaeinae) and the related European firebug (Heteroptera: Pyrrhocoridae: Pyrrhocoris apterus) showing different combinations of the traits "cardenolide resistance" and "cardenolide sequestration," we tested how the two traits affect larval growth upon exposure to dietary cardenolides in an artificial diet system. While cardenolides impaired the growth of P. apterus nymphs neither possessing a resistant Na+/K+-ATPase nor sequestering cardenolides, growth was not affected in the non-sequestering milkweed bug Arocatus longiceps, which possesses a resistant Na+/K+-ATPase. Remarkably, cardenolides increased growth in the sequestering dietary specialists Caenocoris nerii and Oncopeltus fasciatus but not in the sequestering dietary generalist Spilostethus pandurus, which all possess a resistant Na+/K+-ATPase. We furthermore assessed the effect of dietary cardenolides on additional life history parameters, including developmental speed, longevity of adults, and reproductive success in O. fasciatus. Unexpectedly, nymphs under cardenolide exposure developed substantially faster and lived longer as adults. However, fecundity of adults was reduced when maintained on cardenolide-containing diet for their entire lifetime but not when adults were transferred to non-toxic sunflower seeds. We speculate that the resistant Na+/K+-ATPase of milkweed bugs is selected for working optimally in a "toxic environment," that is, when sequestered cardenolides are stored in the body.