Historical Contingency in a Multigene Family Facilitates Adaptive Evolution of Toxin Resistance

Where Does All the Poison Go? Investigating Toxicokinetics of Newt (Taricha) Tetrodotoxin (TTX) in Garter Snakes (Thamnophis)

Animals that consume toxic diets provide models for understanding the molecular and physiological adaptations to ecological challenges. Garter snakes (Thamnophis) in western North America prey on Pacific newts (Taricha), which employ tetrodotoxin (TTX) as an antipredator defense. These snakes possess mutations in voltage-gated sodium channels (Nav), the molecular targets of TTX, that decrease the binding ability of TTX to sodium channels (target-site resistance). However, genetic variation at these loci that cannot explain all the phenotypic variation in TTX resistance in Thamnophis. We explored a separate means of resistance, toxin metabolism, to determine if TTX-resistant snakes either rapidly remove TTX or sequester TTX. We examined the metabolism and distribution of TTX in the body (toxicokinetics), to determine differences between TTX-resistant and TTX-sensitive snakes in the rates at which TTX is eliminated from organs and the whole body (using TTX half-life as our metric). We assayed TTX half-life in snakes from TTX-resistant and TTX-sensitive populations of three garter snake species with a coevolutionary history with newts (T. atratus, T. couchii, T. sirtalis), as well as two non-resistant "outgroup" species (T. elegans, Pituophis catenifer) that seldom (if ever) engage newts. We found TTX half-life varied across species, populations, and tissues. Interestingly, TTX half-life was shortest in T. elegans and P. catenifer compared to all other snakes. Furthermore, TTX-resistant populations of T. couchii and T. sirtalis eliminated TTX faster (shorter TTX half-life) than their TTX-sensitive counterparts, while populations of TTX-resistant and TTX-sensitive T. atratus showed no difference rates of TTX removal (same TTX half-life). The ability to rapidly eliminate TTX may have permitted increased prey consumption, which may have promoted the evolution of additional resistance mechanisms. Finally, snakes still retain substantial amounts of TTX, and we projected that snakes could be dangerous to their own predators days to weeks following the ingestion of a single newt. Thus, aspects of toxin metabolism may have been key in driving predator-prey relationships, and important in determining other ecological interactions.

Coevolution with toxic prey produces functional trade-offs in sodium channels of predatory snakes

Seemingly unrelated traits often share the same underlying molecular mechanisms, potentially generating a pleiotropic relationship whereby selection shaping one trait can simultaneously compromise another. While such functional trade-offs are expected to influence evolutionary outcomes, their actual relevance in nature is masked by obscure links between genotype, phenotype, and fitness. Here, we describe functional trade-offs that likely govern a key adaptation and coevolutionary dynamics in a predator-prey system. Several garter snake (Thamnophis spp.) populations have evolved resistance to tetrodotoxin (TTX), a potent chemical defense in their prey, toxic newts (Taricha spp.). Snakes achieve TTX resistance through mutations occurring at toxin-binding sites in the pore of snake skeletal muscle voltage-gated sodium channels (NaV1.4). We hypothesized that these mutations impair basic NaV functions, producing molecular trade-offs that should ultimately scale up to compromised organismal performance. We investigate biophysical costs in two snake species with unique and independently evolved mutations that confer TTX resistance. We show electrophysiological evidence that skeletal muscle sodium channels encoded by toxin-resistant alleles are functionally compromised. Furthermore, skeletal muscles from snakes with resistance genotypes exhibit reduced mechanical performance. Lastly, modeling the molecular stability of these sodium channel variants partially explains the electrophysiological and muscle impairments. Ultimately, adaptive genetic changes favoring toxin resistance appear to negatively impact sodium channel function, skeletal muscle strength, and organismal performance. These functional trade-offs at the cellular and organ levels appear to underpin locomotor deficits observed in resistant snakes and may explain variation in the population-level success of toxin-resistant alleles across the landscape, ultimately shaping the trajectory of snake-newt coevolution.

Molecular Mechanisms Underlying Vertebrate Adaptive Evolution: A Systematic Review

