Quantitative trait loci for glucosinolate accumulation in Brassica rapa leaves

Construction and Application of an F1-Derived Doubled-Haploid Population and High-Density Genetic Map for Ornamental Kale Breeding

Ornamental kale (Brassica oleracea var. acephala) is an attractive ornamental plant with a range of leaf colors and shapes. Breeding new varieties of ornamental kale has proven challenging due to its lengthy breeding cycle and the limited availability of genetic markers. In this study, a F1DH ornamental kale population comprising 300 DH lines was constructed using microspore culture. A high-density genetic map was developed by conducting whole-genome sequencing on 150 individuals from the F1DH population. The genetic map contained 1696 bin markers with 982,642 single-nucleotide polymorphisms (SNPs) spanning a total distance of 775.81 cM on all nine chromosomes with an average distance between markers of 0.46 cM. The ornamental kale genetic map contained substantially more SNP markers compared with published genetic maps for other B. oleracea crops. Furthermore, utilizing this high-density genetic map, we identified seven quantitative trait loci (QTLs) that significantly influence the leaf shape of ornamental kale. These findings are valuable for understanding the genetic basis of key agronomic traits in ornamental kale. The F1DH progenies provide an excellent resource for germplasm innovation and breeding new varieties of ornamental kale. Additionally, the high-density genetic map provides crucial insights for gene mapping and unraveling the molecular mechanisms behind important agronomic traits in ornamental kale.

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

•

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

•

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

•

Nanocarbon solution obviously reduces 4 indolic glucosinolate productions.

•

Nanocarbon has great great impact on glucosinolate biosynthesis and pathway.

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

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

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

Sustainable Use of Bioactive Compounds from Solanum Tuberosum and Brassicaceae Wastes and by-Products for Crop Protection-A Review

Defatted seed meals of oleaginous Brassicaceae, such as Eruca sativa, and potato peel are excellent plant matrices to recover potentially useful biomolecules from industrial processes in a circular strategy perspective aiming at crop protection. These biomolecules, mainly glycoalkaloids and phenols for potato and glucosinolates for Brassicaceae, have been proven to be effective against microbes, fungi, nematodes, insects, and even parasitic plants. Their role in plant protection is overviewed, together with the molecular basis of their synthesis in plant, and the description of their mechanisms of action. Possible genetic and biotechnological strategies are presented to increase their content in plants. Genetic mapping and identification of closely linked molecular markers are useful to identify the loci/genes responsible for their accumulation and transfer them to elite cultivars in breeding programs. Biotechnological approaches can be used to modify their allelic sequence and enhance the accumulation of the bioactive compounds. How the global challenges, such as reducing agri-food waste and increasing sustainability and food safety, could be addressed through bioprotector applications are discussed here.

Health-promoting phytochemicals and antioxidant capacity in different organs from six varieties of Chinese kale

Chinese kale (Brassica oleracea var. alboglabra) has high nutritional value. This study investigated the contents of glucosinolates, antioxidants (chlorophylls, carotenoids, vitamin C, and total phenolics), and antioxidant capacity in five organs from six varieties of Chinese kale. The highest concentrations of individual and total glucosinolates were in the roots and inflorescences, respectively. The highest levels of antioxidants and antioxidant capacity were in inflorescences and leaves. Plant organs played a predominant role in glucosinolate and antioxidant accumulation. Glucoiberin, glucoraphanin, and glucobrassicin, the main anticarcinogenic glucosinolates, could be enhanced simultaneously because of their high positive correlations. The relationship between glucosinolates and antioxidant capacity indicated that glucobrassicin might contribute to the total antioxidant capacity. These results provide useful information related to consumption, breeding of functional varieties, and use of the non-edible organs of Chinese kale.

Next-Gen Approaches to Flavor-Related Metabolism

Although flavor is an essential element for consumer acceptance of food, breeding programs have focused primarily on yield, leading to significant declines in flavor for many vegetables. The deterioration of flavor quality has concerned breeders; however, the complexity of this trait has hindered efforts to improve or even maintain it. Recently, the integration of flavor-associated metabolic profiling with other omics methodologies derived from big data has become a prominent trend in this research field. Here, we provide an overview of known metabolites contributing to flavor in the major vegetables as well as genetic analyses of the relevant metabolic pathways based on different approaches, especially multi-omics. We present examples demonstrating how omics analyses can help us to understand the accomplishments of historical flavor breeding practices and implement further improvements. The integration of genetics, cultivation, and postharvest practices with genome-scale data analyses will create enormous potential for further flavor quality improvements.

