Concerted evolution reveals co-adapted amino acid substitutions in Na+K+-ATPase of frogs that prey on toxic toads

Evaluating the efficacy of protein quantification methods on membrane proteins

Protein quantification is an important tool for a wide range of biological applications. The most common broadscale methods include the Lowry, bicinchoninic acid (BCA), and Coomassie Bradford assays. Despite their wide applicability, the mechanisms of action imply that these methods may not be ideal for large transmembrane proteins due to the proteins' integration in the plasma membrane. Here, we investigate this problem by assessing the efficacy and applicability of these three common protein quantification methods on a candidate transmembrane protein - the Na,K-ATPase (NKA). We compared these methods to an ELISA, which we newly developed and describe here for the quantification of NKA. The use of a relative standard curve allows this ELISA to be easily adapted to other proteins and across the animal kingdom. Our results revealed that the three conventional methods significantly underestimate the concentration of NKA compared to the ELISA. Further, by applying the protein concentrations determined by the different methods to in vitro assays, we found that variation in the resulting data was consistently low when the assay reactions were prepared based on concentrations determined from the ELISA. Thus, when target protein concentrations vary across samples, the conventional quantification methods cannot produce reliable results in downstream applications. In contrast, the ELISA we describe here consistently provides robust results.

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.

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.

Involvement of the alpha-subunit N-terminus in the mechanism of the Na+,K+-ATPase

Previous studies have shown that cytoplasmic K+ release and the associated E2 → E1 conformational change of the Na+,K+-ATPase is a major rate-determining step of the enzyme's ion pumping cycle and hence a prime site of acute regulatory intervention. From the ionic strength dependence of the enzyme's distribution between the E2 and E1 states, it has also been found that E2 is stabilized by an electrostatic attraction. Any disruption of this electrostatic attraction would, thus, have profound effects on the rate of ion pumping. The aim of this paper is to identify the location of this interaction. Using enhanced-sampling molecular dynamics simulations with a predicted N-terminal structure added to the X-ray crystal structure of the Na+,K+-ATPase, a previously postulated salt bridge between Lys32 and Glu233 (rat sequence numbering) of the enzyme's α-subunit can be excluded. The residues never approach closely enough to form a salt bridge. In contrast, strong interactions with anionic lipid head groups were seen. To investigate the possibility of a protein-lipid interaction experimentally, the surface charge density of Na+,K+-ATPase-containing membrane fragments was estimated from zeta potential measurements to be 0.019 (± 0.001) C m-2. This is in good agreement with the charge density previously determined to be responsible for stabilization of the E2 state of 0.023 (± 0.009) C m-2 and the membrane charge density estimated here from published electron-microscopic images of 0.018C m-2. The results are, therefore, consistent with an interaction of the Na+,K+-ATPase α-subunit N-terminus with negatively-charged lipid head groups of the neighbouring cytoplasmic membrane surface as the origin of the electrostatic interaction stabilising the E2 state.

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.

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.

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.

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.

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

Individual amino acids within a protein work in concert to produce a functionally coherent structure that must be maintained as a protein diverges over time. Given this structure-function relationship, we expect the effects of new mutations to depend on amino acid states at other sites throughout the protein (i.e., background dependence) and that identical mutations will have more similar effects in more closely-related species, for which orthologous proteins will be less diverged. We tested this hypothesis by performing protein-engineering experiments on ATP1A, a protein that mediates resistance to toxins known as cardiotonic steroids (CTS), to reveal the extent of background-dependence across representative tetrapods. We find that, while the effects of mutations at two key sites implicated in CTS-resistance are indeed often background-dependent, the magnitude of these effects does not correlate with overall levels of ATP1A divergence. Our results instead suggest that background-dependent effects are determined by amino acid states at a small number of sites throughout the protein. Evolutionary constraints imposed by relatively few sites may explain the frequent occurrence of identical or similar CTS-resistance substitutions among ATP1A proteins of highly divergent animals (See S1 Text for Spanish translation).

Gene expression phylogenies and ancestral transcriptome reconstruction resolves major transitions in the origins of pregnancy

Structural and physiological changes in the female reproductive system underlie the origins of pregnancy in multiple vertebrate lineages. In mammals, the glandular portion of the lower reproductive tract has transformed into a structure specialized for supporting fetal development. These specializations range from relatively simple maternal nutrient provisioning in egg-laying monotremes to an elaborate suite of traits that support intimate maternal-fetal interactions in Eutherians. Among these traits are the maternal decidua and fetal component of the placenta, but there is considerable uncertainty about how these structures evolved. Previously, we showed that changes in uterine gene expression contributes to several evolutionary innovations during the origins of pregnancy (Mika et al., 2021b). Here, we reconstruct the evolution of entire transcriptomes ('ancestral transcriptome reconstruction') and show that maternal gene expression profiles are correlated with degree of placental invasion. These results indicate that an epitheliochorial-like placenta evolved early in the mammalian stem-lineage and that the ancestor of Eutherians had a hemochorial placenta, and suggest maternal control of placental invasiveness. These data resolve major transitions in the evolution of pregnancy and indicate that ancestral transcriptome reconstruction can be used to study the function of ancestral cell, tissue, and organ systems.

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.

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.

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.

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