Adaptive evolution is a process in which variation that confers an evolutionary advantage in a specific environmental context arises and is propagated through a population. When investigating this process, researchers have mainly focused on describing advantageous phenotypes or putative advantageous genotypes. A recent increase in molecular data accessibility and technological advances has allowed researchers to go beyond description and to make inferences about the mechanisms underlying adaptive evolution. In this systematic review, we discuss articles from 2016 to 2022 that investigated or reviewed the molecular mechanisms underlying adaptive evolution in vertebrates in response to environmental variation. Regulatory elements within the genome and regulatory proteins involved in either gene expression or cellular pathways have been shown to play key roles in adaptive evolution in response to most of the discussed environmental factors. Gene losses were suggested to be associated with an adaptive response in some contexts. Future adaptive evolution research could benefit from more investigations focused on noncoding regions of the genome, gene regulation mechanisms, and gene losses potentially yielding advantageous phenotypes. Investigating how novel advantageous genotypes are conserved could also contribute to our knowledge of adaptive evolution.

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.

In Silico Analysis of Tetrodotoxin Binding in Voltage-Gated Sodium Ion Channels from Toxin-Resistant Animal Lineages

Multiple animal species have evolved resistance to the neurotoxin tetrodotoxin (TTX) through changes in voltage-gated sodium ion channels (VGSCs). Amino acid substitutions in TTX-resistant lineages appear to be positionally convergent with changes in homologous residues associated with reductions in TTX block. We used homology modeling coupled with docking simulations to test whether positionally convergent substitutions generate functional convergence at the level of TTX–channel interactions. We found little evidence that amino acids at convergent positions generated similar patterns among TTX-resistant animal lineages across several metrics, including number of polar contacts, polar contact position, and estimates of binding energy. Though binding energy values calculated for TTX docking were reduced for some TTX-resistant channels, not all TTX-resistant channels and not all of our analyses returned reduced binding energy values for TTX-resistant channels. Our results do not support a simple model of toxin resistance where a reduced number of bonds between TTX and the channel protein prevents blocking. Rather models that incorporate flexibility and movement of the protein overall may better describe how homologous substitutions in the channel cause changes in TTX block.

Historical contingency drives compensatory evolution and rare reversal of phage resistance

Bacteria and lytic viruses (phages) engage in highly dynamic coevolutionary interactions over time, yet we have little idea of how transient selection by phages might shape the future evolutionary trajectories of their host populations. To explore this question, we generated genetically diverse phage-resistant mutants of the bacterium Pseudomonas syringae. We subjected the panel of mutants to prolonged experimental evolution in the absence of phages. Some populations re-evolved phage sensitivity, whereas others acquired compensatory mutations that reduced the costs of resistance without altering resistance levels. To ask whether these outcomes were driven by the initial genetic mechanisms of resistance, we next evolved independent replicates of each individual mutant in the absence of phages. We found a strong signature of historical contingency: some mutations were highly reversible across replicate populations, whereas others were highly entrenched. Through whole-genome sequencing of bacteria over time, we also found that populations with the same resistance gene acquired more parallel sets of mutations than populations with different resistance genes, suggesting that compensatory adaptation is also contingent on how resistance initially evolved. Our study identifies an evolutionary ratchet in bacteria–phage coevolution and may explain previous observations that resistance persists over time in some bacterial populations but is lost in others. We add to a growing body of work describing the key role of phages in the ecological and evolutionary dynamics of their host communities. Beyond this specific trait, our study provides a new insight into the genetic architecture of historical contingency, a crucial component of interpreting and predicting evolution.

Venom resistance mechanisms in centipede show tissue specificity

Venomous animals utilize venom glands to secrete and store powerful toxins for intraspecific and/or interspecific antagonistic interactions, implying that tissue-specific resistance is essential for venom glands to anatomically separate toxins from other tissues. Here, we show the mechanism of tissue-specific resistance in centipedes (Scolopendra subspinipes mutilans), where the splice variant of the receptor repels its own toxin. Unlike the well-known resistance mechanism by mutation in a given exon, we found that the KCNQ1 channel is highly expressed in the venom gland as a unique splice variant in which the pore domain and transmembrane domain six, partially encoded by exon 6 (rather than 7 as found in other tissues), contain eleven mutated residues. Such a splice variant is sufficient to gain resistance to SsTx (a lethal toxin for giant prey capture) in the venom gland due to a partially buried binding site. Therefore, the tissue-specific KCNQ1 modification confers resistance to the toxins, establishing a safe zone in the venom-storing/secreting environment.

The road not taken: Evolution of tetrodotoxin resistance in the Sierra garter snake (Thamnophis couchii) by a path less traveled

