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

A Comprehensive Gene Inventory for Glucosinolate Biosynthetic Pathway in Arabidopsis thaliana

Glucosinolates (GSLs) are plant secondary metabolites comprising sulfur and nitrogen mainly found in plants from the order of Brassicales, such as broccoli, cabbage, and Arabidopsis thaliana. The activated forms of GSL play important roles in fighting against pathogens and have health benefits to humans. The increasing amount of data on A. thaliana generated from various omics technologies can be investigated more deeply in search of new genes or compounds involved in GSL biosynthesis and metabolism. This review describes a comprehensive inventory of A. thaliana GSLs identified from published literature and databases such as KNApSAcK, KEGG, and AraCyc. A total of 113 GSL genes encoding for 23 transcription components, 85 enzymes, and five protein transporters were experimentally characterized in the past two decades. Continuous efforts are still on going to identify all molecules related to the production of GSLs. A manually curated database known as SuCCombase (http://plant-scc.org) was developed to serve as a comprehensive GSL inventory. Realizing lack of information on the regulation of GSL biosynthesis and degradation mechanisms, this review also includes relevant information and their connections with crosstalk among various factors, such as light, sulfur metabolism, and nitrogen metabolism, not only in A. thaliana but also in other crucifers.

Comparative transcriptomics revealed differential regulation of defense related genes in Brassica juncea leading to successful and unsuccessful infestation by aphid species

Productivity of Indian mustard (B. juncea), a major oil yielding crop in rapeseed-mustard group is heavily inflicted by mustard aphid, L. erysimi. Mustard aphid, a specialist aphid species on rapeseed-mustard crops, rapidly multiplies and colonizes the plants leading to successful infestation. In contrary, legume specific cowpea aphid, A. craccivora when released on B. juncea plants fails to build up population and thus remains unsuccessful in infestation. In the present study, differential host response of B. juncea to the two aphid species, one being successful insect-pest and the other being unsuccessful on it has been studied based on transcriptome analysis. Differential feeding efficiency of the two aphid species on mustard plants was evident from the amount of secreted honeydews. Leaf-transcriptomes of healthy and infested plants, treated with the two aphid species, were generated by RNA sequencing on Illumina platform and de novo assembly of the quality reads. A comparative assessment of the differentially expressed genes due to treatments revealed a large extent of overlaps as well as distinctness with respect to the set of genes and their direction of regulation. With respect to host-genes related to transcription factors, oxidative homeostasis, defense hormones and secondary metabolites, L. erysimi led to either suppression or limited activation of the transcript levels compared to A. craccivora. Further, a comprehensive view of the DEGs suggested more potential of successful insect-pests towards transcriptional reprogramming of the host. qRT-PCR based validation of randomly selected up- and down-regulated transcripts authenticated the transcriptome data.

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

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

Drought-Tolerant Brassica rapa Shows Rapid Expression of Gene Networks for General Stress Responses and Programmed Cell Death Under Simulated Drought Stress

Production of oilseed rape/canola (Brassica napus) is increasingly threatened by dry conditions while the demand for vegetable oil is increasing. Brassica rapa is a genetically diverse ancestor of B. napus, and is readily crossed with B. napus. Recently, we reported promising levels of drought tolerance in a wild type of B. rapa which could be a source of drought tolerance for B. napus. We analysed global gene expression by messenger RNA sequencing in seedlings of the drought-tolerant and a drought-sensitive genotype of B. rapa under simulated drought stress and control conditions. A subset of stress-response genes were validated by reverse transcription quantitative PCR. Gene ontology enrichment analysis and pathway enrichment analysis revealed major differences between the two genotypes in the mode and onset of stress responses in the first 12 h of treatment. Drought-tolerant plants reacted uniquely and rapidly by upregulating genes associated with jasmonic acid and salicylic acid metabolism, as well as genes known to cause endoplasmic reticulum stress and induction of programmed cell death. Conversely, active responses in drought-sensitive plants were delayed until 8 or 12 h after stress application. The results may help to identify biomarkers for selection of breeding materials with potentially improved drought tolerance.

Genotypes, Networks, Phenotypes: Moving Toward Plant Systems Genetics

One of the central goals in biology is to understand how and how much of the phenotype of an organism is encoded in its genome. Although many genes that are crucial for organismal processes have been identified, much less is known about the genetic bases underlying quantitative phenotypic differences in natural populations. We discuss the fundamental gap between the large body of knowledge generated over the past decades by experimental genetics in the laboratory and what is needed to understand the genotype-to-phenotype problem on a broader scale. We argue that systems genetics, a combination of systems biology and the study of natural variation using quantitative genetics, will help to address this problem. We present major advances in these two mostly disconnected areas that have increased our understanding of the developmental processes of flowering time control and root growth. We conclude by illustrating and discussing the efforts that have been made toward systems genetics specifically in plants.

Genome-Wide Identification and Characterization of bZIP Transcription Factors in Brassica oleracea under Cold Stress

Cabbages (Brassica oleracea L.) are an important vegetable crop around world, and cold temperature is among the most significant abiotic stresses causing agricultural losses, especially in cabbage crops. Plant bZIP transcription factors play diverse roles in biotic/abiotic stress responses. In this study, 119 putative BolbZIP transcription factors were identified using amino acid sequences from several bZIP domain consensus sequences. The BolbZIP members were classified into 63 categories based on amino acid sequence similarity and were also compared with BrbZIP and AtbZIP transcription factors. Based on this BolbZIP identification and classification, cold stress-responsive BolbZIP genes were screened in inbred lines, BN106 and BN107, using RNA sequencing data and qRT-PCR. The expression level of the 3 genes, Bol008071, Bol033132, and Bol042729, was significantly increased in BN107 under cold conditions and was unchanged in BN106. The upregulation of these genes in BN107, a cold-susceptible inbred line, suggests that they might be significant components in the cold response. Among three identified genes, Bol033132 has 97% sequence similarity to Bra020735, which was identified in a screen for cold-related genes in B. rapa and a protein containing N-rich regions in LCRs. The results obtained in this study provide valuable information for understanding the potential function of BolbZIP transcription factors in cold stress responses.

Selection during crop diversification involves correlated evolution of the circadian clock and ecophysiological traits in Brassica rapa

Crop selection often leads to dramatic morphological diversification, in which allocation to the harvestable component increases. Shifts in allocation are predicted to impact (as well as rely on) physiological traits; yet, little is known about the evolution of gas exchange and related anatomical features during crop diversification. In Brassica rapa, we tested for physiological differentiation among three crop morphotypes (leaf, turnip, and oilseed) and for correlated evolution of circadian, gas exchange, and phenological traits. We also examined internal and surficial leaf anatomical features and biochemical limits to photosynthesis. Crop types differed in gas exchange; oilseed varieties had higher net carbon assimilation and stomatal conductance relative to vegetable types. Phylogenetically independent contrasts indicated correlated evolution between circadian traits and both gas exchange and biomass accumulation; shifts to shorter circadian period (closer to 24 h) between phylogenetic nodes are associated with higher stomatal conductance, lower photosynthetic rate (when CO2 supply is factored out), and lower biomass accumulation. Crop type differences in gas exchange are also associated with stomatal density, epidermal thickness, numbers of palisade layers, and biochemical limits to photosynthesis. Brassica crop diversification involves correlated evolution of circadian and physiological traits, which is potentially relevant to understanding mechanistic targets for crop improvement.

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

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

Molecular Mapping and QTL for Expression Profiles of Flavonoid Genes in Brassica napus

Flavonoids are secondary metabolites that are extensively distributed in the plant kingdom and contribute to seed coat color formation in rapeseed. To decipher the genetic networks underlying flavonoid biosynthesis in rapeseed, we constructed a high-density genetic linkage map with 1089 polymorphic loci (including 464 SSR loci, 97 RAPD loci, 451 SRAP loci, and 75 IBP loci) using recombinant inbred lines (RILs). The map consists of 19 linkage groups and covers 2775 cM of the B. napus genome with an average distance of 2.54 cM between adjacent markers. We then performed expression quantitative trait locus (eQTL) analysis to detect transcript-level variation of 18 flavonoid biosynthesis pathway genes in the seeds of the 94 RILs. In total, 72 eQTLs were detected and found to be distributed among 15 different linkage groups that account for 4.11% to 52.70% of the phenotypic variance atrributed to each eQTL. Using a genetical genomics approach, four eQTL hotspots together harboring 28 eQTLs associated with 18 genes were found on chromosomes A03, A09, and C08 and had high levels of synteny with genome sequences of A. thaliana and Brassica species. Associated with the trans-eQTL hotspots on chromosomes A03, A09, and C08 were 5, 17, and 1 genes encoding transcription factors, suggesting that these genes have essential roles in the flavonoid biosynthesis pathway. Importantly, bZIP25, which is expressed specifically in seeds, MYC1, which controls flavonoid biosynthesis, and the R2R3-type gene MYB51, which is involved in the synthesis of secondary metabolites, were associated with the eQTL hotspots, and these genes might thus be involved in different flavonoid biosynthesis pathways in rapeseed. Hence, further studies of the functions of these genes will provide insight into the regulatory mechanism underlying flavonoid biosynthesis, and lay the foundation for elaborating the molecular mechanism of seed coat color formation in B. napus.

A Systems Genetics Approach Identifies Gene Regulatory Networks Associated with Fatty Acid Composition in Brassica rapa Seed

Fatty acids in seeds affect seed germination and seedling vigor, and fatty acid composition determines the quality of seed oil. In this study, quantitative trait locus (QTL) mapping of fatty acid and transcript abundance was integrated with gene network analysis to unravel the genetic regulation of seed fatty acid composition in a Brassica rapa doubled haploid population from a cross between a yellow sarson oil type and a black-seeded pak choi. The distribution of major QTLs for fatty acids showed a relationship with the fatty acid types: linkage group A03 for monounsaturated fatty acids, A04 for saturated fatty acids, and A05 for polyunsaturated fatty acids. Using a genetical genomics approach, expression quantitative trait locus (eQTL) hotspots were found at major fatty acid QTLs on linkage groups A03, A04, A05, and A09. An eQTL-guided gene coexpression network of lipid metabolism-related genes showed major hubs at the genes BrPLA2-ALPHA, BrWD-40, a number of seed storage protein genes, and the transcription factor BrMD-2, suggesting essential roles for these genes in lipid metabolism. Three subnetworks were extracted for the economically important and most abundant fatty acids erucic, oleic, linoleic, and linolenic acids. Network analysis, combined with comparison of the genome positions of cis- or trans-eQTLs with fatty acid QTLs, allowed the identification of candidate genes for genetic regulation of these fatty acids. The generated insights in the genetic architecture of fatty acid composition and the underlying complex gene regulatory networks in B. rapa seeds are discussed.

Molecular Cloning, Expression Pattern and Genotypic Effects on Glucoraphanin Biosynthetic Related Genes in Chinese Kale (Brassica oleracea var. alboglabra Bailey)

Glucoraphanin is a plant secondary metabolite that is involved in plant defense and imparts health-promoting properties to cruciferous vegetables. In this study, three genes involved in glucoraphanin metabolism, branched-chain aminotransferase 4 (BCAT4), methylthioalkylmalate synthase 1 (MAM1) and dihomomethionine N-hydroxylase (CYP79F1), were cloned from Chinese kale (Brassica oleracea var. alboglabra Bailey). Sequence homology and phylogenetic analysis identified these genes and confirmed the evolutionary status of Chinese kale. The transcript levels of BCAT4, MAM1 and CYP79F1 were higher in cotyledon, leaf and stem compared with flower and silique. BCAT4, MAM1 and CYP79F1 were expressed throughout leaf development with lower transcript levels during the younger stages. Glucoraphanin content varied extensively among different varieties, which ranged from 0.25 to 2.73 µmol·g(-1) DW (dry weight). Expression levels of BCAT4 and MAM1 were high at vegetative-reproductive transition phase, while CYP79F1 was expressed high at reproductive phase. BCAT4, MAM1 and CYP79F1 were expressed significantly high in genotypes with high glucoraphanin content. All the results provided a better understanding of the roles of BCAT4, MAM1 and CYP79F1 in the glucoraphanin biosynthesis of Chinese kale.