Glucosinolate biosynthetic genes in Brassica rapa

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

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

Simulated galactic cosmic ray exposure activates dose-dependent DNA repair response and down regulates glucosinolate pathways in arabidopsis seedlings

Outside the protection of Earth's magnetic field, organisms are constantly exposed to space radiation consisting of energetic protons and other heavier charged particles. With the goal of crewed Mars exploration, the production of fresh food during long duration space missions is critical for meeting astronauts' nutritional and psychological needs. However, the biological effects of space radiation on plants have not been sufficiently investigated and characterized. To that end, 10-day-old Arabidopsis seedlings were exposed to simulated Galactic Cosmic Rays (GCR) and assessed for transcriptomic changes. The simulated GCR irradiation was carried out in the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Lab (BNL). The exposures were conducted acutely for two dose points at 40 cGy or 80 cGy, with sequential delivery of proton, helium, oxygen, silicon, and iron ions. Control and irradiated seedlings were then harvested and preserved in RNAlater at 3 hrs post irradiation. Total RNA was isolated for transcriptomic analyses using RNAseq. The data revealed that the transcriptomic responses were dose-dependent, with significant upregulation of DNA repair pathways and downregulation of glucosinolate biosynthetic pathways. Glucosinolates are important for plant pathogen defense and for the taste of a plant, which are both relevant to growing plants for spaceflight. These findings fill in knowledge gaps of how plants respond to radiation in beyond-Earth environments.

CRISPR/Cas9-mediated BoaAOP2s editing alters aliphatic glucosinolate side-chain metabolic flux and increases the glucoraphanin content in Chinese kale

Glucoraphanin (GRA) is an aliphatic glucosinolate (GSL), and its hydrolysis product has powerful anticancer activity. ALKENYL HYDROXALKYL PRODUCING 2 (AOP2) gene, encodes a 2-oxoglutarate-dependent dioxygenase, which can catalyze GRA to form gluconapin (GNA). However, GRA only present in trace amounts in Chinese kale. To increase the content of GRA in Chinese kale, three copies of BoaAOP2 were isolated and edited using CRISPR/Cas9 system. The content of GRA was 11.71- to 41.29-fold (0.082-0.289 μmol g^-1 FW) higher in T1 generation of boaaop2 mutants than in wild-type plants, and this was accompanied by an increase in the GRA/GNA ratio and reductions in the content of GNA and total aliphatic GSLs. BoaAOP2.1 is an effective gene for the alkenylation of aliphatic GSLs in Chinese kale. Overall, targeted editing of CRISPR/Cas9-mediated BoaAOP2s altered aliphatic GSL side-chain metabolic flux and enhanced the GRA content in Chinese kale, suggesting that metabolic engineering of BoaAOP2s has huge potential in improving nutritional quality of Chinese kale.

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

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

Differences in pseudogene evolution contributed to the contrasting flavors of turnip and Chiifu, two Brassica rapa subspecies

Pseudogenes are important resources for investigation of genome evolution and genomic diversity because they are nonfunctional but have regulatory effects that influence plant adaptation and diversification. However, few systematic comparative analyses of pseudogenes in closely related species have been conducted. Here, we present a turnip (Brassica rapa ssp. rapa) genome sequence and characterize pseudogenes among diploid Brassica species/subspecies. The results revealed that the number of pseudogenes was greatest in Brassica oleracea (CC genome), followed by B. rapa (AA genome) and then Brassica nigra (BB genome), implying that pseudogene differences emerged after species differentiation. In Brassica AA genomes, pseudogenes were distributed asymmetrically on chromosomes because of numerous chromosomal insertions/rearrangements, which contributed to the diversity among subspecies. Pseudogene differences among subspecies were reflected in the flavor-related glucosinolate (GSL) pathway. Specifically, turnip had the highest content of pungent substances, probably because of expansion of the methylthioalkylmalate synthase-encoding gene family in turnips; these genes were converted into pseudogenes in B. rapa ssp. pekinensis (Chiifu). RNA interference-based silencing of the gene encoding 2-oxoglutarate-dependent dioxygenase 2, which is also associated with flavor and anticancer substances in the GSL pathway, resulted in increased abundance of anticancer compounds and decreased pungency of turnip and Chiifu. These findings revealed that pseudogene differences between turnip and Chiifu influenced the evolution of flavor-associated GSL metabolism-related genes, ultimately resulting in the different flavors of turnip and Chiifu.

A high-efficiency PEG-Ca2+-mediated transient transformation system for broccoli protoplasts

Transient transformation of plant protoplasts is an important method for studying gene function, subcellular localization and plant morphological development. In this study, an efficient transient transformation system was established by optimizing the plasmid concentration, PEG4000 mass concentration and genotype selection, key factors that affect transformation efficiency. Meanwhile, an efficient and universal broccoli protoplast isolation system was established. Using 0.5% (w/v) cellulase R-10 and 0.1% (w/v) pectolyase Y-23 to hydrolyze broccoli cotyledons of three different genotypes for 3 h, the yield was more than 5×10⁶/mL/g, and the viability was more than 95%, sufficient to meet the high standards for protoplasts to be used in various experiments. The average transformation efficiency of the two plasmid vectors PHG-eGFP and CP507-YFP in broccoli B1 protoplasts were 61.4% and 41.7%, respectively. Using this system, we successfully performed subcellular localization of the products of three target genes (the clubroot resistance gene CRa and two key genes regulated by glucosinolates, Bol029100 and Bol031350).The results showed that the products of all three genes were localized in the nucleus. The high-efficiency transient transformation system for broccoli protoplasts constructed in this study makes it possible to reliably acquire high-viability protoplasts in high yield. This research provides important technical support for international frontier research fields such as single-cell sequencing, spatial transcriptomics, plant somatic hybridization, gene function analysis and subcellular localization.

Glucosinolates of Sisymbrium officinale and S. orientale

Glucosinolates (GSLs) from Sysimbrium officinale and S. orientale were analyzed qualitatively and quantitatively by their desulfo-counterparts using UHPLC-DAD-MS/MS. Eight GSLs were identified in S. officinale, including Val-derived (glucoputranjivin) and Trp-derived (4-hydroxyglucobrassicin, glucobrassicin, 4-methoxyglucobrassicin, and neoglucobrassicin) as the major ones followed by Leu-derived (Isobutyl GSL), Ile-derived (glucocochlearin) and Phe/Tyr-derived (glucosinalbin). Different S. orientale plant parts contained six GSLs, with Met-derived (progoitrin, epiprogoitrin, and gluconapin) and homoPhe-derived (gluconasturtiin) as the major ones, followed by glucosinalbin and neoglucobrassicin. GSL breakdown products obtained by hydrodistillation (HD) and microwave-assisted distillation from S. officinale, as well as isopropyl isothiocyanate, as the major volatile in both isolates, were tested for their cytotoxic activity using a 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. Generally, all volatile isolates showed similar activity toward the three cancer cell lines. The best activity was shown by isopropyl isothiocyanate at a concentration of 100 µg/mL after 72 h of incubation, with 53.18% for MDA-MB-231, 56.61% for A549, and 60.02% for the T24 cell line.