The repeated evolution of tetrodotoxin (TTX) resistance provides a model for testing hypotheses about the mechanisms of convergent evolution. This poison is broadly employed as a potent antipredator defense, blocking voltage-gated sodium channels (Na_v ) in muscles and nerves, paralyzing and sometimes killing predators. Resistance in taxa bearing this neurotoxin and a few predators appears to come from convergent replacements in specific Na_v residues that interact with TTX. This stereotyped genetic response suggests molecular and phenotypic evolution may be constrained and predictable. Here, we investigate the extent of mechanistic convergence in garter snakes (Thamnophis) that prey on TTX-bearing newts (Taricha) by examining the physiological and genetic basis of TTX resistance in the Sierra garter snake (Th. couchii). We characterize variation in this predatory adaptation across populations at several biological scales: whole-animal TTX resistance; skeletal muscle resistance, functional genetic variation in three Na_v encoding loci; and levels of gene expression for one of these loci. We found Th. couchii possess extensive geographic variation in resistance at the whole-animal and skeletal muscle levels. As in other Thamnophis, resistance at both levels is highly correlated, suggesting convergence across the biological levels linking organism to organ. However, Th. couchii shows no functional variation in Na_v loci among populations or difference in candidate gene expression. Local variation in TTX resistance in Th. couchii cannot be explained by the same relationship between genotype and phenotype seen in other taxa. Thus, historical contingencies may lead different species of Thamnophis down alternative routes to local adaptation.

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.

A review of chemical defense in harlequin toads (Bufonidae: Atelopus)

Toads of the genus Atelopus are chemically defended by a unique combination of endogenously synthesized cardiotoxins (bufadienolides) and neurotoxins which may be sequestered (guanidinium alkaloids). Investigation into Atelopus small-molecule chemical defenses has been primarily concerned with identifying and characterizing various forms of these toxins while largely overlooking their ecological roles and evolutionary implications. In addition to describing the extent of knowledge about Atelopus toxin structures, pharmacology, and biological sources, we review the detection, identification, and quantification methods used in studies of Atelopus toxins to date and conclude that many known toxin profiles are unlikely to be comprehensive because of methodological and sampling limitations. Patterns in existing data suggest that both environmental (toxin availability) and genetic (capacity to synthesize or sequester toxins) factors influence toxin profiles. From an ecological and evolutionary perspective, we summarize the possible selective pressures acting on Atelopus toxicity and toxin profiles, including predation, intraspecies communication, disease, and reproductive status. Ultimately, we intend to provide a basis for future ecological, evolutionary, and biochemical research on Atelopus.

How to approach the study of syndromes in macroevolution and ecology

Syndromes, wherein multiple traits evolve convergently in response to a shared selective driver, form a central concept in ecology and evolution. Recent work has questioned the existence of some classic syndromes, such as pollination and seed dispersal syndromes. Here, we discuss some of the major issues that have afflicted research into syndromes in macroevolution and ecology. First, correlated evolution of traits and hypothesized selective drivers is often relied on as the only evidence for adaptation of those traits to those hypothesized drivers, without supporting evidence. Second, the selective driver is often inferred from a combination of traits without explicit testing. Third, researchers often measure traits that are easy for humans to observe rather than measuring traits that are suited to testing the hypothesis of adaptation. Finally, species are often chosen for study because of their striking phenotypes, which leads to the illusion of syndromes and divergence. We argue that these issues can be avoided by combining studies of trait variation across entire clades or communities with explicit tests of adaptive hypotheses and that taking this approach will lead to a better understanding of syndrome-like evolution and its drivers.

Phylogenomic analyses resolve relationships among garter snakes (Thamnophis: Natricinae: Colubridae) and elucidate biogeographic history and morphological evolution

Garter snakes (Thamnophis) are a successful group of natricines endemic to North America. They have become important natural models for ecological and evolutionary research, yet prior efforts to resolve phylogenetic relationships have resulted in conflicting topologies and weak support for certain relationships. Here, we use genomic data generated with a reduced representation double-digest RADseq approach to reassess evolutionary relationships across Thamnophis. We then use the resulting phylogeny to better understand how biogeography and feeding ecology have influenced lineage diversification and morphological evolution. We recovered highly congruent and strongly supported topologies from maximum likelihood and Bayesian analyses, but some discordance with a multispecies coalescent approach. All phylogenomic estimates split Thamnophis into two clades largely defined by northern and southern North American species. Divergence time estimates and biogeographic analyses indicate a mid-Miocene origin of Thamnophis in Mexico. In addition, historic vicariant events thought to explain biogeographic patterns in other lineages (e.g., Isthmus of Tehuantepec, Rocky Mountain Range, and Trans-Mexican Volcanic Belt) appear to have influenced patterns of diversification in Thamnophis as well. Analyses of morphological traits associated with feeding ecology showed moderate to strong phylogenetic signal. Nevertheless, phylogenetic ANOVA suggested significant differences in certain cranial morphologies between aquatic specialists and garter snakes that are terrestrial-aquatic generalists, independent of evolutionary history. Our new estimate of Thamnophis phylogeny yields an improved understanding of the biogeographic history and morphological evolution of garter snakes, and provides a robust framework for future research on these snakes.