Rocket science: A review of phytochemical & health-related research in Eruca & Diplotaxis species

•

Recent phytochemical research in rocket species is critically reviewed.

•

Glucosinolates and hydrolysis products change over growth and shelf life.

•

Experiments should better consider and account for commercial practices.

•

Research should be focused on providing benefits to the end consumer.

Rocket species (Eruca spp. and Diplotaxis spp.) are becoming increasingly important leafy salad crops across the world. Numerous scientific research papers have been published in recent years surrounding the potential health benefits associated with phytochemicals contained in leaves, such as flavonoids and isothiocyanates. Other research of note has also been conducted into the unique taste and flavour properties of leaves, which can be hot, peppery, bitter, and sweet depending upon the genotype and phytochemical composition. While research into these aspects is increasing, some studies lack cohesion and in-depth knowledge of commercial breeding and cultivation practices that makes interpretation and application of results difficult. This review draws together all significant research findings in these crops over the last five years, and highlights areas that require further exploration and/or resolution. We also advise on experimental considerations for these species to allow for more meaningful utilisation of findings in the commercial sector.

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

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

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

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

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

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

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.

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.

Construction of chromosome segment substitution lines enables QTL mapping for flowering and morphological traits in Brassica rapa

Chromosome segment substitution lines (CSSLs) represent a powerful method for precise quantitative trait loci (QTL) detection of complex agronomical traits in plants. In this study, we used a marker-assisted backcrossing strategy to develop a population consisting of 63 CSSLs, derived from backcrossing of the F1 generated from a cross between two Brassica rapa subspecies: "Chiifu" (ssp. pekinensis), the Brassica "A" genome-represented line used as the donor, and "49caixin" (ssp. parachinensis), a non-heading cultivar used as the recipient. The 63 CSSLs covered 87.95% of the B. rapa genome. Among them, 39 lines carried a single segment; 15 lines, two segments; and nine lines, three or more segments of the donor parent chromosomes. To verify the potential advantage of these CSSL lines, we used them to locate QTL for six morphology-related traits. A total of 58 QTL were located on eight chromosomes for all six traits: 17 for flowering time, 14 each for bolting time and plant height, six for plant diameter, two for leaf width, and five for flowering stalk diameter. Co-localized QTL were mainly distributed on eight genomic regions in A01, A02, A05, A06, A08, A09, and A10, present in the corresponding CSSLs. Moreover, new chromosomal fragments that harbored QTL were identified using the findings of previous studies. The CSSL population constructed in our study paves the way for fine mapping and cloning of candidate genes involved in late bolting, flowering, and plant architecture-related traits in B. rapa. Furthermore, it has great potential for future marker-aided gene/QTL pyramiding of other interesting traits in B. rapa breeding.

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.

Glucosinolate metabolism, functionality and breeding for the improvement of Brassicaceae vegetables

Unique secondary metabolites, glucosinolates (S-glucopyranosyl thiohydroximates), are naturally occurring S-linked glucosides found mainly in Brassicaceae plants. They are enzymatically hydrolyzed to produce sulfate ions, D-glucose, and characteristic degradation products such as isothiocyanates. The functions of glucosinolates in the plants remain unclear, but isothiocyanates possessing a pungent or irritating taste and odor might be associated with plant defense from microbes. Isothiocyanates have been studied extensively in experimental in vitro and in vivo carcinogenesis models for their cancer chemopreventive properties. The beneficial isothiocyanates, glucosinolates that are functional for supporting human health, have received attention from many scientists studying plant breeding, plant physiology, plant genetics, and food functionality. This review presents a summary of recent topics related with glucosinolates in the Brassica family, along with a summary of the chemicals, metabolism, and genes of glucosinolates in Brassicaceae. The bioavailabilities of isothiocyanates from certain functional glucosinolates and the importance of breeding will be described with emphasis on glucosinolates.

Comparative genomics of Brassicaceae crops

