MAM3 catalyzes the formation of all aliphatic glucosinolate chain lengths in Arabidopsis

Characterization of Arabidopsis eskimo1 reveals a metabolic link between xylan O-acetylation and aliphatic glucosinolate metabolism

Glucuronoxylan is present mainly in the dicot of the secondary cell walls, often O-acetylated, which stabilizes cell structure by maintaining interaction with cellulose and other cell wall components. Some members of the Golgi localized Trichome Birefringence-Like (TBL) family function as xylan O-acetyl transferase (XOAT). The primary XOAT in the stem of Arabidopsis is ESKIMO1/TBL29, and its disruption results in decreased xylan acetylation, stunted plant growth, and collapsed xylem vessels. To elucidate the effect on metabolic reprogramming and identify the underlying cause of the stunted growth in eskimo1, we performed transcriptomic, targeted, and untargeted metabolome analysis, mainly in the inflorescence stem tissue. RNA sequencing analysis revealed that the genes involved in the biosynthesis, regulation, and transport of aliphatic glucosinolates (GSLs) were upregulated, whereas those responsible for indolic GSL metabolism were unaffected in the eskimo1 inflorescence stem. Consistently, aliphatic GSLs, such as 4-methylsulfinylbutyl (4MSOB), were increased in stem tissues and seeds. This shift in the profile of aliphatic GSLs in eskimo1 was further supported by the quantification of the soluble acetate, decrease in accumulation of GSL precursor, i.e., methionine, and increase in the level of jasmonic acid.

Bivariate GWA mapping reveals associations between aliphatic glucosinolates and plant responses to thrips and heat stress

Although plants harbor a huge phytochemical diversity, only a fraction of plant metabolites is functionally characterized. In this work, we aimed to identify the genetic basis of metabolite functions during harsh environmental conditions in Arabidopsis thaliana. With machine learning algorithms we predicted stress-specific metabolomes for 23 (a)biotic stress phenotypes of 300 natural Arabidopsis accessions. The prediction models identified several aliphatic glucosinolates (GLSs) and their breakdown products to be implicated in responses to heat stress in siliques and herbivory by Western flower thrips, Frankliniella occidentalis. Bivariate GWA mapping of the metabolome predictions and their respective (a)biotic stress phenotype revealed genetic associations with MAM, AOP, and GS-OH, all three involved in aliphatic GSL biosynthesis. We, therefore, investigated thrips herbivory on AOP, MAM, and GS-OH loss-of-function and/or overexpression lines. Arabidopsis accessions with a combination of MAM2 and AOP3, leading to 3-hydroxypropyl dominance, suffered less from thrips feeding damage. The requirement of MAM2 for this effect could, however, not be confirmed with an introgression line of ecotypes Cvi and Ler, most likely due to other, unknown susceptibility factors in the Ler background. However, AOP2 and GS-OH, adding alkenyl or hydroxy-butenyl groups, respectively, did not have major effects on thrips feeding. Overall, this study illustrates the complex implications of aliphatic GSL diversity in plant responses to heat stress and a cell-content-feeding herbivore.

Kinetic and catalytic mechanisms of the methionine-derived glucosinolate biosynthesis enzyme methylthioalkylmalate synthase

In Brassica plants, methionine-derived aliphatic glucosinolates are chemically diverse natural products that serve as plant defense compounds, as well as molecules with dietary health-promoting effects. During their biosynthesis, methylthioalkylmalate synthase (MAMS) catalyzes the elongation reaction of the aliphatic chain. The MAMS-catalyzed condensation of 4-methylthio-2-oxobutanoic acid and acetyl-CoA generates a 2-malate derivative that either enters the pathway for the synthesis of C3-glucosinolates or undergoes additional extension reactions, which lead to C4- to C9-glucosinolates. Recent determination of the x-ray crystal structure of MAMS from Brassica juncea (Indian mustard) provided insight on the molecular evolution of MAMS, especially substrate specificity changes, from the leucine biosynthesis enzyme α-isopropylmalate synthase but left details of the reaction mechanism unanswered. Here we use the B. juncea MAMS2A (BjMAMS2A) isoform to analyze the kinetic and catalytic mechanisms of this enzyme. Initial velocity studies indicate that MAMS follows an ordered bi bi kinetic mechanism, which based on the x-ray crystal structure, involves binding of 4-methylthio-2-oxobutanoic acid followed by acetyl-CoA. Examination of the pH-dependence of kcat and kcat/Km are consistent with acid/base catalysis. Site-directed mutagenesis of three residues originally proposed to function in the reaction mechanism-Arg89 (R89A, R89K, R89Q), Glu227 (E227A, E227D, E227Q), and His388 (H388A, H388N, H388Q, H388D, and H388E)-showed that only two mutants (E227Q and H388N) retained activity. Based on available structural and biochemical data, a revised reaction mechanism for MAMS-catalyzed elongation of methionine-derived aliphatic glucosinolates is proposed, which is likely also conserved in α-isopropylmalate synthase from leucine biosynthesis in plants and microbes.

Metabolomic and Transcriptomic Profiles in Diverse Brassica oleracea Crops Provide Insights into the Genetic Regulation of Glucosinolate Profiles

Glucosinolates (GSLs) are plant secondary metabolites commonly found in the cruciferous vegetables of the Brassicaceae family, offering health benefits to humans and defense against pathogens and pests to plants. In this study, we investigated 23 GSL compounds' relative abundance in four tissues of five different Brassica oleracea morphotypes. Using the five corresponding high-quality B. oleracea genome assemblies, we identified 183 GSL-related genes and analyzed their expression with mRNA-Seq data. GSL abundance and composition varied strongly, among both tissues and morphotypes, accompanied by different gene expression patterns. Interestingly, broccoli exhibited a nonfunctional AOP2 gene due to a conserved 2OG-FeII_Oxy domain loss, explaining the unique accumulation of two health-promoting GSLs. Additionally, transposable element (TE) insertions were found to affect the gene structure of MAM3 genes. Our findings deepen the understanding of GSL variation and genetic regulation in B. oleracea morphotypes, providing valuable insights for breeding with tailored GSL profiles in these crops.

4-Aldrithiol-based photometric assay for detection of methylthioalkylmalate synthase activity

