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

Developing multifunctional crops by engineering Brassicaceae glucosinolate pathways

Glucosinolates (GSLs), found mainly in species of the Brassicaceae family, are one of the most well-studied classes of secondary metabolites. Produced by the action of myrosinase on GSLs, GSL-derived hydrolysis products (GHPs) primarily defend against biotic stress in planta. They also significantly affect the quality of crop products, with a subset of GHPs contributing unique food flavors and multiple therapeutic benefits or causing disagreeable food odors and health risks. Here, we explore the potential of these bioactive functions, which could be exploited for future sustainable agriculture. We first summarize our accumulated understanding of GSL diversity and distribution across representative Brassicaceae species. We then systematically discuss and evaluate the potential of exploited and unutilized genes involved in GSL biosynthesis, transport, and hydrolysis as candidate GSL engineering targets. Benefiting from available information on GSL and GHP functions, we explore options for multifunctional Brassicaceae crop ideotypes to meet future demand for food diversification and sustainable crop production. An integrated roadmap is subsequently proposed to guide ideotype development, in which maximization of beneficial effects and minimization of detrimental effects of GHPs could be combined and associated with various end uses. Based on several use-case examples, we discuss advantages and limitations of available biotechnological approaches that may contribute to effective deployment and could provide novel insights for optimization of future GSL engineering.

Glucosinolate Biosynthetic Genes of Cabbage: Genome-Wide Identification, Evolution, and Expression Analysis

Cabbage (Brassica oleracea var. capitata) is a vegetable rich in glucosinolates (GSLs) that have proven health benefits. To gain insights into the synthesis of GSLs in cabbage, we systematically analyzed GSLs biosynthetic genes (GBGs) in the entire cabbage genome. In total, 193 cabbage GBGs were identified, which were homologous to 106 GBGs in Arabidopsis thaliana. Most GBGs in cabbage have undergone negative selection. Many homologous GBGs in cabbage and Chinese cabbage differed in expression patterns indicating the unique functions of these homologous GBGs. Spraying five exogenous hormones significantly altered expression levels of GBGs in cabbage. For example, MeJA significantly upregulated side chain extension genes BoIPMILSU1-1 and BoBCAT-3-1, and the expression of core structure construction genes BoCYP83A1 and BoST5C-1, while ETH significantly repressed the expression of side chain extension genes such as BoIPMILSU1-1, BoCYP79B2-1, and BoMAMI-1, and some transcription factors, namely BoMYB28-1, BoMYB34-1, BoMYB76-1, BoCYP79B2-1, and BoMAMI-1. Phylogenetically, the CYP83 family and CYP79B and CYP79F subfamilies may only be involved in GSL synthesis in cruciferous plants. Our unprecedented identification and analysis of GBGs in cabbage at the genome-wide level lays a foundation for the regulation of GSLs synthesis through gene editing and overexpression.

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

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

Seed coat colour of Indian mustard [Brassica juncea (L.) Czern. and Coss.] is associated with Bju.TT8 homologs identifiable by targeted functional markers

Seed coat colour is an important trait in Indian mustard. Breeding for seed coat colour needs precise knowledge of mode of inheritance and markers linked to it. The present study was focussed on genetics and development of functional markers for seed coat colour. F₁s (direct and reciprocal) and F₂ populations were developed by crossing two contrasting parents for seed coat colour (DRMRIJ-31, brown seeded and RLC-3, yellow seeded). Phenotypic results have shown that the seed coat colour trait was under the influence of maternal effect and controlled by digenic-duplicate gene action. Further, Bju.TT8 homologs of both parents (DRMRIJ-31 and RLC-3) were cloned and sequenced. Sequencing results of Bju.TT8 homologs revealed that in RLC-3, gene Bju.ATT8 had an insertion of 1279bp in the 7^th exon; whereas, gene Bju.BTT8 had an SNP (C→T) in the 7^th exon. These two mutations were found to be associated with yellow seed coat colour. Using sequence information, functional markers were developed for both Bju.TT8 homologs, validated on F₂ population and were found highly reliable with no recombination between the markers and the phenotype. Further, these markers were subjected to a germplasm assembly of Indian mustard, and their allelic combination for the seed coat colour genes has been elucidated. The comparative genomics of TT8 genes revealed high degree of similarity between and across the Brassica species, and the respective diploid progenitors in tetraploid Brassica species are the possible donors of TT8 homologs. This study will help in the marker-assisted breeding for seed coat colour, and aid in understanding seed coat colour genetics more precisely.

