Independent evolution of ancestral and novel defenses in a genus of toxic plants (Erysimum, Brassicaceae)

Convergently evolved metabolites are new to me but not to my attackers

This article is a Commentary on Younkin et al. (2024), 242: 2719–2733.

Cardiac glycosides protect wormseed wallflower (Erysimum cheiranthoides) against some, but not all, glucosinolate-adapted herbivores

The chemical arms race between plants and insects is foundational to the generation and maintenance of biological diversity. We asked how the evolution of a novel defensive compound in an already well-defended plant lineage impacts interactions with diverse herbivores. Erysimum cheiranthoides (Brassicaceae), which produces both ancestral glucosinolates and novel cardiac glycosides, served as a model. We analyzed gene expression to identify cardiac glycoside biosynthetic enzymes in E. cheiranthoides and characterized these enzymes via heterologous expression and CRISPR/Cas9 knockout. Using E. cheiranthoides cardiac glycoside-deficient lines, we conducted insect experiments in both the laboratory and field. EcCYP87A126 initiates cardiac glycoside biosynthesis via sterol side-chain cleavage, and EcCYP716A418 has a role in cardiac glycoside hydroxylation. In EcCYP87A126 knockout lines, cardiac glycoside production was eliminated. Laboratory experiments with these lines revealed that cardiac glycosides were highly effective defenses against two species of glucosinolate-tolerant specialist herbivores, but did not protect against all crucifer-feeding specialist herbivores in the field. Cardiac glycosides had lesser to no effect on two broad generalist herbivores. These results begin elucidation of the E. cheiranthoides cardiac glycoside biosynthetic pathway and demonstrate in vivo that cardiac glycoside production allows Erysimum to escape from some, but not all, specialist herbivores.

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

The Madagascar palm genome provides new insights on the evolution of Apocynaceae specialized metabolism

Specialized metabolites possess diverse interesting biological activities and some cardenolides- and monoterpene indole alkaloids- (MIAs) derived pharmaceuticals are currently used to treat human diseases such as cancers or hypertension. While these two families of biocompounds are produced by specific subfamilies of Apocynaceae, one member of this medicinal plant family, the succulent tree Pachypodium lamerei Drake (also known as Madagascar palm), does not produce such specialized metabolites. To explore the evolutionary paths that have led to the emergence and loss of cardenolide and MIA biosynthesis in Apocynaceae, we sequenced and assembled the P. lamerei genome by combining Oxford Nanopore Technologies long-reads and Illumina short-reads. Phylogenomics revealed that, among the Apocynaceae whose genomes have been sequenced, the Madagascar palm is so far the species closest to the common ancestor between MIA producers/non-MIA producers. Transposable elements, constituting 72.48% of the genome, emerge as potential key players in shaping genomic architecture and influencing specialized metabolic pathways. The absence of crucial MIA biosynthetic genes such as strictosidine synthase in P. lamerei and non-Rauvolfioideae species hints at a transposon-mediated mechanism behind gene loss. Phylogenetic analysis not only showcases the evolutionary divergence of specialized metabolite biosynthesis within Apocynaceae but also underscores the role of transposable elements in this intricate process. Moreover, we shed light on the low conservation of enzymes involved in the final stages of MIA biosynthesis in the distinct MIA-producing plant families, inferring independent gains of these specialized enzymes along the evolution of these medicinal plant clades. Overall, this study marks a leap forward in understanding the genomic dynamics underpinning the evolution of specialized metabolites biosynthesis in the Apocynaceae family, with transposons emerging as potential architects of genomics restructuring and gene loss.

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

Aphid Resistance Segregates Independently of Cardenolide and Glucosinolate Content in an Erysimum cheiranthoides (Wormseed Wallflower) F2 Population

Plants in the genus Erysimum produce both glucosinolates and cardenolides as a defense mechanism against herbivory. Two natural isolates of Erysimum cheiranthoides (wormseed wallflower) differed in their glucosinolate content, cardenolide content, and their resistance to Myzus persicae (green peach aphid), a broad generalist herbivore. Both classes of defensive metabolites were produced constitutively and were not further induced by aphid feeding. To investigate the relative importance of glucosinolates and cardenolides in E. cheiranthoides defense, we generated an improved genome assembly, genetic map, and segregating F2 population. The genotypic and phenotypic analysis of the F2 plants identified quantitative trait loci, which affected glucosinolates and cardenolides, but not the aphid resistance. The abundance of most glucosinolates and cardenolides was positively correlated in the F2 population, indicating that similar processes regulate their biosynthesis and accumulation. Aphid reproduction was positively correlated with glucosinolate content. Although the overall cardenolide content had little effect on aphid growth and survival, there was a negative correlation between aphid reproduction and helveticoside abundance. However, this variation in defensive metabolites could not explain the differences in aphid growth on the two parental lines, suggesting that processes other than the abundance of glucosinolates and cardenolides have a predominant effect on aphid resistance in E. cheiranthoides.

Aphid resistance segregates independently of cardiac glycoside and glucosinolate content in an Erysimum cheiranthoides (wormseed wallflower) F2 population

Plants in the genus Erysimum produce both glucosinolates and cardiac glycosides as defense against herbivory. Two natural isolates of Erysimum cheiranthoides (wormseed wallflower) differed in their glucosinolate content, cardiac glycoside content, and resistance to Myzus persicae (green peach aphid), a broad generalist herbivore. Both classes of defensive metabolites were produced constitutively and were not induced further by aphid feeding. To investigate the relative importance of glucosinolates and cardiac glycosides in E. cheiranthoides defense, we generated an improved genome assembly, genetic map, and segregating F2 population. Genotypic and phenotypic analysis of the F2 plants identified quantitative trait loci affecting glucosinolates and cardiac glycosides, but not aphid resistance. The abundance of most glucosinolates and cardiac glycosides was positively correlated in the F2 population, indicating that similar processes regulate their biosynthesis and accumulation. Aphid reproduction was positively correlated with glucosinolate content. Although overall cardiac glycoside content had little effect on aphid growth and survival, there was a negative correlation between aphid reproduction and helveticoside abundance. However, this variation in defensive metabolites could not explain the differences in aphid growth on the two parental lines, suggesting that processes other than the abundance of glucosinolates and cardiac glycosides have a predominant effect on aphid resistance in E. cheiranthoides.

