Genome-wide identification of glucosinolate synthesis genes in Brassica rapa

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.

The fusion of genomes leads to more options: A comparative investigation on the desulfo-glucosinolate sulfotransferases of Brassica napus and homologous proteins of Arabidopsis thaliana

Sulfotransferases (SOTs) (EC 2.8.2.-) play a crucial role in the glucosinolate (Gl) biosynthesis, by catalyzing the final step of the core glucosinolate formation. In Arabidopsis thaliana the three desulfo (ds)-Gl SOTs AtSOT16, AtSOT17 and AtSOT18 were previously characterized, showing different affinities to ds-Gls. But can the knowledge about these SOTs be generally transferred to other Gl-synthesizing plants? It was investigated how many SOTs are present in the economically relevant crop plant Brassica napus L., and if it is possible to predict their characteristics by sequence analysis. The recently sequenced B. napus is a hybrid of Brassica rapa and Brassica oleracea. By database research, 71 putative functional BnSOT family members were identified and at least eleven of those are putative ds-Gl SOTs. Besides the homologs of AtSOT16 - 18, phylogenetic analyses revealed new subfamilies of ds-Gl SOTs, which are not present in A. thaliana. Three of the B. napus ds-Gl SOT proteins were expressed and purified, and characterized by determining the substrate affinities to different ds-Gls. Two of them, BnSOT16-a and BnSOT16-b, showed a significantly higher affinity to an indolic ds-Gl, similarly to AtSOT16. Additionally, BnSOT17-a was characterized and showed a higher affinity to long chained aliphatic Gls, similarly to AtSOT17. Identification of homologs to AtSOT18 was less reliable, because putative SOT18 sequences are more heterogeneous and confirmation of similar characteristics was not possible.

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.

Calcium mitigates the stress caused by ZnSO4 as a sulphur fertilizer and enhances the sulforaphane formation of broccoli sprouts

CaCl₂ improved the growth and sulforaphane formation of broccoli sprouts under ZnSO₄via enhancing key bioactive substances level, antioxidant capacity, myrosinase activity and related genes expression.

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.

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.

Functional identification of genes responsible for the biosynthesis of 1-methoxy-indol-3-ylmethyl-glucosinolate in Brassica rapa ssp. chinensis

BACKGROUND

Brassica vegetables contain a class of secondary metabolites, the glucosinolates (GS), whose specific degradation products determine the characteristic flavor and smell. While some of the respective degradation products of particular GS are recognized as health promoting substances for humans, recent studies also show evidence that namely the 1-methoxy-indol-3-ylmethyl GS might be deleterious by forming characteristic DNA adducts. Therefore, a deeper knowledge of aspects involved in the biosynthesis of indole GS is crucial to design vegetables with an improved secondary metabolite profile.

RESULTS

Initially the leafy Brassica vegetable pak choi (Brassica rapa ssp. chinensis) was established as suitable tool to elicit very high concentrations of 1-methoxy-indol-3-ylmethyl GS by application of methyl jasmonate. Differentially expressed candidate genes were discovered in a comparative microarray analysis using the 2 × 104 K format Brassica Array and compared to available gene expression data from the Arabidopsis AtGenExpress effort. Arabidopsis knock out mutants of the respective candidate gene homologs were subjected to a comprehensive examination of their GS profiles and confirmed the exclusive involvement of polypeptide 4 of the cytochrome P450 monooxygenase subfamily CYP81F in 1-methoxy-indol-3-ylmethyl GS biosynthesis. Functional characterization of the two identified isoforms coding for CYP81F4 in the Brassica rapa genome was performed using expression analysis and heterologous complementation of the respective Arabidopsis mutant.

CONCLUSIONS

Specific differences discovered in a comparative microarray and glucosinolate profiling analysis enables the functional attribution of Brassica rapa ssp. chinensis genes coding for polypeptide 4 of the cytochrome P450 monooxygenase subfamily CYP81F to their metabolic role in indole glucosinolate biosynthesis. These new identified Brassica genes will enable the development of genetic tools for breeding vegetables with improved GS composition in the near future.

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.

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.

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.

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.

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.

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.

Induced production of 1-methoxy-indol-3-ylmethyl glucosinolate by jasmonic acid and methyl jasmonate in sprouts and leaves of pak choi (Brassica rapa ssp. chinensis)

Pak choi plants (Brassica rapa ssp. chinensis) were treated with different signaling molecules methyl jasmonate, jasmonic acid, linolenic acid, and methyl salicylate and were analyzed for specific changes in their glucosinolate profile. Glucosinolate levels were quantified using HPLC-DAD-UV, with focus on induction of indole glucosinolates and special emphasis on 1-methoxy-indol-3-ylmethyl glucosinolate. Furthermore, the effects of the different signaling molecules on indole glucosinolate accumulation were analyzed on the level of gene expression using semi-quantitative realtime RT-PCR of selected genes. The treatments with signaling molecules were performed on sprouts and mature leaves to determine ontogenetic differences in glucosinolate accumulation and related gene expression. The highest increase of indole glucosinolate levels, with considerable enhancement of the 1-methoxy-indol-3-ylmethyl glucosinolate content, was achieved with treatments of sprouts and mature leaves with methyl jasmonate and jasmonic acid. This increase was accompanied by increased expression of genes putatively involved in the indole glucosinolate biosynthetic pathway. The high levels of indole glucosinolates enabled the plant to preferentially produce the respective breakdown products after tissue damage. Thus, pak choi plants treated with methyl jasmonate or jasmonic acid, are a valuable tool to analyze the specific protection functions of 1-methoxy-indole-3-carbinole in the plants defense strategy in the future.