Genome-wide association mapping within a local Arabidopsis thaliana population more fully reveals the genetic architecture for defensive metabolite diversity

A paradoxical finding from genome-wide association studies (GWAS) in plants is that variation in metabolite profiles typically maps to a small number of loci, despite the complexity of underlying biosynthetic pathways. This discrepancy may partially arise from limitations presented by geographically diverse mapping panels. Properties of metabolic pathways that impede GWAS by diluting the additive effect of a causal variant, such as allelic and genetic heterogeneity and epistasis, would be expected to increase in severity with the geographical range of the mapping panel. We hypothesized that a population from a single locality would reveal an expanded set of associated loci. We tested this in a French Arabidopsis thaliana population (less than 1 km transect) by profiling and conducting GWAS for glucosinolates, a suite of defensive metabolites that have been studied in depth through functional and genetic mapping approaches. For two distinct classes of glucosinolates, we discovered more associations at biosynthetic loci than the previous GWAS with continental-scale mapping panels. Candidate genes underlying novel associations were supported by concordance between their observed effects in the TOU-A population and previous functional genetic and biochemical characterization. Local populations complement geographically diverse mapping panels to reveal a more complete genetic architecture for metabolic traits. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

A Comparative Transcriptome and Metabolome Combined Analysis Reveals the Key Genes and Their Regulatory Model Responsible for Glucoraphasatin Accumulation in Radish Fleshy Taproots

Radish (Raphanus sativus L.) is rich in specific glucosinolates (GSLs), which benefit human health and special flavor formation. Although the basic GSLs metabolic pathway in Brassicaceae plants is clear, the regulating mechanism for specific glucosinolates content in radish fleshy taproots is not well understood. In this study, we discovered that there was a significant difference in the GSLs profiles and the content of various GSLs components. Glucoraphasatin (GRH) is the most predominant GSL in radish taproots of different genotypes as assessed by HPLC analysis. Further, we compared the taproot transcriptomes of three radish genotypes with high and low GSLs content by employing RNA-seq. Totally, we identified forty-one differentially expressed genes related to GSLs metabolism. Among them, thirteen genes (RsBCAT4, RsIPMDH1, RsMAM1a, RsMAM1b, RsCYP79F1, RsGSTF9, RsGGP1, RsSUR1, RsUGT74C1, RsST5b, RsAPK1, RsGSL-OH, and RsMYB28) were significantly higher co-expressed in the high content genotypes than in low content genotype. Notably, correlation analysis indicated that the expression level of RsMYB28, as an R2R3 transcription factor directly regulating aliphatic glucosinolate biosynthesis, was positively correlated with the GRH content. Co-expression network showed that RsMYB28 probably positively regulated the expression of the above genes, particularly RsSUR1, and consequently the synthesis of GRH. Moreover, the molecular mechanism of the accumulation of this 4-carbon (4C) GSL in radish taproots was explored. This study provides new perspectives on the GSLs accumulation mechanism and genetic improvements in radish taproots.

Genome-Wide Association Reveals Trait Loci for Seed Glucosinolate Accumulation in Indian Mustard (Brassica juncea L.)

