Omics-based approaches to methionine side chain elongation in Arabidopsis: characterization of the genes encoding methylthioalkylmalate isomerase and methylthioalkylmalate dehydrogenase

Kinetic and Catalytic Mechanism 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 (4-MTOB) and acetyl-CoA generates a 2-malate derivative that either enters the pathway for 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 (IPMS) but left details of the reaction mechanism unanswered. Here we use the B. juncea MAMS2A (BjMAMS2A) isoform to analyze the kinetic and catalytic mechanism 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-MTOB 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 IPMS from leucine biosynthesis in plants and microbes.

Supplementation of Saussurea costus root alleviates sodium nitrite-induced hepatorenal toxicity by modulating metabolic profile, inflammation, and apoptosis

Sodium nitrite (NaNO2) is a widely used food ingredient, although excessive concentrations can pose potential health risks. In the present study, we evaluated the deterioration effects of NaNO2 additives on hematology, metabolic profile, liver function, and kidney function of male Wistar rats. We further explored the therapeutic potential of supplementation with S. costus root ethanolic extract (SCREE) to improve NaNO2-induced hepatorenal toxicity. In this regard, 65 adult male rats were divided into eight groups; Group 1: control, Groups 2, 3, and 4 received SCREE in 200, 400, and 600 mg/kg body weight, respectively, Group 5: NaNO2 (6.5 mg/kg body weight), Groups 6, 7 and 8 received NaNO2 (6.5 mg/kg body weight) in combination with SCREE (200, 400, and 600 mg/kg body weight), respectively. Our results revealed that the NaNO2-treated group shows a significant change in deterioration in body and organ weights, hematological parameters, lipid profile, and hepatorenal dysfunction, as well as immunohistochemical and histopathological alterations. Furthermore, the NaNO2-treated group demonstrated a considerable increase in the expression of TNF-α cytokine and tumor suppressor gene P53 in the kidney and liver, while a significant reduction was detected in the anti-inflammatory cytokine IL-4 and the apoptosis suppressor gene BCL-2, compared to the control group. Interestingly, SCREE administration demonstrated the ability to significantly alleviate the toxic effects of NaNO2 and improve liver function in a dose-dependent manner, including hematological parameters, lipid profile, and modulation of histopathological architecture. Additionally, SCREE exhibited the ability to modulate the expression levels of inflammatory cytokines and apoptotic genes in the liver and kidney. The phytochemical analysis revealed a wide set of primary metabolites in SCREE, including phenolics, flavonoids, vitamins, alkaloids, saponins and tannins, while the untargeted UPLC/T-TOF-MS/MS analysis identified 183 metabolites in both positive and negative ionization modes. Together, our findings establish the potential of SCREE in mitigating the toxic effects of NaNO2 by modulating metabolic, inflammatory, and apoptosis. Together, this study underscores the promise of SCREE as a potential natural food detoxifying additive to counteract the harmful impacts of sodium nitrite.

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.

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.

ATTED-II v11: A Plant Gene Coexpression Database Using a Sample Balancing Technique by Subagging of Principal Components