Promiscuous CYP87A enzyme activity initiates cardenolide biosynthesis in plants

Cardenolides are specialized, steroidal metabolites produced in a wide array of plant families1,2. Cardenolides play protective roles in plants, but these molecules, including digoxin from foxglove (Digitalis spp.), are better known for treatment of congenital heart failure, atrial arrhythmia, various cancers and other chronic diseases3-9. However, it is still unknown how plants synthesize 'high-value', complex cardenolide structures from, presumably, a sterol precursor. Here we identify two cytochrome P450, family 87, subfamily A (CYP87A) enzymes that act on both cholesterol and phytosterols (campesterol and β-sitosterol) to form pregnenolone, the first committed step in cardenolide biosynthesis in the two phylogenetically distant plants Digitalis purpurea and Calotropis procera. Arabidopsis plants overexpressing these CYP87A enzymes ectopically accumulated pregnenolone, whereas silencing of CYP87A in D. purpurea leaves by RNA interference resulted in substantial reduction of pregnenolone and cardenolides. Our work uncovers the key entry point to the cardenolide pathway, and expands the toolbox for sustainable production of high-value plant steroids via synthetic biology.

A pipeline for the rapid collection of color data from photographs

Premise

There are relatively few studies of flower color at landscape scales that can address the relative importance of competing mechanisms (e.g., biotic: pollinators; abiotic: ultraviolet radiation, drought stress) at landscape scales.

Methods

We developed an R shiny pipeline to sample color from images that were automatically downloaded using query results from a search using iNaturalist or the Global Biodiversity Information Facility (GBIF). The pipeline was used to sample ca. 4800 North American wallflower (Erysimum, Brassicaceae) images from iNaturalist. We tested whether flower color was distributed non-randomly across the landscape and whether spatial patterns were correlated with climate. We also used images including ColorCheckers to compare analyses of raw images to color-calibrated images.

Results

Flower color was strongly non-randomly distributed spatially, but did not correlate strongly with climate, with most of the variation explained instead by spatial autocorrelation. However, finer-scale patterns including local correlations between elevation and color were observed. Analyses using color-calibrated and raw images revealed similar results.

Discussion

This pipeline provides users the ability to rapidly capture color data from iNaturalist images and can be a useful tool in detecting spatial or temporal changes in color using citizen science data.

Influences of chemotype and parental genotype on metabolic fingerprints of tansy plants uncovered by predictive metabolomics

Intraspecific plant chemodiversity shapes plant-environment interactions. Within species, chemotypes can be defined according to variation in dominant specialised metabolites belonging to certain classes. Different ecological functions could be assigned to these distinct chemotypes. However, the roles of other metabolic variation and the parental origin (or genotype) of the chemotypes remain poorly explored. Here, we first compared the capacity of terpenoid profiles and metabolic fingerprints to distinguish five chemotypes of common tansy (Tanacetum vulgare) and depict metabolic differences. Metabolic fingerprints captured higher variation in metabolites while preserving the ability to define chemotypes. These differences might influence plant performance and interactions with the environment. Next, to characterise the influence of the maternal origin on chemodiversity, we performed variation partitioning and generalised linear modelling. Our findings revealed that maternal origin was a higher source of chemical variation than chemotype. Predictive metabolomics unveiled 184 markers predicting maternal origin with 89% accuracy. These markers included, among others, phenolics, whose functions in plant-environment interactions are well established. Hence, these findings place parental genotype at the forefront of intraspecific chemodiversity. We recommend considering this factor when comparing the ecology of various chemotypes. Additionally, the combined inclusion of inherited variation in main terpenoids and other metabolites in computational models may help connect chemodiversity and evolutionary principles.

Convergent and divergent evolution of plant chemical defenses

The majority of the several hundred thousand specialized metabolites produced by plants function in defense against insects and other herbivores. Despite this diversity, identical metabolites or structurally distinct metabolites hitting the same targets in herbivorous animals have evolved repeatedly. This convergent evolution may reflect the constraints of plant primary metabolism in providing metabolic precursors, as well as the limited number of readily accessible targets in animals. These restrictions may make it uncommon for plants to develop completely novel toxic and deterrent metabolites, despite the ongoing evolution of resistance mechanisms in insect herbivores. Defensive compounds that are unique to individual genera or species often have long biosynthetic pathways that may complicate the repeated evolution of these metabolites in different plant species.

Cycad phylogeny predicts host plant use of Eumaeus butterflies

Eumaeus butterflies are obligate herbivores of Zamia, the most diverse neotropical genus of cycads. Eumaeus-Zamia interactions have been characterized mainly for species distributed in North and Central America. However, larval host plant use by the southern Eumaeus clade remains largely unknown, precluding a comprehensive study of co-evolution between the genera. Here, we combine fieldwork with museum and literature surveys to expand herbivory records for Eumaeus from 21 to 38 Zamia species. We inferred a time-calibrated phylogeny of Eumaeus to test for distinct macroevolutionary scenarios of larval host plant conservatism and co-evolution. We found a remarkable coincidence between Eumaeus and Zamia diversification, with the butterfly stem group diverging at the same time as the most recent radiation of Zamia in the Miocene. Cophylogenetic reconciliation analyses show a strong cophylogenetic signal between cycads and their butterfly herbivores. Bipartite model-based approaches indicate that this is because closely related Zamia species are used by the same Eumaeus species, suggesting larval host plant resource tracking by the butterfly herbivores. Our results highlight a case of tight evolution between Eumaeus butterflies and cycads, pointing to the generality of correlated evolution and phylogenetic tracking in plant-herbivore interactions across seed plants.