In Brassica plants, glucosinolates are a diverse class of natural products, of which aliphatic methionine-derived glucosinolates are the most abundant form. Their structural diversity comes from the elongation of some side-chains by up to 9 carbons, which, after the formation of the core glucosinolate structure, can undergo further chemical modifications. Methylthioalkylmalate synthase (MAMS) catalyzes the iterative elongation process for aliphatic methionine-derived glucosinolates. Most biochemical studies on MAMS have been performed using liquid chromatography/mass spectrometry (LC/MS)-based assays or high-performance liquid chromatography (HPLC)-based assays. The LC/MS- and HPLC-based methods are endpoint assays, which cannot be monitored in real time and require a laborious process for data collection. These analytical methods are inefficient for performing multiple enzymatic assays needed to determine steady-state kinetic parameters or for mechanistic evaluation of pH-dependence and kinetic isotope effect studies. Although the function of MAMS has long been defined, there is a gap in knowledge as it pertains to biochemical characterization of this plant enzyme. Part of this may be due to the lack of efficient methods that can be used for this type of research. This chapter describes a continuous photometric assay to track MAMS activity in real time using the 4-aldrithiol reagent for reaction detection.

Characterization of unique EDTA-insensitive methylthioalkylmalate synthase from Eutrema japonicum and its potential application in synthetic biology

6-(Methylsulfinyl)hexyl isothiocyanate (6-MSITC), a derivative of glucosinolate with a six-carbon chain, is a compound found in wasabi and has diverse health-promoting properties. The biosynthesis of glucosinolates from methionine depends on a crucial step catalyzed methylthioalkylmalate synthases (MAMs), which are responsible for the generation of glucosinolates with varying chain lengths. In this study, our primary focus was the characterization of two methylthioalkyl malate synthases, MAM1-1 and MAM1-2, derived from Eutrema japonicum, commonly referred to as Japanese wasabi. Eutremajaponicum MAMs (EjMAMs) were expressed in an Escherichiacoli expression system, subsequently purified, and in vitro enzymatic activity was assayed. We explored the kinetic properties, optimal pH conditions, and cofactor preferences of EjMAMs and compared them with those of previously documented MAMs. Surprisingly, EjMAM1-2, categorized as a metallolyase family enzyme, displayed 20% of its maximum activity even in the absence of divalent metal cofactors or under high concentrations of EDTA. Additionally, we utilized AlphaFold2 to generate structural homology models of EjMAMs, and used in silico analysis and mutagenesis studies to investigate the key residues participating in catalytic activity. Moreover, we examined in vivo biosynthesis in E. coli containing Arabidopsis thaliana branched-chain amino acid transferase 3 (AtBCAT3) along with AtMAMs or EjMAMs and demonstrated that EjMAM1-2 exhibited the highest conversion rate among those MAMs, converting l-methionine to 2-(2-methylthio) ethyl malate (2-(2-MT)EM). EjMAM1-2 shows a unique property in vitro and highest activity on converting l-methionine to 2-(2-MT)EM in vivo which displays high potential for isothiocyanate biosynthesis in E. coli platform.

Unravelling Glucoraphanin and Glucoerucin Metabolism across Broccoli Sprout Development: Insights from Metabolite and Transcriptome Analysis

Variations in the concentration of glucoraphanin (GRA) and glucoerucin (GER), as well as the corresponding breakdown products, isothiocyanates (ITCs) and nitriles, were investigated during the growth of broccoli sprouts. The concentrations of GRA and GER decreased sharply from 33.66 µmol/g to 11.48 µmol/g and 12.98 µmol/g to 8.23 µmol/g, respectively, after seed germination. From the third to the seventh day, both GRA and GER were maintained as relatively stable. The highest concentrations of sulforaphane (17.16 µmol/g) and erucin (12.26 µmol/g) were observed on the first day. Hereafter, the concentrations of nitrile hydrolyzed from GRA or GER were higher than those of the corresponding ITCs. Moreover, the ratio of sulforaphane to sulforaphane nitrile decreased from 1.35 to 0.164 from 1 d to 5 d, with a similar trend exhibited for erucin/erucin nitrile after 2 d. RNA-seq analysis showed that BolMYB28 and BolCYP83A1, involved in aliphatic glucosinolate (GSL) biosynthesis, remained largely unexpressed until the third day. In contrast, the genes operating within the GSL-myrosinase hydrolysis pathway were highly expressed right from the beginning, with their expression levels increasing significantly after the third day. Additionally, we identified two BolESPs and six BolNSPs that might play important roles in promoting the production of nitriles during the development of broccoli sprouts.

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.

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.

The domestication-associated L1 gene encodes a eucomic acid synthase pleiotropically modulating pod pigmentation and shattering in soybean

Pod coloration is a domestication-related trait in soybean, with modern cultivars typically displaying brown or tan pods, while their wild relative, Glycine soja, possesses black pods. However, the factors regulating this color variation remain unknown. In this study, we cloned and characterized L1, the classical locus responsible for black pods in soybean. By using map-based cloning and genetic analyses, we identified the causal gene of L1 and revealed that it encodes a hydroxymethylglutaryl-coenzyme A (CoA) lyase-like (HMGL-like) domain protein. Biochemical assays showed that L1 functions as a eucomic acid synthase and facilitates the synthesis of eucomic acid and piscidic acid, both of which contribute to coloration of pods and seed coats in soybean. Interestingly, we found that L1 plants are more prone to pod shattering under light exposure than l1 null mutants because dark pigmentation increases photothermal efficiency. Hence, pleiotropic effects of L1 on pod color and shattering, as well as seed pigmentation, likely contributed to the preference for l1 alleles during soybean domestication and improvement. Collectively, our study provides new insights into the mechanism of pod coloration and identifies a new target for future de novo domestication of legume crops.

Allyl isothiocyanate and 6-(methylsulfinyl) hexyl isothiocyanate contents vary among wild and cultivated wasabi (Eutrema japonium)

