Transcriptome sequencing of rhizome tissue of Sinopodophyllum hexandrum at two temperatures
- Anita Kumari†1, 3,
- Heikham Russiachand Singh†2,
- Ashwani Jha2, 3,
- Mohit Kumar Swarnkar1,
- Ravi Shankar2Email author and
- Sanjay Kumar1Email author
© Kumari et al.; licensee BioMed Central Ltd. 2014
Received: 26 December 2013
Accepted: 17 September 2014
Published: 7 October 2014
Sinopodophyllum hexandrum is an endangered medicinal herb, which is commonly present in elevations ranging between 2,400–4,500 m and is sensitive to temperature. Medicinal property of the species is attributed to the presence of podophyllotoxin in the rhizome tissue. The present work analyzed transcriptome of rhizome tissue of S. hexandrum exposed to 15°C and 25°C to understand the temperature mediated molecular responses including those associated with podophyllotoxin biosynthesis.
Deep sequencing of transcriptome with an average coverage of 88.34X yielded 60,089 assembled transcript sequences representing 20,387 unique genes having homology to known genes. Fragments per kilobase of exon per million fragments mapped (FPKM) based expression analysis revealed genes related to growth and development were over-expressed at 15°C, whereas genes involved in stress response were over-expressed at 25°C. There was a decreasing trend of podophyllotoxin accumulation at 25°C; data was well supported by the expression of corresponding genes of the pathway. FPKM data was validated by quantitative real-time polymerase chain reaction data using a total of thirty four genes and a positive correlation between the two platforms of gene expression was obtained. Also, detailed analyses yielded cytochrome P450s, methyltransferases and glycosyltransferases which could be the potential candidate hitherto unidentified genes of podophyllotoxin biosynthesis pathway.
The present work revealed temperature responsive transcriptome of S. hexandrum on Illumina platform. Data suggested expression of genes for growth and development and podophyllotoxin biosynthesis at 15°C, and prevalence of those associated with stress response at 25°C.
Sinopodophyllum hexandrum (Royle) T.S. Ying, commonly known as Himalayan Mayapple, is an alpine herb of family Berberidaceae that grows on altitudes ranging between 2,400–4,500 m above mean sea level . The species is a perennial, erect, and unbranched herb reaching to a height of 40–50 cm; rhizome gives rise to multiple stems. Rhizome is non-edible, but is the source of a cytotoxic compound podophyllotoxin which is an aryltetralin lignan. Cytotoxicity of the compound is reduced for safe clinical use through its transformation into several derivatives. Some of the derivatives such as etoposide (VP-16) and teniposide (VM-26) have anti-cancerous activities and hence podophyllotoxin is in much demand in pharmaceutical industry . Extensive harvesting of the species from nature to meet the requirement of rhizome for podophyllotoxin extraction has led S. hexandrum to be listed as endangered species (Appendix II of convention for international trade in endangered species; http://www.cites.org/eng/app/appendices.php).
Several studies have reported on propagation and conservation of the species , podophyllotoxin biosynthesis , derivatization of podophyllotoxin, and their mode of action . Also, a few studies have been carried out on plant response to environmental cues. Like other alpine plant species , S. hexandrum has been reported to be sensitive to high temperature . S. hexandrum consumed more starch at 25°C as compared to at 10°C leading to poor root growth at the former temperature . Also, the rate of photosynthesis was reported to be lower at 30°C as compared to at 20°C. Higher temperature was also reported to increase the transpiration rate and reduce water use efficiency indicating the species to be sensitive at high temperature . We conducted a preliminary experiment to observe the growth of S. hexandrum at four different temperatures. S. hexandrum showed partial leaf folding and leaf drooping at 10°C and 35°C, respectively within a span of 3 weeks of temperature exposure, whereas the species did not exhibit such symptoms at 15°C and 25°C; though after 2–3 weeks of temperature, leaves lost shine and developed violet coloured spots at 25°C (Additional file 1). The species did not survive for long at 10°C and 35°C.
The objective of the present work was to study the molecular processes at 15°C and 25°C in commercially important organ rhizome to understand why the species performed better at 15°C. Transcriptome analysis offers a convenient route to understand molecular processes of the tissue in question. For example, analysis of populus (Populus tremula) transcriptome showed the importance of auxin and gibberellins signaling, and starch metabolism during winters ; in sessile oaks [Quercus petraea (Matt.) Liebl.], transcriptome analysis suggested the importance of cell rescue and defense during bud burst after dormancy ; subtracted transcriptome analysis in tea [Camellia sinensis (L.) O. Kuntze] identified the significance of genes involved in cell cycle/cell division and stress-inducible genes including those encoding chaperons during winter dormancy . Next generation sequencing has eased transcriptome analysis, and the technology has been successively used to study the transcriptome of several plant species including Himalayan plant species Picrorhiza kurrooa ( and the references mentioned therein) to discover genes and markers vis-à-vis to gain insight to important biological and molecular processes. More so, temperature responsive transcriptome of S. hexandrum rhizome is not yet reported.
Plants of S. hexandrum were collected from Parashar lake, Mandi (2,730 m above sea level; 30°12′ N 77°47′ E, India) and maintained at the institute in Palampur (1,300 m altitude; 32°06′ N, 76°33′ E, India) for one year in a polyhouse. During the collection procedure, the plants with rhizome of the similar size were selected for various experiments (India’s Biological Diversity Act 2002 permits access of biological resources to bonafide Indians for scientific research ). Plants from the polyhouse were shifted to two separate plant growth chambers (Percival Scientific, USA) maintained at 15°C and 25°C with 16 h photoperiod at a photosynthetic photon flux density of 200 μmol m-2 s-1. Plants were irrigated adequately and the rhizome tissue was harvested at day 0, 14, and 21 of transfer, cut quickly into smaller pieces using a sharp razor, mixed the cut pieces, frozen in liquid nitrogen and stored at -80°C for further analyses.
Extraction and estimation of podophyllotoxin
Frozen samples were ground to fine powder in liquid nitrogen using pestle and mortar followed by addition of 700 μl of 70% methanol for 100 mg tissue with intermittent grinding [14–17]. Extract was transferred to a centrifuge tube and the pestle and mortar was rinsed with 300 μl of 70% methanol to recover the left over sample. Extracts were pooled, vortexed for 5 minutes, sonicated (Ultrasonic Cleaner, MC-109-MP, Oscar, India) for 10 minutes at 25°C and centrifuged at 14,000 rpm for 5 minutes to collect supernatant for podophyllotoxin estimation. Supernatant was filtered through 0.22 μm filter (Millipore, USA) to estimate podophyllotoxin on an Ultra Performance Liquid Chromatography (UPLC) system. The system consisted of Acquity UPLC (Waters, Millford, USA) equipped with binary solvent manager, sample manager, photodiode array detector and a bridged ethyl hybrid workflow shield C18 (1.7 μm particles, 2.1 × 100 mm) analytical column (Waters Corp., Manchester, UK). Mobile phase consisted of methanol and water in a ratio of 60:40. Isocratic elution was carried out at a flow rate of 0.250 ml min-1 with injection volume of 5 μl. Podophyllotoxin was monitored at 240 nm and quantified using standard podophyllotoxin (Sigma, USA). Three separate biological replicates were used for each estimation. Significant difference in podophyllotoxin content at 15°C and 25°C was calculated using Completely Randomized Design to obtain p-value at different sampling times.
RNA extraction, cDNA preparation and transcriptome sequencing
Various protocols were followed essentially as described previously . In brief, total RNA was extracted from rhizome tissue as described by Muoki et al. . Nanodrop 1000 (NanoDrop Technologies, USA) and a Bioanalyzer Chip RNA 7500 series II (Agilent Technologies, USA) were used to determine quality and quantity of RNA. Oligotex mRNA Midi Kit (Qiagen, Germany) was used to purify mRNA from total RNA, and non-directional Illumina RNA sequencing library was prepared with mRNA-Seq 8 Sample Prep Kit (Illumina, USA) using random primers. Library had an average insert size of 200 bp and eight picomoles of the library [quantified using a Bioanalyzer Chip DNA 12000 series II (Agilent Technologies, USA)] was used for cluster generation. Paired end (PE) 36 × 2 bp sequencing was carried out on Illumina Genome Analyzer IIx (Illumina, USA) as per the manufacturer’s instructions.
Expression analysis and validation with quantitative real-time reverse transcriptase mediated polymerase chain reaction (qRT-PCR)
Total RNA was isolated as described previously . RNA was pretreated with RNase-free DNase I (Invitrogen, USA) and reverse transcribed with 1 μg of total RNA using Superscript III (Invitrogen, USA) according to the manufacturer’s instructions. Gene specific expression primers for qRT-PCR were designed using GenScript Real-time PCR Primer Design tool . qRT-PCR was performed with three separate biological replicates on a Stratagene Mx3000P System (Agilent Technologies, USA) using 2× Brilliant II SYBR® Green QPCR Master Mix (Agilent Technologies, USA). PCR was conducted under the following conditions: 10 min at 95°C (enzyme activation), 40 cycles each of 30s at 95°C, 30s at desired Tm (Additional file 2 has details), and 72°C for 30s. A final melting curve analysis was performed (55°C to 95°C) to verify specificity of amplicons. Transcript level of all the genes was normalized to an internal reference actin. Expression was estimated using Relative Expression Software Tool . Expression values were transformed (log2) to generate expression profiles. All primers used in this study are listed in Additional file 2.
De novoassembly and sequence clustering
Information regarding the assembled transcript sequences, associated filtered read data and transcript grouping for unigenes, is available at http://scbb.ihbt.res.in/Podo-12-12-11/. Entire computational analysis was carried out on CentOS Linux based 48 cores 2.2 Ghz AMD processors server with 256 GB random access memory.
To validate assembled transcripts, BLASTN of 1,463 ESTs, downloaded from National Centre for Biotechnology Information (NCBI) dbEST, was performed against transcriptome of S. hexandrum using E-value cut-off of 1e-05. Also, selected unigenes involved in podophyllotoxin biosynthesis were validated after DS clustering against complete coding sequences of corresponding genes of S. hexandrum submitted in GenBank.
Assembled transcripts of S. hexandrum were searched against UniProt database , Gene Ontology (GO) , Kyoto Encyclopedia of Genes and Genomes (KEGG)  and Enzyme Commission number (EC) , using BLASTX with a cut-off E-value of 10-1. Lower cut-off during annotation facilitates annotation even for instances that display only small region matching for some functional domain, which could be missed otherwise with higher cut-off values. Associated GO term was assigned for each assembled transcript by analyzing the GO term for the Uniprot sequence, which returned the highest scoring hit. Similar approach was followed for KEGG and EC classification and annotation for the assembled transcripts. Majority of GO, EC and KEGG based annotations and statistics were done using annotation tool, Annot8r  and an in-house developed scripts. Plant Transcription Factor Database (PlnTFDB)  hosts large number of plant specific transcription factors (TFs), their classification, corresponding nucleic acids and protein sequences. Data for all TFs reported in PlnTFDB, version 3.0, were downloaded for the current study. The assembled transcript sequences were searched against this database using BLASTX with an E-value threshold of 10-5. Top scoring homologous sequences from PlnTFDBs were used to annotate the assembled transcripts displaying highest similarity.
S. hexandrum transcript sequences were searched against MIPS Functional Catalogue (FunCat)  database to classify transcripts based on functional categories classified at MIPS. FunCat hierarchy is based on the Arabidopsis thaliana transcriptome. Therefore, BLASTX search was performed against the Arabidopsis transcriptome (downloaded from TAIR 10 ). The IDs from top scoring Arabidopsis gene hits were extracted and mapped to the S. hexandrum transcript sequences to categorize them into various FunCat categories.
Functional domain search for unknown sequences
The assembled sequences which did not return any homologous sequence hit during BLASTX search were converted into six longest ORFs. These ORFs were scanned against the functional domain databases like Conserved Domain Database (CDD) , using RPS-BLAST UNIX version.
Read mapping and transcript expression abundance measurement
Since a reference genome was not available in the present study, the filtered reads were mapped to the assembled transcripts. This was followed by estimation of total mapped reads. Uniquely mapped reads assigned to each assembled transcript with maximum two mismatches were allowed. Read mapping was done using TopHat which uses Bowtie to map the reads on transcripts  while differential abundance in terms of FPKM was measured using Cuffdiff  which deploys t-test and Benjamini-Hochberg correction to compute significance difference between samples and adjusts the p-value after false discovery rate for multiple hypothesis testing. Using the reads from 15°C and 25°C experiments, the relative abundance of transcripts was estimated for each unigene for the two temperatures for the sequence clusters generated following the DS clustering . DS clustering reduced the inflation of unique genes representation and mapped that information on every part of the study. The read data obtained for the two temperatures were mapped on the assembled transcripts separately. Care was taken to ensure that only those lanes and reads were considered for both the temperatures which had equal sample concentration and produced almost equal percentage of reads qualifying the filtering cut-off, ensuring that experimental conditions do not skew the read representation and analysis.
The associated GO terms (GO terms are derived from dynamic controlled vocabularies or ontologies that can be used to describe the function of genes and gene products) and IDs were mapped to each of these unique transcripts and genes, providing single stop complete information regarding their related expressions at the two temperatures and their corresponding annotations.
Guanine-cytosine (GC) content analysis and simple sequence repeats (SSRs) identification
GC content of the sequences was measured using Emboss GeeCee tool, while sequences were scanned for SSR markers using MISA .
Metabolic networks and pathways analysis
Sequences were scanned against KEGG using BLAST. FPKM was considered for all the assembled sequences. Plant Metabolic Network (PMN) database  was searched for all the significantly up-regulated genes using an in-house developed network running script. Using the same script, pathways network images were generated and corresponding transcript, expression and other associated information were incorporated into the files.
Results and discussion
Rhizome biology for S. hexandrum is of interest not only that rhizome is the source of medicinally important compound podophyllotoxin, but also it allows plant to survive the unfavorable environment of winters when aerial portion of the plant dies. Plant regenerates from the rhizome upon arrival of the growing season. Availability of vast genomic and transcriptomic data in the public domain coupled with next generation sequencing platforms made it possible to understand temperature responsive biology of rhizome of S. hexandrum. Since plants appeared healthier at 15°C as compared to at 25°C (Additional file 1), the present work analyzed transcriptome in the rhizome tissue at the two temperatures vis-à-vis podophyllotoxin content.
Podophyllotoxin accumulation at two temperatures
Reads generation and de novosequence assembly
Summary of transcriptome data generated on Illumina Genome Analyzer IIx for rhizome tissue of S . hexandrum
Total number of paired end reads
Number of reads obtained after quality filtering
Number of primary assembled transcripts of pooled data
Average length of transcripts (bp) of pooled data
Average coverage (X) of pooled data
Homology search, sequence clustering, validation of assembled transcripts and general traits of the transcriptome
The number of unique assembled transcripts reduced from 60,089 to 59,244 sequences after removing the redundancies. These were clustered together under a common Unigene group, and BLASTX  with cut-off E-value of 10-5 showed that a total 25,395 sequences had significant BLAST hit. Remaining 33,849 sequences without significant BLAST hit did show conserved domains in 761 sequences (Additional file 4), suggesting a possible function for these unknown sequences. DS clustering  further reduced the total sequences from 25,395 to 20,387 (Additional file 5; URL: http://scbb.ihbt.res.in/Podo-12-12-11/).
A total of 1,463 ESTs of S. hexandrum were reported in NCBI dbEST, out of which BLAST hits were found for 975 ESTs (66.64%) when searched against the transcriptome. These ESTs showed an average identity of ~93.50% indicating good quality of assembled transcripts. Of the total, 19.69% of ESTs (192/975) had coverage greater than 90% whereas 52.82% ESTs (515/975) showed coverage of 50% or more. To evaluate correct assembly of unigenes, coding sequences of five selected genes (accession numbers GU324975.1, EU240218.2, GU196273.1, KF170372.1 and GU228507.1) involved in podophyllotoxin biosynthesis pathway were searched against unigenes using BLASTN. Analysis yielded average coverage of 95.05% (Additional file 6), suggesting that the unigenes were correctly assembled.
Average GC content of S. hexandrum transcripts was found to be around 44.59% out of which 80% of transcripts had the content in the range of 40-49%, a range comparable to that reported for dicots (44-47%)  (Additional file 7). GC content of genome indicates several features including genome stability and possible repeats dynamics . Assembled transcript sequences of S. hexandrum identified a total of 3,226 SSRs. Most abundant SSR group was of tri-nucleotide (54.40%) with GAA, TTC, TCT, CAC and CCA as predominant repeats, followed by di-nucleotide (36.38%) with prevalent occurrence of poly-CT (261), and mono-nucleotide (19.36%) SSRs represented by Poly-T (272) as a dominant repeat. Only a small fraction of tetra (32) and hexa (3) SSRs contributed to the pool (Additional file 7). SSRs are used in production and control of strains as well as dissemination of important genes apart from being the source of polymorphism and marker of genome .
Functional annotation and classification of the S. hexandrumtranscriptome
Annotation of 20,387 assembled sequences against GO database yielded significant annotation for 15,810 unique transcripts that were classified into biological processes, molecular functions and cellular component using plant specific GO slims. Functional classification in biological process category (Additional file 8) revealed that metabolic process, response to stimulus, cellular process and multicellular organismal development were among the highly represented groups, suggesting tissue to be metabolically active. Genes involved in DNA binding, catalytic and transferase activity were highly represented in molecular function category (Additional file 8), indicating dominance of gene regulation, signal transduction and enzymatically active processes. Among the cellular components (Additional file 8), genes for intracellular location and those encoding for membrane localized protein were most represented, which is a typical character of a eukaryotic cell.
Best EC classification was obtained for 8,172 assembled sequences after DS clustering. An analysis of top 50 abundant enzyme classes showed predominance of serine/threonine protein kinase enzyme (21.73%; Additional file 9A), which is known to be involved in regulation of cell proliferation, programmed cell death (apoptosis) and cell differentiation . Associated KEGG classification could be obtained for 8,759 assembled sequences. An analysis of top 50 KEGG pathways showed that plant hormone signal transduction pathways (5.74%) dominated the group (Additional file 9B), followed by plant-pathogen interaction. GO, EC and KEGG analysis showed rhizome to be a metabolically active tissue.
Global analysis of transcript abundance in response to temperature
PMN analysis using homologous genes exhibiting 2-fold or higher expression at 15°C is mentioned at http://scbb.ihbt.res.in/Podo-12-12-11/podo_pathway_go/Genes_at_15_degree_having_expression_2_fold_or_above/. Every network diagram contains the pathway information, associated homologous gene, transcript ID, and expression/abundance in terms of FPKM; each file name contains the associated transcript’s scaffold/contig ID in its prefix to make the browsing job easier. The file “podophyllum_analysis.doc” in the parent directory of the above link contains tabular representation and details about all such assembled transcripts, associated homologue reported in PMN, associated gene and associated PMN metabolic pathway/network.
Data showed up-regulation of pathways for respiration, photosynthesis, and phenylpropanoid biosynthesis at 15°C. These are fundamental pathways that would be crucial for growth and survival of plants at 15°C. Data showed up-regulation of ABA biosynthesis, ethylene and jasmonic acid, and is suggestive of their roles in plant growth and development at 15°C through functioning as signal molecule . Up-regulation of phenylalanine metabolism, fatty acid biosynthesis, starch and sucrose metabolism and PP pathway at low temperature is also reported in cassava (Manihot esculenta) transcriptome . This global analysis of gene expression provided a comprehensive data-set with each gene represented by its absolute expression level at two temperatures.
Validation of FPKM-based analysis by quantitative real-time polymerase chain reaction (qRT-PCR)
FPKM based expression values were validated by qRT-PCR using sixteen genes related to growth and development and stress response. These genes were NEDD8, CTPS, ATHB-14, ABTB2, AP2, CO, CSLA, CPN60, ZFP, AOX1a, STK, HSP60, HSP, AP2D, FRO and TINY. Cullin-associated NEDD8-dissociated protein 1 regulates cell division, stress responses and hormonal signaling . CTP synthase is involved in embryo development . Homeobox-leucine zipper protein ATHB-14 (HDZip) family of TFs acts in developmental processes and in the mediation of external signals to regulate plant growth . ABTB2, AP2 and CO were reported to play a role in morphogenesis in Arabidopsis [52–54]. CSLA is required for synthesis of cell wall polysaccharide mannan which serves as storage reserve during plant growth and development . CPN60 was essential for development of embryo and seedling .
Transcription factor and expression enrichment
Transcription factor families exhibiting significant difference in expression (based upon FPKM values) at 15°C and 25°C in rhizome of S . hexandrum
Median fold enrichment [15°C/25°C]
Median fold enrichment [25°C/15°C]
Thus, transcriptome analysis including TF, GO, FunCat, EC and PMN data (Figures 3, 4, 5 and 6) suggested that 15°C favored growth and development, whereas 25°C imposed stress on S. hexandrum; and also supported the observation of better growth of plant at 15°C (Additional file 1).
Expression of genes involved in podophyllotoxin biosynthesis at the two temperatures
Studies on biosynthetic pathway of podophyllotoxin in Podophyllum peltatum and Linum flavum established occurrence of dirigent protein mediated coupling of coniferyl alcohol to get pinoresinol by involving the action of dirigent protein oxidase (DPO) . Sequential conversion of pinoresinol to lariciresinol and secoisolariciresinol by the action of pinoresinol–lariciresinol reductase (PLR) was revealed by radiolabeled substrate study [76, 77]. Secoisolariciresinol is dehydrogenated to matairesinol by action of secoisolariciresinol dehydrogenase (SLD) . The enzymatic reactions between matairesinol to deoxypodophyllotoxin are not well characterized. A direct pathway from matairesinol to yatein is proposed by tracer experiments in Anthriscus sylvestris and pathway from yatein to podophyllotoxin via deoxypodophyllotoxin was suggested by feeding experiment , but not yet deciphered. Deoxypodophyllotoxin is considered as the precursor for biosynthesis of podophyllotoxin  as supported by feeding experiment in S. hexandrum and Linum flavum[83, 84]. Hydroxylation of deoxypodophyllotoxin by deoxypodophyllotoxin 7-hydroxylase (DOP7H) involved in the biosynthesis of podophyllotoxin is yet to be characterized  (Figure 8).
Transcriptome data identified cytochrome P450 (CYPs), methyltransferases (MTs) and UDP- glycosyltransferases (UGTs), possibly associated with the biosynthesis of podophyllotoxin
CYPs are membrane bound heme proteins, which catalyze hydroxylation, epioxidation, dealkylation and oxidation reactions. CYPs constitute the third known largest family of plant genes involved in numerous biosynthetic reactions resulting in production of plant hormones, defensive compounds, and secondary metabolites . In podophyllotoxin biosynthetic pathway, hydroxylation at position 7 of deoxypodophyllotoxin by deoxypodophyllotoxin 7-hydroxylase (DOP7H) leads to the synthesis of podophyllotoxin, but the enzyme and the gene are yet to be characterized [82, 91]. Therefore, detailed analysis of CYPs was important.
Transcriptome data identified 116 CYPs based on highest bit score and significant E-value. CYP-39 (scaffold11394_153.2) was found to be significantly up-regulated at 25°C by cuffdiff (Additional file 13). Fifteen CYPs showed more than two fold up-regulation, while 24 showed down-regulation at 15°C as compared to at 25°C (Additional file 14). These 39 CYPs could be the potential candidates for detailed analysis to identify the putative ShDOP7H.
MTs are transferase enzymes that participate in transfer of a methyl group from a donor to an acceptor. O-methylation plays important role in biosynthesis of lignan including podophyllotoxin . Ferulic acid and sinapic acid are methylated compounds and are precursor of monolignols (coniferyl and sinapyl alcohols), the moieties involved in lignin biosynthesis . MTs are essential in PP pathway: activity of CCoAOMT is essential for coniferyl and sinapyl alcohols biosynthesis ; and COMT is the key enzyme involved in methylation using hydroxycinnamates as substrates . Tracer studies in Anthriscus sylvestris revealed that conversion of matairesinol to yatein involves four steps that include hydroxylation, dual methylation and methylenedioxy bridge formation . Of the total 63 MTs identified, 18 MTs exhibited differential modulation by two fold and above (Additional file 13) wherein 2 MTs showed up-regulation, and 16 showed down-regulation at 15°C as compared to at 25°C (Additional file 15). Deeper analysis of these MTs could be promising to identify the genes associated with podophyllotoxin biosynthesis.
Glycosyltransferases (GTs) are super family of enzymes that catalyze transfer of sugars to a wide range of acceptor molecules. Glycosylation of natural products is generally catalyzed by UGTs of family 1 GTs . In higher plants, UGT superfamily glycosylates a broad array of aglycones including plant hormones, plant secondary metabolites and xenobiotics such as herbicides . Glycosylation contributes to increased water solubility, chemical stability and reduced chemical activity and thus plays a role in cell homeostasis, plant growth, development and in response to abiotic stress. Hydroxylated molecules are the common acceptors, while UDP-glucose is the most common donor in the GT catalyzed glycosyl group transferring reaction.
In the lignin biosynthesis pathway, lignin monomers like coumaryl, coniferyl alcohol and sinapyl alcohol are translocated in the form of glucosides from cytosol to cell wall. Involvement of UGTs in the biosynthesis of lignin was also reported . Enzyme activity is expected to be essential for vacuolar storage of otherwise toxic lignans and is shown to be correlated with lignan glucoside accumulation. A podophyllotoxin-glucose glucosidase has been isolated from Podophyllum peltatum. Thus it would be relevant to identify the UGTs involved in enzymatic synthesis of valuable glycoconjugates. BLAST search identified 35 unigenes encoding GTs. A total of 10 UGTs (Additional file 16) exhibited modulation by two fold and above at the two temperatures as analysed through FPKM (Additional file 13), of which 8 UGTs were down-regulated while 2 UGTs showed up-regulation at 15°C as compared to 25°C. In depth analysis of these modulated UGTs would be needed to identify the possible candidates associated with podophyllotoxin biosynthesis.
The present work analysed temperature (15°C and 25°C) responsive transcriptome of rhizomatous tissue of medicinally important endangered species S. hexandrum. A total of 60,089 assembled sequences, representing 20,387 unique genes were obtained on Illumina platform. Assembly of unigenes was validated by full-length sequence data of the available genes of S. hexandrum, which showed correct assembly. Transcriptome had an average coverage of 88.34X, average length of 543.11 bp, average GC content of 44.59% and abundance of trinucleotide SSR (54.40%).
Functional annotation and classification of the S. hexandrum transcriptome revealed metabolic process, response to stimulus, cellular process and multicellular organismal development to be the highly represented groups, suggesting rhizome to be metabolically active tissue.
FPKM data was validated by qRT-PCR data and the two platforms of gene expression showed significantly positive correlation emphasizing the confidence in the two different methods of gene expression. Data showed over-expression of genes and TFs associated with growth and development at 15°C, whereas those associated with stress tolerance over-expressed at 25°C. Also, various genes involved in podophyllotoxin biosynthesis were identified and their expression both by FPKM and the qRT-PCR showed up-regulation of eight out of twelve genes at 15°C. Interestingly, podophyllotoxin accumulation also showed increased trend of accumulation at 15°C. In-depth analyses of CYPs, MTs and UGTs yielded the potential candidate hitherto unidentified genes of podophyllotoxin biosynthesis. Data generated in the present work will contribute to future studies in the field of functional genomics and metabolic engineering for this important plant species. It will also be worthwhile to analyse microRNAs (miRNAs) because of their role in regulating gene expression either by blocking the process of translation or breaking the target transcripts . MiRNAs are core regulatory component defining spatio-temporal states and are involved in stress responses as well. Challenge would be to analyse miRNAs in this species since its genome sequence is not yet reported .
Availability of supporting data
Raw read files used for assembly were submitted to Short Read Archive (SRA) at NCBI under the accession ID SRX682844 (http://www.ncbi.nlm.nih.gov/sra/?term=SRX682844). Transcriptome sequences were submitted to Transcriptome Shotgun Assembly (TSA) project and deposited at DDBJ/EMBL/GenBank under the accession ID GBJY00000000 [http://www.ncbi.nlm.nih.gov/nuccore/GBJY00000000. The version described in this paper is the first version, GBJY01000000 (http://www.ncbi.nlm.nih.gov/Traces/wgs/?val=GBJY01)]. Sequences, which were less than 200 bp long or had a long stretch of “N” (>14 nucleotides) were not accepted by TSA. These sequences along with submitted sequences can be accessed at http://scbb.ihbt.res.in/Podo-12-12-11/.
Authors thank the Director, CSIR-IHBT for providing necessary facilities for research. AK is thankful to CSIR for senior research fellowship, HRS and AJ are thankful to CSIR for assistantship. The work represents IHBT communication number 2375.
The work was supported by the Council of Scientific and Industrial Research (CSIR), India through the projects entitled “Pathway engineering and system biology approach towards homologous and heterologous expression of high-value phytoceuticals (artemisinin, picrosides, morphine, withanolides, podophyllotoxin (NWP-008)”; “Plant diversity: studying adaptation biology and understanding/exploiting medicinally important plants for useful bioactives (SIMPLE-BSC 0109)”; “Genomics of medicinal plants and agronomically important traits (PlaGen-BSC 0107)”; and Computational Biology project, entitled “Genomics and informatics solutions for integrating biology” (GENESIS-BSC 0121).
- Shaw JMH: New combinations for the varieties of Sinopodophyllum hexandrum. Hanburyana. 2009, 4: 33-39.Google Scholar
- Stähelin HF, von Wartburg A: The chemical and biological route from podophyllotoxin glucoside to etoposide: ninth Cain memorial Award lecture. Cancer Res. 1991, 51: 5-15.PubMedGoogle Scholar
- Nadeem M, Palni LMS, Purohit AN, Pandey H, Nandi SK: Propagation and conservation of Podophyllum hexandrum Royle: an important medicinal herb. Biol Conserv. 2000, 92: 121-129. 10.1016/S0006-3207(99)00059-2.View ArticleGoogle Scholar
- Giri A, Lakshmi Narasu M: Production of podophyllotoxin from Podophyllum hexandrum: a potential natural product for clinically useful anticancer drugs. Cytotechnology. 2000, 34: 17-26. 10.1023/A:1008138230896.PubMed CentralPubMedView ArticleGoogle Scholar
- You Y: Podophyllotoxin derivatives: current synthetic approaches for new anticancer agents. Curr Pharm Des. 2005, 11: 1695-1717. 10.2174/1381612053764724.PubMedView ArticleGoogle Scholar
- Kim E, Donohue K: Local adaptation and plasticity of Erysimum capitatum to altitude: its implications for responses to climate change. J Ecol. 2013, 101: 796-805. 10.1111/1365-2745.12077.View ArticleGoogle Scholar
- Kushwaha R, Pandey S, Chanda S, Bhattacharya A, Ahuja PS: Temperature dependent growth and emergence of functional leaves: an adaptive mechanism in the seedlings of the western Himalayan plant Podophyllum hexandrum. J Plant Res. 2008, 121: 299-309. 10.1007/s10265-008-0149-9.PubMedView ArticleGoogle Scholar
- Singh A, Purohit AN: Light and temperature effects on physiological reactions on alpine and temperate populations of Podophyllum hexandrum Royle. J Herbs Spices Med Plants. 1997, 5: 57-66.View ArticleGoogle Scholar
- Schrader J, Moyle R, Bhalerao R, Hertzberg M, Lundeberg J, Nilsson P, Bhalerao RP: Cambial meristem dormancy in trees involves extensive remodelling of the transcriptome. Plant J. 2004, 40: 173-187. 10.1111/j.1365-313X.2004.02199.x.PubMedView ArticleGoogle Scholar
- Derory J, Léger P, Garcia V, Schaeffer J, Hauser MT, Salin F, Luschnig C, Plomion C, Glössl J, Kremer A: Transcriptome analysis of bud burst in sessile oak (Quercus petraea). New Phytol. 2006, 170: 723-738. 10.1111/j.1469-8137.2006.01721.x.PubMedView ArticleGoogle Scholar
- Paul A, Kumar S: Responses to winter dormancy, temperature, and plant hormones share gene networks. Funct Integr Genomics. 2011, 11: 659-664. 10.1007/s10142-011-0233-4.PubMedView ArticleGoogle Scholar
- Gahlan P, Singh HR, Shankar R, Sharma N, Kumari A, Chawla V, Ahuja PS, Kumar S: De novo sequencing and characterization of Picrorhiza kurrooa transcriptome at two temperatures showed major transcriptome adjustments. BMC Genomics. 2012, 13: 126-10.1186/1471-2164-13-126.PubMed CentralPubMedView ArticleGoogle Scholar
- Venkataraman K: India’s Biodiversity Act 2002 and its role in conservation. Trop Ecol. 2009, 50: 23-30.Google Scholar
- Bhattacharyya D, Sinha R, Ghanta S, Chakraborty A, Hazra S, Chattopadhyay S: Proteins differentially expressed in elicited cell suspension culture of Podophyllum hexandrum with enhanced podophyllotoxin content. Proteome Sci. 2012, 10: 34-10.1186/1477-5956-10-34.PubMed CentralPubMedView ArticleGoogle Scholar
- Sultan P, Shawl AS, Abdellah AA, Ramteke PW: Isolation, Characterization and Comparative Study on Podophyllotoxin and Related Glycosides of Podophyllum heaxandrum. Curr Res J Biol Sci. 2010, 2: 345-351.Google Scholar
- Marques JV, Kim KW, Lee C, Costa MA, May GD, Crow JA, Davin LB, Lewis NG: Next generation sequencing in predicting gene function in podophyllotoxin biosynthesis. J Biol Chem. 2013, 288: 466-479. 10.1074/jbc.M112.400689.PubMed CentralPubMedView ArticleGoogle Scholar
- Compound extraction of Podophyllum peltatum. [http://uic.edu/pharmacy/MedPlTranscriptome/d_podo_p_3.html]
- Muoki RC, Paul A, Kumari A, Singh K, Kumar S: An improved protocol for the isolation of RNA from roots of tea (Camellia sinensis (L.) O. Kuntze). Mol Biotechnol. 2012, 52: 82-88. 10.1007/s12033-011-9476-5.PubMedView ArticleGoogle Scholar
- GenScript real-time PCR primer design. [https://www.genscript.com/ssl-bin/app/primer]
- Pfaffl WP, Horgan GW, Dempfle L: Relative expression software tool (REST©) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002, 30: e36-10.1093/nar/30.9.e36.PubMed CentralPubMedView ArticleGoogle Scholar
- Li R, Zhu H, Ruan J, Qian W, Fang X, Shi Z, Li Y, Li S, Shan G, Kristiansen K, Li S, Yang H, Wang J, Wang J: De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2009, 20: 265-272.PubMedView ArticleGoogle Scholar
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22: 1658-1659. 10.1093/bioinformatics/btl158.PubMedView ArticleGoogle Scholar
- UniProt Consortium: Reorganizing the protein space at the Universal Protein Resource (UniProt). Nucleic Acids Res. 2011, 40 (Database issue): D71-D75.Google Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralPubMedView ArticleGoogle Scholar
- Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999, 27: 29-34. 10.1093/nar/27.1.29.PubMed CentralPubMedView ArticleGoogle Scholar
- Webb EC: International Union of Biochemistry and Molecular Biology. Nomenclature Committee. Enzyme nomenclature: Recommendations of the Nomenclature Committee of the International Union of Biochemistry and Molecular Biology on the nomenclature and classification of enzymes. 1992, San Diego: Academic Press, 862-Google Scholar
- Schmid R, Blaxter ML: annot8r: GO, EC and KEGG annotation of EST datasets. BMC Bioinformatics. 2008, 9: 180-10.1186/1471-2105-9-180.PubMed CentralPubMedView ArticleGoogle Scholar
- Pérez-Rodríguez P, Riaño-Pachón DM, Corrêa LG, Rensing SA, Kersten B, Mueller-Roeber B: PlnTFDB: updated content and new features of the plant transcription factor database. Nucleic Acids Res. 2010, 40 (Database issue): D822-D827.View ArticleGoogle Scholar
- Ruepp A, Zollner A, Maier D, Albermann K, Hani J, Mokrejs M, Tetko I, Güldener U, Mannhaupt G, Münsterkötter M, Mewes HW: The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucleic Acids Res. 2004, 32: 5539-5545. 10.1093/nar/gkh894.PubMed CentralPubMedView ArticleGoogle Scholar
- Lamesch P, Berardini TZ, Li D, Swarbreck D, Wilks C, Sasidharan R, Muller R, Dreher K, Alexander DL, Garcia-Hernandez M, Karthikeyan AS, Lee CH, Nelson WD, Ploetz L, Singh S, Wensel A, Huala E: The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 2012, 40 (Database issue): D1202-D1210.PubMed CentralPubMedView ArticleGoogle Scholar
- Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thanki N, Yamashita RA, Zhang D, Zhang N, Zheng C, Bryant SH: CDD: a Conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39: 225-229. 10.1093/nar/gkq769.View ArticleGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-seq. Bioinformatics. 2009, 25: 1105-1111. 10.1093/bioinformatics/btp120.PubMed CentralPubMedView ArticleGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28: 511-515. 10.1038/nbt.1621.PubMed CentralPubMedView ArticleGoogle Scholar
- MISA Program. [http://pgrc.ipk-gatersleben.de/misa]
- Zhang P, Dreher K, Karthikeyan A, Chi A, Pujar A, Caspi R, Karp P, Kirkup V, Latendresse M, Lee C, Mueller LA, Muller R, Rhee SY: Creation of a genome-wide metabolic pathway database for Populus trichocarpa using a new approach for reconstruction and curation of metabolic pathways for plants. Plant Physiol. 2010, 153: 1479-1491. 10.1104/pp.110.157396.PubMed CentralPubMedView ArticleGoogle Scholar
- Salick J, Fangb Z, Byg A: Eastern Himalayan alpine plant ecology, Tibetan ethnobotany, and climate change. Global Environ Chang. 2009, 19: 147-155. 10.1016/j.gloenvcha.2009.01.008.View ArticleGoogle Scholar
- Zobayed SM, Afreen F, Kozai T: Temperature stress can alter the photosynthetic efficiency and secondary metabolite concentrations in St. John’s wort. Plant Physiol Biochem. 2005, 43: 977-984. 10.1016/j.plaphy.2005.07.013.PubMedView ArticleGoogle Scholar
- Alam MA, Gulati P, Aswini KG, Gyan PM, Pradeep KN: Assessment of genetic diversity among Podophyllum hexandrum genotypes of northwestern Himalayan region for podophyllotoxin production. Indian J Biotechnol. 2009, 8: 391-399.Google Scholar
- Chandra S: Effect of altitude on energy exchange characteristics of some alpine medicinal crops from central Himalayas. J Agron Crop Sci. 2004, 190: 13-20. 10.1046/j.0931-2250.2003.00064.x.View ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1016/S0022-2836(05)80360-2.PubMedView ArticleGoogle Scholar
- Vinogradov AE: DNA helix: the importance of being GC-rich. Nucleic Acids Res. 2003, 31: 1838-1844. 10.1093/nar/gkg296.PubMed CentralPubMedView ArticleGoogle Scholar
- Nestor CE, Monckton DG: Correlation of inter-locus polyglutamine toxicity with CAG•CTG triplet repeat expandability and flanking genomic DNA GC content. PLoS One. 2011, 6: e28260-10.1371/journal.pone.0028260.PubMed CentralPubMedView ArticleGoogle Scholar
- Holland JB, Helland SJ, Sharopova N, Rhyne DC: Polymorphism of PCR-based markers targeting exons, introns, promoter regions, and SSRs in maize and introns and repeat sequences in oat. Genome. 2001, 44: 1065-1076. 10.1139/g01-110.PubMedView ArticleGoogle Scholar
- Musashi M, Ota S, Shiroshita N: The role of protein kinase C isoforms in cell proliferation and apoptosis. Int J Hematol. 2000, 72: 12-19.PubMedGoogle Scholar
- Annadurai RS, Jayakumar V, Mugasimangalam RC, Katta MA, Anand S, Gopinathan S, Sarma SP, Fernandes SJ, Mullapudi N, Murugesan S, Rao SN: Next generation sequencing and de novo transcriptome analysis of Costus pictus D. Don, a non-model plant with potent anti-diabetic properties. BMC Genomics. 2012, 13: 663-10.1186/1471-2164-13-663.PubMed CentralPubMedView ArticleGoogle Scholar
- Smith IK, Polle A, Rennenberg H: Glutathione. Stress Responses in Plants: Adaptation and Acclimation Mechanisms. Edited by: Alscher RG, Cumming JR. 1992, New York: Wiley Liss Inc, 201-215.Google Scholar
- Vu JCV, Allen LH, Gesch RW: Up-regulation of photosynthesis and sucrose metabolism enzymes in young expanding leaves of sugarcane under elevated growth CO2. Plant Sci. 2006, 171: 123-131. 10.1016/j.plantsci.2006.03.003.View ArticleGoogle Scholar
- An D, Yang J, Zhang P: Transcriptome profiling of low temperature treated cassava apical shoots showed dynamic responses of tropical plant to cold stress. BMC Genomics. 2012, 13: 64-10.1186/1471-2164-13-64.PubMed CentralPubMedView ArticleGoogle Scholar
- Hotton SK, Callis J: Regulation of cullin RING ligases. Annu Rev Plant Biol. 2008, 59: 467-489. 10.1146/annurev.arplant.58.032806.104011.PubMedView ArticleGoogle Scholar
- TAIR. [https://www.arabidopsis.org/servlets/TairObject?name=AT3G12670&type=locus]
- Henriksson E, Olsson AS, Johannesson H, Johansson H, Hanson J, Engström P, Söderman E: Homeodomain leucine zipper class I genes in Arabidopsis. Expression patterns and phylogenetic relationships. Plant Physiol. 2005, 139: 509-518. 10.1104/pp.105.063461.PubMed CentralPubMedView ArticleGoogle Scholar
- Ha CM, Jun JH, Nam HG, Fletcher JC: BLADE-ON-PETIOLE1 encodes a BTB/POZ domain protein required for leaf morphogenesis in Arabidopsis thaliana. Plant Cell Physiol. 2004, 45: 1361-1370. 10.1093/pcp/pch201.PubMedView ArticleGoogle Scholar
- Jofuku KD, den Boer BG, Van Montagu M, Okamuro JK: Control of Arabidopsis flower and seed development by the homeotic gene APETALA2. Plant Cell. 1994, 6: 1211-1225. 10.1105/tpc.6.9.1211.PubMed CentralPubMedView ArticleGoogle Scholar
- Samach A, Onouchi H, Gold SE, Ditta GS, Schwarz-Sommer Z, Yanofsky MF, Coupland G: Distinct roles of CONSTANS target genes in reproductive development of Arabidopsis. Science. 2000, 288: 1613-1616. 10.1126/science.288.5471.1613.PubMedView ArticleGoogle Scholar
- Kim WC, Reca IB, Kim Y, Park S, Thomashow MF, Keegstra K, Han KH: Transcription factors that directly regulate the expression of CSLA9 encoding mannan synthase in Arabidopsis thaliana. Plant Mol Biol. 2014, 84: 577-587. 10.1007/s11103-013-0154-9.PubMedView ArticleGoogle Scholar
- Apuya NR, Yadegari R, Fischer RL, Harada JJ, Zimmerman JL, Goldberg RB: The Arabidopsis embryo mutant schlepperless has a defect in the chaperonin-60alpha gene. Plant Physiol. 2001, 126: 717-730. 10.1104/pp.126.2.717.PubMed CentralPubMedView ArticleGoogle Scholar
- Davletova S, Schlauch K, Coutu J, Mittler R: The zinc-finger protein Zat12 plays a central role in reactive oxygen and abiotic stress signaling in Arabidopsis. Plant Physiol. 2005, 139: 847-856. 10.1104/pp.105.068254.PubMed CentralPubMedView ArticleGoogle Scholar
- Polidoros AN, Mylona PV, Pasentsis K, Scandalios JG, Tsaftaris AS: The maize alternative oxidase 1a (Aox1a) gene is regulated by signals related to oxidative stress. Redox Rep. 2005, 10: 71-78. 10.1179/135100005X21688.PubMedView ArticleGoogle Scholar
- Diédhiou CJ, Popova OV, Dietz KJ, Golldack D: The SNF1-type serine-threonine protein kinase SAPK4 regulates stress-responsive gene expression in rice. BMC Plant Biol. 2008, 8: 49-10.1186/1471-2229-8-49.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang W, Vinocur B, Shoseyov O, Altman A: Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 2004, 9: 244-252. 10.1016/j.tplants.2004.03.006.PubMedView ArticleGoogle Scholar
- Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K: AP2/ERF family transcription factors in plant abiotic stress responses. Biochim Biophys Acta. 2012, 1819: 86-96. 10.1016/j.bbagrm.2011.08.004.PubMedView ArticleGoogle Scholar
- Yi Y, Guerinot ML: Genetic evidence that induction of root Fe(III) chelate reductase activity is necessary for iron uptake under iron deficiency. Plant J. 1996, 10: 835-844. 10.1046/j.1365-313X.1996.10050835.x.PubMedView ArticleGoogle Scholar
- Sun S, Yu JP, Chen F, Zhao TJ, Fang XH, Li YQ, Sui SF: TINY, a dehydration-responsive element (DRE)-binding protein-like transcription factor connecting the DRE- and ethylene-responsive element-mediated signaling pathways in Arabidopsis. J Biol Chem. 2008, 283: 6261-6271. 10.1074/jbc.M706800200.PubMedView ArticleGoogle Scholar
- Wang J, Chen L, Huang S, Liu J, Ren X, Tian X, Qiao J, Zhang W: RNA-seq based identification and mutant validation of gene targets related to ethanol resistance in cyanobacterial Synechocystis sp. PCC 6803. Biotechnol Biofuels. 2012, 5: 89-10.1186/1754-6834-5-89.PubMed CentralPubMedView ArticleGoogle Scholar
- Xu L, Zhu L, Tu L, Liu L, Yuan D, Jin L, Long L, Zhang X: Lignin metabolism has a central role in the resistance of cotton to the wilt fungus Verticillium dahliae as revealed by RNA-Seq-dependent transcriptional analysis and histochemistry. J Exp Bot. 2011, 62: 5607-5621. 10.1093/jxb/err245.PubMed CentralPubMedView ArticleGoogle Scholar
- Heim MA, Jakoby M, Werber M, Martin C, Weisshaar B, Bailey PC: The basic helix-loop-helix transcription factor family in plants: a genome-wide study of protein structure and functional diversity. Mol Biol Evol. 2003, 20: 735-747. 10.1093/molbev/msg088.PubMedView ArticleGoogle Scholar
- Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L: MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010, 15: 573-581. 10.1016/j.tplants.2010.06.005.PubMedView ArticleGoogle Scholar
- Bienz M: The PHD finger, a nuclear protein-interaction domain. Trends Biochem Sci. 2006, 31: 35-40. 10.1016/j.tibs.2005.11.001.PubMedView ArticleGoogle Scholar
- Tran LS, Nakashima K, Sakuma Y, Simpson SD, Fujita Y, Maruyama K, Fujita M, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Isolation and functional analysis of Arabidopsis stress-inducible NAC transcription factors that bind to a drought-responsive cis-element in the early responsive to dehydration stress 1 promoter. Plant Cell. 2004, 16: 2481-2498. 10.1105/tpc.104.022699.PubMed CentralPubMedView ArticleGoogle Scholar
- Englbrecht CC, Schoof H, Böhm S: Conservation, diversification and expansion of C2H2 zinc finger proteins in the Arabidopsis thaliana genome. BMC Genomics. 2004, 5: 39-10.1186/1471-2164-5-39.PubMed CentralPubMedView ArticleGoogle Scholar
- Jackson DE, Dewick PM: Biosynthesis of Podophyllum lignans-II. Interconversion of aryltetralin lignans in Podophyllum hexandrum. Phytochemistry. 1984, 23: 1037-1042. 10.1016/S0031-9422(00)82604-7.View ArticleGoogle Scholar
- Hoffmann L, Maury S, Martz F, Geoffroy P, Legrand M: Purification, cloning and properties of an acyltransferase controlling shikimate and quinate ester intermediates in phenylpropanoid metabolism. J Biol Chem. 2003, 278: 95-103.PubMedView ArticleGoogle Scholar
- Schoch G, Goepfert S, Morant M, Hehn A, Meyer D, Ullmann P, Werck-Reichhart D: CYP98A3 from Arabidopsis thaliana is a 3′-hydroxylase of phenolic esters, a missing link in the phenylpropanoid pathway. J Biol Chem. 2001, 276: 36566-36574. 10.1074/jbc.M104047200.PubMedView ArticleGoogle Scholar
- Fraser CM, Chapple C: The phenylpropanoid pathway in Arabidopsis. Arabidopsis Book. 2011, 9: e0152-PubMed CentralPubMedView ArticleGoogle Scholar
- Davin LB, Wang HB, Crowell AL, Bedgar DL, Martin DM, Sarkanen S, Lewis NG: Stereoselective bimolecular phenoxy radical coupling by an auxiliary (dirigent) protein without an active center. Science. 1997, 275: 362-366. 10.1126/science.275.5298.362.PubMedView ArticleGoogle Scholar
- Katayama T, Davin LB, Lewis NG: An extraordinary accumulation of (-)-pinoresinol in cell-free extracts of Forsythia intermedia: evidence for enantiospecific reduction of (+)-pinoresinol. Phytochemistry. 1992, 31: 3875-3881. 10.1016/S0031-9422(00)97545-9.PubMedView ArticleGoogle Scholar
- Katayama T, Davin LB, Chu A, Lewis NG: Novel benzylic ether reductions in lignan biogenesis in Forsythia intermedia. Phytochemistry. 1993, 33: 581-591. 10.1016/0031-9422(93)85452-W.View ArticleGoogle Scholar
- Umezawa T, Davin LB, Lewis NG: Formation of lignans, (-)-secoisolariciresinol and (-)-matairesinol with Forsythia intermedia cell-free extracts. J Biol Chem. 1991, 266: 10210-10217.PubMedGoogle Scholar
- Sakakibara N, Suzuki S, Umezawa T, Shimada M: Biosynthesis of yatein in Anthriscus sylvestris. Org Biomol Chem. 2003, 1: 2474-2485. 10.1039/b304411d.PubMedView ArticleGoogle Scholar
- Broomhead AJ, Rahman MA, Dewick PM, Jackson DE, Lucas JA: Matairesinol as precursor of Podophyllum lignans. Phytochemistry. 1991, 30: 1489-1492. 10.1016/0031-9422(91)84194-W.View ArticleGoogle Scholar
- Kuhlmann S, Kranz K, Lücking B: Aspects of cytotoxic lignan biosynthesis in suspension cultures of Linum nodiflorum. Phytochem Rev. 2002, 1: 37-43. 10.1023/A:1015876001812.View ArticleGoogle Scholar
- Kamil WM, Dewick PM: Biosynthesis of the lignans α- and β-peltatin. Phytochemistry. 1986, 25: 2089-2092. 10.1016/0031-9422(86)80071-1.View ArticleGoogle Scholar
- van Uden W, Bouma AS, Bracht Waker JF, Middel O, Wichers HJ, Waard PD, Woerdenbag HJ, Kellogg RM, Pras N: The production of podophyllotoxin and its 5-methoxy derivative through bioconversion of cyclodextrin complexed deoxy-podophyllotoxin by plant cell cultures. Plant Cell Tissue Organ Cult. 1995, 42: 73-79. 10.1007/BF00037684.View ArticleGoogle Scholar
- van Uden W, Bos JA, Boeke GM, Woerdenbag HJ, Pras N: The large-scale isolation of deoxypodophyllotoxin from rhizomes of Anthriscus sylvestris followed by its bioconversion into 5-methoxypodophyllotoxin β-d-glucoside by cell cultures of Linun flavum. J Natural Prod. 1997, 60: 401-403. 10.1021/np960748o.View ArticleGoogle Scholar
- Yousefzadi M, Sharifi M, Behmanesh M, Moyano E, Bonfill M, Cusido RM, Palazon J: Podophyllotoxin: current approaches to its biotechnological production and future challenges. Eng Life Sci. 2010, 10: 281-292. 10.1002/elsc.201000027.View ArticleGoogle Scholar
- Singh K, Rani A, Paul A, Dutt S, Joshi R, Gulati A, Ahuja PS, Kumar S: Differential display mediated cloning of anthocyanidin reductase gene from tea (Camellia sinensis) and its relationship with the concentration of epicatechins. Tree Physiol. 2009, 29: 837-846. 10.1093/treephys/tpp022.PubMedView ArticleGoogle Scholar
- Kawoosa T, Singh H, Kumar A, Sharma SK, Devi K, Dutt S, Vats SK, Sharma M, Ahuja PS, Kumar S: Light and temperature regulated terpene biosynthesis: hepatoprotective monoterpene picroside accumulation in Picrorhiza kurrooa. Funct Integr Genomics. 2010, 10: 393-404. 10.1007/s10142-009-0152-9.PubMedView ArticleGoogle Scholar
- Kumar H, Singh K, Kumar S: 2C-methyl-d-erythritol-2,4-cyclodiphosphate synthase from Stevia rebaudiana Bertoni is a functional gene. Mol Biol Rep. 2012, 12: 10971-10978.View ArticleGoogle Scholar
- Singh RS, Gara RK, Bhardwaj PK, Kaachra A, Malik S, Kumar R, Sharma M, Ahuja PS, Kumar S: Expression of 3-hydroxy-3-methylglutaryl-CoA reductase, p-hydroxybenzoate-m-geranyltransferase and genes of phenylpropanoid pathway exhibits positive correlation with shikonins content in arnebia [Arnebia euchroma (Royle) Johnston]. BMC Mol Biol. 2010, 11: 88-10.1186/1471-2199-11-88.PubMed CentralPubMedView ArticleGoogle Scholar
- Schuler MA: Plant cytochrome P450 monooxygenases. Crit Rev Plant Sci. 1996, 15: 235-284. 10.1080/07352689609701942.View ArticleGoogle Scholar
- Federolf K, Alfermann AW, Fuss E: Aryltetralin-lignan formation in two different cell suspension cultures of Linum album: Deoxypodophyllotoxin 6-hydroxylase, a key-enzyme in the formation of 6-methoxypodophyllotoxin. Phytochemistry. 2007, 68: 1397-1406. 10.1016/j.phytochem.2007.02.031.PubMedView ArticleGoogle Scholar
- Lam KC, Ibrahim RK, Behdad B, Dayanandan S: Structure, function, and evolution of plant O-methyltransferases. Genome. 2007, 50: 1001-1013. 10.1139/G07-077.PubMedView ArticleGoogle Scholar
- Lewis NG, Yamamoto E: Lignin: occurrence, biogenesis and biodegradation. Annu Rev Plant Physiol Plant Mol Biol. 1990, 41: 455-496. 10.1146/annurev.pp.41.060190.002323.PubMedView ArticleGoogle Scholar
- Boerjan W, Ralph J, Baucher M: Lignin biosynthesis. Annu Rev Plant Biol. 2003, 54: 519-546. 10.1146/annurev.arplant.54.031902.134938.PubMedView ArticleGoogle Scholar
- Zhong R, Morrison WH, Himmelsbach DS, Poole FL, Ye ZH: Essential role of caffeoyl coenzyme A O-methyltransferase in lignin biosynthesis in woody poplar plants. Plant Physiol. 2000, 124: 563-578. 10.1104/pp.124.2.563.PubMed CentralPubMedView ArticleGoogle Scholar
- Ross J, Li Y, Lim E, Bowles DJ: Higher plant glycosyltransferases. Genome Biol. 2001, 2: REVIEWS3004-PubMed CentralPubMedView ArticleGoogle Scholar
- Vogt T, Jones P: Glycosyltransferases in plant natural product synthesis: characterization of a supergene family. Trends Plant Sci. 2000, 5: 380-386. 10.1016/S1360-1385(00)01720-9.PubMedView ArticleGoogle Scholar
- Wang J, Hou B: Glycosyltransferases: key players involved in the modification of plant secondary metabolites. Front Biol China. 2009, 4: 39-46. 10.1007/s11515-008-0111-1.View ArticleGoogle Scholar
- Berim A, Ebel R, Schneider B, Petersen M: UDP-glucose:(6-methoxy) podophyllotoxin 7-O-glucosyltransferase from suspension cultures of Linum nodiflorum. Phytochemistry. 2008, 69: 374-381. 10.1016/j.phytochem.2007.07.030.PubMedView ArticleGoogle Scholar
- Jha A, Mehra M, Shankar R: The regulatory epicenter of miRNAs. J Biosci. 2011, 36: 621-638. 10.1007/s12038-011-9109-y.PubMedView ArticleGoogle Scholar
- Jha A, Shankar R: miReader: Discovering novel miRNAs in species without sequenced genome. PLoS One. 2013, 8: e66857-10.1371/journal.pone.0066857.PubMed CentralPubMedView ArticleGoogle 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 credited. 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.