A complex interplay of three R2R3 MYB transcription factors determines the profile of aliphatic glucosinolates in Arabidopsis

Cell-type-aware regulatory landscapes governing monoterpene indole alkaloid biosynthesis in the medicinal plant Catharanthus roseus

In plants, the biosynthetic pathways of some specialized metabolites are partitioned into specialized or rare cell types, as exemplified by the monoterpenoid indole alkaloid (MIA) pathway of Catharanthus roseus (Madagascar Periwinkle), the source of the anticancer compounds vinblastine and vincristine. In the leaf, the C. roseus MIA biosynthetic pathway is partitioned into three cell types with the final known steps of the pathway expressed in the rare cell type termed idioblast. How cell-type specificity of MIA biosynthesis is achieved is poorly understood. We generated single-cell multi-omics data from C. roseus leaves. Integrating gene expression and chromatin accessibility profiles across single cells, as well as transcription factor (TF)-binding site profiles, we constructed a cell-type-aware gene regulatory network for MIA biosynthesis. We showcased cell-type-specific TFs as well as cell-type-specific cis-regulatory elements. Using motif enrichment analysis, co-expression across cell types, and functional validation approaches, we discovered a novel idioblast-specific TF (Idioblast MYB1, CrIDM1) that activates expression of late-stage MIA biosynthetic genes in the idioblast. These analyses not only led to the discovery of the first documented cell-type-specific TF that regulates the expression of two idioblast-specific biosynthetic genes within an idioblast metabolic regulon but also provides insights into cell-type-specific metabolic regulation.

Leaf abaxial and adaxial surfaces differentially affect the interaction of Botrytis cinerea across several eudicots

Eudicot plant species have leaves with two surfaces: the lower abaxial and the upper adaxial surface. Each surface varies in a diversity of components and molecular signals, resulting in potentially different degrees of resistance to pathogens. We tested how Botrytis cinerea, a necrotroph fungal pathogen, interacts with the two different leaf surfaces across 16 crop species and 20 Arabidopsis genotypes. This showed that the abaxial surface is generally more susceptible to the pathogen than the adaxial surface. In Arabidopsis, the differential lesion area between leaf surfaces was associated with jasmonic acid (JA) and salicylic acid (SA) signaling and differential induction of defense chemistry across the two surfaces. When infecting the adaxial surface, leaves mounted stronger defenses by producing more glucosinolates and camalexin defense compounds, partially explaining the differential susceptibility across surfaces. Testing a collection of 96 B. cinerea strains showed the genetic heterogeneity of growth patterns, with a few strains preferring the adaxial surface while most are more virulent on the abaxial surface. Overall, we show that leaf-Botrytis interactions are complex with host-specific, surface-specific, and strain-specific patterns.

NAD+ deficiency primes defense metabolism via 1O2-escalated jasmonate biosynthesis in plants

Nicotinamide adenine dinucleotide (NAD+) is a redox cofactor and signal central to cell metabolisms. Disrupting NAD homeostasis in plant alters growth and stress resistance, yet the underlying mechanisms remain largely unknown. Here, by combining genetics with multi-omics, we discover that NAD+ deficiency in qs-2 caused by mutation in NAD+ biosynthesis gene-Quinolinate Synthase retards growth but induces biosynthesis of defense compounds, notably aliphatic glucosinolates that confer insect resistance. The elevated defense in qs-2 is resulted from activated jasmonate biosynthesis, critically hydroperoxidation of α-linolenic acid by the 13-lipoxygenase (namely LOX2), which is escalated via the burst of chloroplastic ROS-singlet oxygen (1O2). The NAD+ deficiency-mediated JA induction and defense priming sequence in plants is recapitulated upon insect infestation, suggesting such defense mechanism operates in plant stress response. Hence, NAD homeostasis is a pivotal metabolic checkpoint that may be manipulated to navigate plant growth and defense metabolism for stress acclimation.

Gene regulatory networks underlying sulfate deficiency responses in plants

Sulfur (S) is an essential macronutrient for plants and its availability in soils is an important determinant for growth and development. Current regulatory policies aimed at reducing industrial S emissions together with changes in agronomical practices have led to a decline in S contents in soils worldwide. Deficiency of sulfate-the primary form of S accessible to plants in soil-has adverse effects on both crop yield and nutritional quality. Hence, recent research has increasingly focused on unraveling the molecular mechanisms through which plants detect and adapt to a limiting supply of sulfate. A significant part of these studies involves the use of omics technologies and has generated comprehensive catalogs of sulfate deficiency-responsive genes and processes, principally in Arabidopsis together with a few studies centering on crop species such as wheat, rice, or members of the Brassica genus. Although we know that sulfate deficiency elicits an important reprogramming of the transcriptome, the transcriptional regulators orchestrating this response are not yet well understood. In this review, we summarize our current knowledge of gene expression responses to sulfate deficiency and recent efforts towards the identification of the transcription factors that are involved in controlling these responses. We further compare the transcriptional response and putative regulators between Arabidopsis and two important crop species, rice and tomato, to gain insights into common mechanisms of the response to sulfate deficiency.

The Antimicrobial Effects of Myrosinase Hydrolysis Products Derived from Glucosinolates Isolated from Lepidium draba

Lepidium draba (hoary cress) is a perennial plant belonging to the Brassicaceae family that produces two dominant glucosinolates (GLSs): glucoraphanin (GRN) and sinalbin (SBN). They represent the stored form, which is converted upon the myrosinase (Myr) hydrolysis activity to active compounds, mainly isothiocyanates (ITCs) such as sulforaphane (SFN) or p-hydroxybenzyl isothiocyanate (pHBITC). Research on ITCs that have proven anticancer, antimicrobial, and chemoprotective properties is usually conducted with pure commercially available compounds. However, these are chemically reactive, making it difficult to use them directly for preventive purposes in dietary supplements. Efforts are currently being made to prepare dietary supplements enriched with GLS and/or Myr. In this study, we report a simple but efficient chromatographic procedure for the isolation and purification of GLSs from MeOH extract from hoary cress based on a combination of ion exchange and gel permeation chromatography on DEAE-Sephadex A-25 and Sephadex LH-20. To obtain the Myr required for efficient hydrolysis of GLSs into antibacterial ITCs, we developed a rapid method for its extraction from the seeds of Lepidium sativum (garden cress). The yields of GLSs were 22.9 ± 1.2 mg GRN (purity 96%) and 10.4 ± 1.1 mg SBN (purity 92%) from 1 g of dry plant material. Both purified GLSs were used as substrates for the Myr. Analysis of the composition of hydrolysis products (HPs) revealed differences in their hydrolysis rates and in the degree of conversion from GLSs to individual ITCs catalyzed by Myr. When GRNs were cleaved, SFNs were formed in an equimolar ratio, but the formation of pHBITCs was only half that of cleaved SBNs. The decrease in pHBITC content is due to its instability compared to SFN. While SFN is stable in aqueous media during the measurement, pHBITC undergoes non-enzymatic hydrolysis to p-hydroxybenzyl alcohol and thiocyanate ions. Testing of the antimicrobial effects of the HPs formed from GRN by Myr under premix or in situ conditions showed inhibition of the growth of model prokaryotic and eukaryotic microorganisms. This observation could serve as the jumping-off point for the design of a two-component mixture, based on purified GLSs and Myr that is, usable in food or the pharmaceutical industry in the future.

Pseudomonas syringae infectivity correlates to altered transcript and metabolite levels of Arabidopsis mediator mutants

