De novo transcriptome sequencing of radish (Raphanus sativusL.) and analysis of major genes involved in glucosinolate metabolism
- Yan Wang†1, 2, 3,
- Yan Pan†1, 2,
- Zhe Liu1, 2,
- Xianwen Zhu4,
- Lulu Zhai1, 2,
- Liang Xu1, 2,
- Rugang Yu1, 2,
- Yiqin Gong1, 2 and
- Liwang Liu1, 2Email author
© Wang et al.; licensee BioMed Central Ltd. 2013
Received: 17 June 2013
Accepted: 14 November 2013
Published: 27 November 2013
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.
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.
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.
KeywordsRadish De novo assembly RNA-Seq Transcriptome Glucosinolate metabolic pathways
Radish (Raphanus sativus L.) is an annual or biennial herb of the Brassicaceae family, and it is an economically important root vegetable crop produced throughout the world [1, 2]. The edible part of radish is its taproot, which is an excellent source of carbohydrates, dietary fiber, and essential mineral and organic nutrients to human beings [3–5]. Radish roots also contain valuable phytochemicals and have been used for many medicinal purposes [6, 7]. For example, the roots are a rich source of glucosinolates (GS) . GS and their breakdown products such as isothiocyanates (ITC) are secondary metabolites widely present in the Brassicaceae family. The ITC contribute to the flavor and taste of the Brassicaceae vegetables as an important ingredient and have anti-carcinogenic properties [8, 9].
The formation and development of taproot is a complex morphogenetic process controlled by interactions among genetic, environmental and physiological factors [1, 10–12]. Essentially, fleshy root formation is a result of selective expression of related genes. However, the lack of genomic information impedes our understanding of the molecular mechanisms underlying taproot development. Recent analysis of transcript differences between two cDNA libraries from the early and late seedling developmental stages have demonstrated that a set of genes involved in starch and sucrose metabolism, and in phenylpropanoid biosynthesis may be the dominant metabolic pathways during the early stages of taproot formation in radish . This has enabled the mining of genes that are possibly involved in taproot development. However, the molecular mechanisms involved in biosynthesis and metabolism of the related nutritional and flavor components during taproot formation are not well known, especially for many secondary metabolites such as glucosinolates.
Next-generation sequencing (NGS) -based RNA sequencing for transcriptome methods (RNA-seq) allows simultaneous acquisition of sequences for gene discovery as well as transcript identification involved in specific biological processes. This is especially suitable for non-model organisms whose genomic sequences are unknown [14–16]. In recent years, RNA-seq has emerged as a powerful method for discovering and identifying genes involved in biosynthesis of various secondary metabolites, such as, carotenoid biosynthesis in Momordica cochinchinensis, cellulose and lignin biosynthesis in Chinese fir , tea-specific compounds i.e. flavonoid, theanine and caffeine biosynthesis pathways in tea , biosynthesis of flavonoid in Safflower , biosynthesis of active ingredients in Salvia miltiorrhiza  and biosynthesis of capsaicinoid in chili pepper .
Glucosinolate content is a main trait of radish cultivars and is important for flavor formation and nutritional quality of the taproot [8, 9]. Previous studies mainly focused on developing analysis methods to determine GS content in radish, and also to determine variation in GS composition or content in different cultivars, growing conditions, and growth stages [8, 23, 24]. Furthermore, three candidate genes for controlling the GS content in radish roots were identified from single nucleotide polymorphism (SNP) markers developed with GS . However, molecular mechanisms underlying GS metabolism in radish still require elucidation, especially for identification of the full set of genes involved in these related pathways.
In the present study, NGS-based Illumina paired-end solexa sequencing platform was employed to characterize the fleshy taproot de novo transcriptome in radish. A large set of radish transcript sequences were obtained to discover the majority of the activated genes involved in radish taproot. The candidate genes involved in the glucosinolate metabolism and regulation were successfully identified in radish. The sequence of representative genes and expression patterns were further validated. The root de novo transcriptome was comprehensively characterized in radish. This would provide a public information platform for understanding the molecular mechanisms involved in the metabolism of nutritional and flavor components during taproot formation, and facilitate the genetic improvement of quality traits in radish molecular breeding programs.
Results and discussion
Illumina sequencing and de novoassembly of radish root transcriptome
Statistics of output sequencing
Total raw reads
Total clean reads
Total clean nucleotides (nt)
Statistics of assembly quality
Total length (nt)
Mean length (nt)
Total consensus sequences
Functional annotation and classification of the assembled unigenes
Summary statistics of functional annotation for radish root unigenes in public databases
Public protein database
No. of unigene hit
Identification of candidate genes involved in the glucosinolate metabolism of radish
Initially, the parent amino acid is deaminated to form the corresponding 2-oxo acid by a branched-chain amino acid aminotransferase (BCAT, K00826, EC: 126.96.36.199). In A. thaliana, there are seven genes encoding the BCATs, and it is known to be fairly well conserved . In our annotated radish transcriptome unigene dataset, 17 sequences corresponding to five homologous BCAT genes (BCAT 2–5) were successfully identified. Subsequently, methylthioalkylmalate synthase (MAM, K15741, EC: 2.3.3.-) catalyzes 2-oxo acid condense with acetyl-CoA to yield a 2-oxo acid with one more methylene group (−CH2–) than the starting compound. Hereupon, the elongated 2-oxo acid can enter the core glucosinolate structure pathway or proceed through another round of chain elongation. Seven sequences encoding MAM were discovered in our transcriptomic analysis.
The formation of primary glucosinolates involved in core structure biosynthesis is accomplished through five different biochemical steps that synthesize several intermediates. It begins with the oxidation of the precursor amino acids to aldoximes by cytochromes P450 belonging to the CYP79 family, which is composed of a number of catalytic subfamilies. Genome analyses have revealed that Arabidopsis contains seven different CYP79 genes (i.e. CYP79A1, B1, B2, B3, F1, F2 and F3) . In the current study, ten unigene sequences were identified corresponding to the seven different genes with a high homology to CYP79s. All these seven gene members in the Arabidopsis genome were also identified in the radish transcriptome, which further confirmed the close relationship between these two species. Aldoximes are further oxidized to activated compounds (either nitrile oxides or aci-nitro compounds) by cytochromes P450 of the CYP83 family. Based on sequence similarities, four unigenes were identified corresponding to the two CYP83 genes (CYP83A1 and CYP83B1). The activated aldoximes are conjugated with cysteine as a sulfur donor to produce S-alkyl-thiohydroximates; however, it is not clear whether this conjugation is enzyme-mediated. The S-alkylthiohydroximate conjugates are converted to thiohydroximates by the C-S lyase SUPERROOT1 (SUR1, K11819) . In the present study, 12 homolog sequences were discovered encoding SUR1. Thiohydroximates are in turn S-glucosylated by glucosyltransferases of the UGT74 family to form desulfoglucosinolates. Overall, 13 unigenes were identified as UGT74s including UGT74B1, C1, F1 and F2. The final step in the synthesis of the GS core structure was catalyzed by desulfoglucosinolate sulfotransferase (SOT, K11821). There are three close homologous SOT genes (SOT16, 17and18), which were identified in Arabidopsis to catalyze this reaction with a wide variety of desulfoglucosinolate substrates . A total of 11 unigenes from our RNA-seq dataset were identified as SOTs including all three homologies found in Arabidopsis.
The initially produced parent glucosinolate from core structure is subject to a wide range of side chain modifications, which entail various kinds of reactions including oxidations, eliminations, alkylations, and esterifications. Kliebenstein et al. (2011) identified three genes responsible for side chain modification of aliphatic glucosinolates in Arabidopsis by QTL analyses , named GS-OX, GS-AOP and GS-OH; and functionally characterized two genes including AOP2, AOP3 of the GS-AOP cluster. In this study, 20 unigenes ranging from 252 bp to 1,921 bp were homologous to the genes encoding GS-OX; however, the other genes corresponding to the modification of side chain could not be identified.
Upon plant damage, the GS can be degraded to a variety of hydrolysis products such as isothiocyanates, oxazolidine-2-thiones, nitriles, epithionitriles, and thiocyanates. The hydrolytic process is catalyzed by a Beta-thioglucoside glucohydrolase (myrosinase, EC 188.8.131.52, K01188). Until now, myrosinase genes have been isolated from many plant species such as turnip, A. thaliana and mustard , which indicated that these genes are encoded by a multigene family and were classified into four subtypes(MA, MB, MC and TGG) on the basis of amino acid sequences . Additionally, two cDNA clones of myrosinase were isolated from radish seedlings, and both of them were identified as B type myrosinases . In this study, 14 unigenes were found which were homologs of genes encoding myrosinase, and most of them were predicted as MB subtypes.
Identification of genes involved in MYB transcription factors
MYB transcription factors represent a family of proteins that include the conserved MYB DNA-binding domain, which can control diverse pathways and processes corresponding to plant secondary metabolism [47, 48]. It was reported that many members of the MYB family could regulate the expression of related genes at the transcriptional level to control the process of GS metabolism in A. thaliana. For example, MYB28, 29 and 76 exerted a specific and coordinated control on the regulation of aliphatic GS biosynthesis, while MYB34, 51 and 122 could regulate the synthesis of indolic GS [49, 50]. From our radish transcriptome analysis, a total of 257 unigenes were predicted to code MYB proteins including a large number of members (i.e., MYB 2, 3, 4, 25, 28, 29, 43, 47, 52, 56, 58, 65, 69, 73, 78, 95, 103, 108, 121, etc.) (Additional file 4). However, the specific function of the particular MYB member in GS metabolism of radish need to be further verified with functional genomics approach.
Validation and expression analysis of genes involved in GS metabolism
Sequence analyses of the eight putative radish genes involved in glucosinolate metabolism process
In this study, NGS-based Illumina paired-end solexa sequencing platform was employed to characterize the fleshy taproot de novo transcriptome in radish. 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. A total of 67,305 unigenes were successfully annotated by blastx analysis using the publicly available protein database. It was revealed that the main genes activated in radish taproot, were predominately involved in basic physiological and metabolic processes, biosynthesis of secondary metabolites, signal transduction mechanisms, and other cellular components and molecular function related terms based on their matches in the GO, COG and KEGG databases. This study demonstrated that the Illumina paired-end sequencing technology is a fast and cost-effective method for novel gene discovery in non-model plant organisms. Furthermore, radish unigenes provided a comprehensive enough coverage to allow for the discovery of almost all genes known to be involved in GS metabolism and regulation related pathways. Our transcriptome dataset will serve as a valuable public platform to enhance the understanding of molecular mechanisms underlying biosynthesis and metabolism of the nutritional and flavor components during taproot formation. It would further facilitate the genetic improvement of major quality traits in radish breeding programs.
The radish (Raphanus sativus L.) advanced inbred line, ‘NAU-RG’, was used in this study. The surface-sterilized seeds were sown into soil in plastic pots and the seedlings were cultured in a growth chamber with 14 h light at 25°C and 10 h dark at 18°C. For Solexa analysis and T-A cloning sequencing, taproots were sampled at three different developmental stages including seedling, taproot thickening, and mature stages. The subsamples of root, leaf and stem parts were collected at seedling, taproot thickening, and mature stages, respectively for qRT-PCR verification (the skin and flesh at mature stage were separated). All samples were washed with distilled water, immediately frozen in liquid nitrogen and stored at −80°C for RNA extraction.
RNA extraction and Illumina sequencing
Total RNA of the three taproot samples from different stages was isolated using the RNAprep pure Plant Kit (Tiangen Biotech Co., Ltd., China) according to the manufacturer’s protocol. RNA samples were treated with RNase-free DNase I (Takara, Japan) to avoid DNA contamination. cDNA was prepared by equally pooling a total of 10 μg of RNA from each of the taproot sample of three different developmental stages. The mixed root cDNA library named ‘CKA’ was constructed using an mRNA-seq assay for paired-end transcriptome sequencing, which was performed by the Beijing Genomics Institute (BGI, Shenzhen, China).
Poly(A) mRNA was enriched from total RNA by using Sera-mag Magnetic Oligo (dT) Beads (Thermo Fisher Scientific, USA) and then mRNA-enriched RNAs were chemically fragmented to short pieces using 1× fragmentation solution (Ambion, USA) for 2.5 min at 94°C. These short fragments were taken as templates for first-strand cDNA synthesis using random hexamer-primer. The second-strand cDNA was generated using the SuperScript Double-Stranded cDNA Synthesis Kit (Invitrogen, USA). Short fragments were purified with QiaQuick PCR extraction kit and resolved with EB buffer for end repair and tailing A. Thereafter, the short fragments were connected with sequencing adapters, and the suitable fragments were selected for the PCR amplification as templates after agarose gel electrophoresis. Finally, the library was sequenced using Illumina HiSeq™ 2000.
Raw sequence processing and de novoassembly
Raw reads generated by Illumina Hiseq™ 2000 were initially processed to get clean reads. Then, all the clean reads were assembled using a de novo assembly program Trinity . Firstly, clean reads with a certain length of overlap were combined to form longer contiguous sequences (contigs), and then these reads were mapped back to the contigs. The distance and relation among these contigs was calculated based on paired-end reads, which enabled the detection of contigs from the same transcript and also the calculation of distances among these contigs. Finally, the contigs were further assembled using Trinity, and the contigs that could not be extended on either end were defined as unique transcripts. Additionally, the unigenes were divided into two classes by gene family clustering. The prefix CL was given to the clusters following the cluster id. Several unigenes with over 70% similarity were included from one cluster while from the other group the unigenes selected were singletons, for which the prefix unigene was used.
Functional annotation and classification of the assembled transcripts
All of the assembled transcripts were compared with the publicly available protein databases including NCBI non-redundant protein (Nr), Gene Ontology (GO), Clusters of Orthologous Groups (COGs), Swiss-Prot protein and the Kyoto Encyclopedia of Genes and Genomes (KEGG), using the BLASTx analysis with a cut-off E value of 10-5. The best alignments were used to identify sequence direction and to predict the coding regions of the assembled unigenes. If the results from different databases conflicted with each other, a priority order of nr, Swiss-Prot, KEGG and COG was followed. When a unigene happened to be unaligned to none of the above databases, software ESTScan was introduced to decide its sequence direction . For the nr annotations, the BLAST2GO program was used to get GO annotations of unique assembled transcripts for describing biological processes, molecular functions, and cellular components . After getting GO annotations for each transcript, WEGO software  was used to conduct GO functional classification for understanding the distribution of gene functions at the macroscopic level.
Gene validation by T-A cloning and sequencing
Specific PCR primers of the eight selected genes (Additional file 5) were designed corresponding to the conserved region of radish EST sequences from radish cDNA library . PCR was performed in a total volume of 25 μl containing 2.0 mmol/L Mg2+, 0.15 mmol/L dNTPs, 0.4 mmol/L of each primer, 0.8 U Taq DNA polymerase (TAKARA) and 15 ng cDNA with the following conditions: an initial denaturation step at 94°C for 1 min, 35 cycles at 94°C for 50 s, 56°C for 50 s, and 72°C for 90 s, a final extension at 72°C for 10 min and hold at 4°C. The PCR products were separated and ligated into the pMD18-T vector (Takara Bio Inc., China), and then transformed into E. coli DH5α. Positive clones were sequenced with ABI 3730 (Applied Biosystems, USA).
Quantitative real-time PCR (qRT-PCR) analysis
Quantitative real-time PCR was performed on a MyiQ Real-Time PCR Detection System (Bio-Rad) platform using the SYBR Green Master ROX (Roche, Japan) following the manufacturer’s instructions. Primers were designed using Beacon Designer 7.0 software, and Actin2/7 (ACT) (Additional file 6) was selected as the internal control gene . Amplification was achieved by a PCR program having a first denaturation step at 95°C for 5 min, then 40 cycles of denaturation at 95°C for 5 s, followed by annealing and extension at 58°C [30, 59]. The relative expression levels of the selected transcripts were normalized to ACT gene and calculated using the 2-∆∆Ct method. All reactions were performed in three replicates, and the data were analyzed using the Bio-Rad CFX Manager software.
Availability of supporting data
The RNA sequence dataset supporting the results of this article is available in the [NCBI Sequence Read Archive] repository, [SRX316199 and http://www.ncbi.nlm.nih.gov/sra/].
This work was in part supported by grants from the National Key Technology R & D Program of China (2012BAD02B01, 2011GB23600006), NSFC (31171956, 31372064), Key Technology R & D Program of Jiangsu Province (BE2010328), FRFCU(KYZ201209), JASTIF [CX (12)2006] and the PAPD.
- Wang L, He Q: Chinese radish. 2005, Beijing, China: Scientific and Technical Document Publishing HouseGoogle Scholar
- Curtis IS: Genetic transformation of radish (Raphanus sativus L.) by floral-dipping. Transgenic Crops of the World. Edited by: Curtis IS. 2004, Dordrecht: Kluwer Academic Publishers, 271-280.View ArticleGoogle Scholar
- Khanum F, Swamy MS, Krishna KS, Santhanam K, Viswanathan K: Dietary fiber content of commonly fresh and cooked vegetables consumed in India. Plant Foods Hum Nutr. 2000, 55 (3): 207-218. 10.1023/A:1008155732404.View ArticlePubMedGoogle Scholar
- Curtis IS: Genetic engineering of radish: current achievements and future goals. Plant Cell Rep. 2011, 30 (5): 733-744. 10.1007/s00299-010-0978-6.View ArticlePubMedGoogle Scholar
- Vardhini BV, Sujatha E, Rao SSR: Studies on the effect of brassinosteroids on the qualitative changes in the storage roots of radish. Bulg J Agric Sci. 2012, 18 (1): 63-69.Google Scholar
- Hara M, Ito F, Asai T, Kuboi T: Variation in amylase activities in radish (Raphanus sativus) cultivars. Plant Foods Hum Nutr. 2009, 64 (3): 188-192. 10.1007/s11130-009-0129-9.View ArticlePubMedGoogle Scholar
- Hara M, Sawada T, Ito A, Ito F, Kuboi T: A major β-amylase expressed in radish taproots. Food Chem. 2009, 114 (2): 523-528. 10.1016/j.foodchem.2008.09.082.View ArticleGoogle Scholar
- Ishida M, Kakizaki T, Ohara T, Morimitsu Y: Development of a simple and rapid extraction method of glucosinolates from radish roots. Breed Sci. 2011, 61 (2): 208-211. 10.1270/jsbbs.61.208.View ArticleGoogle Scholar
- Holst B, Williamson G: A critical review of the bioavailability of glucosinolates and related compounds. Nat Prod Rep. 2004, 21 (3): 425-447. 10.1039/b204039p.View ArticlePubMedGoogle Scholar
- Usuda H: Effects of growth under elevated CO2 on the capacity of photosynthesis in two radish cultivars differing in capacity of storage root. Plant Prod Sci. 2004, 7 (4): 377-385. 10.1626/pps.7.377.View ArticleGoogle Scholar
- Choi EY, Seo TC, Lee SG, Cho IH, Stangoulis J: Growth and physiological responses of Chinese cabbage and radish to long-term exposure to elevated carbon dioxide and temperature. Hortic Environ Biotechnol. 2011, 52 (4): 376-386. 10.1007/s13580-011-0012-0.View ArticleGoogle Scholar
- Samuolienė G, Sirtautas R, Brazaitytė A, Sakalauskaitė J, Sakalauskienė S, Duchovskis P: The impact of red and blue light-emitting diode illumination on radish physiological indices. Cent Eur J Biol. 2011, 6 (5): 821-828. 10.2478/s11535-011-0059-z.Google Scholar
- Wang S, Wang X, He Q, Liu X, Xu W, Li L, Gao J, Wang F: Transcriptome analysis of the roots at early and late seedling stages using Illumina paired-end sequencing and development of EST-SSR markers in radish. Plant Cell Rep. 2012, 31 (8): 1437-1447. 10.1007/s00299-012-1259-3.View ArticlePubMedGoogle Scholar
- Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37 (1): 1-13. 10.1093/nar/gkn923.PubMed CentralView ArticlePubMedGoogle Scholar
- Varshney RK, Nayak SN, May GD, Jackson SA: Next-generation sequencing technologies and their implications for crop genetics and breeding. Trends Biotechnol. 2009, 27 (9): 522-530. 10.1016/j.tibtech.2009.05.006.View ArticlePubMedGoogle Scholar
- Ward JA, Ponnala L, Weber CA: Strategies for transcriptome analysis in nonmodel plants. Am J Bot. 2012, 99 (2): 267-276. 10.3732/ajb.1100334.View ArticlePubMedGoogle Scholar
- Hyun TK, Rim Y, Jang H-J, Kim CH, Park J, Kumar R, Lee S, Kim BC, Bhak J, Nguyen-Quoc B: De novo transcriptome sequencing of Momordica cochinchinensis to identify genes involved in the carotenoid biosynthesis. Plant Mol Biol. 2012, 79 (4–5): 413-427.View ArticlePubMedGoogle Scholar
- Huang HH, Xu LL, Tong ZK, Lin EP, Liu QP, Cheng LJ, Zhu MY: De novo characterization of the Chinese fir (Cunninghamia lanceolata) transcriptome and analysis of candidate genes involved in cellulose and lignin biosynthesis. BMC Genomics. 2012, 13 (1): 648-10.1186/1471-2164-13-648.PubMed CentralView ArticlePubMedGoogle Scholar
- Shi CY, Yang H, Wei CL, Yu O, Zhang ZZ, Jiang CJ, Sun J, Li YY, Chen Q, Xia T: Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic pathways of tea-specific compounds. BMC Genomics. 2011, 12 (1): 131-10.1186/1471-2164-12-131.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Dong Y, Yang J, Liu X, Wang Y, Yao N, Guan L, Wang N, Wu J, Li X: De novo transcriptome of safflower and the identification of putative genes for oleosin and the biosynthesis of flavonoids. PLoS One. 2012, 7 (2): e30987-10.1371/journal.pone.0030987.PubMed CentralView ArticlePubMedGoogle Scholar
- Hua WP, Zhang Y, Song J, Zhao LJ, Wang ZZ: De novo transcriptome sequencing in Salvia miltiorrhiza to identify genes involved in the biosynthesis of active ingredients. Genomics. 2011, 98 (4): 272-279. 10.1016/j.ygeno.2011.03.012.View ArticleGoogle Scholar
- Liu S, Li W, Wu Y, Chen C, Lei J: De novo transcriptome assembly in chili pepper (Capsicum frutescens) to identify genes involved in the biosynthesis of capsaicinoids. PLoS One. 2013, 8 (1): e48156-10.1371/journal.pone.0048156.PubMed CentralView ArticlePubMedGoogle Scholar
- Ishii G: Glucosinolate in Japanese radish, Raphanus sativus L. Jpn Agric Res Q. 1991, 24: 273-279.Google Scholar
- Ishida M, Nagata M, Ohara T, Kakizaki T, Hatakeyama K, Nishio T: Small variation of glucosinolate composition in Japanese cultivars of radish (Raphanus sativus L.) requires simple quantitative analysis for breeding of glucosinolate component. Breed Sci. 2012, 62 (1): 63-10.1270/jsbbs.62.63.PubMed CentralView ArticlePubMedGoogle Scholar
- Zou Z, Ishida M, Li F, Kakizaki T, Suzuki S, Kitashiba H, Nishio T: 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. PLoS One. 2013, 8 (1): e53541-10.1371/journal.pone.0053541.PubMed CentralView ArticlePubMedGoogle Scholar
- Liang C, Liu X, Yiu S-M, Lim BL: De novo assembly and characterization of Camelina sativa transcriptome by paired-end sequencing. BMC Genomics. 2013, 14 (1): 146-10.1186/1471-2164-14-146.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang LL, Yang X, Sun P, Tong W, Hu SQ: The first Illumina-based de novo transcriptome sequencing and analysis of safflower flowers. PLoS One. 2012, 7 (6): e38653-10.1371/journal.pone.0038653.View ArticleGoogle Scholar
- Xie F, Burklew CE, Yang Y, Liu M, Xiao P, Zhang B, Qiu D: De novo sequencing and a comprehensive analysis of purple sweet potato (Impomoea batatas L.) transcriptome. Planta. 2012, 236 (1): 101-113. 10.1007/s00425-012-1591-4.View ArticlePubMedGoogle Scholar
- Logacheva M, Kasianov A, Vinogradov D, Samigullin T, Gelfand M, Makeev V, Penin A: De novo sequencing and characterization of floral transcriptome in two species of buckwheat (Fagopyrum). BMC Genomics. 2011, 12 (1): 30-10.1186/1471-2164-12-30.PubMed CentralView ArticlePubMedGoogle Scholar
- Dang ZH, Zheng LL, Wang J, Gao Z, Wu SB, Qi Z, Wang YC: Transcriptomic profiling of the salt-stress response in the wild recretohalophyte Reaumuria trigyna. BMC Genomics. 2013, 14 (1): 29-10.1186/1471-2164-14-29.PubMed CentralView ArticlePubMedGoogle Scholar
- Harris M, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C: The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004, 32: 258-261. 10.1093/nar/gkh036.View ArticleGoogle Scholar
- Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.PubMed CentralView ArticlePubMedGoogle Scholar
- Dutkowski J, Tiuryn J: Identification of functional modules from conserved ancestral protein–protein interactions. Bioinformatics. 2007, 23 (13): i149-i158. 10.1093/bioinformatics/btm194.View ArticlePubMedGoogle Scholar
- Zang YX, Kim HU, Kim JA, Lim MH, Jin M, Lee SC, Kwon SJ, Lee SI, Hong JK, Park TH: Genome-wide identification of glucosinolate synthesis genes in Brassica rapa. FEBS J. 2009, 276 (13): 3559-3574. 10.1111/j.1742-4658.2009.07076.x.View ArticlePubMedGoogle Scholar
- Wang H, Wu J, Sun S, Liu B, Cheng F, Sun R, Wang X: Glucosinolate biosynthetic genes in Brassica rapa. Gene. 2011, 487 (2): 135-10.1016/j.gene.2011.07.021.View ArticlePubMedGoogle Scholar
- Wittstock U, Halkier B: Glucosinolate research in the Arabidopsis era. Trends Plant Sci. 2002, 7 (6): 263-270. 10.1016/S1360-1385(02)02273-2.View ArticlePubMedGoogle Scholar
- Yang G, Gao Y-T, Shu X-O, Cai Q, Li G-L, Li H-L, Ji B-T, Rothman N, Dyba M, Xiang Y-B: Isothiocyanate exposure, glutathione S-transferase polymorphisms, and colorectal cancer risk. Am J Clin Nutr. 2010, 91 (3): 704-711. 10.3945/ajcn.2009.28683.PubMed CentralView ArticlePubMedGoogle Scholar
- Sønderby IE, Geu-Flores F, Halkier BA: Biosynthesis of glucosinolates–gene discovery and beyond. Trends Plant Sci. 2010, 15 (5): 283-290. 10.1016/j.tplants.2010.02.005.View ArticlePubMedGoogle Scholar
- Sun C, Li Y, Wu Q, Luo H, Sun Y, Song J, Lui E, Chen S: De novo sequencing and analysis of the American ginseng root transcriptome using a GS FLX Titanium platform to discover putative genes involved in ginsenoside biosynthesis. BMC Genomics. 2010, 11 (1): 262-10.1186/1471-2164-11-262.PubMed CentralView ArticlePubMedGoogle Scholar
- Diebold R, Schuster J, Däschner K, Binder S: The branched-chain amino acid transaminase gene family in Arabidopsis encodes plastid and mitochondrial proteins. Plant Physiol. 2002, 129 (2): 540-550. 10.1104/pp.001602.PubMed CentralView ArticlePubMedGoogle Scholar
- Grubb CD, Abel S: Glucosinolate metabolism and its control. Trends Plant Sci. 2006, 11 (2): 89-100. 10.1016/j.tplants.2005.12.006.View ArticlePubMedGoogle Scholar
- Mikkelsen MD, Naur P, Halkier BA: Arabidopsis mutants in the C-S lyase of glucosinolate biosynthesis establish a critical role for indole-3-acetaldoxime in auxin homeostasis. Plant J. 2004, 37 (5): 770-777. 10.1111/j.1365-313X.2004.02002.x.View ArticlePubMedGoogle Scholar
- Piotrowski M, Schemenewitz A, Lopukhina A, Müller A, Janowitz T, Weiler EW, Oecking C: Desulfoglucosinolate sulfotransferases from Arabidopsis thaliana catalyze the final step in the biosynthesis of the glucosinolate core structure. J Biol Chem. 2004, 279 (49): 50717-50725. 10.1074/jbc.M407681200.View ArticlePubMedGoogle Scholar
- Kliebenstein DJ, Gershenzon J, Mitchell-Olds T: Comparative quantitative trait loci mapping of aliphatic, indolic and benzylic glucosinolate production in Arabidopsis thaliana leaves and seeds. Genetics. 2001, 159 (1): 359-370.PubMed CentralPubMedGoogle Scholar
- Rask L, Andréasson E, Ekbom B, Eriksson S, Pontoppidan B, Meijer J: Myrosinase: gene family evolution and herbivore defense in Brassicaceae. Plant Mol Biol. 2000, 42 (1): 93-113. 10.1023/A:1006380021658.View ArticlePubMedGoogle Scholar
- Hara M, Fujii Y, Sasada Y, Kuboi T: cDNA cloning of radish (Raphanus sativus) myrosinase and tissue-specific expression in root. Plant Cell Physiol. 2000, 41 (10): 1102-1109. 10.1093/pcp/pcd034.View ArticlePubMedGoogle Scholar
- Jin H, Martin C: Multifunctionality and diversity within the plant MYB-gene family. Plant Mol Biol. 1999, 41 (5): 577-585. 10.1023/A:1006319732410.View ArticlePubMedGoogle Scholar
- Laitinen RA, Ainasoja M, Broholm SK, Teeri TH, Elomaa P: Identification of target genes for a MYB-type anthocyanin regulator in Gerbera hybrida. J Exp Bot. 2008, 59 (13): 3691-3703. 10.1093/jxb/ern216.PubMed CentralView ArticlePubMedGoogle Scholar
- Gigolashvili T, Engqvist M, Yatusevich R, Müller C, Flügge UI: HAG2/MYB76 and HAG3/MYB29 exert a specific and coordinated control on the regulation of aliphatic glucosinolate biosynthesis in Arabidopsis thaliana. New Phytol. 2008, 177 (3): 627-642. 10.1111/j.1469-8137.2007.02295.x.View ArticlePubMedGoogle Scholar
- Sønderby IE, Burow M, Rowe HC, Kliebenstein DJ, Halkier BA: A complex interplay of three R2R3 MYB transcription factors determines the profile of aliphatic glucosinolates in Arabidopsis. Plant Physiol. 2010, 153 (1): 348-363. 10.1104/pp.109.149286.PubMed CentralView ArticlePubMedGoogle Scholar
- Schuster J, Binder S: The mitochondrial branched-chain aminotransferase (AtBCAT-1) is capable to initiate degradation of leucine, isoleucine and valine in almost all tissues in Arabidopsis thaliana. Plant Mol Biol. 2005, 57 (2): 241-254. 10.1007/s11103-004-7533-1.View ArticlePubMedGoogle Scholar
- McCully ME, Miller C, Sprague SJ, Huang CX, Kirkegaard JA: Distribution of glucosinolates and sulphur-rich cells in roots of field-grown canola (Brassica napus). New Phytol. 2008, 180 (1): 193-205. 10.1111/j.1469-8137.2008.02520.x.View ArticlePubMedGoogle Scholar
- Li X, Kushad MM: Purification and characterization of myrosinase from horseradish (Armoracia rusticana) roots. Plant Physiol Biochem. 2005, 43 (6): 503-511. 10.1016/j.plaphy.2005.03.015.View ArticlePubMedGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29 (7): 644-652. 10.1038/nbt.1883.PubMed CentralView ArticlePubMedGoogle Scholar
- Iseli C, Jongeneel CV, Bucher P: ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc Int Conf Intell Syst Mol Biol. 1999, 1999: 138-148.Google Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34 (suppl 2): 293-297.View ArticleGoogle Scholar
- Jiang LN, Wang LJ, Liu LW, Zhu XW, Zhai LL, Gong YQ: Development and characterization of cDNA library based novel EST-SSR marker in radish (Raphanus sativus L.). Sci Hortic. 2012, 140: 164-172.View ArticleGoogle Scholar
- Xu Y, Zhu X, Gong Y, Xu L, Wang Y, Liu L: Evaluation of reference genes for gene expression studies in radish (Raphanus sativus L.) using quantitative real-time PCR. Biochem Biophys Res Commun. 2012, 424: 398-403. 10.1016/j.bbrc.2012.06.119.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.