Search

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.

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.

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.

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 Oncopeltu s. 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.

The recurrent evolution of resistance to cardiotonic steroids (CTS) across diverse animals most frequently involves convergent amino-acid substitutions to the H1-H2 extracellular loop of Na+,K+-ATPase (NKA). Previous work established 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 is not known. Using protein engineering, we show that these insertions have distinct effects on CTS resistance of NKA in the two lineages that strongly depend on intramolecular interactions with other residues. Removing the insertion in the chinchilla lineage unexpectedly increases CTS resistance and decreases NKA activity. In the sandgrouse lineage, the insertion works in concert with the substitution Q111R to increase CTS resistance while maintaining wild-type ATPase activity levels. Molecular docking simulations provide additional insight into the biophysical mechanisms responsible for the context-specific 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.

Insects live in a dangerous world and may fall prey to a wide variety of predators, encompassing multiple taxa. As a result, selection may favour defences that are effective against multiple predator types, or target-specific defences that can reduce predation risk from particular groups of predators. Given the variation in sensory systems and hunting tactics, in particular between vertebrate and invertebrate predators, it is not always clear whether defences, such as chemical defences, that are effective against one group will be so against another. Despite this, the majority of research to date has focused on the role of a single predator species when considering the evolution of defended prey. Here we test the effectiveness of the chemical defences of the wood tiger moth, a species previously shown to have defensive chemicals targeted towards ants, against a common invertebrate predator: spiders. We presented both live moths and artificial prey containing their defensive fluids to female Trichonephila senegalensis and recorded their reactions. We found that neither of the moth’s two defensive fluids were able to repel the spiders, and confirmed that methoxypyrazines, a major component of the defences of both the wood tiger moth and many insect species, are ineffective against web-building spiders. Our results highlight the variability between predator taxa in their susceptibility to chemical defences, which can in part explain the vast variation in these chemicals seen in insects, and the existence of multiple defences in a single species.

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.

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.Significance StatementIndividual amino acid residues within a protein work in concert to produce a functionally coherent structure that must be maintained even as orthologous proteins in different species diverge over time. Given this dependence, we expect identical mutations to have more similar effects on protein function in more closely related species. We tested this hypothesis by performing protein-engineering experiments on ATP1A, an enzyme mediating target-site insensitivity to cardiotonic steroids (CTS) in diverse animals. These experiments reveal that the phenotypic effects of substitutions can sometimes be background-dependent, but also that the magnitude of these phenotypic effects does not correlate with overall levels of ATP1A sequence divergence. Our results suggest that epistatic constraints are determined by states at a small number of sites, potentially explaining the frequent convergent CTS resistance substitutions among Na,K-ATPases of highly divergent taxa.

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.

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

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 beetleChrysochus auratusand 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.

BackgroundCardiac glycosides are known to fatally inhibit the Na,K-ATPase throughout the animal kingdom. Several animals, however, evolved target-site insensitivity by substitutions in the otherwise highly conserved cardiac glycoside binding pocket located on the Na,K-ATPase α-subunit. The minimal functional enzyme consist of an α- and a β-subunit, the latter considered mainly as a chaperone responsible for correct folding and membrane integration. We here analyze resistance to cardiac glycosides and kinetic properties of different Na,K-ATPase α/β-combinations of the large milkweed bug, Oncopeltus fasciatus. These insects have adapted to high concentrations of cardiac glycosides in their food plants via several rounds of Na,K-ATPase gene duplications followed by differential resistance conferring substitutions and subfunctionalization of the enzymes.ResultsTo investigate their characteristics we expressed nine combinations of O. fasciatus Na,K-ATPase α/β-sunbunits (three each) in Sf9 cells and tested them with two structurally distinct cardiac glycosides, calotropin, a host plant compound, and ouabain, a commonly used toxin. Differences in the number and identity of amino acid substitutions in the cardiac glycoside binding site resulted in large differences in activity and toxin resistance of the three α-subunits. The enzymes’ kinetics were also influenced by the β-subunits leading to increased activities (αCβ3) or altered resistances. The host plant toxin calotropin proved to be a much more potent inhibitor than ouabain for the phylogenetically oldest αC based enzymes. This effect was compensated for in the αB and αA based enzymes with αAβ1 having higher resistance against calotropin than against ouabain.ConclusionThe originally higher inhibitory potency of the host compound calotropin supports a coevolutionary escalation of plant defenses and herbivore tolerance mechanisms. For the bugs the possession of multiple paralogs improved adaptation to plant toxins in a stepwise manner and mitigates pleiotropic effects by a compromise between ion pumping activity and resistance.