Rapid metabolic responses to pathogens are essential for plant survival and depend on numerous transcription factors. Mediator is the major transcriptional co-regulator for integration and transmission of signals from transcriptional regulators to RNA polymerase II. Using four Arabidopsis Mediator mutants, med16, med18, med25 and cdk8, we studied how differences in regulation of their transcript and metabolite levels correlate to their responses to Pseudomonas syringae infection. We found that med16 and cdk8 were susceptible, while med25 showed increased resistance. Glucosinolate, phytoalexin and carbohydrate levels were reduced already before infection in med16 and cdk8, but increased in med25, which also displayed increased benzenoids levels. Early after infection, wild type plants showed reduced glucosinolate and nucleoside levels, but increases in amino acids, benzenoids, oxylipins and the phytoalexin camalexin. The Mediator mutants showed altered levels of these metabolites and in regulation of genes encoding key enzymes for their metabolism. At later stage, mutants displayed defective levels of specific amino acids, carbohydrates, lipids and jasmonates which correlated to their infection response phenotypes. Our results reveal that MED16, MED25 and CDK8 are required for a proper, coordinated transcriptional response of genes which encode enzymes involved in important metabolic pathways for Arabidopsis responses to Pseudomonas syringae infections.

Analyzing the Diversity of MYB Family Response Strategies to Drought Stress in Different Flax Varieties Based on Transcriptome Data

The MYB transcription factor family has numerous members, and is involved in biological activities, such as ABA signaling, which plays an important role in a plant's resistance to abiotic stresses such as drought. However, the diversity of MYB members that respond to drought stress and their regulatory mechanisms in different flax varieties were unclear. In this study, we obtained 855.69 Gb of clean data from 120 flax root samples from 20 flax (Linum usitatissimum L.) varieties, assembled 92,861 transcripts, and identified 434 MYB family members in each variety. The expression profiles of the MYB transcription factor family from 20 flax varieties under drought stress were analyzed. The results indicated that there are four strategies by which the MYB family responds to drought stress in these 20 flax varieties, each of which has its own specific processes, such as development, reproduction, and localization processes. The four strategies also include common biological processes, such as stimulus responses, metabolic processes, and biological regulation. The WGCNA method was subsequently employed to identify key members of the MYB family involved in response strategies to drought stress. The results demonstrated that a 1R-MYB subfamily gene co-expression network is significantly related to the gibberellin response and cytokinin-activated signaling pathway processes in the 'Strategy 4' for MYB family response to drought, identifying core genes such as Lus.scaffold70.240. Our results showed a diversity of MYB family responses to drought stress within flax varieties, and these results contribute to deciphering the mechanisms of the MYB family regulation of drought resistance. This will promote the more accurate breeding development of flax to adapt to agricultural production under drought conditions.

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.

Transcriptome Analysis Reveals the Mechanism of Exogenous Selenium in Alleviating Cadmium Stress in Purple Flowering Stalks (Brassica campestris var. purpuraria)

In China, cadmium (Cd) stress has a significant role in limiting the development and productivity of purple flowering stalks (Brassica campestris var. purpuraria). Exogenous selenium supplementation has been demonstrated in earlier research to mitigate the effects of Cd stress in a range of plant species; nevertheless, the physiological and molecular processes by which exogenous selenium increases vegetable shoots' resistance to Cd stress remain unclear. Purple flowering stalks (Brassica campestris var. purpuraria) were chosen as the study subject to examine the effects of treatment with sodium selenite (Na2SeO3) on the physiology and transcriptome alterations of cadmium stress. Purple flowering stalk leaves treated with exogenous selenium had higher glutathione content, photosynthetic capacity, and antioxidant enzyme activities compared to the leaves treated with Cd stress alone. Conversely, the contents of proline, soluble proteins, soluble sugars, malondialdehyde, and intercellular CO2 concentration tended to decrease. Transcriptome analysis revealed that 2643 differentially expressed genes (DEGs) were implicated in the response of exogenous selenium treatment to Cd stress. The metabolic pathways associated with flavonoid production, carotenoid synthesis, glutathione metabolism, and glucosinolate biosynthesis were among those enriched in these differentially expressed genes. Furthermore, we discovered DEGs connected to the production route of glucosinolates. This work sheds fresh light on how purple flowering stalks' tolerance to cadmium stress is improved by exogenous selenium.

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.

Is specialized metabolite regulation specialized?

Recent technical and theoretical advances have generated an explosion in the identification of specialized metabolite pathways. In comparison, our understanding of how these pathways are regulated is relatively lagging. This and the relatively young age of specialized metabolite pathways has partly contributed to a default and common paradigm whereby specialized metabolite regulation is theorized as relatively simple with a few key transcription factors and the compounds are non-regulatory end-products. In contrast, studies into model specialized metabolites, such as glucosinolates, are beginning to identify a new understanding whereby specialized metabolites are highly integrated into the plants' core metabolic, physiological, and developmental pathways. This model includes a greatly extended compendium of transcription factors controlling the pathway, key transcription factors that co-evolve with the pathway and simultaneously control core metabolic and developmental components, and finally the compounds themselves evolve regulatory connections to integrate into the plants signaling machinery. In this review, these concepts are illustrated using studies in the glucosinolate pathway within the Brassicales. This suggests that the broader community needs to reconsider how they do or do not integrate specialized metabolism into the regulatory network of their study species.

Phytoalexins of the crucifer Barbarea vulgaris: Structural profile and correlation with glucosinolate turnover

Phytoalexins are antimicrobial plant metabolites elicited by microbial attack or abiotic stress. We investigated phytoalexin profiles after foliar abiotic elicitation in the crucifer Barbarea vulgaris and interactions with the glucosinolate-myrosinase system. The treatment for abiotic elicitation was a foliar spray with CuCl2 solution, a usual eliciting agent, and three independent experiments were carried out. Two genotypes of B. vulgaris (G-type and P-type) accumulated the same three major phytoalexins in rosette leaves after treatment: phenyl-containing nasturlexin D and indole-containing cyclonasturlexin and cyclobrassinin. Phytoalexin levels were investigated daily by UHPLC-QToF MS and tended to differ among plant types and individual phytoalexins. In roots, phytoalexins were low or not detected. In treated leaves, typical total phytoalexin levels were in the range 1-10 nmol/g fresh wt. during three days after treatment while typical total glucosinolate (GSL) levels were three orders of magnitude higher. Levels of some minor GSLs responded to the treatment: phenethylGSL (PE) and 4-substituted indole GSLs. Levels of PE, a suggested nasturlexin D precursor, were lower in treated plants than controls. Another suggested precursor GSL, 3-hydroxyPE, was not detected, suggesting PE hydrolysis to be a key biosynthetic step. Levels of 4-substituted indole GSLs differed markedly between treated and control plants in most experiments, but not in a consistent way. The dominant GSLs, glucobarbarins, are not believed to be phytoalexin precursors. We observed statistically significant linear correlations between total major phytoalexins and the glucobarbarin products barbarin and resedine, suggesting that GSL turnover for phytoalexin biosynthesis was unspecific. In contrast, we did not find correlations between total major phytoalexins and raphanusamic acid or total glucobarbarins and barbarin. In conclusion, two groups of phytoalexins were detected in B. vulgaris, apparently derived from the GSLs PE and indol-3-ylmethylGSL. Phytoalexin biosynthesis was accompanied by depletion of the precursor PE and by turnover of major non-precursor GSLs to resedine. This work paves the way for identifying and characterizing genes and enzymes in the biosyntheses of phytoalexins and resedine.

Brassinosteroids fine-tune secondary and primary sulfur metabolism through BZR1-mediated transcriptional regulation