Wasabi (Japanese horseradish, Eutrema japonicum) is the only cultivated species in the genus Eutrema with functional components that provide a strong pungent flavor. To evaluate genetic resources for wasabi breeding, we surveyed variations in the two most abundant isothiocyanate (ITC) components in wasabi, allyl isothiocyanate (AITC) and 6-methylsulfinyl (hexyl) isothiocyanate (6-MSITC, hexaraphane). We also examined the phylogenetic relationships among 36 accessions of wild and cultivated wasabi in Japan using chloroplast DNA analysis. Our results showed that (i) the 6-MSITC content in currently cultivated wasabi accessions was significantly higher than in escaped cultivars, whereas the AITC content was not significantly different. (ii) Additionally, the 6-MSITC content in cultivated wasabi was significantly lower in the spring than during other seasons. This result suggested that the 6-MSITC content responds to environmental conditions. (iii) The phylogenetic position and the 6-MSITC content of accessions from Rebun, Hokkaido Prefecture had different profiles compared with those from southern Honshu, Japan, indicating heterogeneity of the Rebun populations from other Japanese wasabi accessions. (iv) The total content of AITC and 6-MSITC in cultivated wasabi was significantly higher than that of wild wasabi. In conclusion, old cultivars or landraces of wasabi, "zairai", are the most suitable candidates for immediate use as genetic resources.

Formation, immunomodulatory activities, and enhancement of glucosinolates and sulforaphane in broccoli sprouts: a review for maximizing the health benefits to human

Broccoli sprouts have been considered as functional foods which have received increasing attention because they have been highly prized for glucosinolates, phenolics, and vitamins in particular glucosinolates. One of hydrolysates-sulforaphane from glucoraphanin is positively associated with the attenuation of inflammatory, which could reduce diabetes, cardiovascular and cancer risk. In recent decades, the great interest in natural bioactive components especially for sulforaphane promotes numerous researchers to investigate the methods to enhance glucoraphanin levels in broccoli sprouts and evaluate the immunomodulatory activities of sulforaphane. Therefore, glucosinolates profiles are different in broccoli sprouts varied with genotypes and inducers. Physicochemical, biological elicitors, and storage conditions were widely studied to promote the accumulation of glucosinolates and sulforaphane in broccoli sprouts. These inducers would stimulate the biosynthesis pathway gene expression and enzyme activities of glucosinolates and sulforaphane to increase the concentration in broccoli sprouts. The immunomodulatory activity of sulforaphane was summarized to be a new therapy for diseases with immune dysregulation. The perspective of this review served as a potential reference for customers and industries by application of broccoli sprouts as a functional food and clinical medicine.

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.

Reduced glucosinolate content in oilseed rape (Brassica napus L.) by random mutagenesis of BnMYB28 and BnCYP79F1 genes

The presence of anti-nutritive compounds like glucosinolates (GSLs) in the rapeseed meal severely restricts its utilization as animal feed. Therefore, reducing the GSL content to < 18 µmol/g dry weight in the seeds is a major breeding target. While candidate genes involved in the biosynthesis of GSLs have been described in rapeseed, comprehensive functional analyses are missing. By knocking out the aliphatic GSL biosynthesis genes BnMYB28 and BnCYP79F1 encoding an R2R3 MYB transcription factor and a cytochrome P450 enzyme, respectively, we aimed to reduce the seed GSL content in rapeseed. After expression analyses on single paralogs, we used an ethyl methanesulfonate (EMS) treated population of the inbred winter rapeseed 'Express617' to detect functional mutations in the two gene families. Our results provide the first functional analysis by knock-out for the two GSL biosynthesis genes in winter rapeseed. We demonstrate that independent knock-out mutants of the two genes possessed significantly reduced seed aliphatic GSLs, primarily progoitrin. Compared to the wildtype Express617 control plants (36.3 µmol/g DW), progoitrin levels were decreased by 55.3% and 32.4% in functional mutants of BnMYB28 (16.20 µmol/g DW) and BnCYP79F1 (24.5 µmol/g DW), respectively. Our study provides a strong basis for breeding rapeseed with improved meal quality in the future.

Genetic Improvement of Camelina sativa (L.) Crantz: Opportunities and Challenges

In recent years, a renewed interest in novel crops has been developing due to the environmental issues associated with the sustainability of agricultural practices. In particular, a cover crop, Camelina sativa (L.) Crantz, belonging to the Brassicaceae family, is attracting the scientific community's interest for several desirable features. It is related to the model species Arabidopsis thaliana, and its oil extracted from the seeds can be used either for food and feed, or for industrial uses such as biofuel production. From an agronomic point of view, it can grow in marginal lands with little or no inputs, and is practically resistant to the most important pathogens of Brassicaceae. Although cultivated in the past, particularly in northern Europe and Italy, in the last century, it was abandoned. For this reason, little breeding work has been conducted to improve this plant, also because of the low genetic variability present in this hexaploid species. In this review, we summarize the main works on this crop, focused on genetic improvement with three main objectives: yield, seed oil content and quality, and reduction in glucosinolates content in the seed, which are the main anti-nutritional substances present in camelina. We also report the latest advances in utilising classical plant breeding, transgenic approaches, and CRISPR-Cas9 genome-editing.

Natural variation and artificial selection at the BnaC2.MYB28 locus modulate Brassica napus seed glucosinolate

The degradation products of glucosinolates (GSLs) greatly lower the nutritional value of rapeseed (Brassica napus) meal; thus, reduction of seed GSL content (SGC) has become an important objective of rapeseed breeding. In our previous study, we finely mapped a major QTL (qGSL-C2) for SGC to a 49-kb collinear region on B. rapa chromosome A2. Here, we experimentally validated that BnaC2.MYB28, encoding an R2R3-MYB transcription factor, is the causal gene of qGSL-C2. BnaC2.MYB28 is a nucleus-localized protein mainly expressed in vegetative tissues. Knockout of BnaC2.MYB28 in the high-SGC parent G120 reduced SGC to a value lower than that in the low-SGC parent ZY50, while overexpression of BnaC2.MYB28 in both parental lines (G120 and ZY50) led to extremely high SGC, indicating that BnaC2.MYB28 acts as a positive regulator of SGC in both parents. Molecular characterization revealed that BnaC2.MYB28 forms a homodimer and specifically interacts with BnaMYC3. Moreover, BnaC2.MYB28 can directly activate the expression of GSL biosynthesis genes. Differential expression abundance resulting from the polymorphic promoter sequences, in combination with the different capability in activating downstream genes involved in aliphatic GSL biosynthesis, caused the functional divergence of BnaC2.MYB28 in SGC regulation between the parents. Natural variation of BnaC2.MYB28 was highly associated with SGC in natural germplasm and has undergone artificial selection in modern low-GSL breeding. This study provides important insights into the core function of BnaC2.MYB28 in regulating SGC and a promising strategy for manipulating SGC in rapeseed.

