- Research article
- Open Access
Transcriptome and metabolome analysis in shoot and root of Valeriana fauriei
© Park et al. 2016
- Received: 11 August 2015
- Accepted: 13 April 2016
- Published: 23 April 2016
Valeriana fauriei is commonly used in the treatment of cardiovascular diseases in many countries. Several constituents with various pharmacological properties are present in the roots of Valeriana species. Although many researches on V. fauriei have been done since a long time, further studies in the discipline make a limit due to inadequate genomic information. Hence, Illumina HiSeq 2500 system was conducted to obtain the transcriptome data from shoot and root of V. fauriei.
A total of 97,595 unigenes were noticed from 346,771,454 raw reads after preprocessing and assembly. Of these, 47,760 unigens were annotated with Uniprot BLAST hits and mapped to COG, GO and KEGG pathway. Also, 70,013 and 88,827 transcripts were expressed in root and shoot of V. fauriei, respectively. Among the secondary metabolite biosynthesis, terpenoid backbone and phenylpropanoid biosynthesis were large groups, where transcripts was involved. To characterize the molecular basis of terpenoid, carotenoid, and phenylpropanoid biosynthesis, the levels of transcription were determined by qRT-PCR. Also, secondary metabolites content were measured using GC/MS and HPLC analysis for that gene expression correlated with its accumulation respectively between shoot and root of V. fauriei.
We have identified the transcriptome using Illumina HiSeq system in shoot and root of V. fauriei. Also, we have demonstrated gene expressions associated with secondary metabolism such as terpenoid, carotenoid, and phenylpropanoid.
- Valeriana fauriei
- Illumina HiSeq system
- Digital gene expression
Terpenes, one of the major secondary metabolites in medicinal plants, have many volatile representatives such as isoprenes (C5), monoterpenes (C10), sesquiterpenes (C15), even some diterpenes (C20), and triterepenes (C30) . Members of this family serve many important functions like protecting the plants from many insects, pest, herbivores and microbial pathogens such as bacteria and fungi . All terpenes are assembled from the basic unit of isoprenes, isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP). These universal five carbon precursors are derived from two alternate biosynthetic pathways: mevalonate (MVA) pathway from acetyl-CoA and the 2-C-methylerythritol 4-phosphate (MEP) pathway from glycerol and pyruvate (Additional file 1) . Recent report claimed that even though the operation of this pathway is independent, some metabolic connection were observed between the two pathways .
The MVA pathway initiates with the formation of acetoacetyl-CoA by acetocetyl-CoA thiolase (AACT). The HMG-CoA synthase (HMGS) catalyzes the synthesis of 3-hydroxy-3-methylglutaryl-CoA with one acetyl-CoA and one acetoacetyl-CoA. Next, HMG-CoA reductase (HMGR) synthesizes the MVA. Further, enzymes such as mevalonate kinase (MK), phosphomevalonate kinase (PMK), and mevalonate diphosphate decarboxylase (MVD) catalyze the formation of IPP. The first step in MEP pathway, pyruvate and glyceraldehyde-3-phosphate combined and produce 1-deoxy-D-xylulose 5-phosphate by DOXP synthase (DXS) which is mainly observed in plastids. Then, the conversion of 1-deoxy-D-xylulose-5-phosphate to 2-C-methyl-D-erythritol 4-phosphate is carried out by DOXP reductoisomerase (DXR). MEP is transformed into 1-hydroxy-2-methyl-2-(E)-butenyl 4-phosphate by the action of various catalytic enzymes. Enzyme (E)-4-hydroxy-3-methylbut 2-enyl diphosphate reductase (HDR) catalyzes the synthesis of IPP and DMAPP . IPP and DMAPP are formed by the action of enzymes geranyl diphosphate synthase (GPS) and farnesyl diphosphate synthase (FPS), respectively . IPP and DMAPP are the backbone for the synthesis of all terpenes.
Carotenoids, the natural C40 isoprenoid products, are essential hydrophobic plant compounds which contribute to yellow, orange or red color . Carotenoids play an important role in photosynthesis, photomorphogenesis, and photoprotection. It is also mainly involved in the production of abscisic acid. The main carotenoid biosynthetic pathway was identified in many plants and studied extensively. The pathway begins with condensation of 3 molecules of geranyl-geranylpyrophosphate, a precursors from the upstream MEP pathway for production of 15-cis-phytoene catalyzed by phytoene synthase (PSY) (Additional file 2) . Then, 15-cis-phytoene is converted to lycopene, the first yellow carotenoid, through desaturation reactions which are synthesized by both phytoene desaturase (PDS) and ζ-carotene desaturase (ZDS) . Poly-cis lycopene to trans-lycopene is produced by the action of carotenoid isomerase (CrtISO) . Next, carotenoid biosynthesis is branched to produce α- and β-carotene by enzymatic catalytic activity of two lycopene cyclases, lycopene β-cyclase (LCYB) and lycopene ε-cyclase (LCYE) . α-carotene is hydroxylated into lutein by both β-ring hydroxylase (CHXB) and ε-ring hydroxylase (CHXE). In the other bransch, β-carotene is transformed into zeaxanthin; process is catalyzed by CHXB. Next, zeaxanthin epoxidase (ZEP) allows synthesis of violaxanthin from zeaxanthin. At last, the enzymatic activity of nine-cis-epoxycarotenoid dioxygenases (NCEDs) is responsible for the synthesis of abscisic acid as the final product . Along the pathway, the oxidative activity of the specific enzymes generate apocarotenoids which is further degraded by carotenoid cleavage dioxygenases (CCDs) .
Phenolic compounds, widely distributed in higher plants, belong to one of the major classes of secondary metabolites including lignins, flavonols, isoflavonoids and anthocyanins . These compounds contribute many important functional aspects of plant life such as UV sunscreens, pigments signaling. Additionally, accumulation of phenolic compounds is stimulated by biotic and abiotic responses. Currently, many researchers have been focused on the improvement of phenolic compounds in plants, because of its health promoting properties and curing properties to cancers, neurodegenerative diseases, cardiovascular diseases, osteoporosis and diabetes respectively . Phenolic compounds are synthesized through the phenylpropanoid pathway and its biosynthesis starts with the condensation of the phenylalanine which is end product of shikimate pathway (Additional file 3). In the first step, deamination reaction involved by the action of phenylalanine ammonia-lyase (PAL) for the generation of trans-cinnamic acid . Next, cinnamate 4-hydroxylase (C4H) and 4-coumaroyl CoA ligase (4CL) involved in the production of other intermediate metabolites . In several side branches, p-coumarate acid 3-hydroxylase (C3H) converts p-coumaric acid to caffeic acid which is transformed to ferulic acid, carried out by the enzyme caffeate O-mehtyltransferase (COMT) . Also, the enzyme hydroxycinnamoyl-CoA quinate hydroxyl cinnamoyl transferase (HQT) is involved in synthesis of p-coumaroylquinate from p-coumaroyl CoA which is then converted to chlorogenic acid by C3H . Finally, chalcone synthase (CHS), catalyzes naringenin-chalcone, which is the derivative of the flavanoids. In downstream steps, various enzymes such as isomerases, reductases, and hydroxylases involved in the alternation of the basic flavonoid skeleton, leading to the different flavonoid subclasses .
Recently, whole transcriptome sequencing using next-generation sequencing (NGS) technologies, or RNA sequencing (RNA-Seq) has been widely used for the characterization of genes and their functions in secondary metabolite synthesis . NGS technologies have efficiency due to its much higher levels of sensitivity, accuracy, dymamic-range of gene expression levels and fast run times (ranging from hours to days) compared to traditional low-throughput expressed sequence tag (EST) sequencing by Sanger technology [25, 26]. Also, many information such as quantitative gene expression, development of functional markets, quick insights into the specific gene space, comparative genomic studies and to isolate genes of interest [27, 28]. Among the common NGS platforms, including Illumina, Roche/454, SOLiD, and HelicosHeliScope, the Illumina HiSeq system has been commonly used because of its high throughput sequencing capacity which results in providing a higher coverage and low costs . Despite these advantages, the sequence reads provide sufficient for de novo assemblies of full-length transcripts, except in the case of small classes of RNA . For Valeriana family, only one transcriptome analysis from the Valeriana officinalis have been described in medicinal plant genomics resource (www.medicinalplantgenomics.msu.edu), whereas the transcriptome analysis of V. fauriei has not been investigated yet.
In this current study, we used the IlluminalHiSeq™ 2500 system to obtain the transcriptome in shoot and root of V. fauriei. We present data confirming that new sequencing technology can provide numerous insights into the molecular arrangement of secondary metabolite biosynthesis. In addition, we describe the identification of several full-length and partial-length cDNAs encoding genes related to terpenoid, carotenoid, and phenylpropanoid biosynthetic pathways. The transcripts levels of all genes were determined by real-time PCR and quantify the secondary metabolites with high-performance liquid chromatography (HPLC) to investigate the correlation between the transcriptional regulation of each biosynthetic gene and accumulation of each identified components. Till today no information has been provided the transcriptome characterization of V. fauriei.
Sequencing and transcriptome assembly
Summary of the transcriptome assembly of V. fauriei
346,771,454 (100 %)
35,023,916,854 (100 %)
291,047,352 (83.9 %)
28,761,361,637 (82.1 %)
143,401 (100 %)
70,145,151 (100 %)
97,959 (68.3 %)
61,543,817 (87.7 %)
Summary of the annotations of the V. fauriei transcriptome; GO, Gene Onthology; COG, Cluster Orthologous Groups; KEGG pathway, Kyoto Encyclopedia of Genes and Genomes
Database Ref IDs
97,959 (100 %)
Annotated with Uniprot BLAST hits
47,760 (48.7 %)
Annotated with GO
38,849 (39.6 %)
Annotated with COG
16,532 (16.8 %)
Annotated with KEGG
7041 (7.1 %)
50,199 (51.3 %)
Transcriptome is the full set of transcripts including mRNAs, non-coding RNAs, and small RNAs and their quantity for a specific developmental stage or physiological conditions. Transcriptome profiling is widely performed for clarifying all species of transcripts, confirming the transcriptional structure of genes, and quantifying the changing expression levels of each transcript during development and under different conditions . The development of next-generation sequencing (NGS) technology for whole transcriptome sequencing has offered high-throughput, advances in accuracy and sensitivity and decreased costs compared to traditional EST sequencing by Sanger technology . The improvements in genome sequencing technology have provided a valuable opportunity to sequence increasingly large and complex genome . Large scale sequencing of several non-model plants which are potential to investigate the basis of medicinal properties has been already assessed to several species i.e., Acacia auriculiformis, Acacia mangium , Cajanus cajan L. , Euphorbia fischeriana , Myricarubra , and so on. Although this technology has been extensively practiced to various research areas, genomic information of V. fauriei is still unknown. Therefore, transcriptome profiling using NGS from non-model plants is useful to generate a reference genome and to provide the basis of finding genes associated with particular important functions .
Digital gene expression and Secondary metabolite related gene
Changes in gene expression profile between shoot and root libraries
2 fold changes
Root VS Shoot
Summary for secondary metabolite genes involved in each biosynthetic pathway
No. unique transcripts
Flavone and flavonol biosynthesis
Isoquinoline alkaloid biosynthesis
Terpenoid backbone biosynthesis
Sesquiterpenoid and triterpenoid biosynthesis
Currently, several tools have developed to analyze the gene expression. Especially, digital gene expression (DGE) is getting popular due to its enrichment for a pathway or ontology term by using overlap statistics from variations in the counts of their cognate sequence tags [40, 41]. However, DGE requires a reference sequence to align the relative small read lengths . Indeed, performing DGE and RNA-Seq, which provides reference transcriptome, offers an efficient method to identify the candidate genes encoding enzymes involved in the biosynthesis of secondary metabolites in non-model plants .
Analysis of secondary metabolite biosynthetic genes from V. fauriei unigenes
The sequences of secondary metabolite biosynthetic pathway genes were identified in the NGS of the V. fauriei database. Several full-length cDNAs encoding MCT, HDS, GDS, AACT, HMGS, MK, PMK, IDI, and FDS, and partial-length cDNAs encoding DXS, DXR, CMK, MDS, HDR, and MVD were isolated from V. fauriei in terpenoid biosynthesis (Additional file 4). To confirm this for homology, they were designed as VfDXS (489 aa), VfDXR (473 aa), VfMCT (307 aa), VfCMK (310 aa), VfMDS (222 aa), VfHDS (734 aa), VfHDR (446 aa), VfGDS (418 aa), VfAACT (406 aa), VfHMGS (464 aa), VfHMGR (582 aa), VfMK (389 aa), VfCMK (497 aa), VfMVD (417 aa), VfIDI (235 aa), and VfFDS (345 aa) and showed sequence similarities according to BLAST search. Additional file 5 shows the sequences of carotenoid and phenylpropanoid biosynthetic genes identified from NGS data of V. fauriei. Among carotenoid biosynthetic genes, full-length cDNAs of VfPDS (569 aa), VfZDS (579 aa), VfCHXB (251 aa), and VfNCED (577 aa), and partial-length cDNAs of VfPSY (347 aa), VfCrtISO (401 aa), VfLCYB (363 aa), VfLCYE (249 aa), VfCHXE (403 aa), VfZEP (649 aa), and VfCCD (107 aa) were also exhibited. Additionally, full-length cDNAs of VfCOMT (240 aa), VfCHS (410 aa), VfF3H (353 aa), VfF3’H (326 aa), and VfFLS (332 aa) and partial-length cDNAs of VfPAL (468 aa), VfC4H (407 aa), Vf4CL (406 aa), VfC3H (373 aa), VfHQT (191 aa), VfCHI (172 aa), VfF3’5’H (172 aa), VfFNS (143 aa), VfFNS2 (333 aa), VfGT (175 aa), and VfRT (99 aa) for phenylpropanoid biosynthesis were purified from V. fauriei. A BLAST search at the amino acid level showed that secondary metabolite biosynthetic genes from V. fauriei exhibited high identity to other orthologous genes.
Analysis of terpenoid transcript levels and terpenoid content
Valerenic acid-related compounds content (μg/g dry weight) in root and shoot of V. fauriei. Each value is from three determinations ± SD. n.d, not detected
219.03 ± 6.70
32.22 ± 0.05
Until now, most of the studies have focused on root part of Valeriana species, whereas aerial part has been studied very rarely. Leaf contained valeric acid-related compounds as the major compounds . In addition, patchouli alcohol, α-pinene, and β-pinene were shown mainly in the oil of aerial parts of V. dioschoridis. Patchouli alcohol and isovaleric acid contributed major level of oil in the aerial parts of V. celtica [44, 45]. According to our study, p-cymene and pentanoic acid were detected as major compounds, whereas bornyl isovalerate is present as a little amount in shoot of V. fauriei. The major active compounds present in Valeriana species such as valerenic acid and its derivatives, and oxygenated sesquiterpenoids were accumulated in roots and rhizomes principally . As expected, valerenic acid-related compounds were shown only in root of V. fauriei. Besides, the main compounds of root of this plant were bornyl acetate, cedrol, α-acrorenol. These findings are quite similar to the result of Chen H et al.. They showed that bornyl acetate were the major constituent of the essential oil from the root and rhizomes of V. alternifolia.
Analysis of carotenoid transcript levels and carotenoid content
Carotenoid content (μg/g dry weight) in root and shoot of V. fauriei. Each value is from three determinations ± SD. n.d, not detected
0.42 ± 0.05
0.05 ± 0.01
0.47 ± 0.05
0.16 ± 0.05
10.38 ± 2.08
0.19 ± 0.04
2.99 ± 0.77
0.32 ± 0.11
14.04 ± 3.05
Carotenoids were produced in the photosynthetic organ (shoot) of V. fauriei in abundance. In contrast, the production of carotenoid in the underground organ (root) is rarely occurred. The results of the present study correspond well with those found in the earlier experimental studies in Chinese cabbage  and bitter melon .
Analysis of phenylpropanoid transcript levels and phenolic compound content
Phenolic compounds (μg/mg dry weight) in root and shoot of V. fauriei. Each value is from three determinations ± SD. n.d, not detected
3.819 ± 0.002
0.193 ± 0.001
0.236 ± 0.004
0.836 ± 0.000
0.004 ± 0.001
0.002 ± 0.001
14.726 ± 0.128
0.185 ± 0.001
0.015 ± 0.001
0.288 ± 0.006
0.002 ± 0.000
0.016 ± 0.000
In some studies, phenolic compounds including luteolin, apigenin, quercetin, kaempferol, and ferulic acid have been identified from both above ground biomass and root of Valeriana family . Also, Andres Navarrete et al. have demonstrated that chlorogenic acid quantified in a range of 0.2 to 2 % from V. jatamansi, V. procera, V. edulis, V. sitchensis, and V. officinalis . The well-known phenolic compounds i.e., trans-caffeic acid and rutin were confirmed in V. jatamansi . According to the result of Indra D. Bhatt et al., caffeic acid was highly detected in the aerial part of planted source, whereas chlorogenic acid was found in a higher amount in the root part that in aerial part of V. jatamansi . These findings are in contrast to our results. Here in this study, it was demonstrated that caffeic acid shows quite similarly in both root and shoot of V. fauriei. Besides, chlorogenic acid was more in shoot of V. fauriei where C3H and HQT genes expressed highly. The big differences with regard to chemical constituents indicated that these species are based on morphological .
Metabolic profiles between shoot and root of V. fauriei using GC-TOFMS analysis
The primary core metabolites such as organic acids, amino acids, and sugars determined using GC-TOFMS revealed clear metabolite differentiation between various biological samples. ChromaTOF software was used to support peak findings prior to quantitative analysis and for automated deconvolution of reference mass spectra. The NIST and the in-house libraries for standard chemicals were utilized for the identification of the compounds, which 42 metabolites were detected in the samples.
In further study, we expect that tissue-specific transcript profiling can provide insights into biological, functional differences between independent transcripts of both shoot and root of V. fauriei. In addition, there is an opportunity to produce medicinal substances for industrial purposes in both quality and quantity through genetic engineering.
In present study, we obtained a total of 97,959 unigenes using Illumina HiSeq system from V. fauriei. Among them, there were 47,760 annotated genes with Uniprot BLAST hits and mapped to GO, COG, and KEGG pathway. Individually, 70,013 transcripts and 88,827 transcripts were expressed in root and shoot, respectively of V. fauriei. Transcripts were highly involved in terpenoid backbone biosynthesis and phenylpropanoid biosynthesis compared to other biosynthesis by grouping into KEGG secondary metabolite biosynthetic pathway. Therefore, we found that de novo transciptome sequencing with DGE analysis is a modern novel technique for the identification of specific genes for the candidate enzymes involved in the biosynthesis of secondary metabolites in V. fauriei. We also investigated the expression of genes suggested an association with secondary metabolism in root and shoot of V. fauriei. Most of genes related to terpenoid biosynthesis are highly expressed in root than in shoot of V. fauriei. Also, we confirmed the presence of about 130 volatile compounds those are isolated from V. fauriei using GC/MS and also it is mentionable that valerenic acid-related compounds are shown only in the root of this plant. However, all carotenoid biosynthetic genes except NCED and CCD were expressed in a higher rate in the shoot than that of in the root of this plant. In addition, carotenoids were commonly accumulated in shoot. A couple of phenylpropanoid biosynthetic genes are expressed higher in the shoot and others are expressed in a higher rate in the root of V. fauriei showing in different amounts.
According to transcriptome databases, the comparison of gene expression and metabolite accumulation were achieved in this study. The transcriptome analysis for V. fauriei gives an opportunity to characterize genes leading to the synthesis of secondary metabolites, compounds of interest.
Plant materials and RNA isolation
V. fauriei seeds were collected from Rural Development Administration (RDA, Korea) and V. fauriei plants were established in a greenhouse at the experimental farm of Chungnam National University (Daejeon, Korea). Any temperature condition or additional illumination has not been regulated to the cultures. The plantlets were exposed to outdoor conditions in the greenhouse for 2 months between 19.05.2014 and 20.07.2014. Average temperature in this season was 21.8 °C. Plant materials were excised from 2-month-old young seedling plants and dissected into shoot (leave and stems) and root (roots and rhizomes). The separated plant parts were washed thoroughly with sterile water and frozen in liquid nitrogen immediately and stored at −80 °C. Shoot and root of V. fauriei were ground with liquid nitrogen. Total RNA was extracted from each part of V. fauriei separately using the Total RNA Mini Kit (Geneaid, Taiwan) by following the catalogues instruction and the concentrations were determined by agarose gel electrophoresis and NanoVue Plus Spectrophotometer (GE Healthcare Bio-Science Crop, USA) respectively. All samples were harvested and RNA isolation was performed in triplicate.
Illumina sequencing data was used for the identification of genes involved in terpenoid backbone, carotenoid, and phenylpropanoid metabolic biosynthetic pathway. The individual candidate gene name was searched using the functional annotation file. The selected amino acid sequences of the genes were analyzed for homology using BLAST at the NCBI Genbank database (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/BLAST).
Next generation sequencing of transcriptome
To obtain high-throughput transcriptome data of V. fauriei, we implemented Illumina-based NGS sequencing. Total RNAs were quantified using Nanodrop spectrophotometer (Thermo Scientific) and quality-assessed by RNA 6000 Nano assay kit (Agilent) and Bioanalyser 2100 (Agilent). NGS sequencing libraries were generated from one microgram of total RNA using Truseq RNA Sample Prep Kit (Illumina) according to the manufacturer’s protocol. In brief, the poly-A containing RNA molecules were purified using poly-T oligo attached magnetic beads. After purification, the total poly A+ RNA was fragmented into small pieces using divalent cations under elevated temperature. The cleaved mRNA fragments were reverse transcribed into first strand cDNA using random primers. QiaQuick PCR extraction kit was used for the purification of the shorts fragments and further resolved with EB buffer for end reparation and addition of poly (A). After that, the short fragments with poly (A) tail were interlinked with sequencing adapters. Each library was separated by adjoining distinct MID tag. The resulting cDNA libraries were then paired-end sequenced (2x101bp) with Illumina HiSeq™ 2500 system.
De novo assembly
Complete paired end sequences were obtained as individual fastq files (forward and reverse) from the images by CASAVA (version 1.8.2) base calling software with ASCII Q-score offset 33. Adaptor sequences and low quality bases with PHERD scores (Q) ≤ 20, were removed. Repeat sequences in raw reads were masked by using Repeat Makser against Arabiopdopsis Rep base database. Simple sequence repeats and low complexity sequences were also masked by using SSRIT and DUST software respectively . Finally, masked sequences were subjected to de novo assembly by using CLC Assembly Cell v.4.0 (CLCBio, Inc. Denmark) with customized parameters. To optimize transcriptome assembly, a set of assembly was done with word size of 21 to 63 with increment of 2 (data not shown). From the different assemblies, the best one was selected based on the sequence mapping coverage which was assessed by reference mapping of clean reads to an individual assembly with CLC mapper with length fraction 80 % and similarity 90 %. Based on the reference mapping results, the assembly with word size 63 was chosen as the best assembly. To remove the isoforms and obtain the non-redundant sequence set from the selected assembly contigs, CAP3 was performed with default parameters. Finally, results from CAP3 (singletons and contigs) were merged together to generate reference transcriptome assembly and renamed with unique sequence identifiers.
Non-redundant transcripts (per se reference transcriptome) were subjected to functional annotations by sequence homology search against biological databases such as GO, KEGG, and COG. Sequences were first mapped to UniProt database using BLASTX with e-value cut-off of 1e-5. Biological descriptions, GO terms, KEGG, and COG identifiers were transferred from the best-matched UniProt entry among the mapped sequence for each sequence. Then, the gene ontology functional classifications such as biological process (BP), molecular function (MF) and cellular components (CC), were grouped according to the GO hierarchy level 2 and the distributions were plotted using WEGO . Similarity, e-value, and species distributions were calculated from BLAST results. Transcripts were also grouped according to KEGG map IDs. Finally, sequence descriptions and references were collected from UniProt database using biopython module.
Digital gene expression (DGE) profiling and selection of transcripts involved in secondary metabolic pathway
To characterize the quantitative expression profile of individual sequence, the clean sequence reads from two libraries (root and shoot) were mapped individually to the reference transcriptome using CLC mapper with 90 % similarity and 80 % length fraction respectively. Based on the read count to each transcript, the reads per kilo base per million (RPKM) value was calculated. RPKM = (109 * C) / N * L. C is number of mapped reads per sequence, N is total number of mapped reads and L is length of the sequence . The statistical significance difference between the expression levels of each transcript within each pair of conditions were assessed by Audic and Claveriae’ method and false discovery (FDR) rate control [40, 59]. RPKM values were taken to Gene Spring 12.5 GX to calculate the fold changes (FC) through 2 libraries with default parameter. Fold changes were calculated for category such as root vs Shoot. Finally, sequence were filtered with RPKM ≥0.3, FDR >0.001 and FC ≥2.0 . Annotations obtained from the UniProt database such as gene descriptions, GO terms and sequence descriptions were manually examined from the known keywords (referred by published articles and KEGG pathways) for the candidate gene selections of secondary metabolism related.
cDNA synthesis and quantitative real-time PCR
cDNA was synthesized from 1 μg of total RNA using the ReverTra Ace-α Kit (Toyobo, Osaka, Japan). For all target genes and the 18S gene, an internal reference, primers were designed with GenScript Real-time PCR (TaqMan) Primer Design (www.genscipt.com) to conduct quantitative real-time PCR (Additional file 8). Quantitative real-time PCR was performed in a BIO-RAD CFX96 Real-time PCR system (Bio-Rad Laboratories, Hercules, CA) with the 2X Real-Time PCR smart mix (Solgent Co., Ltd. Daejeon, Korea) under the following conditions: pre-denaturation at 95 °C for 15 min, denaturation for 20 s at 95 °C, and reaction cycle was repeated for 39 cycles at annealing for 40 s, 20 s at 72 °C and final extension at 72 °C for 10 min. The differences between treatment means were evaluated from tree independent replicates for each sample. Then, average mean value and standard deviation value from three replicates of respective samples were analyzed.
Analysis of GC and GC-mass spectrometry
Ten gram of the fresh shoot and root of V. fauriei were weighed and transferred into 25 ml headspace vials. A fused-silica fiber covered with a 75 μm layer of carboxen/polydimethylsiloxane (CAR/PDMS) was used for absorption of the volatile compounds in the fresh shoot and root. The vials containing the samples and the solvents were kept at 25 °C for 20 min and was then removed from the vial and injected into the GC where analysis was performed at 250 °C for 3 min.
GC analysis was performed using an Agilent 6890 N GC mainframe equipped with an HP-5 (30 m × 0.32mmID, film thickness 0.25 μm) fused-silica capillary column (Agilent, USA) and a flame ionization detector. The injector and detector temperatures for each analysis were 250 and 280 °C, respectively. The carrier gas was nitrogen at a flow rate of 1.0 mL · min-1. The column temperature was maintained at 50 °C for 5 min and afterward programmed as follow : increase from 50 to 260 °C at a rate of 3 °C · min−1, increase from 260 to 280 °C at a rate of 10 °C · min−1 and hold at 280 °C for 5 min.
GC-MS analysis was performed on a GC/MSD Polaris Q (thermoFinnigan, USA) with an HP-5 (30 m × 0.32 mm ID, film thickness 0.25 μm) fused-silica capillary column (Agilent, USA). Helium was used as the carrier gas at a flow rate of 1.0 mL · min−1. For GC-MS detection, an electron ionization system with system energy 70 eV, trap current 250 μA, and ion source temperature 200 °C was used. The oven temperature program was the same as that described for GC, and injections were used in the split less mode. To identify of samples, components compared of mass spectra with the NIST and WILLY library data of the GC-MS system and with data from the literature. Total ion current chromatograms were recorded in a mass rage of 40–400 amu.
High performance liquid chromatography analysis
Collected samples were dried in the freeze-dryer at −80 °C for 3 days. Dried samples were ground into a fine powder using a mortar and pestle.
Quantification of valerenic acid and its derivatives
One gram of powdered samples was extracted with 10 mL of 90 % (v/v) methanol at room temperature for 30 min and the extracts were centrifuged at 12,000 rpm for 10 min. This step was repeated for three times. Thereafter, the final extract were evaporated using LABOROTA 4000 (Heidolph, Germany) and filtered with a 0.45-μm Acrodisc syringe filter (Pall Corp.; Port Washington, NY), for HPLC analysis. HPLC analysis was performed with a C18 column (μBondapak™ C18 10 μm 125Å 3.9 × 300 nm column). The mobile phase was a gradient prepared from mixtures of acetonitrile and 0.25 % phosphoric acid and the column temperature was maintained at 30 °C. The flow rate was maintained at 0.7 mL/min. Injection volume of 20 μL and 221 nm wavelengths were used for detection. The compounds of standard were determined by using a standard curve.
For carotenoids quantification, 300 mg of V. fauriei samples were mixed with 3 ml of ethanol containing 0.1 % ascorbic acid (w/v). This mixture was mixed throughly for a while, then incubated in a water bate at 85 °C for 5 min. In subsequent step, separately 120 μl of potassium hydroxide (80 % w/v) was added to saponify any potential interfering oils. After vortex, incubated at 85 °C for 10 min again and immediately keep it ice for 5 min. Separately, 1.5 ml of cold deionized water was mixed with 0.05 ml of b-Apo-80-carotenal (1.25 μg) was added to mixture as an internal standard. 1.5 ml of hexane was used for the extraction of extraction of complete carotenoids. The isolated carotenoids layers were separated by centrifugation at 1200 rpm for 10 min. The resulting extracts were freeze-dried by passing nitrogen gas and resuspended in 50:50 (v/v) dichloromethane/methanol. The carotenoids were separated on an Agilent 1260 HPLC system. The separating solvent consists of mixture of solvent A (methanol/water (92 % v/v) including 10 mM ammonium acetate) and solvent B consisted of 100 % methyl tert-butyl ether (MTBE). The separating solvents were flow at 1.0 ml/min and the injection sling was 20 μl. Samples were eluted with the following gradient: 0 min, 90 % A/10 % B; 20 min, 83 % A/17 % B; 29 min, 75 % A/25 % B, 35 min, 30 % A/70 % B; 40 min, 30 % A/ 70 % B; 42 min, 25 % A/75 % B; 45 min, 90 % A/10 % B; 55 min, 90 % A/10 % B. With the help of previous guidelines, based on the retention time and UV-visible spectrum data, carotenoids were identified.
Phenolic compounds extracted by mixing 100 mg samples with 80 % methanol at room temperature for 1 h in a 10 ml sterile test tube. After proper mixing, the crude mixture was centrifuged and the supernatant was filtered through 0.45 μm filter before HPLC analysis. Individual phenolic compounds were eluted, separated and quantified in a Futecs model NS-4000 HPLC apparatus (Daejeon, Korea). The HPLC analysis was performed in C18 column (maintained at 30 °C) and the phenolic compounds were detected at 280 nm. Gradient mobile phase system [water: acetic acid (98:2 v/v)] was used for the elution of the phenolic compounds. The gradient solvents were flown at 1.0 ml/min and 20 μl of the sample was used for injection. The contents of individual phenol compounds were calculated using a standard calibration curve. All samples were run in triplicate.
GC-TOFMS analysis of polar metabolites
The extraction of polar metabolite was performed as described by Kim et al. (2013) . The sample preparation methods, chemicals and reagents, esterification procedures, GC-TOFMS instrument operating conditions, analytical methods for the separation of the samples and the scanning and the detection of the compounds rang were as implemented as described in our previous research paper .
In V. fauriei samples, 42 metabolites were identified by GC-TOFMS. The quantitative calculations of all analytes were based on the peak area ratios relative to that of the IS. The relative quantification data acquired from GC-TOFMS were subjected to PCA (SIMCA-P version 13.0; Umetrics, Umeå, Sweden) to evaluate the relationships in terms of similarity or dissimilarity among groups of multivariate data. The PCA output consisted of score plots to visualize the contrast between different samples and loading plots to explain the cluster separation. The data file was scaled with unit variance scaling before all the variables were subjected to the PCA.
Ethics approval and consent to participate
Consent for publication
Availability of data and material
The datasets supporting the conclusions of this article are included within the article and its additional files.
This work was carried out with the support of “Cooperative Research Program for Agriculture Science & Technology Development (Project No. PJ011582022015)” Rural Development Administration, Republic of Korea and was supported by the grant entitled “The Genetic and Genomic Evaluation of Indigenous Biological Resources” funded by the National Institute of BIological Resources (NIBR201503101). The authors would like to extend their sincere appreciation to the Deanship of Scientific Research at King Saud University for its funding of this research through the Research Group project NO (RG-1435-071).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Osima Y, Matsuoka S, Ohizumi Y. Antidepressant principles of Valeriana fauriei roots. Chem Pharm Bull. 1995;43(1):169–70.View ArticleGoogle Scholar
- Sharma M, Jain UK, Patel A, Gupta N. A Comprehensive Pharmacognostic Report on Valerian. Int J Pharma Sci Res. 2010;1(7):6–40.Google Scholar
- Raal A, Orav A, Arak E, Kailas T, Müürisepp M. Variation in the composition of the essential oil of Valeriana officinalis L. roots from Estonia. Proc Estonian Acad Sci Chem. 2007;56(2):67–74.Google Scholar
- Liu X-G, Zhang W-C, Gao P-Y, Wang G-S, Li L-Z, Song S-J, Zhang X, Yao X-S, Liu K, Zhang Z-X. Two New Sesquiterpenes from the Roots of Valeriana fauriei Briq. Helvetica Chimica Acta. 2013;96:651–5.View ArticleGoogle Scholar
- Houghton PJ. The Scientific Basis for the Reputed Activity of Valerian. J Pharm Pharmacol. 1999;51:505–12.View ArticlePubMedGoogle Scholar
- Zhenxue Z, Xinsheng Y. Progress on chemical studies of Valeriana. Chin J Med Chem. 2000;10:226–30.Google Scholar
- Mathela CS, Chanotiya CS, Sammal SS, Pant AK, Pandey S. Compositional Diversity of Terpenoids in the Himalayan Valeriana Genera. Chem Biodivers. 2005;2(9):1174–82.View ArticlePubMedGoogle Scholar
- Bos R, Woerdenbag HJ, van Putten FM, Hendriks H, Scheffer JJ. Seasonal Variation of the Essential Oil, Valerenic Acid and Derivatives, and Valepotriates in < EM EMTYPE=. Planta Med. 1998;64(02):143–7.View ArticlePubMedGoogle Scholar
- Yuan C-S, Mehendale S, Xiao Y, Aung HH, Xie J-T, Ang-Lee MK. The Gamma-Aminobutyric Acidergic Effects of Valerian and Valerenic Acid on Rat Brainstem Neuronal Activity. Anesth Analg. 2004;98(2):353–8.View ArticlePubMedGoogle Scholar
- Dudareva N, Pichersky E, Gershenzon J. Biochemistry of Plant Volatiles. Plant Physiol. 2004;135:1893–902.View ArticlePubMedPubMed CentralGoogle Scholar
- Dudareva N, Negre F, Nagegowda DA, Orlova I. Plant Volatiles: Recent Advances and Future Perspectives. Crit Rev Plant Sci. 2006;25(5):417–40.View ArticleGoogle Scholar
- Nagegowda DA. Plant volatile terpenoid metabolism: Biosynthetic genes, transcriptional regulation and subcellular compartmentation. Fed Eur Biochem Soc. 2010;584:2965–73.View ArticleGoogle Scholar
- Hemmerlin A, Hoeffler J-F, Meyer O, Tritsch D, Kagan IA, Grosdemange-Billiard C, et al. Cross-talk between the Cytosolic Mevalonate and the Plastidial Methylerythritol Phosphate Pathways in Tobacco Bright Yellow-2 Cells. J Biol Chem. 2003;278:26666–76.Google Scholar
- Howitt CA, Pogson BJ. Carotenoid accumulation and function in seeds and non-green tissues. Plant, Cell Environ. 2006;29:435–45.View ArticleGoogle Scholar
- Shumskaya M, Wurtzel ET. The carotenoid biosynthetic pathway: Thinking in all dimensions. Plant Sci. 2013;208:58–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Sandmann G. Evolution of carotene desaturation: The complication of a simple pathway. Arch Biochem Biophys. 2009;483(2):169–74.View ArticlePubMedGoogle Scholar
- Clotault J, Peltier D, Berruyer R, Thomas M, Briard M, Geoffriau E. Expression of carotenoid biosynthesis genes during carrot root development. J Exp Bot. 2008;59(13):3563–73.View ArticlePubMedGoogle Scholar
- Cazzonelli CI. Carotenoids in nature: insights from plants and beyond. Funct Plant Biol. 2011;38:833–47.View ArticleGoogle Scholar
- Auldridge ME, McCarty DR, Klee HJ. Plant carotenoid cleavage oxygenases and their apocarotenoid products. Curr Opin Plant Biol. 2006;9(3):315–21.View ArticlePubMedGoogle Scholar
- Vogt T. Phenylpropanoid Biosynthesis. Mol Plant. 2010;3(1):2–20.View ArticlePubMedGoogle Scholar
- Christelle M, Andréa B, Schafleitnerc R, Legaya S, Lefèvrea I, Aliagac CAA, et al. Gene expression changes related to the production of phenolic compounds in potato tubers grown under drought stress. Phytochemistry. 2009;70(9):1107–16.Google Scholar
- Fraser CM, Chapple C. The Phenylpropanoid Pathway in Arabidopsis. Arabidopsis book 2011;e1052. doi:10.1199/tab.0152.
- Martens S, Preuß A, Matern U. Multifunctional flavonoid dioxygenases: Flavonol and anthocyanin biosynthesis in Arabidopsis thaliana L. Phytochemistry. 2010;71(10):1040–9.View ArticlePubMedGoogle Scholar
- Angeloni F, Wagemaker CAM, Jetten MSM, Camp HJMOD, Janssen-Megens EM, Francoijs K-J, et al. De novo transcriptome characterization and development of genomic tools for Scabiosa columbaria L. using next-generation sequencing techniques. Mol Ecol Resour. 2011;11:662–74.Google Scholar
- Glenn TC. Field guide to next-generation DNA sequencers. Mol Ecol Resour. 2011;11:759–69.View ArticlePubMedGoogle Scholar
- Xia J, Wang Q, Jia P, Wang B, Pao W, Zhao Z. NGS Catalog: A Database of Next Generation Sequencing Studies in Humans. Human Mutation Database in Brief. 2012;33(6):E2341–55.View ArticleGoogle Scholar
- Morozova O, Hirst M, Marra MA. Applications of New Sequencing Technologies for Transcriptome Analysis. Annu Rev Genomics Hum Genet. 2009;10(1):135–51.View ArticlePubMedGoogle Scholar
- Garg R, Jain M. RNA-Seq for transcriptome analysis in non-model plants. Methods Mol Biol. 2013;1069:43–58.View ArticlePubMedGoogle Scholar
- Strausberg RL, Levy S, Rogers Y-H. Emerging DNA sequencing technologies for human genomic medicine. Drug Discov Today. 2008;13:569–77.View ArticlePubMedGoogle Scholar
- Li R, Fan W, Tian G, Zhu H, He L, Cai J, Huang Q, Cai Q, Li B, Bai Y. The sequence and de novo assembly of the giant panda genome. Nature. 2009;463(7279):311–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Guo Q, Ma X, Wei S, Qiu D, Wilson IW, Wu P, Tang Q, Liu L, Dong S, Zu W. De novo transcriptome sequencing and digital gene expression analysis predict biosynthetic pathway of rhynchophylline and isorhynchophylline from Uncaria rhynchophylla, a non-model plant with potent anti-alzheimer’s properties. BMC Genomics. 2014;15(1):676.View ArticlePubMedPubMed CentralGoogle Scholar
- Imelfort M, Edwards D. De novo sequencing of plant genomes using second-generation technologies. Brief Bioinform. 2009;10(6):609–18.View ArticlePubMedGoogle Scholar
- Wong MM, Cannon CH, Wickneswari R. Identification of lignin genes and regulatory sequences involved in secondary cell wall formation in Acacia auriculiformis and Acacia mangium via de novo transcriptome sequencing. BMC Genomics. 2011;12(1):342.View ArticlePubMedPubMed CentralGoogle Scholar
- Kudapa H, Bharti AK, Cannon SB, Farmer AD, Mulaosmanovic B, Kramer R, Bohra A, Weeks NT, Crow JA, Tuteja R. A comprehensive transcriptome assembly of pigeonpea (Cajanus cajan L.) using Sanger and second-generation sequencing platforms. Mol Plant. 2012;5(5):1020–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Barrero RA, Chapman B, Yang Y, Moolhuijzen P, Keeble-Gagnère G, Zhang N, Tang Q, Bellgard MI, Qiu D. De novo assembly of Euphorbia fischeriana root transcriptome identifies prostratin pathway related genes. BMC Genomics. 2011;12(1):600.View ArticlePubMedPubMed CentralGoogle Scholar
- Feng C, Chen M, Xu C-J, Bai L, Yin X-r, Li X, Allan AC, Ferguson IB, Chen K-s. Transcriptomic analysis of Chinese bayberry (Myrica rubra) fruit development and ripening using RNA-Seq. BMC Genomics. 2012;13(1):19.Google Scholar
- Zhou Y, Gao F, Liu R, Feng J, Li H. De novo sequencing and analysis of root transcriptome using 454 pyrosequencing to discover putative genes associated with drought tolerance in Ammopiptanthus mongolicus. BMC Genomics. 2012;13(1):266.View ArticlePubMedPubMed CentralGoogle Scholar
- Tang Q, Ma X, Mo C, Wilson IW, Song C, Zhao H, Yang Y, Fu W, Qiu D. An efficient approach to finding Siraitia grosvenorii triterpene biosynthetic genes by RNA-seq and digital gene expression analysis. BMC Genomics. 2011;12:343.View ArticlePubMedPubMed CentralGoogle Scholar
- Audic S, Claverie J-M. The significance of digital gene expression profiles. Genome Res. 1997;7(10):986–95.PubMedGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Wei Z, Wanxia X, Jiewen L, Shujing S, Lei L, Shujing Z. Transcriptome profiling and digital gene expression analysis of Fallopia multiflora to discover putative genes involved in the biosynthesis of 2,3,5,4'-tetrahydroxy stilbene-2-O-β-D-glucoside. Gene. 2014;547(1):126–35.View ArticleGoogle Scholar
- Coassini Lokar L, Moneghini M. Geographical variation in the monoterpenes of Valeriana officinalis Leaf. Biochem Syst Ecol. 1989;17(7):563–7.View ArticleGoogle Scholar
- Bicchi C, Sandra P, Schelfaut M, Verzele M. Studies on the essential oil of Valeriana celtica L. J High Resolut Chromatogr. 1983;6(4):213–5.View ArticleGoogle Scholar
- Tzakou O, Couladis M, Pavlovic M, Soković M. Composition and antifungal activity of the oil from aerial parts and rhizomes of Valeriana dioscoridis from Greece. J Essent Oil Res. 2004;16(5):500–3.View ArticleGoogle Scholar
- Wang Y, Jin L, Yu S, Shi Q, Gu Y, Kiyota H. Chemical constituents of plants from the genus Valeriana. Mini-Reviews in Organic Chemistry. 2010;7(2):161–72.View ArticleGoogle Scholar
- Chen H, Xu X, Zhao V. Chemical constituents of essential oil from roots and rhizomes of Valeriana alernifolia Bunge. Zhongguo Zhong Yao Za Zhi. 1998;23(7):418–9. 447.PubMedGoogle Scholar
- Tuan PA, Park NI, Park WT, Kim YB, Kim JK, Lee J, Lee S-h, Yang T-J, Park SU. Carotenoids accumulation and expression of carotenogenesis genes during seedling and leaf development in Chinese cabbage (‘Brassica rapa subsp. pekinensis’). 2012.Google Scholar
- Tuan PA, Kim JK, Park NI, Lee SY, Park SU. Carotenoid content and expression of phytoene synthase and phytoene desaturase genes in bitter melon (Momordica charantia). Food Chem. 2011;126(4):1686–92.View ArticlePubMedGoogle Scholar
- Navarrete A, Avula B, Choi Y-W, Khan IA. Chemical fingerprinting of valeriana species: simultaneous determination of valerenic acids, flavonoids, and phenylpropanoids using liquid chromatography with ultraviolet detection. J AOAC Int. 2006;89(1):8–15.PubMedGoogle Scholar
- Tang Y-P, Liu X, Yu B. Two new flavone glycosides from Valeriana jatamansi. J Asian Nat Prod Res. 2003;5(4):257–61.View ArticlePubMedGoogle Scholar
- Bhatt ID, Dauthal P, Rawat S, Gaira KS, Jugran A, Rawal RS, Dhar U. Characterization of essential oil composition, phenolic content, and antioxidant properties in wild and planted individuals of Valeriana jatamansi Jones. Sci Hortic. 2012;136:61–8.View ArticleGoogle Scholar
- Yamazaki M, Mochida K, Asano T, Nakabayashi R, Chiba M, Udomson N, Yamazaki Y, Goodenowe DB, Sankawa U, Yoshida T. Coupling deep transcriptome analysis with untargeted metabolic profiling in Ophiorrhiza pumila to further the understanding of the biosynthesis of the anti-cancer alkaloid camptothecin and anthraquinones. Plant Cell Physiol. 2013;54(5):686–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee JM, Joung JG, McQuinn R, Chung MY, Fei Z, Tieman D, Klee H, Giovannoni J. Combined transcriptome, genetic diversity and metabolite profiling in tomato fruit reveals that the ethylene response factor SlERF6 plays an important role in ripening and carotenoid accumulation. The Plant Journal. 2012;70(2):191–204.View ArticlePubMedGoogle Scholar
- Temnykh S, DeClerck G, Lukashova A, Lipovich L, Cartinhour S, McCouch S. Computational and Experimental Analysis of Microsatellites in Rice (Oryza sativa L.): Frequency, Length Variation, Transposon Associations, and Genetic Marker Potential. Genome Res. 2001;11(8):1441–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang X, Madan A. CAP3: A DNA Sequence Assembly Program. Genome Res. 1999;9(9):868–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Shengting Li, Ruiqiang Li, Bolund L et al. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34:293–7.View ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.Google Scholar
- Kim JK, Park S-Y, Lim S-H, Yeo Y, Cho HS, Ha S-H. Comparative metabolic profiling of pigmented rice (Oryza sativa L.) cultivars reveals primary metabolites are correlated with secondary metabolites. J Cereal Sci. 2013;57(1):14–20.View ArticleGoogle Scholar