Evidence that toxin resistance in poison birds and frogs is not rooted in sodium channel mutations and may rely on "toxin sponge" proteins

Many poisonous organisms carry small-molecule toxins that alter voltage-gated sodium channel (NaV) function. Among these, batrachotoxin (BTX) from Pitohui poison birds and Phyllobates poison frogs stands out because of its lethality and unusual effects on NaV function. How these toxin-bearing organisms avoid autointoxication remains poorly understood. In poison frogs, a NaV DIVS6 pore-forming helix N-to-T mutation has been proposed as the BTX resistance mechanism. Here, we show that this variant is absent from Pitohui and poison frog NaVs, incurs a strong cost compromising channel function, and fails to produce BTX-resistant channels in poison frog NaVs. We also show that captivity-raised poison frogs are resistant to two NaV-directed toxins, BTX and saxitoxin (STX), even though they bear NaVs sensitive to both. Moreover, we demonstrate that the amphibian STX "toxin sponge" protein saxiphilin is able to protect and rescue NaVs from block by STX. Taken together, our data contradict the hypothesis that BTX autoresistance is rooted in the DIVS6 N→T mutation, challenge the idea that ion channel mutations are a primary driver of toxin resistance, and suggest the possibility that toxin sequestration mechanisms may be key for protecting poisonous species from the action of small-molecule toxins.

From Poison to Promise: The Evolution of Tetrodotoxin and Its Potential as a Therapeutic

Tetrodotoxin (TTX) is a potent neurotoxin that was first identified in pufferfish but has since been isolated from an array of taxa that host TTX-producing bacteria. However, determining its origin, ecosystem roles, and biomedical applications has challenged researchers for decades. Recognized as a poison and for its lethal effects on humans when ingested, TTX is primarily a powerful sodium channel inhibitor that targets voltage-gated sodium channels, including six of the nine mammalian isoforms. Although lethal doses for humans range from 1.5-2.0 mg TTX (blood level 9 ng/mL), when it is administered at levels far below LD₅₀, TTX exhibits therapeutic properties, especially to treat cancer-related pain, neuropathic pain, and visceral pain. Furthermore, TTX can potentially treat a variety of medical ailments, including heroin and cocaine withdrawal symptoms, spinal cord injuries, brain trauma, and some kinds of tumors. Here, we (i) describe the perplexing evolution and ecology of tetrodotoxin, (ii) review its mechanisms and modes of action, and (iii) offer an overview of the numerous ways it may be applied as a therapeutic. There is much to be explored in these three areas, and we offer ideas for future research that combine evolutionary biology with therapeutics. The TTX system holds great promise as a therapeutic and understanding the origin and chemical ecology of TTX as a poison will only improve its general benefit to humanity.

Every breath you take: assessing metabolic costs of toxin resistance in garter snakes (Thamnophis)

Trait specialization often comes at the expense of original trait function, potentially causing evolutionary tradeoffs that may render specialist populations vulnerable to extinction. However, many specialized adaptations evolve repeatedly, suggesting selection favors specialization in specific environments. Some garter snake (Thamnophis) populations possess specialized mutations in voltage-gated sodium channels that allow them to consume Pacific newts (Taricha) defended by a highly potent neurotoxin (tetrodotoxin). These mutations, however, also decrease protein and muscle function, suggesting garter snakes may suffer evolutionary tradeoffs. We measured a key physiological process, standard metabolic rate (SMR), to investigate whether specialized adaptations in toxin-resistant garter snakes affect baseline energy expenditure. In snakes, skeletal muscles influence metabolism and power ventilation, so inefficiencies of sodium channels in these muscles might impact whole-animal energy expenditure. Further, because sodium channels are membrane-bound proteins, inefficiencies of channel kinetics and performance might be exacerbated at suboptimal temperatures. We measured SMR in 2 species, Thamnophis atratus and Thamnophis sirtalis, that independently evolved tetrodotoxin resistance through unique mutations, providing replicate experiments with distinct underlying genetics and potential physiological costs. Despite our expectations, neither resistance phenotype nor sodium channel genotype affected metabolism and resistant snakes did not perform worse under suboptimal body temperature. Instead, T. atratus and T. sirtalis show nearly identical rates of mass-adjusted energy expenditure at both temperatures, despite differing eco-morphologies, life histories, and distant phylogenetic positions. These findings suggest SMR may be a conserved feature of Thamnophis, and that any organismal tradeoffs may be compensated to retain whole-animal function.