Identification of key genes controlling soluble sugar and glucosinolate biosynthesis in Chinese cabbage by integrating metabolome and genome-wide transcriptome analysis

INTRODUCTION

Soluble sugar and glucosinolate are essential components that determine the flavor of Chinese cabbage and consumer preferences. However, the underlying regulatory networks that modulate the biosynthesis of soluble sugar and glucosinolate in Chinese cabbage remain largely unknown.

METHODS

The glucosinolate and carotene content in yellow inner-leaf Chinese cabbage were observed, followed by the combination of metabolome and transcriptome analysis to explore the metabolic basis of glucosinolate and soluble sugar.

RESULTS

This study observed high glucosinolate and carotene content in yellow inner-leaf Chinese cabbage, which showed a lower soluble sugar content. The differences between the yellow and the white inner-leaf Chinese cabbage were compared using the untargeted metabonomic and transcriptomic analyses in six cultivars of Chinese cabbage to explore the metabolic basis of glucosinolate and soluble sugar. Aliphatic glucosinolate and two soluble sugars (fructose and glucose) were the key metabolites that caused the difference in Chinese cabbage's glucosinolate and soluble sugar. By integrating soluble sugar and glucosinolate-associated metabolism and transcriptome data, we indicated BraA05gAOP1 and BraA04gAOP4, BraA03gHT7 and BraA01gHT4 were the glucosinolates and soluble sugar biosynthesis structural genes. Moreover, BraA01gCHR11 and BraA07gSCL1 were two vital transcription factors that regulate soluble sugar and glucosinolate biosynthesis.

DISCUSSION

These findings provide novel insights into glucosinolate and soluble sugar biosynthesis and a possible explanation for the significant difference in nutrients between yellow and white inner-leaf Chinese cabbage. Moreover, it will facilitate genetic modification to improve the Chinese cabbage's nutritional and health values.

Effects of nanocarbon solution treatment on the nutrients and glucosinolate metabolism in broccoli

•

Nanocarbon application could enhance the total protein in broccoli, directly.

•

18.75 L·ha⁻¹ nanocarbon solution greatly increase 22.9 % of glucoraphanin.

•

Nanocarbon solution obviously reduces 4 indolic glucosinolate productions.

•

Nanocarbon has great great impact on glucosinolate biosynthesis and pathway.

The effects of a nanocarbon solution on the nutrients, glucosinolate metabolism and glucoraphanin pathway in broccoli were investigated. Significant positive linear relationships were observed between the nanocarbon solution and total protein yield, although effects on the soluble sugars, vitamin C and dry matter production were not observed. All nanocarbon solutions significantly increased the glucoraphanin content (p < 0.05), and the 18.75 L·ha⁻¹ nanocarbon solution maximally increased the glucoraphanin content by 22.9 %. However, these treatments also significantly reduced the contents of glucobrassicin, 4-methoxyglucobrassicin, 4-hydroxyglucobrassicin and neoglucobrassicin. Further research demonstrated that the 18.75 L·ha⁻¹ nanocarbon solution significantly upregulated the MAM1, IPMI2, CYP79F1, FMOgs-ox2, AOP2, and TGG1 expression levels, which directly resulted in the accumulation of glucoraphanin and glucoerucin. This study provides insights into the prospective nanotechnological approaches for developing efficient and environmentally friendly nanocarbon solution for use on crops.

TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies

