Further insight into decreases in seed glucosinolate content based on QTL mapping and RNA-seq in Brassica napus L

The QTL hotspots determining seed glucosinolate content instead of only four HAG1 loci and elucidation of a potential regulatory model for rapeseed SGC variation. Glucosinolates (GSLs) are amino acid-derived, sulfur-rich secondary metabolites that function as biopesticides and flavor compounds, but the high seed glucosinolate content (SGC) reduces seed quality for rapeseed meal. To dissect the genetic mechanism and further reduce SGC in rapeseed, QTL mapping was performed using an updated high-density genetic map based on a doubled haploid (DH) population derived from two parents that showed significant differences in SGC. In 15 environments, a total of 162 significant QTLs were identified for SGC and then integrated into 59 consensus QTLs, of which 32 were novel QTLs. Four QTL hotspot regions (QTL-HRs) for SGC variation were discovered on chromosomes A09, C02, C07 and C09, including seven major QTLs that have previously been reported and four novel major QTLs in addition to HAG1 loci. SGC was largely determined by superimposition of advantage allele in the four QTL-HRs. Important candidate genes directly related to GSL pathways were identified underlying the four QTL-HRs, including BnaC09.MYB28, BnaA09.APK1, BnaC09.SUR1 and BnaC02.GTR2a. Related differentially expressed candidates identified in the minor but environment stable QTLs indicated that sulfur assimilation plays an important rather than dominant role in SGC variation. A potential regulatory model for rapeseed SGC variation constructed by combining candidate GSL gene identification and differentially expressed gene analysis based on RNA-seq contributed to a better understanding of the GSL accumulation mechanism. This study provides insights to further understand the genetic regulatory mechanism of GSLs, as well as the potential loci and a new route to further diminish the SGC in rapeseed.

Metabolite analysis of Arabidopsis CYP79A2 overexpression lines reveals turnover of benzyl glucosinolate and an additive effect of different aldoximes on phenylpropanoid repression

Indole-3-acetaldoxime (IAOx) and phenylacetaldoxime (PAOx) are precursors for the growth hormones indole-3-acetic acid (IAA) and phenylacetic acid (PAA) and the defense compounds glucosinolates in Brassicales. Our recent work has shown that Arabidopsis transgenic lines overexpressing AtCYP79A2, a PAOx-production enzyme, accumulate the PAOx-derived compounds benzyl glucosinolate and PAA. Here we report that they also accumulate the benzyl glucosinolate hydrolysis products benzyl isothiocyanate and benzyl cyanide, which indicates that the turnover of benzyl glucosinolate can occur in intact tissues. Myrosinases or β-glucosidases are known to catalyze glucosinolate breakdown. However, transcriptomics analysis detected no substantial increase in expression of known myrosinases or putative β-glucosidases in AtCYP79A2 overexpressing lines. It was previously shown that accumulation of aldoximes or their derivatives represses the phenylpropanoid pathway. For instance, ref2 mutant having a defect in one of the aldoxime catabolic enzymes decreases phenylpropanoid production. Considering that AtCYP79A2 is not expressed in most organs under optimal growth condition, ref2 accumulates aliphatic aldoximes but not PAOx. Interestingly, overexpression of AtCYP79A2 in ref2 resulted in a further decrease in sinapoylmalate content compared to ref2. This indicates that accumulation of PAOx has an additive effect on phenylpropanoid pathway suppression mediated by other aldoximes.

Transcriptomics analysis of genes induced by melatonin related to glucosinolates synthesis in broccoli hairy roots

Glucoraphanin (GRA) is found in the seeds and vegetative organs of broccoli (Brassica oleracea L. var. italica Planch) as the precursor of anti-carcinogen sulforaphane (SF). The yield of GRA obtained from these materials is weak and the cost is high. In recent years, the production of plant secondary metabolites by large-scale hairy roots culture in vitro has succeeded in some species. Melatonin (MT) is a natural hormone which existed in numerous organisms. Studies have demonstrated that MT can improve the synthesis of secondary metabolites in plants. At present, it has not been reported that MT regulates the biosynthesis of glucoraphanin in broccoli hairy roots. In this study, the broccoli hairy roots that grew for 20 d were respectively treated by 500 µM MT for 0, 6, 12, 20 and 32. To explore the reason of changes in secondary metabolites and reveal the biosynthetic pathway of glucoraphanin at transcriptional level. Compared with 0 h, the yield of GRA under other treatments was increased, and the overall trend was firstly increased and then decreased. The total yield of GRA reached the highest at 12 h, which was 1.22-fold of 0 h. Then, the genome of broccoli as the reference, a total of 13234 differentially expressed genes (DEGs) were identified in broccoli hairy roots under treatment with 500 µM MT for 0, 6, 12, 20 and 32 h, respectively. Among these DEGs, 6266 (47.35%) were upregulated and 6968 (52.65%) were downregulated. It was found that the pathway of 'Glucosinolates biosynthesis (ko00966)' was enriched in the 16th place by Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of the upregulated DEGs. The expression of key genes in the GRA biosynthesis pathway was upregulated at all time points, and a deduced GRA biosynthesis pathway map was constructed for reference.

Arabidopsis assemble distinct root-associated microbiomes through the synthesis of an array of defense metabolites

Plant associated microbiomes are known to confer fitness advantages to the host. Understanding how plant factors including biochemical traits influence host associated microbiome assembly could facilitate the development of microbiome-mediated solutions for sustainable plant production. Here, we examined microbial community structures of a set of well-characterized Arabidopsis thaliana mutants disrupted in metabolic pathways for the production of glucosinolates, flavonoids, or a number of defense signalling molecules. A. thaliana lines were grown in a natural soil and maintained under greenhouse conditions for 4 weeks before collection of roots for bacterial and fungal community profiling. We found distinct relative abundances and diversities of bacterial and fungal communities assembled in the individual A. thaliana mutants compared to their parental lines. Bacterial and fungal genera were mostly enriched than depleted in secondary metabolite and defense signaling mutants, except for flavonoid mutations on fungi communities. Bacterial genera Azospirillum and Flavobacterium were significantly enriched in most of the glucosinolate, flavonoid and signalling mutants while the fungal taxa Sporobolomyces and Emericellopsis were enriched in several glucosinolates and signalling mutants. Whilst the present study revealed marked differences in microbiomes of Arabidopsis mutants and their parental lines, it is suggestive that unknown enzymatic and pleiotropic activities of the mutated genes could contribute to the identified host-associated microbiomes. Notwithstanding, this study revealed interesting gene-microbiota links, and thus represents valuable resource data for selecting candidate A. thaliana mutants for analyzing the links between host genetics and the associated microbiome.

The BrGI Circadian Clock Gene Is Involved in the Regulation of Glucosinolates in Chinese Cabbage

Circadian clocks integrate environmental cues with endogenous signals to coordinate physiological outputs. Clock genes in plants are involved in many physiological and developmental processes, such as photosynthesis, stomata opening, stem elongation, light signaling, and floral induction. Many Brassicaceae family plants, including Chinese cabbage (Brassica rapa ssp. pekinensis), produce a unique glucosinolate (GSL) secondary metabolite, which enhances plant protection, facilitates the design of functional foods, and has potential medical applications (e.g., as antidiabetic and anticancer agents). The levels of GSLs change diurnally, suggesting a connection to the circadian clock system. We investigated whether circadian clock genes affect the biosynthesis of GSLs in Brassica rapa using RNAi-mediated suppressed transgenic Brassica rapa GIGENTEA homolog (BrGI knockdown; hereafter GK1) Chinese cabbage. GIGANTEA plays an important role in the plant circadian clock system and is related to various developmental and metabolic processes. Using a validated GK1 transgenic line, we performed RNA sequencing and high-performance liquid chromatography analyses. The transcript levels of many GSL pathway genes were significantly altered in GK1 transgenic plants. In addition, GSL contents were substantially reduced in GK1 transgenic plants. We report that the BrGI circadian clock gene is required for the biosynthesis of GSLs in Chinese cabbage plants.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Influence of Genotype on High Glucosinolate Synthesis Lines of Brassica rapa

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

The Arabidopsis Iron-Sulfur (Fe-S) Cluster Gene MFDX1 Plays a Role in Host and Nonhost Disease Resistance by Accumulation of Defense-Related Metabolites