Gene duplication is an important source of evolutionary innovation, but the adaptive division-of-labor between duplicates can be opposed by ongoing gene conversion between them. Here we document a tandem duplication of Na+,K+-ATPase subunit α1 (ATP1A1) shared by frogs in the genusLeptodactylus, a group of species that feeds on toxic toads. One ATP1A1 paralog evolved resistance to toad toxins while the other paralog retained ancestral susceptibility. We show that the twoLeptodactylusparalogs are distinguished by 12 amino acid substitutions that were maintained by strong selection that counteracted the homogenizing effect of gene conversion. Protein-engineering experiments show that two major-effect substitutions confer toxin resistance, whereas the 10 additional substitutions mitigate deleterious pleiotropic effects on enzyme function. Our results highlight how trans-specific, neofunctionalized gene duplicates can provide unique insights into interactions between adaptive substitutions and the genetic backgrounds on which they arise.One Sentence SummarySelection counteracts gene conversion to maintain an adaptive division-of-labor between tandemly duplicated genes.

The larvae of several species in the hawk moth genus Hyles, including H. euphorbiae, feed on plants of the genus Euphorbia containing phorbol esters and are insensitive to addition of the standard phorbol ester, tetradecanoyl-phorbol-13-acetate (TPA) to their artificial diet. Specialised non-Euphorbia feeding larvae were sensitive, whereas polyphagous ones were insensitive if their natural food plant spectrum also includes Euphorbia. Larvae of Hippotion celerio, an out-group species with polyphagous larvae not using Euphorbia as food plants, were sensitive. A highly conserved sequence of the TPA binding site of the protein kinase C in H. euphorbiae and Hippotion celerio demonstrates that intoxication by phorbol esters is not avoided by preventing target binding. H. euphorbiae larvae showed no vitality loss after chemical destruction of their peritrophic matrix and subsequent TPA treatment. TPA fed larvae that had putatively piperonyl butoxide (PBO)-inhibited cytochrome P450 enzymes also showed no deficits, indicating that this is not the only detoxification pathway in H. euphorbiae. Based on these qualitative results, we postulate that proto-Hyles was polyphagous with the ability to use Euphorbia as food plants. The most ancestral Hyles species presumably remained polyphagous, with the ability to switch to the toxic food plants if necessary. Within the youngest, Palearctic radiation, the species specialised on non-Euphorbia plants subsequently lost detoxification abilities, albeit not to 100%. The species in South America, the origin of Hyles, can detoxify TPA, indicating this ability to be a plesiomorphic character state within the genus that enabled its adaptive radiation in the Palearctic, the distribution and diversity centre of Euphorbia.

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.

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.