Vegetables are an indispensable part of the daily diet of humans. Therefore, it is vital to systematically study the genomic data of vegetables and build a platform for data sharing and analysis. In this study, a comprehensive platform for vegetables with a user-friendly Web interface-The Vegetable Information Resource (TVIR, http://tvir.bio2db.com)-was built based on the genomes of 59 vegetables. TVIR database contains numerous important functional genes, including 5215 auxin genes, 2437 anthocyanin genes, 15 002 flowering genes, 79 830 resistance genes, and 2639 glucosinolate genes of 59 vegetables. In addition, 2597 N6-methyladenosine (m6A) genes were identified, including 513 writers, 1058 erasers, and 1026 readers. A total of 2 101 501 specific clustered regularly interspaced short palindromic repeat (CRISPR) guide sequences and 17 377 miRNAs were detected and deposited in TVIR database. Information on gene synteny, duplication, and orthologs is also provided for 59 vegetable species. TVIR database contains 2 346 850 gene annotations by the Swiss-Prot, TrEMBL, Gene Ontology (GO), Pfam, and Non-redundant (Nr) databases. Synteny, Primer Design, Blast, and JBrowse tools are provided to facilitate users in conducting comparative genomic analyses. This is the first large-scale collection of vegetable genomic data and bioinformatic analysis. All genome and gene sequences, annotations, and bioinformatic results can be easily downloaded from TVIR. Furthermore, transcriptome data of 98 vegetables have been collected and collated, and can be searched by species, tissues, or different growth stages. TVIR is expected to become a key hub for vegetable research globally. The database will be updated with newly assembled vegetable genomes and comparative genomic studies in the future.

The Brassicaceae genome resource (TBGR): A comprehensive genome platform for Brassicaceae plants

The Brassicaceae is an important plant family. We built a user-friendly, web-based, comparative, and functional genomic database, The Brassicaceae Genome Resource (TBGR, http://www.tbgr.org.cn), based on 82 released genomes from 27 Brassicaceae species. The TBGR database contains a large number of important functional genes, including 4,096 glucosinolate genes, 6,625 auxin genes, 13,805 flowering genes, 36,632 resistance genes, 1,939 anthocyanin genes, and 1,231 m6A genes. A total of 1,174,049 specific guide sequences for clustered regularly interspaced short palindromic repeats and 5,856,479 transposable elements were detected in Brassicaceae. TBGR also provides information on synteny, duplication, and orthologs for 27 Brassicaceae species. The TBGR database contains 1,183,851 gene annotations obtained using the TrEMBL, Swiss-Prot, Nr, GO, and Pfam databases. The BLAST, Synteny, Primer Design, Seq_fetch, and JBrowse tools are provided to help users perform comparative genomic analyses. All the genome assemblies, gene models, annotations, and bioinformatics results can be easily downloaded from the TBGR database. We plan to improve and continuously update the database with newly assembled genomes and comparative genomic studies. We expect the TBGR database to become a key resource for the study of the Brassicaceae.

Insights into glucosinolate accumulation and metabolic pathways in Isatis indigotica Fort

BACKGROUND

Glucosinolates (GSLs) play important roles in defending against exogenous damage and regulating physiological activities in plants. However, GSL accumulation patterns and molecular regulation mechanisms are largely unknown in Isatis indigotica Fort.

RESULTS

Ten GSLs were identified in I. indigotica, and the dominant GSLs were epiprogoitrin (EPI) and indole-3-methyl GSL (I3M), followed by progoitrin (PRO) and gluconapin (GNA). The total GSL content was highest (over 20 μmol/g) in reproductive organs, lowest (less than 1.0 μmol/g) in mature organs, and medium in fresh leaves (2.6 μmol/g) and stems (1.5 μmol/g). In the seed germination process, the total GSL content decreased from 27.2 μmol/g (of seeds) to 2.7 μmol/g (on the 120th day) and then increased to 4.0 μmol/g (180th day). However, the content of indole GSL increased rapidly in the first week after germination and fluctuated between 1.13 μmol/g (28th day) and 2.82 μmol/g (150th day). Under the different elicitor treatments, the total GSL content increased significantly, ranging from 2.9-fold (mechanical damage, 3 h) to 10.7-fold (MeJA, 6 h). Moreover, 132 genes were involved in GSL metabolic pathways. Among them, no homologs of AtCYP79F2 and AtMAM3 were identified, leading to a distinctive GSL profile in I. indigotica. Furthermore, most genes involved in the GSL metabolic pathway were derived from tandem duplication, followed by dispersed duplication and segmental duplication. Purifying selection was observed, although some genes underwent relaxed selection. In addition, three tandem-arrayed GSL-OH genes showed different expression patterns, suggesting possible subfunctionalization during evolution.

CONCLUSIONS

Ten different GSLs with their accumulation patterns and 132 genes involved in the GSL metabolic pathway were explored, which laid a foundation for the study of GSL metabolism and regulatory mechanisms in I. indigotica.

Metabolite Characteristics Analysis of Siliques and Effects of Lights on the Accumulation of Glucosinolates in Siliques of Rapeseed

Glucosinolates (GSLs) are naturally occurring secondary metabolites found in the Brassicaceae family, which mainly synthesize in the siliques with a wide range of functions. In this study, we investigated the effects of lights on metabolites in siliques of rapeseed through ultra high-performance liquid chromatography (UPLC)-heated electrospray ionization (HESI)-tandem mass spectrometry (MS/MS). A total of 249 metabolites, including 29 phenolic acids, 38 flavonoids, 22 GSLs, 93 uncalculated and 67 unknown compounds, were identified in siliques of rapeseed. Meanwhile, 62 metabolites showed significant differences after shading treatment, which were mainly GSLs and unknown compounds. Interestingly, the amounts of 10 GSLs had high accumulation levels in siliques, while the expression levels of their corresponding biosynthetic genes (AOP, GSL-OH, IGMT, and ST5a) were obviously reduced after shading treatment. Further evidence showed that the amounts of GSLs were significantly reduced in seeds, in accordance with the expression profiles of transporter genes (BnaGTRs). Our findings indicated that lights could affect the accumulation and transportation of GSLs from siliques to seeds in rapeseed. Therefore, this study facilitates a better understanding of metabolic characteristics of siliques and provides insight into the importance of light for GSLs accumulation and transportation in siliques and seeds of rapeseed.

Jasmonate regulates the FAMA/mediator complex subunit 8-THIOGLUCOSIDE GLUCOHYDROLASE 1 cascade and myrosinase activity

Myrosinases are β-thioglucoside glucosidases that are unique to the Brassicales order. These enzymes hydrolyze glucosinolates to produce compounds that have direct antibiotic effects or that function as signaling molecules in the plant immune system, protecting plants from pathogens and insect pests. However, the effects of jasmonic acid (JA), a plant hormone that is crucial for plant disease resistance, on myrosinase activity remain unclear. Here, we systematically studied the effects of JA on myrosinase activity and explored the associated internal transcriptional regulation mechanisms. Exogenous application of JA significantly increased myrosinase activity, while the inhibition of endogenous JA biosynthesis and signaling reduced myrosinase activity. In addition, some myrosinase genes in Arabidopsis (Arabidopsis thaliana) were upregulated by JA. Further genetic and biochemical evidence showed that transcription factor FAMA interacted with a series of JASMONATE ZIM-DOMAIN proteins and affected JA-mediated myrosinase activity. However, among the JA-upregulated myrosinase genes, only THIOGLUCOSIDE GLUCOHYDROLASE 1 (TGG1) was positively regulated by FAMA. Further biochemical analysis showed that FAMA bound to the TGG1 promoter to directly mediate TGG1 expression in conjunction with Mediator complex subunit 8 (MED8). Together, our results provide evidence that JA acts as an important signal upstream of the FAMA/MED8-TGG1 pathway to positively regulate myrosinase activity in Arabidopsis.

Influence of Genotype on High Glucosinolate Synthesis Lines of Brassica rapa

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

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

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

The Molecular Basis of Kale Domestication: Transcriptional Profiling of Developing Leaves Provides New Insights Into the Evolution of a Brassica oleracea Vegetative Morphotype

Morphotypes of Brassica oleracea are the result of a dynamic interaction between genes that regulate the transition between vegetative and reproductive stages and those that regulate leaf morphology and plant architecture. In kales, ornate leaves, extended vegetative phase, and nutritional quality are some of the characters potentially selected by humans during domestication. We used a combination of developmental studies and transcriptomics to understand the vegetative domestication syndrome of kale. To identify candidate genes that are responsible for the evolution of domestic kale, we searched for transcriptome-wide differences among three vegetative B. oleracea morphotypes. RNA-seq experiments were used to understand the global pattern of expressed genes during a mixture of stages at one time in kale, cabbage, and the rapid cycling kale line TO1000. We identified gene expression patterns that differ among morphotypes and estimate the contribution of morphotype-specific gene expression that sets kale apart (3958 differentially expressed genes). Differentially expressed genes that regulate the vegetative to reproductive transition were abundant in all morphotypes. Genes involved in leaf morphology, plant architecture, defense, and nutrition were differentially expressed in kale. This allowed us to identify a set of candidate genes we suggest may be important in the kale domestication syndrome. Understanding candidate genes responsible for kale domestication is of importance to ultimately improve Cole crop production.

Regulation of glucosinolate biosynthesis

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

Perception of unrelated microbe-associated molecular patterns triggers conserved yet variable physiological and transcriptional changes in Brassica rapa ssp. pekinensis

Pattern-triggered immunity (PTI) includes the different transcriptional and physiological responses that enable plants to ward off microbial invasion. Surface-localized pattern-recognition receptors (PRRs) recognize conserved microbe-associated molecular patterns (MAMPs) and initiate a branched signaling cascade that culminate in an effective restriction of pathogen growth. In the model species Arabidopsis thaliana, early PTI events triggered by different PRRs are broadly conserved although their nature or intensity is dependent on the origin and features of the detected MAMP. In order to provide a functional basis for disease resistance in leafy vegetable crops, we surveyed the conservation of PTI events in Brassica rapa ssp. pekinensis. We identified the PRR homologs present in B. rapa genome and found that only one of the two copies of the bacterial Elongation factor-Tu receptor (EFR) might function. We also characterized the extent and unexpected specificity of the transcriptional changes occurring when B. rapa seedlings are treated with two unrelated MAMPs, the bacterial flagellin flg22 peptide and the fungal cell wall component chitin. Finally, using a MAMP-induced protection assay, we could show that bacterial and fungal MAMPs elicit a robust immunity in B. rapa, despite significant differences in the kinetic and amplitude of the early signaling events. Our data support the relevance of PTI for crop protection and highlight specific functional target for disease resistance breeding in Brassica crops.

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

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

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

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

Independent evolution of ancestral and novel defenses in a genus of toxic plants (Erysimum, Brassicaceae)

Phytochemical diversity is thought to result from coevolutionary cycles as specialization in herbivores imposes diversifying selection on plant chemical defenses. Plants in the speciose genus Erysimum (Brassicaceae) produce both ancestral glucosinolates and evolutionarily novel cardenolides as defenses. Here we test macroevolutionary hypotheses on co-expression, co-regulation, and diversification of these potentially redundant defenses across this genus. We sequenced and assembled the genome of E. cheiranthoides and foliar transcriptomes of 47 additional Erysimum species to construct a phylogeny from 9868 orthologous genes, revealing several geographic clades but also high levels of gene discordance. Concentrations, inducibility, and diversity of the two defenses varied independently among species, with no evidence for trade-offs. Closely related, geographically co-occurring species shared similar cardenolide traits, but not glucosinolate traits, likely as a result of specific selective pressures acting on each defense. Ancestral and novel chemical defenses in Erysimum thus appear to provide complementary rather than redundant functions.

Plants are often attacked by insects and other herbivores. As a result, they have evolved to defend themselves by producing many different chemicals that are toxic to these pests. As producing each chemical costs energy, individual plants often only produce one type of chemical that is targeted towards their main herbivore. Related species of plants often use the same type of chemical defense so, if a particular herbivore gains the ability to cope with this chemical, it may rapidly become an important pest for the whole plant family.

To escape this threat, some plants have gained the ability to produce more than one type of chemical defense. Wallflowers, for example, are a group of plants in the mustard family that produce two types of toxic chemicals: mustard oils, which are common in most plants in this family; and cardenolides, which are an innovation of the wallflowers, and which are otherwise found only in distantly related plants such as foxglove and milkweed. The combination of these two chemical defenses within the same plant may have allowed the wallflowers to escape attacks from their main herbivores and may explain why the number of wallflower species rapidly increased within the last two million years.

Züst et al. have now studied the diversity of mustard oils and cardenolides present in many different species of wallflower. This analysis revealed that almost all of the tested wallflower species produced high amounts of both chemical defenses, while only one species lacked the ability to produce cardenolides. The levels of mustard oils had no relation to the levels of cardenolides in the tested species, which suggests that the regulation of these two defenses is not linked. Furthermore, Züst et al. found that closely related wallflower species produced more similar cardenolides, but less similar mustard oils, to each other. This suggests that mustard oils and cardenolides have evolved independently in wallflowers and have distinct roles in the defense against different herbivores.

The evolution of insect resistance to pesticides and other toxins is an important concern for agriculture. Applying multiple toxins to crops at the same time is an important strategy to slow the evolution of resistance in the pests. The findings of Züst et al. describe a system in which plants have naturally evolved an equivalent strategy to escape their main herbivores. Understanding how plants produce multiple chemical defenses, and the costs involved, may help efforts to breed crop species that are more resistant to herbivores and require fewer applications of pesticides.

Expression profiles of glucosinolate biosynthetic genes in turnip (Brassica rapa var. rapa) at different developmental stages and effect of transformed flavin-containing monooxygenase genes on hairy root glucosinolate content

BACKGROUND

Glucosinolates (GSLs) are secondary metabolites, mainly existing in Brassica vegetables. Their breakdown products have health benefits and contribute to the distinctive taste of these vegetables. Because of their high value, there is a lot of interest in developing breeding strategies to increase the content of beneficial GSLs in Brassica species. GSLs are synthesized from certain amino acids and their biological roles depend largely on the structure of their side chains. Flavin-containing monooxygenase (FMO_GS-OX ) genes are involved in the synthesis of these side chains. To better understand GSL biosynthesis, we sequenced the transcriptomes of turnip (Brassica rapa var. rapa) tubers at four developmental stages (S1-S4) and determined their GSL content.

RESULTS

The total GSL content was high at the early stage (S1) of tuber development and increased up to S3, then decreased at S4. We detected 61 differentially expressed genes, including five FMO_GS-OX genes, that were related for GSL biosynthesis among the four developmental stages. Most of these genes were highly expressed at stages S1 to S3, but their expression was much lower at S4. We estimated the effect of the five FMO_GS-OX genes on GSL content by overexpressing them in turnip hairy roots and found that the amount of aliphatic GSLs increased significantly in the transgenic plants.

CONCLUSION

The transcriptome data and characterization of genes involved in GSL biosynthesis, particularly the FMO_GS-OX genes, will be valuable for improving the yield of beneficial GSLs in turnip and other Brassica crops. © 2019 Society of Chemical Industry.

Natural Variation of Glucosinolates and Their Breakdown Products in Broccoli (Brassica oleracea var. italica) Seeds

Seeds of 32 pure lines and 6 commercial broccoli cultivars were used to investigate variation in glucosinolates and their breakdown products. The aliphatic glucosinolate content was 54.5-218.7 μmol/g fresh weight, accounting for >90% of the total glucosinolates. The major glucosinolates found were glucoraphanin and glucoerucin in 27 samples and progoitrin in 7 samples. A gas chromatography-flame ionization detector (GC-FID) method was used to identify glucosinolate breakdown products; nine products were directly determined using standards. Using Arabidopsis thaliana lines myb28myb29 and Landsberg erecta to hydrolyze each reference glucosinolate, seven products were tentatively identified. 4-(Methylsulfinyl)butyl isothiocyanate and 5-(methylsulfinyl)pentanenitrile contents were 2.6-91.1 μmol/g fresh weight and 0-35.4 μmol/g fresh weight, respectively, with epithionitriles being more common than nitriles in accessions rich in alkenyl glucosinolate. Additionally, (S)-5-vinyl-1,3-oxazolidine-2-thione was detected in accessions rich in progoitrin. Specific lines with altered glucosinolate profiles and breakdown products were obtained and discussed according to the putative glucosinolate metabolism pathway.

Overexpression of the MYB29 transcription factor affects aliphatic glucosinolate synthesis in Brassica oleracea

KEY MESSAGE

Overexpression of BoMYB29 gene up-regulates the aliphatic glucosinolate pathway in Brassica oleracea plants increasing the production of the anti-cancer metabolite glucoraphanin, and the toxic and pungent sinigrin. Isothiocyanates, the bio-active hydrolysis products of glucosinolates, naturally produced by several Brassicaceae species, play an important role in human health and agriculture. This study aims at correlating the content of aliphatic glucosinolates to the expression of genes involved in their synthesis in Brassica oleracea, and perform functional analysis of BoMYB29 gene. To this purpose, three genotypes were used: a sprouting broccoli, a cabbage, and a wild genotype (Winspit), a high glucosinolate containing accession. Winspit showed the highest transcript level of BoMYB28, BoMYB29 and BoAOP2 genes, and BoAOP2 expression was positively correlated with that of the two MYB genes. Further analyses of the aliphatic glucosinolates also showed a positive correlation between the expression of BoAOP2 and the production of sinigrin and gluconapin in Winspit. The Winspit BoMYB29 CDS was cloned and overexpressed in Winspit and in the DH AG1012 line. Overexpressing Winspit plants produced higher quantities of alkenyl glucosinolates, such as sinigrin. Conversely, the DH AG1012 transformants showed a higher production of methylsulphinylalkyl glucosinolates, including glucoraphanin, and, despite an up-regulation of the aliphatic glucosinolate genes, no increase in alkenyl glucosinolates. The latter may be explained by the absence of a functional AOP2 gene in DH AG1012. Nevertheless, an extract of DH AG1012 lines overexpressing BoMYB29 provided a chemoprotective effect on human colon cells. This work exemplifies how the genetic diversity of B. oleracea may be used by breeders to select for higher expression of transcription factors for glucosinolate biosynthesis to improve its natural, health-promoting properties.

Glucosinolate variability between turnip organs during development

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

Genome-wide identification of loci affecting seed glucosinolate contents in Brassica napus L

Glucosinolates are amino acid-derived secondary metabolites that act as chemical defense agents against pests. However, the presence of high levels of glucosinolates severely diminishes the nutritional value of seed meals made from rapeseed (Brassica napus L.). To identify the loci affecting seed glucosinolate content (SGC), we conducted genome-wide resequencing in a population of 307 diverse B. napus accessions from the three B. napus ecotype groups, namely, spring, winter, and semi-winter. These resequencing data were used for a genome-wide association study (GWAS) to identify the loci affecting SGC. In the three ecotype groups, four common and four ecotype-specific haplotype blocks (HBs) were significantly associated with SGC. To identify candidate genes controlling SGC, transcriptome analysis was carried out in 36 accessions showing extreme SGC values. Analyses of haplotypes, genomic variation, and candidate gene expression pointed to five and three candidate genes in the common and spring group-specific HBs, respectively. Our expression analyses demonstrated that additive effects of the three candidate genes in the spring group-specific HB play important roles in the SGC of B. napus.

Enhancement of Glucosinolate Production in Watercress ( Nasturtium officinale) Hairy Roots by Overexpressing Cabbage Transcription Factors

Glucosinolates are secondary metabolites that play important roles in plant defense and human health, as their production in plants is enhanced by overexpressing transcription factors. Here, four cabbage transcription factors (IQD1-1, IQD1-2, MYB29-1, and MYB29-2) affecting genes in both aliphatic and indolic glucosinolates biosynthetic pathways and increasing glucosinolates accumulation were overexpressed in watercress. Five IQD1-1, six IQD1-2, five MYB29-1, six MYB29-2, and one GUS hairy root lines were created. The expression of all genes involved in glucosinolates biosynthesis was higher in transgenic lines than in the GUS hairy root line, in agreement with total glucosinolates contents, determined by high-performance liquid chromatography. In transgenic IQD1-1 (1), IQD1-2 (4), MYB29-1 (2), and MYB29-2 (1) hairy root lines, total glucosinolates were 3.39-, 3.04-, 2.58-, and 4.69-fold higher than those in the GUS hairy root lines, respectively. These results suggest a central regulatory function for IQD1-1, IQD1-2, MYB29-1, and MYB29-2 transcription factors in glucosinolates biosynthesis in watercress hairy roots.

Glucosinolates: Molecular structure, breakdown, genetic, bioavailability, properties and healthy and adverse effects

Glucosinolates are a large group of plant secondary metabolites with nutritional effects and biologically active compounds. Glucosinolates are mainly found in cruciferous plants such as Brassicaceae family, including common edible plants such as broccoli (Brassica oleracea var. italica), cabbage (B. oleracea var. capitata f. alba), cauliflower (B. oleracea var. botrytis), rapeseed (Brassica napus), mustard (Brassica nigra), and horseradish (Armoracia rusticana). If cruciferous plants are consumed without processing, myrosinase enzyme will hydrolyze the glucosinolates to various metabolites, such as isothiocyanates, nitriles, oxazolidine-2-thiones, and indole-3-carbinols. On the other hand, when cruciferous are cooked before consumption, myrosinase is inactivated and glucosinolates could be partially absorbed in their intact form through the gastrointestinal mucosa. This review paper summarizes the glucosinolate molecular breakdown, their genetic aspects from biosynthesis to precursors, their bioavailability (assimilation, absorption, and elimination of these molecules), their sensory properties, identified healthy and adverse effects, as well as the impact of processing on their bioavailability.

Effect of heat root stress and high salinity on glucosinolates metabolism in wild rocket

Wild rocket (Diplotaxis tenuifolia L.) is a leafy vegetable appreciated for its characteristic sensory properties which are mainly due to the presence of glucosinolates (GSLs). Short-term exposure to abiotic stresses can induce physiological responses and transcriptional changes which involve GSLs. For this reason, the aim of this work was to study the mechanisms of regulation of GSLs metabolism in rocket subjected to heat stress (40 °C) and high salinity (200 mM NaCl) imposed for up to 48 h. GSLs levels and the expression of methylthioalkylmalate synthase1 (DtMAM1), cytochromeP79F1 (DtCYP79F1), cytochromeP45083A1 (DtCYP83A1), cytosolic-sulfotransferase5b (DtST5b), cytosolic-sulfotransferase5c (DtST5c), flavinmono-oxygenase (DtFMO), myrosinase (DtMyro) and thio-methyl transferase (DtTMT) were analyzed under stress conditions. In addition, the effect on chlorophyll and glucose levels, as well as on chlorophyll a fluorescence were evaluated. Chlorophyll and chlorophyll fluorescence were not affected by the short-term application of stresses. Glucose levels in roots were doubled in response to high salinity, while, in the same organ, GSLs were three fold lower in response to both stresses. The relative content of several aliphatic GSLs was significantly reduced in leaves as a response to both stresses. A key role in GSLs metabolism and in the response to salinity is hypothesized for the gene DtTMT, as it showed an increment in transcripts accumulation (three-fold) consistent with the decrement in the GSLs levels found in salt-exposed leaves and roots. The results obtained in this study can be used in breeding programmes aiming to enhance rocket sensory quality and to improve the resistance to abiotic stresses.

Statistical modeling for estimating glucosinolate content in Chinese cabbage by growth conditions

BACKGROUND

Glucosinolate in Chinese cabbage (Brassica campestris L. ssp. pekinensis (Lour.) Rupr) has potential benefits for human health, and its content is affected by growth conditions. In this study, we used a statistical model to identify the relationship between glucosinolate content and growth conditions, and to predict glucosinolate content in Chinese cabbage.

RESULT

Multiple regression analysis was employed to develop the model's growth condition parameters of growing period, temperature, humidity and glucosinolate content measured in Chinese cabbage grown in a plant factory. The developed model was represented by a second-order multi-polynomial equation with two independent parameters: growth duration and temperature (adjusted R²  = 0.81), and accurately predicted glucosinolate content after 14 days of seeding.

CONCLUSION

To our knowledge, this study presents the first statistical model for evaluating glucosinolate content, suggesting a useful methodology for designing glucosinolate-related experiments, and optimizing glucosinolate content in Chinese cabbage cultivation. © 2018 Society of Chemical Industry.

Genotypic Variation of Glucosinolates and Their Breakdown Products in Leaves of Brassica rapa

An in-depth glucosinolate (GLS) profiling was performed on a core collection of 91 Brassica rapa accessions, representing diverse morphotypes of heterogeneous geographical origin, to better understand the natural variation in GLS accumulation and GLS breakdown product formation. Leaves of the 91 B. rapa accessions were analyzed for their GLS composition by UHPLC-DAD and the corresponding breakdown products by GC-MS. Fifteen different GLSs were identified, and aliphatic GLSs prevailed regarding diversity and concentration. Twenty-three GLS breakdown products were identified, among them nine isothiocyanates, ten nitriles, and four epithionitriles. Epithionitriles were the prevailing breakdown products due to the high abundance of alkenyl GLSs. The large scale data set allowed the identification of correlations in abundance of specific GLSs or of GLS breakdown products. Discriminant function analysis identified subspecies with high levels of similarity in the acquired metabolite profiles. In general, the five main subspecies grouped significantly in terms of their GLS profiles.

A naturally occurring variation in the BrMAM-3 gene is associated with aliphatic glucosinolate accumulation in Brassica rapa leaves

Glucosinolate profiles significantly vary among Brassica rapa genotypes. However, the molecular basis of these variations is largely unknown. In this study, we investigated a major quantitative trait locus (QTL) controlling aliphatic glucosinolate accumulation in B. rapa leaves. The QTL, which encompasses three tandem MAM genes and two MYB genes, was detected in two BC₂DH populations. Among the five-candidate genes, only the expression level of BrMAM-3 (Bra013007) was significantly correlated with the accumulation of aliphatic glucosinolates in B. rapa leaves. We identified a naturally occurring insertion within exon 1 of BrMAM-3, which is predicted to be a loss-of-function mutation, as confirmed by qRT-PCR. We determined that the loss of function was associated with the low glucosinolate content in B. rapa accessions. Furthermore, overexpressing the BrMAM-3 gene resulted in an increase in total aliphatic glucosinolates in Arabidopsis transgenic lines. Our study provides insights into the molecular mechanism underlying the accumulation of aliphatic glucosinolates in B. rapa leaves, thereby facilitating in the manipulation of total aliphatic glucosinolate content in Brassica crops.

Insights into the species-specific metabolic engineering of glucosinolates in radish (Raphanus sativus L.) based on comparative genomic analysis

Glucosinolates (GSLs) and their hydrolysis products present in Brassicales play important roles in plants against herbivores and pathogens as well as in the protection of human health. To elucidate the molecular mechanisms underlying the formation of species-specific GSLs and their hydrolysed products in Raphanus sativus L., we performed a comparative genomics analysis between R. sativus and Arabidopsis thaliana. In total, 144 GSL metabolism genes were identified, and most of these GSL genes have expanded through whole-genome and tandem duplication in R. sativus. Crucially, the differential expression of FMOGS-OX2 in the root and silique correlates with the differential distribution of major aliphatic GSL components in these organs. Moreover, MYB118 expression specifically in the silique suggests that aliphatic GSL accumulation occurs predominantly in seeds. Furthermore, the absence of the expression of a putative non-functional epithiospecifier (ESP) gene in any tissue and the nitrile-specifier (NSP) gene in roots facilitates the accumulation of distinctive beneficial isothiocyanates in R. sativus. Elucidating the evolution of the GSL metabolic pathway in R. sativus is important for fully understanding GSL metabolic engineering and the precise genetic improvement of GSL components and their catabolites in R. sativus and other Brassicaceae crops.

Expression profiles of BrMYB transcription factors related to glucosinolate biosynthesis and stress response in eight subspecies of Brassica rapa

Brassica rapa is a polyploid species with phenotypically diverse cultivated subspecies. Glucosinolates (GSLs) are secondary metabolites that contribute to anticarcinogenic activity and plant defense in Brassicaceae. Previously, complete coding sequences of 13 BrMYB transcription factors (TFs) related to GSL biosynthesis were identified in the B. rapa genome. In the present study, we investigated GSL content and expression levels of these BrMYBTFs in 38 accessions belonging to eight subspecies of B. rapa. Twelve identified GSLs were detected and were classified into three chemical groups based on patterns of GSL content and expression profiles of the BrMYBTFs. GSL content and BrMYBTF expression levels differed among genotypes, including B. rapa subspecies pekinensis, chinensis and rapa. BrMYB28.3, BrMYB51.1 and BrMYB122.2 positively regulated GSL content in 38 accessions. Furthermore, expression levels of BrMYB28s and BrMYB34.3 increased under most abiotic and biotic stress treatments. The three BrMYB51 paralogs also showed drastically increased expression levels after infection with Pectobacterium carotovorum. The results of the present study improve our understanding of the functional diversity of these 13 BrMYBTFs during the evolution of polyploid B. rapa.

Using RNA-Seq for Genomic Scaffold Placement, Correcting Assemblies, and Genetic Map Creation in a Common Brassica rapa Mapping Population

Brassica rapa is a model species for agronomic, ecological, evolutionary, and translational studies. Here, we describe high-density SNP discovery and genetic map construction for a B. rapa recombinant inbred line (RIL) population derived from field collected RNA sequencing (RNA-Seq) data. This high-density genotype data enables the detection and correction of putative genome misassemblies and accurate assignment of scaffold sequences to their likely genomic locations. These assembly improvements represent 7.1-8.0% of the annotated B. rapa genome. We demonstrate how using this new resource leads to a significant improvement for QTL analysis over the current low-density genetic map. Improvements are achieved by the increased mapping resolution and by having known genomic coordinates to anchor the markers for candidate gene discovery. These new molecular resources and improvements in the genome annotation will benefit the Brassicaceae genomics community and may help guide other communities in fine-tuning genome annotations.

RapaNet: A Web Tool for the Co-Expression Analysis of Brassica rapa Genes

Accumulated microarray data are used for assessing gene function by providing statistical values for co-expressed genes; however, only a limited number of Web tools are available for analyzing the co-expression of genes of Brassica rapa. We have developed a Web tool called RapaNet (http://bioinfo.mju.ac.kr/arraynet/brassica300k/query/), which is based on a data set of 143 B rapa microarrays compiled from various organs and at different developmental stages during exposure to biotic or abiotic stress. RapaNet visualizes correlated gene expression information via correlational networks and phylogenetic trees using Pearson correlation coefficient (r). In addition, RapaNet provides hierarchical clustering diagrams, scatterplots of log ratio intensities, related pathway maps, and cis-element lists of promoter regions. To ascertain the functionality of RapaNet, the correlated genes encoding ribosomal protein (L7Ae), photosystem II protein D1 (psbA), and cytochrome P450 monooxygenase in glucosinolate biosynthesis (CYP79F1) were retrieved from RapaNet and compared with their Arabidopsis homologues. An analysis of the co-expressed genes revealed their shared and unique features.

De novo transcriptome analysis and glucosinolate profiling in watercress (Nasturtium officinale R. Br.)

Background

Watercress (Nasturtium officinale R. Br.) is an aquatic herb species that is a rich source of secondary metabolites such as glucosinolates. Among these glucosinolates, watercress contains high amounts of gluconasturtiin (2-phenethyl glucosinolate) and its hydrolysis product, 2-phennethyl isothiocyanate, which plays a role in suppressing tumor growth. However, the use of N. officinale as a source of herbal medicines is currently limited due to insufficient genomic and physiological information.

Results

To acquire precise information on glucosinolate biosynthesis in N. officinale, we performed a comprehensive analysis of the transcriptome and metabolome of different organs of N. officinale. Transcriptome analysis of N. officinale seedlings yielded 69,570,892 raw reads. These reads were assembled into 69,635 transcripts, 64,876 of which were annotated to transcripts in public databases. On the basis of the functional annotation of N. officinale, we identified 33 candidate genes encoding enzymes related to glucosinolate biosynthetic pathways and analyzed the expression of these genes in the leaves, stems, roots, flowers, and seeds of N. officinale. The expression of NoMYB28 and NoMYB29, the main regulators of aliphatic glucosinolate biosynthesis, was highest in the stems, whereas the key regulators of indolic glucosinolate biosynthesis, such as NoDof1.1, NoMYB34, NoMYB51, and NoMYB122, were strongly expressed in the roots. Most glucosinolate biosynthetic genes were highly expressed in the flowers. HPLC analysis enabled us to detect eight glucosinolates in the different organs of N. officinale. Among these glucosinolates, the level of gluconasturtiin was considerably higher than any other glucosinolate in individual organs, and the amount of total glucosinolates was highest in the flower.

Conclusions

This study has enhanced our understanding of functional genomics of N. officinale, including the glucosinolate biosynthetic pathways of this plant. Ultimately, our data will be helpful for further research on watercress bio-engineering and better strategies for exploiting its anti-carcinogenic properties.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3792-5) contains supplementary material, which is available to authorized users.

De novo Transcriptome Assembly of Chinese Kale and Global Expression Analysis of Genes Involved in Glucosinolate Metabolism in Multiple Tissues

Chinese kale, a vegetable of the cruciferous family, is a popular crop in southern China and Southeast Asia due to its high glucosinolate content and nutritional qualities. However, there is little research on the molecular genetics and genes involved in glucosinolate metabolism and its regulation in Chinese kale. In this study, we sequenced and characterized the transcriptomes and expression profiles of genes expressed in 11 tissues of Chinese kale. A total of 216 million 150-bp clean reads were generated using RNA-sequencing technology. From the sequences, 98,180 unigenes were assembled for the whole plant, and 49,582~98,423 unigenes were assembled for each tissue. Blast analysis indicated that a total of 80,688 (82.18%) unigenes exhibited similarity to known proteins. The functional annotation and classification tools used in this study suggested that genes principally expressed in Chinese kale, were mostly involved in fundamental processes, such as cellular and molecular functions, the signal transduction, and biosynthesis of secondary metabolites. The expression levels of all unigenes were analyzed in various tissues of Chinese kale. A large number of candidate genes involved in glucosinolate metabolism and its regulation were identified, and the expression patterns of these genes were analyzed. We found that most of the genes involved in glucosinolate biosynthesis were highly expressed in the root, petiole, and in senescent leaves. The expression patterns of ten glucosinolate biosynthetic genes from RNA-seq were validated by quantitative RT-PCR in different tissues. These results provided an initial and global overview of Chinese kale gene functions and expression activities in different tissues.

Enriching Glucoraphanin in Brassica rapa Through Replacement of BrAOP2.2/BrAOP2.3 with Non-functional Genes

Sulforaphane, the hydrolytic product of glucoraphanin glucosinolate, is a potent anticarcinogen that reduces the risk of several human cancers. However, in most B. rapa vegetables, glucoraphanin is undetectable or only present in trace amounts, since the glucoraphanin that is present is converted to gluconapin by three functional BrAOP2 genes. In this study, to enrich beneficial glucoraphanin content in B. rapa, the functional BrAOP2 alleles were replaced by non-functional counterparts through marker-assisted backcrossing (MAB). We identified non-functional mutations of two BrAOP2 genes from B. rapa. The backcross progenies with introgression of both non-functional braop2.2 and braop2.3 alleles significantly increased the glucoraphanin content by 18 times relative to the recurrent parent. In contrast, replacement or introgression of single non-functional braop2.2 or braop2.3 locus did not change glucoraphanin content. Our results suggest that replacement of these two functional BrAOP2 genes with non-functional alleles has the potential for producing improved Brassica crops with enriched beneficial glucoraphanin content.

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.

How specialized volatiles respond to chronic and short-term physiological and shock heat stress in Brassica nigra

Brassicales release volatile glucosinolate breakdown products upon tissue mechanical damage, but it is unclear how the release of glucosinolate volatiles responds to abiotic stresses such as heat stress. We used three different heat treatments, simulating different dynamic temperature conditions in the field to gain insight into stress-dependent changes in volatile blends and photosynthetic characteristics in the annual herb Brassica nigra (L.) Koch. Heat stress was applied by either heating leaves through temperature response curve measurements from 20 to 40 °C (mild stress), exposing plants for 4 h to temperatures 25-44 °C (long-term stress) or shock-heating leaves to 45-50 °C. Photosynthetic reduction through temperature response curves was associated with decreased stomatal conductance, while the reduction due to long-term stress and collapse of photosynthetic activity after heat shock stress were associated with non-stomatal processes. Mild stress decreased constitutive monoterpene emissions, while long-term stress and shock stress resulted in emissions of the lipoxygenase pathway and glucosinolate volatiles. Glucosinolate volatile release was more strongly elicited by long-term stress and lipoxygenase product released by heat shock. These results demonstrate that glucosinolate volatiles constitute a major part of emission blend in heat-stressed B. nigra plants, especially upon chronic stress that leads to induction responses.

Characterization of the watercress (Nasturtium officinale R. Br.; Brassicaceae) transcriptome using RNASeq and identification of candidate genes for important phytonutrient traits linked to human health

BACKGROUND

Consuming watercress is thought to provide health benefits as a consequence of its phytonutrient composition. However, for watercress there are currently limited genetic resources underpinning breeding efforts for either yield or phytonutritional traits. In this paper, we use RNASeq data from twelve watercress accessions to characterize the transcriptome, perform candidate gene mining and conduct differential expression analysis for two key phytonutritional traits: antioxidant (AO) capacity and glucosinolate (GLS) content.

RESULTS

The watercress transcriptome was assembled to 80,800 transcripts (48,732 unigenes); 71 % of which were annotated based on orthology to Arabidopsis. Differential expression analysis comparing watercress accessions with 'high' and 'low' AO and GLS resulted in 145 and 94 differentially expressed loci for AO capacity and GLS respectively. Differentially expressed loci between high and low AO watercress were significantly enriched for genes involved in plant defence and response to stimuli, in line with the observation that AO are involved in plant stress-response. Differential expression between the high and low GLS watercress identified links to GLS regulation and also novel transcripts warranting further investigation. Additionally, we successfully identified watercress orthologs for Arabidopsis phenylpropanoid, GLS and shikimate biosynthesis pathway genes, and compiled a catalogue of polymorphic markers for future applications.

CONCLUSIONS

Our work describes the first transcriptome of watercress and establishes the foundation for further molecular study by providing valuable resources, including sequence data, annotated transcripts, candidate genes and markers.

De novo Transcriptome Analysis of Sinapis alba in Revealing the Glucosinolate and Phytochelatin Pathways

Sinapis alba is an important condiment crop and can also be used as a phytoremediation plant. Though it has important economic and agronomic values, sequence data, and the genetic tools are still rare in this plant. In the present study, a de novo transcriptome based on the transcriptions of leaves, stems, and roots was assembled for S. alba for the first time. The transcriptome contains 47,972 unigenes with a mean length of 1185 nt and an N50 of 1672 nt. Among these unigenes, 46,535 (97%) unigenes were annotated by at least one of the following databases: NCBI non-redundant (Nr), Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, Gene Ontology (GO), and Clusters of Orthologous Groups of proteins (COGs). The tissue expression pattern profiles revealed that 3489, 1361, and 8482 unigenes were predominantly expressed in the leaves, stems, and roots of S. alba, respectively. Genes predominantly expressed in the leaf were enriched in photosynthesis- and carbon fixation-related pathways. Genes predominantly expressed in the stem were enriched in not only pathways related to sugar, ether lipid, and amino acid metabolisms but also plant hormone signal transduction and circadian rhythm pathways, while the root-dominant genes were enriched in pathways related to lignin and cellulose syntheses, involved in plant-pathogen interactions, and potentially responsible for heavy metal chelating, and detoxification. Based on this transcriptome, 14,727 simple sequence repeats (SSRs) were identified, and 12,830 pairs of primers were developed for 2522 SSR-containing unigenes. Additionally, the glucosinolate (GSL) and phytochelatin metabolic pathways, which give the characteristic flavor and the heavy metal tolerance of this plant, were intensively analyzed. The genes of aliphatic GSLs pathway were predominantly expressed in roots. The absence of aliphatic GSLs in leaf tissues was due to the shutdown of BCAT4, MAM1, and CYP79F1 expressions. Glutathione was extensively converted into phytochelatin in roots, but it was actively converted to the oxidized form in leaves, indicating the different mechanisms in the two tissues. This transcriptome will not only benefit basic research and molecular breeding of S. alba but also be useful for the molecular-assisted transfer of beneficial traits to other crops.

BjuB.CYP79F1 Regulates Synthesis of Propyl Fraction of Aliphatic Glucosinolates in Oilseed Mustard Brassica juncea: Functional Validation through Genetic and Transgenic Approaches

Among the different types of methionine-derived aliphatic glucosinolates (GS), sinigrin (2-propenyl), the final product in 3C GS biosynthetic pathway is considered very important as it has many pharmacological and therapeutic properties. In Brassica species, the candidate gene regulating synthesis of 3C GS remains ambiguous. Earlier reports of GSL-PRO, an ortholog of Arabidopsis thaliana gene At1g18500 as a probable candidate gene responsible for 3C GS biosynthesis in B. napus and B. oleracea could not be validated in B. juncea through genetic analysis. In this communication, we report the isolation and characterization of the gene CYP79F1, an ortholog of A. thaliana gene At1g16410 that is involved in the first step of core GS biosynthesis. The gene CYP79F1 in B. juncea showed presence-absence polymorphism between lines Varuna that synthesizes sinigrin and Heera virtually free from sinigrin. Using this presence-absence polymorphism, CYP79F1 was mapped to the previously mapped 3C GS QTL region (J16Gsl4) in the LG B4 of B. juncea. In Heera, the gene was observed to be truncated due to an insertion of a ~4.7 kb TE like element leading to the loss of function of the gene. Functional validation of the gene was carried out through both genetic and transgenic approaches. An F2 population segregating only for the gene CYP79F1 and the sinigrin phenotype showed perfect co-segregation. Finally, genetic transformation of a B. juncea line (QTL-NIL J16Gsl4) having high seed GS but lacking sinigrin with the wild type CYP79F1 showed the synthesis of sinigrin validating the role of CYP79F1 in regulating the synthesis of 3C GS in B. juncea.

Identification of candidate genes involved in Witches' broom disease resistance in a segregating mapping population of Theobroma cacao L. in Brazil

Background

Witches’ broom disease (WBD) caused by the fungus Moniliophthora perniciosa is responsible for considerable economic losses for cacao producers. One of the ways to combat WBD is to plant resistant cultivars. Resistance may be governed by a few genetic factors, mainly found in wild germplasm.

Results

We developed a dense genetic linkage map with a length of 852.8 cM that contains 3,526 SNPs and is based on the MP01 mapping population, which counts 459 trees from a cross between the resistant ‘TSH 1188’ and the tolerant ‘CCN 51’ at the Mars Center for Cocoa Science in Barro Preto, Bahia, Brazil. Seven quantitative trait loci (QTL) that are associated with WBD were identified on five different chromosomes using a multi-trait QTL analysis for outbreeders. Phasing of the haplotypes at the major QTL region on chromosome IX on a diversity panel of genotypes clearly indicates that the major resistance locus comes from a well-known source of WBD resistance, the clone ‘SCAVINA 6’. Various potential candidate genes identified within all QTL may be involved in different steps leading to disease resistance. Preliminary expression data indicate that at least three of these candidate genes may play a role during the first 12 h after infection, with clear differences between ‘CCN 51’ and ‘TSH 1188’.

Conclusions

We combined the information from a large mapping population with very distinct parents that segregate for WBD, a dense set of mapped markers, rigorous phenotyping capabilities and the availability of a sequenced genome to identify several genomic regions that are involved in WBD resistance. We also identified a novel source of resistance that most likely comes from the ‘CCN 51’ parent. Thanks to the large population size of the MP01 population, we were able to pick up QTL and markers with relatively small effects that can contribute to the creation and selection of more tolerant/resistant plant material.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2415-x) contains supplementary material, which is available to authorized users.