The small subunit 1 of the Arabidopsis isopropylmalate isomerase is required for normal growth and development and the early stages of glucosinolate formation

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.

Proteomic and metabolomic analysis of Nicotiana benthamiana under dark stress

Exposure to extended periods of darkness is a common source of abiotic stress that significantly affects plant growth and development. To understand how Nicotiana benthamiana responds to dark stress, the proteomes and metabolomes of leaves treated with darkness were studied. In total, 5763 proteins and 165 primary metabolites were identified following dark treatment. Additionally, the expression of autophagy-related gene (ATG) proteins was transiently upregulated. Weighted gene coexpression network analysis (WGCNA) was utilized to find the protein modules associated with the response to dark stress. A total of four coexpression modules were obtained. The results indicated that heat-shock protein (HSP70), SnRK1-interacting protein 1, 2A phosphatase-associated protein of 46 kDa (Tap46), and glutamate dehydrogenase (GDH) might play crucial roles in N. benthamiana's response to dark stress. Furthermore, a protein-protein interaction (PPI) network was constructed and top-degreed proteins were predicted to identify potential key factors in the response to dark stress. These proteins include isopropylmalate isomerase (IPMI), eukaryotic elongation factor 5A (ELF5A), and ribosomal protein 5A (RPS5A). Finally, metabolic analysis suggested that some amino acids and sugars were involved in the dark-responsive pathways. Thus, these results provide a new avenue for understanding the defensive mechanism against dark stress at the protein and metabolic levels in N. benthamiana.

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.

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.

Homeostasis of branched-chain amino acids is critical for the activity of TOR signaling in Arabidopsis

The target of rapamycin (TOR) kinase is an evolutionarily conserved hub of nutrient sensing and metabolic signaling. In plants, a functional connection of TOR activation with glucose availability was demonstrated, while it is yet unclear whether branched-chain amino acids (BCAAs) are a primary input of TOR signaling as they are in yeast and mammalian cells. Here, we report on the characterization of an Arabidopsis mutant over-accumulating BCAAs. Through chemical interventions targeting TOR and by examining mutants of BCAA biosynthesis and TOR signaling, we found that BCAA over-accumulation leads to up-regulation of TOR activity, which causes reorganization of the actin cytoskeleton and actin-associated endomembranes. Finally, we show that activation of TOR is concomitant with alteration of cell expansion, proliferation and specialized metabolism, leading to pleiotropic effects on plant growth and development. These results demonstrate that BCAAs contribute to plant TOR activation and reveal previously uncharted downstream subcellular processes of TOR signaling.

Exploring the basis of 2-propenyl and 3-butenyl glucosinolate synthesis by QTL mapping and RNA-sequencing in Brassica juncea

Brassica juncea is used as a condiment, as vegetables and as an oilseed crop, especially in semiarid areas. In the present study, we constructed a genetic map using one recombinant inbred line (RIL) of B. juncea. A total of 304 ILP (intron length polymorphism) markers were mapped to 18 linkage groups designated LG01-LG18 in B. juncea. The constructed map covered a total genetic length of 1671.13 cM with an average marker interval of 5.50 cM. The QTLs for 2-propenyl glucosinolates (GSLs) colocalized with the QTLs for 3-butenyl GSLs between At1g26180 and BnapPIP1580 on LG08 in the field experiments of 2016 and 2017. These QTLs accounted for an average of 42.3% and 42.6% phenotypic variation for 2-propenyl and 3-butenyl GSLs, respectively. Furthermore, the Illumina RNA-sequencing technique was used to excavate the genes responsible for the synthesis of GSLs in the siliques of the parental lines of the RIL mapping population, because the bulk of the seed GSLs might originate from the siliques. Comparative analysis and annotation by gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) revealed that 324 genes were involved in GSL metabolism, among which only 24 transcripts were differentially expressed genes (DEGs). Among those DEGs, 15 genes were involved in the biosynthesis and transport of aliphatic GSLs, and their expression patterns were further validated by qRT-PCR analysis. Joint QTL mapping and RNA-sequencing analyses reveal one candidate gene of IIL1 (LOC106416451) for GSL metabolism in B. juncea. These results will be helpful for further fine mapping, gene cloning and genetic mechanisms of 2-propenyl and 3-butenyl GSLs in B. juncea.