In the struggle against dietary toxins, insects are known to employ target site insensitivity, metabolic detoxification, and transporters that shunt away toxins. Specialized insects across six taxonomic orders feeding on cardenolide-containing plants have convergently evolved target site insensitivity via specific amino acid substitutions in the Na/K-ATPase. Nonetheless, in vitro pharmacological experiments have suggested a role for multidrug transporters (Mdrs) and organic anion transporting polypeptides (Oatps), which may provide a basal level of protection in both specialized and non-adapted insects. Because the genes coding for these proteins are evolutionarily conserved and in vivo genetic evidence in support of this hypothesis is lacking, here we used wildtype and mutant Drosophila melanogaster (Drosophila) in capillary feeder (CAFE) assays to quantify toxicity of three chemically diverse, medically relevant cardenolides. We examined multiple components of fitness, including mortality, longevity, and LD50, and found that, while the three cardenolides each stimulated feeding (i.e., no deterrence to the toxin), all decreased lifespan, with the most apolar cardenolide having the lowest LD50 value. Flies showed a clear non-monotonic dose response and experienced high levels of toxicity at the cardenolide concentration found in plants. At this concentration, both Mdr and Oatp knockout mutant flies died more rapidly than wildtype flies, and the mutants also experienced more adverse neurological effects on high-cardenolide-level diets. Our study further establishes Drosophila as a model for the study of cardenolide pharmacology and solidifies support for the hypothesis that multidrug and organic anion transporters are key players in insect protection against dietary cardenolides.

Mullerian mimicry is a classic example of adaptation, yet Muller's original theory does not account for the diversity often observed in mimicry rings. Here, we aimed to assess how well classical Mullerian 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 Mullerian 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.

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.

Herbivorous insects and their adaptations against plant toxins provide striking opportunities to investigate the genetic basis of traits involved in coevolutionary interactions. Target site insensitivity to cardenolides has evolved convergently across six orders of insects, involving identical substitutions in the Na,K-ATPase gene and repeated convergent gene duplications. The large milkweed bug, Oncopeltus fasciatus, has three copies of the Na,K-ATPase α-subunit gene that bear differing numbers of amino acid substitutions in the binding pocket for cardenolides. To analyze the effect of these substitutions on cardenolide resistance and to infer possible trade-offs in gene function, we expressed the cardenolide-sensitive Na,K-ATPase of Drosophila melanogaster in vitro and introduced four distinct combinations of substitutions observed in the three gene copies of O. fasciatus. With an increasing number of substitutions, the sensitivity of the Na,K-ATPase to a standard cardenolide decreased in a stepwise manner. At the same time, the enzyme's overall activity decreased significantly with increasing cardenolide resistance and only the least substituted mimic of the Na,K-ATPase α1C copy maintained activity similar to the wild-type enzyme. Our results suggest that the Na,K-ATPase copies in O. fasciatus have diverged in function, enabling specific adaptations to dietary cardenolides while maintaining the functionality of this critical ion carrier.

Species of the heteropteran subfamily Lygaeinae possess special subcuticular compartments to store cardiac glycosides, plant-derived defensive compounds, which they release upon predator attack. In adults of the large milkweed bug, Oncopeltus fasciatus, these storage compartments have previously been described as a modified integument, forming a fluid-filled dorsolateral space. Here we use three-dimensional imaging of serial histological sections and synchrotron radiation-based micro-computed tomography data to reveal the morphology of these storage compartments and the mechanisms used for the release of a cardiac glycoside-rich fluid upon attack. Our comparative analysis revealed that the morphology and release mechanism vary among the species investigated. By reconstructing these traits on a recent molecular phylogeny of the Lygaeinae, we demonstrate that the adaptations for the storage and release of cardiac glycosides have evolved in a stepwise manner.