Quantifying chemodiversity considering biochemical and structural properties of compounds with the R package chemodiv

Plants produce large numbers of phytochemical compounds affecting plant physiology and interactions with their biotic and abiotic environment. Recently, chemodiversity has attracted considerable attention as an ecologically and evolutionary meaningful way to characterize the phenotype of a mixture of phytochemical compounds. Currently used measures of phytochemical diversity, and related measures of phytochemical dissimilarity, generally do not take structural or biosynthetic properties of compounds into account. Such properties can be indicative of the compounds' function and inform about their biosynthetic (in)dependence, and should therefore be included in calculations of these measures. We introduce the R package chemodiv, which retrieves biochemical and structural properties of compounds from databases and provides functions for calculating and visualizing chemical diversity and dissimilarity for phytochemicals and other types of compounds. Our package enables calculations of diversity that takes the richness, relative abundance and - most importantly - structural and/or biosynthetic dissimilarity of compounds into account. We illustrate the use of the package with examples on simulated and real datasets. By providing the R package chemodiv for quantifying multiple aspects of chemodiversity, we hope to facilitate investigations of how chemodiversity varies across levels of biological organization, and its importance for the ecology and evolution of plants and other organisms.

Hybridization and introgression are prevalent in Southern European Erysimum (Brassicaceae) species

BACKGROUND AND AIMS

Hybridization is a common and important force in plant evolution. One of its outcomes is introgression - the transfer of small genomic regions from one taxon to another by hybridization and repeated backcrossing. This process is believed to be common in glacial refugia, where range expansions and contractions can lead to cycles of sympatry and isolation, creating conditions for extensive hybridization and introgression. Polyploidization is another genome-wide process with a major influence on plant evolution. Both hybridization and polyploidization can have complex effects on plant evolution. However, these effects are often difficult to understand in recently evolved species complexes.

METHODS

We combined flow cytometry, analyses of transcriptomic sequences and pollen tube growth assays to investigate the consequences of polyploidization, hybridization and introgression on the recent evolution of several Erysimum (Brassicaceae) species from the South of the Iberian Peninsula, a well-known glacial refugium. This species complex differentiated in the last 2 million years, and its evolution has been hypothesized to be determined mainly by polyploidization, interspecific hybridization and introgression.

KEY RESULTS

Our results support a scenario of widespread hybridization involving both extant and 'ghost' taxa. Several taxa studied here, most notably those with purple corollas, are polyploids, probably of allopolyploid origin. Moreover, hybridization in this group might be an ongoing phenomenon, as pre-zygotic barriers appeared weak in many cases.

CONCLUSIONS

The evolution of Erysimum spp. has been determined by hybridization to a large extent. Species with purple (polyploids) and yellow flowers (mostly diploid) exhibit a strong signature of introgression in their genomes, indicating that hybridization occurred regardless of colour and across ploidy levels. Although the adaptive value of such genomic exchanges remains unclear, our results demonstrate the significance of hybridization for plant diversification, which should be taken into account when studying plant evolution.

Quantification of plant cardenolides by HPLC, measurement of Na(+)/K(+)-ATPase inhibition activity, and characterization of target enzymes

The biosynthesis of cardiac glycosides, broadly classified as cardenolides and bufadienolides, has evolved repeatedly among flowering plants. Individual species can produce dozens or even hundreds of structurally distinct cardiac glycosides. Although all cardiac glycosides exhibit biological activity by inhibiting the function of the essential Na⁺/K⁺-ATPase in animal cells, they differ in their level of inhibitory activity. For within- and between-species comparisons of cardiac glycosides to address ecological and evolutionary questions, it is necessary to not only quantify their relative abundance, but also their effectiveness in inhibiting the activity of different animal Na⁺/K⁺-ATPases. Here we describe protocols for characterizing the amount and toxicity of cardenolides from plant samples and the degree of insect Na⁺/K⁺-ATPase tolerance to inhibition: (1) an HPLC-based assay to quantify the abundance of individual cardenolides in plant extracts, (2) an assay to quantify inhibition of Na⁺/K⁺-ATPase activity by plant extracts, and (3) extraction of insect Na⁺/K⁺-ATPases for inhibition assays.

Testing hypotheses of a coevolutionary key innovation reveals a complex suite of traits involved in defusing the mustard oil bomb

Coevolutionary interactions are responsible for much of the Earth's biodiversity, with key innovations driving speciation bursts on both sides of the interaction. One persistent question is whether macroevolutionary traits identified as key innovations accurately predict functional performance and selection dynamics within species, as this necessitates characterizing their function, investigating their fitness consequences, and exploring the selection dynamics acting upon them. Here, we used CRISPR-Cas9 mediating nonhomologous end joining (NHEJ) in the butterfly species Pieris brassicae to knock out and directly assess the function and fitness impacts of nitrile specifier protein (NSP) and major allergen (MA). These are two closely related genes that facilitate glucosinolate (GSL) detoxification capacity, which is a key innovation in mustard feeding Pierinae butterflies. We find NSP and MA are both required for survival on plants containing GSLs, with expression differences arising in response to variable GSL profiles, concordant with detoxification performance. Importantly, this concordance was only observed when using natural host plants, likely reflecting the complexity of how these enzymes interact with natural plant variation in GSLs and myrosinases. Finally, signatures of positive selection for NSP and MA were detected across Pieris species, consistent with these genes' importance in recent coevolutionary interactions. Thus, the war between these butterflies and their host plants involves more than the mere presence of chemical defenses and detoxification mechanisms, as their regulation and activation represent key components of complex interactions. We find that inclusion of these dynamics, in ecologically relevant assays, is necessary for coevolutionary insights in this system and likely others.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Lack of ITS sequence homogenization in Erysimum species (Brassicaceae) with different ploidy levels