Transcriptomic and Metabolomic Analysis of a Pseudomonas-Resistant versus a Susceptible Arabidopsis Accession

Accessions of one plant species may show significantly different levels of susceptibility to stresses. The Arabidopsis thaliana accessions Col-0 and C24 differ significantly in their resistance to the pathogen Pseudomonas syringae pv. tomato (Pst). To help unravel the underlying mechanisms contributing to this naturally occurring variance in resistance to Pst, we analyzed changes in transcripts and compounds from primary and secondary metabolism of Col-0 and C24 at different time points after infection with Pst. Our results show that the differences in the resistance of Col-0 and C24 mainly involve mechanisms of salicylic-acid-dependent systemic acquired resistance, while responses of jasmonic-acid-dependent mechanisms are shared between the two accessions. In addition, arginine metabolism and differential activity of the biosynthesis pathways of aliphatic glucosinolates and indole glucosinolates may also contribute to the resistance. Thus, this study highlights the difference in the defense response strategies utilized by different genotypes.

Genetic disruption of Arabidopsis secondary metabolite synthesis leads to microbiome-mediated modulation of nematode invasion

In-depth understanding of metabolite-mediated plant-nematode interactions can guide us towards novel nematode management strategies. To improve our understanding of the effects of secondary metabolites on soil nematode communities, we grew Arabidopsis thaliana genetically altered in glucosinolate, camalexin, or flavonoid synthesis pathways, and analyzed their root-associated nematode communities using metabarcoding. To test for any modulating effects of the associated microbiota on the nematode responses, we characterized the bacterial and fungal communities. Finally, as a proxy of microbiome-modulating effects on nematode invasion, we isolated the root-associated microbiomes from the mutants and tested their effect on the ability of the plant parasitic nematode Meloidogyne incognita to penetrate tomato roots. Most mutants had altered relative abundances of several nematode taxa with stronger effects on the plant parasitic Meloidogyne hapla than on other root feeding taxa. This probably reflects that M. hapla invades and remains embedded within root tissues and is thus intimately associated with the host. When transferred to tomato, microbiomes from the flavonoid over-producing pap1-D enhanced M. incognita root-invasion, whereas microbiomes from flavonoid-deficient mutants reduced invasion. This suggests microbiome-mediated effect of flavonoids on Meloidogyne infectivity plausibly mediated by the alteration of the abundances of specific microbial taxa in the transferred microbiomes, although we could not conclusively pinpoint such causative microbial taxa.

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'.

Optimization of Heterologous Glucoraphanin Production In Planta

Glucoraphanin is a plant specialized metabolite found in cruciferous vegetables that has long been a target for production in a heterologous host because it can subsequently be hydrolyzed to form the chemopreventive compound sulforaphane before and during consumption. However, previous studies have only been able to produce small amounts of glucoraphanin in heterologous plant and microbial systems compared to the levels found in glucoraphanin-producing plants, suggesting that there may be missing auxiliary genes that play a role in improving production in planta. In an effort to identify auxiliary genes required for high glucoraphanin production, we leveraged transient expression in Nicotiana benthamiana to screen a combination of previously uncharacterized coexpressed genes and rationally selected genes alongside the glucoraphanin biosynthetic pathway. This strategy alleviated metabolic bottlenecks, which improved glucoraphanin production by 4.74-fold. Our optimized glucoraphanin biosynthetic pathway provides a pathway amenable for high glucoraphanin production.

Metabolic Profiling of Organic Acids Reveals the Involvement of HuIPMS2 in Citramalic Acid Synthesis in Pitaya

Pitayas are rich in organic acids, especially citramalic acid, which is significantly higher than the plants. However, the mechanism of citramalic acid biosynthesis remains to be fully elucidated. In this study, organic acid compositions and contents, as well as expression patterns of key genes related to organic acid metabolism were analyzed during fruit maturation of four different pitaya cultivars i.e., ‘Guanhuabai’ (GHB), ‘Guanhuahong’ (GHH), ‘Wucihuanglong’ (WCHL), and ‘Youcihuanglong’ (YCHL). The total organic acid contents increased first and then declined during fruit maturation. The main organic acids were citramalic acid during the early stages of GHB, GHH, and WCHL pitayas, and dominated by malic acid as fruit maturation. In comparison, citric acid and malic acid were main organic acid for ‘YCHL’ pitaya. Citramalate synthase (IPMS) was involved in the synthesis of citramalic acid, and three types of HuIPMS i.e., HuIPMS1, HuIPMS2, and HuIPMS3, were obtained in our study. Highest expression levels of HuIPMS1 were detected in sepals, while HuIPMS2 and HuIPMS3 exhibited preferential expression in tender stems and ovaries. The expression levels of HuIPMS2 and HuIPMS3 were positively correlated with the content of citramalic acid in the four pitaya cultivars. HuIPMS2 was a chloroplast-localized protein, while HuIPMS3 presented a cytoplasmic-like and nuclear subcellular localization. These findings provide an important basis for further understanding of the molecular mechanism that leads to citramalic acid metabolism during pitaya fruit maturation.

Arabidopsis Glutathione-S-Transferases GSTF11 and GSTU20 Function in Aliphatic Glucosinolate Biosynthesis

Glutathione (GSH) conjugation with intermediates is required for the biosynthesis of glucosinolate (GSL) by serving as a sulfur supply. Glutathione-S-transferases (GSTs) primarily work on GSH conjugation, suggesting their involvement in GSL metabolism. Although several GSTs, including GSTF11 and GSTU20, have been recently postulated to act in GSL biosynthesis, molecular evidence is lacking. Here, we demonstrated that GSTF11 and GSTU20 play non-redundant, although partially overlapping, roles in aliphatic GSL biosynthesis. In addition, GSTU20 plays a more important role than GSTF11, which is manifested by the greater loss of aliphatic GSLs associated with GSTU20 mutant and a greater number of differentially expressed genes in GSTU20 mutant compared to GSTF11 mutant. Moreover, a double mutation leads to a greater aggregate loss of aliphatic GSLs, suggesting that GSTU20 and GSTF11 may function in GSL biosynthesis in a dosage-dependent manner. Together, our results provide direct evidence that GSTU20 and GSTF11 are critically involved in aliphatic GSL biosynthesis, filling the knowledge gap that has been speculated in recent decades.