Changing substrate specificity and iteration of amino acid chain elongation in glucosinolate biosynthesis through targeted mutagenesis of Arabidopsis methylthioalkylmalate synthase 1

Methylthioalkylmalate synthases catalyse the committing step of amino acid chain elongation in glucosinolate biosynthesis. As such, this group of enzymes plays an important role in determining the glucosinolate composition of Brassicaceae species, including Arabidopsis thaliana. Based on protein structure modelling of MAM1 from A. thaliana and analysis of 57 MAM sequences from Brassicaceae species, we identified four polymorphic residues likely to interact with the 2-oxo acid substrate. Through site-directed mutagenesis, the natural variation in these residues and the effect on product composition were investigated. Fifteen MAM1 variants as well as the native MAM1 and MAM3 from A. thaliana were characterised by heterologous expression of the glucosinolate chain elongation pathway in Escherichia coli. Detected products derived from leucine, methionine or phenylalanine were elongated with up to six methylene groups. Product profile and accumulation were changed in 14 of the variants, demonstrating the relevance of the identified residues. The majority of the single amino acid substitutions decreased the length of methionine-derived products, while approximately half of the substitutions increased the phenylalanine-derived products. Combining two substitutions enabled the MAM1 variant to increase the number of elongation rounds of methionine from three to four. Notably, characterisation of the native MAMs indicated that MAM1 and not MAM3 is responsible for homophenylalanine production. This hypothesis was confirmed by glucosinolate analysis in mam1 and mam3 mutants of A. thaliana.

Molecular Basis of the Evolution of Methylthioalkylmalate Synthase and the Diversity of Methionine-Derived Glucosinolates

The globally cultivated Brassica species possess diverse aliphatic glucosinolates, which are important for plant defense and animal nutrition. The committed step in the side chain elongation of methionine-derived aliphatic glucosinolates is catalyzed by methylthioalkylmalate synthase, which likely evolved from the isopropylmalate synthases of leucine biosynthesis. However, the molecular basis for the evolution of methylthioalkylmalate synthase and its generation of natural product diversity in Brassica is poorly understood. Here, we show that Brassica genomes encode multiple methylthioalkylmalate synthases that have differences in expression profiles and 2-oxo substrate preferences, which account for the diversity of aliphatic glucosinolates across Brassica accessions. Analysis of the 2.1 Å resolution x-ray crystal structure of Brassica juncea methylthioalkylmalate synthase identified key active site residues responsible for controlling the specificity for different 2-oxo substrates and the determinants of side chain length in aliphatic glucosinolates. Overall, these results provide the evolutionary and biochemical foundation for the diversification of glucosinolate profiles across globally cultivated Brassica species, which could be used with ongoing breeding strategies toward the manipulation of beneficial glucosinolate compounds for animal health and plant protection.

ARSENATE INDUCED CHLOROSIS 1/ TRANSLOCON AT THE OUTER ENVOLOPE MEMBRANE OF CHLOROPLASTS 132 Protects Chloroplasts from Arsenic Toxicity

Arsenic (As) is highly toxic to plants and detoxified primarily through complexation with phytochelatins (PCs) and other thiol compounds. To understand the mechanisms of As toxicity and detoxification beyond PCs, we isolated an arsenate-sensitive mutant of Arabidopsis (Arabidopsis thaliana), arsenate induced chlorosis1 (aic1), in the background of the PC synthase-defective mutant cadmium-sensitive1-3 (cad1-3). Under arsenate stress, aic1 cad1-3 showed larger decreases in chlorophyll content and the number and size of chloroplasts than cad1-3 and a severely distorted chloroplast structure. The aic1 single mutant also was more sensitive to arsenate than the wild type (Columbia-0). As concentrations in the roots, shoots, and chloroplasts were similar between aic1 cad1-3 and cad1-3 Using genome resequencing and complementation, TRANSLOCON AT THE OUTER ENVOLOPE MEMBRANE OF CHLOROPLAST132 (TOC132) was identified as the mutant gene, which encodes a translocon protein involved in the import of preproteins from the cytoplasm into the chloroplasts. Proteomic analysis showed that the proteome of aic1 cad1-3 chloroplasts was more affected by arsenate stress than that of cad1-3 A number of proteins related to chloroplast ribosomes, photosynthesis, compound synthesis, and thioredoxin systems were less abundant in aic1 cad1-3 than in cad1-3 under arsenate stress. Our results indicate that chloroplasts are a sensitive target of As toxicity and that AIC1/Toc132 plays an important role in protecting chloroplasts from As toxicity.