For adaptation to ever-changing environments, plants have evolved elaborate metabolic systems coupled to a regulatory network for optimal growth and defense. Regulation of plant secondary metabolic pathways such as glucosinolates (GSLs) by defense phytohormones in response to different stresses and nutrient deficiency has been intensively investigated, while how growth-promoting hormone balances plant secondary and primary metabolism has been largely unexplored. Here, we found that growth-promoting hormone brassinosteroid (BR) inhibits GSLs accumulation while enhancing biosynthesis of primary sulfur metabolites, including cysteine (Cys) and glutathione (GSH) both in Arabidopsis and Brassica crops, fine-tuning secondary and primary sulfur metabolism to promote plant growth. Furthermore, we demonstrate that of BRASSINAZOLE RESISTANT 1 (BZR1), the central component of BR signaling, exerts distinct transcriptional inhibition regulation on indolic and aliphatic GSL via direct MYB51 dependent repression of indolic GSL biosynthesis, while exerting partial MYB29 dependent repression of aliphatic GSL biosynthesis. Additionally, BZR1 directly activates the transcription of APR1 and APR2 which encodes rate-limiting enzyme adenosine 5'-phosphosulfate reductases in the primary sulfur metabolic pathway. In summary, our findings indicate that BR inhibits the biosynthesis of GSLs to prioritize sulfur usage for primary metabolites under normal growth conditions. These findings expand our understanding of BR promoting plant growth from a metabolism perspective.

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.

Transcriptomic and epigenomic analyses revealed that polycomb repressive complex 2 regulates not only developmental but also stress responsive metabolism in Brassica rapa

Polycomb group proteins (PcG) play a crucial role in developmental programs in eukaryotic organisms, including plants. PcG-mediated gene repression is achieved by epigenetic histone modification on target chromatins. Loss of PcG components leads to severe developmental defects. CURLY LEAF (CLF), a PcG component in Arabidopsis, catalyzes the trimethylation of histone H3 on lysine 27 (H3K27me3), a repressive histone mark in numerous genes in Arabidopsis. In this study, we isolated a single homolog of Arabidopsis CLF, namely, BrCLF, in Brassica rapa ssp. trilocularis. Transcriptomic analysis revealed that BrCLF participated in B. rapa developmental processes, such as seed dormancy, leaf and flower organ development, and floral transition. BrCLF was also involved in stress signaling and stress-responsive metabolism, such as aliphatic and indolic glucosinolate metabolism in B. rapa. Epigenome analysis showed that H3K27me3 was substantially enriched in genes related to these developmental and stress-responsive processes. Thus, this study provided a basis for elucidating the molecular mechanism of the PcG-mediated regulation of development and stress responses in B. rapa.

Loss of MYB34 Transcription Factor Supports the Backward Evolution of Indole Glucosinolate Biosynthesis in a Subclade of the Camelineae Tribe and Releases the Feedback Loop in This Pathway in Arabidopsis

Glucosinolates are specialized defensive metabolites characteristic of the Brassicales order. Among them, aliphatic and indolic glucosinolates (IGs) are usually highly abundant in species from the Brassicaceae family. The exceptions this trend are species representing a subclade of the Camelineae tribe, including Capsella and Camelina genera, which have reduced capacity to produce and metabolize IGs. Our study addresses the contribution of specific glucosinolate-related myeloblastosis (MYB) transcription factors to this unprecedented backward evolution of IG biosynthesis. To this end, we performed phylogenomic and functional studies of respective MYB proteins. The obtained results revealed weakened conservation of glucosinolate-related MYB transcription factors, including loss of functional MYB34 protein, in the investigated species. We showed that the introduction of functional MYB34 from Arabidopsis thaliana partially restores IG biosynthesis in Capsella rubella, indicating that the loss of this transcription factor contributes to the backward evolution of this metabolic pathway. Finally, we performed an analysis of the impact of particular myb mutations on the feedback loop in IG biosynthesis, which drives auxin overproduction, metabolic dysregulation and strong growth retardation caused by mutations in IG biosynthetic genes. This uncovered the unique function of MYB34 among IG-related MYBs in this feedback regulation and consequently in IG conservation in Brassicaceae plants.

Glucosinolate Biosynthetic Genes of Cabbage: Genome-Wide Identification, Evolution, and Expression Analysis

Cabbage (Brassica oleracea var. capitata) is a vegetable rich in glucosinolates (GSLs) that have proven health benefits. To gain insights into the synthesis of GSLs in cabbage, we systematically analyzed GSLs biosynthetic genes (GBGs) in the entire cabbage genome. In total, 193 cabbage GBGs were identified, which were homologous to 106 GBGs in Arabidopsis thaliana. Most GBGs in cabbage have undergone negative selection. Many homologous GBGs in cabbage and Chinese cabbage differed in expression patterns indicating the unique functions of these homologous GBGs. Spraying five exogenous hormones significantly altered expression levels of GBGs in cabbage. For example, MeJA significantly upregulated side chain extension genes BoIPMILSU1-1 and BoBCAT-3-1, and the expression of core structure construction genes BoCYP83A1 and BoST5C-1, while ETH significantly repressed the expression of side chain extension genes such as BoIPMILSU1-1, BoCYP79B2-1, and BoMAMI-1, and some transcription factors, namely BoMYB28-1, BoMYB34-1, BoMYB76-1, BoCYP79B2-1, and BoMAMI-1. Phylogenetically, the CYP83 family and CYP79B and CYP79F subfamilies may only be involved in GSL synthesis in cruciferous plants. Our unprecedented identification and analysis of GBGs in cabbage at the genome-wide level lays a foundation for the regulation of GSLs synthesis through gene editing and overexpression.

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.

Ablation of glucosinolate accumulation in the oil crop Camelina sativa by targeted mutagenesis of genes encoding the transporters GTR1 and GTR2 and regulators of biosynthesis MYB28 and MYB29

Camelina sativa is an oil crop with low input costs and resistance to abiotic and biotic stresses. The presence of glucosinolates, plant metabolites with adverse health effects, restricts the use of camelina for human and animal nutrition. Cas9 endonuclease-based targeted mutagenesis of the three homeologs of each of the glucosinolate transporters CsGTR1 and CsGTR2 caused a strong decrease in glucosinolate amounts, highlighting the power of this approach for inactivating multiple genes in a hexaploid crop. Mutagenesis of the three homeologs of each of the transcription factors CsMYB28 and CsMYB29 resulted in the complete loss of glucosinolates, representing the first glucosinolate-free Brassicaceae crop. The oil and protein contents and the fatty acid composition of the csgtr1csgtr2 and csmyb28csmyb29 mutant seeds were not affected. The decrease and elimination of glucosinolates improves the quality of the oil and press cake of camelina, which thus complies with international standards regulating glucosinolate levels for human consumption and animal feeding.

Anticancer Effects of High Glucosinolate Synthesis Lines of Brassica rapa on Colorectal Cancer Cells

Chemoprevention is a method of health control in modern industrialized societies. Traditional breeding (hybridization) has been widely used to produce new (sub)species with beneficial phenotypes. Previously, we produced a number of doubled haploid (DH) lines of Brassica rapa with a high glucosinolate (GSL) content. In this study, we evaluated the anticancer activities of extracts from three selected high-GSL (HGSL)-containing DH lines (DHLs) of Brassica rapa in human colorectal cancer (CRC) cells. The three HGSL DHL extracts showed anti-proliferative activities in the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium assay and pro-apoptotic activities in the cell cycle or annexin V analysis with the induction of pro-apoptotic protein expression in CRC cells. Mechanistically, HGSL DHL extracts inhibited the NF-κB and ERK pathways, leading to a reduction in the nuclear localization of NF-κB p65. In addition, reactive oxygen species were induced by HGSL DHL extract treatment in CRC cells. In conclusion, our data suggest that the newly developed HGSL DHLs possess enhanced anticancer activities and are potentially helpful as a daily vegetable supplement with chemopreventive activities.

