A novel 2-oxoacid-dependent dioxygenase involved in the formation of the goiterogenic 2-hydroxybut-3-enyl glucosinolate and generalist insect resistance in Arabidopsis,

A 2-oxoglutarate-dependent dioxygenase, GLUCORAPHASATIN SYNTHASE 1 (GRS1) is a major determinant for different aliphatic glucosinolates between radish and Chinese cabbage

Glucosinolates (GSLs) are secondary metabolites in Brassicaceae plants and play a defensive role against a variety of abiotic and biotic stresses. Also, it exhibits anti-cancer activity against cancer cell in human. Different profiles of aliphatic GSL compounds between radish and Chinese cabbage were previously reported. However, molecular details underlying the divergent profile between two species were not clearly understood. In this study, we found that major difference of aliphatic GSLs profiles between two species is determined by the dominantly expressed genes in first step of the secondary modification phase, which are responsible for enzymatic catalysis of methylthioalkyl-glucosinolate. For instance, active expression of GLUCORAPHASATIN SYNTHASE 1 (GRS1) gene in radish play an important role in the production of glucoraphasatin (GRH) and glucoraphenin (GRE), a major aliphatic GSLs in radish. Meanwhile, Chinese cabbage was found to merely produce glucoraphasatin (GRH), instead producing glucoraphanin (GRA) and gluconapin (GNP) due to the mere expression of GRS1 homologs and abundant expressions of FLAVIN-CONTAINING MONOOXYGENASES (FMO GS-OX) homologs in Chinese cabbage. In addition, we noticed that wounding treatment on leaf tissues substantially enhanced the production of aliphatic and benzenic GSLs in both Chinese cabbage and radish, indicating that GSLs are wound-induced defensive compounds in both Chinese cabbage and radish plants.

ELONGATED HYPOCOTYL 5 interacts with HISTONE DEACETYLASE 9 to suppress glucosinolate biosynthesis in Arabidopsis

Glucosinolates (GSLs) are defensive secondary metabolites produced by Brassicaceae species in response to abiotic and biotic stresses. The biosynthesis of GSL compounds and the expression of GSL-related genes are highly modulated by endogenous signals (i.e. circadian clocks) and environmental cues, such as temperature, light, and pathogens. However, the detailed mechanism by which light signaling influences GSL metabolism remains poorly understood. In this study, we found that a light-signaling factor, ELONGATED HYPOCOTYL 5 (HY5), was involved in the regulation of GSL content under light conditions in Arabidopsis (Arabidopsis thaliana). In hy5-215 mutants, the transcript levels of GSL pathway genes were substantially upregulated compared with those in wild-type (WT) plants. The content of GSL compounds was also substantially increased in hy5-215 mutants, whereas 35S::HY5-GFP/hy5-215 transgenic lines exhibited comparable levels of GSL-related transcripts and GSL content to those in WT plants. HY5 physically interacts with HISTONE DEACETYLASE9 and binds to the proximal promoter region of MYB29 and IMD1 to suppress aliphatic GSL biosynthetic processes. These results demonstrate that HY5 suppresses GSL accumulation during the daytime, thus properly modulating GSL content daily in Arabidopsis plants.

Allelic variation of TaABI5-A4 significantly affects seed dormancy in bread wheat

KEY MESSAGE

We identified a pivotal transcription factor TaABI5-A4 that is significantly associated with pre-harvest sprouting in wheat; its function in regulating seed dormancy was confirmed in transgenic rice. ABI5 is a critical transcription factor in regulation of crop seed maturation, dormancy, germination, and post-germination. Sixteen copies of homologous sequences of ABI5 were identified in Chinese wheat line Zhou 8425B. Cultivars of two haplotypes TaABI5-A4a and TaABI5-A4b showed significantly different seed dormancies. Based on two SNPs between the sequences of TaABI5-A4a and TaABI5-A4b, two complementary dominant sequence-tagged site (STS) markers were developed and validated in a natural population of 103 Chinese wheat cultivars and advanced lines and 200 recombinant inbred lines (RILs) derived from the Yangxiaomai/Zhongyou 9507 cross; the STS markers can be used efficiently and reliably to evaluate the dormancy of wheat seeds. The transcription level of TaABI5-A4b was significantly increased in TaABI5-A4a-GFP transgenic rice lines compared with that in TaABI5-A4b-GFP. The average seed germination index of TaABI5-A4a-GFP transgenic rice lines was significantly lower than those of TaABI5-A4b-GFP. In addition, seeds of TaABI5-A4a-GFP transgenic lines had higher ABA sensitivity and endogenous ABA content, lower endogenous GA content and plant height, and thicker stem internodes than those of TaABI5-A4b-GFP. Allelic variation of TaABI5-A4-affected wheat seed dormancy and the gene function was confirmed in transgenic rice. The transgenic rice lines of TaABI5-A4a and TaABI5-A4b had significantly different sensitivities to ABA and contents of endogenous ABA and GA in mature seeds, thereby influencing the seed dormancy, plant height, and stem internode length and diameter.

Adaptive modifications in plant sulfur metabolism over evolutionary time

Sulfur (S) is an essential element for life on Earth. Plants are able to take up and utilize sulfate (SO42-), the most oxidized inorganic form of S compounds on Earth, through the reductive S assimilatory pathway that couples with photosynthetic energy conversion. Organic S compounds are subsequently synthesized in plants and made accessible to animals, primarily as the amino acid methionine. Thus, plant S metabolism clearly has nutritional importance in the global food chain. S metabolites may be part of redox regulation and drivers of essential metabolic pathways as cofactors and prosthetic groups, such as Fe-S centers, CoA, thiamine, and lipoic acid. The evolution of the S metabolic pathways and enzymes reflects the critical importance of functional innovation and diversifications. Here we review the major evolutionary alterations that took place in S metabolism across different scales and outline research directions that may take advantage of understanding the evolutionary adaptations.

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.

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.

Tomato root specialized metabolites evolved through gene duplication and regulatory divergence within a biosynthetic gene cluster

Tremendous plant metabolic diversity arises from phylogenetically restricted specialized metabolic pathways. Specialized metabolites are synthesized in dedicated cells or tissues, with pathway genes sometimes colocalizing in biosynthetic gene clusters (BGCs). However, the mechanisms by which spatial expression patterns arise and the role of BGCs in pathway evolution remain underappreciated. In this study, we investigated the mechanisms driving acylsugar evolution in the Solanaceae. Previously thought to be restricted to glandular trichomes, acylsugars were recently found in cultivated tomato roots. We demonstrated that acylsugars in cultivated tomato roots and trichomes have different sugar cores, identified root-enriched paralogs of trichome acylsugar pathway genes, and characterized a key paralog required for root acylsugar biosynthesis, SlASAT1-LIKE (SlASAT1-L), which is nested within a previously reported trichome acylsugar BGC. Last, we provided evidence that ASAT1-L arose through duplication of its paralog, ASAT1, and was trichome-expressed before acquiring root-specific expression in the Solanum genus. Our results illuminate the genomic context and molecular mechanisms underpinning metabolic diversity in plants.

2-oxoglutarate-dependent dioxygenases and BAHD acyltransferases drive the structural diversification of orobanchol in Fabaceae plants