Cardiac glycosides are a prime example of highly toxic plant secondary compounds, which block an essential transmembrane carrier in animals, the Na,K-ATPase. Nevertheless, over 100 insect species from diverse orders are known to feed on plants containing these compounds and in many cases these toxins are additionally sequestered without ill effect. We investigated whether the insects’ adaptations for handling cardiac glycosides are based on a single physiological mechanism or whether various strategies have evolved across groups. We analyzed gene sequences of the Na,K-ATPase α-subunit from cardiac glycoside-adapted insects and screened for amino-acid substitutions which could alter the affinity of the enzyme toward cardiac glycosides. In representatives from five insect orders, separated by over 300 million years of evolutionary divergence, we uncovered amino-acid substitutions at identical positions. Especially striking is the convergent substitution of a histidine for the conserved asparagine at position 122, which we report here for the first time in a sawfly, Monophadnus latus Costa (Hymenoptera: Tenthredinidae), and which was previously observed in the orders Lepidoptera, Coleoptera, Hemiptera, and Diptera. Prior in vitro expression and enzyme assays indicated that this substitution as well as combined substitutions with other residues result in a strongly increased cardenolide resistance of the Na,K-ATPase. The substitutions to threonine111 and histidine122 observed in M. latus are highly effective and were previously known only in lygaeid bugs. However, not all insects dealing with dietary cardenolides rely on target-site insensitivity as a mechanism of resistance. An impermeable gut or the exclusion of cardenolides from the nervous tissue with the greatest expression of Na,K-ATPase by the perineurium, the insect blood brain barrier, apparently represent alternative strategies. Immuno-histochemical data presented here support the existence of P-glycoprotein-like efflux transporters in insect gut membranes that might prevent the uptake of allelochemicals like cardenolides.

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. Electronic supplementary material The online version of this article (10.1186/s12862-017-1097-6) contains supplementary material, which is available to authorized users.

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.

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.

Specular reflection appears as a bright spot or highlight on any smooth glossy convex surface and is caused by a near mirror-like reflectance off the surface. Convex shapes always provide the ideal geometry for highlights, areas of very strong reflectance, regardless of the orientation of the surface or position of the receiver. Despite highlights and glossy appearance being common in chemically defended insects, their potential signalling function is unknown. We tested the role of highlights in warning colouration of a chemically defended, alpine leaf beetle, Oreina cacaliae. We reduced the beetles’ glossiness, hence their highlights, by applying a clear matt finish varnish on their elytra. We used blue tits as predators to examine whether the manipulation affected their initial latency to attack, avoidance learning and generalization of warning colouration. The birds learned to avoid both dull and glossy beetles but they initially avoided glossy prey more than dull prey. Interestingly, avoidance learning was generalized asymmetrically: birds that initially learned to avoid dull beetles avoided glossy beetles equally strongly, but not vice versa. We conclude that specular reflectance and glossiness can amplify the warning signal of O. cacaliae, augmenting avoidance learning, even if it is not critical for it.

The blowfly Lucilia bufonivora shows high host specificity for toads despite the host’s toxic skin secretion, which consists mainly of bufadienolides. These toxins are effective blockers of the Na+, K+-ATPase, an enzyme that is essential for many physiological processes in animals. Whereas common toad (Bufo bufo) toxins were identified in the larvae of the fly, few toxins were found in the pupae and empty puparia as trace amounts, while adult flies were entirely free of these toxic compounds. Similar results were obtained when larvae of generalist necrophagous blowflies (L. sericata, Calliphora vicina) fed on toad carcasses. Analysis of the Na+, K+-ATPase gene revealed no amino acid substitution at positions known to mediate resistance to bufadienolides in other systems. Alternative mechanisms of resistance such as efficient excretion of the compounds may enable the flies to use this poisonous food source.

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. CRISPR–Cas9 engineering of the Drosophila Atpα gene (encoding the α-subunit of the sodium pump) is used to study the ability of mutations that evolved independently in several insect orders to confer resistance to keystone plant toxins.

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.

Infections with maternally inherited Wolbachia bacteria may have dramatic influences on reproductive traits and speciation patterns of their hosts. We here show that in the beetle genus Altica, infection has influenced phylogenetic patterns of the host's mtDNA and different strains led to repeated selective sweeps. By comparing a COI/II-based phylogeny of the hosts with a phylogeny of the bacteria based on ftsZ, we show that cospeciation is rare and restricted to few recently diverged species. While in general each species apparently harbours a single Wolbachia strain, Altica lythri presents a strikingly different pattern: in the polyphyletic species, three highly divergent mtDNA haplotypes (2.1-4.6% p-distance) are coupled with three different Wolbachia strains (wLytA1, wLytA2 and wLytB). These haplotypes and Wolbachia strains are widely distributed and mostly found in sympatry. A phylogeny based on microsatellite data supports the monophyly of A. lythri. The discrepancy between mtDNA and nuclear phylogeny may best be explained by interspecific hybridization that led to introgression of mtDNA coupled with a different Wolbachia strain. Selective sweeps apparently drove the introgressed haplotypes to widespread distribution. As for effects of Wolbachia on reproduction, infection with wLytA1 appears to be correlated with a substantial sex ratio distortion, which was most prominent in A. lythri.

