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.

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.

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.

When does the female bias arise? Insights from the sex determination cascade of a flea beetle with a strongly skewed sex ratio.

Reproduction-manipulating bacteria like Wolbachia can shift sex ratios in insects towards females, but skewed sex ratios may also arise from genetic conflicts. The flea beetle Altica lythri harbors three main mtDNA strains that are coupled to three different Wolbachia infections. Depending on the mtDNA types, the females produce either offspring with a balanced sex ratio or exclusively daughters. To obtain markers that can monitor when sex bias arises in the beetle's ontogeny, we elucidated the sex determination cascade of A. lythri. We established a RT-PCR method based on length variants of dsx (doublesex) transcripts to determine the sex of morphologically indistinguishable eggs and larvae. In females of one mtDNA type (HT1/HT1*) known to produce only daughters, male offspring were already missing at the egg stage while for females of another type (HT2), the dsx splice variants revealed a balanced sex ratio among eggs and larvae. Our data suggest that the sex determination cascade in A. lythri is initiated by maternally transmitted female-specific tra (transformer) mRNA as primary signal. This tra mRNA seems to be involved in a positive feedback loop that maintains the production of the female splice variant, as known for female offspring in Tribolium castaneum. The translation of the maternally transmitted female tra mRNA must be inhibited in male offspring, but the underlying primary genetic signal remains to be identified. We discuss which differences between the mtDNA types can influence sex determination and lead to the skewed sex ratio of HT1.

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.

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.

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

Functional evidence supports adaptive plant chemical defense along a geographical cline.

Environmental clines in organismal defensive traits are usually attributed to stronger selection by enemies at lower latitudes or near the host's range center. Nonetheless, little functional evidence has supported this hypothesis, especially for coevolving plants and herbivores. We quantified cardenolide toxins in seeds of 24 populations of common milkweed (Asclepias syriaca) across 13 degrees of latitude, revealing a pattern of increasing cardenolide concentrations toward the host's range center. The unusual nitrogen-containing cardenolide labriformin was an exception and peaked at higher latitudes. Milkweed seeds are eaten by specialist lygaeid bugs that are even more tolerant of cardenolides than the monarch butterfly, concentrating most cardenolides (but not labriformin) from seeds into their bodies. Accordingly, whether cardenolides defend seeds against these specialist bugs is unclear. We demonstrate that Oncopeltus fasciatus (Lygaeidae) metabolized two major compounds (glycosylated aspecioside and labriformin) into distinct products that were sequestered without impairing growth. We next tested several isolated cardenolides in vitro on the physiological target of cardenolides (Na+/K+-ATPase); there was little variation among compounds in inhibition of an unadapted Na+/K+-ATPase, but tremendous variation in impacts on that of monarchs and Oncopeltus. Labriformin was the most inhibitive compound tested for both insects, but Oncopeltus had the greater advantage over monarchs in tolerating labriformin compared to other compounds. Three metabolized (and stored) cardenolides were less toxic than their parent compounds found in seeds. Our results suggest that a potent plant defense is evolving by natural selection along a geographical cline and targets specialist herbivores, but is met by insect tolerance, detoxification, and sequestration.

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.

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

Although gene duplication is an important source of evolutionary innovation, the functional divergence of duplicates can be opposed by ongoing gene conversion between them. Here, we report on the evolution of a tandem duplication of Na+,K+-ATPase subunit α1 (ATP1A1) shared by frogs in the genus Leptodactylus, a group of species that feeds on toxic toads. One ATP1A1 paralog evolved resistance to toad toxins although the other retained ancestral susceptibility. Within species, frequent non-allelic gene conversion homogenized most of the sequence between the two copies but was counteracted by strong selection on 12 amino acid substitutions that distinguish the two paralogs. Protein-engineering experiments show that two of these substitutions substantially increase toxin resistance, whereas the additional 10 mitigate their deleterious effects on ATPase activity. Our results reveal how examination of neo-functionalized gene duplicate evolution can help pinpoint key functional substitutions and interactions with the genetic backgrounds on which they arise.

ABCB transporters in a leaf beetle respond to sequestered plant toxins.