The internal transcribed spacers (ITS) exhibit concerted evolution by the fast homogenization of these sequences at the intragenomic level. However, the rate and extension of this process are unclear and might be conditioned by the number and divergence of the different ITS copies. In some cases, such as hybrid species and polyploids, ITS sequence homogenization appears incomplete, resulting in multiple haplotypes within the same organism. Here, we studied the dynamics of concerted evolution in 85 individuals of seven plant species of the genus Erysimum (Brassicaceae) with multiple ploidy levels. We estimated the rate of concerted evolution and the degree of sequence homogenization separately for ITS1 and ITS2 and whether these varied with ploidy. Our results showed incomplete sequence homogenization, especially for polyploid samples, indicating a lack of concerted evolution in these taxa. Homogenization was usually higher in ITS2 than in ITS1, suggesting that concerted evolution operates more efficiently on the former. Furthermore, the hybrid origin of several species appears to contribute to the maintenance of high haplotype diversity, regardless of the level of ploidy. These findings indicate that sequence homogenization of ITS is a dynamic and complex process that might result in varying intra- and inter-genomic diversity levels.

Genomic insights into the evolution of plant chemical defense

Plant trait evolution can be impacted by common mechanisms of genome evolution, including whole-genome and small-scale duplication, rearrangement, and selective pressures. With the increasing accessibility of genome sequencing for non-model species, comparative studies of trait evolution among closely related or divergent lineages have supported investigations into plant chemical defense. Plant defensive compounds include major chemical classes, such as terpenoids, alkaloids, and phenolics, and are used in primary and secondary plant functions. These include the promotion of plant health, facilitation of pollination, defense against pathogens, and responses to a rapidly changing climate. We discuss mechanisms of genome evolution and use examples from recent studies to impress a stronger understanding of the link between genotype and phenotype as it relates to the evolution of plant chemical defense. We conclude with considerations for how to leverage genomics, transcriptomics, metabolomics, and functional assays for studying the emergence and evolution of chemical defense systems.

Using interdisciplinary, phylogeny-guided approaches to understand the evolution of plant metabolism

To cope with relentless environmental pressures, plants produce an arsenal of structurally diverse chemicals, often called specialized metabolites. These lineage-specific compounds are derived from the simple building blocks made by ubiquitous core metabolic pathways. Although the structures of many specialized metabolites are known, the underlying metabolic pathways and the evolutionary events that have shaped the plant chemical diversity landscape are only beginning to be understood. However, with the advent of multi-omics data sets and the relative ease of studying pathways in previously intractable non-model species, plant specialized metabolic pathways are now being systematically identified. These large datasets also provide a foundation for comparative, phylogeny-guided studies of plant metabolism. Comparisons of metabolic traits and features like chemical abundances, enzyme activities, or gene sequences from phylogenetically diverse plants provide insights into how metabolic pathways evolved. This review highlights the power of studying evolution through the lens of comparative biochemistry, particularly how placing metabolism into a phylogenetic context can help a researcher identify the metabolic innovations enabling the evolution of structurally diverse plant metabolites.

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.

Ethylene modulates translation dynamics in Arabidopsis under submergence via GCN2 and EIN2

General translational repression is a key process that reduces energy consumption under hypoxia. Here, we show that plant stress-activated general control nonderepressible 2 (GCN2) was activated to regulate the reduction in polysome loading during submergence in Arabidopsis. GCN2 signaling was activated by ethylene under submergence. GCN2 activity was reduced in etr1-1, but not in ein2-5 or eil1ein3, under submergence, suggesting that GCN2 activity is regulated by a noncanonical ethylene signaling pathway. Polysome loading was not reduced in ein2-5 under submergence, implying that ethylene modulates translation via both EIN2 and GCN2. Transcriptomic analysis demonstrated that EIN2 and GCN2 regulate not only general translational repression but also translational enhancement of specific mRNAs under submergence. Together, these results demonstrate that during submergence, entrapped ethylene triggers GCN2 and EIN2 to regulate translation dynamics and ensure the translation of stress response proteins.

Orchid Phylotranscriptomics: The Prospects of Repurposing Multi-Tissue Transcriptomes for Phylogenetic Analysis and Beyond

The Orchidaceae is rivaled only by the Asteraceae as the largest plant family, with the estimated number of species exceeding 25,000 and encompassing more than 700 genera. To gain insights into the mechanisms driving species diversity across both global and local scales, well-supported phylogenies targeting different taxonomic groups and/or geographical regions will be crucial. High-throughput sequencing technologies have revolutionized the field of molecular phylogenetics by simplifying the process of obtaining genome-scale sequence data. Consequently, there has been an explosive growth of such data in public repositories. Here we took advantage of this unprecedented access to transcriptome data from predominantly non-phylogenetic studies to assess if it can be repurposed to gain rapid and accurate phylogenetic insights across the orchids. Exhaustive searches revealed transcriptomic data for more than 100 orchid species spanning 5 subfamilies, 13 tribes, 21 subtribes, and 50 genera that were amendable for exploratory phylotranscriptomic analysis. Next, we performed re-assembly of the transcriptomes before strategic selection of the final samples based on a gene completeness evaluation. Drawing on these data, we report phylogenetic analyses at both deep and shallow evolutionary scales via maximum likelihood and shortcut coalescent species tree methods. In this perspective, we discuss some key outcomes of this study and conclude by highlighting other complementary, albeit rarely explored, insights beyond phylogenetic analysis that repurposed multi-tissue transcriptome can offer.

Evidence for tissue-specific defense-offense interactions between milkweed and its community of specialized herbivores