Until recently, genes from the iron-sulfur (Fe-S) cluster pathway were not known to have a role in plant disease resistance. The Nitrogen Fixation S (NIFS)-like 1 (NFS1) and Mitochondrial Ferredoxin-1 (MFDX1) genes are part of a set of 27 Fe-S cluster genes induced after infection with host and nonhost pathogens in Arabidopsis. A role for AtNFS1 in plant immunity was recently demonstrated. In this work, we showed that MFDX1 is also involved in plant defense. More specifically, Arabidopsis mfdx1 mutants were compromised for nonhost resistance against Pseudomonas syringae pv. tabaci, and showed increased susceptibility to the host pathogen P. syringae pv. tomato DC3000. Arabidopsis AtMFDX1 overexpression lines were less susceptible to P. syringae pv. tomato DC3000. Metabolic profiling revealed a reduction of several defense-related primary and secondary metabolites, such as asparagine and glucosinolates in the Arabidopsis mfdx1-1 mutant when compared to Col-0. A reduction of 5-oxoproline and ornithine metabolites that are involved in proline synthesis in mitochondria and affect abiotic stresses was also observed in the mfdx1-1 mutant. In contrast, an accumulation of defense-related metabolites such as glucosinolates was observed in the Arabidopsis NFS1 overexpressor when compared to wild-type Col-0. Additionally, mfdx1-1 plants displayed shorter primary root length and reduced number of lateral roots compared to the Col-0. Taken together, these results provide additional evidence for a new role of Fe-S cluster pathway in plant defense responses.

LATE ELONGATED HYPOCOTYL potentiates resistance conferred by CIRCADIAN CLOCK ASSOCIATED1 to aphid by co-regulating the expression of indole glucosinolate biosynthetic genes

CIRCADIAN CLOCK ASSOCIATED1 (CCA1) and LATE ELONGATED HYPOCOTYL (LHY) are core components of the circadian clock in Arabidopsis thaliana that impacts plant response to biotic stresses. Their clock-regulating functions are believed to be partially redundant, and mutation of either gene leads to shortened periods of the circadian cycle. Our recent study has demonstrated that CCA1 promotes plant resistance to the green peach aphid (Myzus persicae) through modulation of indole glucosinolate biosynthesis, but the role of LHY remains to be elucidated. Here we showed that, similar to cca1-11, single mutant lhy-21 became more susceptible to aphid infestation. Damage to the cca1-11 lhy-21 double mutant by aphids was most pronounced, indicating that the defensive roles of CCA1 and LHY were not entirely redundant. Also, the cyclic expression pattern of key indole glucosinolate biosynthetic genes was considerably disturbed in both single mutants and this was more severe in the double mutant. Apparently, both CCA1 and LHY were necessary for circadian-regulated indole glucosinolate biosynthesis. Taken together, LHY-CCA1 coordination in transcriptional regulation of indole glucosinolate biosynthetic genes most likely contributed to plant defensive capacity against aphids.

RcTGA1 and glucosinolate biosynthesis pathway involvement in the defence of rose against the necrotrophic fungus Botrytis cinerea

BACKGROUND

Rose is an important economic crop in horticulture. However, its field growth and postharvest quality are negatively affected by grey mould disease caused by Botrytis c. However, it is unclear how rose plants defend themselves against this fungal pathogen. Here, we used transcriptomic, metabolomic and VIGS analyses to explore the mechanism of resistance to Botrytis c.

RESULT

In this study, a protein activity analysis revealed a significant increase in defence enzyme activities in infected plants. RNA-Seq of plants infected for 0 h, 36 h, 60 h and 72 h produced a total of 54 GB of clean reads. Among these reads, 3990, 5995 and 8683 differentially expressed genes (DEGs) were found in CK vs. T36, CK vs. T60 and CK vs. T72, respectively. Gene annotation and cluster analysis of the DEGs revealed a variety of defence responses to Botrytis c. infection, including resistance (R) proteins, MAPK cascade reactions, plant hormone signal transduction pathways, plant-pathogen interaction pathways, Ca²⁺ and disease resistance-related genes. qPCR verification showed the reliability of the transcriptome data. The PTRV2-RcTGA1-infected plant material showed improved susceptibility of rose to Botrytis c. A total of 635 metabolites were detected in all samples, which could be divided into 29 groups. Metabonomic data showed that a total of 59, 78 and 74 DEMs were obtained for T36, T60 and T72 (T36: Botrytis c. inoculated rose flowers at 36 h; T60: Botrytis c. inoculated rose flowers at 60 h; T72: Botrytis c. inoculated rose flowers at 72 h) compared to CK, respectively. A variety of secondary metabolites are related to biological disease resistance, including tannins, amino acids and derivatives, and alkaloids, among others; they were significantly increased and enriched in phenylpropanoid biosynthesis, glucosinolates and other disease resistance pathways. This study provides a theoretical basis for breeding new cultivars that are resistant to Botrytis c.

CONCLUSION

Fifty-four GB of clean reads were generated through RNA-Seq. R proteins, ROS signalling, Ca²⁺ signalling, MAPK signalling, and SA signalling were activated in the Old Blush response to Botrytis c. RcTGA1 positively regulates rose resistance to Botrytis c. A total of 635 metabolites were detected in all samples. DEMs were enriched in phenylpropanoid biosynthesis, glucosinolates and other disease resistance pathways.

Plant surface metabolites as potent antifungal agents

Triunsaturated fatty acids are substrates for the synthesis of the defense hormone jasmonate which plays roles in resistance to numerous fungal pathogens. However, relatively little is known about other potential roles of di-unsaturated and triunsaturated fatty acids in resistance to fungal pathogens - in particular those that can attack plants at the seedling stage. We examined the roles of polyunsaturated fatty acids (PUFAs) in Arabidopsis thaliana during attack by the necrotrophic pathogen, Botrytis cinerea. We found that PUFA-deficient Arabidopsis mutants (fad2-1, fad2-3 and fad3-2 fad7-2 fad8 [fad trip]) displayed an unexpectedly strong resistance to B. cinerea at the cotyledon stage. Preliminary analyses revealed no changes in the expression of defense genes, however cuticle permeability defects were detected in both fad2-1 and fad trip mutants. Analysis of B. cinerea development on the surface of cotyledons revealed arrested hyphal growth on fad2-3 and fad trip mutants and 28% reduction in fungal adhesion on fad2-3 cotyledons. Surface metabolite analysis from the cotyledons of PUFA mutants led to the identification of 7-methylsulfonylheptyl glucosinolate (7MSOHG), which over-accumulated on the plant surface. We linked the appearance of 7MSOHG to defects in cuticle composition and permeability of mutants and show that its appearance correlates with resistance to B. cinerea.

Epistatic Transcription Factor Networks Differentially Modulate Arabidopsis Growth and Defense

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

Fine mapping and candidate gene analysis of a seed glucosinolate content QTL, qGSL-C2, in rapeseed (Brassica napus L.)

QTL mapping and candidate gene analysis indicate that allelic variations in BnaC2.MYB28 resulted from homeologous exchange and determine difference in seed glucosinolate content. A low seed glucosinolate content has long been an important breeding objective in rapeseed improvement. However, the molecular mechanisms underlying seed GSL content variations remain to be elucidated in allotetraploid Brassica napus. Here, we developed a double haploid population from a cross between two B. napus accessions that possess relatively low, but significantly different seed GSL contents and identified a major QTL, qGSL-C2, on chromosome C02 that explains 30.88-72.87% of the phenotypic variation observed in five environments. Using near-isogenic lines, we further delimited qGSL-C2 to a physical region of 49 kb on the B. rapa chromosome A02 which is highly homologous to the target C02 interval. Among five candidate genes, BnaC2.MYB28, a homologue of the Arabidopsis MYB28 encoding a putative R2R3-MYB-type transcription factor functioning in aliphatic methionine-derived GSL synthesis, was most likely to be the target gene underlying the QTL. Sequence analysis revealed multiple insertion/deletion and SNP variations in the genomic region between the alleles of the NILs. Furthermore, the allelic variations in BnaC2.MYB28 in the natural B. napus population were significantly associated with seed GSL content. Remarkably, the phylogenetic analysis and sequence comparison suggested that while the BnaC2.MYB28 allele from the parental line G120 was inherited from B. oleracea BolC2.MYB28, its counterpart from the other parent, 9172, most likely evolved from B. rapa BraA2.MYB28 via possible homeologous exchange. Our study promotes greater understanding of the molecular regulation of seed GSL content and provides useful molecular markers for seed GSL improvement in B. napus.

