Divergence of defensive cucurbitacins in independent Cucurbita pepo domestication events leads to differences in specialist herbivore preference

Crop domestication and improvement often concurrently affect plant resistance to pests and production of secondary metabolites, creating challenges for isolating the ecological implications of selection for specific metabolites. Cucurbitacins are bitter triterpenoids with extreme phenotypic differences between Cucurbitaceae lineages, yet we lack integrated models of herbivore preference, cucurbitacin accumulation, and underlying genetic mechanisms. In Cucurbita pepo, we dissected the effect of cotyledon cucurbitacins on preference of a specialist insect pest (Acalymma vittatum) for multiple tissues, assessed genetic loci underlying cucurbitacin accumulation in diverse germplasm and a biparental F₂ population (from a cross between two independent domesticates), and characterized quantitative associations between gene expression and metabolites during seedling development. Acalymma vittatum affinity for cotyledons is mediated by cucurbitacins, but other traits contribute to whole-plant resistance. Cotyledon cucurbitacin accumulation was associated with population structure, and our genetic mapping identified a single locus, Bi-4, containing genes relevant to transport and regulation - not biosynthesis - that diverged between lineages. These candidate genes were expressed during seedling development, most prominently a putative secondary metabolite transporter. Taken together, these findings support the testable hypothesis that breeding for plant resistance to insects involves targeting genes for regulation and transport of defensive metabolites, in addition to core biosynthesis genes.

The role of toxic nectar secondary compounds in driving differential bumble bee preferences for milkweed flowers

While morphological differences such as tongue length are often featured as drivers of pollinator floral preferences, differences in chemical detection and tolerance to secondary compounds may also play a role. We sought to better understand the role of secondary compounds in floral preference by examining visitation of milkweed flowers, which can contain toxic cardenolides in their nectar, by bumble bees (Bombus spp.), some of their most abundant and important pollinators. We examine bumble bee species visitation of common milkweed (Asclepias syriaca) compared to other flowers in the field and test whether observed preferences may be influenced by avoidance and tolerance of cardenolides, as measured by the cardenolide ouabain, in the lab. We reveal that common milkweed is visited predominantly by one bumble bee species, Bombus griseocollis, in a ratio much higher than the abundance of this species in the community. We confirmed the presence and toxicity of cardenolides in A. syriaca nectar. Lab experiments revealed that B. griseocollis, compared to the common bumble bees B. impatiens and B. bimaculatus, exhibit greater avoidance of cardenolides, but only at levels that start to induce illness, whereas the other species exhibit either no or reduced avoidance of cardenolides, resulting in illness and mortality in these bees. Toxicity experiments reveal that B. griseocollis also has a substantially higher tolerance for cardenolides than B. impatiens. Together, these results support a potential evolutionary association between B. griseocollis and milkweed that may involve increased ability to both detect and tolerate milkweed cardenolides.

Butterflies and the people who love them Wendy Williams Simon and Schuster, 2020. 240 pp

A meandering investigation hints at how much is left to learn about these charismatic insects

Agrobacterium tumefaciens-Mediated Transformation of Three Milkweed Species (Asclepias hallii, A. syriaca, and A. tuberosa: Apocynaceae)

