Glucosinolate biosynthetic genes in Brassica rapa

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.

User Guidelines for the Brassica Database: BRAD

The genome sequence of Brassica rapa was first released in 2011. Since then, further Brassica genomes have been sequenced or are undergoing sequencing. It is therefore necessary to develop tools that help users to mine information from genomic data efficiently. This will greatly aid scientific exploration and breeding application, especially for those with low levels of bioinformatic training. Therefore, the Brassica database (BRAD) was built to collect, integrate, illustrate, and visualize Brassica genomic datasets. BRAD provides useful searching and data mining tools, and facilitates the search of gene annotation datasets, syntenic or non-syntenic orthologs, and flanking regions of functional genomic elements. It also includes genome-analysis tools such as BLAST and GBrowse. One of the important aims of BRAD is to build a bridge between Brassica crop genomes with the genome of the model species Arabidopsis thaliana, thus transferring the bulk of A. thaliana gene study information for use with newly sequenced Brassica crops.

Aromatic Glucosinolate Biosynthesis Pathway in Barbarea vulgaris and its Response to Plutella xylostella Infestation

The inducibility of the glucosinolate resistance mechanism is an energy-saving strategy for plants, but whether induction would still be triggered by glucosinolate-tolerant Plutella xylostella (diamondback moth, DBM) after a plant had evolved a new resistance mechanism (e.g., saponins in Barbara vulgaris) was unknown. In B. vulgaris, aromatic glucosinolates derived from homo-phenylalanine are the dominant glucosinolates, but their biosynthesis pathway was unclear. In this study, we used G-type (pest-resistant) and P-type (pest-susceptible) B. vulgaris to compare glucosinolate levels and the expression profiles of their biosynthesis genes before and after infestation by DBM larvae. Two different stereoisomers of hydroxylated aromatic glucosinolates are dominant in G- and P-type B. vulgaris, respectively, and are induced by DBM. The transcripts of genes in the glucosinolate biosynthesis pathway and their corresponding transcription factors were identified from an Illumina dataset of G- and P-type B. vulgaris. Many genes involved or potentially involved in glucosinolate biosynthesis were induced in both plant types. The expression patterns of six DBM induced genes were validated by quantitative PCR (qPCR), while six long-fragment genes were validated by molecular cloning. The core structure biosynthetic genes showed high sequence similarities between the two genotypes. In contrast, the sequence identity of two apparent side chain modification genes, the SHO gene in the G-type and the RHO in P-type plants, showed only 77.50% identity in coding DNA sequences and 65.48% identity in deduced amino acid sequences. The homology to GS-OH in Arabidopsis, DBM induction of the transcript and a series of qPCR and glucosinolate analyses of G-type, P-type and F1 plants indicated that these genes control the production of S and R isomers of 2-hydroxy-2-phenylethyl glucosinolate. These glucosinolates were significantly induced by P. xylostella larvae in both the susceptiple P-type and the resistant G-type, even though saponins are the main DBM-resistance causing metabolites in G-type plants. Indol-3-ylmethylglucosinolate was induced in the G-type only. These data will aid our understanding of the biosynthesis and induction of aromatic glucosinolates at the molecular level and also increase our knowledge of the complex mechanisms underpinning defense induction in plants.

Biofortification of oilseed Brassica juncea with the anti-cancer compound glucoraphanin by suppressing GSL-ALK gene family

Glucosinolates are amino acids derived secondary metabolites, invariably present in Brassicales, which have huge health and agricultural benefits. Sulphoraphane, the breakdown product of glucosinolate glucoraphanin is known to posses anti-cancer properties. AOP (2-oxoglutarate-dependent dioxygenases) or GSL-ALK enzyme catalyzes the conversion of desirable glucoraphanin to deleterious gluconapin and progoitrin, which are present in very high amounts in most of the cultivable Brassica species including Brassica juncea. In this study we showed that B. juncea encodes four functional homologs of GSL-ALK gene and constitutive silencing of GSL-ALK homologs resulted in accumulation of glucoraphanin up to 43.11 μmoles g(-1) DW in the seeds with a concomitant reduction in the anti-nutritional glucosinolates. Glucoraphanin content was found remarkably high in leaves as well as sprouts of the transgenic lines. Transcript quantification of high glucoraphanin lines confirmed significant down-regulation of GSL-ALK homologs. Growth and other seed quality parameters of the transgenic lines did not show drastic difference, compared to the untransformed control. High glucoraphanin lines also showed higher resistance towards stem rot pathogen Sclerotinia sclerotiorum. Our results suggest that metabolic engineering of GSL-ALK has huge potential for enriching glucoraphanin content, and improve the oil quality and vegetable value of Brassica crops.

Three genes encoding AOP2, a protein involved in aliphatic glucosinolate biosynthesis, are differentially expressed in Brassica rapa

The glucosinolate biosynthetic gene AOP2 encodes an enzyme that plays a crucial role in catalysing the conversion of beneficial glucosinolates into anti-nutritional ones. In Brassica rapa, three copies of BrAOP2 have been identified, but their function in establishing the glucosinolate content of B. rapa is poorly understood. Here, we used phylogenetic and gene structure analyses to show that BrAOP2 proteins have evolved via a duplication process retaining two highly conserved domains at the N-terminal and C-terminal regions, while the middle part has experienced structural divergence. Heterologous expression and in vitro enzyme assays and Arabidopsis mutant complementation studies showed that all three BrAOP2 genes encode functional BrAOP2 proteins that convert the precursor methylsulfinyl alkyl glucosinolate to the alkenyl form. Site-directed mutagenesis showed that His356, Asp310, and Arg376 residues are required for the catalytic activity of one of the BrAOP2 proteins (BrAOP2.1). Promoter-β-glucuronidase lines revealed that the BrAOP2.3 gene displayed an overlapping but distinct tissue- and cell-specific expression profile compared with that of the BrAOP2.1 and BrAOP2.2 genes. Quantitative real-time reverse transcription-PCR assays demonstrated that BrAOP2.1 showed a slightly different pattern of expression in below-ground tissue at the seedling stage and in the silique at the reproductive stage compared with BrAOP2.2 and BrAOP2.3 genes in B. rapa. Taken together, our results revealed that all three BrAOP2 paralogues are active in B. rapa but have functionally diverged.