Strigolactones (SLs), a class of plant apocarotenoids, serve dual roles as rhizosphere-signaling molecules and plant hormones. Orobanchol, a major naturally occurring SL, along with its various derivatives, has been detected in the root exudates of plants of the Fabaceae family. Medicaol, fabacyl acetate, and orobanchyl acetate were identified in the root exudates of barrel medic (Medicago truncatula), pea (Pisum sativum), and cowpea (Vigna unguiculata), respectively. Although the biosynthetic pathway leading to orobanchol production has been elucidated, the biosynthetic pathways of the orobanchol derivatives have not yet been fully elucidated. Here, we report the identification of 2-oxoglutarate-dependent dioxygenases (DOXs) and BAHD acyltransferases responsible for converting orobanchol to these derivatives in Fabaceae plants. First, the metabolic pathways downstream of orobanchol were analyzed using substrate feeding experiments. Prohexadione, an inhibitor of DOX inhibits the conversion of orobanchol to medicaol in barrel medic. The DOX inhibitor also reduced the formation of fabacyl acetate and fabacol, a precursor of fabacyl acetate, in pea. Subsequently, we utilized a dataset based on comparative transcriptome analysis to select a candidate gene encoding DOX for medicaol synthase in barrel medic. Recombinant proteins of the gene converted orobanchol to medicaol. The candidate genes encoding DOX and BAHD acyltransferase for fabacol synthase and fabacol acetyltransferase, respectively, were selected by co-expression analysis in pea. The recombinant proteins of the candidate genes converted orobanchol to fabacol and acetylated fabacol. Furthermore, fabacol acetyltransferase and its homolog in cowpea acetylated orobanchol. The kinetics and substrate specificity analyses revealed high affinity and strict recognition of the substrates of the identified enzymes. These findings shed light on the molecular mechanisms underlying the structural diversity of SLs.

Genome-wide association studies identify loci controlling specialized seed metabolites in Arabidopsis

Plants synthesize specialized metabolites to facilitate environmental and ecological interactions. During evolution, plants diversified in their potential to synthesize these metabolites. Quantitative differences in metabolite levels of natural Arabidopsis (Arabidopsis thaliana) accessions can be employed to unravel the genetic basis for metabolic traits using genome-wide association studies (GWAS). Here, we performed metabolic GWAS on seeds of a panel of 315 A. thaliana natural accessions, including the reference genotypes C24 and Col-0, for polar and semi-polar seed metabolites using untargeted ultra-performance liquid chromatography-mass spectrometry. As a complementary approach, we performed quantitative trait locus (QTL) mapping of near-isogenic introgression lines between C24 and Col-0 for specific seed specialized metabolites. Besides common QTL between seeds and leaves, GWAS revealed seed-specific QTL for specialized metabolites, indicating differences in the genetic architecture of seeds and leaves. In seeds, aliphatic methylsulfinylalkyl and methylthioalkyl glucosinolates associated with the ALKENYL HYDROXYALKYL PRODUCING loci (GS-ALK and GS-OHP) on chromosome 4 containing alkenyl hydroxyalkyl producing 2 (AOP2) and 3 (AOP3) or with the GS-ELONG locus on chromosome 5 containing methylthioalkyl malate synthase (MAM1) and MAM3. We detected two unknown sulfur-containing compounds that were also mapped to these loci. In GWAS, some of the annotated flavonoids (kaempferol 3-O-rhamnoside-7-O-rhamnoside, quercetin 3-O-rhamnoside-7-O-rhamnoside) were mapped to transparent testa 7 (AT5G07990), encoding a cytochrome P450 75B1 monooxygenase. Three additional mass signals corresponding to quercetin-containing flavonols were mapped to UGT78D2 (AT5G17050). The association of the loci and associating metabolic features were functionally verified in knockdown mutant lines. By performing GWAS and QTL mapping, we were able to leverage variation of natural populations and parental lines to study seed specialized metabolism. The GWAS data set generated here is a high-quality resource that can be investigated in further studies.

Specificity and breadth of plant specialized metabolite-microbe interactions

Plant specialized metabolites shape plant interactions with the environment including plant-microbe interactions. While we often group compounds into generic classes, it is the precise structure of a compound that creates a specific role in plant-microbe or-pathogen interactions. Critically, the structure guides definitive targets in individual interactions, yet single compounds are not limited to singular mechanistic targets allowing them to influence interactions across broad ranges of attackers, from bacteria to fungi to animals. Further, the direction of the effect can be altered by counter evolution within the interacting organism leading to single compounds being both beneficial and detrimental. Thus, the benefit of a single compound to a host needs to be assessed by measuring the net benefit across all interactions while in each specific interaction. Factoring this complexity for single compounds in plant-microbe interactions with the massive expansion in our identification of specialized metabolite pathways means that we need systematic studies to classify the full breadth of activities. Only with this full biological knowledge we can develop mechanistic, ecological, and evolutionary models to understand how plant specialized metabolites fully influence plant-microbe and plant-biotic interactions more broadly.

Large-scale identification of novel transcriptional regulators of the aliphatic glucosinolate pathway in Arabidopsis

Aliphatic glucosinolates are a large group of plant secondary metabolites characteristic of Brassicaceae, including the model plant Arabidopsis. The diverse and complex degradation products of aliphatic glucosinolates contribute to plant responses to herbivory, pathogen attack, and environmental stresses. Most of the biosynthesis genes in the aliphatic glucosinolate pathway have been cloned in Arabidopsis, and the research focus has recently shifted to the regulatory mechanisms controlling aliphatic glucosinolate accumulation. Up till now, more than 40 transcriptional regulators have been identified as regulating the aliphatic glucosinolate pathway, but many more novel regulators likely remain to be discovered based on research evidence over the past decade. In the current study, we took a systemic approach to functionally test 155 candidate transcription factors in Arabidopsis identified by yeast one-hybrid assay, and successfully validated at least 30 novel regulators that could significantly influence the accumulation of aliphatic glucosinolates in our experimental set-up. We also showed that the regulators of the aliphatic glucosinolate pathway have balanced positive and negative effects, and glucosinolate metabolism and plant development can be coordinated. Our work is the largest scale effort so far to validate transcriptional regulators of a plant secondary metabolism pathway, and provides new insights into how the highly diverse plant secondary metabolism is regulated at the transcriptional level.

Genome-Wide Identification of 2-Oxoglutarate and Fe (II)-Dependent Dioxygenase (2ODD-C) Family Genes and Expression Profiles under Different Abiotic Stresses in Camellia sinensis (L.)

The 2-oxoglutarate and Fe (II)-dependent dioxygenase (2ODD-C) family of 2-oxoglutarate-dependent dioxygenases potentially participates in the biosynthesis of various metabolites under various abiotic stresses. However, there is scarce information on the expression profiles and roles of 2ODD-C genes in Camellia sinensis. We identified 153 Cs2ODD-C genes from C. sinensis, and they were distributed unevenly on 15 chromosomes. According to the phylogenetic tree topology, these genes were divided into 21 groups distinguished by conserved motifs and an intron/exon structure. Gene-duplication analyses revealed that 75 Cs2ODD-C genes were expanded and retained after WGD/segmental and tandem duplications. The expression profiles of Cs2ODD-C genes were explored under methyl jasmonate (MeJA), polyethylene glycol (PEG), and salt (NaCl) stress treatments. The expression analysis showed that 14, 13, and 49 Cs2ODD-C genes displayed the same expression pattern under MeJA and PEG treatments, MeJA and NaCl treatments, and PEG and NaCl treatments, respectively. A further analysis showed that two genes, Cs2ODD-C36 and Cs2ODD-C21, were significantly upregulated and downregulated after MeJA, PEG, and NaCl treatments, indicating that these two genes played positive and negative roles in enhancing the multi-stress tolerance. These results provide candidate genes for the use of genetic engineering technology to modify plants by enhancing multi-stress tolerance to promote phytoremediation efficiency.

Developing multifunctional crops by engineering Brassicaceae glucosinolate pathways