Analysis of Glucosinolate Content and Metabolism Related Genes in Different Parts of Chinese Flowering Cabbage

Glucosinolates (GSLs) are important secondary metabolites that play important defensive roles in cruciferous plants. Chinese flowering cabbage, one of the most common vegetable crops, is rich in GSLs and thus has the potential to reduce the risk of cancer in humans. Many genes that are involved in GSL biosynthesis and metabolism have been identified in the model plant Arabidopsis thaliana; however, few studies investigated the genes related to GSL biosynthesis and metabolism in Chinese flowering cabbage. In the present study, the GSL composition and content in three different organs of Chinese flowering cabbage (leaf, stalk, and flower bud) were determined. Our results showed that the total GSL content in flower buds was significantly higher than in stalks and leaves, and aliphatic GSLs were the most abundant GSL type. To understand the molecular mechanisms underlying the variations of GSL content, we analyzed the expression of genes encoding enzymes involved in GSL biosynthesis and transport in different tissues of Chinese flowering cabbage using RNA sequencing; the expression levels of most genes were found to be consistent with the pattern of total GSL content. Correlation and consistency analysis of differentially expressed genes from different organs with the GSL content revealed that seven genes (Bra029966, Bra012640, Bra016787, Bra011761, Bra006830, Bra011759, and Bra029248) were positively correlated with GSL content. These findings provide a molecular basis for further elucidating GSL biosynthesis and transport in Chinese flowering cabbage.

Genome- and transcriptome-wide association studies reveal the genetic basis and the breeding history of seed glucosinolate content in Brassica napus

A high content of seed glucosinolates and their degradation products imposes anti-nutritional effects on livestock; therefore, persistent efforts are made to reduce the seed GSL content to increase the commercial value of rapeseed meal. Here, we dissected the genetic structure of SGC by genome-wide association studies (GWAS) combined with transcriptome-wide association studies (TWAS). Fifteen reliable quantitative trait loci (QTLs) were identified to be associated with the reduced SGC in modern B. napus cultivars by GWAS. Analysis of the selection strength and haplotypes at these QTLs revealed that low SGC was predominantly generated by the co-selection of qGSL.A02.2, qGSL.C02.1, qGSL.A09.2, and qGSL.C09.1. Integration of the results from TWAS, comprehensive bioinformatics, and POCKET algorithm analyses indicated that BnaC02.GTR2 (BnaC02g42260D) is a candidate gene underlying qGSL.C02.1. Using CRISPR/Cas9-derived Bna.gtr2s knockout mutants, we experimentally verified that both BnaC02.GTR2 and its three paralogs positively regulate seed GSL accumulation but negatively regulated vegetative tissue GSL contents. In addition, we observed smaller seeds with higher seed oil content in these Bna.gtr2 mutants. Furthermore, both RNA-seq and correlation analyses suggested that Bna.GTR2s might play a comprehensive role in seed development, such as amino acid accumulation, GSL synthesis, sugar assimilation, and oil accumulation. This study unravels the breeding selection history of low-SGC improvement and provides new insights into the molecular function of Bna.GTR2s in both seed GSL accumulation and seed development in B. napus.

Herbivore feeding preference corroborates optimal defense theory for specialized metabolites within plants

Numerous plants protect themselves from attackers by using specialized metabolites. The biosynthesis of these deterrent, often toxic metabolites is costly, as their synthesis diverts energy and resources on account of growth and development. How plants diversify investments into growth and defense is explained by the optimal defense theory. The central prediction of the optimal defense theory is that plants maximize growth and defense by concentrating specialized metabolites in tissues that are decisive for fitness. To date, supporting physiological evidence relies on the correlation between plant metabolite presence and animal feeding preference. Here, we use glucosinolates as a model to examine the effect of changes in chemical defense distribution on feeding preference. Taking advantage of the uniform glucosinolate distribution in transporter mutants, we show that high glucosinolate accumulation in tissues important to fitness protects them by guiding larvae of a generalist herbivore to feed on other tissues. Moreover, we show that the mature leaves of Arabidopsis thaliana supply young leaves with glucosinolates to optimize defense against herbivores. Our study provides physiological evidence for the central hypothesis of the optimal defense theory and sheds light on the importance of integrating glucosinolate biosynthesis and transport for optimizing plant defense.

Effects of Sulfur Assimilation in Pseudomonas fluorescens SS101 on Growth, Defense, and Metabolome of Different Brassicaceae

Genome-wide analysis of plant-growth-promoting Pseudomonas fluorescens strain SS101 (PfSS101) followed by site-directed mutagenesis previously suggested that sulfur assimilation may play an important role in growth promotion and induced systemic resistance in Arabidopsis. Here, we investigated the effects of sulfur metabolism in PfSS101 on growth, defense, and shoot metabolomes of Arabidopsis and the Brassica crop, Broccoli. Root tips of seedlings of Arabidopsis and two Broccoli cultivars were treated with PfSS101 or with a mutant disrupted in the adenylsulfate reductase cysH, a key gene in cysteine and methionine biosynthesis. Phenotyping of plants treated with wild-type PfSS101 or its cysH mutant revealed that sulfur assimilation in PfSS101 was associated with enhanced growth of Arabidopsis but with a reduction in shoot biomass of two Broccoli cultivars. Untargeted metabolomics revealed that cysH-mediated sulfur assimilation in PfSS101 had significant effects on shoot chemistry of Arabidopsis, in particular on chain elongation of aliphatic glucosinolates (GLSs) and on indole metabolites, including camalexin and the growth hormone indole-3-acetic acid. In Broccoli, PfSS101 sulfur assimilation significantly upregulated the relative abundance of several shoot metabolites, in particular, indolic GLSs and phenylpropanoids. These metabolome changes in Broccoli plants coincided with PfSS101-mediated suppression of leaf infections by Xanthomonas campestris. Our study showed the metabolic interconnectedness of plants and their root-associated microbiota.