Glucosinolates (GSLs) are sulphur- and nitrogen-containing secondary metabolites implicated in the fitness of Brassicaceae and appreciated for their pungency and health-conferring properties. In Indian mustard (Brassica juncea L.), GSL content and composition are seed-quality-determining traits affecting its economic value. Depending on the end use, i.e., condiment or oil, different GSL levels constitute breeding targets. The genetic control of GSL accumulation in Indian mustard, however, is poorly understood, and current knowledge of GSL biosynthesis and regulation is largely based on Arabidopsis thaliana. A genome-wide association study was carried out to dissect the genetic architecture of total GSL content and the content of two major GSLs, sinigrin and gluconapin, in a diverse panel of 158 Indian mustard lines, which broadly grouped into a South Asia cluster and outside-South-Asia cluster. Using 14,125 single-nucleotide polymorphisms (SNPs) as genotyping input, seven distinct significant associations were discovered for total GSL content, eight associations for sinigrin content and 19 for gluconapin. Close homologues of known GSL structural and regulatory genes were identified as candidate genes in proximity to peak SNPs. Our results provide a comprehensive map of the genetic control of GLS biosynthesis in Indian mustard, including priority targets for further investigation and molecular marker development.

Current achievements and future prospects of genetic engineering in Indian mustard (Brassica juncea L. Czern & Coss.)

Transgenic technology in Indian mustard has expedited crop improvement programs. Further, there is a need to optimize gene editing protocols and find out the suitable target genes to harvest the benefits of gene editing technology in this important edible oilseed crop. Brassica juncea is an economically and industrially important oilseed crop being grown mainly in India and in some parts of Canada, Russia, China and Australia. Besides being consumed as edible oil, it also has numerous applications in food and paint industry. However, its overall production and productivity are being hampered by a number of biotic and abiotic stress factors. Further, its oil and seedmeal quality needs to be improved for increasing food as well as feed value. However, the lack of resistant crossable germplasm or varieties necessitated the use of genetic engineering interventions in Indian mustard crop improvement. A number of genes conferring resistance to biotic stresses including lectins for aphids' control, chitinase, glucanase and osmotin for disease control and for abiotic stresses, CODA, LEA and ion antiporter genes have been transferred to Indian mustard. Both antisense and RNAi technologies have been employed for improving oil and seedmeal quality. Efforts have been made to improve the phytoremediation potential of this crop through genetic engineering approach. The deployment of barnase/barstar gene system for developing male sterile and restorer lines has really expedited hybrid development programs in Indian mustard. Further, there is a need to optimize gene editing protocols and to find out suitable target genes for gene editing in this crop. In this review paper, authors have attempted to review various genetic transformation efforts carried out in Indian mustard for its improvement to combat biotic and abiotic stress challenges, quality improvement and hybrid development.

Comparative analysis of glucosinolates and metabolite profiling of green and red mustard (brassica juncea) hairy roots

Here, accumulation of glucosinolates and expression of glucosinolates biosynthesis genes in green and red mustard hairy roots were identified and quantified by HPLC and qRT-PCR analyses. The total glucosinolates content of green mustard hairy root (10.09 µg/g dry weight) was 3.88 times higher than that of red mustard hairy root. Indolic glucosinolates (glucobrassicin, 4-methoxyglucobrassicin, and neoglucobrassicin) in green mustard were found at 30.92, 6.95, and 5.29 times higher than in red mustard hairy root, respectively. Conversely, levels of glucotropaeolin (aromatic glucosinolate) was significantly higher in red mustard than in green mustard. Accumulation of glucoraphasatin, an aliphatic glucosinolate, was only observed only in red mustard hairy roots. Quantitative real-time PCR analysis showed that the expression level of genes related to aliphatic and aromatic glucosinolate biosynthesis were higher in red mustard, exception BjCYP83B. The expression of BjCYP79B2, which encodes a key enzyme involved in the indolic glucosinolate biosynthetic pathway, was higher in green mustard than in red mustard. Additionally, to further distinguish between green mustard and red mustard hairy roots, hydrophilic and lipophilic compounds were identified by gas chromatography-mass spectrometry and subjected to principal component analysis. The results indicated that core primary metabolites and glucosinolate levels were higher in the hairy roots of green mustard than in those of red mustard.