Differential expression of major genes involved in the biosynthesis of aliphatic glucosinolates in intergeneric Baemoochae (Brassicaceae) and its parents during development

Thus study found the temporal and spatial relationship between production of aliphatic glucosinolate compounds and the expression profile of glucosinolate-related genes during growth and development in radish, Chinese cabbage, and their intergeneric hybrid baemoochae plants. Glucosinolates (GSLs) are one of major bioactive compounds in Brassicaceae plants. GSLs play a role in defense against microbes as well as chemo-preventative activity against cancer, which draw attentions from plant scientists. We investigated the temporal relationship between production of aliphatic Glucosinolate (GSLs) compounds and the expression profile of GSL related genes during growth and development in radish, Chinese cabbage, and their intergeneric hybrid, baemoochae. Over the complete life cycle, Glucoraphasatin (GRH) and glucoraphanin (GRE) predominated in radish, whereas gluconapin (GNP), glucobrassicanapin (GBN), and glucoraphanin (GRA) abounded in Chinese cabbage. Baemoochae contained intermediate levels of all GSLs studied, indicating inheritance from both radish and Chinese cabbage. Expression patterns of BCAT4, CYP79F1, CYP83A1, UGT74B1, GRS1, FMOgs-ox1, and AOP2 genes showed a correlation to their corresponding encoded proteins in radish, Chinese cabbage, and baemoochae. Interestingly, there is a sharp change in gene expression pattern involved in side chain modification, particularly GRS1, FMOgs-ox1, and AOP2, among these plants during the vegetative and reproductive stage. For instance, the GRS1 was strongly expressed during leaf development, while both of FMOgs-ox1 and AOP2 was manifested high in floral tissues. Furthermore, expression of GRS1 gene which is responsible for GRH production was predominantly expressed in leaf tissues of radish and baemoochae, whereas it was only slightly detected in Chinese cabbage root tissue, explaining why radish has an abundance of GRH compared to other Brassica plants. Altogether, our comprehensive and comparative data proved that aliphatic GSLs biosynthesis is dynamically and precisely regulated in a tissue- and development-dependent manner in Brassicaceae family members.

Genetic architecture of glucosinolate variation in Brassica napus

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

Differential regulation of host plant adaptive genes in Pieris butterflies exposed to a range of glucosinolate profiles in their host plants

Specialist herbivores have often evolved highly sophisticated mechanisms to counteract defenses mediated by major plant secondary-metabolites. Plant species of the herbivore host range often display high chemical diversity and it is not well understood how specialist herbivores respond to this chemical diversity. Pieris larvae overcome toxic products from glucosinolate hydrolysis, the major chemical defense of their Brassicaceae hosts, by expressing nitrile-specifier proteins (NSP) in their gut. Furthermore, Pieris butterflies possess so-called major allergen (MA) proteins, which are multi-domain variants of a single domain major allergen (SDMA) protein expressed in the guts of Lepidopteran larvae. Here we show that Pieris larvae fine-tune NSP and MA gene expression depending on the glucosinolate profiles of their Brassicaceae hosts. Although the role of MA is not yet fully understood, the expression levels of NSP and MA in larvae that fed on plants whose glucosinolate composition varied was dramatically changed, whereas levels of SDMA expression remained unchanged. In addition, we found a similar regulation pattern among these genes in larvae feeding on Arabidopsis mutants with different glucosinolate profiles. Our results demonstrate that Pieris larvae appear to use different host plant adaptive genes to overcome a wide range of glucosinolate profiles in their host plants.

Altered Glucosinolate Profiles and Expression of Glucosinolate Biosynthesis Genes in Ringspot-Resistant and Susceptible Cabbage Lines

Ringspot, caused by the fungus Mycosphaerella brassicicola, is a serious disease of Brassica crops worldwide. Despite noteworthy progress to reveal the role of glucosinolates in pathogen defense, the host⁻pathogen interaction between cabbage (Brassica oleracea) and M. brassicicola has not been fully explored. Here, we investigated the glucosinolate profiles and expression of glucosinolate biosynthesis genes in the ringspot-resistant (R) and susceptible (S) lines of cabbage after infection with M. brassicicola. The concomitant rise of aliphatic glucoiberverin (GIV) and indolic glucobrassicin (GBS) and methoxyglucobrassicin (MGBS) was linked with ringspot resistance in cabbage. Pearson's correlation and principle component analysis showed a significant positive association between GIV contents and the expression of the glucosinolate biosynthesis gene ST5b-Bol026202 and between GBS contents and the expression of the glucosinolate biosynthesis gene MYB34-Bol017062. Our results confirmed that M. brassicicola infection induces the expression of glucosinolate biosynthesis genes in cabbage, which alters the content of individual glucosinolates. This link between the expression of glucosinolate biosynthesis genes and the accumulation of their respective glucosinolates with the resistance to ringspot extends our molecular sense of glucosinolate-negotiated defense against M. brassicicola in cabbage.

Brassinosteroids regulate glucosinolate biosynthesis in Arabidopsis thaliana

Plants must constantly adjust their growth and defense responses to deal with the wide variety of stresses they encounter in their environment. Among phytohormones, brassinosteroids (BRs) are an important group of plant steroid hormones involved in numerous aspects of the plant lifecycle including growth, development and responses to various stresses including insect attacks. Here, we show that BRs regulate glucosinolate (GS) biosynthesis and function in insect herbivory. Preference tests and larval feeding experiments using the generalist herbivore, diamondback moth (Plutella xylostella), revealed that the larvae prefer to feed on Arabidopsis thaliana brassinosteroid insensitive 1 (bri1-5) plants over wild-type Ws-2 or BRI1-Flag (bri1-5 background) transgenic plants, which results in an increase in larval weight. Analysis of GS contents showed that 3-(methylsulfinyl) propyl GS (C3) levels were higher in bri1-5 than in Ws2 and BRI1-Flag transgenic plants, whereas sinigrin (2-propenylglucosinolate), glucoerucin (4-methylthiobutylglucosinolate) and glucobrassicin (indol-3-ylmethylglucosinolate) levels were lower in this mutant. We investigated the effect of brassinolide (BL) on GS biosynthesis in Arabidopsis and radish (Raphanus sativus L.) by monitoring the expression levels of GS biosynthetic genes, including MAM1, MAM3, BCAT4 and AOP2, which increased in a BL-dependent manner. These results suggest that BRs regulate GS profiles in higher plants, which function in defense responses against insects.

Differences in the enzymatic hydrolysis of glucosinolates increase the defense metabolite diversity in 19 Arabidopsis thaliana accessions

Plants of the order Brassicales produce glucosinolates (GS), a group of secondary metabolites that are part of an elaborate defense system. But it is not the GS itself rather its enzymatic hydrolysis products that cause the bioactive effects protecting the plants against pests and pathogens. Thus the enzymatic hydrolysis and a variety of additional influential factors determine the structural outcome of the GS degradation process. To evaluate the possible diversity of defense metabolites a range of 19 Arabidopsis thaliana accessions were selected showing divergence in their geographical origin, in their phenotype, and in their GS profile. These particular accessions accumulate several alkenyl GS, hydroxyalkyl GS, methylthioalkyl GS, and methylsulfinylalkyl GS in their rosette leaves whereas the indole GS contents are relatively invariant, as analyzed by UHPLC-DAD. After tissue disruption the enzymatic formation of GS hydrolysis products was examined and breakdown products were identified and quantified by GC-MS. Great differences in the amount and structure of volatile enzymatic degradation products could be observed in the different accessions, with strong variation in formation of epithionitriles, nitriles, and isothiocyanates. The occurrence of specific GS hydrolysis products was put in relation to relative gene expression profiles of myrosinases and specifier proteins as measured by RT-qPCR, and in relation to relative protein abundance of epithiospecifier protein. Dependent on the different GS profiles and reliant on degradation protein abundance and composition the ecotypes strongly varied in their ability to form isothiocyanates, nitriles and epithionitriles, thus increasing the plants' equipment of defense metabolites.