Glucosinolates (GSLs), found mainly in species of the Brassicaceae family, are one of the most well-studied classes of secondary metabolites. Produced by the action of myrosinase on GSLs, GSL-derived hydrolysis products (GHPs) primarily defend against biotic stress in planta. They also significantly affect the quality of crop products, with a subset of GHPs contributing unique food flavors and multiple therapeutic benefits or causing disagreeable food odors and health risks. Here, we explore the potential of these bioactive functions, which could be exploited for future sustainable agriculture. We first summarize our accumulated understanding of GSL diversity and distribution across representative Brassicaceae species. We then systematically discuss and evaluate the potential of exploited and unutilized genes involved in GSL biosynthesis, transport, and hydrolysis as candidate GSL engineering targets. Benefiting from available information on GSL and GHP functions, we explore options for multifunctional Brassicaceae crop ideotypes to meet future demand for food diversification and sustainable crop production. An integrated roadmap is subsequently proposed to guide ideotype development, in which maximization of beneficial effects and minimization of detrimental effects of GHPs could be combined and associated with various end uses. Based on several use-case examples, we discuss advantages and limitations of available biotechnological approaches that may contribute to effective deployment and could provide novel insights for optimization of future GSL engineering.

Identification and in vitro enzymatic activity analysis of the AOP2 gene family associated with glucosinolate biosynthesis in Tumorous stem mustard (Brassica juncea var. tumida)

The major enzyme encoded by the glucosinolate biosynthetic gene AOP2 is involved in catalyzing the conversion of glucoiberin (GIB) into sinigrin (SIN) in Brassicaceae crops. The AOP2 proteins have previously been identified in several Brassicaceae species, but not in Tumorous stem mustard. As per this research, the five identified members of the AOP2 family from the whole genome of Brassica juncea named BjuAOP2.1-BjuAOP2.5 were found to be evenly distributed on five chromosomes. The subcellular localization results implied that BjuAOP2 proteins were mainly concentrated in the cytoplasm. Phylogenetic analysis of the AOP2 proteins from the sequenced Brassicaceae species in BRAD showed that BjuAOP2 genes were more closely linked to Brassica carinata and Brassica rapa than Arabidopsis. In comparison with other Brassicaceae plants, the BjuAOP2 members were conserved in terms of gene structures, protein sequences, and motifs. The light response and hormone response elements were included in the BjuAOP2 genes' cis-regulatory elements. The expression pattern of BjuAOP2 genes was influenced by the different stages of development and the type of tissue being examined. The BjuAOP2 proteins were used to perform the heterologous expression experiment. The results showed that all the five BjuAOP2 proteins can catalyze the conversion of GIB to SIN with different catalytic activity. These results provide the basis for further investigation of the functional study of BjuAOP2 in Tumorous stem mustard glucosinolate biosynthesis.

Glucosinolates and Cytotoxic Activity of Collard Volatiles Obtained Using Microwave-Assisted Extraction

Glucosinolates (GSLs) in Brassica oleracea L. convar. acephala var. viridis (collard) flower, leaf, stem, and root were analyzed qualitatively and quantitatively via their desulfo-counterparts using UHPLC-DAD-MS/MS. Twelve GSLs were identified, including Met-derived GSLs (sinigrin, glucoibervirin, glucoerucin, glucoiberin, glucoraphanin, progoitrin), Trp-derived GSLs (4-hydroxyglucobrassicin, glucobrassicin, 4-methoxyglucobrassicin, and neoglucobrassicin), and Phe-derived GSLs (glucotropaeolin and gluconasturtiin). Total GSL content was highest in the root, having 63.40 μmol/g dried weight (DW), with gluconasturtiin (34.02 μmol/g DW) as the major GSL, followed by sinigrin and glucoibervirin (12.43 and 7.65 μmol/g DW, respectively). Total GSL contents in the flower, leaf, and stem were lower than in root, having 6.27, 2.64, and 1.84 μmol/g DW, respectively, with Trp and/or Met-derived GSLs as the predominant ones. GSL breakdown products were obtained via microwave hydrodiffusion and gravity (MHG) and volatile breakdown products were analyzed using GC-MS techniques. Volatile isolates were tested for their cytotoxic activity using MTT assay. MHG volatile extract from the root demonstrated the best cytotoxic activity against human bladder cancer cell line T24 and breast cancer cell line MDA-MB-231 during an incubation time of 72 h (IC₅₀ 21.58, and 11.62 μg/mL, respectively). The activity of the root extract can be attributed to its major volatile, 2-phenylethyl isothiocyanate (gluconasturtiin breakdown product).

The rhizosphere microbiome and host plant glucosinolates exhibit feedback cycles in Brassica rapa

The rhizosphere microbiome influences many aspects of plant fitness, including production of secondary compounds and defence against insect herbivores. Plants also modulate the composition of the microbial community in the rhizosphere via secretion of root exudates. We tested both the effect of the rhizosphere microbiome on plant traits, and host plant effects on rhizosphere microbes using recombinant inbred lines (RILs) of Brassica rapa that differ in production of glucosinolates (GLS), secondary metabolites that contribute to defence against insect herbivores. First, we investigated the effect of genetic variation in GLS production on the composition of the rhizosphere microbiome. Using a Bayesian Dirichlet-multinomial regression model (DMBVS), we identified both negative and positive associations between bacteria from six genera and the concentration of five GLS compounds produced in plant roots. Additionally, we tested the effects of microbial inoculation (an intact vs. disrupted soil microbiome) on GLS production and insect damage in these RILs. We found a significant microbial treatment × genotype interaction, in which total GLS was higher in the intact relative to the disrupted microbiome treatment in some RILs. However, despite differences in GLS production between microbial treatments, we observed no difference in insect damage between treatments. Together, these results provide evidence for a full feedback cycle of plant-microbe interactions mediated by GLS; that is, GLS compounds produced by the host plant "feed-down" to influence rhizosphere microbial community and rhizosphere microbes "feed-up" to influence GLS production.

Transcriptome and QTL mapping analyses of major QTL genes controlling glucosinolate contents in vegetable- and oilseed-type Brassica rapa plants

Glucosinolates (GSLs) are secondary metabolites providing defense against pathogens and herbivores in plants, and anti-carcinogenic activity against human cancer cells. Profiles of GSLs vary greatly among members of genus Brassica. In this study, we found that a reference line of Chinese cabbage (B. rapa ssp. pekinensis), 'Chiifu' contains significantly lower amounts of total GSLs than the oilseed-type B. rapa (B. rapa ssp. trilocularis) line 'LP08'. This study aimed to identify the key regulators of the high accumulation of GSLs in Brassica rapa plants using transcriptomic and linkage mapping approaches. Comparative transcriptome analysis showed that, in total, 8,276 and 9,878 genes were differentially expressed between 'Chiifu' and 'LP08' under light and dark conditions, respectively. Among 162 B. rapa GSL pathway genes, 79 were related to GSL metabolism under light conditions. We also performed QTL analysis using a single nucleotide polymorphism-based linkage map constructed using 151 F₅ individuals derived from a cross between the 'Chiifu' and 'LP08' inbred lines. Two major QTL peaks were successfully identified on chromosome 3 using high-performance liquid chromatography to obtain GSL profiles from 97 F₅ recombinant inbred lines. The MYB-domain transcription factor gene BrMYB28.1 (Bra012961) was found in the highest QTL peak region. The second highest peak was located near the 2-oxoacid-dependent dioxygenase gene BrGSL-OH.1 (Bra022920). This study identified major genes responsible for differing profiles of GSLs between 'Chiifu' and 'LP08'. Thus, our study provides molecular insights into differences in GSL profiles between vegetative- and oilseed-type B. rapa plants.

Genetic variation underlying differential ammonium and nitrate responses in Arabidopsis thaliana