Jasmonate regulates plant resistance to Pectobacterium brasiliense by inducing indole glucosinolate biosynthesis

Pectobacterium brasiliense (P. brasiliense) is a necrotrophic bacterium that causes the soft rot disease in Brassica rapa. However, the mechanisms underlying plant immune responses against necrotrophic bacterial pathogens with a broad host range are still not well understood. Using a flg22-triggered seedling growth inhibition (SGI) assay with 455 Brassica rapa inbred lines, we selected six B. rapa flagellin-insensitive lines (Brfin2-7) and three B. rapa flagellin-sensitive lines (Brfs1-3). Brfin lines showed compromised flg22-induced immune responses (oxidative burst, mitogen-activated protein kinase (MAPK) activation, and seedling growth inhibition) compared to the control line R-o-18; nevertheless, they were resistant to P. brasiliense. To explain this, we analyzed the phytohormone content and found that most Brfin lines had higher P. brasiliense-induced jasmonic acid (JA) than Brfs lines. Moreover, MeJA pretreatment enhanced the resistance of B. rapa to P. brasiliense. To explain the correlation between the resistance of Brfin lines to P. brasiliense and activated JA signaling, we analyzed pathogen-induced glucosinolate (GS) content in B. rapa. Notably, in Brfin7, the neoglucobrassicin (NGBS) content among indole glucosinolates (IGS) was significantly higher than that in Brfs2 following P. brasiliense inoculation, and genes involved in IGSs biosynthesis were also highly expressed. Furthermore, almost all Brfin lines with high JA levels and resistance to P. brasiliense had higher P. brasiliense-induced NGBS levels than Brfs lines. Thus, our results show that activated JA-mediated signaling attenuates flg22-triggered immunity but enhances resistance to P. brasiliense by inducing indole glucosinolate biosynthesis in Brassica rapa. This study provides novel insights into the role of JA-mediated defense against necrotrophic bacterial pathogens within a broad host range.

Glucosinolates and Biotic Stress Tolerance in Brassicaceae with Emphasis on Cabbage: A Review

Glucosinolates (GSLs) and GSL-associated genes are receiving increasing attention from molecular biologists due to their multifunctional properties. GSLs are secondary metabolites considered to be highly active in most Brassica species. Their importance has motivated the discovery and functional analysis of the GSLs and GSL hydrolysis products involved in disease development in brassicas and other plants. Comprehensive knowledge of the GSL content of Brassica species and the molecular details of GSL-related genes will help elucidate the molecular control of this plant defense system. This report provides an overview of the current status of knowledge on GSLs, GSL biosynthesis, as well as hydrolysis related genes, and GSL hydrolysis products that regulate fungal, bacterial, and insect resistance in cabbage and other brassicas.

Structural insights into partner selection for MYB and bHLH transcription factor complexes

MYB and basic helix-loop-helix (bHLH) transcription factors form complexes to regulate diverse metabolic and developmental processes in plants. However, the molecular mechanisms responsible for MYB-bHLH interaction and partner selection remain unclear. Here, we report the crystal structures of three MYB-bHLH complexes (WER-EGL3, CPC-EGL3 and MYB29-MYC3), uncovering two MYB-bHLH interaction modes. WER and CPC are R2R3- and R3-type MYBs, respectively, but interact with EGL3 through their N-terminal R3 domain in a similar mode. A single amino acid of CPC, Met49, is crucial for competition with WER to interact with EGL3. MYB29, a R2R3-type MYB transcription factor, interacts with MYC3 by its C-terminal MYC-interaction motif. The WER-EGL3 and MYB29-MYC3 binding modes are widely applied among MYB-bHLH complexes in Arabidopsis and evolve independently in plants.

Effects of Plant Hormones, Metal Ions, Salinity, Sugar, and Chemicals Pollution on Glucosinolate Biosynthesis in Cruciferous Plant

Cruciferous vegetable crops are grown widely around the world, which supply a multitude of health-related micronutrients, phytochemicals, and antioxidant compounds. Glucosinolates (GSLs) are specialized metabolites found widely in cruciferous vegetables, which are not only related to flavor formation but also have anti-cancer, disease-resistance, and insect-resistance properties. The content and components of GSLs in the Cruciferae are not only related to genotypes and environmental factors but also are influenced by hormones, plant growth regulators, and mineral elements. This review discusses the effects of different exogenous substances on the GSL content and composition, and analyzes the molecular mechanism by which these substances regulate the biosynthesis of GSLs. Based on the current research status, future research directions are also proposed.

Sulfur compounds: From plants to humans and their role in chronic disease prevention

Sulfur is essential for the health of plants and is an indispensable dietary component for human health and disease prevention. Its incorporation into our food supply is heavily reliant upon the uptake of sulfur into plant tissue and our subsequent intake. Dietary requirements for sulfur are largely calculated based upon requirements for the sulfur-containing amino acids (SAA), cysteine and methionine, to meet the demands for synthesis of proteins, enzymes, co-enzymes, vitamins, and hormones. SAA are found in abundance in animal sources and are relatively low in plants. However, some plants, particularly cruciferous and allium vegetables, produce many protective sulfur-containing secondary metabolites, such as glucosinolates and cysteine sulfoxides. The variety and quantity of these sulfur-containing metabolites are extensive and their effects on human health are wide-reaching. Many benefits appear to be related to sulfur's role in redox biochemistry, protecting against uncontrolled oxidative stress and inflammation; features consistent within cardiometabolic dysfunction and many chronic metabolic diseases of aging. This narrative explores the origins and importance of sulfur, its incorporation into our food supply and dietary sources. It also explores the overarching potential of sulfur for human health, particularly around the amelioration of oxidative stress and chronic inflammation, and subsequent chronic disease prevention.

Characterisation of the Introgression of Brassica villosa Genome Into Broccoli to Enhance Methionine-Derived Glucosinolates and Associated Health Benefits

Broccoli cultivars that have enhanced accumulation of methionine-derived glucosinolates have been developed through the introgression of a novel allele of the MYB28 transcription factor from the wild species Brassica villosa. Through a novel k-mer approach, we characterised the extent of the introgression of unique B. villosa genome sequences into high glucosinolate broccoli genotypes. RNAseq analyses indicated that the introgression of the B. villosa MYB28 C2 allele resulted in the enhanced expression of the MYB28 transcription factor, and modified expression of genes associated with sulphate absorption and reduction, and methionine and glucosinolate biosynthesis when compared to standard broccoli. A adenine-thymine (AT) short tandem repeat (STR) was identified within the 5' untranslated region (UTR) B. villosa MYB28 allele that was absent from two divergent cultivated forms of Brassica oleracea, which may underpin the enhanced expression of B. villosa MYB28.

A Comparative Transcriptome and Metabolome Combined Analysis Reveals the Key Genes and Their Regulatory Model Responsible for Glucoraphasatin Accumulation in Radish Fleshy Taproots