Exogenous Methyl Jasmonate and Salicylic Acid Induce Subspecies-Specific Patterns of Glucosinolate Accumulation and Gene Expression in Brassica oleracea L

Glucosinolates have anti-carcinogenic properties. In the recent decades, the genetics of glucosinolate biosynthesis has been widely studied, however, the expression of specific genes involved in glucosinolate biosynthesis under exogenous phytohormone treatment has not been explored at the subspecies level in Brassica oleracea. Such data are vital for strategies aimed at selective exploitation of glucosinolate profiles. This study quantified the expression of 38 glucosinolate biosynthesis-related genes in three B. oleracea subspecies, namely cabbage, broccoli and kale, and catalogued associations between gene expression and increased contents of individual glucosinolates under methyl jasmonate (MeJA) and salicylic acid (SA) treatments. Glucosinolate accumulation and gene expression in response to phytohormone elicitation was subspecies specific. For instance, cabbage leaves showed enhanced accumulation of the aliphatic glucoiberin, progoitrin, sinigrin and indolic neoglucobrassicin under both MeJA and SA treatment. MeJA treatment induced strikingly higher accumulation of glucobrassicin (GBS) in cabbage and kale and of neoglucobrassicin (NGBS) in broccoli compared to controls. Notably higher expression of ST5a (Bol026200), CYP81F1 (Bol028913, Bol028914) and CYP81F4 genes was associated with significantly higher GBS accumulation under MeJA treatment compared to controls in all three subspecies. CYP81F4 genes, trans-activated by MYB34 genes, were expressed at remarkably high levels in all three subspecies under MeJA treatment, which also induced in higher indolic NGBS accumulation in all three subspecies. Remarkably higher expression of MYB28 (Bol036286), ST5b, ST5c, AOP2, FMOGS-OX5 (Bol031350) and GSL-OH (Bol033373) was associated with much higher contents of aliphatic glucosinolates in kale leaves compared to the other two subspecies. The genes expressed highly could be utilized in strategies to selectively increase glucosinolate compounds in B. oleracea subspecies. These results promote efforts to develop genotypes of B. oleracea and other species with enhanced levels of desired glucosinolates.

Comparative Transcriptome Analyses Reveal a Special Glucosinolate Metabolism Mechanism in Brassica alboglabra Sprouts

Brassica sprouts contain abundant phytochemicals, especially glucosinolates (GSs). Various methods have been used to enhance GS content in sprouts. However, the molecular basis of GS metabolism in sprouts remains an open question. Here we employed RNA-seq analysis to compare the transcriptomes of high-GS (JL-08) and low-GS (JL-09) Brassica alboglabra sprouts. Paired-end Illumina RNA-seq reads were generated and mapped to the Brassica oleracea reference genome. The differentially expressed genes were analyzed between JL-08 and JL-09. Among these, 1477 genes were up-regulated and 1239 down-regulated in JL-09 compared with JL-08. Enrichment analysis of these differentially expressed genes showed that the GS biosynthesis had the smallest enrichment factor and the highest Q-value of all metabolic pathways in Kyoto Encyclopedia of Genes and Genomes database, indicating the main metabolic difference between JL-08 and JL-09 is the GS biosynthetic pathway. Thirty-seven genes of the sequenced data were annotated as putatively involved in GS biosynthesis, degradation, and regulation, of which 11 were differentially expressed in JL-08 and JL-09. The expression level of GS degradation enzyme myrosinase in high-GS JL-08 was lower compared with low-GS JL-09. Surprisingly, in high-GS JL-08, the expression levels of GS biosynthesis genes were also lower than those in low-GS JL-09. As the GS contents in sprouts are determined by dynamic equilibrium of seed stored GS mobilization, de novo synthesis, degradation, and extra transport, the result of this study leads us to suggest that efforts to increase GS content should focus on either raising GS content in seeds or decreasing myrosinase activity, rather than improving the expression level of GS biosynthesis genes in sprouts.