Phytophagous insects can tolerate and detoxify toxic compounds present in their host plants and have evolved intricate adaptations to this end. Some insects even sequester the toxins for their defence. This necessitates specific mechanisms, especially carrier proteins that regulate uptake and transport to specific storage sites or protect sensitive tissues from noxious compounds. We identified three ATP-binding cassette subfamily B (ABCB) transporters from the transcriptome of the cardenolide-sequestering leaf beetle Chrysochus auratus and analysed their functional role in the sequestration process. These were heterologously expressed and tested for their ability to interact with various potential substrates: verapamil (standard ABCB substrate), the cardenolides digoxin (commonly used), cymarin (present in the species's host plant) and calotropin (present in the ancestral host plants). Verapamil stimulated all three ABCBs and each was activated by at least one cardenolide, however, they differed as to which they were activated by. While the expression of the most versatile transporter fits with a protective role in the blood-brain barrier, the one specific for cymarin shows an extreme abundance in the elytra, coinciding with the location of the defensive glands. Our data thus suggest a key role of ABCBs in the transport network needed for cardenolide sequestration.

Biased predation could promote convergence yet maintain diversity within Müllerian mimicry rings of Oreina leaf beetles.

Müllerian mimicry is a classic example of adaptation, yet Müller's original theory does not account for the diversity often observed in mimicry rings. Here, we aimed to assess how well classical Müllerian mimicry can account for the colour polymorphism found in chemically defended Oreina leaf beetles by using field data and laboratory assays of predator behaviour. We also evaluated the hypothesis that thermoregulation can explain diversity between Oreina mimicry rings. We found that frequencies of each colour morph were positively correlated among species, a critical prediction of Müllerian mimicry. Predators learned to associate colour with chemical defences. Learned avoidance of the green morph of one species protected green morphs of another species. Avoidance of blue morphs was completely generalized to green morphs, but surprisingly, avoidance of green morphs was less generalized to blue morphs. This asymmetrical generalization should favour green morphs: indeed, green morphs persist in blue communities, whereas blue morphs are entirely excluded from green communities. We did not find a correlation between elevation and coloration, rejecting thermoregulation as an explanation for diversity between mimicry rings. Biased predation could explain within-community diversity in warning coloration, providing a solution to a long-standing puzzle. We propose testable hypotheses for why asymmetric generalization occurs, and how predators maintain the predominance of blue morphs in a community, despite asymmetric generalization.

Genome editing retraces the evolution of toxin resistance in the monarch butterfly.

Identifying the genetic mechanisms of adaptation requires the elucidation of links between the evolution of DNA sequence, phenotype, and fitness1. Convergent evolution can be used as a guide to identify candidate mutations that underlie adaptive traits2-4, and new genome editing technology is facilitating functional validation of these mutations in whole organisms1,5. We combined these approaches to study a classic case of convergence in insects from six orders, including the monarch butterfly (Danaus plexippus), that have independently evolved to colonize plants that produce cardiac glycoside toxins6-11. Many of these insects evolved parallel amino acid substitutions in the α-subunit (ATPα) of the sodium pump (Na+/K+-ATPase)7-11, the physiological target of cardiac glycosides12. Here we describe mutational paths involving three repeatedly changing amino acid sites (111, 119 and 122) in ATPα that are associated with cardiac glycoside specialization13,14. We then performed CRISPR-Cas9 base editing on the native Atpα gene in Drosophila melanogaster flies and retraced the mutational path taken across the monarch lineage11,15. We show in vivo, in vitro and in silico that the path conferred resistance and target-site insensitivity to cardiac glycosides16, culminating in triple mutant 'monarch flies' that were as insensitive to cardiac glycosides as monarch butterflies. 'Monarch flies' retained small amounts of cardiac glycosides through metamorphosis, a trait that has been optimized in monarch butterflies to deter predators17-19. The order in which the substitutions evolved was explained by amelioration of antagonistic pleiotropy through epistasis13,14,20-22. Our study illuminates how the monarch butterfly evolved resistance to a class of plant toxins, eventually becoming unpalatable, and changing the nature of species interactions within ecological communities2,6-11,15,17-19.

New ways to acquire resistance: imperfect convergence in insect adaptations to a potent plant toxin.

Evolution of insensitivity to the toxic effects of cardiac glycosides has become a model in the study of convergent evolution, as five taxonomic orders of insects use the same few similar amino acid substitutions in the otherwise highly conserved Na,K-ATPase α. We show here that insensitivity in pyrgomorphid grasshoppers evolved along a slightly divergent path. As in other lineages, duplication of the Na,K-ATPase α gene paved the way for subfunctionalization: one copy maintains the ancestral, sensitive state, while the other copy is resistant. Nonetheless, in contrast with all other investigated insects, the grasshoppers' resistant copy shows length variation by two amino acids in the first extracellular loop, the main part of the cardiac glycoside-binding pocket. RT-qPCR analyses confirmed that this copy is predominantly expressed in tissues exposed to the toxins, while the ancestral copy predominates in the nervous tissue. Functional tests with genetically engineered Drosophila Na,K-ATPases bearing the first extracellular loop of the pyrgomorphid genes showed the derived form to be highly resistant, while the ancestral state is sensitive. Thus, we report convergence in gene duplication and in the gene targets for toxin insensitivity; however, the means to the phenotypic end have been novel in pyrgomorphid grasshoppers.

Robust reference gene design and validation for expression studies in the large milkweed bug, Oncopeltus fasciatus, upon cardiac glycoside stress.

Despite its history as a developmental and evolutionary model organism, gene expression analysis in the large milkweed bug, Oncopeltus fasciatus, has rarely been explored using quantitative real-time PCR. The strength of this method depends greatly on the endogenous controls used for normalization, which are lacking for the milkweed bug system. Here, to fill in this gap in our knowledge, we validated the stability of a set of ten candidate reference genes identified from the O. fasciatus transcriptome, and did so upon exposure to a dietary toxin, a cardiac glycoside, and across four different exposure periods. To increase robustness against gDNA contaminants, genome resources were used to design intron-bridging primers. A comprehensive stability validation by the Bestkeeper, Normfinder, geNorm and comparative ΔCt methods identified ef1a and tubulin as the most stable genes across treatments and time points, whereas 18S rRNA was the most unstable. However, accounting for the temporal scale indicated that time point confined normalizers might enable higher quantification accuracy for treatment comparison. Overall this study demonstrates: (i) a robust RT-qPCR primer design approach is possible for non-model organisms where genome annotation is often incomplete, and (ii) the importance of detailed reference gene stability exploration in multifactorial experimental designs.

Ontogeny of Defensive Chemistry in Longitarsus Flea Beetles (Coleoptera, Chrysomelidae): More Protection for the Vulnerable Stages?

Several species of the flea beetles genus Longitarsus sequester pyrrolizidine alkaloids (PAs) from their host plants. Previous data demonstrated that PAs may be transferred from root-feeding larvae into the adult beetles. Here we compared the patterns and concentrations found in larvae and pupae of L. anchusae and L. echii with those of the roots of their respective hosts, Symphytum officinale and Echium vulgare (Boraginaceae). PA patterns and concentrations in the roots were complex and variable, whereas those in the larvae and pupae were simpler and more constant. In L. anchusae, intermedine and lycopsamine were the dominant PAs even if they could not be detected in the roots. In L. echii simpler, hydrolized PAs prevailed. Overall, the concentrations of total PAs of larvae and pupae were significantly higher than those of the roots the larvae had been feeding on. Larvae and pupae of both species also had considerably higher PA concentrations than determined previously for field collected beetles. Possibly the rather immobile juvenile stages enjoy a better protection by higher PA concentrations. On the other hand, we could not detect PAs in eggs of either species, indicating that transmission of appreciable amounts of PAs from mother to offspring does not occur.

The function and evolutionary significance of a triplicated Na,K-ATPase gene in a toxin-specialized insect.

BACKGROUND: The Na,K-ATPase is a vital animal cell-membrane protein that maintains the cell's resting potential, among other functions. Cardenolides, a group of potent plant toxins, bind to and inhibit this pump. The gene encoding the α-subunit of the pump has undergone duplication events in some insect species known to feed on plants containing cardenolides. Here we test the function of these duplicated gene copies in the cardenolide-adapted milkweed bug, Oncopeltus fasciatus, which has three known copies of the gene: α1A, α1B and α1C. RESULTS: Using RT-qPCR analyses we demonstrate that the α1C is highly expressed in neural tissue, where the pump is generally thought to be most important for neuron excitability. With the use of in vivo RNAi in adult bugs we found that α1C knockdowns suffered high mortality, where as α1A and α1B did not, supporting that α1C is most important for effective ion pumping. Next we show a role for α1A and α1B in the handling of cardenolides: expression results find that both copies are primarily expressed in the Malpighian tubules, the primary insect organ responsible for excretion, and when we injected either α1A or α1B knockdowns with cardenolides this proved fatal (whereas not in controls). CONCLUSIONS: These results show that the Na,K-ATPα gene-copies have taken on diverse functions. Having multiple copies of this gene appears to have allowed the newly arisen duplicates to specialize on resistance to cardenolides, whereas the ancestral copy of the pump remains comparatively sensitive, but acts as a more efficient ion carrier. Interestingly both the α1A and α1B were required for cardenolide handling, suggesting that these two copies have separate and vital functions. Gene duplications of the Na,K-ATPase thus represent an excellent example of subfunctionalization in response to a new environmental challenge.

Substitutions in the cardenolide binding site and interaction of subunits affect kinetics besides cardenolide sensitivity of insect Na,K-ATPase.

Substitutions within the cardenolide target site of several insects' Na,K-ATPase α-subunits may confer resistance against toxic cardenolides. However, to which extent these substitutions alter the Na,K-ATPase's kinetic properties and how they interact with different ß-subunits is not clear. The cardenolide-adapted milkweed bug Oncopeltus fasciatus possesses three paralogs of the α-subunit (A, B, and C) that differ in number and identity of resistance-conferring substitutions. We introduced these substitutions into the α-subunit of Drosophila melanogaster and combined them with the ß-subunits Nrv2.2 and Nrv3. The substitutions Q111T-N122H-F786N-T797A (A-copy mimic) and Q111T-N122H-F786N (B-copy mimic) mediated high insensitivity to ouabain, yet they drastically lowered ATPase activity. Remarkably, the identity of the ß-subunit was decisive and all α-subunits were less active when combined with Nrv3 than when combined with Nrv2.2. Both the substitutions and the co-expressed ß-subunit strongly affected the enyzme's affinity for Na+ and K+. Na+ affinity was considerably higher for all enzymes expressed with nrv3 while expression with nrv2.2 mostly increased K+ affinity. Our results provide the first evidence that resistance against cardenolides comes at the cost of significantly altered kinetic properties of the Na,K-ATPase. The ß-subunit can strongly modulate these properties but cannot fully compensate for the effect of the substitutions.

Convergently Evolved Toxic Secondary Metabolites in Plants Drive the Parallel Molecular Evolution of Insect Resistance.

Natural selection imposed by natural toxins has led to striking levels of convergent evolution at the molecular level. Cardiac glycosides represent a group of plant toxins that block the Na,K-ATPase, a vital membrane protein in animals. Several herbivorous insects have convergently evolved resistant Na,K-ATPases, and in some species, convergent gene duplications have also arisen, likely to cope with pleiotropic costs of resistance. To understand the genetic basis and predictability of these adaptations, we studied five independent lineages of leaf-mining flies (Diptera: Agromyzidae). These flies have colonized host plants in four botanical families that convergently evolved cardiac glycosides of two structural types: cardenolides and bufadienolides. We compared each of six fly species feeding on such plants to a phylogenetically related but nonadapted species. Irrespective of the type of cardiac glycoside in the host plant, five out of six exposed species displayed substitutions in the cardiac glycoside-binding site of the Na,K-ATPase that were previously described in other insect orders; in only one species was the gene duplicated. In vitro assays of nervous tissue extractions confirmed that the substitutions lead to increased resistance of the Na,K-ATPase. Our results demonstrate that target site insensitivity of Na,K-ATPase is a common response to dietary cardiac glycosides leading to highly predictable amino acid changes; nonetheless, convergent evolution of gene duplication for this multifunctional enzyme appears more constrained.

Toad toxin-resistant snake (Thamnophis elegans) expresses high levels of mutant Na+/K+-ATPase mRNA in cardiac muscle.

Toads are chemically defended by bufadienolides, which are lethal to most predators. These toxins exert their lethal effects by binding to and disabling the Na+/K+-ATPases of cell membranes. Many species of snakes exhibit resistance to the effects of bufadienolides due to target-site insensitivity of the Na+/K+-ATPase. Mutations that confer resistance have previously been identified in ATP1a3, the gene that codes for the Na+/K+-ATPase α-3 paralog. We have found that this mutant gene is expressed at a significantly elevated level in heart tissue compared to gut, kidney, and liver of the bufadienolide-resistant snake, Thamnophis elegans. Furthermore, we found that exposure to bufadienolides elicits a significant increase in the expression levels of ATP1a3 in the heart, but not in the kidneys, liver, or gut 1h after exposure. We suggest that upregulation of ATP1a3 in the heart plays an important role in the physiological processes involved in tolerance of bufadienolides among genetically resistant snakes.