The family Brassicaceae is one of the major groups of the plant kingdom and comprises diverse species of great economic, agronomic and scientific importance, including the model plant Arabidopsis. The sequencing of the Arabidopsis genome has revolutionized our knowledge in the field of plant biology and provides a foundation in genomics and comparative biology. Genomic resources have been utilized in Brassica for diversity analyses, construction of genetic maps and identification of agronomic traits. In Brassicaceae, comparative sequence analysis across the species has been utilized to understand genome structure, evolution and the detection of conserved genomic segments. In this review, we focus on the progress made in genetic resource development, genome sequencing and comparative mapping in Brassica and related species. The utilization of genomic resources and next-generation sequencing approaches in improvement of Brassica crops is also discussed.

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.

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.

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.

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.

Quantitative trait loci mapping in Brassica rapa revealed the structural and functional conservation of genetic loci governing morphological and yield component traits in the A, B, and C subgenomes of Brassica species

Brassica rapa is an important crop species that produces vegetables, oilseed, and fodder. Although many studies reported quantitative trait loci (QTL) mapping, the genes governing most of its economically important traits are still unknown. In this study, we report QTL mapping for morphological and yield component traits in B. rapa and comparative map alignment between B. rapa, B. napus, B. juncea, and Arabidopsis thaliana to identify candidate genes and conserved QTL blocks between them. A total of 95 QTL were identified in different crucifer blocks of the B. rapa genome. Through synteny analysis with A. thaliana, B. rapa candidate genes and intronic and exonic single nucleotide polymorphisms in the parental lines were detected from whole genome resequenced data, a few of which were validated by mapping them to the QTL regions. Semi-quantitative reverse transcriptase PCR analysis showed differences in the expression levels of a few genes in parental lines. Comparative mapping identified five key major evolutionarily conserved crucifer blocks (R, J, F, E, and W) harbouring QTL for morphological and yield components traits between the A, B, and C subgenomes of B. rapa, B. juncea, and B. napus. The information of the identified candidate genes could be used for breeding B. rapa and other related Brassica species.

Genetic analysis of metabolome-phenotype interactions: from model to crop species

The past decade has seen increased interest from the scientific community, and particularly plant biologists, in integrating metabolic approaches into research aimed at unraveling phenotypic diversity and its underlying genetic variation. Advances in plant metabolomics have enabled large-scale analyses that have identified qualitative and quantitative variation in the metabolic content of various species, and this variation has been linked to genetic factors through genetic-mapping approaches, providing a glimpse of the genetic architecture of the plant metabolome. Parallel analyses of morphological phenotypes and physiological performance characteristics have further enhanced our understanding of the complex molecular mechanisms regulating these quantitative traits. This review aims to illustrate the advantages of including assessments of phenotypic and metabolic diversity in investigations of the genetic basis of complex traits, and the value of this approach in studying agriculturally important crops. We highlight the ground-breaking work on model species and discuss recent achievements in important crop species.

Reducing progoitrin and enriching glucoraphanin in Brassica napus seeds through silencing of the GSL-ALK gene family

The hydrolytic products of glucosinolates in brassica crops are bioactive compounds. Some glucosinolate derivatives such as oxazolidine-2-thione from progoitrin in brassica oilseed meal are toxic and detrimental to animals, but some isothiocyanates such as sulforaphane are potent anti-carcinogens that have preventive effects on several human cancers. In most B. rapa, B. napus and B. juncea vegetables and oilseeds, there is no or only trace amount of glucoraphanin that is the precursor to sulforaphane. In this paper, RNA interference (RNAi) of the GSL-ALK gene family was used to down-regulate the expression of GSL-ALK genes in B. napus. The detrimental glucosinolate progoitrin was reduced by 65 %, and the beneficial glucosinolate glucoraphanin was increased to a relatively high concentration (42.6 μmol g(-1) seed) in seeds of B. napus transgenic plants through silencing of the GSL-ALK gene family. Therefore, there is potential application of the new germplasm with reduced detrimental glucosinolates and increased beneficial glucosinolates for producing improved brassica vegetables.

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

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

Glucosinolate biosynthetic genes in Brassica rapa

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

Genic microsatellite markers in Brassica rapa: development, characterization, mapping, and their utility in other cultivated and wild Brassica relatives