Volatile compounds emitted by diverse phytopathogenic microorganisms promote plant growth and flowering through cytokinin action

It is known that volatile emissions from some beneficial rhizosphere microorganisms promote plant growth. Here we show that volatile compounds (VCs) emitted by phylogenetically diverse rhizosphere and non-rhizhosphere bacteria and fungi (including plant pathogens and microbes that do not normally interact mutualistically with plants) promote growth and flowering of various plant species, including crops. In Arabidopsis plants exposed to VCs emitted by the phytopathogen Alternaria alternata, changes included enhancement of photosynthesis and accumulation of high levels of cytokinins (CKs) and sugars. Evidence obtained using transgenic Arabidopsis plants with altered CK status show that CKs play essential roles in this phenomenon, because growth and flowering responses to the VCs were reduced in mutants with CK-deficiency (35S:AtCKX1) or low receptor sensitivity (ahk2/3). Further, we demonstrate that the plant responses to fungal VCs are light-dependent. Transcriptomic analyses of Arabidopsis leaves exposed to A. alternata VCs revealed changes in the expression of light- and CK-responsive genes involved in photosynthesis, growth and flowering. Notably, many genes differentially expressed in plants treated with fungal VCs were also differentially expressed in plants exposed to VCs emitted by the plant growth promoting rhizobacterium Bacillus subtilis GB03, suggesting that plants react to microbial VCs through highly conserved regulatory mechanisms.

Arabidopsis Responds to Alternaria alternata Volatiles by Triggering Plastid Phosphoglucose Isomerase-Independent Mechanisms

Volatile compounds (VCs) emitted by phylogenetically diverse microorganisms (including plant pathogens and microbes that do not normally interact mutualistically with plants) promote photosynthesis, growth, and the accumulation of high levels of starch in leaves through cytokinin (CK)-regulated processes. In Arabidopsis (Arabidopsis thaliana) plants not exposed to VCs, plastidic phosphoglucose isomerase (pPGI) acts as an important determinant of photosynthesis and growth, likely as a consequence of its involvement in the synthesis of plastidic CKs in roots. Moreover, this enzyme plays an important role in connecting the Calvin-Benson cycle with the starch biosynthetic pathway in leaves. To elucidate the mechanisms involved in the responses of plants to microbial VCs and to investigate the extent of pPGI involvement, we characterized pPGI-null pgi1-2 Arabidopsis plants cultured in the presence or absence of VCs emitted by Alternaria alternata We found that volatile emissions from this fungal phytopathogen promote growth, photosynthesis, and the accumulation of plastidic CKs in pgi1-2 leaves. Notably, the mesophyll cells of pgi1-2 leaves accumulated exceptionally high levels of starch following VC exposure. Proteomic analyses revealed that VCs promote global changes in the expression of proteins involved in photosynthesis, starch metabolism, and growth that can account for the observed responses in pgi1-2 plants. The overall data show that Arabidopsis plants can respond to VCs emitted by phytopathogenic microorganisms by triggering pPGI-independent mechanisms.

Optimization of Engineered Production of the Glucoraphanin Precursor Dihomomethionine in Nicotiana benthamiana