Structural Studies of Aliphatic Glucosinolate Chain-Elongation Enzymes

Plants evolved specialized metabolic pathways through gene duplication and functional divergence of enzymes involved in primary metabolism. The results of this process are varied pathways that produce an array of natural products useful to both plants and humans. In plants, glucosinolates are a diverse class of natural products. Glucosinolate function stems from their hydrolysis products, which are responsible for the strong flavors of Brassicales plants, such as mustard, and serve as plant defense molecules by repelling insects, fighting fungal infections, and discouraging herbivory. Additionally, certain hydrolysis products such as isothiocyanates can potentially serve as cancer prevention agents in humans. The breadth of glucosinolate function is a result of its great structural diversity, which comes from the use of aliphatic, aromatic and indole amino acids as precursors and elongation of some side chains by up to nine carbons, which, after the formation of the core glucosinolate structure, can undergo further chemical modifications. Aliphatic methionine-derived glucosinolates are the most abundant form of these compounds. Although both elongation and chemical modification of amino acid side chains are important for aliphatic glucosinolate diversity, its elongation process has not been well described at the molecular level. Here, we summarize new insights on the iterative chain-elongation enzymes methylthioalkylmalate synthase (MAMS) and isopropylmalate dehydrogenase (IPMDH).

Comparative transcriptomic analyses of glucosinolate metabolic genes during the formation of Chinese kale seeds

BACKGROUND

To understand the mechanism of glucosinolates (GSs) accumulation in the specific organs, combined analysis of physiological change and transcriptome sequencing were applied in the current study. Taking Chinese kale as material, seeds and silique walls were divided into different stages based on the development of the embryo in seeds and then subjected to GS analysis and transcriptome sequencing.

RESULTS

The main GS in seeds of Chinese kale were glucoiberin and gluconapin and their content changed with the development of the seed. During the transition of the embryo from torpedo- to the early cotyledonary-embryo stage, the accumulation of GS in the seed was accompanied by the salient decline of GS in the corresponding silique wall. Thus, the seed and corresponding silique wall at these two stages were subjected to transcriptomic sequencing analysis. 135 genes related to GS metabolism were identified, of which 24 genes were transcription factors, 81 genes were related to biosynthetic pathway, 25 genes encoded catabolic enzymes, and 5 genes matched with transporters. The expression of GS biosynthetic genes was detected both in seeds and silique walls. The high expression of FMOGS-OX and AOP2, which is related to the production of gluconapin by side modification, was noted in seeds at both stages. Interestingly, the expression of GS biosynthetic genes was higher in the silique wall compared with that in the seed albeit lower content of GS existed in the silique wall than in the seed. Combined with the higher expression of transporter genes GTRs in silique walls than in seeds, it was proposed that the transportation of GS from the silique wall to the seed is an important source for seed GS accumulation. In addition, genes related to GS degradation expressed abundantly in the seed at the early cotyledonary-embryo stage indicating its potential role in balancing seed GS content.

CONCLUSIONS

Two stages including the torpedo-embryo and the early cotyledonary-embryo stage were identified as crucial in GS accumulation during seed development. Moreover, we confirmed the transportation of GS from the silique wall to the seed and proposed possible sidechain modification of GS biosynthesis may exist during seed formation.

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.

Brassica carinata genome characterization clarifies U's triangle model of evolution and polyploidy in Brassica

Ethiopian mustard (Brassica carinata) in the Brassicaceae family possesses many excellent agronomic traits. Here, the high-quality genome sequence of B. carinata is reported. Characterization revealed a genome anchored to 17 chromosomes with a total length of 1.087 Gb and an N50 scaffold length of 60 Mb. Repetitive sequences account for approximately 634 Mb or 58.34% of the B. carinata genome. Notably, 51.91% of 97,149 genes are confined to the terminal 20% of chromosomes as a result of the expansion of repeats in pericentromeric regions. Brassica carinata shares one whole-genome triplication event with the five other species in U's triangle, a classic model of evolution and polyploidy in Brassica. Brassica carinata was deduced to have formed ∼0.047 Mya, which is slightly earlier than B. napus but later than B. juncea. Our analysis indicated that the relationship between the two subgenomes (BcaB and BcaC) is greater than that between other two tetraploid subgenomes (BjuB and BnaC) and their respective diploid parents. RNA-seq datasets and comparative genomic analysis were used to identify several key genes in pathways regulating disease resistance and glucosinolate metabolism. Further analyses revealed that genome triplication and tandem duplication played important roles in the expansion of those genes in Brassica species. With the genome sequencing of B. carinata completed, the genomes of all six Brassica species in U's triangle are now resolved. The data obtained from genome sequencing, transcriptome analysis, and comparative genomic efforts in this study provide valuable insights into the genome evolution of the six Brassica species in U's triangle.

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.

Engineering and optimization of the 2-phenylethylglucosinolate production in Nicotiana benthamiana by combining biosynthetic genes from Barbarea vulgaris and Arabidopsis thaliana

2-Phenylethylglucosinolate (2PE) derived from homophenylalanine is present in plants of the Brassicales order as a defense compound. It is associated with multiple biological properties, including deterrent effects on pests and antimicrobial and health-promoting functions, due to its hydrolysis product 2-phenylethyl isothiocyanate, which confers 2PE as a potential application in agriculture and industry. In this study, we characterized the putative key genes for 2PE biosynthesis from Barbarea vulgaris W.T. Aiton and demonstrated the feasibility of engineering 2PE production in Nicotiana benthamiana Domin. We used different combinations of genes from B. vulgaris and Arabidopsis thaliana (L.) Heynh. to demonstrate that: (i) BvBCAT4 performed more efficiently than AtBCAT4 in biosynthesis of both homophenylalanine and dihomomethionine; (ii) MAM1 enzymes were critical for the chain-elongated profile, while CYP79F enzymes accepted both chain-elongated methionine and homophenylalanine; (iii) aliphatic but not aromatic core structure pathway catalyzed the 2PE biosynthesis; (iv) a chimeric pathway containing BvBCAT4, BvMAM1, AtIPMI and AtIPMDH1 resulted in a two-fold increase in 2PE production compared with the B. vulgaris-specific chain elongation pathway; and (v) profiles of chain-elongated products and glucosinolates partially mirrored the profiles in the gene donor plant, but were wider in N. benthamiana than in the native plants. Our study provides a strategy to produce the important homophenylalanine and 2PE in a heterologous host. Furthermore, chimeric engineering of the complex 2PE biosynthetic pathway enabled detailed understanding of catalytic properties of individual enzymes - a prerequisite for understanding biochemical evolution. The new-to-nature gene combinations have the potential for application in biotechnological and plant breeding.

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.