Genic microsatellite markers, also known as functional markers, are preferred over anonymous markers as they reveal the variation in transcribed genes among individuals. In this study, we developed a total of 707 expressed sequence tag-derived simple sequence repeat markers (EST-SSRs) and used for development of a high-density integrated map using four individual mapping populations of B. rapa. This map contains a total of 1426 markers, consisting of 306 EST-SSRs, 153 intron polymorphic markers, 395 bacterial artificial chromosome-derived SSRs (BAC-SSRs), and 572 public SSRs and other markers covering a total distance of 1245.9 cM of the B. rapa genome. Analysis of allelic diversity in 24 B. rapa germplasm using 234 mapped EST-SSR markers showed amplification of 2 alleles by majority of EST-SSRs, although amplification of alleles ranging from 2 to 8 was found. Transferability analysis of 167 EST-SSRs in 35 species belonging to cultivated and wild brassica relatives showed 42.51% (Sysimprium leteum) to 100% (B. carinata, B. juncea, and B. napus) amplification. Our newly developed EST-SSRs and high-density linkage map based on highly transferable genic markers would facilitate the molecular mapping of quantitative trait loci and the positional cloning of specific genes, in addition to marker-assisted selection and comparative genomic studies of B. rapa with other related species.

Design of a Brassica rapa core collection for association mapping studies

A Brassica rapa collection of 239 accessions, based on two core collections representing different morphotypes from different geographical origins, is presented and its use for association mapping is illustrated for flowering time. We analyzed phenotypic variation of leaf and seed pod traits, plant architecture, and flowering time using data collected from three field experiments and evaluated the genetic diversity with a set of SSR markers. The Wageningen University and Research Centre (WUR) and the Vavilov Research Institute of Plant Industry (VIR) core collections had similar representations of most morphotypes, as illustrated by the phenotypic and genetic variation within these groups. The analysis of population structure revealed five subgroups in the collection, whereas previous studies of the WUR core collection indicated four subgroups; the fifth group identified consisted mainly of oil accessions from the VIR core collection, winter oils from Pakistan, and a number of other types. A very small group of summer oils is described, that is not related to other oil accessions. A candidate gene approach was chosen for association mapping of flowering time with a BrFLC1 biallelic CAPS marker and a BrFLC2 multiallelic SSR marker. The two markers were significantly associated with flowering time, but their effects were confined to certain morphotypes and (or) alleles. Based on these results, we discuss the optimal design for an association mapping population and the need to fix the heterogeneous accessions to facilitate phenotyping and genotyping.

Comparative methods for association studies: a case study on metabolite variation in a Brassica rapa core collection

Background

Association mapping is a statistical approach combining phenotypic traits and genetic diversity in natural populations with the goal of correlating the variation present at phenotypic and allelic levels. It is essential to separate the true effect of genetic variation from other confounding factors, such as adaptation to different uses and geographical locations. The rapid availability of large datasets makes it necessary to explore statistical methods that can be computationally less intensive and more flexible for data exploration.

Methodology/Principal Findings

A core collection of 168 Brassica rapa accessions of different morphotypes and origins was explored to find genetic association between markers and metabolites: tocopherols, carotenoids, chlorophylls and folate. A widely used linear model with modifications to account for population structure and kinship was followed for association mapping. In addition, a machine learning algorithm called Random Forest (RF) was used as a comparison. Comparison of results across methods resulted in the selection of a set of significant markers as promising candidates for further work. This set of markers associated to the metabolites can potentially be applied for the selection of genotypes with elevated levels of these metabolites.

Conclusions/Significance

The incorporation of the kinship correction into the association model did not reduce the number of significantly associated markers. However incorporation of the STRUCTURE correction (Q matrix) in the linear regression model greatly reduced the number of significantly associated markers. Additionally, our results demonstrate that RF is an interesting complementary method with added value in association studies in plants, which is illustrated by the overlap in markers identified using RF and a linear mixed model with correction for kinship and population structure. Several markers that were selected in RF and in the models with correction for kinship, but not for population structure, were also identified as QTLs in two bi-parental DH populations.

Phytochemicals of Brassicaceae in plant protection and human health--influences of climate, environment and agronomic practice