Radish (Raphanus sativus L.) is rich in specific glucosinolates (GSLs), which benefit human health and special flavor formation. Although the basic GSLs metabolic pathway in Brassicaceae plants is clear, the regulating mechanism for specific glucosinolates content in radish fleshy taproots is not well understood. In this study, we discovered that there was a significant difference in the GSLs profiles and the content of various GSLs components. Glucoraphasatin (GRH) is the most predominant GSL in radish taproots of different genotypes as assessed by HPLC analysis. Further, we compared the taproot transcriptomes of three radish genotypes with high and low GSLs content by employing RNA-seq. Totally, we identified forty-one differentially expressed genes related to GSLs metabolism. Among them, thirteen genes (RsBCAT4, RsIPMDH1, RsMAM1a, RsMAM1b, RsCYP79F1, RsGSTF9, RsGGP1, RsSUR1, RsUGT74C1, RsST5b, RsAPK1, RsGSL-OH, and RsMYB28) were significantly higher co-expressed in the high content genotypes than in low content genotype. Notably, correlation analysis indicated that the expression level of RsMYB28, as an R2R3 transcription factor directly regulating aliphatic glucosinolate biosynthesis, was positively correlated with the GRH content. Co-expression network showed that RsMYB28 probably positively regulated the expression of the above genes, particularly RsSUR1, and consequently the synthesis of GRH. Moreover, the molecular mechanism of the accumulation of this 4-carbon (4C) GSL in radish taproots was explored. This study provides new perspectives on the GSLs accumulation mechanism and genetic improvements in radish taproots.

Transcriptomic and metabolic analyses revealed the modulatory effect of vernalization on glucosinolate metabolism in radish (Raphanus sativus L.)

Vernalization is the process by which long-term cold like winter triggers transition to flowering in plants. Many biennial and perennial plants including Brassicaceae family plants require vernalization for floral transition. Not only floral transition, but dynamic physiological and metabolic changes might also take place during vernalization. However, vernalization-mediated metabolic change is merely investigated so far. One of secondary metabolites found in Brassiceceae family plants is glucosinolates (GSLs). GSLs provides defense against pathogens and herbivores attack in plants and also exhibits inhibitory activity against human cancer cell. Profiles of GSLs are highly modulated by different environmental stresses in Brassciaceae family plants. To grasp the effect of vernalization on GSLs metabolic dynamics in radish (Raphanus sativus L.), we performed transcriptomic and metabolic analysis during vernalization in radish. Through transcriptome analysis, we found many GSLs metabolic genes were significantly down-regulated by vernalization in radish plants. Ultra-High Performance Liquid Chromatography analysis also revealed that GSLs compounds were substantially reduced in vernalized radish samples compared to non-vernalized radish samples. Furthermore, we found that repressive histone modification (i.e. H3K27me3) is involved in the modulation of GSLs metabolism via epigenetic suppression of Glucoraphasatin Synthase 1 (GRS1) during vernalization in radish. This study revealed that GSLs metabolism is modulated by vernalization, suggestive of a newly identified target of vernalization in radish.

A genome-scale TF-DNA interaction network of transcriptional regulation of Arabidopsis primary and specialized metabolism

Plant metabolism is more complex relative to individual microbes. In single-celled microbes, transcriptional regulation by single transcription factors (TFs) is sufficient to shift primary metabolism. Corresponding genome-level transcriptional regulatory maps of metabolism reveal the underlying design principles responsible for these shifts as a model in which master regulators largely coordinate specific metabolic pathways. Plant primary and specialized metabolism occur within innumerable cell types, and their reactions shift depending on internal and external cues. Given the importance of plants and their metabolites in providing humanity with food, fiber, and medicine, we set out to develop a genome-scale transcriptional regulatory map of Arabidopsis metabolic genes. A comprehensive set of protein-DNA interactions between Arabidopsis thaliana TFs and gene promoters in primary and specialized metabolic pathways were mapped. To demonstrate the utility of this resource, we identified and functionally validated regulators of the tricarboxylic acid (TCA) cycle. The resulting network suggests that plant metabolic design principles are distinct from those of microbes. Instead, metabolism appears to be transcriptionally coordinated via developmental- and stress-conditional processes that can coordinate across primary and specialized metabolism. These data represent the most comprehensive resource of interactions between TFs and metabolic genes in plants.

Silver nanoparticles and Chlorella treatments induced glucosinolates and kaempferol key biosynthetic genes in Eruca sativa

Background

Microalgae and nanoparticles are currently considered promising tools for numerous agricultural and biotechnological applications. The green microalga Chlorella sp. MF1 and its biosynthesized silver nanoparticles (AgNPs) were used in this study as biofortification agents to enhance glucosinolate and kaempferol levels in Eruca sativa. UV–visible spectroscopy, XRD, FTIR and TEM were comprehensively used for characterizing Chlorella-based AgNPs.

Results

The biosynthesized AgNPs were found to be spherical in shape, with size ranging from 1.45 to 5.08 nm. According to FTIR measurements, silver ions were reduced to AgNPs by functional groups such as amide, hydroxyl and carboxylate. Different experimental treatments were conducted, including either soaking seeds of E. sativa or foliar spray with various concentrations of Chlorella suspension (1, 2, 3 and 4 g L−1) and AgNPs (5, 10, 20 and 40 mg L−1). Expression levels of five key genes in the biosynthetic pathway of glucosinolates (MAM1, SUR1, MYB34 and MYB51) and kaempferol (CHS) were assessed using qRT-PCR. The results indicated an upregulation in the gene expression levels in all treatments compared to control, recording the highest level at 40 mg L−1 AgNPs and 4 g L−1Chlorella suspension. In addition, high glucosinolates and kaempferol content was detected in plants whose leaves were sprayed with AgNPs and Chlorella suspension (40 mg L−1 and 4 g L−1) based on HPLC analysis. Sequence analysis of amplified CHS fragments from E. sativa plants treated with AgNPs (40 mg L−1) showed high sequence similarity to A. thaliana CHS gene. However, there were several CHS regions with sequence polymorphism (SNPs and Indels) in foliar sprayed plants.

Conclusions

Results of this study evidenced that the application of AgNPs and Chlorella suspension increased glucosinolates and kaempferol content in E. sativa through upregulation of key genes in their biosynthetic pathway.

Multiple indole glucosinolates and myrosinases defend Arabidopsis against Tetranychus urticae herbivory

Arabidopsis (Arabidopsis thaliana) defenses against herbivores are regulated by the jasmonate (JA) hormonal signaling pathway, which leads to the production of a plethora of defense compounds. Arabidopsis defense compounds include tryptophan-derived metabolites, which limit Arabidopsis infestation by the generalist herbivore two-spotted spider mite, Tetranychus urticae. However, the phytochemicals responsible for Arabidopsis protection against T. urticae are unknown. Here, we used Arabidopsis mutants disrupted in the synthesis of tryptophan-derived secondary metabolites to identify phytochemicals involved in the defense against T. urticae. We show that of the three tryptophan-dependent pathways found in Arabidopsis, the indole glucosinolate (IG) pathway is necessary and sufficient to assure tryptophan-mediated defense against T. urticae. We demonstrate that all three IGs can limit T. urticae herbivory, but that they must be processed by myrosinases to hinder T. urticae oviposition. Putative IG breakdown products were detected in mite-infested leaves, suggesting in planta processing by myrosinases. Finally, we demonstrate that besides IGs, there are additional JA-regulated defenses that control T. urticae herbivory. Together, our results reveal the complexity of Arabidopsis defenses against T. urticae that rely on multiple IGs, specific myrosinases, and additional JA-dependent defenses.

Influence of Genotype on High Glucosinolate Synthesis Lines of Brassica rapa