Coevolution between plants and herbivores often involves escalation of defense-offense strategies, but attack by multiple herbivores may obscure the match of plant defense to any one attacker. As herbivores often specialize on distinct plant parts, we hypothesized that defense-offense interactions in coevolved systems may become physiologically and evolutionarily compartmentalized between plant tissues. We report that roots, leaves, flower buds and seeds of the tropical milkweed (Asclepias curassavica) show increasing concentrations of cardenolide toxins acropetally, with latex showing the highest concentration. In vitro assays of the physiological target of cardenolides, the Na⁺ /K⁺ -ATPase (hereafter 'sodium pump'), of three specialized milkweed herbivores (root-feeding Tetraopes tetrophthalmus, leaf-feeding Danaus plexippus, and seed-feeding Oncopeltus fasciatus) show that they are proportionally tolerant to the cardenolide concentrations of the tissues they eat. Indeed, molecular substitutions in the insects' sodium pumps predicted their tolerance to toxins from their target tissues. Nonetheless, the relative inhibition of the sodium pumps of these specialists by the concentration vs. composition (inhibition controlled for concentration, what we term 'potency') of cardenolides from their target vs. non-target plant tissues revealed different degrees of insect adaptation to tissue-specific toxins. In addition, a trade-off between toxin concentration and potency emerged across plant tissues, potentially reflecting coevolutionary history or plant physiological constraints. Our findings suggest that tissue-specific coevolutionary dynamics may be proceeding between the plant and its specialized community of herbivores. This novel finding may be common in nature, contributing to ways in which coevolution proceeds in multi-species communities.

Insights into glucosinolate accumulation and metabolic pathways in Isatis indigotica Fort

BACKGROUND

Glucosinolates (GSLs) play important roles in defending against exogenous damage and regulating physiological activities in plants. However, GSL accumulation patterns and molecular regulation mechanisms are largely unknown in Isatis indigotica Fort.

RESULTS

Ten GSLs were identified in I. indigotica, and the dominant GSLs were epiprogoitrin (EPI) and indole-3-methyl GSL (I3M), followed by progoitrin (PRO) and gluconapin (GNA). The total GSL content was highest (over 20 μmol/g) in reproductive organs, lowest (less than 1.0 μmol/g) in mature organs, and medium in fresh leaves (2.6 μmol/g) and stems (1.5 μmol/g). In the seed germination process, the total GSL content decreased from 27.2 μmol/g (of seeds) to 2.7 μmol/g (on the 120th day) and then increased to 4.0 μmol/g (180th day). However, the content of indole GSL increased rapidly in the first week after germination and fluctuated between 1.13 μmol/g (28th day) and 2.82 μmol/g (150th day). Under the different elicitor treatments, the total GSL content increased significantly, ranging from 2.9-fold (mechanical damage, 3 h) to 10.7-fold (MeJA, 6 h). Moreover, 132 genes were involved in GSL metabolic pathways. Among them, no homologs of AtCYP79F2 and AtMAM3 were identified, leading to a distinctive GSL profile in I. indigotica. Furthermore, most genes involved in the GSL metabolic pathway were derived from tandem duplication, followed by dispersed duplication and segmental duplication. Purifying selection was observed, although some genes underwent relaxed selection. In addition, three tandem-arrayed GSL-OH genes showed different expression patterns, suggesting possible subfunctionalization during evolution.

CONCLUSIONS

Ten different GSLs with their accumulation patterns and 132 genes involved in the GSL metabolic pathway were explored, which laid a foundation for the study of GSL metabolism and regulatory mechanisms in I. indigotica.

Acropetal and basipetal cardenolide transport in Erysimum cheiranthoides (wormseed wallflower)

Plant specialized metabolites are often subject to within-plant transport and have tissue-specific distribution patterns. Among plants in the Brassicaceae, the genus Erysimum is unique in producing not only glucosinolates but also cardenolides. Ten cardenolides were detected with varying abundance in different tissues of Erysimum cheiranthoides L (Brassicaceae; wormseed wallflower). As is predicted by the optimal defense theory, cardenolides were most abundant in young leaves and reproductive tissues. The lowest concentrations were observed in senescing leaves and roots. Crosses between wildtype E. cheiranthoides and a mutant line with an altered cardenolide profile showed that the seed cardenolide phenotype is determined entirely by the maternal genotype. Prior to the development of the first true leaves, seedling cotyledons also had the maternal cardenolide profile. Hypocotyl grafting experiments showed that the root cardenolide profile is determined entirely by the aboveground plant genotype. In further grafting experiments, there was no evidence of cardenolide transport into the leaves, but a mixed cardenolide profile was observed in the stems and inflorescences of plants that had been grafted at vegetative and flowering growth stages, respectively. Together, these results indicate that E. cheiranthoides leaves are a site of cardenolide biosynthesis.

Intraspecific variation in plant-associated herbivore communities is phylogenetically structured in Brassicaceae

As a result of co‐evolution between plants and herbivores, related plants often interact with similar herbivore communities. Variation in plant–herbivore interactions is determined by variation in underlying functional traits and by ecological and stochastic processes. Hence, typically, only a subset of possible interactions is realised on individual plants. We show that insect herbivore communities assembling on individual plants are structured by plant phylogeny among 12 species in two phylogenetic lineages of Brassicaceae. This community sorting to plant phylogeny was retained when splitting the community according to herbivore feeding guilds. Relative abundance of herbivores as well as the size of the community structured community dissimilarity among plant species. Importantly, the amount of intraspecific variation in realised plant–herbivore interactions is also phylogenetically structured. We argue that variability in realised interactions that are not directly structured by plant traits is ecologically relevant and must be considered in the evolution of plant defences.

Selenium mediated host plant-mite conflict: defense and adaptation

BACKGROUND

Selenium has shown effectiveness in protecting plants from herbivores. However, some insects have evolved adaptability to selenium.

RESULTS

Selenium accumulation in host plants protected them against spider mite feeding. Selenium showed toxic effects on spider mites by reducing growth and interfering with reproduction. After 40 generations on selenium-rich plants, a Tetranychus cinnabarinus strain (Tc-Se) developed adaptability to selenium, with an increased rate of population growth and enhanced ability for selenium metabolism. The high expression of two genes (GSTd07 and SPS1) in the selenium metabolism pathway might be involved in selenium metabolism in spider mites. After GSTd07 and SPS1 were silenced, the selenium adaptability decreased. Recombinant GSTd07 protein promoted the reaction between sodium selenite and glutathione (GSH) and increased the production of sodium selenite metabolites. The results indicated that GSTd07 was involved in the first step of selenium metabolism.