In this review, we provide an overview of the role of glucosinolates and other phytochemical compounds present in the Brassicaceae in relation to plant protection and human health. Current knowledge of the factors that influence phytochemical content and profile in the Brassicaceae is also summarized and multi-factorial approaches are briefly discussed. Variation in agronomic conditions (plant species, cultivar, developmental stage, plant organ, plant competition, fertilization, pH), season, climatic factors, water availability, light (intensity, quality, duration) and CO(2) are known to significantly affect content and profile of phytochemicals. Phytochemicals such as the glucosinolates and leaf surface waxes play an important role in interactions with pests and pathogens. Factors that affect production of phytochemicals are important when designing plant protection strategies that exploit these compounds to minimize crop damage caused by plant pests and pathogens. Brassicaceous plants are consumed increasingly for possible health benefits, for example, glucosinolate-derived effects on degenerative diseases such as cancer, cardiovascular and neurodegenerative diseases. Thus, factors influencing phytochemical content and profile in the production of brassicaceous plants are worth considering both for plant and human health. Even though it is known that factors that influence phytochemical content and profile may interact, studies of plant compounds were, until recently, restricted by methods allowing only a reductionistic approach. It is now possible to design multi-factorial experiments that simulate their combined effects. This will provide important information to ecologists, plant breeders and agronomists.

Ecological genetics and genomics of plant defenses: Evidence and approaches

Herbivores exert significant selection on plants, and plants have evolved a variety of constitutive and inducible defenses to resist and tolerate herbivory. Assessing the genetic mechanisms that influence defenses against herbivores will deepen our understanding of the evolution of essential phenotypic traits.Ecogenomics is a powerful interdisciplinary approach that can address fundamental questions about the ecology and evolutionary biology of species, such as: which evolutionary forces maintain variation within a population? and What is the genetic architecture of adaptation? This field seeks to identify gene regions that influence ecologically-important traits, assess the fitness consequences under natural conditions of alleles at key quantitative trait loci (QTLs), and test how the abiotic and biotic environment affects gene expression.Here, we review ecogenomics techniques and emphasize how this framework can address long-standing and emerging questions relating to anti-herbivore defenses in plants. For example, ecogenomics tools can be used to investigate: inducible vs. constitutive defenses; tradeoffs between resistance and tolerance; adaptation to the local herbivore community; selection on alleles that confer resistance and tolerance in natural populations; and whether different genes are activated in response to specialist vs. generalist herbivores and to different types of damage.Ecogenomic studies can be conducted with model species, such as Arabidopsis, or their relatives, in which case myriad molecular tools are already available. Burgeoning sequence data will also facilitate ecogenomic studies of non-model species. Throughout this paper, we highlight approaches that are particularly suitable for ecological studies of non-model organisms, discuss the benefits and disadvantages of specific techniques, and review bioinformatic tools for analyzing data.We focus on established and promising techniques, such as QTL mapping with pedigreed populations, genome wide association studies, transcription profiling strategies, population genomics, and transgenic methodologies. Many of these techniques are complementary and can be used jointly to investigate the genetic architecture of defense traits and selection on alleles in nature.

BrFLC2 (FLOWERING LOCUS C) as a candidate gene for a vernalization response QTL in Brassica rapa

Flowering time is an important agronomic trait, and wide variation exists among Brassica rapa. In Arabidopsis, FLOWERING LOCUS C (FLC) plays an important role in modulating flowering time and the response to vernalization. Brassica rapa contains several paralogues of FLC at syntenic regions. BrFLC2 maps under a major flowering time and vernalization response quantitative trait locus (QTL) at the top of A02. Here the effects of vernalization on flowering time in a double haploid (DH) population and on BrFLC2 expression in selected lines of a DH population in B. rapa are descibed. The effect of the major flowering time QTL on the top of A02 where BrFLC2 maps clearly decreases upon vernalization, which points to a role for BrFLC2 underlying the QTL. In all developmental stages and tissues (seedlings, cotyledons, and leaves), BrFLC2 transcript levels are higher in late flowering pools of DH lines than in pools of early flowering DH lines. BrFLC2 expression diminished after different durations of seedling vernalization in both early and late DH lines. The reduction of BrFLC2 expression upon seedling vernalization of both early and late flowering DH lines was strongest at the seedling stage and diminished in subsequent growth stages, which suggests that the commitment to flowering is already set at very early developmental stages. Taken together, these data support the hypothesis that BrFLC2 is a candidate gene for the flowering time and vernalization response QTL in B. rapa.