Novel glucosinolate composition lacking 4-methylthio-3-butenyl glucosinolate in Japanese white radish (Raphanus sativus L.)

Genetic analysis and gene mapping of the 4-methylthio-3-butenyl glucosinolate-less trait of white radish were performed and a white radish cultivar with new glucosinolate composition was developed. A spontaneous mutant having significantly low 4-methylthio-3-butenyl glucosinolate (4MTB-GSL) content was identified from a landrace of Japanese white radish (Raphanus sativus L.) through intensive evaluation of glucosinolate profiles of 632 lines including genetic resources and commercial cultivars using high-performance liquid chromatography (HPLC) analysis. A line lacking 4MTB-GSL was developed using the selected mutant as a gene source. Genetic analyses of F1, F2, and BC1F1 populations of this line suggested that the 4MTB-GSL-less trait is controlled by a single recessive allele. Using SNP and SCAR markers, 96 F2 plants were genotyped, and a linkage map having nine linkage groups with a total map distance of 808.3 cM was constructed. A gene responsible for the 4MTB-GSL-less trait was mapped between CL1753 and CL5895 at the end of linkage group 1. The genetic distance between these markers was 4.2 cM. By selfing and selection of plants lacking 4MTB-GSL, a new cultivar, 'Daikon parental line No. 5', was successfully developed. This cultivar was characterized by glucoerucin, which accounted for more than 90% of the total glucosinolates (GSLs). The total GSL content in roots was ca. 12 μmol/g DW, significantly lower than those in common white radish cultivars. Significance of this line in radish breeding is discussed.

Detoxification function of the Arabidopsis sulphotransferase AtSOT12 by sulphonation of xenobiotics

Cytosolic sulphotransferases have been implicated in inactivation of endogenous steroid hormones and detoxification of xenobiotics in human and animals. Yet, the function of plant sulphotransferases in xenobiotic sulphonation and detoxification has not been reported. In this study, we show that the Arabidopsis sulphotransferase AtSOT12 could sulphonate the bacterial-produced toxin cycloheximide. Loss-of-function mutant sot12 exhibited hypersensitive phenotype to cycloheximide, and expression of AtSOT12 protein in yeast cells conferred resistance to this toxic compound. AtSOT12 exhibited broad specificity and could sulphonate a variety of xenobiotics including phenolic and polycyclic compounds. Enzyme kinetics analysis indicated that AtSOT12 has different selectivity for simple phenolics with different side chains, and the position of the side chain in the simple phenolic compounds affects substrate binding affinity and catalytic efficiency. We proposed that the broad specificity and induced production of AtSOT12 may have rendered this enzyme to not only modify endogenous molecules such as salicylic acid as we previously reported, but also sulphonate pathogen-produced toxic small molecules to protect them from infection. Sulphonation of small molecules in plants may constitute a rapid way to inactivate or change the physiochemical properties of biologically active molecules that could have profound effects on plant growth, development and defence.

Identification and expression analysis of glucosinolate biosynthetic genes and estimation of glucosinolate contents in edible organs of Brassica oleracea subspecies

Glucosinolates are anti-carcinogenic, anti-oxidative biochemical compounds that defend plants from insect and microbial attack. Glucosinolates are abundant in all cruciferous crops, including all vegetable and oilseed Brassica species. Here, we studied the expression of glucosinolate biosynthesis genes and determined glucosinolate contents in the edible organs of a total of 12 genotypes of Brassica oleracea: three genotypes each from cabbage, kale, kohlrabi and cauliflower subspecies. Among the 81 genes analyzed by RT-PCR, 19 are transcription factor-related, two different sets of 25 genes are involved in aliphatic and indolic biosynthesis pathways and the rest are breakdown-related. The expression of glucosinolate-related genes in the stems of kohlrabi was remarkably different compared to leaves of cabbage and kale and florets of cauliflower as only eight genes out of 81 were expressed in the stem tissues of kohlrabi. In the stem tissue of kohlrabi, only one aliphatic transcription factor-related gene, Bol036286 (MYB28) and one indolic transcription factor-related gene, Bol030761 (MYB51), were expressed. The results indicated the expression of all genes is not essential for glucosinolate biosynthesis. Using HPLC analysis, a total of 16 different types of glucosinolates were identified in four subspecies, nine of them were aliphatic, four of them were indolic and one was aromatic. Cauliflower florets measured the highest number of 14 glucosinolates. Among the aliphatic glucosinolates, only gluconapin was found in the florets of cauliflower. Glucoiberverin and glucobrassicanapin contents were the highest in the stems of kohlrabi. The indolic methoxyglucobrassicin and aromatic gluconasturtiin accounted for the highest content in the florets of cauliflower. A further detailed investigation and analyses is required to discern the precise roles of each of the genes for aliphatic and indolic glucosinolate biosynthesis in the edible organs.