Milkweeds have ecological significance for insect herbivores that rely on them as hosts for either part of or the entirety of their life cycles. Interesting interactions, some of which are not completely understood, have evolved over time. To develop these species as models to elucidate the interplay with insect herbivores, we established Agrobacterium tumefaciens-mediated transformation approaches for Asclepias hallii (Hall's milkweed), A. syriaca (common milkweed), and A. tuberosa (butterflyweed). The method is based on infection of stem internodal explants, which were more amenable to transformation than leaf explants. We found that addition of freshly prepared dithiothreitol was critical to prevent browning of stem explants. Depending on the species, the time from infection to the regeneration of transgenic lines ranges from 2 to 4 months. Transformation efficiency for A. hallii was 9%, whereas efficiencies for A. syriaca and A. tuberosa were 6% and 13%, respectively. © 2020 by John Wiley & Sons, Inc. Basic Protocol 1: Agrobacterium tumefaciens-mediated transformation of Asclepias internodal stem explants Basic Protocol 2: Preparation of Agrobacterium glycerol stocks containing gene constructs.

Evolution of phenotypic plasticity: Genetic differentiation and additive genetic variation for induced plant defence in wild arugula Eruca sativa

Phenotypic plasticity is the primary mechanism of organismal resilience to abiotic and biotic stress, and genetic differentiation in plasticity can evolve if stresses differ among populations. Inducible defence is a common form of adaptive phenotypic plasticity, and long-standing theory predicts that its evolution is shaped by costs of the defensive traits, costs of plasticity and a trade-off in allocation to constitutive versus induced traits. We used a common garden to study the evolution of defence in two native populations of wild arugula Eruca sativa (Brassicaceae) from contrasting desert and Mediterranean habitats that differ in attack by caterpillars and aphids. We report genetic differentiation and additive genetic variance for phenology, growth and three defensive traits (toxic glucosinolates, anti-nutritive protease inhibitors and physical trichome barriers) as well their inducibility in response to the plant hormone jasmonic acid. The two populations were strongly differentiated for plasticity in nearly all traits. There was little evidence for costs of defence or plasticity, but constitutive and induced traits showed a consistent additive genetic trade-off within each population for the three defensive traits. We conclude that these populations have evolutionarily diverged in inducible defence and retain ample potential for the future evolution of phenotypic plasticity in defence.

Plant Defense by Latex: Ecological Genetics of Inducibility in the Milkweeds and a General Review of Mechanisms, Evolution, and Implications for Agriculture

Latex occurs in 10% of plant families, has evolved independently many times, and is the most effective defense of milkweeds against its chewing herbivores. Here we report on new experiments on the heritability and inducibility of latex in several milkweed species. In addition, we review what is known about the genetic and environmental determinants of latex exudation, hormonal regulation, evolution within and among species, and the role and frequency of latex in agricultural crops. We first evaluated genotype-by-environment interactions using ~20 full-sibling genetic families in each of seven Asclepias species treated as controls or attacked by monarch butterfly caterpillars. All species showed substantial genetic variation for latex exudation and six of seven species responded to monarch herbivory (two species increased latex, two species decreased, and two showed variation among genetic families). Exogenous application of jasmonic acid (JA) to three species induced a consistent increase in latex (including species which showed a decline following caterpillar herbivory). We next evaluated three hypotheses for what determines genetic variation for induced latex in A. syriaca: 1) a trade-off with constitutive investment, 2) differential endogenous JA induction, or 3) variation in responsiveness to JA. We only found support for the second hypothesis: genetic families with a stronger JA-burst showed the greatest latex exudation following herbivory. We conclude that most species exhibit a genetic and inducible basis for latex, although genetic variation in inducibility is not pervasive. Finally, we summarized studies across 22 species of Asclepias and found that neither a species' latitude nor its phylogenetic position predicted latex inducibility. Nonetheless, a negative association between constitutive and induced latex across species indicates a macroevolutionary trade-off in allocation to this defense. Our review indicates that jasmonic acid is a key regulator of latex exudation, laticifer morphology, and defensive metabolites within latex. Biotic and abiotic factors strongly modulate latex expression. A survey of latex in food crops revealed that latex and analogous exudates (gums, resins, mucilage) are more common than expected based on their distribution across all plants. In conclusion, despite its widespread occurrence, the literature on latex is currently dominated by rubber trees and milkweeds, and we look forward to the broadening of ecological, agricultural, and mechanistic research into other systems.

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 fitness¹. Convergent evolution can be used as a guide to identify candidate mutations that underlie adaptive traits^2-4, and new genome editing technology is facilitating functional validation of these mutations in whole organisms^1,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 toxins^6-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 glycosides¹². Here we describe mutational paths involving three repeatedly changing amino acid sites (111, 119 and 122) in ATPα that are associated with cardiac glycoside specialization^13,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 lineage^11,15. We show in vivo, in vitro and in silico that the path conferred resistance and target-site insensitivity to cardiac glycosides¹⁶, 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 predators^17-19. The order in which the substitutions evolved was explained by amelioration of antagonistic pleiotropy through epistasis^13,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 communities^{2,6-11,15,17-19}.