Gene Conversion Facilitates the Adaptive Evolution of Self-Resistance in Highly Toxic Newts

Reconstructing the histories of complex adaptations and identifying the evolutionary mechanisms underlying their origins are two of the primary goals of evolutionary biology. Taricha newts, which contain high concentrations of the deadly toxin tetrodotoxin (TTX) as an antipredator defense, have evolved resistance to self-intoxication, which is a complex adaptation requiring changes in six paralogs of the voltage-gated sodium channel (Na_v) gene family, the physiological target of TTX. Here, we reconstruct the origins of TTX self-resistance by sequencing the entire Na_v gene family in newts and related salamanders. We show that moderate TTX resistance evolved early in the salamander lineage in three of the six Na_v paralogs, preceding the proposed appearance of tetrodotoxic newts by ∼100 My. TTX-bearing newts possess additional unique substitutions across the entire Na_v gene family that provide physiological TTX resistance. These substitutions coincide with signatures of positive selection and relaxed purifying selection, as well as gene conversion events, that together likely facilitated their evolution. We also identify a novel exon duplication within Na_v1.4 encoding an expressed TTX-binding site. Two resistance-conferring changes within newts appear to have spread via nonallelic gene conversion: in one case, one codon was copied between paralogs, and in the second, multiple substitutions were homogenized between the duplicate exons of Na_v1.4. Our results demonstrate that gene conversion can accelerate the coordinated evolution of gene families in response to a common selection pressure.

The chemistry and biology of guanidine secondary metabolites

Covering: 2017-2019Guanidine natural products isolated from microorganisms, marine invertebrates and terrestrial plants, amphibians and spiders, represented by non-ribosomal peptides, guanidine-bearing polyketides, alkaloids, terpenoids and shikimic acid derived, are the subject of this review. The topics include the discovery of new metabolites, total synthesis of natural guanidine compounds, biological activity and mechanism-of-action, biosynthesis and ecological functions.

The geographic mosaic in parallel: Matching patterns of newt tetrodotoxin levels and snake resistance in multiple predator-prey pairs

The Geographic Mosaic Theory of Coevolution predicts that coevolutionary arms races will vary over time and space because of the diverse ecological settings and population histories of interacting species across the landscape. Thus, understanding coevolution may require investigating broad sets of populations sampled across the range of the interaction. In addition, comparing coevolutionary dynamics between similar systems may reveal the importance of specific factors that structure coevolution. Here, we examine geographic patterns of prey traits and predator traits in the relatively unstudied interaction between the Sierra garter snake (Thamnophis couchii) and sympatric prey, the rough-skinned newt (Taricha granulosa), Sierra newt (Ta. sierrae) and California newt (Ta. torosa). This system parallels, in space and phenotypes, a classic example of coevolution between predatory common garter snakes (Th. sirtalis) and their toxic newt prey exhibiting hotspots of newt tetrodotoxin (TTX) levels and matching snake TTX resistance. We quantified prey and predator traits from hundreds of individuals across their distributions, and functional trait matching at sympatric sites. We show strong regional patterns of trait covariation across the shared ranges of Th. couchii and newt prey. Traits differ significantly among localities, with lower newt TTX levels and snake TTX resistance at the northern latitudes, and higher TTX levels and snake resistance at southern latitudes. Newts and snakes in northern populations show the highest degree of functional trait matching despite possessing the least extreme traits. Conversely, newts and snakes in southern populations show the greatest mismatch despite possessing exaggerated traits, with some snakes so resistant to TTX they would be unaffected by any sympatric newt. Nevertheless, individual variation was substantial, and appears to offer the opportunity for continued reciprocal selection in most populations. Overall, the three species of newts appear to be engaged in a TTX-mediated arms race with Th. couchii. These patterns are congruent with those seen between newts and Th. sirtalis, including the same latitudinal gradient in trait covariation, and the potential 'escape' from the arms race by snake predators. Such concordance in broad scale patterns across two distinct systems suggests common phenomena might structure geographic mosaics in similar ways.

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.