Nitrogen is an essential element required for plant growth and productivity. Understanding the mechanisms and natural genetic variation underlying nitrogen use in plants will facilitate the engineering of plant nitrogen use to maximize crop productivity while minimizing environmental costs. To understand the scope of natural variation that may influence nitrogen use, we grew 1,135 Arabidopsis thaliana natural genotypes on two nitrogen sources, nitrate and ammonium, and measured both developmental and defense metabolite traits. By using different environments and focusing on multiple traits, we identified a wide array of different nitrogen responses. These responses are associated with numerous genes, most of which were not previously associated with nitrogen responses. Only a small portion of these genes appear to be shared between environments or traits, while most are predominantly specific to a developmental or defense trait under a specific nitrogen source. Finally, by using a large population, we were able to identify unique nitrogen responses, such as preferring ammonium or nitrate, which appear to be generated by combinations of loci rather than a few large-effect loci. This suggests that it may be possible to obtain novel phenotypes in complex nitrogen responses by manipulating sets of genes with small effects rather than solely focusing on large-effect single gene manipulations.

BocODD1 and BocODD2 Regulate the Biosynthesis of Progoitrin Glucosinolate in Chinese Kale

Progoitrin (2-hydroxy-3-butenyl glucosinolate, PRO) is the main source of bitterness of Brassica plants. Research on the biosynthesis of PRO glucosinolate can aid the understanding of the nutritional value in Brassica plants. In this study, four ODD genes likely involved in PRO biosynthesis were cloned from Chinese kale. These four genes, designated as BocODD1-4, shared 75-82% similarities with the ODD sequence of Arabidopsis. The sequences of these four BocODDs were analyzed, and BocODD1 and BocODD2 were chosen for further study. The gene BocODD1,2 showed the highest expression levels in the roots, followed by the leaves, flowers, and stems, which is in accordance with the trend of the PRO content in the same tissues. Both the expression levels of BocODD1,2 and the content of PRO were significantly induced by high- and low-temperature treatments. The function of BocODDs involved in PRO biosynthesis was identified. Compared with the wild type, the content of PRO was increased twofold in the over-expressing BocODD1 or BocODD2 plants. Meanwhile, the content of PRO was decreased in the BocODD1 or BocODD2 RNAi lines more than twofold compared to the wildtype plants. These results suggested that BocODD1 and BocODD2 may play important roles in the biosynthesis of PRO glucosinolate in Chinese kale.

Genome-wide association mapping within a local Arabidopsis thaliana population more fully reveals the genetic architecture for defensive metabolite diversity

A paradoxical finding from genome-wide association studies (GWAS) in plants is that variation in metabolite profiles typically maps to a small number of loci, despite the complexity of underlying biosynthetic pathways. This discrepancy may partially arise from limitations presented by geographically diverse mapping panels. Properties of metabolic pathways that impede GWAS by diluting the additive effect of a causal variant, such as allelic and genetic heterogeneity and epistasis, would be expected to increase in severity with the geographical range of the mapping panel. We hypothesized that a population from a single locality would reveal an expanded set of associated loci. We tested this in a French Arabidopsis thaliana population (less than 1 km transect) by profiling and conducting GWAS for glucosinolates, a suite of defensive metabolites that have been studied in depth through functional and genetic mapping approaches. For two distinct classes of glucosinolates, we discovered more associations at biosynthetic loci than the previous GWAS with continental-scale mapping panels. Candidate genes underlying novel associations were supported by concordance between their observed effects in the TOU-A population and previous functional genetic and biochemical characterization. Local populations complement geographically diverse mapping panels to reveal a more complete genetic architecture for metabolic traits. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

A keystone gene underlies the persistence of an experimental food web

Genes encode information that determines an organism's fitness. Yet we know little about whether genes of one species influence the persistence of interacting species in an ecological community. Here, we experimentally tested the effect of three plant defense genes on the persistence of an insect food web and found that a single allele at a single gene promoted coexistence by increasing plant growth rate, which in turn increased the intrinsic growth rates of species across multiple trophic levels. Our discovery of a "keystone gene" illustrates the need to bridge between biological scales, from genes to ecosystems, to understand community persistence.

The ease and complexity of identifying and using specialized metabolites for crop engineering

Plants produce a broad variety of specialized metabolites with distinct biological activities and potential applications. Despite this potential, most biosynthetic pathways governing specialized metabolite production remain largely unresolved across the plant kingdom. The rapid advancement of genetics and biochemical tools has enhanced our ability to identify plant specialized metabolic pathways. Further advancements in transgenic technology and synthetic biology approaches have extended this to a desire to design new pathways or move existing pathways into new systems to address long-running difficulties in crop systems. This includes improving abiotic and biotic stress resistance, boosting nutritional content, etc. In this review, we assess the potential and limitations for (1) identifying specialized metabolic pathways in plants with multi-omics tools and (2) using these enzymes in synthetic biology or crop engineering. The goal of these topics is to highlight areas of research that may need further investment to enhance the successful application of synthetic biology for exploiting the myriad of specialized metabolic pathways.

RNA Sequencing of Arabidopsis thaliana Seedlings after Non-Thermal Plasma-Seed Treatment Reveals Upregulation in Plant Stress and Defense Pathways

Not all agricultural practices are sustainable; however, non-thermal plasma treatment of seeds may be an eco-friendly alternative to improve macroscopic plant growth parameters. Despite the numerous successful results of plasma-seed treatments reported in the literature, there is a large gap in our understanding of how non-thermal plasma treatments affect seeds, especially due to the plethora of physical, chemical, and biological variables. This study uses RNA sequencing to characterize the changes in gene transcription in Arabidopsis thaliana (L.) Heynh. seeds 6 days after exposure to surface dielectric barrier discharge plasma treatment. Here, we provide an overview of all pathways that are differentially expressed where few genes are upregulated and many genes are downregulated. Our results reveal that plasma treatment time is a parameter that can activate different pathways in plant defense. An 80 s treatment upregulates the glucosinolate pathway, a defense response to insects and herbivores to deter feeding, whereas a shorter treatment of 60 s upregulates the phenylpropanoid pathway, which reinforces the cell wall with lignin and produces antimicrobial compounds, a defense response to bacterial or fungal plant pathogens. It seems that plasma elicits a wounding response from the seed in addition to redox changes. This suggests that plasma treatment can be potentially applied in agriculture to protect plants against abiotic and biotic stresses without discharging residues into the environment.

A 2-oxoglutarate-dependent dioxygenase converts dihydrofuran to furan in Salvia diterpenoids

Tanshinone ⅡA (TⅡA), a diterpene quinone with a furan ring, is a bioactive compound found in the medicinal herb redroot sage (Salvia miltiorrhiza Bunge), in which both furan and dihydrofuran analogs are present in abundance. Progress has been made recently in elucidating the tanshinone biosynthetic pathway, including heterocyclization of the dihydrofuran D-ring by cytochrome P450s; however, dehydrogenation of dihydrofuran to furan, a key step of furan ring formation, remains uncharacterized. Here, by differential transcriptome mining, we identified six 2-oxoglutarate-dependent dioxygenase (2-ODD) genes whose expressions corresponded to tanshinone biosynthesis. We showed that Sm2-ODD14 acts as a dehydrogenase catalyzing the furan ring aromatization. In vitro Sm2-ODD14 converted cryptotanshinone to TⅡA and thus was designated TⅡA synthase (SmTⅡAS). Furthermore, SmTⅡAS showed a strict substrate specificity, and repression of SmTⅡAS expression in hairy root by RNAi led to increased accumulation of total dihydrofuran-tanshinones and decreased production of furan-tanshinones. We conclude that SmTⅡAS controls the metabolite flux from dihydrofuran- to furan-tanshinones, which influences medicinal properties of S. miltiorrhiza.

Glucosinolate Profiles in Different Organs of 111 Radish Accessions and Candidate Genes Involved in Converting Glucobrassicin to 4-Hydroxyglucobrassicin