Carotenoid biosynthetic genes in Brassica rapa: comparative genomic analysis, phylogenetic analysis, and expression profiling

Background

Carotenoids are isoprenoid compounds synthesized by all photosynthetic organisms. Despite much research on carotenoid biosynthesis in the model plant Arabidopsis thaliana, there is a lack of information on the carotenoid pathway in Brassica rapa. To better understand its carotenoid biosynthetic pathway, we performed a systematic analysis of carotenoid biosynthetic genes at the genome level in B. rapa.

Results

We identified 67 carotenoid biosynthetic genes in B. rapa, which were orthologs of the 47 carotenoid genes in A. thaliana. A high level of synteny was observed for carotenoid biosynthetic genes between A. thaliana and B. rapa. Out of 47 carotenoid biosynthetic genes in A. thaliana, 46 were successfully mapped to the 10 B. rapa chromosomes, and most of the genes retained more than one copy in B. rapa. The gene expansion was caused by the whole-genome triplication (WGT) event experienced by Brassica species. An expression analysis of the carotenoid biosynthetic genes suggested that their expression levels differed in root, stem, leaf, flower, callus, and silique tissues. Additionally, the paralogs of each carotenoid biosynthetic gene, which were generated from the WGT in B. rapa, showed significantly different expression levels among tissues, suggesting differentiated functions for these multi-copy genes in the carotenoid pathway.

Conclusions

This first systematic study of carotenoid biosynthetic genes in B. rapa provides insights into the carotenoid metabolic mechanisms of Brassica crops. In addition, a better understanding of carotenoid biosynthetic genes in B. rapa will contribute to the development of conventional and transgenic B. rapa cultivars with enriched carotenoid levels in the future.

Electronic supplementary material

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

The radish genome and comprehensive gene expression profile of tuberous root formation and development

Understanding the processes that regulate plant sink formation and development at the molecular level will contribute to the areas of crop breeding, food production and plant evolutionary studies. We report the annotation and analysis of the draft genome sequence of the radish Raphanus sativus var. hortensis (long and thick root radish) and transcriptome analysis during root development. Based on the hybrid assembly approach of next-generation sequencing, a total of 383 Mb (N50 scaffold: 138.17 kb) of sequences of the radish genome was constructed containing 54,357 genes. Syntenic and phylogenetic analyses indicated that divergence between Raphanus and Brassica coincide with the time of whole genome triplication (WGT), suggesting that WGT triggered diversification of Brassiceae crop plants. Further transcriptome analysis showed that the gene functions and pathways related to carbohydrate metabolism were prominently activated in thickening roots, particularly in cell proliferating tissues. Notably, the expression levels of sucrose synthase 1 (SUS1) were correlated with root thickening rates. We also identified the genes involved in pungency synthesis and their transcription factors.

Genome-wide identification of SSR and SNP markers from the non-heading Chinese cabbage for comparative genomic analyses

BACKGROUND

Non-heading Chinese cabbage (NHCC), belonging to Brassica, is an important leaf vegetable in Asia. Although genetic analyses have been performed through conventional selection and breeding efforts, the domestication history of NHCC and the genetics underlying its morphological diversity remain unclear. Thus, the reliable molecular markers representative of the whole genome are required for molecular-assisted selection in NHCC.

RESULTS

A total of 20,836 simple sequence repeats (SSRs) were detected in NHCC, containing repeat types from mononucleotide to nonanucleotide. The average density was 62.93 SSRs/Mb. In gene regions, 5,435 SSRs were identified in 4,569 genes. A total of 5,008 primer pairs were designed, and 74 were randomly selected for validation. Among these, 60 (81.08%) were polymorphic in 18 Cruciferae. The number of polymorphic bands ranged from two to five, with an average of 2.70 for each primer. The average values of the polymorphism information content, observed heterozygosity, Hardy-Weinberg equilibrium, and Shannon's information index were 0.2970, 0.4136, 0.5706, and 0.5885, respectively. Four clusters were classified according to the unweighted pair-group method with arithmetic average cluster analysis of 18 genotypes. In addition, a total of 1,228,979 single nucleotide polymorphisms (SNPs) were identified in the NHCC through a comparison with the genome of Chinese cabbage, and the average SNP density in the whole genome was 4.33/Kb. The number of SNPs ranged from 341,939 to 591,586 in the 10 accessions, and the average heterozygous SNPs ratio was ~42.53%. All analyses showed these markers were high quality and reliable. Therefore, they could be used in the construction of a linkage map and for genetic diversity studies for NHCC in future.

CONCLUSIONS

This is the first systematic and comprehensive analysis and identification of SSRs in NHCC and 17 species. The development of a large number of SNP and SSR markers was successfully achieved for NHCC. These novel markers are valuable for constructing genetic linkage maps, comparative genome analysis, quantitative trait locus (QTL) mapping, genome-wide association studies, and marker-assisted selection in NHCC breeding system research.

Deciphering allelic variations for seed glucosinolate traits in oilseed mustard (Brassica juncea) using two bi-parental mapping populations