Chemopreventive glucosinolate accumulation in various broccoli and collard tissues: Microfluidic-based targeted transcriptomics for by-product valorization

Floret, leaf, and root tissues were harvested from broccoli and collard cultivars and extracted to determine their glucosinolate and hydrolysis product profiles using high performance liquid chromatography and gas chromotography. Quinone reductase inducing bioactivity, an estimate of anti-cancer chemopreventive potential, of the extracts was measured using a hepa1c1c7 murine cell line. Extracts from root tissues were significantly different from other tissues and contained high levels of gluconasturtiin and glucoerucin. Targeted gene expression analysis on glucosinolate biosynthesis revealed that broccoli root tissue has elevated gene expression of AOP2 and low expression of FMOGS-OX homologs, essentially the opposite of what was observed in broccoli florets, which accumulated high levels of glucoraphanin. Broccoli floret tissue has significantly higher nitrile formation (%) and epithionitrile specifier protein gene expression than other tissues. This study provides basic information of the glucosinolate metabolome and transcriptome for various tissues of Brassica oleracea that maybe utilized as potential byproducts for the nutraceutical market.

Understanding of MYB Transcription Factors Involved in Glucosinolate Biosynthesis in Brassicaceae

Glucosinolates (GSLs) are widely known secondary metabolites that have anticarcinogenic and antioxidative activities in humans and defense roles in plants of the Brassicaceae family. Some R2R3-type MYB (myeloblastosis) transcription factors (TFs) control GSL biosynthesis in Arabidopsis. However, studies on the MYB TFs involved in GSL biosynthesis in Brassica species are limited because of the complexity of the genome, which includes an increased number of paralog genes as a result of genome duplication. The recent completion of the genome sequencing of the Brassica species permits the identification of MYB TFs involved in GSL biosynthesis by comparative genome analysis with A. thaliana. In this review, we describe various findings on the regulation of GSL biosynthesis in Brassicaceae. Furthermore, we identify 63 orthologous copies corresponding to five MYB TFs from Arabidopsis, except MYB76 in Brassica species. Fifty-five MYB TFs from the Brassica species possess a conserved amino acid sequence in their R2R3 MYB DNA-binding domain, and share close evolutionary relationships. Our analysis will provide useful information on the 55 MYB TFs involved in the regulation of GSL biosynthesis in Brassica species, which have a polyploid genome.

Epistasis × environment interactions among Arabidopsis thaliana glucosinolate genes impact complex traits and fitness in the field

Despite the growing number of studies showing that genotype × environment and epistatic interactions control fitness, the influences of epistasis × environment interactions on adaptive trait evolution remain largely uncharacterized. Across three field trials, we quantified aliphatic glucosinolate (GSL) defense chemistry, leaf damage, and relative fitness using mutant lines of Arabidopsis thaliana varying at pairs of causal aliphatic GSL defense genes to test the impact of epistatic and epistasis × environment interactions on adaptive trait variation. We found that aliphatic GSL accumulation was primarily influenced by additive and epistatic genetic variation, leaf damage was primarily influenced by environmental variation and relative fitness was primarily influenced by epistasis and epistasis × environment interactions. Epistasis × environment interactions accounted for up to 48% of the relative fitness variation in the field. At a single field site, the impact of epistasis on relative fitness varied significantly over 2 yr, showing that epistasis × environment interactions within a location can be temporally dynamic. These results suggest that the environmental dependency of epistasis can profoundly influence the response to selection, shaping the adaptive trajectories of natural populations in complex ways, and deserves further consideration in future evolutionary studies.

How does a plant orchestrate defense in time and space? Using glucosinolates in Arabidopsis as case study

The sessile nature of plants has caused plants to develop means to defend themselves against attacking organisms. Multiple strategies range from physical barriers to chemical warfare including pre-formed anticipins as well as phytoalexins produced only upon attack. While phytoalexins require rapid induction, constitutive defenses can impose ecological costs if they deter pollinators or attract specialized herbivores. In the model Arabidopsis thaliana, the well-characterized glucosinolate anticipins are categorized into different classes, aliphatic and indole glucosinolates, depending on their amino acid precursor. Using glucosinolates in Arabidopsis as case study, we will discuss how plants orchestrate synthesis, storage and activation of pre-formed defense compounds spatially and temporally.

A novel Filamentous Flower mutant suppresses brevipedicellus developmental defects and modulates glucosinolate and auxin levels

BREVIPEDICELLUS (BP) encodes a class-I KNOTTED1-like homeobox (KNOX) transcription factor that plays a critical role in conditioning a replication competent state in the apical meristem, and it also governs growth and cellular differentiation in internodes and pedicels. To search for factors that modify BP signaling, we conducted a suppressor screen on bp er (erecta) plants and identified a mutant that ameliorates many of the pleiotropic defects of the parent line. Map based cloning and complementation studies revealed that the defect lies in the FILAMENTOUS FLOWER (FIL) gene, a member of the YABBY family of transcriptional regulators that contribute to meristem organization and function, phyllotaxy, leaf and floral organ growth and polarity, and are also known to repress KNOX gene expression. Genetic and cytological analyses of the fil-10 suppressor line indicate that the role of FIL in promoting growth is independent of its previously characterized influences on meristem identity and lateral organ polarity, and likely occurs non-cell-autonomously from superior floral organs. Transcription profiling of inflorescences revealed that FIL downregulates numerous transcription factors which in turn may subordinately regulate inflorescence architecture. In addition, FIL, directly or indirectly, activates over a dozen genes involved in glucosinolate production in part by activating MYB28, a known activator of many aliphatic glucosinolate biosynthesis genes. In the bp er fil-10 suppressor mutant background, enhanced expression of CYP71A13, AMIDASE1 (AMI) and NITRILASE genes suggest that auxin levels can be modulated by shunting glucosinolate metabolites into the IAA biosynthetic pathway, and increased IAA levels in the bp er fil-10 suppressor accompany enhanced internode and pedicel elongation. We propose that FIL acts to oppose KNOX1 gene function through a complex regulatory network that involves changes in secondary metabolites and auxin.

Mild osmotic stress promotes 4-methoxy indolyl-3-methyl glucosinolate biosynthesis mediated by the MKK9-MPK3/MPK6 cascade in Arabidopsis

MKK9-MPK3/MPK6 cascade positively regulates IGSs' biosynthetic genes. Glucosinolates (GSs), secondary metabolites well known for their roles in plant defense, have been implicated to play an important role in plant abiotic stress response; however, the exact role in these processes and the underlying regulatory mechanisms remain elusive. Mitogen-activated protein kinase (MAPK) cascades are extensively involved in plant abiotic stress response. In this study, we examined the levels of four indolic glucosinolates (IGSs) in the shoots of Arabidopsis seedlings under mild osmotic stress conditions and found that 4-methoxy indolyl-3-methyl glucosinolate (4MI3G) accumulated and that MPK3 and MPK6 were activated. Loss of MPK3 or MPK6 function led to a reduction in mild osmotic stress-induced 4MI3G. Further analyses revealed that MKK9 acts upstream of MPK3 and MPK6 to promote 4MI3G accumulation. The level of 4MI3G induced by mild osmotic stress was reduced in the mkk9 mutant. Conversely, 4MI3G increased in MKK9 ^DD , a constitutively activate mutant of MKK9. Gene expression analyses indicated that the activated MKK9-MPK3/MPK6 cascade upregulates the IGS biosynthetic genes. Moreover, the lack of MYB51, the transcription factor controlling biosynthetic genes responsible for synthesizing the IGS core structure, or CYP81F2, the enzyme catalyzing core structure modification to 4MI3G, significantly reduced mild osmotic stress- and MKK9 ^DD -induced 4MI3G. Thus, our study demonstrates that mild osmotic stress promotes 4MI3G biosynthesis and the accumulation in Arabidopsis through activation of the MKK9-MPK3/MPK6 cascade and provides an MAPK-mediated signaling pathway for the IGS response to abiotic stress in plants.