Glucosinolate (GSL) not only has highly physiological function for plants but also has considerable human interest. We analyzed the GSL compositions and levels in four organs of 111 radish accessions. Seven major GSLs were detected (approximately 5-245 μmol g^-1 DW), among which 4-(methylsulfinyl)but-3-enyl GSL and 4-methylsulfanyl-3-butenyl GSL were the dominant GSLs. GSL levels varied substantially among species and groups, and some genotypes/groups with special GSL profiles were identified. The total GSL level was higher in seeds than in sprouts, taproots, and leaves. Additionally, a correlation analysis revealed that seed 4-(methylsulfinyl)but-3-enyl GSL levels were highly correlated with sprout GSL levels. Moreover, a candidate gene (RsCYP81F2.3) encoding an enzyme that catalyzes the conversion of indol-3-ylmethyl GSL to 4-hydroxyindol-3-ylmethyl GSL was identified based on the detection and analysis of three radish accessions with relatively high indol-3-ylmethyl GSL, low 4-hydroxyindol-3-ylmethyl GSL, and 4-methoxyindol-3-ylmethyl GSL levels in their seeds. Our results provide some insights for finding materials and genes relevant for breeding new varieties with ideal GSL compositions and levels.

Tomato E8 Encodes a C-27 Hydroxylase in Metabolic Detoxification of α-Tomatine during Fruit Ripening

Tomato (Solanum lycopersicum) contains α-tomatine, a steroidal glycoalkaloid that contributes to the plant defense against pathogens and herbivores through its bitter taste and toxicity. It accumulates at high levels in all the plant tissues, especially in leaves and immature green fruits, whereas it decreases during fruit ripening through metabolic conversion to the nontoxic esculeoside A, which accumulates in the mature red fruit. This study aimed to identify the gene encoding a C-27 hydroxylase that is a key enzyme in the metabolic conversion of α-tomatine to esculeoside A. The E8 gene, encoding a 2-oxoglutalate-dependent dioxygenase, is well known as an inducible gene in response to ethylene during fruit ripening. The recombinant E8 was found to catalyze the C-27 hydroxylation of lycoperoside C to produce prosapogenin A and is designated as Sl27DOX. The ripe fruit of E8/Sl27DOX-silenced transgenic tomato plants accumulated lycoperoside C and exhibited decreased esculeoside A levels compared with the wild-type (WT) plants. Furthermore, E8/Sl27DOX deletion in tomato accessions resulted in higher lycoperoside C levels in ripe fruits than in WT plants. Thus, E8/Sl27DOX functions as a C-27 hydroxylase of lycoperoside C in the metabolic detoxification of α-tomatine during tomato fruit ripening, and the efficient detoxification by E8/27DOX may provide an advantage in the domestication of cultivated tomatoes.

Genome-Wide Identification and Expression Profiling of 2OGD Superfamily Genes from Three Brassica Plants

The 2-oxoglutarate and Fe(II)-dependent dioxygenase (2OGD) superfamily is the second largest enzyme family in the plant genome, and its members are involved in various oxygenation and hydroxylation reactions. Due to their important biochemical significance in metabolism, a systematic analysis of the plant 2OGD genes family is necessary. Here, we identified 160, 179, and 337 putative 2OGDs from Brassica rapa, Brassica oleracea, and Brassica napus. According to their gene structure, domain, phylogenetic features, function, and previous studies, we also divided 676 2OGDs into three subfamilies: DOXA, DOXB, and DOXC. Additionally, homologous and phylogenetic comparisons of three subfamily genes provided valuable insight into the evolutionary characteristics of the 2OGD genes from Brassica plants. Expression profiles derived from the transcriptome and Genevestigator database exhibited distinct expression patterns of the At2OGD, Br2OGD, and Bo2OGD genes in different developmental stages, tissues, or anatomical parts. Some 2OGD genes showed high expression levels in various tissues, such as callus, seed, silique, and root tissues, while other 2OGD genes were expressed at very low levels in other tissues. Analysis of six Bo2OGD genes in different tissues by qRT-PCR indicated that these genes are involved in the metabolism of gibberellin, which in turn regulates plant growth and development. Our working system analysed 2OGD gene families of three Brassica plants and laid the foundation for further study of their functional characterization.

Effects of the Developmental Regulator BOLITA on the Plant Metabolome

Transcription factors are important regulators of gene expression. They can orchestrate the activation or repression of hundreds or thousands of genes and control diverse processes in a coordinated way. This work explores the effect of a master regulator of plant development, BOLITA (BOL), in plant metabolism, with a special focus on specialized metabolism. For this, we used an Arabidopsis thaliana line in which the transcription factor activity can be induced. Fingerprinting metabolomic analyses of whole plantlets were performed at different times after induction. After 96 h, all induced replicas clustered as a single group, in contrast with all controls which did not cluster. Metabolomic analyses of shoot and root tissues enabled the putative identification of differentially accumulated metabolites in each tissue. Finally, the analysis of global gene expression in induced vs. non-induced root samples, together with enrichment analyses, allowed the identification of enriched metabolic pathways among the differentially expressed genes and accumulated metabolites after the induction. We concluded that the induction of BOL activity can modify the Arabidopsis metabolome. Future work should investigate whether its action is direct or indirect, and the implications of the metabolic changes for development regulation and bioprospection.

Genetic variation, environment and demography intersect to shape Arabidopsis defense metabolite variation across Europe

Plants produce diverse metabolites to cope with the challenges presented by complex and ever-changing environments. These challenges drive the diversification of specialized metabolites within and between plant species. However, we are just beginning to understand how frequently new alleles arise controlling specialized metabolite diversity and how the geographic distribution of these alleles may be structured by ecological and demographic pressures. Here, we measure the variation in specialized metabolites across a population of 797 natural Arabidopsis thaliana accessions. We show that a combination of geography, environmental parameters, demography and different genetic processes all combine to influence the specific chemotypes and their distribution. This showed that causal loci in specialized metabolism contain frequent independently generated alleles with patterns suggesting potential within-species convergence. This provides a new perspective about the complexity of the selective forces and mechanisms that shape the generation and distribution of allelic variation that may influence local adaptation.

Since plants cannot move, they have evolved chemical defenses to help them respond to changes in their surroundings. For example, where animals run from predators, plants may produce toxins to put predators off. This approach is why plants are such a rich source of drugs, poisons, dyes and other useful substances. The chemicals plants produce are known as specialized metabolites, and they can change a lot between, and even within, plant species. The variety of specialized metabolites is a result of genetic changes and evolution over millions of years.

Evolution is a slow process, yet plants are able to rapidly develop new specialized metabolites to protect them from new threats. Even different populations of the same species produce many distinct metabolites that help them survive in their surroundings. However, the factors that lead plants to produce new metabolites are not well understood, and it is not known how this affects genetic variation.

To gain a better understanding of this process, Katz et al. studied 797 European variants of a common weed species called Arabidopsis thaliana, which is widely studied. The investigation found that many factors affect the range of specialized metabolites in each variant. These included local geography and environment, as well as genetics and population history (demography). Katz et al. revealed a pattern of relationships between the variants that could mirror their evolutionary history as the species spread and adapted to new locations.

These results highlight the complex network of factors that affect plant evolution. Rapid diversification is key to plant survival in new and changing environments and has resulted in a wide range of specialized metabolites. As such they are of interest both for studying plant evolution and for understanding their ecology. Expanding similar work to more populations and other species will broaden the scope of our ability to understand how plants adapt to their surroundings.

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.

Abscisic acid-mediated induction of FLAVIN-CONTAINING MONOOXYGENASE 2 leads to reduced accumulation of methylthioalkyl glucosinolates in Arabidopsis thaliana