Glucosinolates are natural products characteristic of the Brassicales order, which include vegetables such as cabbages and the model plant Arabidopsis thaliana. Glucoraphanin is the major glucosinolate in broccoli and associated with the health-promoting effects of broccoli consumption. Toward our goal of creating a rich source of glucoraphanin for dietary supplements, we have previously reported the feasibility of engineering glucoraphanin in Nicotiana benthamiana through transient expression of glucoraphanin biosynthetic genes from A. thaliana (Mikkelsen et al., 2010). As side-products, we obtained fivefold to eightfold higher levels of chain-elongated leucine-derived glucosinolates, not found in the native plant. Here, we investigated two different strategies to improve engineering of the methionine chain elongation part of the glucoraphanin pathway in N. benthamiana: (1) coexpression of the large subunit (LSU1) of the heterodimeric isopropylmalate isomerase and (2) coexpression of BAT5 transporter for efficient transfer of intermediates across the chloroplast membrane. We succeeded in raising dihomomethionine (DHM) levels to a maximum of 432 nmol g⁻¹ fresh weight that is equivalent to a ninefold increase compared to the highest production of this intermediate, as previously reported (Mikkelsen et al., 2010). The increased DHM production without increasing leucine-derived side-product levels provides new metabolic engineering strategies for improved glucoraphanin production in a heterologous host.

The cytosolic branched-chain aminotransferases of Arabidopsis thaliana influence methionine supply, salvage and glucosinolate metabolism

Arabidopsis thaliana possesses six branched-chain aminotransferases (BCAT1-6). Previous studies revealed that some members of this protein family are involved in the biosynthesis of branched-chain amino acids and/or in the Met chain elongation pathway, the initial steps towards the biosynthesis of Met-derived glucosinolates. We now analyzed branched-chain aminotransferase 6 (BCAT6). In vivo GFP-tagging experiments strongly suggest this enzyme to be localized to the cytosol. Substrate specificity assays performed with recombinant enzyme revealed that BCAT6 transaminates Val, Leu and Ile as well as the corresponding 2-oxo acids but also transaminates Met and its cognate ketoacid 4-methyl-2-oxobutanoate. We established single (bcat6-1), double (bcat4-2/bcat6-1) and triple (bcat3-1/bcat4-2/bcat6-1) mutants involving BCAT6 with the latter exhibiting a clear macroscopic phenotype with smaller plants and abnormal leaf morphology. Metabolite profiling of these mutants demonstrated that BCAT6 can contribute to Met chain elongation with the triple mutant line lacking BCAT3, 4 and 6 showing a dramatic reduction of Met-derived glucosinolate species down to 32 and 14% of wild-type levels in plant foliage and seeds, respectively. This drop in glucosinolate levels is accompanied by a 46-fold increase of free Met, demonstrating the important role of the three branched-chain aminotransferases in converting Met to its 2-oxo acid for glucosinolate chain elongation. In addition, we determined the relative amounts of 5'-deoxy-5'-methylthioadenosine, an intermediate of the Met recycling pathway. This metabolite accumulated to relative high amounts in the absence of the cytosolic BCAT4 and BCAT6, suggesting that cytosolic Met salvage also contributes to the biosynthesis of glucosinolates.

Distinct patterns of the histone marks associated with recruitment of the methionine chain-elongation pathway from leucine biosynthesis

Aliphatic glucosinolates (GLSs) are derived from chain-elongated methionine produced by an iterative three-step process, known to be evolutionarily recruited from leucine biosynthesis. The divergence of homologous genes between two pathways is mainly linked to the alterations in biochemical features. In this study, it was discovered that a distinct pattern of histone modifications is associated with and/or contributes to the divergence of the two pathways. In general, genes involved in leucine biosynthesis were robustly associated with H3k4me2 and H3K4me3. In contrast, despite the considerable abundances of H3K4me2 observed in some of genes involved in methionine chain elongation, H3K4me3 was completely missing. This H3K4m3-depleted pattern had no effect on gene transcription, whereas it seemingly co-evolved with the entire pathway of aliphatic GLS biosynthesis. The results reveal a novel association of the epigenetic marks with plant secondary metabolism, and may help to understand the recruitment of the methionine chain-elongation pathway from leucine biosynthesis.