Proposed Method for Estimating Health-Promoting Glucosinolates and Hydrolysis Products in Broccoli (Brassica oleracea var. italica) Using Relative Transcript Abundance

Due to the importance of glucosinolates and their hydrolysis products in human nutrition and plant defense, optimizing the content of these compounds is a frequent breeding objective for Brassica crops. Toward this goal, we investigated the feasibility of using models built from relative transcript abundance data for the prediction of glucosinolate and hydrolysis product concentrations in broccoli. We report that predictive models explaining at least 50% of the variation for a number of glucosinolates and their hydrolysis products can be built for prediction within the same season, but prediction accuracy decreased when using models built from one season's data for prediction of an opposing season. This method of phytochemical profile prediction could potentially allow for lower phytochemical phenotyping costs and larger breeding populations. This, in turn, could improve selection efficiency for phase II induction potential, a type of chemopreventive bioactivity, by allowing for the quick and relatively cheap content estimation of phytochemicals known to influence the trait.

Arabidopsis myrosinases link the glucosinolate-myrosinase system and the cuticle

Both physical barriers and reactive phytochemicals represent two important components of a plant's defence system against environmental stress. However, these two defence systems have generally been studied independently. Here, we have taken an exclusive opportunity to investigate the connection between a chemical-based plant defence system, represented by the glucosinolate-myrosinase system, and a physical barrier, represented by the cuticle, using Arabidopsis myrosinase (thioglucosidase; TGG) mutants. The tgg1, single and tgg1 tgg2 double mutants showed morphological changes compared to wild-type plants visible as changes in pavement cells, stomatal cells and the ultrastructure of the cuticle. Extensive metabolite analyses of leaves from tgg mutants and wild-type Arabidopsis plants showed altered levels of cuticular fatty acids, fatty acid phytyl esters, glucosinolates, and indole compounds in tgg single and double mutants as compared to wild-type plants. These results point to a close and novel association between chemical defence systems and physical defence barriers.

Glucosinolate biosynthesis in Eruca sativa

Glucosinolates (GSLs) are a highly important group of secondary metabolites in the Caparalles order, both due to their significance in plant-biome interactions and to their chemoprotective properties. This study identified genes involved in all steps of aliphatic and indolic GSL biosynthesis in Eruca sativa, a cultivated plant closely related to Arabidopsis thaliana with agronomic and nutritional value. The impact of nitrogen (N) and sulfur (S) availability on GSL biosynthetic pathways at a transcriptional level, and on the final GSL content of plant leaf and root tissues, was investigated. N and S supply had a significant and interactive effect on the GSL content of leaves, in a structure-specific and tissue-dependent manner; the metabolites levels were significantly correlated with the relative expression of the genes involved in their biosynthesis. A more complex effect was observed in roots, where aliphatic and indolic GSLs and related biosynthetic genes responded differently to the various nutritional treatments suggesting that nitrogen and sulfur availability are important factors that control plant GSL content at a transcriptional level. The biological activity of extracts derived from these plants grown under the specific nutritional schemes was examined. N and S availability were found to significantly affect the cytotoxicity of E. sativa extracts on human cancer cells, supporting the notion that carefully designed nutritional schemes can promote the accumulation of chemoprotective substances in edible plants.

Sulfur deficiency-induced repressor proteins optimize glucosinolate biosynthesis in plants

Glucosinolates (GSLs) in the plant order of the Brassicales are sulfur-rich secondary metabolites that harbor antipathogenic and antiherbivory plant-protective functions and have medicinal properties, such as carcinopreventive and antibiotic activities. Plants repress GSL biosynthesis upon sulfur deficiency (-S); hence, field performance and medicinal quality are impaired by inadequate sulfate supply. The molecular mechanism that links -S to GSL biosynthesis has remained understudied. We report here the identification of the -S marker genes sulfur deficiency induced 1 (SDI1) and SDI2 acting as major repressors controlling GSL biosynthesis in Arabidopsis under -S condition. SDI1 and SDI2 expression negatively correlated with GSL biosynthesis in both transcript and metabolite levels. Principal components analysis of transcriptome data indicated that SDI1 regulates aliphatic GSL biosynthesis as part of -S response. SDI1 was localized to the nucleus and interacted with MYB28, a major transcription factor that promotes aliphatic GSL biosynthesis, in both yeast and plant cells. SDI1 inhibited the transcription of aliphatic GSL biosynthetic genes by maintaining the DNA binding composition in the form of an SDI1-MYB28 complex, leading to down-regulation of GSL biosynthesis and prioritization of sulfate usage for primary metabolites under sulfur-deprived conditions.

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.

Modification of oil and glucosinolate content in canola seeds with altered expression of Brassica napus LEAFY COTYLEDON1

Over the last few decades, research focusing on canola (Brassica napus L.) seed oil content and composition has expanded. Oil production and accumulation are influenced by genes participating in embryo and seed development. The Arabidopsis LEAFY COTYLEDON1 (LEC1) is a well characterized regulator of embryo development that also enhances the expression of genes involved in fatty acid (FA) synthesis. B. napus lines over-expressing or down-regulating BnLEC1 were successfully generated by Agrobacterium-mediated transformation. The constitutive expression of BnLEC1 in B. napus var. Polo, increased seed oil content by 7-16%, while the down-regulation of BnLEC1 in B. napus var. Topas reduced oil content by 9-12%. Experimental manipulation of BnLEC1 caused transcriptional changes in enzymes participating in sucrose metabolism, glycolysis, and FA biosynthesis, suggesting an enhanced carbon flux towards FA biosynthesis in tissues over-expressing BnLEC1. The increase in oil content induced by BnLEC1 was not accompanied by alterations in FA composition, oil nutritional value or glucosinolate (GLS) levels. Suppression of BnLEC1 reduced seed oil accumulation and elevated the level of GLS possibly through the transcriptional regulation of BnST5a (Sulphotransferase5a), the last GLS biosynthetic enzyme. Collectively, these findings demonstrate that experimental alterations of BnLEC1 expression can be used to influence oil production and quality in B. napus.

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.

Genome survey sequencing provides clues into glucosinolate biosynthesis and flowering pathway evolution in allotetrapolyploid Brassica juncea

BACKGROUND

Brassica juncea is an economically important vegetable crop in China, oil crop in India, condiment crop in Europe and selected for canola quality recently in Canada and Australia. B. juncea (2n = 36, AABB) is an allotetraploid derived from interspecific hybridization between B. rapa (2n = 20, AA) and B. nigra (2n = 16, BB), followed by spontaneous chromosome doubling.

RESULTS

Comparative genome analysis by genome survey sequence (GSS) of allopolyploid B. juncea with B. rapa was carried out based on high-throughput sequencing approaches. Over 28.35 Gb of GSS data were used for comparative analysis of B. juncea and B. rapa, producing 45.93% reads mapping to the B. rapa genome with a high ratio of single-end reads. Mapping data suggested more structure variation (SV) in the B. juncea genome than in B. rapa. We detected 2,921,310 single nucleotide polymorphisms (SNPs) with high heterozygosity and 113,368 SVs, including 1-3 bp Indels, between B. juncea and B. rapa. Non-synonymous polymorphisms in glucosinolate biosynthesis genes may account for differences in glucosinolate biosynthesis and glucosinolate components between B. juncea and B. rapa. Furthermore, we identified distinctive vernalization-dependent and photoperiod-dependent flowering pathways coexisting in allopolyploid B. juncea, suggesting contribution of these pathways to adaptation for survival during polyploidization.

CONCLUSIONS

Taken together, we proposed that polyploidization has allowed for accelerated evolution of the glucosinolate biosynthesis and flowering pathways in B. juncea that likely permit the phenotypic variation observed in the crop.

Novel bioresources for studies of Brassica oleracea: identification of a kale MYB transcription factor responsible for glucosinolate production