This study was conducted to investigate doubled haploid (DH) lines produced between high GSL (HGSL) Brassica rapa ssp. trilocularis (yellow sarson) and low GSL (LGSL) B. rapa ssp. chinensis (pak choi) parents. In total, 161 DH lines were generated. GSL content of HGSL DH lines ranged from 44.12 to 57.04 μmol·g⁻¹·dry weight (dw), which is within the level of high GSL B. rapa ssp. trilocularis (47.46 to 59.56 μmol g⁻¹ dw). We resequenced five of the HGSL DH lines and three of the LGSL DH lines. Recombination blocks were formed between the parental and DH lines with 108,328 single-nucleotide polymorphisms in all chromosomes. In the measured GSL, gluconapin occurred as the major substrate in HGSL DH lines. Among the HGSL DH lines, BrYSP_DH005 had glucoraphanin levels approximately 12-fold higher than those of the HGSL mother plant. The hydrolysis capacity of GSL was analyzed in HGSL DH lines with a Korean pak choi cultivar as a control. Bioactive compounds, such as 3-butenyl isothiocyanate, 4-pentenyl isothiocyanate, 2-phenethyl isothiocyanate, and sulforaphane, were present in the HGSL DH lines at 3-fold to 6.3-fold higher levels compared to the commercial cultivar. The selected HGSL DH lines, resequencing data, and SNP identification were utilized for genome-assisted selection to develop elite GSL-enriched cultivars and the industrial production of potential anti-cancerous metabolites such as gluconapin and glucoraphanin.

New Insights on the Regulation of Glucosinolate Biosynthesis via COP1 and DELLA Proteins in Arabidopsis Thaliana

The biosynthesis of defensive secondary metabolites, such as glucosinolates (GSLs), is a costly process, which requires nutrients, ATP, and reduction equivalents, and, therefore, needs well-orchestrated machinery while coordinating defense and growth. We discovered that the key repressor of light signaling, the CONSTITUTIVE PHOTOMORPHOGENIC 1/SUPPRESSOR OF PHYTOCHROME A-105 (COP1/SPA) complex, is a crucial component of GSL biosynthesis regulation. Various mutants in this COP1/SPA complex exhibited a strongly reduced level of GSL and a low expression of jasmonate (JA)-dependent genes. Furthermore, cop1, which is known to accumulate DELLA proteins in the dark, shows reduced gibberellin (GA) and JA signaling, thereby phenocopying other DELLA-accumulating mutants. This phenotype can be complemented by a dominant gain-of-function allele of MYC3 and by crossing with a mutant having low DELLA protein levels. Hence, SPA1 interacts with DELLA proteins in a yeast two-hybrid screen, whereas high levels of DELLA inhibit MYC function and suppress JA signaling. DELLA accumulation leads to reduced synthesis of GSL and inhibited growth. Thus, the COP1/SPA-mediated degradation of DELLA not only affects growth but also regulates the biosynthesis of GSLs.

An Arabidopsis expression predictor enables inference of transcriptional regulators for gene modules

The regulation of gene expression by transcription factors (TFs) has been studied for a long time, but no model that can accurately predict transcriptome profiles based on TF activities currently exists. Here, we developed a computational approach, named EXPLICIT (Expression Prediction via Log-linear Combination of Transcription Factors), to construct a universal predictor for Arabidopsis to predict the expression of 29 182 non-TF genes using 1678 TFs. When applied to RNA-Seq samples from diverse tissues, EXPLICIT generated accurate predicted transcriptomes correlating well with actual expression, with an average correlation coefficient of 0.986. After recapitulating the quantitative relationships between TFs and their target genes, EXPLICIT enabled downstream inference of TF regulators for genes and gene modules functioning in diverse plant pathways, including those involved in suberin, flavonoid, glucosinolate metabolism, lateral root, xylem, secondary cell wall development or endoplasmic reticulum stress response. Our approach showed a better ability to recover the correct TF regulators when compared with existing plant tools, and provides an innovative way to study transcriptional regulation.

Transcriptome, proteome, and metabolome reveal the mechanism of tolerance to selenate toxicity in Cardamine violifolia

Cardamine violifolia was found here to accumulate selenium (Se) to over 9000 mg kg^-1 dry weight. To investigate the mechanism of Se accumulation and tolerance in C. violifolia, metabolome, transcriptome, and proteome technologies were applied to C. violifolia seedlings treated with selenate. Several sulfate transporter (Sultr) genes (Sultr1;1, Sultr1;2, and Sultr2;1) and sulfur assimilatory enzyme genes showed high expression levels in response to selenate. Many calcium protein and cysteine-rich kinase genes of C. violifolia were downregulated, whereas selenium-binding protein 1 (SBP1) and protein sulfur deficiency-induced 2 (SDI2) of C. violifolia were upregulated by selenate. The expression of genes involved in the ribosome and posttranslational modifications and chaperones in C. violifolia were also detected in response to selenate. Based on the results of this study and previous findings, we suggest that the downregulated expression of calcium proteins and cysteine-rich kinases, and the upregulated expression of SBP1 and SDI2, were important contributors to the Se tolerance of C. violifolia. The downregulation of cysteine-rich kinases and calcium proteins would enhance Se tolerance of C. violifolia is a novel proposition that has not been reported on other Se hyperaccumulators. This study provides us novel insights to understand Se accumulation and tolerance in plants.

ADAP is a possible negative regulator of glucosinolate biosynthesis in Arabidopsis thaliana based on clustering and gene expression analyses

Glucosinolates (GSLs) are plant secondary metabolites consisting of sulfur and nitrogen, commonly found in Brassicaceae crops, such as Arabidopsis thaliana. These compounds are known for their roles in plant defense mechanisms against pests and pathogens. 'Guilt-by-association' (GBA) approach predicts genes encoding proteins with similar function tend to share gene expression pattern generated from high throughput sequencing data. Recent studies have successfully identified GSL genes using GBA approach, followed by targeted verification of gene expression and metabolite data. Therefore, a GSL co-expression network was constructed using known GSL genes obtained from our in-house database, SuCComBase. DPClusO was used to identify subnetworks of the GSL co-expression network followed by Fisher's exact test leading to the discovery of a potential gene that encodes the ARIA-interacting double AP2-domain protein (ADAP) transcription factor (TF). Further functional analysis was performed using an effective gene silencing system known as CRES-T. By applying CRES-T, ADAP TF gene was fused to a plant-specific EAR-motif repressor domain (SRDX), which suppresses the expression of ADAP target genes. In this study, ADAP was proposed as a negative regulator in aliphatic GSL biosynthesis due to the over-expression of downstream aliphatic GSL genes (UGT74C1 and IPMI1) in ADAP-SRDX line. The significant over-expression of ADAP gene in the ADAP-SRDX line also suggests the behavior of the TF that negatively affects the expression of UGT74C1 and IPMI1 via a feedback mechanism in A. thaliana.

Fluorescent protein-based imaging and tissue-specific RNA-seq analysis of Arabidopsis hydathodes

Hydathodes are typically found at leaf teeth in vascular plants and are involved in water release to the outside. Although morphological and physiological analysis of hydathodes has been performed in various plants, little is known about the genes involved in hydathode function. In this study, we performed fluorescent protein-based imaging and tissue-specific RNA-seq analysis in Arabidopsis hydathodes. We used the enhancer trap line E325, which has been reported to express green fluorescent protein (GFP) at its hydathodes. We found that E325-GFP was expressed in small cells found inside the hydathodes (named E cells) that were distributed between the water pores and xylem ends. No fluorescence of the phloem markers pSUC2:GFP and pSEOR1:SEOR1-YFP was observed in the hydathodes. These observations indicate that Arabidopsis hydathodes are composed of three major components: water pores, xylem ends, and E cells. In addition, we performed transcriptome analysis of the hydathode using the E325-GFP line. Microsamples were collected from GFP-positive or -negative regions of E325 leaf margins with a needle-based device (~130 µm in diameter). RNA-seq was performed with each single microsample using a high-throughput library preparation method called Lasy-Seq. We identified 72 differentially expressed genes. Among them, 68 genes showed significantly higher and four genes showed significantly lower expression in the hydathode. Our results provide new insights into the molecular basis for hydathode physiology and development.

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.