Sex linkage of the skeletal muscle sodium channel gene (SCN4A) explains apparent deviations from Hardy-Weinberg equilibrium of tetrodotoxin-resistance alleles in garter snakes (Thamnophis sirtalis)

The arms race between tetrodotoxin-bearing Pacific newts (Taricha) and their garter snake predators (Thamnophis) in western North America has become a classic example of coevolution, shedding light on predator-prey dynamics, the molecular basis of adaptation, and patterns of convergent evolution. Newts are defended by tetrodotoxin (TTX), a neurotoxin that binds to voltage-gated sodium channels (Na_v proteins), arresting electrical activity in nerves and muscles and paralyzing would-be predators. However, populations of the common garter snake (T. sirtalis) have overcome this defense, largely through polymorphism at the locus SCN4A, which renders the encoded protein (Na_v1.4) less vulnerable to TTX. Previous work suggests that SCN4A commonly shows extreme deviations from Hardy-Weinberg equilibrium (HWE) in these populations, which has been interpreted as the result of intense selection imposed by newts. Here we show that much of this apparent deviation can be attributed to sex linkage of SCN4A. Using genomic data and quantitative PCR, we show that SCN4A is on the Z chromosome in Thamnophis and other advanced snakes. Taking Z-linkage into account, we find that most apparent deviations from HWE can be explained by female hemizygosity rather than low heterozygosity. Sex linkage can affect mutation rates, selection, and drift, and our results suggest that Z-linkage of SCN4A may make significant contributions to the overall dynamics of the coevolutionary arms race between newts and snakes.

Tolerance of Novel Toxins through Generalized Mechanisms: Simulating Gradual Host Shifts of Butterflies

Organisms encounter a wide range of toxic compounds in their environments, from chemicals that serve anticonsumption or anticompetition functions to pollutants and pesticides. Although we understand many detoxification mechanisms that allow organisms to consume toxins typical of their diet, we know little about why organisms vary in their ability to tolerate entirely novel toxins. We tested whether variation in generalized stress responses, such as antioxidant pathways, may underlie variation in reactions to novel toxins and, if so, their associated costs. We used an artificial diet to present cabbage white butterfly caterpillars (Pieris rapae) with plant material containing toxins not experienced in their evolutionary history. Families that maintained high performance (e.g., high survival, fast development time, large body size) on diets containing one novel toxic plant also performed well when exposed to two other novel toxic plants, consistent with a generalized response. Variation in constitutive (but not induced) expression of genes involved in oxidative stress responses was positively related to performance on the novel diets. While we did not detect reproductive trade-offs of this generalized response, there was a tendency to have less melanin investment in the wings, consistent with the role of melanin in oxidative stress responses. Taken together, our results support the hypothesis that variation in generalized stress responses, such as genes involved in oxidative stress responses, may explain the variation in tolerance to entirely novel toxins and may facilitate colonization of novel hosts and environments.

Duplication history and molecular evolution of the rbcS multigene family in angiosperms

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) is considered to be the main enzyme determining the rate of photosynthesis. The small subunit of the protein, encoded by the rbcS gene, has been shown to influence the catalytic efficiency, CO2 specificity, assembly, activity, and stability of RuBisCO. However, the evolution of the rbcS gene remains poorly studied. We inferred the phylogenetic tree of the rbcS gene in angiosperms using the nucleotide sequences and found that it is composed of two lineages that may have existed before the divergence of land plants. Although almost all species sampled carry at least one copy of lineage 1, genes of lineage 2 were lost in most angiosperm species. We found the specific residues that have undergone positive selection during the evolution of the rbcS gene. We detected intensive coevolution between each rbcS gene copy and the rbcL gene encoding the large subunit of RuBisCO. We tested the role played by each rbcS gene copy on the stability of the RuBisCO protein through homology modelling. Our results showed that this evolutionary constraint could limit the level of divergence seen in the rbcS gene, which leads to the similarity among the rbcS gene copies of lineage 1 within species.

Adaptive substitutions underlying cardiac glycoside insensitivity in insects exhibit epistasis in vivo