Plants belonging to the Brassicaceae family exhibit species-specific profiles of glucosinolates (GSLs), a class of defence compounds against pathogens and insects. GSLs also exhibit various human health-promoting properties. Among them, glucoraphanin (aliphatic 4-methylsulphinylbutyl GSL) has attracted the most attention because it hydrolyses to form a potent anticancer compound. Increased interest in developing commercial varieties of Brassicaceae crops with desirable GSL profiles has led to attempts to identify genes that are potentially valuable for controlling GSL biosynthesis. However, little attention has been focused on genes of kale (Brassica oleracea var. acephala). In this study, we established full-length kale cDNA libraries containing 59 904 clones, which were used to generate an expressed sequence tag (EST) data set with 119 204 entries. The EST data set clarified genes related to the GSL biosynthesis pathway in kale. We specifically focused on BoMYB29, a homolog of Arabidopsis MYB29/PMG2/HAG3, not only to characterize its function but also to demonstrate its usability as a biological resource. BoMYB29 overexpression in wild-type Arabidopsis enhanced the expression of aliphatic GSL biosynthetic genes and the accumulation of aliphatic GSLs. When expressed in the myb28myb29 mutant, which exhibited no detectable aliphatic GSLs, BoMYB29 restored the expression of biosynthetic genes and aliphatic GSL accumulation. Interestingly, the ratio of methylsulphinyl GSL content, including glucoraphanin, to that of methylthio GSLs was greatly increased, indicating the suitability of BoMYB29 as a regulator for increasing methylsulphinyl GSL content. Our results indicate that these biological resources can facilitate further identification of genes useful for modifications of GSL profiles and accumulation in kale.

AtMYB44 regulates resistance to the green peach aphid and diamondback moth by activating EIN2-affected defences in Arabidopsis

Recently we showed that the transcription activator AtMYB44 regulates expression of EIN2, a gene essential for ethylene signalling and insect resistance, in Arabidopsis thaliana (Arabidopsis). To link the transactivation with insect resistance, we investigated the wild-type and atmyb44 mutant plants, genetically Complemented atmyb44 (Catmyb44) and AtMYB44-Overexpression Transgenic Arabidopsis (MYB44OTA). We found that AtMYB44 played a critical role in Arabidopsis resistance to the phloem-feeding generalist green peach aphid (Myzus persicae Sulzer) and leaf-chewing specialist caterpillar diamondback moth (Plutella xylostella L.). AtMYB44 was required not only for the development of constitutive resistance but also for the induction of resistance by both herbivorous insects. Levels of constitutive and herbivore-induced resistance were consistent with corresponding amounts of the AtMYB44 protein constitutively produced in MYB44OTA and induced by herbivory in Catmyb44. In both cases, AtMYB44 promoted EIN2 expression to a greater extent in MYB44OTA than in Catmyb44. However, AtMYB44-promoted EIN2 expression was arrested with reduced resistance levels in the EIN2-deficient Arabidopsis mutant ein2-1 and the MYB44OTA ein2-1 hybrid. In the different plant genotypes, only MYB44OTA constitutively displayed phloem-based defences, which are specific to phloem-feeding insects, and robust expression of genes involved in the biosynthesis of glucosinolates, which are the secondary plant metabolites known as deterrents to generalist herbivores. Phloem-based defences and glucosinolate-related gene expression were not detected in ein2-1 and MYB44OTA ein2-1. These results establish a genetic connection between the regulatory role of AtMYB44 in EIN2 expression and the development of Arabidopsis resistance to insects.

Targeted silencing of BjMYB28 transcription factor gene directs development of low glucosinolate lines in oilseed Brassica juncea

Brassica juncea (Indian mustard), a globally important oilseed crop, contains relatively high amount of seed glucosinolates ranging from 80 to 120 μmol/g dry weight (DW). One of the major breeding objectives in oilseed Brassicas is to improve the seed-meal quality through the development of low-seed-glucosinolate lines (<30 μmol/g DW), as high amounts of certain seed glucosinolates are known to be anti-nutritional and reduce the meal palatability. Here, we report the development of transgenic B. juncea lines having seed glucosinolates as low as 11.26 μmol/g DW, through RNAi-based targeted suppression of BjMYB28, a R2R3-MYB transcription factor family gene involved in aliphatic glucosinolate biosynthesis. Targeted silencing of BjMYB28 homologs provided significant reduction in the anti-nutritional aliphatic glucosinolates fractions, without altering the desirable nonaliphatic glucosinolate pool, both in leaves and seeds of transgenic plants. Molecular characterization of single-copy, low glucosinolate homozygous lines confirmed significant down-regulation of BjMYB28 homologs vis-à-vis enhanced accumulation of BjMYB28-specific siRNA pool. Consequently, these low glucosinolate lines also showed significant suppression of genes involved in aliphatic glucosinolate biosynthesis. The low glucosinolate trait was stable in subsequent generations of the transgenic lines with no visible off-target effects on plant growth and development. Various seed quality parameters including fatty acid composition, oil content, protein content and seed weight of the low glucosinolate lines also remained unaltered, when tested under containment conditions in the field. Our results indicate that targeted silencing of a key glucosinolate transcriptional regulator MYB28 has huge potential for reducing the glucosinolates content and improving the seed-meal quality of oilseed Brassica crops.

Evaluation of glucosinolate variation in a collection of turnip (Brassica rapa) germplasm by the analysis of intact and desulfo glucosinolates

Glucosinolates (GLS) are secondary metabolites occurring in cruciferous species. These compounds are important for plant defense, human health, and the characteristic flavor of Brassica vegetables. In this study, the GLS in tubers from a collection of 48 turnip ( Brassica rapa ) accessions from different geographic origin were analyzed. Two different methods were used: desulfo GLS were analyzed by high-performance liquid chromatography with a photodiode array detector, and intact GLS were analyzed by accurate mass liquid chromatography-mass spectrometry. For most GLS, desulfo and intact signals correlated well, and the analytical reproducibility for individual GLS was similar for both methods. A total of 11 different GLS was monitored in the turnip tubers, through both intact and desulfo GLS analysis methods. Four clusters of accessions could be clearly distinguished based on GLS composition of the turnip tuber. Clustering based on tuber GLS differed markedly from a previously published clustering based on leaf GLS.

Glucosinolate variation in leaves of Brassica rapa crops

Total and individual glucosinolate (GSL) content of leaves of vegetable turnip rape (Brassica rapa L. var. rapa) was determined in a set of 45 varieties consisting in early, medium and late types grown at two locations in northwestern Spain. The objectives were to determine the diversity among varieties in GSL content and to relate that variation with earliness and plant habit. Eight GSL were identified, being two aliphatic GSL, gluconapin (84.4 % of the total GSL) and glucobrassicanapin (7.2 % of the total GSL) the most abundant. Indolic and aromatic GSL content were low but also showed significant differences among varieties. Differences in total and individual GSL content were found among varieties, plant habit groups, and earliness groups. Total GSL content ranged from 19 to 37.3 μmol g(-1) dw in early and extra-late groups, respectively, and from 19.5 to 36.3 μmol g(-1) dw for turnips and turnip greens groups, respectively. These differences were consistent to values found for gluconapin content where the turnip group had the highest values (31.8 μmol g(-1) dw) and the turnip top group had the lowest (15.7 μmol g(-1) dw). Two varieties, MBG-BRS0429 and MBG-BRS0550 (from turnip greens and extra-late groups) and MBG-BRS0438 (from turnips and late groups), stood out as they had the highest total GSL content and could be used as a good source of these beneficial bioactive compounds. Elucidation of genetic diversity among crops can provide useful information to assist plant breeders to design improved breeding strategies in order to obtain varieties rich on GSL.

Isolation and expression of glucosinolate synthesis genes CYP83A1 and CYP83B1 in Pak Choi (Brassica rapa L. ssp. chinensis var. communis (N. Tsen & S.H. Lee) Hanelt)