Despite the monarch butterfly (Danaus plexippus) being famous for its adaptations to the defensive traits of its milkweed host plants, little is known about the macroevolution of these traits. Unlike most other animal species, monarchs are largely insensitive to cardenolides, because their target site, the sodium pump (Na(+)/K(+) -ATPase), has evolved amino acid substitutions that reduce cardenolide binding (so-called target site insensitivity, TSI). Because many, but not all, species of milkweed butterflies (Danaini) are associated with cardenolide-containing host plants, we analyzed 16 species, representing all phylogenetic lineages of milkweed butterflies, for the occurrence of TSI by sequence analyses of the Na(+)/K(+) -ATPase gene and by enzymatic assays with extracted Na(+)/K(+) -ATPase. Here we report that sensitivity to cardenolides was reduced in a stepwise manner during the macroevolution of milkweed butterflies. Strikingly, not all Danaini typically consuming cardenolides showed TSI, but rather TSI was more strongly associated with sequestration of toxic cardenolides. Thus, the interplay between bottom-up selection by plant compounds and top-down selection by natural enemies can explain the evolutionary sequence of adaptations to these toxins.

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

The extent of convergent molecular evolution is largely unknown, yet is critical to understanding the genetics of adaptation. Target site insensitivity to cardenolides is a prime candidate for studying molecular convergence because herbivores in six orders of insects have specialized on these plant poisons, which gain their toxicity by blocking an essential transmembrane carrier, the sodium pump (Na,K-ATPase). We investigated gene sequences of the Na,K-ATPase α-subunit in 18 insects feeding on cardenolide-containing plants (spanning 15 genera and four orders) to screen for amino acid substitutions that might lower sensitivity to cardenolides. The replacement N122H that was previously shown to confer resistance in the monarch butterfly ( Danaus plexippus ) and Chrysochus leaf beetles was found in four additional species, Oncopeltus fasciatus and Lygaeus kalmii (Heteroptera, Lygaeidae), Labidomera clivicollis (Coleoptera, Chrysomelidae), and Liriomyza asclepiadis (Diptera, Agromyzidae). Thus, across 300 Myr of insect divergence, specialization on cardenolide-containing plants resulted in molecular convergence for an adaptation likely involved in coevolution. Our screen revealed a number of other substitutions connected to cardenolide binding in mammals. We confirmed that some of the particular substitutions provide resistance to cardenolides by introducing five distinct constructs of the Drosophila melanogaster gene into susceptible eucaryotic cells under an ouabain selection regime. These functional assays demonstrate that combined substitutions of Q 111 and N 122 are synergistic, with greater than twofold higher resistance than either substitution alone and >12-fold resistance over the wild type. Thus, even across deep phylogenetic branches, evolutionary degrees of freedom seem to be limited by physiological constraints, such that the same molecular substitutions confer adaptation.

We analyzed several species of the weevil family Mecininae (Coleoptera, Curculionidae) that all feed on iridoid glycoside (IG) containing plants of the Plantaginaceae to investigate whether the beetles sequester these deterrent substances from their host plants. Within the Mecininae two genera of the tribe Cionini were found to sequester aucubin and catalpol: Cionus Clairvilleand Schellenberg and Cleopus Dejean. Both analyzed genera of the Mecinini, Mecinus Germar and Rhinusa Stephens, do not sequester IGs although the compounds are present in their food plants. They thus represent the first case of specialists on IG plants that have not evolved adaptations to use the compounds. However, in contrast to the Cionini these genera have a hidden lifestyle, so that their need for defence might be lower. Both Cionus and Cleopus, sequester catalpol with a higher efficiency than aucubin. However, in contrast to Cionus species, Cleopus species only sequester catalpol. In species feeding on Scrophularia, the aucubin concentration is higher while in beetles on Verbascum catalpol is usually dominating. This pattern can also be detected in the only species living on both plants, Cionus hortulanus. The ability to sequester IGs must have a single origin at the base of the sister genera Cionus and Cleopus.