Side-chain modification contributes to the structural diversity of aliphatic glucosinolates (GSLs), a class of sulfur-containing secondary metabolites found in Brassicales. The first step in side-chain modification of aliphatic GSLs is the S-oxygenation of the methylthioalkyl (MT) moiety to the methylsulfinylalkyl (MS) moiety. This reaction is catalyzed by flavin-containing monooxygenase (FMO_GS-OX), which is encoded by seven genes in Arabidopsis thaliana. Therefore, the regulation of FMO_GS-OX gene expression is key to controlling side-chain structural diversity. In this study, we demonstrated that the expression of FMO_GS-OX2 and FMO_GS-OX4 was induced by glucose treatment, independent of MYB28/29 and MYC2/3/4, the transcription factors that positively regulate aliphatic GSL biosynthesis. Glucose treatment of the abscisic acid (ABA)-related mutants indicated that glucose-triggered upregulation of FMO_GS-OX2 and FMO_GS-OX4 was partially regulated by ABA through the key negative regulators ABI1 and ABI2, and the positive regulator SnRK2, but not via the transcription factor ABI5. In wild-type plants, glucose treatment drastically reduced the accumulation of 4-methylthiobutyl (4MT) GSL, whereas a decrease in 4MT GSL was not observed in the fmo_gs-ox2, abi1-1, abi2-1, aba2-1, or aba3-1 mutants. This result indicated that the decreased accumulation of 4MT GSL by glucose treatment was attributed to upregulation of FMO_GS-OX2 via the ABA signaling pathway.

Glucosinolate Biosynthesis and the Glucosinolate–Myrosinase System in Plant Defense

Insect pests represent a major global challenge to important agricultural crops. Insecticides are often applied to combat such pests, but their use has caused additional challenges such as environmental contamination and human health issues. Over millions of years, plants have evolved natural defense mechanisms to overcome insect pests and pathogens. One such mechanism is the production of natural repellents or specialized metabolites like glucosinolates. There are three types of glucosinolates produced in the order Brassicales: aliphatic, indole, and benzenic glucosinolates. Upon insect herbivory, a “mustard oil bomb” consisting of glucosinolates and their hydrolyzing enzymes (myrosinases) is triggered to release toxic degradation products that act as insect deterrents. This review aims to provide a comprehensive summary of glucosinolate biosynthesis, the “mustard oil bomb”, and how these metabolites function in plant defense against pathogens and insects. Understanding these defense mechanisms will not only allow us to harness the benefits of this group of natural metabolites for enhancing pest control in Brassicales crops but also to transfer the “mustard oil bomb” to non-glucosinolate producing crops to boost their defense and thereby reduce the use of chemical pesticides.

Greenbug (Schizaphis graminum) herbivory significantly impacts protein and phosphorylation abundance in switchgrass (Panicum virgatum)

Switchgrass (Panicum virgatum L.) is an important crop for biofuel production but it also serves as host for greenbugs (Schizaphis graminum Rondani; GB). Although transcriptomic studies have been done to infer the molecular mechanisms of plant defense against GB, little is known about the effect of GB infestation on the switchgrass protein expression and phosphorylation regulation. The global response of the switchgrass cultivar Summer proteome and phosphoproteome was monitored by label-free proteomics shotgun in GB-infested and uninfested control plants at 10 days post infestation. Peptides matching a total of 3,594 proteins were identified and 429 were differentially expressed proteins in GB-infested plants relative to uninfested control plants. Among these, 291 and 138 were up and downregulated by GB infestation, respectively. Phosphoproteome analysis identified 310 differentially phosphorylated proteins (DP) from 350 phosphopeptides with a total of 399 phosphorylated sites. These phosphopeptides had more serine phosphorylated residues (79%), compared to threonine phosphorylated sites (21%). Overall, KEGG pathway analysis revealed that GB feeding led to the enriched accumulation of proteins important for biosynthesis of plant defense secondary metabolites and repressed the accumulation of proteins involved in photosynthesis. Interestingly, defense modulators such as terpene synthase, papain-like cysteine protease, serine carboxypeptidase, and lipoxygenase2 were upregulated at the proteome level, corroborating previously published transcriptomic data.

A Comprehensive Gene Inventory for Glucosinolate Biosynthetic Pathway in Arabidopsis thaliana

