De novo assembly of Euphorbia fischeriana root transcriptome identifies prostratin pathway related genes
- Roberto A Barrero†1,
- Brett Chapman†1,
- Yanfang Yang2,
- Paula Moolhuijzen1,
- Gabriel Keeble-Gagnère1,
- Nan Zhang2,
- Qi Tang2, 3,
- Matthew I Bellgard1Email author and
- Deyou Qiu2Email author
© Barrero et al; licensee BioMed Central Ltd. 2011
Received: 16 May 2011
Accepted: 13 December 2011
Published: 13 December 2011
Euphorbia fischeriana is an important medicinal plant found in Northeast China. The plant roots contain many medicinal compounds including 12-deoxyphorbol-13-acetate, commonly known as prostratin that is a phorbol ester from the tigliane diterpene series. Prostratin is a protein kinase C activator and is effective in the treatment of Human Immunodeficiency Virus (HIV) by acting as a latent HIV activator. Latent HIV is currently the biggest limitation for viral eradication. The aim of this study was to sequence, assemble and annotate the E. fischeriana transcriptome to better understand the potential biochemical pathways leading to the synthesis of prostratin and other related diterpene compounds.
In this study we conducted a high throughput RNA-seq approach to sequence the root transcriptome of E. fischeriana. We assembled 18,180 transcripts, of these the majority encoded protein-coding genes and only 17 transcripts corresponded to known RNA genes. Interestingly, we identified 5,956 protein-coding transcripts with high similarity (> = 75%) to Ricinus communis, a close relative to E. fischeriana. We also evaluated the conservation of E. fischeriana genes against EST datasets from the Euphorbeacea family, which included R. communis, Hevea brasiliensis and Euphorbia esula. We identified a core set of 1,145 gene clusters conserved in all four species and 1,487 E. fischeriana paralogous genes. Furthermore, we screened E. fischeriana transcripts against an in-house reference database for genes implicated in the biosynthesis of upstream precursors to prostratin. This identified 24 and 9 candidate transcripts involved in the terpenoid and diterpenoid biosyntehsis pathways, respectively. The majority of the candidate genes in these pathways presented relatively low expression levels except for 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase (HDS) and isopentenyl diphosphate/dimethylallyl diphosphate synthase (IDS), which are required for multiple downstream pathways including synthesis of casbene, a proposed precursor to prostratin.
The resources generated in this study provide new insights into the upstream pathways to the synthesis of prostratin and will likely facilitate functional studies aiming to produce larger quantities of this compound for HIV research and/or treatment of patients.
Euphorbia fischeriana is a perennial herbaceous flowering plant, located in northeast of China. The dried plant roots of E. fischeriana and other Euphorbia species are called "Lang-Du", which means extremely toxic in Chinese. The roots have often been used in traditional Chinese medicine to treat a wide range of ailments, such as edema, ascites, and cancer . Previous studies of the plant have revealed it mainly contains diterpenoids, triterpenoids and steroids . One particular tigliane diterpene and phorbol ester found in the roots of E. fischeriana is 12-deoxyphorbol-13-acetate, commonly known as prostratin.
Prostratin was first discovered during ethno botanical studies carried out on the island of Savai'i, Samoa, where native healers used the bark of Homalanthus nutans, a small rain-forest tree from the Euphorbiaceae family, called "mamala" to treat hepatitis infection [3, 4]. Prostratin has been shown in previous studies to be a protein kinase C activator  and effective in the treatment of HIV infected patients, by acting as a latent HIV1 activator, enabling anti-viral drugs to eliminate a large proportion of viral reservoirs in the body [6, 7]. Latent HIV is one of the biggest obstacles HIV researchers face in regards to the treatment of patients . Prostratin has also been found to be an anti-tumour promoter, unlike other diterpene compounds that are known to promote tumorigenesis [9, 10].
Very little is known about the biosynthetic pathway of prostratin, although a lot is known about its biochemical structure. The main pathways which are likely to play a role in the synthesis of prostratin and which were chosen for this study were the Terpenoid Backbone Biosynthesis (TBB) pathway (KEGG map00900) and the Diterpenoid Biosynthesis (DB) pathway (KEGG map00904). The TBB pathway is an important pathway for the synthesis of terpenoid or isoprene (5-carbon) compounds that are the building blocks for many important complex compounds synthesised further down the pathway (down stream) in other important biosynthetic pathways. The TBB pathway consists of two separate parallel biosynthetic pathways, one corresponding to the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway that takes place in the plastids, and another comprising the mevalonic acid (MVA) pathway that occurs in the cytosol. Both pathways lead to the synthesis of terpenoid building blocks: isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP). These terpenoid compounds are then used to generate compounds geranyl diphosphate (GPP), farsenyl diphosphate (FPP) and geranylgeranyl diphosphate (GGPP), which are all precursor compounds to pathways further down stream, such as the DB pathway. The DB pathway is an important pathway, which begins with GGPP and leads to the synthesis of many important diterpenoids (20-carbon) compounds. The synthesis of one particular diterpenoid compound that is investigated in this study is casbene, which has a very similar skeleton structure to prostratin. It is highly likely that casbene is a precursor to prostratin .
Similar properties to prostratin have been reported for 12-deoxyphorbol 13-phenylacetate (DPP), which has been derived from Euphorbia resinifera and very close relative of E. fischeriana . DPP has very similar structure to prostratin, with the addition of an ester group at C13, concluding that many more phorbol esters may have anti-HIV properties and would be worth investigating in future studies [13, 14].
Cross talk and/or interaction between distinct metabolic pathways are important to determine the potential metabolite profile of a cell under specific conditions. Distinct pathways may share common intermediate compounds for their downstream processing. An associated pathway to the DB pathway is the Kaurenol and the Zeatin Biosynthesis (ZB) pathway (KEGG map00908). Kaurenol is derived from GGPP and serves as a precursor for the synthesis of various gibberellin compounds that act as plant hormones modulating growth and development . The intermediate DMAPP compound of the TBB pathway is utilized to initiate the ZB pathway that leads to the synthesis of Zeatin, a member of the cytokinin family, a class of plant hormones involved in various processes of plant growth and development (KEGG map00908).
In this study we present the results of next generation sequencing, de novo assembly and annotation of the E. fischeriana root transcriptome. Based on literature review and KEGG pathway information we identified candidate genes involved in the synthesis of upstream precursors to prostratin and estimated the expression levels of these enzymes.
Results and Discussion
Sequencing and de novo transcriptome assembly
Summary of Euphorbia fischeriana transcriptome assembly
Total number of mate-pair reads (before trimming)
Total number of read base pairs (bp)
Average read length (before trimming; bp)
Total number of read mate-pairs (after trimming)
Total number of read singletons (after trimming)
Average read length (after trimming; bp)
Total number of ESTs
Total number of EST base pairs (bp)
Average EST length
Total number of transcripts assembled (pre-isoform filtering)
Total number of transcripts assembled (post-isoform filtering)
Average length of all transcripts (bp)
Transcripts with E-value > = 1e-05 against nr
Average length of transcripts (bp)
To determine the best parameters for transcriptome de novo assembly using Oases  multiple k-mers were compared (Additional file 3). Our analysis determined that de novo assembly using a k-mer of 25 provided the best compromise between high and low abundant transcripts (see Methods). We also determined that a minimum k-mer coverage threshold of two is suitable for de novo assembly as this removes the majority of sequencing errors (Additional file 4). Thus, a k-mer of 25 and a k-mer coverage cut-off of two were used to assemble the E. fischeriana root transcriptome using Oases  by combining both high quality Illumina short reads and ESTs in a hybrid assembly approach.
Another feature of Oases is that attempts to cluster assembled transcripts into 'gene clusters' containing two or more putative alternative splicing (AS) isoforms. Our preliminary manual inspection of randomly selected gene clusters revealed that the majority of the predicted AS isoforms corresponded to spurious calls including RNA degradation products, sequence gaps denoted by Ns that were introduced in the scaffolding step and clustering of unrelated sense-antisense transcripts among others (Additional files 5 and 6). These assembly artefacts occur in part due to the extreme variability in coverage depth between genes, isoforms and along each isoform increasing the complexity of the de Brujin graph structure. To remove spurious isoforms from downstream analyses we selected a single transcript from each gene cluster based on various filtering criteria: i) the transcript has the highest Oases confidence score (ratio of nodes in transcript to nodes in gene cluster) that represents the transcripts with the largest number of exons, ii) encodes the longest ORF, iii) corresponds to the longest nucleotide transcript, and iv) in cases where two or more transcripts have the same length then the one with highest sequence coverage was selected. This generated a reference E. fischeriana root transcriptome of 18,180 transcripts (Table 1).
Statistics of functional annotation of transcripts
% Total transcripts
a) Protein-coding genes
Similar to known proteins
Conserved hypothetical proteins
Hypothetical proteins (ORF > = 80aa)
b) Non-coding RNA genes
Putative long ncRNAs
Pseudo tRNA genes
Total number of assembled transcripts
The species distribution showed that the majority of top matches were to Ricinus communis (68%), Populus trichocarpa (18%) and Vitis vinifera (8%) (Figure 2C). The top matches to R. communis were further evaluated and identified that 5,956 transcripts were highly similar (> = 75%) to the E. fischeriana transcriptome. The BLASTx species distribution showed a bias towards R. communis owing to the over representation of this species within the database compared to other species such as Euphorbia esula, a closer relative of E. fischeriana.
We then evaluated the relative expression levels of E. fischeriana transcripts as described in methods and categorized these into three expression ranges 1-125, 125-250 and more than 250 mean coverage expression level. The majority of the transcripts fall into the first category (14,365; 96%), while only 339 (2%) and 261 (2%) were assigned to the 125-250 and 250 plus categories, respectively.
Functional and pathway annotations
Among the 2,428 transcripts assigned to the response to stimulus category, 370 and 257 were assigned to 'defence response' and 'response to cadmium ion', respectively. Interestingly, 191 of the transcripts assigned to defence response were associated with defence against bacteria, which correlates with the known anti-microbial properties of root extracts from Euphorbia species, such as Euphorbia fusiformis .
Under Molecular Function GO the two most abundant categories were 'binding' and 'catalytic activity' accounting for 33.2% and 11.8% of the transcripts, respectively (Figure 3). The 6,034 transcripts assigned to the binding category were classified into various categories including protein binding (2,552; 14%), nucleotide binding (1,934; 10.6%), ATP binding (1,708; 9.4%), metal ion binding (1,562; 8.6%), DNA binding (868; 4.8%) and RNA binding (596; 3.3%). The large fraction of transcripts associated with protein, nucleotide and/or metabolite binding suggest the presence of an intricate interactome network in E. fischeriana roots.
Next we evaluated the expression levels for each GO code by averaging the coverage of all transcripts matching to each GO code (see Methods). We found that the GO Molecular Function codes 'Structural molecular activity' and 'Antioxidant activity' showed the highest overall expression levels with an average mean coverage of 291 and 136, respectively (Figure 3). Interestingly the majority of the GO Biological Process and Cellular Component codes presented an overall expression level close to 50 (Figure 3).
We next determined the overall expression levels for each KEGG pathway as described above for GO codes. The pathways with the highest expression levels were 'Translation' and 'Signaling molecules and interaction' with mean coverage values of 327 and 148, respectively (Figure 4A). The majority of all other pathways displayed similar expression levels with an overall mean coverage close to 50.
Comparison of E. fischeriana transcriptome with related species
Candidate genes upstream of prostratin biosynthesis pathway
Transcripts matching to genes encoding enzymes involved in the TBB, DB and ZB pathways were found through BLASTx nr searches using an E-value threshold of 1e-05 (Figure 6A). We found two transcripts, encoding for geranylgeranyl diphosphate synthase (GGPPS) and 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) contained transcripts on the opposite strands encoding different genes, not classified as TBB genes. These transcripts were further evaluated to determine if these corresponded to Naturally occurring Antisense Transcripts (NATs), but these were not NATs (Additional file 6). Furthermore, this highlighted another artefact produced using the Oases tool when unrelated transcripts may be clustered into the same cluster.
The DB pathway is of particular interest, as one of the downstream products is the phorbol ester prostratin. All diterpenes begin synthesis from GGPP, a 20-carbon isoprenoid diphosphate, derived primarily from the MEP pathway (Figure 6A). GGPP is also the precursor to many other compounds in plants, such as chlorophylls, prenylated proteins and gibberellins [22, 23]. On entry to the DB pathway, GGPP is converted to a variety of diterpenes by diterpene synthases. Only a hand full of diterpene synthases have been identified thus far, including casbene and neocembrene synthases from R. communis and other well known Euphorbiaceae species [24, 25].
Casbene is considered the most likely precursor to prostratin , although other diterpenes with similar structures could possibly go through a further series of structural changes with the same end product, so other possible paths within the DB pathway were considered such as the initial structural changes of GGPP, to ent-copalyl diphosphate and then ent-kaurene. A number of transcripts coding for these enzymes involved in this path were found with BLASTx nr searches and their expression levels were evaluated based on the mean read coverage of these transcripts as described previously.
The expression levels of enzymes in the TBB and DB pathways were compared to those of enzymes in the ZB pathway. We found that the TBB pathway showed an overall higher expression level as compared to other downstream pathways. Interestingly 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase (HDS) and isopentenyl diphosphate/dimethylallyl diphosphate synthase (IDS) from the MEP pathway presented the highest expression levels (Figure 6A). Comparison of the MEP and MVA pathways within the TBB pathway revealed that genes in the MVA pathway had lower overall expression levels, although enzymes such as 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) and mevalonate diphosphate decarboxylase (MDD) had higher relative expression levels. These findings suggest that overall within the TBB pathway, there is a preference for the synthesis of larger amounts of dimethylallyl diphosphate (DMAPP) and isopentenyl diphosphate (IPP) that is required to drive various downstream pathways including diterpene synthesis (Figure 6A).
Real time PCR primers used for expression validations of selected enzymes.
5'-Forward Primer Sequence-3'
5'-Reverse Primer Sequence-3'
There were no transcripts with sequence similarity to geranyl diphosphate synthase (GPPS) identified in this study (Figure 6A and 6B). This may be due to the low expression level of GPPS and the inefficient assembly of poorly expressed genes. GPPS enzyme is critical for the synthesis of geranyl diphosphate (GPP), which is essential for synthesis of farnesyl diphosphate (FPP) via farnesyl diphosphate synthase (FPPS). FPPS was detected in the transcriptome and its expression was found to be relatively low (Figure 6A and 6B). Additionally, mevalonate kinase (MK) transcripts were not found in the transcriptome even though downstream enzyme transcripts were present (Figure 6B). This may be due to similar reasons for the absence of GPPS. As both GPPS and MK enzymes are present in the related H. brasiliensis species and stored in public databases under the accession numbers AB294710 and AB294693, respectively. These findings emphasis the need to increase the amount of sequencing data and/or the evaluation of lower k-mer values to de novo assemble low expressed genes. Overall we identified 26 enzymes involved in terpenoid and diterpenoid biosynthesis, including two casbene synthases that are a valuable resource for further biochemical and functional studies leading to increase the production of prostratin.
The de novo assembly of the E. fischeriana root transcriptome identified 18,180 transcripts, of these 15,191 encoded genes with sequence similarity in other species and 1,487 represent paralogous genes. This study identified 26 transcripts encoding enzymes involved in various pathways upstream of the casbene biosynthesis pathway, which is a proposed precursor for prostratin. Furthermore we revealed the high expression of HDS and IDS enzymes in the TBB pathway. Critically we found a significant higher expression level of the ent-Kaurene oxidase (ent-Kox) and tRNA Dimethylallyltransferase (tRNA-DMAT) enzymes driving the synthesis of kaurenol and cis-zeatin-O-glucoside + UDP, which compete for available GGPP and DMAPP, respectively. DMAPP is essential for the synthesis of GGPP further downstream, while GGPP is directly essential for the synthesis of casbene. The resources generated in this study will likely facilitate further functional studies aiming to increase the production of prostratin, DPP and other phorbol esters of interest for the advancement of HIV research and/or treatment of patients.
Plant material, RNA isolation and deep sequencing
Live plants of E. fischeriana were collected in June 2008 from Jiagedaqi of Hei longjiang province of China. The plants were then grown in the green house of Chinese Academy of Forestry, Beijing. The root of E. fischeriana was washed with tap water and cut into small pieces. The root materials were immediately frozen in liquid nitrogen and were stored at -80°C until further processing.
Total RNA was isolated according to the method described by Chang et al. . After the RNA pellets were dried, RNA was dissolved in 500 μL of RNase free water. Total RNA purity was checked with Agilent 2100 Nano drop machine. The RNA was stored in a -80°C freezer before being sent to the Beijing Genome Institute (BGI) at Shenzhen with dry ice for mRNA purification and cDNA construction.
The library for transcriptome sequencing was constructed with Illumina's kit following manufacturer's protocol. The mRNA was purified from 10 μg of total RNA using oligo (dT) magnetic beads. After purification, the mRNA was fragmented into small pieces using divalent cations under elevated temperatures. The RNA fragments were used for first strand cDNA synthesis with random primers. Second strand cDNA synthesis was done by using DNA polymerase I and RNaseH. The cDNA fragments then went through an end repair process and were ligated to adapters. The products were purified and enriched with PCR before sequencing on the Illumina GAII sequencing platform. Image deconvolution and quality value calculations were performed using the Illumina GA pipeline 1.3.
RNA isolation and EST sequencing
Frozen root samples stored at -80°C were sent to the Beijing Genome Institute (BGI) at Beijing on dry ice. Total RNA was isolated as described above. The RNA was stored in a -80°C freezer until further processing.
Approximately 1 μg of total RNA was used for preparing a cDNA library using the Creator Smart cDNA Library construction kit (Clontech, Mountain View, CA, USA) following manufacture's instructions. The resulting second cDNA strand products were then run on an agarose gel and those with a size between 1-3 kbp were excised and purified using the QIAquick PCR Purification kit (QIAGEN, Shanghai, China) according to the manufacturer's protocol. The products were transformed into DH10B competent cells. Library was checked with a titer of 2 × 105 pfu/mL and a capacity of 1.2 × 106 clones. A total of 2,099 ESTs were sequenced using capillary sequencing. Vector sequences were removed and 1,884 good EST sequences with an average length of 677 bp and a minimum length of 101 bp were submitted to dbEST at GenBank. The assigned accession numbers are the following: [GenBank:JK523124] to [GenBank:JK525007].
We evaluated several assemblers for the de novo assembly of the E. fischeriana root transcriptome, including Oases , Velvet , QSRA , Euler-SR , Edena  and SOAPdenovo . Preliminary assembled contigs by each tool were blasted against NCBI non-redundant protein database. We found that Oases (Velvet version 1.0.12 and Oases version 0.1.15) was the tool with the largest number of database hits (data not shown) and was selected for downstream analyses.
The reads were first trimmed using the adaptive trimming function of a trimming perl script implemented by Nik Joshi at The Bioinformatics Core, UC Davis Genome Centre http://wiki.bioinformatics.ucdavis.edu/index.php/Trim.pl. Additional files 1 and 2 show the results of quality assessment using FastQC prior and after trimming of poor bases and/or removal of poor reads , respectively. To assess the best parameters to use for this assembly, multiple assemblies from k-mer 17 to 47 were compared based on N50, the number of transcripts and the number of gene clusters (Additional file 3). A k-mer of 25 was determined to be the best k-mer, with the highest N50, highest number of transcripts and the highest number of gene clusters. A minimum transcript size of 100 bp was also compared to 300 bp for all assemblies in the comparison. The appropriate k-mer coverage cut off was determined using an R package plotrix (Additional file 4). All assemblies used a minimum k-mer coverage of 2× and a pair-end insert size of 200 bp was used and the assembly was assisted using 1,884 E. fischeriana ESTs. All other parameters were used on default settings. The final assembly used for annotation used a minimum transcript length of 300 bp.
In many instances Oases predicted spurious isoforms and to increase confidence in the assembled sequences a single transcript from each gene cluster was selected based on the following criteria: i) the transcript has the highest Oases confidence score (ratio of nodes in transcript to nodes in a gene cluster) that represents the transcript with the largest number of exons, ii) encodes the longest ORF, iii) corresponds to the longest nucleotide transcript, and iv) in cases where two or more transcripts have the same length then the one with highest sequence coverage is selected. This generated a dataset of 18,180 transcripts, of these 9,883 transcripts were submitted to GenBank/DDBJ/EMBL following their submission guidelines under the project ID 66759 and Locus Tag EFI. The remaining 8,297 sequences contained gaps denoted by 'Ns' that were introduced during the scaffolding step using pair-end short reads with an expected insert size of 200 bp. A fasta file of all 18,180 transcripts is presented in the Additional file 10.
Annotation by peptide sequence was done by searching transcripts against the NCBI non-redundant (nr) peptide database which includes all non-redundant GenBank CDS translations, RefSeq Proteins, PDB, SwissProt, PIR and PRF, excluding environmental samples from WGS projects. The search was conducted using BLASTx with an E-value cut-off of 1e-05 and matching to the top hits.
The assembled transcripts were scanned for the presence of tRNA and rRNA sequences using the programs tRNAscan-SE  and RNAmmer , respectively. The tRNA transcripts were predicted from the original assembly using a k-mer of 25 and a minimum transcript size of 300 bp. To identify additional tRNAs we conducted a new assembly using a shorter k-mer of 17 and a minimum transcript length of 100 bp. Newly assembled transcripts were then screened for tRNAs as described above.
Gene Ontology (GO) and KEGG pathways
The GO and KEGG annotations were carried out using the annotation program Annot8r , which assigned GO and KEGG pathway terms to the transcripts. The program requires a prepared MySQL database and the transcripts in a fasta file. The user progresses through a series of menus, selecting fasta file name, database name, E-value cut off and number of top hits. The assembled transcripts were annotated using an E-value cut off of 1e-20 and the top 5 hits were used for the annotation of each sequence.
Related species comparison analysis
The EST datasets of closely related species, namely E. esula, H. brasiliensis and R. communis were downloaded from the NCBI EST database. Non-redundant datasets were then generated using CD-HIT-EST as previously described . This yielded non-redundant sequence datasets for E. fischeriana (18,179 sequences), E. esula (32,532 sequences), H. brasiliensis (5,457 sequences), and R. communis (19,893 sequences). Sequence similarity comparisons and clustering were performed using tBLASTx in conjunction with OrthoMCL  using a defined E-value cut off of 1e-20.
Expression analysis and prostratin candidate genes
To determine the relative expression levels of E. fischeriana transcripts high quality trimmed short reads were mapped onto these transcripts using the Burrows-Wheeler Aligner (BWA)  and coverage for each nucleotide was determined using SAMtools . The mean coverage for each transcript was then calculated by averaging the coverage for each nucleotide within the transcript. The expression levels of transcripts were then categorized into various expression ranges.
An in-house database of prostratin pathway related candidate genes was created by interrogating the literature and KEGG pathways for genes matching to the TBB, DB and the comparative downstream pathway, the ZB pathway. We then screened E. fischeriana transcripts against this in-house database using BLASTx to identify significant matches to enzymes in the TBB (including MEP and MVA pathways), DB and ZB pathways. The ZB pathway, which has little relevance to the synthesis of prostratin and other diterpenes, was chosen for use as a comparison to the DB pathway, to compare other possible competing downstream pathways. The mean coverage values for all transcripts were plotted to determine the changes in expression levels through the pathways.
RNA isolation, reverse transcription and Real-time PCR
E. fischeriana total RNA was isolated from the roots using Column Plant RNAout kit (TIANDZ, Beijing, China), according to the manufacturer's protocol. RNA was treated with DNase I (Takara, Dalian, China) to remove residual genomic DNA. The concentration of the isolated RNA and the 260/280-absorbance ratio was measured in triplicates with Nanodrop ND-8000 (Thermo, USA). The quality of RNA samples was confirmed by electrophoresis on a 1.2% agarose. Total RNA was reverse-transcribed to cDNA using PrimeScript RT reagent Kit (Perfect Real Time) (Takara, Dalian, China) in a total volume of 10 μl, according to the manufacturer's instruction. About 600 ng of total RNA, 2 μl 5 × PrimeScript buffer (for Real Time), 0.5 μl PrimeScript RT Enzyme Mix I, 0.5 μl Oligo dT Primer (50 μM) and 2 μl Random 6 mers (100 μM) were mixed. The reaction was carried out at 37°C for 15 min and 85°C for 5 s.
Several enzymes from the Terpenoid Biosynthesis pathway (Figure 6) that showed various levels of expression were selected for validation using real time PCR. Forward (Fwd) and reverse (Rev) primers were designed using Primer3 as described previously . Table 3 shows the primers for the selected enzymes and controls. The Real-time PCR assays were conducted in an optional 96-well plate with ABI7500 system (ABI, USA) and a commercial SRBR-Green master mix kit (Takara, Dalian, China), according to the manufacturer's protocol. The reactions were performed with 1 μl cDNA in 20 μl reaction mix containing 10 μl 2 × SYBR Premix Ex Taq and 1.0 μl primers. The conditions were as follows: initial holding at 95°C for 3 min, followed by a two-step program of 95°C for 15 s and 58°C for 33 s for 40 cycles. Each sample was analyzed in three technical replicates and mean Ct values were calculated. Reverse transcriptase negative controls and "no template controls" (without cDNA in PCR) were included.
Human Immunodeficiency Virus
Expressed Sequence Tag
Kyoto Encyclopedia of Genes and Genomes
Open Reading Frame
acetoacetyl-coenzyme A (CoA) thiloase
2-C-methyl-erythritol 4-phosphate cytidyl transferase
1-deoxy-D-xylulose 5-phosphate reductoisomerase
farnesyl diphosphate synthase
geranylgeranyl diphosphate synthase
geranyl diphosphate synthase
3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) reductase
isopentenyl diphosphate isomerase
4-(cytidine 5'-diphospho)-2-C-methyl-D-erythritol kinase
mevalonate diphosphate decarboxylase
isopentenyl diphosphate/dimethylallyl diphosphate synthase
2-C-methyl-D-erythritol 2,4-cyclodiphosphate synthase
1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase
3S-hydroxy-3-methylglutaryl coenzyme A
4-(cytidine 5'-diphospho)-2C-methyl-D-erythritol 2-phosphate
- ent-KSA ent:
-Kaurene synthase A
- ent-KSB ent:
-Kaurene synthase B
- ent-Kox ent:
- cis-ZOG cis-:
- ent-CPP ent:
- ent-K ent:
Terpenoid Backbone Biosynthesis
We acknowledge the National Collaborative Research Infrastructure (NCRIS) of Australia for providing access to bioinformatics computational pipelines used in this study and The Research Institute of Forestry, Chinese Academy of Forestry for financial support (ZD200902).
- Wang YB, Huang R, Wang HB, Jin HZ, Lou LG, Qin GW: Diterpenoids from the roots of Euphorbia fischeriana. J Nat Prod. 2006, 69 (6): 967-970. 10.1021/np0600088.PubMedView ArticleGoogle Scholar
- Qin GW, Xu RS: Recent advances on bioactive natural products from Chinese medicinal plants. Med Res Rev. 1998, 18 (6): 375-382. 10.1002/(SICI)1098-1128(199811)18:6<375::AID-MED2>3.0.CO;2-8.PubMedView ArticleGoogle Scholar
- Cox PA: Saving the Ethnopharmacological Heritage of Samoa. Journal of Ethnopharmacology. 1993, 38 (2-3): 181-188. 10.1016/0378-8741(93)90015-W.PubMedView ArticleGoogle Scholar
- Cox PA: Ensuring equitable benefits: The Falealupo covenant and the isolation of anti-viral drug prostratin from a Samoan medicinal plant. Pharm Biol. 2001, 39: 33-40.PubMedGoogle Scholar
- Gustafson KR, Cardellina JH, McMahon JB, Gulakowski RJ, Ishitoya J, Szallasi Z, Lewin NE, Blumberg PM, Weislow OS, Beutler JA, Buckheit RW, Cragg GM, Cox PA, Bader JP, Boyd MR: A nonpromoting phorbol from the samoan medicinal plant Homalanthus nutans inhibits cell killing by HIV-1. J Med Chem. 1992, 35 (11): 1978-1986. 10.1021/jm00089a006.PubMedView ArticleGoogle Scholar
- Hezareh M: Prostratin as a new therapeutic agent targeting HIV viral reservoirs. Drug News Perspect. 2005, 18 (8): 496-500. 10.1358/dnp.2005.18.8.944543.PubMedView ArticleGoogle Scholar
- Kulkosky J, Culnan DM, Roman J, Dornadula G, Schnell M, Boyd MR, Pomerantz RJ: Prostratin: activation of latent HIV-1 expression suggests a potential inductive adjuvant therapy for HAART. Blood. 2001, 98 (10): 3006-3015. 10.1182/blood.V98.10.3006.PubMedView ArticleGoogle Scholar
- Brown HJ, McBride WH, Zack JA, Sun R: Prostratin and bortezomib are novel inducers of latent Kaposi's sarcoma-associated herpesvirus. Antivir Ther. 2005, 10 (6): 745-751.PubMedGoogle Scholar
- Szallasi Z, Krsmanovic L, Blumberg PM: Nonpromoting 12-deoxyphorbol 13-esters inhibit phorbol 12-myristate 13-acetate induced tumor promotion in CD-1 mouse skin. Cancer Res. 1993, 53 (11): 2507-2512.PubMedGoogle Scholar
- Gulakowski RJ, McMahon JB, Buckheit RW, Gustafson KR, Boyd MR: Antireplicative and anticytopathic activities of prostratin, a non-tumor-promoting phorbol ester, against human immunodeficiency virus (HIV). Antiviral Res. 1997, 33 (2): 87-97. 10.1016/S0166-3542(96)01004-2.PubMedView ArticleGoogle Scholar
- Kirby J, Nishimoto M, Park JG, Withers ST, Nowroozi F, Behrendt D, Rutledge EJ, Fortman JL, Johnson HE, Anderson JV, Keasling JD: Cloning of casbene and neocembrene synthases from Euphorbiaceae plants and expression in Saccharomyces cerevisiae. Phytochemistry. 2010, 71 (13): 1466-1473. 10.1016/j.phytochem.2010.06.001.PubMedView ArticleGoogle Scholar
- Bocklandt S, Blumberg PM, Hamer DH: Activation of latent HIV-1 expression by the potent anti-tumor promoter 12-deoxyphorbol 13-phenylacetate. Antiviral Res. 2003, 59 (2): 89-98. 10.1016/S0166-3542(03)00034-2.PubMedView ArticleGoogle Scholar
- Brooks G, Morrice NA, Ellis C, Aitken A, Evans AT, Evans FJ: Toxic phorbol esters from Chinese tallow stimulate protein kinase C. Toxicon. 1987, 25 (11): 1229-1233. 10.1016/0041-0101(87)90141-3.PubMedView ArticleGoogle Scholar
- Marquez N, Calzado MA, Sanchez-Duffhues G, Perez M, Minassi A, Pagani A, Appendino G, Diaz L, Munoz-Fernandez MA, Munoz E: Differential effects of phorbol-13-monoesters on human immunodeficiency virus reactivation. Biochem Pharmacol. 2008, 75 (6): 1370-1380. 10.1016/j.bcp.2007.12.004.PubMedView ArticleGoogle Scholar
- Plant-Hormones. [http://www.plant-hormones.info/gibberellins.htm]
- Oases Transcriptome Assembler. [http://www.ebi.ac.uk/~zerbino/oases/]
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25 (5): 955-964. 10.1093/nar/25.5.955.PubMed CentralPubMedView ArticleGoogle Scholar
- Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW: RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007, 35 (9): 3100-3108. 10.1093/nar/gkm160.PubMed CentralPubMedView ArticleGoogle 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
- DeLong A, Mockaitis K, Christensen S: Protein phosphorylation in the delivery of and response to auxin signals. Plant Mol Biol. 2002, 49 (3-4): 285-303.PubMedView ArticleGoogle Scholar
- Natarajan D, Britto SJ, Srinivasan K, Nagamurugan N, Mohanasundari C, Perumal G: Anti-bacterial activity of Euphorbia fusiformis--a rare medicinal herb. J Ethnopharmacol. 2005, 102: 4-View ArticleGoogle Scholar
- Gerber E, Hemmerlin A, Hartmann M, Heintz D, Hartmann MA, Mutterer J, Rodriguez-Concepcion M, Boronat A, Van Dorsselaer A, Rohmer M, Crowell DN, Bach TJ: The plastidial 2-C-methyl-D-erythritol 4-phosphate pathway provides the isoprenyl moiety for protein geranylgeranylation in tobacco BY-2 cells. Plant Cell. 2009, 21 (1): 285-300. 10.1105/tpc.108.063248.PubMed CentralPubMedView ArticleGoogle Scholar
- Kasahara H, Hanada A, Kuzuyama T, Takagi M, Kamiya Y, Yamaguchi S: Contribution of the mevalonate and methylerythritol phosphate pathways to the biosynthesis of gibberellins in Arabidopsis. J Biol Chem. 2002, 277 (47): 45188-45194. 10.1074/jbc.M208659200.PubMedView ArticleGoogle Scholar
- Mau CJ, West CA: Cloning of casbene synthase cDNA: evidence for conserved structural features among terpenoid cyclases in plants. Proc Natl Acad Sci USA. 1994, 91 (18): 8497-8501. 10.1073/pnas.91.18.8497.PubMed CentralPubMedView ArticleGoogle Scholar
- Wildung MR, Croteau R: A cDNA clone for taxadiene synthase, the diterpene cyclase that catalyzes the committed step of taxol biosynthesis. J Biol Chem. 1996, 271 (16): 9201-9204. 10.1074/jbc.271.16.9201.PubMedView ArticleGoogle Scholar
- Chang S, Puryear J, Cairney J: A simple and efficient method for isolating RNA from pine trees. Plant Molecular Biology Reporter. 1993, 113-116. 11
- Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18 (5): 821-829. 10.1101/gr.074492.107.PubMed CentralPubMedView ArticleGoogle Scholar
- Bryant DW, Wong WK, Mockler TC: QSRA: a quality-value guided de novo short read assembler. BMC Bioinformatics. 2009, 10: 69-10.1186/1471-2105-10-69.PubMed CentralPubMedView ArticleGoogle Scholar
- Chaisson MJ, Pevzner PA: Short read fragment assembly of bacterial genomes. Genome Res. 2008, 18 (2): 324-330. 10.1101/gr.7088808.PubMed CentralPubMedView ArticleGoogle Scholar
- Hernandez D, Francois P, Farinelli L, Osteras M, Schrenzel J: De novo bacterial genome sequencing: millions of very short reads assembled on a desktop computer. Genome Res. 2008, 18 (5): 802-809. 10.1101/gr.072033.107.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. 2010, 20 (2): 265-272. 10.1101/gr.097261.109.PubMed CentralPubMedView ArticleGoogle Scholar
- Krueger F, Andrews SR, Osborne CS: Large scale loss of data in low-diversity illumina sequencing libraries can be recovered by deferred cluster calling. PLoS One. 6 (1): e16607-
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22 (13): 1658-1659. 10.1093/bioinformatics/btl158.PubMedView ArticleGoogle Scholar
- Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003, 13 (9): 2178-2189. 10.1101/gr.1224503.PubMed CentralPubMedView ArticleGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralPubMedView ArticleGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.PubMed CentralPubMedView ArticleGoogle Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.