De novo assembly and characterization of Camelina sativa transcriptome by paired-end sequencing

Background

Biofuels extracted from the seeds of Camelina sativa have recently been used successfully as environmentally friendly jet-fuel to reduce greenhouse gas emissions. Camelina sativa is genetically very close to Arabidopsis thaliana, and both are members of the Brassicaceae. Although public databases are currently available for some members of the Brassicaceae, such as A. thaliana, A. lyrata, Brassica napus, B. juncea and B. rapa, there are no public Expressed Sequence Tags (EST) or genomic data for Camelina sativa. In this study, a high-throughput, large-scale RNA sequencing (RNA-seq) of the Camelina sativa transcriptome was carried out to generate a database that will be useful for further functional analyses.

Results

Approximately 27 million clean “reads” filtered from raw reads by removal of adaptors, ambiguous reads and low-quality reads (2.42 gigabase pairs) were generated by Illumina paired-end RNA-seq technology. All of these clean reads were assembled de novo into 83,493 unigenes and 103,196 transcripts using SOAPdenovo and Trinity, respectively. The average length of the transcripts generated by Trinity was 697 bp (N50 = 976), which was longer than the average length of unigenes (319 bp, N50 = 346 bp). Nonetheless, the assembly generated by SOAPdenovo produced similar number of non-redundant hits (22,435) with that of Trinity (22,433) in BLASTN searches of the Arabidopsis thaliana CDS sequence database (TAIR). Four public databases, the Kyoto Encyclopedia of Genes and Genomes (KEGG), Swiss-prot, NCBI non-redundant protein (NR), and the Cluster of Orthologous Groups (COG), were used for unigene annotation; 67,791 of 83,493 unigenes (81.2%) were finally annotated with gene descriptions or conserved protein domains that were mapped to 25,329 non-redundant protein sequences. We mapped 27,042 of 83,493 unigenes (32.4%) to 119 KEGG metabolic pathways.

Conclusions

This is the first report of a transcriptome database for Camelina sativa, an environmentally important member of the Brassicaceae. We showed that C. savita is closely related to Arabidopsis spp. and more distantly related to Brassica spp. Although the majority of annotated genes had high sequence identity to those of A. thaliana, a substantial proportion of disease-resistance genes (NBS-encoding LRR genes) were instead more closely similar to the genes of other Brassicaceae; these genes included BrCN, BrCNL, BrNL, BrTN, BrTNL in B. rapa. As plant genomes are under long-term selection pressure from environmental stressors, conservation of these disease-resistance genes in C. sativa and B. rapa genomes implies that they are exposed to the threats from closely-related pathogens in their natural habitats.

Genotypic variation of the glucosinolate profile in pak choi (Brassica rapa ssp. chinensis)

Thirteen different pak choi (Brassica rapa ssp. chinensis) cultivars were characterized regarding their glucosinolate profile analyzed by HPLC-DAD-MS. The identified glucosinolates were subjected to principal component analysis, and three distinct groups of pak choi sprouts were identified. Group differences were marked mainly by variations in the aliphatic glucosinolate profile such as differing levels of 3-butenyl glucosinolate and 2-hydroxy-3-butenyl glucosinolate as well as by their varying proportional ratios. In addition, the three groups of pak choi sprouts varied by the presence or absence of 2-hydroxy-4-pentenyl glucosinolate and in level and composition of butyl glucosinolates. This classification is reflected by relative mRNA expression level of 2-oxoacid-dependent dioxygenase. As in sprouts, the major glucosinolates in mature leaves were found to be the aliphatic glucosinolates. However, unlike in sprouts, an additional aliphatic glucosinolate, 5-methylsulfinylpentyl glucosinolate, was detected as characteristic ontogenetic variation in mature leaves in 12 of the 13 pak choi cultivars analyzed.

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.

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

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

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.

Identification of novel soybean microRNAs involved in abiotic and biotic stresses

BACKGROUND

Small RNAs (19-24 nt) are key regulators of gene expression that guide both transcriptional and post-transcriptional silencing mechanisms in eukaryotes. Current studies have demonstrated that microRNAs (miRNAs) act in several plant pathways associated with tissue proliferation, differentiation, and development and in response to abiotic and biotic stresses. In order to identify new miRNAs in soybean and to verify those that are possibly water deficit and rust-stress regulated, eight libraries of small RNAs were constructed and submitted to Solexa sequencing.

RESULTS

The libraries were developed from drought-sensitive and tolerant seedlings and rust-susceptible and resistant soybeans with or without stressors. Sequencing the library and subsequent analyses detected 256 miRNAs. From this total, we identified 24 families of novel miRNAs that had not been reported before, six families of conserved miRNAs that exist in other plants species, and 22 families previously reported in soybean. We also observed the presence of several isomiRNAs during our analyses. To validate novel miRNAs, we performed RT-qPCR across the eight different libraries. Among the 11 miRNAs analyzed, all showed different expression profiles during biotic and abiotic stresses to soybean. The majority of miRNAs were up-regulated during water deficit stress in the sensitive plants. However, for the tolerant genotype, most of the miRNAs were down regulated. The pattern of miRNAs expression was also different for the distinct genotypes submitted to the pathogen stress. Most miRNAs were down regulated during the fungus infection in the susceptible genotype; however, in the resistant genotype, most miRNAs did not vary during rust attack. A prediction of the putative targets was carried out for conserved and novel miRNAs families.

CONCLUSIONS

Validation of our results with quantitative RT-qPCR revealed that Solexa sequencing is a powerful tool for miRNA discovery. The identification of differentially expressed plant miRNAs provides molecular evidence for the possible involvement of miRNAs in the process of water deficit- and rust-stress responses.