Trade-offs constrain the evolution of an inducible defense within but not between plant species

Inducible defense is a common form of phenotypic plasticity, and inducibility (change in defense after herbivore attack) has long been predicted to trade off with constitutive (or baseline) defense to manage resource allocation. Although such trade-offs likely constrain evolution within species, the extent to which they influence divergence among species is unresolved. We studied cardenolide toxins among genetic families in eight North American Asclepias species, spanning the full range of defense in the genus. Using common environment experiments and chemical assays, we report a consistent trade-off (negative genetic correlation) between concentrations of constitutive cardenolides and their inducibility within each species. However, no trade-off was found in a phylogenetic analysis across species. To investigate factors driving differences in defense allocation among species we used latitude as a proxy for growing season and herbivore pressure and found that divergence into lower latitudes resulted in evolution of higher cardenolides overall. Next we used an enzymatic assay of the cellular target of cardenolides (sodium-potassium ATPase) and confirm that higher cardenolides resulted in stronger toxicity to a sensitive species, but not to specialized monarch butterflies. Thus, plant speciation into biogeographic regions with alternative resources or pest pressure resulted in the macroevolution of cardenolide defense, especially against unspecialized herbivores. Nonetheless, trade-offs persist in the extent to which this defense is allocated constitutively or is inducible among genotypes within each species.

Ontogenetic strategies in insect herbivores and their impact on tri-trophic interactions

Insect herbivores express tremendous ontogenetic variation in traits related to growth and maturation, but also as an evolutionary consequence of ecological interactions with plants and predators. These selective pressures can either reinforce or restrict expression of particular ontogenetic strategies, allowing herbivores to simultaneously cope with plant resistance and risk of predation through ontogenetic change. For example, whereas an increase in defense-sabotaging behavior, aposematism and sequestration along herbivore ontogeny seems to be reinforced by both bottom-up and top-down forces, some ontogenetic trends in anti-predator behavior can be limited by plant resistance. Communication among plants, herbivores and their natural enemies is also influenced by insect ontogenies. The study of ontogenetic strategies of herbivores requires the assessment of the genetic variation, heritability and adaptive value across herbivore development, considering the variation in plant quality and predation risk.

Cardenolide Intake, Sequestration, and Excretion by the Monarch Butterfly along Gradients of Plant Toxicity and Larval Ontogeny

Monarch butterflies, Danaus plexippus, migrate long distances over which they encounter host plants that vary broadly in toxic cardenolides. Remarkably little is understood about the mechanisms of sequestration in Lepidoptera that lay eggs on host plants ranging in such toxins. Using closely-related milkweed host plants that differ more than ten-fold in cardenolide concentrations, we mechanistically address the intake, sequestration, and excretion of cardenolides by monarchs. We show that on high cardenolide plant species, adult butterflies saturate in cardenolides, resulting in lower concentrations than in leaves, while on low cardenolide plants, butterflies concentrate toxins. Butterflies appear to focus their sequestration on particular compounds, as the diversity of cardenolides is highest in plant leaves, lower in frass, and least in adult butterflies. Among the variety of cardenolides produced by the plant, sequestered compounds may be less toxic to the butterflies themselves, as they are more polar on average than those in leaves. In accordance with this, results from an in vitro assay based on inhibition of Na⁺/K⁺ ATPase (the physiological target of cardenolides) showed that on two milkweed species, including the high cardenolide A. perennis, extracts from butterflies have lower inhibitory effects than leaves when standardized by cardenolide concentration, indicating selective sequestration of less toxic compounds from these host plants. To understand how ontogeny shapes sequestration, we examined cardenolide concentrations in caterpillar body tissues and hemolymph over the course of development. Caterpillars sequestered higher concentrations of cardenolides as early instars than as late instars, but within the fifth instar, concentration increased with body mass. Although it appears that large amounts of sequestration occurs in early instars, a host switching experiment revealed that caterpillars can compensate for feeding on low cardenolide host plants with substantial sequestration in the fifth instar. We highlight commonalities and striking differences in the mechanisms of sequestration depending on host plant chemistry and developmental stage, which have important implications for monarch defense.