Regulation of glucosinolate biosynthesis

Glucosinolates are secondary defense metabolites produced by plants of the order Brassicales, which includes the model species Arabidopsis and many crop species. In the past 13 years, the regulation of glucosinolate synthesis in plants has been intensively studied, with recent research revealing complex molecular mechanisms that connect glucosinolate production with responses to other central pathways. In this review, we discuss how the regulation of glucosinolate biosynthesis is ecologically relevant for plants, how it is controlled by transcription factors, and how this transcriptional machinery interacts with hormonal, environmental, and epigenetic mechanisms. We present the central players in glucosinolate regulation, MYB and basic helix-loop-helix transcription factors, as well as the plant hormone jasmonate, which together with other hormones and environmental signals allow the coordinated and rapid regulation of glucosinolate genes. Furthermore, we highlight the regulatory connections between glucosinolates, auxin, and sulfur metabolism and discuss emerging insights and open questions on the regulation of glucosinolate biosynthesis.

Arabidopsis thaliana transcription factors MYB28 and MYB29 shape ammonium stress responses by regulating Fe homeostasis

Although ammonium (NH₄⁺ ) is a key intermediate of plant nitrogen metabolism, high concentrations of NH₄⁺ in the soil provoke physiological disorders that lead to the development of stress symptoms. Ammonium nutrition was shown to induce the accumulation of glucosinolates (GSLs) in leaves of different Brassicaceae species. To further understand the link between ammonium nutrition and GSLs, we analysed the ammonium stress response of Arabidopsis mutants impaired in GSL metabolic pathway. We showed that the MYB28 and MYB29 double mutant (myb28myb29), which is almost deprived of aliphatic GSLs, is highly hypersensitive to ammonium nutrition. Moreover, we evidenced that the stress symptoms developed were not a consequence of the lack of aliphatic GSLs. Transcriptomic analysis highlighted the induction of an iron (Fe) deficiency response in myb28myb29 under ammonium nutrition. Consistently, ammonium-grown myb28myb29 plants showed altered Fe accumulation and homeostasis. Interestingly, we showed overall that growing Arabidopsis with increased Fe availability relieved ammonium stress symptoms and that this was associated with MYB28 and MYB29 expression. Taken together, our data indicated that the control of Fe homeostasis was crucial for the Arabidopsis response to ammonium nutrition and evidenced that MYB28 and MYB29 play a role in this control.

A bipartite transcription factor module controlling expression in the bundle sheath of Arabidopsis thaliana

C₄ photosynthesis evolved repeatedly from the ancestral C₃ state, improving photosynthetic efficiency by ~50%. In most C₄ lineages, photosynthesis is compartmented between mesophyll and bundle sheath cells, but how gene expression is restricted to these cell types is poorly understood. Using the C₃ model Arabidopsis thaliana, we identified cis-elements and transcription factors driving expression in bundle sheath strands. Upstream of the bundle sheath preferentially expressed MYB76 gene, we identified a region necessary and sufficient for expression containing two cis-elements associated with the MYC and MYB families of transcription factors. MYB76 expression is reduced in mutant alleles for these transcription factors. Moreover, downregulated genes shared by both mutants are preferentially expressed in the bundle sheath. Our findings are broadly relevant for understanding the spatial patterning of gene expression, provide specific insights into mechanisms associated with the evolution of C₄ photosynthesis and identify a short tuneable sequence for manipulating gene expression in the bundle sheath.

The Eruca sativa Genome and Transcriptome: A Targeted Analysis of Sulfur Metabolism and Glucosinolate Biosynthesis Pre and Postharvest

Rocket (Eruca sativa) is a source of health-related metabolites called glucosinolates (GSLs) and isothiocyanates (ITCs) but little is known of the genetic and transcriptomic mechanisms responsible for regulating pre and postharvest accumulations. We present the first de novo reference genome assembly and annotation, with ontogenic and postharvest transcriptome data relating to sulfur assimilation, transport, and utilization. Diverse gene expression patterns related to sulfur metabolism, GSL biosynthesis, and glutathione biosynthesis are present between inbred lines of rocket. A clear pattern of differential expression determines GSL abundance and the formation of hydrolysis products. One breeding line sustained GSL accumulation and hydrolysis product formation throughout storage. Multiple copies of MYB28, SLIM1, SDI1, and ESM1 have increased and differential expression postharvest, and are associated with GSLs and hydrolysis product formation. Two glucosinolate transporter gene (GTR2) copies were found to be associated with increased GSL accumulations in leaves. Monosaccharides (which are essential for primary metabolism and GSL biosynthesis, and contribute to the taste of rocket) were also quantified in leaves, with glucose concentrations significantly correlated with the expression of numerous GSL-related genes. Significant negative correlations were observed between the expression of glutathione synthetase (GSH) genes and those involved in GSL metabolism. Breeding line "B" showed increased GSH gene expression and low GSL content compared to two other lines where the opposite was observed. Co-expression analysis revealed senescence (SEN1) and oxidative stress-related (OXS3) genes have higher expression in line B, suggesting that postharvest deterioration is associated with low GSL concentrations.

Association Mapping of Seed Quality Traits Under Varying Conditions of Nitrogen Application in Brassica juncea L. Czern & Coss

Indian mustard (Brassica juncea) is a major source of vegetable oil in the Indian subcontinent. The seed cake left after the oil extraction is used as livestock feed. We examined the genetic architecture of oil, protein, and glucosinolates by conducting a genome-wide association study (GWAS), using an association panel comprising 92 diverse genotypes. We conducted trait phenotyping over 2 years at two levels of nitrogen (N) application. Genotyping by sequencing was used to identify 66,835 loci, covering 18 chromosomes. Genetic diversity and phenotypic variations were high for the studied traits. Trait performances were stable when averaged over years and N levels. However, individual performances differed. General and mixed linear models were used to estimate the association between the SNP markers and the seed quality traits. Population structure, principal components (PCs) analysis, and discriminant analysis of principal components (DAPCs) were included as covariates to overcome the bias due to the population stratification. We identified 16, 23, and 27 loci associated with oil, protein, and glucosinolates, respectively. We also established LD patterns and haplotype structures for the candidate genes. The average block sizes were larger on A-genome chromosomes as compared to the B- genome chromosomes. Genetic associations differed over N levels. However, meta-analysis of GWAS datasets not only improved the power to recognize associations but also helped to identify common SNPs for oil and protein contents. Annotation of the genomic region around the identified SNPs led to the prediction of 21 orthologs of the functional candidate genes related to the biosynthesis of oil, protein, and glucosinolates. Notable among these are: LACS5 (A09), FAD6 (B05), ASN1 (A06), GTR2 (A06), CYP81G1 (B06), and MYB44 (B06). The identified loci will be very useful for marker-aided breeding for seed quality modifications in B. juncea.

Induction of Glucoraphasatin Biosynthesis Genes by MYB29 in Radish (Raphanus sativus L.) Roots