Protein complex formation in methionine chain-elongation and leucine biosynthesis

During the past two decades, glucosinolate (GLS) metabolic pathways have been under extensive studies because of the importance of the specialized metabolites in plant defense against herbivores and pathogens. The studies have led to a nearly complete characterization of biosynthetic genes in the reference plant Arabidopsis thaliana. Before methionine incorporation into the core structure of aliphatic GLS, it undergoes chain-elongation through an iterative three-step process recruited from leucine biosynthesis. Although enzymes catalyzing each step of the reaction have been characterized, the regulatory mode is largely unknown. In this study, using three independent approaches, yeast two-hybrid (Y2H), coimmunoprecipitation (Co-IP) and bimolecular fluorescence complementation (BiFC), we uncovered the presence of protein complexes consisting of isopropylmalate isomerase (IPMI) and isopropylmalate dehydrogenase (IPMDH). In addition, simultaneous decreases in both IPMI and IPMDH activities in a leuc:ipmdh1 double mutants resulted in aggregated changes of GLS profiles compared to either leuc or ipmdh1 single mutants. Although the biological importance of the formation of IPMI and IPMDH protein complexes has not been documented in any organisms, these complexes may represent a new regulatory mechanism of substrate channeling in GLS and/or leucine biosynthesis. Since genes encoding the two enzymes are widely distributed in eukaryotic and prokaryotic genomes, such complexes may have universal significance in the regulation of leucine biosynthesis.

Citramalate synthase yields a biosynthetic pathway for isoleucine and straight- and branched-chain ester formation in ripening apple fruit

A plant pathway that initiates with the formation of citramalate from pyruvate and acetyl-CoA by citramalate synthase (CMS) is shown to contribute to the synthesis of α-ketoacids and important odor-active esters in apple (Malus × domestica) fruit. Microarray screening led to the discovery of a gene with high amino acid similarity to 2-isopropylmalate synthase (IPMS). However, functional analysis of recombinant protein revealed its substrate preference differed substantially from IPMS and was more typical of CMS. MdCMS also lacked the regulatory region present in MdIPMS and was not sensitive to feedback inhibition. 13C-acetate feeding of apple tissue labeled citramalate and α-ketoacids in a manner consistent with the presence of the citramalate pathway, labeling both straight- and branched-chain esters. Analysis of genomic DNA (gDNA) revealed the presence of two nearly identical alleles in "Jonagold" fruit (MdCMS_1 and MdCMS_2), differing by two nonsynonymous single-nucleotide polymorphisms (SNPs). The mature proteins differed only at amino acid 387, possessing either glutamine387 (MdCMS_1) or glutamate387 (MdCMS_2). Glutamate387 was associated with near complete loss of activity. MdCMS expression was fruit-specific, increasing severalfold during ripening. The translated protein product was detected in ripe fruit. Transient expression of MdCMS_1 in Nicotiana benthamiana induced the accumulation of high levels of citramalate, whereas MdCMS_2 did not. Domesticated apple lines with MdCMS isozymes containing only glutamate387 produced a very low proportion of 2-methylbutanol- and 2-methylbutanoate (2MB) and 1-propanol and propanoate (PROP) esters. The citramalate pathway, previously only described in microorganisms, is shown to function in ripening apple and contribute to isoleucine and 2MB and PROP ester biosynthesis without feedback regulation.

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.

Effect of selenium-sulfur interaction on the anabolism of sulforaphane in broccoli

The effects of S (as sulphate) and Se (as selenite) treatment (S mM/Se μM: 1/0, 1/50, 1/100, 1/150, 4/0, 4/50, 4/100, and 4/150) on the production of sulforaphane (an anticancer compound), the accumulation of its precursor substance, and the expression of genes related to glucoraphanin biosynthesis in broccoli were examined. Sulforaphane yield and myrosinase activity increased significantly with the combined application of 4 mM S and 100 μM Se on broccoli. Furthermore, the concentrations of glucoraphanin (a sulforaphane precursor) and methionine (a glucoraphanin substrate) slightly changed after Se application. And the strong anticancer activity of compound Se-SMC was further improved. Analysis of related gene expression showed that MY, which encodes myrosinase, was strongly induced by Se treatment. Thus, the myrosinase activity induced by Se treatment is the dominant factor affecting sulforaphane yield from glucoraphanin hydrolyzation. Selenium-sulfur biofortification provides a technical support for the cultivation of broccoli with high sulforaphane and high anti-cancer selenium compounds.

Multiomic Analysis Elucidates the Reasons Underlying the Differential Metabolite Accumulation in Citrus Mature Leaves and Fruit Juice Sacs

Fruit and leaf possess distinctly different metabolites. Here, metabolites and transcriptome were compared between mature leaves (ML) and juice sacs (JS) of Citrus grandis "Hirado Buntan" to investigate the possible reasons. Results indicated that the remarkable difference in starch, total flavonoids and carotenoids, l-ascorbate, and jasmonic acid between ML and JS was tightly related to the expression levels of their biosynthesis-related genes, while the significant difference in abscisic acid and citrate was mainly related to the gene expression level(s) of 9-cis-epoxycarotenoid dioxygenase and proton pump genes, respectively. In addition, ATP citrate lyase probably plays a key role in the levels of flavonoids between ML and JS via providing different levels of acetyl-CoA. Taken together, these results identified some key candidate genes responsible for the content of a given metabolite and will contribute to research in regulating such metabolite content in citrus fruits.

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.

In Arabidopsis thaliana Substrate Recognition and Tissue- as Well as Plastid Type-Specific Expression Define the Roles of Distinct Small Subunits of Isopropylmalate Isomerase