Glucosinolates (GSLs) are plant secondary metabolites comprising sulfur and nitrogen mainly found in plants from the order of Brassicales, such as broccoli, cabbage, and Arabidopsis thaliana. The activated forms of GSL play important roles in fighting against pathogens and have health benefits to humans. The increasing amount of data on A. thaliana generated from various omics technologies can be investigated more deeply in search of new genes or compounds involved in GSL biosynthesis and metabolism. This review describes a comprehensive inventory of A. thaliana GSLs identified from published literature and databases such as KNApSAcK, KEGG, and AraCyc. A total of 113 GSL genes encoding for 23 transcription components, 85 enzymes, and five protein transporters were experimentally characterized in the past two decades. Continuous efforts are still on going to identify all molecules related to the production of GSLs. A manually curated database known as SuCCombase (http://plant-scc.org) was developed to serve as a comprehensive GSL inventory. Realizing lack of information on the regulation of GSL biosynthesis and degradation mechanisms, this review also includes relevant information and their connections with crosstalk among various factors, such as light, sulfur metabolism, and nitrogen metabolism, not only in A. thaliana but also in other crucifers.

Dissection of genetic architecture for glucosinolate accumulations in leaves and seeds of Brassica napus by genome-wide association study

Glucosinolates (GSLs), whose degradation products have been shown to be increasingly important for human health and plant defence, compose important secondary metabolites found in the order Brassicales. It is highly desired to enhance pest and disease resistance by increasing the leaf GSL content while keeping the content low in seeds of Brassica napus, one of the most important oil crops worldwide. Little is known about the regulation of GSL accumulation in the leaves. We quantified the levels of 9 different GSLs and 15 related traits in the leaves of 366 accessions and found that the seed and leaf GSL content were highly correlated (r = 0.79). A total of 78 loci were associated with GSL traits, and five common and eleven tissue-specific associated loci were related to total leaf and seed GSL content. Thirty-six candidate genes were inferred to be involved in GSL biosynthesis. The candidate gene BnaA03g40190D (BnaA3.MYB28) was validated by DNA polymorphisms and gene expression analysis. This gene was responsible for high leaf/low seed GSL content and could explain 30.62% of the total leaf GSL variation in the low seed GSL panel and was not fixed during double-low rapeseed breeding. Our results provide new insights into the genetic basis of GSL variation in leaves and seeds and may facilitate the metabolic engineering of GSLs and the breeding of high leaf/low seed GSL content in B. napus.

Exploring glucosinolates diversity in Brassicaceae: a genomic and chemical assessment for deciphering abiotic stress tolerance

Brassica is one of the most economically important genus of the Brassicaceae family, encompassing several key crops like Brassica napus (cabbage) and broccoli (Brassica oleraceae var. italica). This family is well known for their high content of characteristic secondary metabolites such as glucosinolates (GLS) compounds, recognize for their beneficial health properties and role in plants defense. In this work, we have looked through gene clusters involved in the biosynthesis of GLS, by combining genomic analysis with biochemical pathways and chemical diversity assessment. A total of 101 Brassicaceae genes involved in GLS biosynthesis were identified, using a multi-database approach. Through a UPGMA and PCA analysis on the 101 GLS genes recorded, revealed a separation between the genes mainly involved in GLS core structure synthesis and genes belonging to the CYP450s and MYBs gene families. After, a detailed phylogenetic analysis was conducted to better understand the disjunction of the aliphatic and indolic genes, by focusing on CYP79F1-F2 and CYP81F1-F4, respectively. Our results point to a recent diversification of the aliphatic CYP79F1 and F2 genes in Brassica crops, while for indolic genes an earliest diversification is observed for CYP81F1-F4 genes. Chemical diversity revealed that Brassica crops have distinct GLS chemo-profiles from other Brassicaceae genera; being highlighted the high contents of GLS found among the Diplotaxis species. Also, we have explored GLS-rich species as a new source of taxa with great agronomic potential, particularly in abiotic stress tolerance, namely Diplotaxis, the closest wild relatives of Brassica crops.

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.

Identification of α-Tomatine 23-Hydroxylase Involved in the Detoxification of a Bitter Glycoalkaloid

Tomato plants (Solanum lycopersicum) contain steroidal glycoalkaloid α-tomatine, which functions as a chemical barrier to pathogens and predators. α-Tomatine accumulates in all tissues and at particularly high levels in leaves and immature green fruits. The compound is toxic and causes a bitter taste, but its presence decreases through metabolic conversion to nontoxic esculeoside A during fruit ripening. This study identifies the gene encoding a 23-hydroxylase of α-tomatine, which is a key to this process. Some 2-oxoglutarate-dependent dioxygenases were selected as candidates for the metabolic enzyme, and Solyc02g062460, designated Sl23DOX, was found to encode α-tomatine 23-hydroxylase. Biochemical analysis of the recombinant Sl23DOX protein demonstrated that it catalyzes the 23-hydroxylation of α-tomatine and the product spontaneously isomerizes to neorickiioside B, which is an intermediate in α-tomatine metabolism that appears during ripening. Leaves of transgenic tomato plants overexpressing Sl23DOX accumulated not only neorickiioside B but also another intermediate, lycoperoside C (23-O-acetylated neorickiioside B). Furthermore, the ripe fruits of Sl23DOX-silenced transgenic tomato plants contained lower levels of esculeoside A but substantially accumulated α-tomatine. Thus, Sl23DOX functions as α-tomatine 23-hydroxylase during the metabolic processing of toxic α-tomatine in tomato fruit ripening and is a key enzyme in the domestication of cultivated tomatoes.

Plant Networks as Traits and Hypotheses: Moving Beyond Description

Biology relies on the central thesis that the genes in an organism encode molecular mechanisms that combine with stimuli and raw materials from the environment to create a final phenotypic expression representative of the genomic programming. While conceptually simple, the genotype-to-phenotype linkage in a eukaryotic organism relies on the interactions of thousands of genes and an environment with a potentially unknowable level of complexity. Modern biology has moved to the use of networks in systems biology to try to simplify this complexity to decode how an organism's genome works. Previously, biological networks were basic ways to organize, simplify, and analyze data. However, recent advances are allowing networks to move beyond description and become phenotypes or hypotheses in their own right. This review discusses these efforts, like mapping responses across biological scales, including relationships among cellular entities, and the direct use of networks as traits or hypotheses.

Glucosinolate variability between turnip organs during development

Turnip (Brassica rapa spp. rapa) is an important vegetable species, with a unique physiology. Several plant parts, including both the turnip tubers and leaves, are important for human consumption. During the development of turnip plants, the leaves function as metabolic source tissues, while the tuber first functions as a sink, while later the tuber turns into a source for development of flowers and seeds. In the present study, chemical changes were determined for two genotypes with different genetic background, and included seedling, young leaves, mature leaves, tuber surface, tuber core, stalk, flower and seed tissues, at seven different time points during plant development. As a basis for understanding changes in glucosinolates during plant development, the profile of glucosinolates was analysed using liquid chromatography (LC) coupled to mass spectrometry (MS). This analysis was complemented by a gene expression analysis, focussed on GLS biosynthesis, which could explain part of the observed variation, pointing to important roles of specific gene orthologues for defining the chemical differences. Substantial differences in glucosinolate profiles were observed between above-ground tissues and turnip tuber, reflecting the differences in physiological role. In addition, differences between the two genotypes and between tissues that were harvested early or late during the plant lifecycle. The importance of the observed differences in glucosinolate profile for the ecophysiology of the turnip and for breeding turnips with optimal chemical profiles is discussed.

A high-density genetic map and QTL mapping of leaf traits and glucosinolates in Barbarea vulgaris

Background

Barbarea vulgaris is a wild cruciferous plant and include two distinct types: the G- and P-types named after their glabrous and pubescent leaves, respectively. The types differ significantly in resistance to a range of insects and diseases as well as glucosinolates and other chemical defenses. A high-density linkage map was needed for further progress to be made in the molecular research of this plant.

Results

We performed restriction site-associated DNA sequencing (RAD-Seq) on an F₂ population generated from G- and P-type B. vulgaris. A total of 1545 SNP markers were mapped and ordered in eight linkage groups, which represents the highest density linkage map to date for the crucifer tribe Cardamineae. A total of 722 previously published genome contigs (50.2 Mb, 30% of the total length) can be anchored to this high density genetic map, an improvement compared to a previously published map (431 anchored contigs, 38.7 Mb, 23% of the assembly genome). Most of these (572 contigs, 31.2 Mb) were newly anchored to the map, representing a significant improvement. On the basis of the present high-density genetic map, 37 QTL were detected for eleven traits, each QTL explaining 2.9–71.3% of the phenotype variation. QTL of glucosinolates, leaf size and color traits were in most cases overlapping, possibly implying a functional connection.

Conclusions

This high-density linkage map and the QTL obtained in this study will be useful for further understanding of the genetic of the B. vulgaris and molecular basis of these traits, many of which are shared in the related crop watercress.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5769-z) contains supplementary material, which is available to authorized users.

Transcriptome reveals the gene expression patterns of sulforaphane metabolism in broccoli florets

Sulforaphane is a new and effective anti-cancer component that is abundant in broccoli. In the past few years, the patterns of variability in glucosinolate content and its regulation in A. thaliana have been described in detail. However, the diversity of glucosinolate and sulforaphane contents in different organs during vegetative and reproductive stages has not been clearly explained. In this paper, we firstly investigated the transcriptome profiles of the developing buds and leaves at bolting stage of broccoli (B52) to further assess the gene expression patterns involved in sulforaphane synthesis. The CYP79F1 gene, as well as nine other genes related to glucorahpanin biosynthesis, MAM1, MAM3, St5b-2, FMO GS-OX1, MY, AOP2, AOP3, ESP and ESM1 were selected by digital gene expression analysis and were validated by quantitative real-time PCR (qRT-PCR). Meanwhile, the compositions of glucosinolates and sulforaphane were detected for correlation analysis with related genes. Finally the RNA sequencing libraries generated 147 957 344 clean reads, and 8 539 unigene assemblies were produced. In digital result, only CYP79F1, in the glucoraphanin pathway, was up-regulated in young buds but absent from the other organs, which was consistent with the highest level of sulforaphane content being in this organ compared to mature buds, buds one day before flowering, flowers and leaves. The sequencing results also presented that auxin and cytokinin might affect glucoraphanin accumulation. The study revealed that up-regulated expression of CYP79F1 plays a fundamental and direct role in sulforaphane production in inflorescences. Two genes of MAM1 and St5b-2 could up-regulated glucoraphanin generation. Synergistic expression of MAM1, MAM3, St5b-2, FMO GS-OX1, MY, ESP and ESM1 was found in sulforaphane metabolism. This study will be beneficial for understanding the diversity of sulforaphane in broccoli organs.

The Role of the Glucosinolate-Myrosinase System in Mediating Greater Resistance of Barbarea verna than B. vulgaris to Mamestra brassicae Larvae