Glucoraphasatin (GRH) is a specific aliphatic glucosinolate (GSL) that is only abundant in radish (Raphanus sativus L.). The gene expression regulating GRH biosynthesis in radish is still poorly understood. We employed a total of 59 radish accessions to analyze GSL profiles and showed that GRH was specific and predominant among the aliphatic GSLs in radish roots. We selected five accessions roots with high, moderate and low GSL biosynthesis, respectively, to conduct a comparative transcriptome analysis and the qRT-PCR of the biosynthesis genes for aliphatic GSLs. In this study, among all the accessions tested, roots with the accession RA157-74 had a high GRH content and showed a significant expression of the aliphatic GSL biosynthesis genes. We defined the genes involved in the GRH biosynthesis process and found that they were regulated by a transcription factor (RSG00789) at the MYB29 locus in radish roots. We found 13 aliphatic GSL biosynthesis genes regulated by the RSG00789 gene in the GRH biosynthesis pathway.

FRS7 and FRS12 recruit NINJA to regulate expression of glucosinolate biosynthesis genes

The sessile lifestyle of plants requires accurate physiology adjustments to be able to thrive in a changing environment. Plants integrate environmental timing signals to control developmental and stress responses. Here, we identified Far1 Related Sequence (FRS) 7 and FRS12, two transcriptional repressors that accumulate in short-day conditions, as regulators of Arabidopsis glucosinolate (GSL) biosynthesis. Loss of function of FRS7 and FRS12 results in plants with increased amplitudes of diurnal expression of GSL pathway genes. Protein interaction analyses revealed that FRS7 and FRS12 recruit the NOVEL INTERACTOR OF JAZ (NINJA) to assemble a transcriptional repressor complex. Genetic and molecular evidence demonstrated that FRS7, FRS12 and NINJA jointly regulate the expression of GSL biosynthetic genes, and thus constitute a molecular mechanism that modulates specialized metabolite accumulation.

Involvement of BGLU30 in Glucosinolate Catabolism in the Arabidopsis Leaf under Dark Conditions

Glucosinolates (GSLs) are secondary metabolites that play important roles in plant defense and are suggested to act as storage compounds. Despite their important roles, metabolic dynamics of GSLs under various growth conditions remain poorly understood. To determine how light conditions influence the levels of different GSLs and their distribution in Arabidopsis leaves, we visualized the GSLs under different light conditions using matrix-assisted laser desorption/ionization mass spectrometry imaging. We observed the unique distribution patterns of each GSL in the inner regions of leaves and marked decreases under darkness, indicating light conditions influenced GSL metabolism. GSLs are hydrolyzed by a group of ß-glucosidase (BGLU) called myrosinase. Previous transcriptome data for GSL metabolism under light and dark conditions have revealed the highly induced expression of BGLU30, one of the putative myrosinases, which is also annotated as Dark INducible2, under darkness. Impairment of the darkness-induced GSL decrease in the disruption mutants of BGLU30, bglu30, indicated that BGLU30 mediated GSL hydrolysis under darkness. Based on the GSL profiles in the wild-type and bglu30 leaves under both conditions, short-chain GSLs were potentially preferable substrates for BGLU30. Our findings provide an effective way of visualizing GSL distribution in plants and highlighted the carbon storage GSL function.

mGWAS Uncovers Gln-Glucosinolate Seed-Specific Interaction and its Role in Metabolic Homeostasis

Gln is a key player in plant metabolism. It is one of the major free amino acids that is transported into the developing seed and is central for nitrogen metabolism. However, Gln natural variation and its regulation and interaction with other metabolic processes in seeds remain poorly understood. To investigate the latter, we performed a metabolic genome-wide association study (mGWAS) of Gln-related traits measured from the dry seeds of the Arabidopsis (Arabidopsis thaliana) diversity panel using all potential ratios between Gln and the other members of the Glu family as traits. This semicombinatorial approach yielded multiple candidate genes that, upon further analysis, revealed an unexpected association between the aliphatic glucosinolates (GLS) and the Gln-related traits. This finding was confirmed by an independent quantitative trait loci mapping and statistical analysis of the relationships between the Gln-related traits and the presence of specific GLS in seeds. Moreover, an analysis of Arabidopsis mutants lacking GLS showed an extensive seed-specific impact on Gln levels and composition that manifested early in seed development. The elimination of GLS in seeds was associated with a large effect on seed nitrogen and sulfur homeostasis, which conceivably led to the Gln response. This finding indicates that both Gln and GLS play key roles in shaping the seed metabolic homeostasis. It also implies that select secondary metabolites might have key functions in primary seed metabolism. Finally, our study shows that an mGWAS performed on dry seeds can uncover key metabolic interactions that occur early in seed development.

Epistatic Transcription Factor Networks Differentially Modulate Arabidopsis Growth and Defense

Plants integrate internal and external signals to finely coordinate growth and defense for maximal fitness within a complex environment. A common model suggests that growth and defense show a trade-offs relationship driven by energy costs. However, recent studies suggest that the coordination of growth and defense likely involves more conditional and intricate connections than implied by the trade-off model. To explore how a transcription factor (TF) network may coordinate growth and defense, we used a high-throughput phenotyping approach to measure growth and flowering in a set of single and pairwise mutants previously linked to the aliphatic glucosinolate (GLS) defense pathway. Supporting a link between growth and defense, 17 of the 20 tested defense-associated TFs significantly influenced plant growth and/or flowering time. The TFs' effects were conditional upon the environment and age of the plant, and more critically varied across the growth and defense phenotypes for a given genotype. In support of the coordination model of growth and defense, the TF mutant's effects on short-chain aliphatic GLS and growth did not display a simple correlation. We propose that large TF networks integrate internal and external signals and separately modulate growth and the accumulation of the defensive aliphatic GLS.

Glucosinolate structural diversity, identification, chemical synthesis and metabolism in plants

The glucosinolates (GSLs) is a well-defined group of plant metabolites characterized by having an S-β-d-glucopyrano unit anomerically connected to an O-sulfated (Z)-thiohydroximate function. After enzymatic hydrolysis, the sulfated aglucone can undergo rearrangement to an isothiocyanate, or form a nitrile or other products. The number of GSLs known from plants, satisfactorily characterized by modern spectroscopic methods (NMR and MS) by mid-2018, is 88. In addition, a group of partially characterized structures with highly variable evidence counts for approximately a further 49. This means that the total number of characterized GSLs from plants is somewhere between 88 and 137. The diversity of GSLs in plants is critically reviewed here, resulting in significant discrepancies with previous reviews. In general, the well-characterized GSLs show resemblance to C-skeletons of the amino acids Ala, Val, Leu, Trp, Ile, Phe/Tyr and Met, or to homologs of Ile, Phe/Tyr or Met. Insufficiently characterized, still hypothetic GSLs include straight-chain alkyl GSLs and chain-elongated GSLs derived from Leu. Additional reports (since 2011) of insufficiently characterized GSLs are reviewed. Usually the crucial missing information is correctly interpreted NMR, which is the most effective tool for GSL identification. Hence, modern use of NMR for GSL identification is also reviewed and exemplified. Apart from isolation, GSLs may be obtained by organic synthesis, allowing isotopically labeled GSLs and any kind of side chain. Enzymatic turnover of GSLs in plants depends on a considerable number of enzymes and other protein factors and furthermore depends on GSL structure. Identification of GSLs must be presented transparently and live up to standard requirements in natural product chemistry. Unfortunately, many recent reports fail in these respects, including reports based on chromatography hyphenated to MS. In particular, the possibility of isomers and isobaric structures is frequently ignored. Recent reports are re-evaluated and interpreted as evidence of the existence of "isoGSLs", i.e. non-GSL isomers of GSLs in plants. For GSL analysis, also with MS-detection, we stress the importance of using authentic standards.