In Arabidopsis thaliana, the heterodimeric isopropylmalate isomerase (IPMI) is composed of a single large (IPMI LSU1) and one of three different small subunits (IPMI SSU1 to 3). The function of IPMI is defined by the small subunits. IPMI SSU1 is required for Leu biosynthesis and has previously also been proposed to be involved in the first cycle of Met chain elongation, the first phase of the synthesis of Met-derived glucosinolates. IPMI SSU2 and IPMI SSU3 participate in the Met chain elongation pathway. Here, we investigate the role of the three IPMI SSUs through the analysis of the role of the substrate recognition region spanning five amino acids on the substrate specificity of IPMI SSU1. Furthermore, we analyze in detail the expression pattern of fluorophore-tagged IPMI SSUs throughout plant development. Our study shows that the substrate recognition region that differs between IPMI SSU1 and the other two IMPI SSUs determines the substrate preference of IPMI. Expression of IPMI SSU1 is spatially separated from the expression of IPMI SSU2 and IPMI SSU3, and IPMI SSU1 is found in small plastids, whereas IMPI SSU2 and SSU3 are found in chloroplasts. Our data show a distinct role for IMPI SSU1 in Leu biosynthesis and for IMPI SSU2 and SSU3 in the Met chain elongation pathway.

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.

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.

Quantitative proteomics reveals that tea tree oil effects Botrytis cinerea mitochondria function

Tea tree oil (TTO) inhibits the spore germination and mycelial growth of Botrytis cinerea, and induces mitochondrial dysfunction of B. cinerea. To further determine the effects of TTO on mitochondria in B. cinerea, label-free quantitative proteomics analysis was performed. A total of 85 differentially expression proteins (DEPs) were identified; Among them 51 were more abundant in TTO-treated samples, and 34 were less abundant. DEPs were then annotated and classified into 34 functional groups based on Gene Ontology analysis. Subsequent Kyoto Encyclopedia of Genes and Genomes analysis linked identified DEPs to 83 different pathways. This study suggests that TTO inhibits the tricarboxylic acid cycle, pyruvate metabolism, amino acid metabolism, and membrane-related pathways in mitochondria, and also promotes sphingolipid metabolism, which may accelerate cell death in B. cinerea.

Genomic Origin and Diversification of the Glucosinolate MAM Locus

Glucosinolates are a diverse group of plant metabolites that characterize the order Brassicales. The MAM locus is one of the most significant QTLs for glucosinolate diversity. However, most of what we understand about evolution at the locus is focused on only a few species and not within a phylogenetic context. In this study, we utilize a micro-synteny network and phylogenetic inference to investigate the origin and diversification of the MAM/IPMS gene family. We uncover unique MAM-like genes found at the orthologous locus in the Cleomaceae that shed light on the transition from IPMS to MAM. In the Brassicaceae, we identify six distinct MAM clades across Lineages I, II, and III. We characterize the evolutionary impact and consequences of local duplications, transpositions, whole genome duplications, and gene fusion events, generating several new hypothesizes on the function and diversity of the MAM locus.

Overexpression of HMG-CoA synthase promotes Arabidopsis root growth and adversely affects glucosinolate biosynthesis

3-Hydroxy-3-methylglutaryl-CoA synthase (HMGS) catalyses the second step of the mevalonate (MVA) pathway. An HMGS inhibitor (F-244) has been reported to retard growth in wheat, tobacco, and Brassica juncea, but the mechanism remains unknown. Although the effects of HMGS on downstream isoprenoid metabolites have been extensively reported, not much is known on how it might affect non-isoprenoid metabolic pathways. Here, the mechanism of F-244-mediated inhibition of primary root growth in Arabidopsis and the relationship between HMGS and non-isoprenoid metabolic pathways were investigated by untargeted SWATH-MS quantitative proteomics, quantitative real-time PCR, and target metabolite analysis. Our results revealed that the inhibition of primary root growth caused by F-244 was a consequence of reduced stigmasterol, auxin, and cytokinin levels. Interestingly, proteomic analyses identified a relationship between HMGS and glucosinolate biosynthesis. Inhibition of HMGS activated glucosinolate biosynthesis, resulting from the induction of glucosinolate biosynthesis-related genes, suppression of sterol biosynthesis-related genes, and reduction in sterol levels. In contrast, HMGS overexpression inhibited glucosinolate biosynthesis, due to down-regulation of glucosinolate biosynthesis-related genes, up-regulation of sterol biosynthesis-related genes, and increase in sterol content. Thus, HMGS might represent a target for the manipulation of glucosinolate biosynthesis, given the regulatory relationship between HMGS in the MVA pathway and glucosinolate biosynthesis.

Genetic architecture of glucosinolate variation in Brassica napus

The diverse biological activities of glucosinolate (GSL) hydrolysis products play significant biological and economical roles in the defense system and nutritional qualities of Brassica napus (oilseed rape). Yet, genomic-based study of the B. napus GSL regulatory mechanisms are scarce due to the complexity of working with polyploid species. To address these challenges, we used transcriptome-based GWAS approach, Associative Transcriptomics (AT), across a diversity panel of 288 B. napus genotypes to uncover the underlying genetic basis controlling quantitative variation of GSLs in B. napus vegetative tissues. Single nucleotide polymorphism (SNP) markers and gene expression markers (GEMs) associations identify orthologues of MYB28/HAG1 (AT5G61420), specifically the copies on chromosome A9 and C2, to be the key regulators of aliphatic GSL variation in leaves. We show that the positive correlation observed between aliphatic GSLs in seed and leaf is due to the amount synthesized, as controlled by Bna.HAG1.A9 and Bna.HAG1.C2, rather than by variation in the transport processes. In addition, AT and differential expression analysis in root tissues implicate an orthologue of MYB29/HAG3 (AT5G07690), Bna.HAG3.A3, as controlling root aromatic GSL variation. Based on the root expression data we also propose Bna.MAM3.A3 to have a role in controlling phenylalanine chain elongation for aromatic GSL biosynthesis. This work uncovers a regulator of homophenylalanine-derived aromatic GSLs and implicates the shared biosynthetic pathways between aliphatic and aromatic GSLs.