Specializing on host plants with toxic secondary compounds enforces specific adaptation in insect herbivores. In this review, we focus on two compound classes, iridoid glycosides and cardenolides, which can be found in the food plants of a large number of insect species that display various degrees of adaptation to them. These secondary compounds have very different modes of action: Iridoid glycosides are usually activated in the gut of the herbivores by β-glucosidases that may either stem from the food plant or be present in the gut as standard digestive enzymes. Upon cleaving, the unstable aglycone is released that unspecifically acts by crosslinking proteins and inhibiting enzymes. Cardenolides, on the other hand, are highly specific inhibitors of an essential ion carrier, the sodium pump. In insects exposed to both kinds of toxins, carriers either enabling the safe storage of the compounds away from the activating enzymes or excluding the toxins from sensitive tissues, play an important role that deserves further analysis. To avoid toxicity of iridoid glycosides, repression of activating enzymes emerges as a possible alternative strategy. Cardenolides, on the other hand, may lose their toxicity if their target site is modified and this strategy has evolved multiple times independently in cardenolide-adapted insects.

Sequestration of toxic plant compounds by herbivorous insects as a defence against predators has been observed in many tritrophic systems. In this study, we focused on the potential benefit of sequestered cardenolides, potent toxins that block the Na+K+-ATPase against an invertebrate predator. To evaluate the effect of cardenolides we used three different types of cardenolide-containing prey, which were offered to the orb-weaving spider Nephila senegalensis. As a natural system spiders were presented with the cardenolide-sequestering moth Empyreuma pugione. To reduce the complexity of the predator–prey interaction spiders were furthermore fed with ouabain-laced non-toxic flies. Lastly, we developed a new reductionistic assay involving gelatine capsules as insect dummies. Our results show that ouabain elicits significant rejection behaviour and suggest cardenolides to be the basis for rejection of E. pugione by N. senegalensis.

Iridoid glycosides are secondary plant compounds that have deterrent, growth reducing or even toxic effects on non-adapted herbivorous insects. To investigate the effects of iridoid glycoside containing plants on the digestive metabolism of a generalist herbivore, larvae of Spilosoma virginica (Lepidoptera: Arctiidae) were reared on three plant species that differ in their secondary plant chemistry: Taraxacum officinale (no iridoid glycosides), Plantago major (low iridoid glycoside content), and P. lanceolata (high iridoid glycoside content). Midguts of fifth instar larvae were assayed for the activity and kinetic properties of β-glucosidase using different substrates. Compared to the larvae on T. officinale, the β-glucosidase activity of larvae feeding on P. lanceolata was significantly lower measured with 4-nitrophenyl-β-d-glucopyranoside. Using the iridoid glycoside aucubin as a substrate, we did not find differences in the β-glucosidase activity of the larvae reared on the three plants. Heat inactivation experiments revealed the existence of a heat-labile and a more heat-stable β-glucosidase with similar Michaelis constants for 4-nitrophenyl-β-d-glucopyranoside. We discuss possible mechanisms leading to the observed decrease of β-glucosidase activity for larvae reared on P. lanceolata and its relevance for generalist herbivores in adapting to iridoid glycoside containing plant species and their use as potential host plants.