CONCLUSION

Plants can resist spider mite feeding by accumulating selenium. Spider mites subjected to long-term selenium exposure can adapt to selenium by increasing the expression of key genes involved in selenium metabolism. These results elucidate the mechanism of the interaction between mites and host plants mediated by selenium. This study of the interaction between selenium-mediated host plants and spider mites may lead to the development of new and less toxic methods for the prevention and control of spider mites. © 2021 Society of Chemical Industry.

Glucosinolate profiles and phylogeny in Barbarea compared to other tribe Cardamineae (Brassicaceae) and Reseda (Resedaceae), based on a library of ion trap HPLC-MS/MS data of reference desulfoglucosinolates

A library of ion trap MS2 spectra and HPLC retention times reported here allowed distinction in plants of at least 70 known glucosinolates (GSLs) and some additional proposed GSLs. We determined GSL profiles of selected members of the tribe Cardamineae (Brassicaceae) as well as Reseda (Resedaceae) used as outgroup in evolutionary studies. We included several accessions of each species and a range of organs, and paid attention to minor peaks and GSLs not detected. In this way, we obtained GSL profiles of Barbarea australis, Barbarea grayi, Planodes virginica selected for its apparent intermediacy between Barbarea and the remaining tribe and family, and Rorippa sylvestris and Nasturtium officinale, for which the presence of acyl derivatives of GSLs was previously untested. We also screened Armoracia rusticana, with a remarkably diverse GSL profile, the emerging model species Cardamine hirsuta, for which we discovered a GSL polymorphism, and Reseda luteola and Reseda odorata. The potential for aliphatic GSL biosynthesis in Barbarea vulgaris was of interest, and we subjected P-type and G-type B. vulgaris to several induction regimes in an attempt to induce aliphatic GSL. However, aliphatic GSLs were not detected in any of the B. vulgaris types. We characterized the investigated chemotypes phylogenetically, based on nuclear rDNA internal transcribed spacer (ITS) sequences, in order to understand their relation to the species B. vulgaris in general, and found them to be representative of the species as it occurs in Europe, as far as documented in available ITS-sequence repositories. In short, we provide GSL profiles of a wide variety of tribe Cardamineae plants and conclude aliphatic GSLs to be absent or below our limit of detection in two major evolutionary lines of B. vulgaris. Concerning analytical chemistry, we conclude that availability of authentic reference compounds or reference materials is critical for reliable GSL analysis and characterize two publicly available reference materials: seeds of P. virginica and N. officinale.

Comparison of glucosinolate diversity in the crucifer tribe Cardamineae and the remaining order Brassicales highlights repetitive evolutionary loss and gain of biosynthetic steps

We review glucosinolate (GSL) diversity and analyze phylogeny in the crucifer tribe Cardamineae as well as selected species from Brassicaceae (tribe Brassiceae) and Resedaceae. Some GSLs occur widely, while there is a scattered distribution of many less common GSLs, tentatively sorted into three classes: ancient, intermediate and more recently evolved. The number of conclusively identified GSLs in the tribe (53 GSLs) constitute 60% of all GSLs known with certainty from any plant (89 GSLs) and apparently unique GSLs in the tribe constitute 10 of those GSLs conclusively identified (19%). Intraspecific, qualitative GSL polymorphism is known from at least four species in the tribe. The most ancient GSL biosynthesis in Brassicales probably involved biosynthesis from Phe, Val, Leu, Ile and possibly Trp, and hydroxylation at the β-position. From a broad comparison of families in Brassicales and tribes in Brassicaceae, we estimate that a common ancestor of the tribe Cardamineae and the family Brassicaceae exhibited GSL biosynthesis from Phe, Val, Ile, Leu, possibly Tyr, Trp and homoPhe (ancient GSLs), as well as homologs of Met and possibly homoIle (intermediate age GSLs). From the comparison of phylogeny and GSL diversity, we also suggest that hydroxylation and subsequent methylation of indole GSLs and usual modifications of Met-derived GSLs (formation of sulfinyls, sulfonyls and alkenyls) occur due to conserved biochemical mechanisms and was present in a common ancestor of the family. Apparent loss of homologs of Met as biosynthetic precursors was deduced in the entire genus Barbarea and was frequent in Cardamine (e.g. C. pratensis, C. diphylla, C. concatenata, possibly C. amara). The loss was often associated with appearance of significant levels of unique or rare GSLs as well as recapitulation of ancient types of GSLs. Biosynthetic traits interpreted as de novo evolution included hydroxylation at rare positions, acylation at the thioglucose and use of dihomoIle and possibly homoIle as biosynthetic precursors. Biochemical aspects of the deduced evolution are discussed and testable hypotheses proposed. Biosyntheses from Val, Leu, Ile, Phe, Trp, homoPhe and homologs of Met are increasingly well understood, while GSL biosynthesis from mono- and dihomoIle is poorly understood. Overall, interpretation of known diversity suggests that evolution of GSL biosynthesis often seems to recapitulate ancient biosynthesis. In contrast, unprecedented GSL biosynthetic innovation seems to be rare.

Cardenolides, toxicity, and the costs of sequestration in the coevolutionary interaction between monarchs and milkweeds