CYP83A1 and CYP83B1 are two key synthesis genes in the glucosinolate biosynthesis pathway. CYP83A1 mainly metabolizes the aliphatic oximes to form aliphatic glucosinolate and CYP83B1 mostly catalyzes aromatic oximes to synthesis corresponding substrates for aromatic and indolic glucosinolates. In this study, two CYP83A1 genes named BcCYP83A1-1 (JQ289997), BcCYP83A1-2 (JQ289996) respectively and one CYP83B1 (BcCYP83B1, HM347235) gene were cloned from the leaves of pak choi (Brassica rapa L. ssp. chinensis var. communis (N. Tsen & S.H. Lee) Hanelt) “Hangzhou You Dong Er” cultivar. Their ORFs were 1506, 1509 and 1500 bp in length, encoding 501, 502 and 499 amino acids, respectively. The predicted amino acid sequences of CYP83A1-1, CYP83A1-2 and CYP83B1 shared high sequence identity of 87.65, 86.48 and 95.59% to the corresponding ones in Arabidopsis, and 98.80, 98.61 and 98.80% to the corresponding ones in Brassica pekinensis (Chinese cabbage), respectively. Quantitative real-time PCR analysis indicated that both CYP83A1 and CYP83B1 expressed in roots, leaves and petioles of pak choi, while the transcript abundances of CYP83A1 were higher in leaves than in petioles and roots, whereas CYP83B1 showed higher abundances in roots. The expression levels of glucosinolate biosynthetic genes were consistent with the glucosinolate profile accumulation in shoots of seven cultivars and three organs. The isolation and characterization of the glucosinolate synthesis genes in pak choi would promote the way for further development of agronomic traits via genetic engineering.

Expression pattern of the glucosinolate side chain biosynthetic genes MAM1 and MAM3 of Arabidopsis thaliana in different organs and developmental stages

Aliphatic glucosinolates, secondary metabolites known to be involved in plant defence, make up the majority of the glucosinolate content of Arabidopsis thaliana, and their structural diversity arises in part from chain elongations of methionine before the formation of the glucosinolate core structure. The key enzymatic step in determining the length of the chain is the condensation of acetyl-coenzyme A with a series of ω-methylthio-2-oxoalkanoic acids, catalyzed by methylthioalkylmalate (MAM) synthases. The existence of two MAM synthases has been previously reported in A. thaliana, ecotype Columbia-0. MAM1 catalyses the condensation step of the first three elongation cycles while MAM3 catalyzes the condensation step of all six elongation cycles. We studied the expression patterns of MAM1 and MAM3 genes in different organs and developmental stages using promoter-GUS fusion lines and qRT-PCR. The promoter-GUS lines revealed MAM1 and MAM3 expression in varying degrees in all organs, but this was generally restricted to the phloem, except in wounded tissue where expression was general. No difference was found between the two genes. The qRT-PCR measurements showed that expression was generally highest in seedlings and vegetative parts at the reproductive phase, but low in flowers and fruits. Since high amounts of glucosinolates accumulate in flowers and fruits, these data indicate possible transport from vegetative to reproductive organs. The expression of MAM1 was different than that of MAM3 with MAM3 having relative more expression in seedlings and roots than MAM1.

Engineering of glucosinolate biosynthesis: candidate gene identification and validation

The diverse biological roles of glucosinolates as plant defense metabolites and anticancer compounds have spurred a strong interest in their biosynthetic pathways. Since the completion of the Arabidopsis genome, functional genomics approaches have enabled significant progress on the elucidation of glucosinolate biosynthesis, although in planta validation of candidate gene function often is hampered by time-consuming generation of knockout and overexpression lines in Arabidopsis. To better exploit the increasing amount of data available from genomic sequencing, microarray database and RNAseq, time-efficient methods for identification and validation of candidate genes are needed. This chapter covers the methodology we are using for gene discovery in glucosinolate engineering, namely, guilt-by-association-based in silico methods and fast proof-of-function screens by transient expression in Nicotiana benthamiana. Moreover, the lessons learned in the rapid, transient tobacco system are readily translated to our robust, versatile yeast expression platform, where additional genes critical for large-scale microbial production of glucosinolates can be identified. We anticipate that the methodology presented here will be beneficial to elucidate and engineer other plant biosynthetic pathways.

Varied response of Spodoptera littoralis against Arabidopsis thaliana with metabolically engineered glucosinolate profiles

Upon herbivory glucosinolates are known to be degraded into a cascade of secondary products that can be detrimental for certain herbivores. We performed herbivory bioassays using first and second instar generalist Lepidoptera larvae Spodoptera littoralis on Arabidopsis thaliana engineered to overexpress novel glucosinolates. A differential response in larval feeding patterns was observed on the plants engineered with novel glucosinolates. Larvae fed on plants overexpressing 4-hydroxybenzyl glucosinolate and isopropyl glucosinolate showed little response. Larvae fed on 35S:CYP79A2 plants engineered to overexpress benzyl glucosinolates, however, showed reduced larval and pupal weights. Upon herbivory a high expression of JA signalling gene LOX2 was observed on the 35S:CYP79A2 plants compared to the PR1a and VSP2 expression. To confirm the role of benzyl isothiocyanate (BITC), a degradation product of benzyl glucosinolate overexpressing plants, in the retarded larval growth we used Virus Induced Gene Silencing (VIGS) approach to silence LOX2 expression in the 35S:CYP79A2 plants. S. littoralis larvae fed on LOX2 silenced 35S:CYP79A2 plants exhibited a retarded larval growth thus indicating that BITC played a pivotal role in anti-herbivory and not only the JA signalling pathway.

Genomic analysis of QTLs and genes altering natural variation in stochastic noise

Quantitative genetic analysis has long been used to study how natural variation of genotype can influence an organism's phenotype. While most studies have focused on genetic determinants of phenotypic average, it is rapidly becoming understood that stochastic noise is genetically determined. However, it is not known how many traits display genetic control of stochastic noise nor how broadly these stochastic loci are distributed within the genome. Understanding these questions is critical to our understanding of quantitative traits and how they relate to the underlying causal loci, especially since stochastic noise may be directly influenced by underlying changes in the wiring of regulatory networks. We identified QTLs controlling natural variation in stochastic noise of glucosinolates, plant defense metabolites, as well as QTLs for stochastic noise of related transcripts. These loci included stochastic noise QTLs unique for either transcript or metabolite variation. Validation of these loci showed that genetic polymorphism within the regulatory network alters stochastic noise independent of effects on corresponding average levels. We examined this phenomenon more globally, using transcriptomic datasets, and found that the Arabidopsis transcriptome exhibits significant, heritable differences in stochastic noise. Further analysis allowed us to identify QTLs that control genomic stochastic noise. Some genomic QTL were in common with those altering average transcript abundance, while others were unique to stochastic noise. Using a single isogenic population, we confirmed that natural variation at ELF3 alters stochastic noise in the circadian clock and metabolism. Since polymorphisms controlling stochastic noise in genomic phenotypes exist within wild germplasm for naturally selected phenotypes, this suggests that analysis of Arabidopsis evolution should account for genetic control of stochastic variance and average phenotypes. It remains to be determined if natural genetic variation controlling stochasticity is equally distributed across the genomes of other multi-cellular eukaryotes.

Identification of enzymatic and regulatory genes of plant metabolism through QTL analysis in Arabidopsis