QTL mapping by two DH mapping populations deciphered allelic variations for five different seed glucosinolate traits in B. juncea. Allelic variations for five different seed glucosinolate (GS) traits, namely % propyl, % butyl, % pentyl, aliphatics and total GS content were studied through QTL analysis using two doubled haploid (DH) mapping populations. While the high GS parents in two populations differed in their profiles of seed aliphatic GS, the low GS parents were similar. Phenotypic data of seed GS traits from three environments of the two populations were subjected to QTL analysis. The first population (referred to as DE population) detected a total of 60 QTL from three environments which upon intra-population meta-QTL analysis were merged to 17 S-QTL (Stable QTL) and 15 E-QTL (Environment QTL). The second population (referred to as VH population) detected 58 QTL from the three environments that were merged to 15S-QTL and 16E-QTL. In both the populations, majority of S-QTL were detected as major QTL. Inter-population meta-analysis identified three C-QTL (consensus QTL) formed by merging major QTL from the two populations. Candidate genes of GS pathway were co-localized to the QTL regions either through genetic mapping or through in silico comparative analysis. Parental allelic variants of QTL or of the co-mapped candidate gene(s) were determined on the basis of the significantly different R (2) values of the component QTL from the two populations which were merged to form C-QTL. The results of the study are significant for marker-assisted transfer of the low GS trait and also for developing lines with lower GS than are present in Brassica juncea.

Genomic origin, expression differentiation and regulation of multiple genes encoding CYP83A1, a key enzyme for core glucosinolate biosynthesis, from the allotetraploid Brassica juncea

The multiple BjuCYP83A1 genes formed as a result of polyploidy have retained cell-, tissue-, and condition-specific transcriptional sub-functionalization to control the complex aliphatic glucosinolates biosynthesis in the allotetraploid Brassica juncea. Glucosinolates along with their breakdown products are associated with diverse roles in plant metabolism, plant defense and animal nutrition. CYP83A1 is a key enzyme that oxidizes aliphatic aldoximes to aci-nitro compounds in the complex aliphatic glucosinolate biosynthetic pathway. In this study, we reported the isolation of four CYP83A1 genes named BjuCYP83A1-1, -2, -3, and -4 from allotetraploid Brassica juncea (AABB genome), an economically important oilseed crop of Brassica genus. The deduced BjuCYP83A1 proteins shared 85.7-88.4 % of sequence identity with A. thaliana AtCYP83A1 and 84.2-95.8 % among themselves. Phylogenetic and divergence analysis revealed that the four BjuCYP83A1 proteins are evolutionary conserved and have evolved via duplication and hybridization of two relatively simpler diploid Brassica genomes namely B. rapa (AA genome) and B. nigra (BB genome), and have retained high level of sequence conservation following allopolyploidization. Ectopic over-expression of BjuCYP83A1-1 in A. thaliana showed that it is involved mainly in the synthesis of C4 aliphatic glucosinolates. Detailed expression analysis using real-time qRT-PCR in B. juncea and PromoterBjuCYP83A1-GUS lines in A. thaliana confirmed that the four BjuCYP83A1 genes have retained ubiquitous, overlapping but distinct expression profiles in different tissue and cell types of B. juncea, and in response to various elicitor treatments and environmental conditions. Taken together, this study demonstrated that transcriptional sub-functionalization and coordinated roles of multiple BjuCYP83A1 genes control the biosynthesis of aliphatic glucosinolates in the allotetraploid B. juncea, and provide a framework for metabolic engineering of aliphatic glucosinolates in economically important Brassica species.

Lineage-specific evolution of Methylthioalkylmalate synthases (MAMs) involved in glucosinolates biosynthesis

Methylthioalkylmalate synthases (MAMs) encoded by MAM genes are central to the diversification of the glucosinolates, which are important secondary metabolites in Brassicaceae species. However, the evolutionary pathway of MAM genes is poorly understood. We analyzed the phylogenetic and synteny relationships of MAM genes from 13 sequenced Brassicaceae species. Based on these analyses, we propose that the syntenic loci of MAM genes, which underwent frequent tandem duplications, divided into two independent lineage-specific evolution routes and were driven by positive selection after the divergence from Aethionema arabicum. In the lineage I species Capsella rubella, Camelina sativa, Arabidopsis lyrata, and A. thaliana, the MAM loci evolved three tandem genes encoding enzymes responsible for the biosynthesis of aliphatic glucosinolates with different carbon chain-lengths. In lineage II species, the MAM loci encode enzymes responsible for the biosynthesis of short-chain aliphatic glucosinolates. Our proposed model of the evolutionary pathway of MAM genes will be useful for understanding the specific function of these genes in Brassicaceae species.

The multi-protein family of sulfotransferases in plants: composition, occurrence, substrate specificity, and functions

All members of the sulfotransferase (SOT, EC 2.8.2.-) protein family transfer a sulfuryl group from the donor 3'-phosphoadenosine 5'-phosphosulfate (PAPS) to an appropriate hydroxyl group of several classes of substrates. The primary structure of these enzymes is characterized by a histidine residue in the active site, defined PAPS binding sites and a longer SOT domain. Proteins with this SOT domain occur in all organisms from all three domains, usually as a multi-protein family. Arabidopsis thaliana SOTs, the best characterized SOT multi-protein family, contains 21 members. The substrates for several plant enzymes have already been identified, such as glucosinolates, brassinosteroids, jasmonates, flavonoids, and salicylic acid. Much information has been gathered on desulfo-glucosinolate (dsGl) SOTs in A. thaliana. The three cytosolic dsGl SOTs show slightly different expression patterns. The recombinant proteins reveal differences in their affinity to indolic and aliphatic dsGls. Also the respective recombinant dsGl SOTs from different A. thaliana ecotypes differ in their kinetic properties. However, determinants of substrate specificity and the exact reaction mechanism still need to be clarified. Probably, the three-dimensional structures of more plant proteins need to be solved to analyze the mode of action and the responsible amino acids for substrate binding. In addition to A. thaliana, more plant species from several families need to be investigated to fully elucidate the diversity of sulfated molecules and the way of biosynthesis catalyzed by SOT enzymes.

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

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