For highly specialized insect herbivores, plant chemical defenses are often co-opted as cues for oviposition and sequestration. In such interactions, can plants evolve novel defenses, pushing herbivores to trade off benefits of specialization with costs of coping with toxins? We tested how variation in milkweed toxins (cardenolides) impacted monarch butterfly (Danaus plexippus) growth, sequestration, and oviposition when consuming tropical milkweed (Asclepias curassavica), one of two critical host plants worldwide. The most abundant leaf toxin, highly apolar and thiazolidine ring-containing voruscharin, accounted for 40% of leaf cardenolides, negatively predicted caterpillar growth, and was not sequestered. Using whole plants and purified voruscharin, we show that monarch caterpillars convert voruscharin to calotropin and calactin in vivo, imposing a burden on growth. As shown by in vitro experiments, this conversion is facilitated by temperature and alkaline pH. We next employed toxin-target site experiments with isolated cardenolides and the monarch's neural Na⁺/K⁺-ATPase, revealing that voruscharin is highly inhibitory compared with several standards and sequestered cardenolides. The monarch's typical >50-fold enhanced resistance to cardenolides compared with sensitive animals was absent for voruscharin, suggesting highly specific plant defense. Finally, oviposition was greatest on intermediate cardenolide plants, supporting the notion of a trade-off between benefits and costs of sequestration for this highly specialized herbivore. There is apparently ample opportunity for continued coevolution between monarchs and milkweeds, although the diffuse nature of the interaction, due to migration and interaction with multiple milkweeds, may limit the ability of monarchs to counteradapt.

Engineering insect resistance using plant specialized metabolites

Plants in nature are protected against insect herbivory by a wide variety of specialized metabolites. Although insect herbivores generally tolerate the defensive metabolites of their preferred host plants, the presence of additional chemical defenses in otherwise closely related plant species can nevertheless provide resistance. This chemical resistance to insect herbivory can be enhanced by genetic engineering to increase the production of endogenous defensive metabolites, modify existing biochemical pathways, or move the biosynthesis of entirely new classes of specialized metabolites into recipient plants. However, current plant genetic engineering strategies are limited by insufficient knowledge of the biosynthetic pathways of plant specialized metabolism, unintended side-effects that result from redirecting plant metabolism, inadequate transgene construction and delivery methods, and requirements for tissue-specific production of defensive metabolites to enhance herbivore resistance.

Sustainable Use of Bioactive Compounds from Solanum Tuberosum and Brassicaceae Wastes and by-Products for Crop Protection-A Review

Defatted seed meals of oleaginous Brassicaceae, such as Eruca sativa, and potato peel are excellent plant matrices to recover potentially useful biomolecules from industrial processes in a circular strategy perspective aiming at crop protection. These biomolecules, mainly glycoalkaloids and phenols for potato and glucosinolates for Brassicaceae, have been proven to be effective against microbes, fungi, nematodes, insects, and even parasitic plants. Their role in plant protection is overviewed, together with the molecular basis of their synthesis in plant, and the description of their mechanisms of action. Possible genetic and biotechnological strategies are presented to increase their content in plants. Genetic mapping and identification of closely linked molecular markers are useful to identify the loci/genes responsible for their accumulation and transfer them to elite cultivars in breeding programs. Biotechnological approaches can be used to modify their allelic sequence and enhance the accumulation of the bioactive compounds. How the global challenges, such as reducing agri-food waste and increasing sustainability and food safety, could be addressed through bioprotector applications are discussed here.

A story from the Miocene: Clock-dated phylogeny of Sisymbrium L. (Sisymbrieae, Brassicaceae)

Morphological variability and imprecise generic boundaries have hindered systematic, taxonomical, and nomenclatural studies of Sisymbrium L. (Brassicaceae, Sisymbrieae DC.). The members of this almost exclusively Old-World genus grow mostly on highly porous substrates across open steppe, semidesert, or ruderal habitats in the temperate zone of the Northern Hemisphere and African subtropics. The present study placed the biological history of Sisymbrium L. into time and space and rendered the tribus Sisymbrieae as monotypic. Five nuclear-encoded and three chloroplast-encoded loci of approximately 85% of all currently accepted species were investigated. Several accessions per species covering their whole distribution range allowed for a more representative assessment of intraspecific genetic diversity. In the light of fossil absence, the impact of different secondary calibration methods and taxon sets on time spans was tested, and we showed that such a combinatorial nested dating approach is beneficial. Multigene phylogeny accompanied with a time divergence estimation analysis placed the onset and development of this tribus into the western Irano-Turanian floristic region during the Miocene. Continuous increase in continentality and decrease in temperatures promoted the diversity of the Sisymbrieae, which invaded the open grasslands habitats in Eurasia, Mediterranean, and South Africa throughout the Pliocene and Pleistocene. Our results support the assumption of the Irano-Turanian region as a biodiversity reservoir for adjacent regions.

Evolutionary history and ecology shape the diversity and abundance of phytochemical arsenals across monkeyflowers

We examine the extent to which phylogenetic effects and ecology are associated with macroevolutionary patterns of phytochemical defence production across the Mimulus phylogeny. We grew plants from 21 species representing the five major sections of the Mimulus phylogeny in a common garden to assess how the arsenals (NMDS groupings) and abundances (concentrations) of a phytochemical defence, phenylpropanoid glycosides (PPGs), vary across the phylogeny. Very few PPGs are widespread across the genus, but many are common to multiple sections of the genus. Phytochemical arsenals cluster among sections in an NMDS and are not associated with total concentration of PPGs. There is a strong phylogenetic signal for phytochemical arsenal composition across the Mimulus genus, whereas ecological variables such as growing season length, latitude, and elevation do not significantly influence arsenal. In contrast, there is little phylogenetic signal for total PPG concentration, and this trait is significantly influenced by several ecological factors. Phytochemical arsenals and abundances are influenced by plant life history form. Both phylogenetic effects and ecology are related to phytochemical patterns across species, albeit in different ways. The independence of phytochemical defence concentrations from arsenal compositions indicates that these aspects of defence may continue to evolve independently of one another.

Evolutionary changes in the glucosinolate biosynthetic capacity in species representing Capsella, Camelina and Neslia genera