Population Variation, Environmental Gradients, and the Evolutionary Ecology of Plant Defense against Herbivory

A central tenet of plant defense theory is that adaptation to the abiotic environment sets the template for defense strategies, imposing a trade-off between plant growth and defense. Yet this trade-off, commonly found among species occupying divergent resource environments, may not occur across populations of single species. We hypothesized that more favorable climates and higher levels of herbivory would lead to increases in growth and defense across plant populations. We evaluated whether plant growth and defense traits covaried across 18 populations of showy milkweed (Asclepias speciosa) inhabiting an east-west climate gradient spanning 25° of longitude. A suite of traits impacting defense (e.g., latex, cardenolides), growth (e.g., size), or both (e.g., specific leaf area [SLA], trichomes) were measured in natural populations and in a common garden, allowing us to evaluate plastic and genetically based variation in these traits. In natural populations, herbivore pressure increased toward warmer sites with longer growing seasons. Growth and defense traits showed strong clinal patterns and were positively correlated. In a common garden, clines with climatic origin were recapitulated only for defense traits. Correlations between growth and defense traits were also weaker and more negative in the common garden than in the natural populations. Thus, our data suggest that climatically favorable sites likely facilitate the evolution of greater defense at minimal costs to growth, likely because of increased resource acquisition.

Toxicity of Milkweed Leaves and Latex: Chromatographic Quantification Versus Biological Activity of Cardenolides in 16 Asclepias Species

Cardenolides are classically studied steroidal defenses in chemical ecology and plant-herbivore coevolution. Although milkweed plants (Asclepias spp.) produce up to 200 structurally different cardenolides, all compounds seemingly share the same well-characterized mode of action, inhibition of the ubiquitous Na⁺/K⁺ ATPase in animal cells. Over their evolutionary radiation, milkweeds show a quantitative decline of cardenolide production and diversity. This reduction is contrary to coevolutionary predictions and could represent a cost-saving strategy, i.e. production of fewer but more toxic cardenolides. Here we test this hypothesis by tandem cardenolide quantification using HPLC (UV absorption of the unsaturated lactone) and a pharmacological assay (in vitro inhibition of a sensitive Na⁺/K⁺ ATPase) in a comparative study of 16 species of Asclepias. We contrast cardenolide concentrations in leaf tissue to the subset of cardenolides present in exuding latex. Results from the two quantification methods were strongly correlated, but the enzymatic assay revealed that milkweed cardenolide mixtures often cause stronger inhibition than equal amounts of a non-milkweed reference cardenolide, ouabain. Cardenolide concentrations in latex and leaves were positively correlated across species, yet latex caused 27% stronger enzyme inhibition than equimolar amounts of leaf cardenolides. Using a novel multiple regression approach, we found three highly potent cardenolides (identified as calactin, calotropin, and voruscharin) to be primarily responsible for the increased pharmacological activity of milkweed cardenolide mixtures. However, contrary to an expected trade-off between concentration and toxicity, later-diverging milkweeds had the lowest amounts of these potent cardenolides, perhaps indicating an evolutionary response to milkweed's diverse community of specialist cardenolide-sequestering insect herbivores.