We investigated the influences of two structurally similar glucosinolates, phenethylglucosinolate (gluconasturtiin, NAS) and its (S)-2-hydroxyl derivative glucobarbarin (BAR), as well as their hydrolysis products on larvae of the generalist Mamestra brassicae (Lepidoptera: Noctuidae). Previous results suggested a higher defensive activity of BAR than NAS based on resistance toward M. brassicae larvae of natural plant genotypes of Barbarea vulgaris R. Br. (Brassicaceae) dominated by BAR. In the present study, the hypothesis of a higher defensive activity of BAR than NAS was tested by comparing two Barbarea species similarly dominated either by BAR or by NAS and by testing effects of isolated BAR and NAS on larval survival and feeding preferences. Larvae reared on leaf disks of B. verna (Mill.) Asch. had a lower survival than those reared on B. vulgaris P- and G-chemotypes. Leaves of B. verna were dominated by NAS, whereas B. vulgaris chemotypes were dominated by BAR or its epimer. In addition, B. verna leaves showed a threefold higher activity of the glucosinolate-activating myrosinase enzymes. The main product of NAS from breakdown by endogenous enzymes including myrosinases ("autolysis") in B. verna leaves was phenethyl isothiocyanate, while the main products of BAR in autolyzed B. vulgaris leaves were a cyclized isothiocyanate product, namely an oxazolidine-2-thione, and a downstream metabolite, an oxazolidin-2-one. The glucosinolates BAR and NAS were isolated and offered to larvae on disks of cabbage. Both glucosinolates exerted similar negative effects on larval survival but effects of NAS tended to be more detrimental. Low concentrations of BAR, but not of NAS, stimulated larval feeding, whereas high BAR concentrations acted deterrent. NAS only tended to be deterrent at the highest concentration, but the difference was not significant. Recoveries of NAS and BAR on cabbage leaf disks were similar, and when hydrolyzed by mechanical leaf damage, the same isothiocyanate-type products as in Barbarea plants were formed with further conversion of BAR to cyclic products, (R)-5-phenyloxazolidine-2-thione [(R)-barbarin] and (R)-5-phenyloxazolidin-2-one [(R)-resedine]. We conclude that a previously proposed generally higher defensive activity of BAR than NAS to M. brassicae larvae could not be confirmed. Indeed, the higher resistance of NAS-containing B. verna plants may be due to a combined effect of rather high concentrations of NAS and a relatively high myrosinase activity or other plant traits not investigated yet.

Glucosinolate turnover in Brassicales species to an oxazolidin-2-one, formed via the 2-thione and without formation of thioamide

Glucosinolates are found in plants of the order Brassicales and hydrolyzed to different breakdown products, particularly after tissue damage. In Barbarea vulgaris R.Br. (Brassicaceae), the dominant glucosinolate in the investigated "G-type" is glucobarbarin, (S)-2-hydroxy-2-phenylethylglucosinolate. Formation of the nitrile from glucobarbarin was observed in vitro, while a previously suggested thioamide (synonym thionamide) was not confirmed. Resedine (5-phenyl-1,3-oxazolidin-2-one) was detected after glucobarbarin hydrolysis in crushed B. vulgaris leaves and siliques, but not in intact parts. The abundance increased for several hours after completion of hydrolysis. The corresponding 1,3-oxazolidine-2-thione (OAT), with the common name barbarin, was also formed, and appeared to be the precursor of resedine. Addition of each of two non-endogenous OATs, (S)-5-ethyl-5-methylOAT and (R)-5-vinylOAT (R-goitrin), to a leaf homogenate resulted in formation of the corresponding 1,3-oxazolidin-2-ones (OAOs), confirming the metabolic connection of OAT to OAO. Formation of OAOs was inhibited by prior brief heating of the homogenate, suggesting enzyme involvement. We suggest the conversion of OATs to OAOs to be catalyzed by an enzyme ("oxazolidinethionase") responsible for turnover of OAT formed in intact plants. Resedine had been reported as an alkaloid from another species - Reseda luteola L. (Resedaceae) - naturally containing the glucosinolate glucobarbarin. However, resedine was not detected in intact R. luteola plants, but formed after tissue damage. The formation of resedine in two families suggests a broad distribution of putative OATases in the Brassicales; potentially involved in glucosinolate turnover that needs myrosinase activity as the committed step. In agreement with the proposed function of OATase, several candidate genes for myrosinases in glucosinolate turnover in intact plants were discovered in the B. vulgaris genome. We also suggest that biotechnological conversion of OATs to OAOs might improve the nutritional value of Brassicales protein. HPLC-MS/MS methods for detection of these glucobarbarin products are described.

Comparison of the Relative Potential for Epigenetic and Genetic Variation To Contribute to Trait Stability

The theoretical ability of epigenetic variation to influence the heritable variation of complex traits is gaining traction in the study of adaptation. This theory posits that epigenetic marks can control adaptive phenotypes but the relative potential of epigenetic variation in comparison to genetic variation in these traits is not presently understood. To compare the potential of epigenetic and genetic variation in adaptive traits, we analyzed the influence of DNA methylation variation on the accumulation of chemical defense compounds glucosinolates from the order Brassicales. Several decades of work on glucosinolates has generated extensive knowledge about their synthesis, regulation, genetic variation and contribution to fitness establishing this pathway as a model pathway for complex adaptive traits. Using high-throughput phenotyping with a randomized block design of ddm1 derived Arabidopsis thaliana epigenetic Recombinant Inbred Lines, we measured the correlation between DNA methylation variation and mean glucosinolate variation and within line stochastic variation. Using this information, we identified epigenetic Quantitative Trait Loci that contained specific Differentially Methylated Regions associated with glucosinolate traits. This showed that variation in DNA methylation correlates both with levels and variance of glucosinolates and flowering time with trait-specific loci. By conducting a meta-analysis comparing the results to different genetically variable populations, we conclude that the influence of DNA methylation variation on these adaptive traits is much lower than the corresponding impact of standing genetic variation. As such, selective pressure on these traits should mainly affect standing genetic variation to lead to adaptation.

Mapping the Arabidopsis Metabolic Landscape by Untargeted Metabolomics at Different Environmental Conditions

Metabolic genome-wide association studies (mGWAS), whereupon metabolite levels are regarded as traits, can help unravel the genetic basis of metabolic networks. A total of 309 Arabidopsis accessions were grown under two independent environmental conditions (control and stress) and subjected to untargeted LC-MS-based metabolomic profiling; levels of the obtained hydrophilic metabolites were used in GWAS. Our two-condition-based GWAS for more than 3000 semi-polar metabolites resulted in the detection of 123 highly resolved metabolite quantitative trait loci (p ≤ 1.0E-08), 24.39% of which were environment-specific. Interestingly, differently from natural variation in Arabidopsis primary metabolites, which tends to be controlled by a large number of small-effect loci, we found several major large-effect loci alongside a vast number of small-effect loci controlling variation of secondary metabolites. The two-condition-based GWAS was followed by integration with network-derived metabolite-transcript correlations using a time-course stress experiment. Through this integrative approach, we selected 70 key candidate associations between structural genes and metabolites, and experimentally validated eight novel associations, two of them showing differential genetic regulation in the two environments studied. We demonstrate the power of combining large-scale untargeted metabolomics-based GWAS with time-course-derived networks both performed under different abiotic environments for identifying metabolite-gene associations, providing novel global insights into the metabolic landscape of Arabidopsis.

Network-Guided Discovery of Extensive Epistasis between Transcription Factors Involved in Aliphatic Glucosinolate Biosynthesis

Plants use diverse mechanisms influenced by vast regulatory networks of indefinite scale to adapt to their environment. These regulatory networks have an unknown potential for epistasis between genes within and across networks. To test for epistasis within an adaptive trait genetic network, we generated and tested 47 Arabidopsis thaliana double mutant combinations for 20 transcription factors, which all influence the accumulation of aliphatic glucosinolates, the defense metabolites that control fitness. The epistatic combinations were used to test if there is more or less epistasis depending on gene membership within the same or different phenotypic subnetworks. Extensive epistasis was observed between the transcription factors, regardless of subnetwork membership. Metabolite accumulation displayed antagonistic epistasis, suggesting the presence of a buffering mechanism. Epistasis affecting enzymatic estimated activity was highly conditional on the tissue and environment and shifted between both antagonistic and synergistic forms. Transcriptional analysis showed that epistasis shifts depend on how the trait is measured. Because the 47 combinations described here represent a small sampling of the potential epistatic combinations in this genetic network, there is potential for significantly more epistasis. Additionally, the main effect of the individual gene was not predictive of the epistatic effects, suggesting that there is a need for further studies.