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

Heat waves induce milkweed resistance to a specialist herbivore via increased toxicity and reduced nutrient content

Over the last decade, a large effort has been made to understand how extreme climate events disrupt species interactions. Yet, it is unclear how these events affect plants and herbivores directly, via metabolic changes, and indirectly, via their subsequent altered interaction. We exposed common milkweed (Asclepias syriaca) and monarch caterpillars (Danaus plexippus) to control (26:14°C, day:night) or heat wave (HW) conditions (36:24°C, day:night) for 4 days and then moved each organism to a new control or HW partner to disentangle the direct and indirect effects of heat exposure on each organism. We found that the HW directly benefited plants in terms of growth and defence expression (increased latex exudation and total cardenolides) and insect her'bivores through faster larval development. Conversely, indirect HW effects caused both plant latex and total cardenolides to decrease after subsequent herbivory. Nonetheless, increasing trends of more toxic cardenolides and lower leaf nutritional quality after herbivory by HW caterpillars likely led to reduced plant damage compared to controls. Our findings reveal that indirect impacts of HWs may play a greater role in shaping plant-herbivore interactions via changes in key physiological traits, providing valuable understanding of how ecological interactions may proceed in a changing world.

Four enzymes control natural variation in the steroid core of Erysimum cardenolides

Plants commonly produce families of structurally related metabolites with similar defensive functions. This apparent redundancy raises the question of underlying molecular mechanisms and adaptive benefits of such chemical variation. Cardenolides, a class defensive compounds found in the wallflower genus Erysimum (L., Brassicaceae) and scattered across other plant families, show substantial structural variation, with glycosylation and hydroxylation being common modifications of a steroid core, which itself may vary in terms of stereochemistry and saturation. Through a combination of chemical mutagenesis and analysis of gene coexpression networks, we identified four enzymes involved in cardenolide biosynthesis in Erysimum that work together to determine stereochemistry at carbon 5 of the steroid core: Ec3βHSD, a 3β-hydroxysteroid dehydrogenase, Ec3KSI, a ketosteroid isomerase, EcP5βR2, a progesterone 5β-reductase, and EcDET2, a steroid 5α-reductase. We biochemically characterized the activity of these enzymes in vitro and generated CRISPR/Cas9 knockout lines to confirm activity in vivo. Cardenolide biosynthesis was not eliminated in any of the knockouts. Instead, mutant plants accumulated cardenolides with altered saturation and stereochemistry of the steroid core. Furthermore, we found variation in carbon 5 configuration among the cardenolides of 44 species of Erysimum, where the occurrence of some 5β-cardenolides is associated with the expression and sequence of P5βR2. This may have allowed Erysimum species to fine-tune their defensive profiles to target specific herbivore populations over the course of evolution.

SIGNIFICANCE STATEMENT

Plants use an array of toxic compounds to defend themselves from attack against insects and other herbivores. One mechanism through which plants may evolve more toxic compounds is through modifications to the structure of compounds they already produce. In this study, we show how plants in the wallflower genus Erysimum use four enzymes to fine-tune the structure of toxic metabolites called cardenolides. Natural variation in the sequence and expression of a single enzyme called progesterone 5β-reductase 2 partly explains the variation in cardenolides observed across the Erysimum genus. These alterations to cardenolide structure over the course of evolution suggests that there may be context-dependent benefits to Erysimum to invest in one cardenolide variant over another.

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.

Overexpression and RNAi-mediated Knockdown of Two 3β-hydroxy-Δ5-steroid dehydrogenase Genes in Digitalis lanata Shoot Cultures Reveal Their Role in Cardenolide Biosynthesis

3β-hydroxy-Δ5-steroid dehydrogenases (3βHSDs) are supposed to be involved in 5β-cardenolide biosynthesis. Here, a novel 3βHSD (Dl3βHSD2) was isolated from Digitalis lanata shoot cultures and expressed in E. coli. Recombinant Dl3βHSD1 and Dl3βHSD2 shared 70% amino acid identity, reduced various 3-oxopregnanes and oxidised 3-hydroxypregnanes, but only rDl3βHSD2 converted small ketones and secondary alcohols efficiently. To explain these differences in substrate specificity, we established homology models using borneol dehydrogenase of Salvia rosmarinus (6zyz) as the template. Hydrophobicity and amino acid residues in the binding pocket may explain the difference in enzyme activities and substrate preferences. Compared to Dl3βHSD1, Dl3βHSD2 is weakly expressed in D. lanata shoots. High constitutive expression of Dl3βHSDs was realised by Agrobacterium-mediated transfer of Dl3βHSD genes fused to the CaMV-35S promotor into the genome of D. lanata wild type shoot cultures. Transformed shoots (35S:Dl3βHSD1 and 35S:Dl3βHSD2) accumulated less cardenolides than controls. The levels of reduced glutathione (GSH), which is known to inhibit cardenolide formation, were higher in the 35S:Dl3βHSD1 lines than in the controls. In the 35S:Dl3βHSD1 lines cardenolide levels were restored after adding of the substrate pregnane-3,20-dione in combination with buthionine-sulfoximine (BSO), an inhibitor of GSH formation. RNAi-mediated knockdown of the Dl3βHSD1 yielded several shoot culture lines with strongly reduced cardenolide levels. In these lines, cardenolide biosynthesis was fully restored after addition of the downstream precursor pregnan-3β-ol-20-one, whereas upstream precursors such as progesterone had no effect, indicating that no shunt pathway could overcome the Dl3βHSD1 knockdown. These results can be taken as the first direct proof that Dl3βHSD1 is indeed involved in 5β-cardenolide biosynthesis.

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.

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.