Myrosin idioblast cell fate and development are regulated by the Arabidopsis transcription factor FAMA, the auxin pathway, and vesicular trafficking

Crucifer shoots harbor a glucosinolate-myrosinase system that defends against insect predation. Arabidopsis thaliana myrosinase (thioglucoside glucohydrolase [TGG]) accumulates in stomata and in myrosin idioblasts (MIs). This work reports that the basic helix-loop-helix transcription factor FAMA that is key to stomatal development is also expressed in MIs. The loss of FAMA function abolishes MI fate as well as the expression of the myrosinase genes TGG1 and TGG2. MI cells have previously been reported to be located in the phloem. Instead, we found that MIs arise from the ground meristem rather than provascular tissues and thus are not homologous with phloem. Moreover, MI patterning and morphogenesis are abnormal when the function of the ARF-GEF gene GNOM is lost as well as when auxin efflux and vesicular trafficking are chemically disrupted. Stomata and MI cells constitute part of a wider system that reduces plant predation, the so-called "mustard oil bomb," in which vacuole breakage in cells harboring myrosinase and glucosinolate yields a brew toxic to many animals, especially insects. This identification of the gene that confers the fate of MIs, as well as stomata, might facilitate the development of strategies for engineering crops to mitigate predation.

Regulatory network of secondary metabolism in Brassica rapa: insight into the glucosinolate pathway

Brassica rapa studies towards metabolic variation have largely been focused on the profiling of the diversity of metabolic compounds in specific crop types or regional varieties, but none aimed to identify genes with regulatory function in metabolite composition. Here we followed a genetical genomics approach to identify regulatory genes for six biosynthetic pathways of health-related phytochemicals, i.e carotenoids, tocopherols, folates, glucosinolates, flavonoids and phenylpropanoids. Leaves from six weeks-old plants of a Brassica rapa doubled haploid population, consisting of 92 genotypes, were profiled for their secondary metabolite composition, using both targeted and LC-MS-based untargeted metabolomics approaches. Furthermore, the same population was profiled for transcript variation using a microarray containing EST sequences mainly derived from three Brassica species: B. napus, B. rapa and B. oleracea. The biochemical pathway analysis was based on the network analyses of both metabolite QTLs (mQTLs) and transcript QTLs (eQTLs). Co-localization of mQTLs and eQTLs lead to the identification of candidate regulatory genes involved in the biosynthesis of carotenoids, tocopherols and glucosinolates. We subsequently focused on the well-characterized glucosinolate pathway and revealed two hotspots of co-localization of eQTLs with mQTLs in linkage groups A03 and A09. Our results indicate that such a large-scale genetical genomics approach combining transcriptomics and metabolomics data can provide new insights into the genetic regulation of metabolite composition of Brassica vegetables.

A comprehensive set of transcript sequences of the heavy metal hyperaccumulator Noccaea caerulescens

Noccaea caerulescens is an extremophile plant species belonging to the Brassicaceae family. It has adapted to grow on soils containing high, normally toxic, concentrations of metals such as nickel, zinc, and cadmium. Next to being extremely tolerant to these metals, it is one of the few species known to hyperaccumulate these metals to extremely high concentrations in their aboveground biomass. In order to provide additional molecular resources for this model metal hyperaccumulator species to study and understand the mechanism of adaptation to heavy metal exposure, we aimed to provide a comprehensive database of transcript sequences for N. caerulescens. In this study, 23,830 transcript sequences (isotigs) with an average length of 1025 bp were determined for roots, shoots and inflorescences of N. caerulescens accession "Ganges" by Roche GS-FLEX 454 pyrosequencing. These isotigs were grouped into 20,378 isogroups, representing potential genes. This is a large expansion of the existing N. caerulescens transcriptome set consisting of 3705 unigenes. When translated and compared to a Brassicaceae proteome set, 22,232 (93.2%) of the N. caerulescens isotigs (corresponding to 19,191 isogroups) had a significant match and could be annotated accordingly. Of the remaining sequences, 98 isotigs resembled non-plant sequences and 1386 had no significant similarity to any sequence in the GenBank database. Among the annotated set there were many isotigs with similarity to metal homeostasis genes or genes for glucosinolate biosynthesis. Only for transcripts similar to Metallothionein3 (MT3), clear evidence for an additional copy was found. This comprehensive set of transcripts is expected to further contribute to the discovery of mechanisms used by N. caerulescens to adapt to heavy metal exposure.

The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes

Polyploidization has provided much genetic variation for plant adaptive evolution, but the mechanisms by which the molecular evolution of polyploid genomes establishes genetic architecture underlying species differentiation are unclear. Brassica is an ideal model to increase knowledge of polyploid evolution. Here we describe a draft genome sequence of Brassica oleracea, comparing it with that of its sister species B. rapa to reveal numerous chromosome rearrangements and asymmetrical gene loss in duplicated genomic blocks, asymmetrical amplification of transposable elements, differential gene co-retention for specific pathways and variation in gene expression, including alternative splicing, among a large number of paralogous and orthologous genes. Genes related to the production of anticancer phytochemicals and morphological variations illustrate consequences of genome duplication and gene divergence, imparting biochemical and morphological variation to B. oleracea. This study provides insights into Brassica genome evolution and will underpin research into the many important crops in this genus.

Nitrates and glucosinolates as strong determinants of the nutritional quality in rocket leafy salads