Glucosinolates are unique thioglucosides that evolved in the order Brassicales. These compounds function in plant adaptation to the environment, including combating plant pathogens, herbivore deterrence and abiotic stress tolerance. In line with their defensive functions glucosinolates usually accumulate constitutively in relatively high amounts in all tissues of Brassicaceae plants. Here we performed glucosinolate analysis in different organs of selected species representing Capsella, Camelina and Neslia genera, which similarly as the model plant Arabidopsis thaliana belong to the Camelineae tribe. We also identified orthologs of A. thaliana glucosinolate biosynthetic genes in the published genomes of some of the investigated species. Subsequent gene expression and phylogenetic analyses enabled us an insight into the evolutionary changes in the transcription of these genes and in the sequences of respective proteins that occurred within the Camelineae tribe. Our results indicated that glucosinolates are highly abundant in siliques and roots of the investigated species but hardly, if at all, produced in leaves. In addition to this unusual tissular distribution we revealed reduced structural diversity of methionine-derived aliphatic glucosinolates (AGs) with elevated accumulation of rare long chain AGs. This preference seems to correlate with evolutionary changes in genes encoding methylthioalkylmalate synthases that are responsible for the elongation of AG side chains. Finally, our results indicate that the biosynthetic pathway for tryptophan-derived indolic glucosinolates likely lost its main functions in immunity and resistance towards sucking insects and is on its evolutionary route to be shut off in the investigated species.

Less Is More: a Mutation in the Chemical Defense Pathway of Erysimum cheiranthoides (Brassicaceae) Reduces Total Cardenolide Abundance but Increases Resistance to Insect Herbivores

Erysimum cheiranthoides L (Brassicaceae; wormseed wallflower) accumulates not only glucosinolates, which are characteristic of the Brassicaceae, but also abundant and diverse cardenolides. These steroid toxins, primarily glycosylated forms of digitoxigenin, cannogenol, and strophanthidin, inhibit the function of essential Na⁺/K⁺-ATPases in animal cells. We screened a population of 659 ethylmethanesulfonate-mutagenized E. cheiranthoides plants to identify isolates with altered cardenolide profiles. One mutant line exhibited 66% lower cardenolide content, resulting from greatly decreased cannogenol and strophanthidin glycosides, partially compensated for by increases in digitoxigenin glycosides. This phenotype was likely caused by a single-locus recessive mutation, as evidenced by a wildtype phenotype of F1 plants from a backcross, a 3:1 wildtype:mutant segregation in the F2 generation, and genetic mapping of the altered cardenolide phenotype to one position in the genome. The mutation created a more even cardenolide distribution, decreased the average cardenolide polarity, but did not impact most glucosinolates. Growth of generalist herbivores from two feeding guilds, Myzus persicae Sulzer (Hemiptera: Aphididae; green peach aphid) and Trichoplusia ni Hübner (Lepidoptera: Noctuidae; cabbage looper), was decreased on the mutant line compared to wildtype. Both herbivores accumulated cardenolides in proportion to the plant content, with T. ni accumulating higher total concentrations than M. persicae. Helveticoside, a relatively abundant cardenolide in E. cheiranthoides, was not detected in M. persicae feeding on these plants. Our results support the hypothesis that increased digitoxigenin glycosides provide improved protection against M. persicae and T. ni, despite an overall decrease in cardenolide content of the mutant line.

Independent evolution of ancestral and novel defenses in a genus of toxic plants (Erysimum, Brassicaceae)

Phytochemical diversity is thought to result from coevolutionary cycles as specialization in herbivores imposes diversifying selection on plant chemical defenses. Plants in the speciose genus Erysimum (Brassicaceae) produce both ancestral glucosinolates and evolutionarily novel cardenolides as defenses. Here we test macroevolutionary hypotheses on co-expression, co-regulation, and diversification of these potentially redundant defenses across this genus. We sequenced and assembled the genome of E. cheiranthoides and foliar transcriptomes of 47 additional Erysimum species to construct a phylogeny from 9868 orthologous genes, revealing several geographic clades but also high levels of gene discordance. Concentrations, inducibility, and diversity of the two defenses varied independently among species, with no evidence for trade-offs. Closely related, geographically co-occurring species shared similar cardenolide traits, but not glucosinolate traits, likely as a result of specific selective pressures acting on each defense. Ancestral and novel chemical defenses in Erysimum thus appear to provide complementary rather than redundant functions.

Plants are often attacked by insects and other herbivores. As a result, they have evolved to defend themselves by producing many different chemicals that are toxic to these pests. As producing each chemical costs energy, individual plants often only produce one type of chemical that is targeted towards their main herbivore. Related species of plants often use the same type of chemical defense so, if a particular herbivore gains the ability to cope with this chemical, it may rapidly become an important pest for the whole plant family.

To escape this threat, some plants have gained the ability to produce more than one type of chemical defense. Wallflowers, for example, are a group of plants in the mustard family that produce two types of toxic chemicals: mustard oils, which are common in most plants in this family; and cardenolides, which are an innovation of the wallflowers, and which are otherwise found only in distantly related plants such as foxglove and milkweed. The combination of these two chemical defenses within the same plant may have allowed the wallflowers to escape attacks from their main herbivores and may explain why the number of wallflower species rapidly increased within the last two million years.

Züst et al. have now studied the diversity of mustard oils and cardenolides present in many different species of wallflower. This analysis revealed that almost all of the tested wallflower species produced high amounts of both chemical defenses, while only one species lacked the ability to produce cardenolides. The levels of mustard oils had no relation to the levels of cardenolides in the tested species, which suggests that the regulation of these two defenses is not linked. Furthermore, Züst et al. found that closely related wallflower species produced more similar cardenolides, but less similar mustard oils, to each other. This suggests that mustard oils and cardenolides have evolved independently in wallflowers and have distinct roles in the defense against different herbivores.

The evolution of insect resistance to pesticides and other toxins is an important concern for agriculture. Applying multiple toxins to crops at the same time is an important strategy to slow the evolution of resistance in the pests. The findings of Züst et al. describe a system in which plants have naturally evolved an equivalent strategy to escape their main herbivores. Understanding how plants produce multiple chemical defenses, and the costs involved, may help efforts to breed crop species that are more resistant to herbivores and require fewer applications of pesticides.