The biochemical diversity in the plant kingdom is estimated to well exceed 100,000 distinct compounds (Weckwerth, 2003) and 4000 to 20,000 metabolites per species seem likely (Fernie et al., 2004). In recent years extensive progress has been made towards the identification of enzymes and regulatory genes working in a complex network to generate this large arsenal of metabolites. Genetic loci influencing quantitative traits, e.g. metabolites or biomass, may be mapped to associated molecular markers, a method called quantitative trait locus mapping (QTL mapping), which may facilitate the identification of novel genes in biochemical pathways. Arabidopsis thaliana, as a model organism for seed plants, is a suitable target for metabolic QTL (mQTL) studies due to the availability of highly developed molecular and genetic tools, and the extensive knowledge accumulated on the metabolite profile. While intensely studied, in particular since the availability of its complete sequence, the genome of Arabidopsis still comprises a large proportion of genes with only tentative function based on sequence homology. From a total number of 33,518 genes currently listed (TAIR 9, http://www.arabidopsis.org), only about 25% have direct experimental evidence for their molecular function and biological process, while for more than 30% no biological data are available. Modern metabolomics approaches together with continually extended genomic resources will facilitate the task of assigning functions to those genes. In our previous study we reported on the identification of mQTL (Lisec et al., 2008). In this paper, we summarize the current status of mQTL analyses and causal gene identification in Arabidopsis and present evidence that a candidate gene located within the confidence interval of a fumarate mQTL (AT5G50950) encoding a putative fumarase is likely to be the causal gene of this QTL. The total number of genes molecularly identified based on mQTL studies is still limited, but the advent of multi-parallel analysis techniques for measurement of gene expression, as well as protein and metabolite abundances and for rapid gene identification will assist in the important task of assigning enzymes and regulatory genes to the growing network of known metabolic reactions.

Combining genome-wide association mapping and transcriptional networks to identify novel genes controlling glucosinolates in Arabidopsis thaliana

BACKGROUND

Genome-wide association (GWA) is gaining popularity as a means to study the architecture of complex quantitative traits, partially due to the improvement of high-throughput low-cost genotyping and phenotyping technologies. Glucosinolate (GSL) secondary metabolites within Arabidopsis spp. can serve as a model system to understand the genomic architecture of adaptive quantitative traits. GSL are key anti-herbivory defenses that impart adaptive advantages within field trials. While little is known about how variation in the external or internal environment of an organism may influence the efficiency of GWA, GSL variation is known to be highly dependent upon the external stresses and developmental processes of the plant lending it to be an excellent model for studying conditional GWA.

METHODOLOGY/PRINCIPAL FINDINGS

To understand how development and environment can influence GWA, we conducted a study using 96 Arabidopsis thaliana accessions, >40 GSL phenotypes across three conditions (one developmental comparison and one environmental comparison) and ∼230,000 SNPs. Developmental stage had dramatic effects on the outcome of GWA, with each stage identifying different loci associated with GSL traits. Further, while the molecular bases of numerous quantitative trait loci (QTL) controlling GSL traits have been identified, there is currently no estimate of how many additional genes may control natural variation in these traits. We developed a novel co-expression network approach to prioritize the thousands of GWA candidates and successfully validated a large number of these genes as influencing GSL accumulation within A. thaliana using single gene isogenic lines.

CONCLUSIONS/SIGNIFICANCE

Together, these results suggest that complex traits imparting environmentally contingent adaptive advantages are likely influenced by up to thousands of loci that are sensitive to fluctuations in the environment or developmental state of the organism. Additionally, while GWA is highly conditional upon genetics, the use of additional genomic information can rapidly identify causal loci en masse.

Role of camalexin, indole glucosinolates, and side chain modification of glucosinolate-derived isothiocyanates in defense of Arabidopsis against Sclerotinia sclerotiorum

Plant secondary metabolites are known to facilitate interactions with a variety of beneficial and detrimental organisms, yet the contribution of specific metabolites to interactions with fungal pathogens is poorly understood. Here we show that, with respect to aliphatic glucosinolate-derived isothiocyanates, toxicity against the pathogenic ascomycete Sclerotinia sclerotiorum depends on side chain structure. Genes associated with the formation of the secondary metabolites camalexin and glucosinolate were induced in Arabidopsis thaliana leaves challenged with the necrotrophic pathogen S. sclerotiorum. Unlike S. sclerotiorum, the closely related ascomycete Botrytis cinerea was not identified to induce genes associated with aliphatic glucosinolate biosynthesis in pathogen-challenged leaves. Mutant plant lines deficient in camalexin, indole, or aliphatic glucosinolate biosynthesis were hypersusceptible to S. sclerotiorum, among them the myb28 mutant, which has a regulatory defect resulting in decreased production of long-chained aliphatic glucosinolates. The antimicrobial activity of aliphatic glucosinolate-derived isothiocyanates was dependent on side chain elongation and modification, with 8-methylsulfinyloctyl isothiocyanate being most toxic to S. sclerotiorum. This information is important for microbial associations with cruciferous host plants and for metabolic engineering of pathogen defenses in cruciferous plants that produce short-chained aliphatic glucosinolates.

Evidence for a nitrate-independent function of the nitrate sensor NRT1.1 in Arabidopsis thaliana

NRT1.1 is a putative nitrate sensor and is involved in many nitrate-dependent responses. On the other hand, a nitrate-independent function of NRT1.1 has been implied, but the clear-cut evidence is unknown. We found that NRT1.1 mutants showed enhanced tolerance to concentrated ammonium as sole N source in Arabidopsis thaliana. This unique phenotype was not observed in mutants of NLP7, which has been suggested to play a role in the nitrate-dependent signaling pathway. Our real-time PCR analysis, and evidence from a literature survey revealed that several genes relevant to the aliphatic glucosinolate-biosynthetic pathway were regulated via a nitrate-independent signal from NRT1.1. When taken together, the present study strongly suggests the existence of a nitrate-independent function of NRT1.1.

Barbarea vulgaris linkage map and quantitative trait loci for saponins, glucosinolates, hairiness and resistance to the herbivore Phyllotreta nemorum

Combined genomics and metabolomics approaches were used to unravel molecular mechanisms behind interactions between winter cress (Barbarea vulgaris) and flea beetle (Phyllotreta nemorum). B. vulgaris comprises two morphologically, biochemically and cytologically deviating types, which differ in flea beetle resistance, saponin and glucosinolate profiles, as well as leaf pubescence. An F2 population generated from a cross between the two B. vulgaris types was used to construct a B. vulgaris genetic map based on 100 AFLP and 31 microsatellite markers. The map was divided into eight linkage groups. QTL (quantitative trait loci) analysis revealed a total of 15 QTL affecting eight traits, including nine QTL for four saponins, two QTL for two glucosinolates, two QTL for hairiness, and two QTL for flea beetle resistance. The two QTL for resistance towards flea beetles in B. vulgaris co-localized with QTL for the four saponins associated with resistance. Furthermore, global QTL analysis of B. vulgaris metabolites identified QTL for a number of flavonoid glycosides and additional saponins from both resistant and susceptible types. The transcriptome of the resistant B. vulgaris type was sequenced by pyrosequencing, and sequences containing microsatellites were identified. Microsatellite types in B. vulgaris were similar to Arabidopsis thaliana but different from Oryza sativa. Comparative analysis between B. vulgaris and A. thaliana revealed a remarkable degree of synteny between a large part of linkage groups 1 and 4 of B. vulgaris harboring the two QTL for flea beetle resistance and Arabidopsis chromosomes 3 and 1. Gene candidates that may underlie QTL for resistance and saponin biosynthesis are discussed.

Glucosinolate biochemical diversity and innovation in the Brassicales

Glucosinolates were analysed from herbarium specimens and living tissues from representative of all families of the Brassicales, following the phylogenetic schemes of Rodman et al. (1998) and Hall et al. (2002, 2004), including specimens of Akania, Setchellanthus, Emblingia, Stixis, Forchhammeria and members of the Capparaceae for which glucosinolate content had not previously been reported. The results are reviewed along with additional published data on glucosinolate content of members of the Brassicales. In addition to providing an overview of the evolution of glucosinolate biochemical diversity within the core Brassicales, there were three main findings. Firstly, the glucosinolate content of some 'orphan' taxa of the Brassicales, such as Setchellanthus and Emblingia were consistent with recent phylogentic analyses based upon DNA sequence comparisons, while further analyses of Tirania and Stixis is required. Secondly, methyl glucosinolate is found within the Capparaceae and Cleomaceae, but also, unexpectedly, within Forchhammeria, with implications for the biochemical and evolutionary origin of methyl glucosinolate and the phylogenetic relationships of Forchhammeria. Thirdly, whereas Old World Capparaceae contain methyl glucosinolate, New World Capparaceae, including New World Capparis, either contain methyl glucosinolates or glucosinolates of complex and unresolved structures, indicative of continued innovation in glucosinolate biosynthesis. These taxa may be productive sources of glucosinolate biosynthetic genes and alleles that are not found in the model plant Arabidopsis thaliana.

Regulatory networks of glucosinolates shape Arabidopsis thaliana fitness

Systems biology approaches address higher levels of complex, but dynamic metabolic regulatory networks utilizing single accessions of a species. This contrasts with the likelihood that plants utilize genetic diversity of both individual genes and regulatory networks as a solution to surviving in a complex environment. This would require systems biology to begin a more inclusive search for 'all' networks within a species. In this review, we will highlight how natural genetic diversity within particularly aliphatic glucosinolates in Arabidopsis thaliana and related species has resulted in highly complex, dynamic regulatory networks enabling the plant to adapt to a highly changing environment. We will discuss how this diversity is essential for the fitness performance of A. thaliana.