Rocket is an important leafy vegetable crop and a good source of antioxidants and anticancer molecules such as glucosinolates and other sulfur compounds. Rocket is also a hyper-accumulator of nitrates which have been considered for long time the main factors that cause gastro-intestinal cancer. In this review, the content of these compounds in rocket tissues and their levels at harvest and during storage are discussed. Moreover, the effect of these compounds in preventing or inducing human diseases is also highlighted. This review provides an update to all the most recent studies carried out on rocket encouraging the consumption of this leafy vegetable to reduce the risk of contracting cancer and other cardiovascular diseases.

Identification of metabolic QTLs and candidate genes for glucosinolate synthesis in Brassica oleracea leaves, seeds and flower buds

Glucosinolates are major secondary metabolites found in the Brassicaceae family. These compounds play an essential role in plant defense against biotic and abiotic stresses, but more interestingly they have beneficial effects on human health. We performed a genetic analysis in order to identify the genome regions regulating glucosinolates biosynthesis in a DH mapping population of Brassica oleracea. In order to obtain a general overview of regulation in the whole plant, analyses were performed in the three major organs where glucosinolates are synthesized (leaves, seeds and flower buds). Eighty two significant QTLs were detected, which explained a broad range of variability in terms of individual and total glucosinolate (GSL) content. A meta-analysis rendered eighteen consensus QTLs. Thirteen of them regulated more than one glucosinolate and its content. In spite of the considerable variability of glucosinolate content and profiles across the organ, some of these consensus QTLs were identified in more than one tissue. Consensus QTLs control the GSL content by interacting epistatically in complex networks. Based on in silico analysis within the B. oleracea genome along with synteny with Arabidopsis, we propose seven major candidate loci that regulate GSL biosynthesis in the Brassicaceae family. Three of these loci control the content of aliphatic GSL and four of them control the content of indolic glucosinolates. GSL-ALK plays a central role in determining aliphatic GSL variation directly and by interacting epistatically with other loci, thus suggesting its regulatory effect.

A Multi-Layered Screening Method to Identify Plant Regulatory Genes

We used a seven-step process to identify genes involved in glucosinolate biosynthesis and metabolism in the Chinese cabbage (Brassica rapa). We constructed an annotated data set with 34,570 unigenes from B. rapa and predicted 11,526 glucosinolate-related candidate genes using expression profiles generated across nine stages of development on a 47k-gene microarray. Using our multi-layered screening method, we screened 392 transcription factors, 843 pathway genes, and 4,162 ortholog genes associated with glucosinolate-related biosynthesis. Finally, we identified five genes by comparison of the pathway-network genes including the transcription-factor genes and the ortholog-ontology genes. The five genes were anchored to the chromosomes of B. rapa to characterize their genetic-map positions, and phylogenetic reconstruction with homologous genes was performed. These anchored genes were verified by reverse-transcription polymerase chain reaction. While the five genes identified by our multi-layered screen require further characterization and validation, our study demonstrates the power of multi-layered screening after initial identification of genes on microarrays.

Genome triplication drove the diversification of Brassica plants

The genus Brassica belongs to the plant family Brassicaceae, which includes many important crop species that are used as oilseed, condiments, or vegetables throughout the world. Brassica plants comprise many diverse species, and each species contains rich morphotypes showing extreme traits. Brassica species experienced an extra whole genome triplication (WGT) event compared with the model plant Arabidopsis thaliana. Whole genome sequencing of the Brassica species Brassica rapa, Brassica oleracea and others demonstrated that WGT plays an important role in the speciation and morphotype diversification of Brassica plants. Comparative genomic analysis based on the genome sequences of B. rapa and A. thaliana clearly identified the WGT event and further demonstrated that the translocated Proto-Calepine Karyotype (tPCK, n=7) was the diploid ancestor of the three subgenomes in B. rapa. Following WGT, subsequent extensive genome fractionation, block reshuffling and chromosome reduction accompanied by paleocentromere descent from the three tPCK subgenomes during the rediploidization process produced stable diploid species. Genomic rearrangement of the diploid species and their hybridization then contributed to Brassica speciation. The subgenome dominance effect and biased gene retention, such as the over-retention of auxin-related genes after WGT, promoted functional gene evolution and thus propelled the expansion of rich morphotypes in the Brassica species. In conclusion, the WGT event initiated subsequent genomic and gene-level evolution, which further drove Brassica speciation and created rich morphotypes in each species.

De novo transcriptome sequencing of radish (Raphanus sativus L.) and analysis of major genes involved in glucosinolate metabolism

Background

Radish (Raphanus sativus L.), is an important root vegetable crop worldwide. Glucosinolates in the fleshy taproot significantly affect the flavor and nutritional quality of radish. However, little is known about the molecular mechanisms underlying glucosinolate metabolism in radish taproots. The limited availability of radish genomic information has greatly hindered functional genomic analysis and molecular breeding in radish.

Results

In this study, a high-throughput, large-scale RNA sequencing technology was employed to characterize the de novo transcriptome of radish roots at different stages of development. Approximately 66.11 million paired-end reads representing 73,084 unigenes with a N50 length of 1,095 bp, and a total length of 55.73 Mb were obtained. Comparison with the publicly available protein database indicates that a total of 67,305 (about 92.09% of the assembled unigenes) unigenes exhibit similarity (e –value ≤ 1.0e^-5) to known proteins. The functional annotation and classification including Gene Ontology (GO), Clusters of Orthologous Group (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the main activated genes in radish taproots are predominately involved in basic physiological and metabolic processes, biosynthesis of secondary metabolite pathways, signal transduction mechanisms and other cellular components and molecular function related terms. The majority of the genes encoding enzymes involved in glucosinolate (GS) metabolism and regulation pathways were identified in the unigene dataset by targeted searches of their annotations. A number of candidate radish genes in the glucosinolate metabolism related pathways were also discovered, from which, eight genes were validated by T-A cloning and sequencing while four were validated by quantitative RT-PCR expression profiling.

Conclusions

The ensuing transcriptome dataset provides a comprehensive sequence resource for molecular genetics research in radish. It will serve as an important public information platform to further understanding of the molecular mechanisms involved in biosynthesis and metabolism of the related nutritional and flavor components during taproot formation in radish.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-836) contains supplementary material, which is available to authorized users.