GERMINATOR: a software package for high-throughput scoring and curve fitting of Arabidopsis seed germination

Over the past few decades seed physiology research has contributed to many important scientific discoveries and has provided valuable tools for the production of high quality seeds. An important instrument for this type of research is the accurate quantification of germination; however gathering cumulative germination data is a very laborious task that is often prohibitive to the execution of large experiments. In this paper we present the germinator package: a simple, highly cost-efficient and flexible procedure for high-throughput automatic scoring and evaluation of germination that can be implemented without the use of complex robotics. The germinator package contains three modules: (i) design of experimental setup with various options to replicate and randomize samples; (ii) automatic scoring of germination based on the color contrast between the protruding radicle and seed coat on a single image; and (iii) curve fitting of cumulative germination data and the extraction, recap and visualization of the various germination parameters. The curve-fitting module enables analysis of general cumulative germination data and can be used for all plant species. We show that the automatic scoring system works for Arabidopsis thaliana and Brassica spp. seeds, but is likely to be applicable to other species, as well. In this paper we show the accuracy, reproducibility and flexibility of the germinator package. We have successfully applied it to evaluate natural variation for salt tolerance in a large population of recombinant inbred lines and were able to identify several quantitative trait loci for salt tolerance. Germinator is a low-cost package that allows the monitoring of several thousands of germination tests, several times a day by a single person.

A Brassica rapa linkage map of EST-based SNP markers for identification of candidate genes controlling flowering time and leaf morphological traits

For identification of genes responsible for varietal differences in flowering time and leaf morphological traits, we constructed a linkage map of Brassica rapa DNA markers including 170 EST-based markers, 12 SSR markers, and 59 BAC sequence-based markers, of which 151 are single nucleotide polymorphism (SNP) markers. By BLASTN, 223 markers were shown to have homologous regions in Arabidopsis thaliana, and these homologous loci covered nearly the whole genome of A. thaliana. Synteny analysis between B. rapa and A. thaliana revealed 33 large syntenic regions. Three quantitative trait loci (QTLs) for flowering time were detected. BrFLC1 and BrFLC2 were linked to the QTLs for bolting time, budding time, and flowering time. Three SNPs in the promoter, which may be the cause of low expression of BrFLC2 in the early-flowering parental line, were identified. For leaf lobe depth and leaf hairiness, one major QTL corresponding to a syntenic region containing GIBBERELLIN 20 OXIDASE 3 and one major QTL containing BrGL1, respectively, were detected. Analysis of nucleotide sequences and expression of these genes suggested possible involvement of these genes in leaf morphological traits.

Genome-wide identification of glucosinolate synthesis genes in Brassica rapa

Glucosinolates play important roles in plant defense against herbivores and microbes, as well as in human nutrition. Some glucosinolate-derived isothiocyanate and nitrile compounds have been clinically proven for their anticarcinogenic activity. To better understand glucosinolate biosynthesis in Brassica rapa, we conducted a comparative genomics study with Arabidopsis thaliana and identified total 56 putative biosynthetic and regulator genes. This established a high colinearity in the glucosinolate biosynthesis pathway between Arabidopsis and B. rapa. Glucosinolate genes in B. rapa share 72-94% nucleotide sequence identity with the Arabidopsis orthologs and exist in different copy numbers. The exon/intron split pattern of B. rapa is almost identical to that of Arabidopsis, although inversion, insertion, deletion and intron size variations commonly occur. Four genes appear to be nonfunctional as a result of the presence of a frame shift mutation and retrotransposon insertion. At least 12 paralogs of desulfoglucosinolate sulfotransferase were found in B. rapa, whereas only three were found in Arabidopsis. The expression of those paralogs was not tissue-specific but varied greatly depending on B. rapa tissue types. Expression was also developmentally regulated in some paralogs but not in other paralogs. Most of the regulator genes are present as triple copies. Accordingly, glucosinolate synthesis and regulation in B. rapa appears to be more complex than that of Arabidopsis. With the isolation and further characterization of the endogenous genes, health-beneficial vegetables or desirable animal feed crops could be developed by metabolically engineering the glucosinolate pathway.