The caterpillars of the oleander hawk moth, Daphnis nerii (Linnaeus, 1758) (Lepidoptera: Sphingidae) feed primarily on oleander (Nerium oleander). This plant is rich in cardenolides, which specifically inhibit the Na+K+-ATPase. Since some insects feeding on cardenolide plants possess cardenolide-resistant Na+K+-ATPases, we tested whether D. nerii also possesses this strategy for circumventing cardenolide toxicity. To do so, we established a physiological assay, which allowed direct measurement of Na+K+-ATPase cardenolide sensitivity. Using Schistocerca gregaria, as a cardenolide-sensitive reference species, we showed that D. nerii Na+K+-ATPase was extremely sensitive to the cardenolide ouabain. Surprisingly, its sensitivity is even higher than that of the cardenolide-sensitive generalist, S. gregaria. The presence or absence of cardenolides in the diet of D. nerii did not influence the enzyme’s cardenolide sensitivity, indicating that target-site insensitivity is not inducible in this species. However, despite the sensitivity of their Na+K+-ATPase, caterpillars of D. nerii quickly recovered from an injection of an excessive amount of ouabain into their haemocoel. We conclude that D. nerii possesses adaptations, which enable it to feed on a cardenolide-rich diet other than that previously described in cardenolide specialized insects, and discuss other potential resistance mechanisms.

Whenever potentially noxious plant compounds are taken up and recycled by herbivorous insects, a protective function of these sequestered compounds is assumed. The flea beetle Longitarsus melanocephalus sequesters iridoid glycosides from its host plant up to a concentration of 2% DW, yet so far it remained unknown whether the insects gain protection from natural enemies by sequestering plant compounds at these comparatively low concentrations. Here we tested whether iridoid glycosides might deter or inhibit various soil and litter-dwelling potential enemies and pathogens. In choice experiments presenting L. melanocephalus pupae together with Tribolium castaneum pupae, the predator Lithobius forficatus was deterred by the iridoid glycoside containing pupae, while Forficula auricularia as well as the nematode Heterorhabditis bacteriophora were not deterred. L. forficatus also avoided artificial baits doted with 2% iridoid glycosides while F. auricularia showed no aversion to iridoid glycosides at these concentrations and H. bacteriophora did not suffer any toxic effect. Of the pathogens tested, the entomopathogenic fungi Beauveria bassiana and Metarhizium anisopliae were not inhibited in their growth by iridoid glycosides ranging up to 2%. However, an inhibitory effect could be observed against the entomopathogenic bacterium Bacillus thuringiensis , even at the relatively small concentrations that are common in L. melanocephalus . The antibacterial effect might thus be another important selective value of iridoid glycoside sequestration in this species.

Selection pressure by natural enemies on phytophagous insect larvae is intense and has frequently triggered the evolution of chemical defence as an effective counterstrategy. In the chrysomelid subfamily Galerucinae, glandular structures and defensive fluids have been described for the tribe Sermylini Wilcox, 1965. Previous morphological and ultrastructural studies raised doubts that these defensive devices in Sermylini can be traced back to a common origin. The taxonomy of the Galerucinae cannot clear these doubts because the phylogeny of this taxon is a matter of current debate. We therefore investigated the phylogeny of the Galerucinae based on approximately 1740 bp of the mitochondrial 12S and 16S rRNA and the nuclear elongation factor 1 alpha genes. Our data support the hypothesized close relationship between the subfamilies Galerucinae and Alticinae, yet, by contrast to other recent analyses, the two groups are mostly resolved as monophyletic sister groups or, in some analyses, with the Galerucinae nested paraphyletically within the Alticinae. Within the subfamily Galerucinae, only the tribe Galerucini formed a monophyletic taxon, except for one species, Cerochroa brachialis Stal, 1858. In none of our analyses were the Sermylini recovered as a monophyletic tribe. However, our data support monophyly of each of the three groups within the Sermylini that have morphologically distinguishable larval defensive openings. We conclude that the defensive structures in larvae of Sermylini have no common origin, but evolved independently. Our data suggest that the tremendous selection pressure by natural enemies led to the recurrent evolution of similar chemical defensive devices in Sermylini larvae. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 93, 165–175.