Testing local host adaptation and phenotypic plasticity in a herbivore when alternative related host plants occur sympatrically

Host race formation in phytophagous insects can be an early stage of adaptive speciation. However, the evolution of phenotypic plasticity in host use is another possible outcome. Using a reciprocal transplant experiment we tested the hypothesis of local adaptation in the aphid Brevicoryne brassicae. Aphid genotypes derived from two sympatric host plants, Brassica oleracea and B. campestris, were assessed in order to measure the extent of phenotypic plasticity in morphological and life history traits in relation to the host plants. We obtained an index of phenotypic plasticity for each genotype. Morphological variation of aphids was summarized by principal components analysis. Significant effects of recipient host on morphological variation and life history traits (establishment, age at first reproduction, number of nymphs, and intrinsic growth rate) were detected. We did not detected genotype × host plant interaction; in general the genotypes developed better on B. campestris, independent of the host plant species from which they were collected. Therefore, there was no evidence to suggest local adaptation. Regarding plasticity, significant differences among genotypes in the index of plasticity were detected. Furthermore, significant selection on PC1 (general aphid body size) on B. campestris, and on PC1 and PC2 (body length relative to body size) on B. oleracea was detected. The elevation of the reaction norm of PC1 and the slope of the reaction norm for PC2 (i.e., plasticity) were under directional selection. Thus, host plant species constitute distinct selective environments for B. brassicae. Aphid genotypes expressed different phenotypes in response to the host plant with low or nil fitness costs. Phenotypic plasticity and gene flow limits natural selection for host specialization promoting the maintenance of genetic variation in host exploitation.

Four genes encoding MYB28, a major transcriptional regulator of the aliphatic glucosinolate pathway, are differentially expressed in the allopolyploid Brassica juncea

Glucosinolates are Capparales-specific secondary metabolites that have immense potential in human health and agriculture. Unlike Arabidopsis thaliana, our knowledge about glucosinolate regulators in the Brassica crops is sparse. In the current study, four MYB28 homologues were identified (BjuMYB28-1,-2,-3,-4) from the polyploid Brassica juncea, and the effects of allopolyploidization on the divergence of gene sequence, structure, function, and expression were assessed. The deduced protein sequences of the four BjuMYB28 genes showed 76.1-83.1% identity with the Arabidopsis MYB28. Phylogenetic analysis revealed that the four BjuMYB28 proteins have evolved via the hybridization and duplication processes forming the B. juncea genome (AABB) from B. rapa (AA) and B. nigra (BB), while retaining high levels of sequence conservation. Mutant complementation and over-expression studies in A. thaliana showed that all four BjuMYB28 genes encode functional MYB28 proteins and resulted in similar aliphatic glucosinolate composition and content. Detailed expression analysis using qRT-PCR assays and promoter-GUS lines revealed that the BjuMYB28 genes have both tissue- and cell-specific expression partitioning in B. juncea. The two B-genome origin BjuMYB28 genes had more abundant transcripts during the early stages of plant development than the A-genome origin genes. However, with the onset of the reproductive phase, expression levels of all four BjuMYB28 increased significantly, which may be necessary for producing and maintaining high amounts of aliphatic glucosinolates during the later stages of plant development. Taken together, our results suggest that the four MYB28 genes are differentially expressed and regulated in B. juncea to play discrete though overlapping roles in controlling aliphatic glucosinolate biosynthesis.

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.

Quantitative trait loci analysis of non-enzymatic glucosinolate degradation rates in Brassica oleracea during food processing

Epidemiological and mechanistic studies show health-promoting effects of glucosinolates and their breakdown products. In literature, differences in non-enzymatic glucosinolate degradation rates during food processing between different vegetables are described, which provide the basis for studying the genetic effects of this trait and breeding vegetables with high glucosinolate retention during food processing. Non-enzymatic glucosinolate degradation, induced by heat, was studied in a publicly available Brassica oleracea doubled haploid population. Data were modeled to obtain degradation rate constants that were used as phenotypic traits to perform quantitative trait loci (QTL) mapping. Glucosinolate degradation rate constants were determined for five aliphatic and two indolic glucosinolates. Degradation rates were independent of the initial glucosinolate concentration. Two QTL were identified for the degradation rate of the indolic glucobrassicin and one QTL for the degradation of the aliphatic glucoraphanin, which co-localized with one of the QTL for glucobrassicin. Factors within the plant matrix might influence the degradation of different glucosinolates in different genotypes. In addition to genotypic effects, we demonstrated that growing conditions influenced glucosinolate degradation as well. The study identified QTL for glucosinolate degradation, giving the opportunity to breed vegetables with a high retention of glucosinolates during food processing, although the underlying mechanisms remain unknown.

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.

Genetic analysis of health-related secondary metabolites in a Brassica rapa recombinant inbred line population