Population Variation, Environmental Gradients, and the Evolutionary Ecology of Plant Defense against Herbivory

A central tenet of plant defense theory is that adaptation to the abiotic environment sets the template for defense strategies, imposing a trade-off between plant growth and defense. Yet this trade-off, commonly found among species occupying divergent resource environments, may not occur across populations of single species. We hypothesized that more favorable climates and higher levels of herbivory would lead to increases in growth and defense across plant populations. We evaluated whether plant growth and defense traits covaried across 18 populations of showy milkweed (Asclepias speciosa) inhabiting an east-west climate gradient spanning 25° of longitude. A suite of traits impacting defense (e.g., latex, cardenolides), growth (e.g., size), or both (e.g., specific leaf area [SLA], trichomes) were measured in natural populations and in a common garden, allowing us to evaluate plastic and genetically based variation in these traits. In natural populations, herbivore pressure increased toward warmer sites with longer growing seasons. Growth and defense traits showed strong clinal patterns and were positively correlated. In a common garden, clines with climatic origin were recapitulated only for defense traits. Correlations between growth and defense traits were also weaker and more negative in the common garden than in the natural populations. Thus, our data suggest that climatically favorable sites likely facilitate the evolution of greater defense at minimal costs to growth, likely because of increased resource acquisition.

Plant chemical defense indirectly mediates aphid performance via interactions with tending ants

The benefits of mutualistic interactions are often highly context dependent. We studied the interaction between the milkweed aphid Aphis asclepiadis and a tending ant, Formica podzolica. Although this interaction is generally considered beneficial, variation in plant genotype may alter it from mutualistic to antagonistic. Here we link the shift in strength and relative benefit of the ant-aphid interaction to plant genotypic variation in the production of cardenolides, a class of toxic defensive chemicals. In a field experiment with highly variable genotypes of the common milkweed (Asclepias syriaca), we show that plant cardenolides, especially polar forms, are ingested by aphids and excreted in honeydew proportionally to plant concentrations without directly affecting aphid performance. Ants consume honeydew, and aphids that excreted high amounts of cardenolides received fewer ant visits, which in turn reduced aphid survival. On at least some plant genotypes, aphid numbers per plant were reduced in the presence of ants to levels lower than in corresponding ant-exclusion treatments, suggesting antagonistic ant behavior. Although cardenolides appear ineffective as direct plant defenses against aphids, the multi-trophic context reveals an ant-mediated negative indirect effect on aphid performance and population dynamics.

Relative Selectivity of Plant Cardenolides for Na+/K+-ATPases From the Monarch Butterfly and Non-resistant Insects

A major prediction of coevolutionary theory is that plants may target particular herbivores with secondary compounds that are selectively defensive. The highly specialized monarch butterfly (Danaus plexippus) copes well with cardiac glycosides (inhibitors of animal Na+/K+-ATPases) from its milkweed host plants, but selective inhibition of its Na+/K+-ATPase by different compounds has not been previously tested. We applied 17 cardiac glycosides to the D. plexippus-Na+/K+-ATPase and to the more susceptible Na+/K+-ATPases of two non-adapted insects (Euploea core and Schistocerca gregaria). Structural features (e.g., sugar residues) predicted in vitro inhibitory activity and comparison of insect Na+/K+-ATPases revealed that the monarch has evolved a highly resistant enzyme overall. Nonetheless, we found evidence for relative selectivity of individual cardiac glycosides reaching from 4- to 94-fold differences of inhibition between non-adapted Na+/K+-ATPase and D. plexippus-Na+/K+-ATPase. This toxin receptor specificity suggests a mechanism how plants could target herbivores selectively and thus provides a strong basis for pairwise coevolutionary interactions between plants and herbivorous insects.

Toxicity of the spiny thick-foot Pachypodium

PREMISE OF THE STUDY