ATTED-II (https://atted.jp) is a gene coexpression database for nine plant species based on publicly available RNAseq and microarray data. One of the challenges in constructing condition-independent coexpression data based on publicly available gene expression data is managing the inherent sampling bias. Here, we report ATTED-II version 11, wherein we adopted a coexpression calculation methodology to balance the samples using principal component analysis and ensemble calculation. This approach has two advantages. First, omitting principal components with low contribution rates reduces the main contributors of noise. Second, balancing large differences in contribution rates enables considering various sample conditions entirely. In addition, based on RNAseq- and microarray-based coexpression data, we provide species-representative, integrated coexpression information to enhance the efficiency of interspecies comparison of the coexpression data. These coexpression data are provided as a standardized z-score to facilitate integrated analysis with different data sources. We believe that with these improvements, ATTED-II is more valuable and powerful for supporting interspecies comparative studies and integrated analyses using heterogeneous data.

Genome-Wide Association Study of Phenylalanine Derived Glucosinolates in Brassica rapa

Glucosinolates (GSLs) are sulfur-containing bioactive compounds usually present in Brassicaceae plants and are usually responsible for a pungent flavor and reduction of the nutritional values of seeds. Therefore, breeding rapeseed varieties with low GSL levels is an important breeding objective. Most GSLs in Brassica rapa are derived from methionine or tryptophan, but two are derived from phenylalanine, one directly (benzylGSL) and one after a round of chain elongation (phenethylGSL). In the present study, two phenylalanine (Phe)-derived GSLs (benzylGSL and phenethylGSL) were identified and quantified in seeds by liquid chromatography and mass spectrometry (LC-MS) analysis. Levels of benzylGSL were low but differed among investigated low and high GSL genotypes. Levels of phenethylGSL (also known as 2-phenylethylGSL) were high but did not differ among GSL genotypes. Subsequently, a genome-wide association study (GWAS) was conducted using 159 B. rapa accessions to demarcate candidate regions underlying 43 and 59 QTNs associated with benzylGSL and phenethylGSL that were distributed on 10 chromosomes and 9 scaffolds, explaining 0.56% to 70.86% of phenotypic variations, respectively. Furthermore, we find that 15 and 18 known or novel candidate genes were identified for the biosynthesis of benzylGSL and phenethylGSL, including known regulators of GSL biosynthesis, such as BrMYB34, BrMYB51, BrMYB28, BrMYB29 and BrMYB122, and novel regulators or structural genes, such as BrMYB44/BrMYB77 and BrMYB60 for benzylGSL and BrCYP79B2 for phenethylGSL. Finally, we investigate the expression profiles of the biosynthetic genes for two Phe-derived GSLs by transcriptomic analysis. Our findings provide new insight into the complex machinery of Phe-derived GSLs in seeds of B. rapa and help to improve the quality of Brassicaceae plant breeding.

Identification of Potential Genes Encoding Protein Transporters in Arabidopsis thaliana Glucosinolate (GSL) Metabolism

Several species in Brassicaceae produce glucosinolates (GSLs) to protect themselves against pests. As demonstrated in A. thaliana, the reallocation of defence compounds, of which GSLs are a major part, is highly dependent on transport processes and serves to protect high-value tissues such as reproductive tissues. This study aimed to identify potential GSL-transporter proteins (TPs) using a network-biology approach. The known A. thaliana GSL genes were retrieved from the literature and pathway databases and searched against several co-expression databases to generate a gene network consisting of 1267 nodes and 14,308 edges. In addition, 1151 co-expressed genes were annotated, integrated, and visualised using relevant bioinformatic tools. Based on three criteria, 21 potential GSL genes encoding TPs were selected. The AST68 and ABCG40 potential GSL TPs were chosen for further investigation because their subcellular localisation is similar to that of known GSL TPs (SULTR1;1 and SULTR1;2) and ABCG36, respectively. However, AST68 was selected for a molecular-docking analysis using AutoDOCK Vina and AutoDOCK 4.2 with the generated 3D model, showing that both domains were well superimposed on the homologs. Both molecular-docking tools calculated good binding-energy values between the sulphate ion and Ser419 and Val172, with the formation of hydrogen bonds and van der Waals interactions, respectively, suggesting that AST68 was one of the sulphate transporters involved in GSL biosynthesis. This finding illustrates the ability to use computational analysis on gene co-expression data to screen and characterise plant TPs on a large scale to comprehensively elucidate GSL metabolism in A. thaliana. Most importantly, newly identified potential GSL transporters can serve as molecular tools in improving the nutritional value of crops.

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

Potential Arabidopsis thaliana glucosinolate genes identified from the co-expression modules using graph clustering approach

Background

Glucosinolates (GSLs) are plant secondary metabolites that contain nitrogen-containing compounds. They are important in the plant defense system and known to provide protection against cancer in humans. Currently, increasing the amount of data generated from various omics technologies serves as a hotspot for new gene discovery. However, sometimes sequence similarity searching approach is not sufficiently effective to find these genes; hence, we adapted a network clustering approach to search for potential GSLs genes from the Arabidopsis thaliana co-expression dataset.

Methods

We used known GSL genes to construct a comprehensive GSL co-expression network. This network was analyzed with the DPClusOST algorithm using a density of 0.5. 0.6. 0.7, 0.8, and 0.9. Generating clusters were evaluated using Fisher’s exact test to identify GSL gene co-expression clusters. A significance score (SScore) was calculated for each gene based on the generated p-value of Fisher’s exact test. SScore was used to perform a receiver operating characteristic (ROC) study to classify possible GSL genes using the ROCR package. ROCR was used in determining the AUC that measured the suitable density value of the cluster for further analysis. Finally, pathway enrichment analysis was conducted using ClueGO to identify significant pathways associated with the GSL clusters.

Results

The density value of 0.8 showed the highest area under the curve (AUC) leading to the selection of thirteen potential GSL genes from the top six significant clusters that include IMDH3, MVP1, T19K24.17, MRSA2, SIR, ASP4, MTO1, At1g21440, HMT3, At3g47420, PS1, SAL1, and At3g14220. A total of Four potential genes (MTO1, SIR, SAL1, and IMDH3) were identified from the pathway enrichment analysis on the significant clusters. These genes are directly related to GSL-associated pathways such as sulfur metabolism and valine, leucine, and isoleucine biosynthesis. This approach demonstrates the ability of the network clustering approach in identifying potential GSL genes which cannot be found from the standard similarity search.

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.

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

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

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.

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.

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.

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.

[Biosynthesis of indole glucosinolates and ecological role of secondary modification pathways]

Indole glucosinolates are plant secondary metabolites derived from the amino acid tryptophan. They are part of a large group of sulfur-containing molecules almost exclusively found among Brassicales, which include the mustard family (Brassicaceae) with many edible plant species of major nutritional importance. These compounds mediate numerous interactions between these plants and their natural enemies and are therefore of major biological and economical interest. This literature review aims at taking stock of recent advances of our knowledge about the biosynthetic pathways of indole glucosinolates, but also about the defense strategies and ecological processes involving these metabolites.

[Development of Plant Metabolomics and Medicinal Plant Genomics]

 A variety of chemicals produced by plants, often referred to as 'phytochemicals', have been used as medicines, food, fuels and industrial raw materials. Recent advances in the study of genomics and metabolomics in plant science have accelerated our understanding of the mechanisms, regulation and evolution of the biosynthesis of specialized plant products. We can now address such questions as how the metabolomic diversity of plants is originated at the levels of genome, and how we should apply this knowledge to drug discovery, industry and agriculture. Our research group has focused on metabolomics-based functional genomics over the last 15 years and we have developed a new research area called 'Phytochemical Genomics'. In this review, the development of a research platform for plant metabolomics is discussed first, to provide a better understanding of the chemical diversity of plants. Then, representative applications of metabolomics to functional genomics in a model plant, Arabidopsis thaliana, are described. The extension of integrated multi-omics analyses to non-model specialized plants, e.g., medicinal plants, is presented, including the identification of novel genes, metabolites and networks for the biosynthesis of flavonoids, alkaloids, sulfur-containing metabolites and terpenoids. Further, functional genomics studies on a variety of medicinal plants is presented. I also discuss future trends in pharmacognosy and related sciences.

Targeted LC-MS Analysis for Plant Secondary Metabolites

Recent technological developments and methodological advances of both liquid chromatography (LC) and mass spectrometry (MS) have allowed LC-MS-based plant metabolomics to become a common tool for investigating quantity, quality, and chemical diversity of plant metabolites. Targeted LC-MS metabolite analysis focuses on the detection and quantitation of the researcher's target metabolites. Whilst the word "target analysis" has been used for the analytical measurement to obtain the absolute concentrations evaluated by authentic and/or stable-isotope-labeled standards, over time the phrase came to be also used in a broad sense for the measurement of annotatable metabolites by structural information obtained from the combination of different strategies such as MS/MS analysis, reference extracts, mutant analysis and database search. Here, we describe a general protocol for targeted LC-MS metabolite profiling of plant secondary metabolites. Additionally, we introduce some examples of peak annotation using the combination approach.

Fire and Brimstone: Molecular Interactions between Sulfur and Glucosinolate Biosynthesis in Model and Crop Brassicaceae

Glucosinolates (GSLs) represent one of the most widely studied classes of plant secondary metabolite, and have a wide range of biological activities. Their unique properties also affect livestock and human health, and have been harnessed for food and other end-uses. Since GSLs are sulfur (S)-rich there are many lines of evidence suggesting that plant S status plays a key role in determining plant GSL content. However, there is still a need to establish a detailed knowledge of the distribution and remobilization of S and GSLs throughout the development of Brassica crops, and to represent this in terms of primary and secondary sources and sinks. The increased genome complexity, gene duplication and divergence within brassicas, together with their ontogenetic plasticity during crop development, appear to have a marked effect on the regulation of S and GSLs. Here, we review the current understanding of inorganic S (sulfate) assimilation into organic S forms, including GSLs and their precursors, the intracellular and inter-organ transport of inorganic and organic S forms, and the accumulation of GSLs in specific tissues. We present this in the context of overlapping sources and sinks, transport processes, signaling molecules and their associated molecular interactions. Our analysis builds on recent insights into the molecular regulation of sulfate uptake and transport by different transporters, transcription factors and miRNAs, and the role that these may play in GSL biosynthesis. We develop a provisional model describing the key processes that could be targeted in crop breeding programs focused on modifying GSL content.

Flower power and the mustard bomb: Comparative analysis of gene and genome duplications in glucosinolate biosynthetic pathway evolution in Cleomaceae and Brassicaceae

PREMISE OF THE STUDY

Glucosinolates (GS) are a class of plant secondary metabolites that provide defense against herbivores and may play an important role in pollinator attraction. Through coevolution with plant-interacting organisms, glucosinolates have diversified into a variety of chemotypes through gene sub- and neofunctionalization. Polyploidy has been of major importance in the evolutionary history of these gene families and the development of chemically separate GS types. Here we study the effects of polyploidy in Tarenaya hassleriana (Cleomaceae) on the genes underlying GS biosynthesis.

METHODS

We established putative orthologs of all gene families involved in GS biosynthesis through sequence comparison and their duplication method through calculation of synonymous substitution ratios, phylogenetic gene trees, and synteny comparison. We drew expression data from previously published work of the identified genes and compared expression in several tissues.

KEY RESULTS

We show that the majority of gene family expansion in T. hassleriana has taken place through the retention of polyploid duplicates, together with tandem and transpositional duplicates. We also show that the large majority (>75%) is actively expressed either globally or in specific tissues. We show that MAM and CYP83 gene families, which are crucial to GS diversification in Brassicaceae, are also recruited into specific tissue expression pathways in Cleomaceae.

CONCLUSIONS

Many GS genes have expanded through polyploidy, gene transposition duplication, and tandem duplication in Cleomaceae. Duplicate retention through these mechanisms is similar to A. thaliana, but based on the expression of GS genes, Cleomaceae-specific diversification of GS genes has taken place.

The Regulation of Essential Amino Acid Synthesis and Accumulation in Plants

Although amino acids are critical for all forms of life, only proteogenic amino acids that humans and animals cannot synthesize de novo and therefore must acquire in their diets are classified as essential. Nine amino acids-lysine, methionine, threonine, phenylalanine, tryptophan, valine, isoleucine, leucine, and histidine-fit this definition. Despite their nutritional importance, several of these amino acids are present in limiting quantities in many of the world's major crops. In recent years, a combination of reverse genetic and biochemical approaches has been used to define the genes encoding the enzymes responsible for synthesizing, degrading, and regulating these amino acids. In this review, we describe recent advances in our understanding of the metabolism of the essential amino acids, discuss approaches for enhancing their levels in plants, and appraise efforts toward their biofortification in crop plants.

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.

Proteasome targeting of proteins in Arabidopsis leaf mesophyll, epidermal and vascular tissues

Protein and transcript levels are partly decoupled as a function of translation efficiency and protein degradation. Selective protein degradation via the Ubiquitin-26S proteasome system (UPS) ensures protein homeostasis and facilitates adjustment of protein abundance during changing environmental conditions. Since individual leaf tissues have specialized functions, their protein composition is different and hence also protein level regulation is expected to differ. To understand UPS function in a tissue-specific context we developed a method termed Meselect to effectively and rapidly separate Arabidopsis thaliana leaf epidermal, vascular and mesophyll tissues. Epidermal and vascular tissue cells are separated mechanically, while mesophyll cells are obtained after rapid protoplasting. The high yield of proteins was sufficient for tissue-specific proteome analyses after inhibition of the proteasome with the specific inhibitor Syringolin A (SylA) and affinity enrichment of ubiquitylated proteins. SylA treatment of leaves resulted in the accumulation of 225 proteins and identification of 519 ubiquitylated proteins. Proteins that were exclusively identified in the three different tissue types are consistent with specific cellular functions. Mesophyll cell proteins were enriched for plastid membrane translocation complexes as targets of the UPS. Epidermis enzymes of the TCA cycle and cell wall biosynthesis specifically accumulated after proteasome inhibition, and in the vascular tissue several enzymes involved in glucosinolate biosynthesis were found to be ubiquitylated. Our results demonstrate that protein level changes and UPS protein targets are characteristic of the individual leaf tissues and that the proteasome is relevant for tissue-specific functions.

Associative transcriptomics study dissects the genetic architecture of seed glucosinolate content in Brassica napus

Breeding new varieties with low seed glucosinolate (GS) concentrations has long been a prime target in Brassica napus. In this study, a novel association mapping methodology termed 'associative transcriptomics' (AT) was applied to a panel of 101 B. napus lines to define genetic regions and also candidate genes controlling total seed GS contents. Over 100,000 informative single-nucleotide polymorphisms (SNPs) and gene expression markers (GEMs) were developed for AT analysis, which led to the identification of 10 SNP and 7 GEM association peaks. Within these peaks, 26 genes were inferred to be involved in GS biosynthesis. A weighted gene co-expression network analysis provided additional 40 candidate genes. The transcript abundance in leaves of two candidate genes, BnaA.GTR2a located on chromosome A2 and BnaC.HAG3b on C9, was correlated with seed GS content, explaining 18.8 and 16.8% of phenotypic variation, respectively. Resequencing of genomic regions revealed six new SNPs in BnaA.GTR2a and four insertions or deletions in BnaC.HAG3b. These deletion polymorphisms were then successfully converted into polymerase chain reaction-based diagnostic markers that can, due to high linkage disequilibrium observed in these regions of the genome, be used for marker-assisted breeding for low seed GS lines.

Whole genome and tandem duplicate retention facilitated glucosinolate pathway diversification in the mustard family

Plants share a common history of successive whole-genome duplication (WGD) events retaining genomic patterns of duplicate gene copies (ohnologs) organized in conserved syntenic blocks. Duplication was often proposed to affect the origin of novel traits during evolution. However, genetic evidence linking WGD to pathway diversification is scarce. We show that WGD and tandem duplication (TD) accelerated genetic versatility of plant secondary metabolism, exemplified with the glucosinolate (GS) pathway in the mustard family. GS biosynthesis is a well-studied trait, employing at least 52 biosynthetic and regulatory genes in the model plant Arabidopsis. In a phylogenomics approach, we identified 67 GS loci in Aethionema arabicum of the tribe Aethionemae, sister group to all mustard family members. All but one of the Arabidopsis GS gene families evolved orthologs in Aethionema and all but one of the orthologous sequence pairs exhibit synteny. The 45% fraction of duplicates among all protein-coding genes in Arabidopsis was increased to 95% and 97% for Arabidopsis and Aethionema GS pathway inventory, respectively. Compared with the 22% average for all protein-coding genes in Arabidopsis, 52% and 56% of Aethionema and Arabidopsis GS loci align to ohnolog copies dating back to the last common WGD event. Although 15% of all Arabidopsis genes are organized in tandem arrays, 45% and 48% of GS loci in Arabidopsis and Aethionema descend from TD, respectively. We describe a sequential combination of TD and WGD events driving gene family extension, thereby expanding the evolutionary playground for functional diversification and thus potential novelty and success.

How to discover a metabolic pathway? An update on gene identification in aliphatic glucosinolate biosynthesis, regulation and transport

Identification of enzymes, regulators and transporters involved in different metabolic processes is the foundation to understand how organisms function. There are, however, many difficulties in identifying candidate genes as well as in proving their in vivo roles. In this review, we describe different approaches utilized in Arabidopsis thaliana to identify gene candidates and experiments required to prove the function of a given candidate. For example, we use the production of methionine-derived aliphatic glucosinolates that represent major defence compounds in A. thaliana. Nearly all biosynthetic genes, as well as the first sets of regulators and transporters, have been identified. An array of approaches, i.e. classical mapping, quantitative trait loci (QTL) mapping, eQTL mapping, co-expression, genome wide association studies (GWAS), mutant screens and phylogenetic analyses, has been exploited to increase the number of identified genes. Here we summarize the lessons learned from the different approaches used over the years with the aim to help designing and combining new approaches in the future.

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

In Arabidopsis thaliana the evolutionary and functional relationship between Leu biosynthesis and the Met chain elongation pathway, the first part of glucosinolate formation, is well documented. Nevertheless the exact functions of some pathway components are still unclear. Isopropylmalate isomerase (IPMI), an enzyme usually involved in Leu biosynthesis, is a heterodimer consisting of a large and a small subunit. While the large protein is encoded by a single gene (ISOPROPYLMALATE ISOMERASE LARGE SUBUNIT1), three genes encode small subunits (ISOPROPYLMALATE ISOMERASE SMALL SUBUNIT1 to 3). We have now analyzed small subunit 1 (ISOPROPYLMALATE ISOMERASE SMALL SUBUNIT1) employing artificial microRNA for a targeted knockdown of the encoding gene. Strong reduction of corresponding mRNA levels to less than 5% of wild-type levels resulted in a severe phenotype with stunted growth, narrow pale leaf blades with green vasculature and abnormal adaxial-abaxial patterning as well as anomalous flower morphology. Supplementation of the knockdown plants with leucine could only partially compensate for the morphological and developmental abnormalities. Detailed metabolite profiling of the knockdown plants revealed changes in the steady state levels of isopropylmalate and glucosinolates as well as their intermediates demonstrating a function of IPMI SSU1 in both leucine biosynthesis and the first cycle of Met chain elongation. Surprisingly the levels of free leucine slightly increased suggesting an imbalanced distribution of leucine within cells and/or within plant tissues.

Genomic signature of adaptation to climate in Medicago truncatula

Local adaptation and adaptive clines are pervasive in natural plant populations, yet the effects of these types of adaptation on genomic diversity are not well understood. With a data set of 202 accessions of Medicago truncatula genotyped at almost 2 million single nucleotide polymorphisms, we used mixed linear models to identify candidate loci responsible for adaptation to three climatic gradients-annual mean temperature (AMT), precipitation in the wettest month (PWM), and isothermality (ITH)-representing the major axes of climate variation across the species' range. Loci with the strongest association to these climate gradients tagged genome regions with high sequence similarity to genes with functional roles in thermal tolerance, drought tolerance, or resistance to herbivores of pathogens. Genotypes at these candidate loci also predicted the performance of an independent sample of plant accessions grown in climate-controlled conditions. Compared to a genome-wide sample of randomly drawn reference SNPs, candidates for two climate gradients, AMT and PWM, were significantly enriched for genic regions, and genome segments flanking genic AMT and PWM candidates harbored less nucleotide diversity, elevated differentiation between haplotypes carrying alternate alleles, and an overrepresentation of the most common haplotypes. These patterns of diversity are consistent with a history of soft selective sweeps acting on loci underlying adaptation to climate, but not with a history of long-term balancing selection.

Transporters in plant sulfur metabolism

Sulfur is an essential nutrient, necessary for synthesis of many metabolites. The uptake of sulfate, primary and secondary assimilation, the biosynthesis, storage, and final utilization of sulfur (S) containing compounds requires a lot of movement between organs, cells, and organelles. Efficient transport systems of S-containing compounds across the internal barriers or the plasma membrane and organellar membranes are therefore required. Here, we review a current state of knowledge of the transport of a range of S-containing metabolites within and between the cells as well as of their long distance transport. An improved understanding of mechanisms and regulation of transport will facilitate successful engineering of the respective pathways, to improve the plant yield, biotic interaction and nutritional properties of crops.

Genome-wide analysis of branched-chain amino acid levels in Arabidopsis seeds

Branched-chain amino acids (BCAAs) are three of the nine essential amino acids in human and animal diets and are important for numerous processes in development and growth. However, seed BCAA levels in major crops are insufficient to meet dietary requirements, making genetic improvement for increased and balanced seed BCAAs an important nutritional target. Addressing this issue requires a better understanding of the genetics underlying seed BCAA content and composition. Here, a genome-wide association study and haplotype analysis for seed BCAA traits in Arabidopsis thaliana revealed a strong association with a chromosomal interval containing two branched-chain amino acid transferases, BCAT1 and BCAT2. Linkage analysis, reverse genetic approaches, and molecular complementation analysis demonstrated that allelic variation at BCAT2 is responsible for the natural variation of seed BCAAs in this interval. Complementation analysis of a bcat2 null mutant with two significantly different alleles from accessions Bayreuth-0 and Shahdara is consistent with BCAT2 contributing to natural variation in BCAA levels, glutamate recycling, and free amino acid homeostasis in seeds in an allele-dependent manner. The seed-specific phenotype of bcat2 null alleles, its strong transcription induction during late seed development, and its subcellular localization to the mitochondria are consistent with a unique, catabolic role for BCAT2 in BCAA metabolism in seeds.

Novel insights into the function of Arabidopsis R2R3-MYB transcription factors regulating aliphatic glucosinolate biosynthesis

Arabidopsis transcription factors, MYB28, MYB29 and MYB76, positively regulate aliphatic glucosinolate (AGSL) biosynthesis. Mutual transcriptional regulation among these MYB genes makes it difficult to elucidate their individual function simply by analyzing knock-out mutants or ectopically overexpressing lines of these genes. In this study, we constructed transgenic lines expressing each MYB gene driven by its own promoter in the myb28myb29 background, where the expression of the endogenous MYB28, MYB29 and MYB76 was repressed with no AGSL accumulation. In leaves, transgenic MYB28 expression activated AGSL biosynthetic genes and restored accumulation of AGSLs with short side chains. Transgenic MYB29 expression activated the same biosynthetic pathway, but induction of the genes involved in side chain elongation was weaker than that by MYB28, resulting in a weaker recovery of AGSLs. Neither MYB28 nor MYB29 recovered long-chain AGSL accumulation. MYB76 was considered to require both MYB28 and MYB29 for its normal level of expression in leaves, and could not activate AGSL biosynthesis on its own. Interestingly, the accumulation in seeds of long- and short-chain AGSLs was restored by transgenic expression of MYB28 and MYB76, respectively. A sulfur stress experiment indicated that MYB28 expression was induced by sulfur deficiency, while the expression levels of MYB29 and MYB76 were positively correlated with sulfur concentration. This study illustrated how the individual MYBs work in regulating AGSL biosynthesis when expressed alone under normal transcriptional regulation.

Integrated LC-MS/MS system for plant metabolomics

Liquid chromatography-tandem mass spectrometry (LC-MS/MS) is highly sensitive, selective, and enables extensive detection of metabolites within a sample. The result allows us to characterize comprehensive metabolite accumulation patterns without dependence on authentic standard compounds and isolation of the individual metabolites. A reference database search is essential for the structural assignment process of un-targeted MS and MS/MS data. Moreover, the characterization of unknown metabolites is challenging, since these cannot be assigned a candidate structure by using a reference database. In this case study, integrated LC-MS/MS based plant metabolomics allows us to detect several hundred metabolites in a sample; and integrated omics analyses, e.g., large-scale reverse genetics, linkage mapping, and association mapping, provides a powerful tool for candidate structure selection or rejection. We also examine emerging technology and applications for LC-MS/MS-based un-targeted plant metabolomics. These activities promote the characterization of massive extended detectable metabolites.

Integrated proteomics and metabolomics of Arabidopsis acclimation to gene-dosage dependent perturbation of isopropylmalate dehydrogenases

Maintaining metabolic homeostasis is critical for plant growth and development. Here we report proteome and metabolome changes when the metabolic homeostasis is perturbed due to gene-dosage dependent mutation of Arabidopsis isopropylmalate dehydrogenases (IPMDHs). By integrating complementary quantitative proteomics and metabolomics approaches, we discovered that gradual ablation of the oxidative decarboxylation step in leucine biosynthesis caused imbalance of amino acid homeostasis, redox changes and oxidative stress, increased protein synthesis, as well as a decline in photosynthesis, which led to rearrangement of central metabolism and growth retardation. Disruption of IPMDHs involved in aliphatic glucosinolate biosynthesis led to synchronized increase of both upstream and downstream biosynthetic enzymes, and concomitant repression of the degradation pathway, indicating metabolic regulatory mechanisms in controlling glucosinolate biosynthesis.

Current metabolomics: practical applications

The field of metabolomics continues to grow rapidly over the last decade and has been proven to be a powerful technology in predicting and explaining complex phenotypes in diverse biological systems. Metabolomics complements other omics, such as transcriptomics and proteomics and since it is a 'downstream' result of gene expression, changes in the metabolome is considered to best reflect the activities of the cell at a functional level. Thus far, metabolomics might be the sole technology capable of detecting complex, biologically essential changes. As one of the omics technology, metabolomics has exciting applications in varied fields, including medical science, synthetic biology, medicine, and predictive modeling of plant, animal and microbial systems. In addition, integrated applications with genomics, transcriptomics, and proteomics provide greater understanding of global system biology. In this review, we discuss recent applications of metabolomics in microbiology, plant, animal, food, and medical science.

QTL analysis using SNP markers developed by next-generation sequencing for identification of candidate genes controlling 4-methylthio-3-butenyl glucosinolate contents in roots of radish, Raphanus sativus L

SNP markers for QTL analysis of 4-MTB-GSL contents in radish roots were developed by determining nucleotide sequences of bulked PCR products using a next-generation sequencer. DNA fragments were amplified from two radish lines by multiplex PCR with six primer pairs, and those amplified by 2,880 primer pairs were mixed and sequenced. By assembling sequence data, 1,953 SNPs in 750 DNA fragments, 437 of which have been previously mapped in a linkage map, were identified. A linkage map of nine linkage groups was constructed with 188 markers, and five QTLs were detected in two F(2) populations, three of them accounting for more than 50% of the total phenotypic variance being repeatedly detected. In the identified QTL regions, nine SNP markers were newly produced. By synteny analysis of the QTLs regions with Arabidopsis thaliana and Brassica rapa genome sequences, three candidate genes were selected, i.e., RsMAM3 for production of aliphatic glucosinolates linked to GSL-QTL-4, RsIPMDH1 for leucine biosynthesis showing strong co-expression with glucosinolate biosynthesis genes linked to GSL-QTL-2, and RsBCAT4 for branched-chain amino acid aminotransferase linked to GSL-QTL-1. Nucleotide sequences and expression of these genes suggested their possible function in 4MTB-GSL biosynthesis in radish roots.

Partial deficiency of isoleucine impairs root development and alters transcript levels of the genes involved in branched-chain amino acid and glucosinolate metabolism in Arabidopsis

Isoleucine is one of the branched-chain amino acids (BCAAs) that are essential substrates for protein synthesis in all organisms. Although the metabolic pathway for isoleucine has been well characterized in higher plants, it is not known whether it plays a specific role in plant development. In this study, an Arabidopsis mutant, lib (low isoleucine biosynthesis), that has defects in both cell proliferation and cell expansion processes during root development, was characterized. The lib mutant carries a T-DNA insertion in the last exon of the OMR1 gene that encodes a threonine deaminase/dehydratase (TD). TD catalyses the deamination and dehydration of threonine, which is the first and also the committed step in the biosynthesis of isoleucine. This T-DNA insertion results in a partial deficiency of isoleucine in lib root tissues but it does not affect its total protein content. Application of exogenous isoleucine or introduction of a wild-type OMR1 gene into the lib mutant can completely rescue the mutant phenotypes. These results reveal an important role for isoleucine in plant development. In addition, microarray analysis indicated that the partial deficiency of isoleucine in the lib mutant triggers a decrease in transcript levels of the genes encoding the major enzymes involved in the BCAA degradation pathway; the analysis also indicated that many genes involved in the biosynthesis of methionine-derived glucosinolates are up-regulated.

RIKEN tandem mass spectral database (ReSpect) for phytochemicals: a plant-specific MS/MS-based data resource and database

The fragment pattern analysis of tandem mass spectrometry (MS/MS) has long been used for the structural characterization of metabolites. The construction of a plant-specific MS/MS data resource and database will enable complex phytochemical structures to be narrowed down to candidate structures. Therefore, a web-based database of MS/MS data pertaining to phytochemicals was developed and named ReSpect (RIKEN tandem mass spectral database). Of the 3595 metabolites in ReSpect, 76% were derived from 163 literature reports, whereas the rest was obtained from authentic standards. As a main web application of ReSpect, a fragment search was established based on only the m/z values of query data and records. The confidence levels of the annotations were managed using the MS/MS fragmentation association rule, which is an algorithm for discovering common fragmentations in MS/MS data. Using this data resource and database, a case study was conducted for the annotation of untargeted MS/MS data that were selected after quantitative trait locus analysis of the accessions (Gifu and Miyakojima) of a model legume Lotus japonicus. In the case study, unknown metabolites were successfully narrowed down to putative structures in the website.

Rapid profiling of intact glucosinolates in Arabidopsis leaves by UHPLC-QTOFMS using a charged surface hybrid column

INTRODUCTION

The analysis of glucosinolates (GS) is traditionally performed by reverse-phase liquid chromatography coupled to ultraviolet detection after a time-consuming desulphation step, which is required for increased retention. Simpler and more efficient alternative methods that can shorten both sample preparation and analysis are much needed.

OBJECTIVE

To evaluate the feasibility of using ultrahigh-pressure liquid chromatography coupled to quadrupole time-of-flight mass spectrometry (UHPLC-QTOFMS) for the rapid profiling of intact GS.

METHODOLOGY

A simple and short extraction of GS from Arabidopsis thaliana leaves was developed. Four sub-2 µm reverse-phase columns were tested for the rapid separation of these polar compounds using formic acid as the chromatographic additive. High-resolution QTOFMS was used to detect and identify GS.

RESULTS

A novel charged surface hybrid (CSH) column was found to provide excellent retention and separation of GS within a total running time of 11 min. Twenty-one GS could be identified based on their accurate mass as well as isotopic and fragmentation patterns. The method was applied to determine the changes in GS content that occur after herbivory in Arabidopsis. In addition, we evaluated its applicability to the profiling of other Brassicaceae species.

CONCLUSION

The method developed can profile the full range of GS, including the most polar ones, in a shorter time than previous methods, and is highly compatible with mass spectrometric detection.

Glucosinolate structures in evolution

By 2000, around 106 natural glucosinolates (GSLs) were probably documented. In the past decade, 26 additional natural GSL structures have been elucidated and documented. Hence, the total number of documented GSLs from nature by 2011 can be estimated to around 132. A considerable number of additional suggested structures are concluded not to be sufficiently documented. In many cases, NMR spectroscopy would have provided the missing structural information. Of the GSLs documented in the past decade, several are of previously unexpected structures and occur at considerable levels. Most originate from just four species: Barbarea vulgaris, Arabidopsis thaliana, Eruca sativa and Isatis tinctoria. Acyl derivatives of known GSLs comprised 15 of the 26 newly documented structures, while the remaining exhibited new substitution patterns or chain length, or contained a mercapto group or related thio-functionality. GSL identification methods are reviewed, and the importance of using authentic references and structure-sensitive detection methods such as MS and NMR is stressed, especially when species with relatively unknown chemistry are analyzed. An example of qualitative GSL analysis is presented with experimental details (group separation and HPLC of both intact and desulfated GSLs, detection and structure determination by UV, MS, NMR and susceptibility to myrosinase) with emphasis on the use of NMR for structure elucidation of even minor GSLs and GSL hydrolysis products. The example includes identification of a novel GSL, (R)-2-hydroxy-2-(3-hydroxyphenyl)ethylglucosinolate. Recent investigations of GSL evolution, based on investigations of species with well established phylogeny, are reviewed. From the relatively few such investigations, it is already clear that GSL profiles are regularly subject to evolution. This result is compatible with natural selection for specific GSL side chains. The probable existence of structure-specific GSL catabolism in intact plants suggests that biochemical evolution of GSLs has more complex implications than the mere liberation of a different hydrolysis product upon tissue disruption.

Integrative systems biology: an attempt to describe a simple weed

Genome-scale studies hold great promise for revealing novel plant biology. Because of the complexity of these techniques, numerous considerations need to be made before embarking on a study. Here we focus on the Arabidopsis model system because of the wealth of available genome-scale data. Many approaches are available that provide genome-scale information regarding the state of a given organism (e.g. genomics, epigenomics, transcriptomics, proteomics, metabolomics interactomics, ionomics, phenomics, etc.). Integration of all of these types of data will be necessary for a comprehensive description of Arabidopsis. In this review we propose that 'triangulation' among transcriptomics, proteomics and metabolomics is a meaningful approach for beginning this integrative analysis and uncovering a systems level perspective of Arabidopsis biology.

Characterization of metabolite quantitative trait loci and metabolic networks that control glucosinolate concentration in the seeds and leaves of Brassica napus

• Glucosinolates are a major class of secondary metabolites found in the Brassicaceae, whose degradation products are proving to be increasingly important for human health and in crop protection. • The genetic and metabolic basis of glucosinolate accumulation was dissected through analysis of total glucosinolate concentration and its individual components in both leaves and seeds of a doubled-haploid (DH) mapping population of oilseed rape/canola (Brassica napus). • The quantitative trait loci (QTL) that had an effect on glucosinolate concentration in either or both of the organs were integrated, resulting in 105 metabolite QTL (mQTL). Pairwise correlations between individual glucosinolates and prior knowledge of the metabolic pathways involved in the biosynthesis of different glucosinolates allowed us to predict the function of genes underlying the mQTL. Moreover, this information allowed us to construct an advanced metabolic network and associated epistatic interactions responsible for the glucosinolate composition in both leaves and seeds of B. napus. • A number of previously unknown potential regulatory relationships involved in glucosinolate synthesis were identified and this study illustrates how genetic variation can affect a biochemical pathway.

Advances in omics and bioinformatics tools for systems analyses of plant functions

Omics and bioinformatics are essential to understanding the molecular systems that underlie various plant functions. Recent game-changing sequencing technologies have revitalized sequencing approaches in genomics and have produced opportunities for various emerging analytical applications. Driven by technological advances, several new omics layers such as the interactome, epigenome and hormonome have emerged. Furthermore, in several plant species, the development of omics resources has progressed to address particular biological properties of individual species. Integration of knowledge from omics-based research is an emerging issue as researchers seek to identify significance, gain biological insights and promote translational research. From these perspectives, we provide this review of the emerging aspects of plant systems research based on omics and bioinformatics analyses together with their associated resources and technological advances.

Advances in omics and bioinformatics tools for systems analyses of plant functions

Omics and bioinformatics are essential to understanding the molecular systems that underlie various plant functions. Recent game-changing sequencing technologies have revitalized sequencing approaches in genomics and have produced opportunities for various emerging analytical applications. Driven by technological advances, several new omics layers such as the interactome, epigenome and hormonome have emerged. Furthermore, in several plant species, the development of omics resources has progressed to address particular biological properties of individual species. Integration of knowledge from omics-based research is an emerging issue as researchers seek to identify significance, gain biological insights and promote translational research. From these perspectives, we provide this review of the emerging aspects of plant systems research based on omics and bioinformatics analyses together with their associated resources and technological advances.

Glucosinolate biosynthetic genes in Brassica rapa

Glucosinolates (GS) are a group of amino acid-derived secondary metabolites found throughout the Cruciferae family. Glucosinolates and their degradation products play important roles in pathogen and insect interactions, as well as in human health. In order to elucidate the glucosinolate biosynthetic pathway in Brassica rapa, we conducted comparative genomic analyses of Arabidopsis thaliana and B. rapa on a genome-wide level. We identified 102 putative genes in B. rapa as the orthologs of 52 GS genes in A. thaliana. All but one gene was successfully mapped on 10 chromosomes. Most GS genes exist in more than one copy in B. rapa. A high co-linearity in the glucosinolate biosynthetic pathway between A. thaliana and B. rapa was also established. The homologous GS genes in B. rapa and A. thaliana share 59-91% nucleotide sequence identity and 93% of the GS genes exhibit synteny between B. rapa and A. thaliana. Moreover, the structure and arrangement of the B. rapa GS (BrGS) genes correspond with the known evolutionary divergence of B. rapa, and may help explain the profiles and accumulation of GS in B. rapa.

Structural and functional evolution of isopropylmalate dehydrogenases in the leucine and glucosinolate pathways of Arabidopsis thaliana

The methionine chain-elongation pathway is required for aliphatic glucosinolate biosynthesis in plants and evolved from leucine biosynthesis. In Arabidopsis thaliana, three 3-isopropylmalate dehydrogenases (AtIPMDHs) play key roles in methionine chain-elongation for the synthesis of aliphatic glucosinolates (e.g. AtIPMDH1) and leucine (e.g. AtIPMDH2 and AtIPMDH3). Here we elucidate the molecular basis underlying the metabolic specialization of these enzymes. The 2.25 Å resolution crystal structure of AtIPMDH2 was solved to provide the first detailed molecular architecture of a plant IPMDH. Modeling of 3-isopropylmalate binding in the AtIPMDH2 active site and sequence comparisons of prokaryotic and eukaryotic IPMDH suggest that substitution of one active site residue may lead to altered substrate specificity and metabolic function. Site-directed mutagenesis of Phe-137 to a leucine in AtIPMDH1 (AtIPMDH1-F137L) reduced activity toward 3-(2'-methylthio)ethylmalate by 200-fold, but enhanced catalytic efficiency with 3-isopropylmalate to levels observed with AtIPMDH2 and AtIPMDH3. Conversely, the AtIPMDH2-L134F and AtIPMDH3-L133F mutants enhanced catalytic efficiency with 3-(2'-methylthio)ethylmalate ∼100-fold and reduced activity for 3-isopropylmalate. Furthermore, the altered in vivo glucosinolate profile of an Arabidopsis ipmdh1 T-DNA knock-out mutant could be restored to wild-type levels by constructs expressing AtIPMDH1, AtIPMDH2-L134F, or AtIPMDH3-L133F, but not by AtIPMDH1-F137L. These results indicate that a single amino acid substitution results in functional divergence of IPMDH in planta to affect substrate specificity and contributes to the evolution of specialized glucosinolate biosynthesis from the ancestral leucine pathway.

Natural diversity and adaptation in plant secondary metabolism

Technological advances in metabolomics, transcriptomics and genomics have facilitated the detection of genes that contribute to diversification in plant secondary metabolism. Statistical tools from molecular population genetics may help in evaluating whether the corresponding genes or genomic regions carry a signature of selection and answering the question of whether novel compounds are 'adaptive'. Gene duplication fuels diversification in plant secondary metabolism and the evolutionary mechanism for adaptation may follow a path of neofunctionalization subsequent to gene duplication, or gene duplication may occur subsequent to--and resolve--an adaptive conflict present in a single ancestral gene sequence.

Metabolic engineering in Nicotiana benthamiana reveals key enzyme functions in Arabidopsis indole glucosinolate modification

Indole glucosinolates, derived from the amino acid Trp, are plant secondary metabolites that mediate numerous biological interactions between cruciferous plants and their natural enemies, such as herbivorous insects, pathogens, and other pests. While the genes and enzymes involved in the Arabidopsis thaliana core biosynthetic pathway, leading to indol-3-yl-methyl glucosinolate (I3M), have been identified and characterized, the genes and gene products responsible for modification reactions of the indole ring are largely unknown. Here, we combine the analysis of Arabidopsis mutant lines with a bioengineering approach to clarify which genes are involved in the remaining biosynthetic steps in indole glucosinolate modification. We engineered the indole glucosinolate biosynthesis pathway into Nicotiana benthamiana, showing that it is possible to produce indole glucosinolates in a noncruciferous plant. Building upon this setup, we demonstrate that all members of a small gene subfamily of cytochrome P450 monooxygenases, CYP81Fs, are capable of carrying out hydroxylation reactions of the glucosinolate indole ring, leading from I3M to 4-hydroxy-indol-3-yl-methyl and/or 1-hydroxy-indol-3-yl-methyl glucosinolate intermediates, and that these hydroxy intermediates are converted to 4-methoxy-indol-3-yl-methyl and 1-methoxy-indol-3-yl-methyl glucosinolates by either of two family 2 O-methyltransferases, termed indole glucosinolate methyltransferase 1 (IGMT1) and IGMT2.