Predicting how species will respond to selection pressures requires understanding the factors that constrain their evolution. We use genome engineering of Drosophila to investigate constraints on the repeated evolution of unrelated herbivorous insects to toxic cardiac glycosides, which primarily occurs via a small subset of possible functionally-relevant substitutions to Na⁺,K⁺-ATPase. Surprisingly, we find that frequently observed adaptive substitutions at two sites, 111 and 122, are lethal when homozygous and adult heterozygotes exhibit dominant neural dysfunction. We identify a phylogenetically correlated substitution, A119S, that partially ameliorates the deleterious effects of substitutions at 111 and 122. Despite contributing little to cardiac glycoside-insensitivity in vitro, A119S, like substitutions at 111 and 122, substantially increases adult survivorship upon cardiac glycoside exposure. Our results demonstrate the importance of epistasis in constraining adaptive paths. Moreover, by revealing distinct effects of substitutions in vitro and in vivo, our results underscore the importance of evaluating the fitness of adaptive substitutions and their interactions in whole organisms.

Species-level predation network uncovers high prey specificity in a Neotropical army ant community

Army ants are among the top arthropod predators and considered keystone species in tropical ecosystems. During daily mass raids with many thousand workers, army ants hunt live prey, likely exerting strong top-down control on prey species. Many tropical sites exhibit a high army ant species diversity (>20 species), suggesting that sympatric species partition the available prey niches. However, whether and to what extent this is achieved has not been intensively studied yet. We therefore conducted a large-scale diet survey of a community of surface-raiding army ants at La Selva Biological Station in Costa Rica. We systematically collected 3,262 prey items from eleven army ant species (genera Eciton, Nomamyrmex and Neivamyrmex). Prey items were classified as ant prey or non-ant prey. The prey nearly exclusively consisted of other ants (98%), and most booty was ant brood (87%). Using morphological characters and DNA barcoding, we identified a total of 1,103 ant prey specimens to the species level. One hundred twenty-nine ant species were detected among the army ant prey, representing about 30% of the known local ant diversity. Using weighted bipartite network analyses, we show that prey specialization in army ants is unexpectedly high and prey niche overlap very small. Besides food niche differentiation, we uncovered a spatiotemporal niche differentiation in army ant raid activity. We discuss competition-driven multidimensional niche differentiation and predator-prey arms races as possible mechanisms underlying prey specialization in army ants. By combining systematic prey sampling with species-level prey identification and network analyses, our integrative approach can guide future research by portraying how predator-prey interactions in complex communities can be reliably studied, even in cases where morphological prey identification is infeasible.

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

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

Large-effect mutations generate trade-off between predatory and locomotor ability during arms race coevolution with deadly prey

Adaptive evolution in response to one selective challenge may disrupt other important aspects of performance. Such evolutionary trade‐offs are predicted to arise in the process of local adaptation, but it is unclear if these phenotypic compromises result from the antagonistic effects of simple amino acid substitutions. We tested for trade‐offs associated with beneficial mutations that confer tetrodotoxin (TTX) resistance in the voltage‐gated sodium channel (Na_V1.4) in skeletal muscle of the common garter snake (Thamnophis sirtalis). Separate lineages in California and the Pacific Northwest independently evolved TTX‐resistant changes to the pore of Na_V1.4 as a result of arms race coevolution with toxic prey, newts of the genus Taricha. Snakes from the California lineage that were homozygous for an allele known to confer large increases in toxin resistance (Na_V1.4^LVNV) had significantly reduced crawl speed compared to individuals with the ancestral TTX‐sensitive channel. Heterologous expression of native snake Na_V1.4 proteins demonstrated that the same Na_V1.4^LVNV allele confers a dramatic increase in TTX resistance and a correlated decrease in overall channel excitability. Our results suggest the same mutations that accumulate during arms race coevolution and beneficially interfere with toxin‐binding also cause changes in electrophysiological function of the channel that may affect organismal performance. This trade‐off was only evident in the predator lineage where coevolution has led to the most extreme resistance phenotype, determined by four critical amino acid substitutions. If these biophysical changes also translate to a fitness cost—for example, through the inability of T. sirtalis to quickly escape predators—then pleiotropy at this single locus could contribute to observed variation in levels of TTX resistance across the mosaic landscape of coevolution.

Comparison of Strategies to Overcome Drug Resistance: Learning from Various Kingdoms