Pachypodium (Apocynaceae) is a genus of iconic stem-succulent and poisonous plants endemic to Madagascar and southern Africa. We tested hypotheses about the mode of action and macroevolution of toxicity in this group. We further hypothesized that while monarch butterflies are highly resistant to cardenolide toxins (a type of cardiac glycoside) from American Asclepias, they may be negatively affected by Pachypodium defenses, which evolved independently.

METHODS

We grew 16 of 21 known Pachypodium spp. and quantified putative cardenolides by HPLC and also by inhibition of animal Na⁺ /K⁺ -ATPase (the physiological target of cardiac glycosides) using an in vitro assay. Pachypodium extracts were tested against monarch caterpillars in a feeding bioassay. We also tested four Asclepias spp. and five Pachypodium spp. extracts, contrasting inhibition of the cardenolide-sensitive porcine Na⁺ /K⁺ -ATPase to the monarch's resistant form.

KEY RESULTS

We found evidence for low cardenolides by HPLC, but substantial toxicity when extracts were assayed on Na⁺ /K⁺ -ATPases. Toxicity showed phylogenetic signal, and taller species showed greater toxicity (this was marginal after phylogenetic correction). Application of Pachypodium extracts to milkweed leaves reduced monarch growth, and this was predicted by inhibition of the sensitive Na⁺ /K⁺ -ATPase in phylogenetic analyses. Asclepias extracts were 100-fold less potent against the monarch compared to the porcine Na⁺ /K⁺ -ATPase, but this difference was absent for Pachypodium extracts.

CONCLUSIONS

Pachypodium contains potent toxicity capable of inhibiting sensitive and cardenolide-adapted Na⁺ /K⁺ -ATPases. Given the monarch's sensitivity to Pachypodium, we suggest that these plants contain novel cardiac glycosides or other compounds that facilitate toxicity by binding to Na⁺ /K⁺ -ATPases.

Fitness consequences of occasional outcrossing in a functionally asexual plant (Oenothera biennis)

Many clonal organisms occasionally outcross, but the long-term consequences of such infrequent events are often unknown. During five years, representing three to five plant generations, we followed 16 experimental field populations of the forb, Oenothera biennis, originally planted with the same 18 original genotypes. Oenothera biennis usually self fertilizes, which, due to its genetic system (permanent translocation heterozygosity), results in seeds that are clones of the maternal plant. However, rare outcrossing produces genetically novel offspring (but without recombination or increased heterozygosity). We sought to understand whether novel genotypes produced through natural outcrossing had greater fecundity or different multigenerational dynamics compared to our original genotypes. We further assessed whether any differences in fitness or abundances through time between original and novel genotypes were exaggerated in the presence vs. absence of insect herbivores. Over the course of the experiment, we genotyped >12,500 plants using microsatellite DNA markers to identify and track the frequency of specific genotypes and estimated fecundity on a subset (>3,000) of plants. The effective outcrossing rate was 7.3% in the first year and ultimately 50% of the plants were of outcrossed origin by the final year of the experiment. Lifetime fruit production per plant was on average 32% higher across all novel genotypes produced via outcrossing compared to the original genotypes, and this fecundity advantage was significantly enhanced in populations lacking herbivores. Among 43 novel genotypes that were abundant enough to phenotype with replication, plants produced nearly 30% more fruits than the average of their specific two parental genotypes, and marginally more fruits (8%) than their most fecund parent. Mean per capita fecundity of novel genotypes predicted their relative frequencies at the end of the experiment. Novel genotypes increased more dramatically in herbivore-present compared to suppressed populations (45% vs. 27% of all plants), countering the increased competition from dandelions (Taraxacum officinale) that resulted from herbivore suppression. Increased interspecific competition likely also lead to the lower realized fitness of novel vs. original genotypes in herbivore-suppressed populations. These results demonstrate that rare outcrossing and the generation of novel genotypes can create high-fecundity progeny, with the biotic environment influencing the dynamical outcome of such advantages.