The genetic basis of the wide variation for nutritional traits in Brassica rapa is largely unknown. A new Recombinant Inbred Line (RIL) population was profiled using High Performance Liquid Chromatography (HPLC) and Nuclear Magnetic Resonance (NMR) analysis to detect quantitative trait loci (QTLs) controlling seed tocopherol and seedling metabolite concentrations. RIL population parent L58 had a higher level of glucosinolates and phenylpropanoids, whereas levels of sucrose, glucose and glutamate were higher in the other RIL population parent, R-o-18. QTL related to seed tocopherol (α-, β-, γ-, δ-, α-/γ- and total tocopherol) concentrations were detected on chromosomes A3, A6, A9 and A10, explaining 11%–35% of the respective variation. The locus on A3 co-locates with the BrVTE1gene, encoding tocopherol cyclase. NMR spectroscopy identified the presence of organic/amino acid, sugar/glucosinolate and aromatic compounds in seedlings. QTL positions were obtained for most of the identified compounds. Compared to previous studies, novel loci were found for glucosinolate concentrations. This work can be used to design markers for marker-assisted selection of nutritional compounds in B. rapa.

MYB transcription factors regulate glucosinolate biosynthesis in different organs of Chinese cabbage (Brassica rapa ssp. pekinensis)

In this study, we investigated the expression of seven MYB transcription factors (a total of 17 genes that included Dof1.1, IQD1-1, MYB28, MYB29, MYB34, MYB51, and MYB122 and their isoforms) involved in aliphatic and indolic glucosinolate (GSL) biosynthesis and analyzed the aliphatic and indolic GSL content in different organs of Chinese cabbage (Brassica rapassp. Pekinensis). MYB28 and MYB29 expression in the stem was dramatically different when compared with the levels in the other organs. MYB34, MYB122, MYB51, Dof1.1, and IQD1-1 showed very low transcript levels among different organs. HPLC analysis showed that the glucosinolates (GSLs) consisted of five aliphatic GSLs (progoitrin, sinigrin, glucoalyssin, gluconapin, and glucobrassicanapin) and four indolic GSLs (4-hydroxyglucobrassicin, glucobrassicin, 4-methoxygluco-brassicin, and neoglucobrassicin). Aliphatic GSLs exhibited 63.3% of the total GSLs content, followed by aromatic GSL (19.0%), indolic GSLs (10%), and unknown GSLs (7.7%) in different organs of Chinese cabbage. The total GSL content of different parts (ranked in descending order) was as follows: seed > flower > young leaves > stem > root > old leaves. The relationship between GSLs accumulation and expression of GSLs biosynthesis MYB TFs genes in different organs may be helpful to understand the mechanism of MYB TFs regulating GSL biosynthesis in Chinese cabbage.

Genetic regulation of glucoraphanin accumulation in Beneforté broccoli

· Diets rich in broccoli (Brassica oleracea var italica) have been associated with maintenance of cardiovascular health and reduction in risk of cancer. These health benefits have been attributed to glucoraphanin that specifically accumulates in broccoli. The development of broccoli with enhanced concentrations of glucoraphanin may deliver greater health benefits. · Three high-glucoraphanin F1 broccoli hybrids were developed in independent programmes through genome introgression from the wild species Brassica villosa. Glucoraphanin and other metabolites were quantified in experimental field trials. Global SNP analyses quantified the differential extent of B. villosa introgression · The high-glucoraphanin broccoli hybrids contained 2.5-3 times the glucoraphanin content of standard hybrids due to enhanced sulphate assimilation and modifications in sulphur partitioning between sulphur-containing metabolites. All of the high-glucoraphanin hybrids possessed an introgressed B. villosa segment which contained a B. villosa Myb28 allele. Myb28 expression was increased in all of the high-glucoraphanin hybrids. Two high-glucoraphanin hybrids have been commercialised as Beneforté broccoli. · The study illustrates the translation of research on glucosinolate genetics from Arabidopsis to broccoli, the use of wild Brassica species to develop cultivars with potential consumer benefits, and the development of cultivars with contrasting concentrations of glucoraphanin for use in blinded human intervention studies.

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

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

Flavin-dependent enzymes in cancer prevention

Statistical studies have demonstrated that various agents may reduce the risk of cancer's development. One of them is activity of flavin-dependent enzymes such as flavin-containing monooxygenase (FMO)(GS-OX1), FAD-dependent 5,10-methylenetetrahydrofolate reductase and flavin-dependent monoamine oxidase. In the last decade, many papers concerning their structure, reaction mechanism and role in the cancer prevention were published. In our work, we provide a more in-depth analysis of flavin-dependent enzymes and their contribution to the cancer prevention. We present the actual knowledge about the glucosinolate synthesized by flavin-containing monooxygenase (FMO)(GS-OX1) and its role in cancer prevention, discuss the influence of mutations in FAD-dependent 5,10-methylenetetrahydrofolate reductase on the cancer risk, and describe FAD as an important cofactor for the demethylation of histons. We also present our views on the role of riboflavin supplements in the prevention against cancer.

Engineering glucosinolates in plants: current knowledge and potential uses

Glucosinolates (GSL) and their derivatives are well known for the characteristic roles they play in plant defense as signaling molecules and as bioactive compounds for human health. More than 130 GSLs have been reported so far, and most of them belong to the Brassicaceae family. Several enzymes and transcription factors involved in the GSL biosynthesis have been studied in the model plant, Arabidopsis, and in a few other Brassica crop species. Recent studies in GSL research have defined the regulation, distribution, and degradation of GSL biosynthetic pathways; however, the underlying mechanism behind transportation of GSLs in plants is still largely unknown. This review highlights the recent advances in the metabolic engineering of GSLs in plants and discusses their potential applications.