Drug resistance, especially antibiotic resistance, is a growing threat to human health. To overcome this problem, it is significant to know precisely the mechanisms of drug resistance and/or self-resistance in various kingdoms, from bacteria through plants to animals, once more. This review compares the molecular mechanisms of the resistance against phycotoxins, toxins from marine and terrestrial animals, plants and fungi, and antibiotics. The results reveal that each kingdom possesses the characteristic features. The main mechanisms in each kingdom are transporters/efflux pumps in phycotoxins, mutation and modification of targets and sequestration in marine and terrestrial animal toxins, ABC transporters and sequestration in plant toxins, transporters in fungal toxins, and various or mixed mechanisms in antibiotics. Antibiotic producers in particular make tremendous efforts for avoiding suicide, and are more flexible and adaptable to the changes of environments. With these features in mind, potential alternative strategies to overcome these resistance problems are discussed. This paper will provide clues for solving the issues of drug resistance.

Interacting amino acid replacements allow poison frogs to evolve epibatidine resistance

Animals that wield toxins face self-intoxication. Poison frogs have a diverse arsenal of defensive alkaloids that target the nervous system. Among them is epibatidine, a nicotinic acetylcholine receptor (nAChR) agonist that is lethal at microgram doses. Epibatidine shares a highly conserved binding site with acetylcholine, making it difficult to evolve resistance yet maintain nAChR function. Electrophysiological assays of human and frog nAChR revealed that one amino acid replacement, which evolved three times in poison frogs, decreased epibatidine sensitivity but at a cost of acetylcholine sensitivity. However, receptor functionality was rescued by additional amino acid replacements that differed among poison frog lineages. Our results demonstrate how resistance to agonist toxins can evolve and that such genetic changes propel organisms toward an adaptive peak of chemical defense.

Step-wise evolution of complex chemical defenses in millipedes: a phylogenomic approach

With fossil representatives from the Silurian capable of respiring atmospheric oxygen, millipedes are among the oldest terrestrial animals, and likely the first to acquire diverse and complex chemical defenses against predators. Exploring the origin of complex adaptive traits is critical for understanding the evolution of Earth's biological complexity, and chemical defense evolution serves as an ideal study system. The classic explanation for the evolution of complexity is by gradual increase from simple to complex, passing through intermediate "stepping stone" states. Here we present the first phylogenetic-based study of the evolution of complex chemical defenses in millipedes by generating the largest genomic-based phylogenetic dataset ever assembled for the group. Our phylogenomic results demonstrate that chemical complexity shows a clear pattern of escalation through time. New pathways are added in a stepwise pattern, leading to greater chemical complexity, independently in a number of derived lineages. This complexity gradually increased through time, leading to the advent of three distantly related chemically complex evolutionary lineages, each uniquely characteristic of each of the respective millipede groups.

Venomics: integrative venom proteomics and beyond*

Venoms are integrated phenotypes that evolved independently in, and are used for predatory and defensive purposes by, a wide phylogenetic range of organisms. The same principles that contribute to the evolutionary success of venoms, contribute to making the study of venoms of great interest in such diverse fields as evolutionary ecology and biotechnology. Evolution is profoundly contingent, and nature also reinvents itself continuosly. Changes in a complex phenotypic trait, such as venom, reflect the influences of prior evolutionary history, chance events, and selection. Reconstructing the natural history of venoms, particularly those of snakes, which will be dealt with in more detail in this review, requires the integration of different levels of knowledge into a meaningful and comprehensive evolutionary framework for separating stochastic changes from adaptive evolution. The application of omics technologies and other disciplines have contributed to a qualitative and quantitative advance in the road map towards this goal. In this review we will make a foray into the world of animal venoms, discuss synergies and complementarities of the different approaches used in their study, and identify current bottlenecks that prevent inferring the evolutionary mechanisms and ecological constraints that molded snake venoms to their present-day variability landscape.

Snake Genome Sequencing: Results and Future Prospects

Snake genome sequencing is in its infancy—very much behind the progress made in sequencing the genomes of humans, model organisms and pathogens relevant to biomedical research, and agricultural species. We provide here an overview of some of the snake genome projects in progress, and discuss the biological findings, with special emphasis on toxinology, from the small number of draft snake genomes already published. We discuss the future of snake genomics, pointing out that new sequencing technologies will help overcome the problem of repetitive sequences in assembling snake genomes. Genome sequences are also likely to be valuable in examining the clustering of toxin genes on the chromosomes, in designing recombinant antivenoms and in studying the epigenetic regulation of toxin gene expression.