Analysis of Dendrobium huoshanense transcriptome unveils putative genes associated with active ingredients synthesis
BMC Genomics volume 19, Article number: 978 (2018)
Dendrobium huoshanense C.Z. Tang et S.J. Cheng is a traditional Chinese herbal medicine with high medicinal value in China. Polysaccharides and alkaloids are its main active ingredients. To understand the difference of main active ingredients in different tissues, we determined the contents of polysaccharides and alkaloids in the roots, stems and leaves of D. huoshanense. In order to explore the reasons for the differences of active ingredients at the level of transcription, we selected roots, stems and leaves of D. huoshanenese for transcriptome sequencing and pathway mining.
The contents of polysaccharides and alkaloids of D. huoshanense were determined and it was found that there were significant differences in different tissues. A total of 716,634,006 clean reads were obtained and 478,361 unigenes were assembled by the Illumina platform sequencing. We identified 1407 carbohydrate-active related unigenes against CAZy database including 447 glycosyltransferase genes (GTs), 818 glycoside hydrolases (GHs), 60 carbohydrate esterases (CEs), 62 carbohydrate-binding modules (CBMs), and 20 polysaccharide lyases (PLs). In the glycosyltransferases (GTs) family, 315 differential expression genes (DEGs) were identified. In total, 124 and 58 DEGs were associated with the biosynthesis of alkaloids in Dh_L vs. Dh_S and Dh_R vs. Dh_L, respectively. A total of 62 DEGs associated with the terpenoid pathway were identified between Dh_R and Dh_S. Five key enzyme genes involved in the terpenoids pathway were identified, and their expression patterns in different tissues was validated using quantitative real-time PCR.
In summary, our study presents a transcriptome profile of D. huoshanense. These data contribute to our deeper relevant researches on active ingredients and provide useful insights into the molecular mechanisms regulating polysaccharides and alkaloids in Dendrobium.
Dendrobium is the second largest genus of the Orchidaceae widely distributed in tropical and subtropical regions of Asia, Oceania, and other areas . Dendrobium huoshanense C.Z. Tang et S.J. Cheng is a perennial epiphytic herb, which belongs to Orchidaceae Dendrobium Sw. The wild D. huoshanense is endangered because of over-exploitation and habitat deterioration. It is distributed only in the Ta-pieh Mountains of China including Huoshan, Jinzhai, Yuexi and Shucheng County in Anhui Province and Yingshan County in Hubei Province. D. huoshanense has been used for nourishing the stomach, promoting the secretion of body fluids as well as treating throat inflammation and enhancing immunity in China . The active medicinal components of D. huoshanense are very complex, including polysaccharides, alkaloids, amino acids, phenols, coumarins, terpenes, flavonoids, benzyl compounds, and several trace mineral elements [3, 4], and its major components are polysaccharides and alkaloids [5, 6].
Previous work focused on structural identification and characterization of polysaccharides and alkaloid components of Dendrobium species. Dendrobium polysaccharides are mainly composed of glucose, galactose, xylose, rhamnose, mannose and other monosaccharides. However, monosaccharide components and contents of polysaccharides differ among Dendrobium species [7, 8]. Polysaccharides are natural macromolecules with anti-tumor, anti-oxidation, anti-aging, antibacterial, antiviral, hypoglycemic, blood lipid metabolism, anti-radiation, anticoagulant and other biological activities [9,10,11]. The polysaccharide HPS-1B23 was isolated from D. huoshanense by chemical methods and the nuclear magnetic resonance (NMR) technique. It was composed of glucose, mannose and galactose, and the molar ratio of each monosaccharide was 31:10:8 . A new Dendrobium polysaccharide DHPD1 with a molecular weight of 3.2 × 103 Da was isolated from the protocorm of D. huoshanense, and the results showed that it was mainly composed of glucose, arabinose and galactose at a molar ratio of 1.023: 0.023: 0.021 and small amounts of mannose and xylose . A water-soluble polysaccharide was isolated and purified from the fresh stems of D. huoshanense, which mainly consisted of glucose, xylose, and galactose in molar ratio of 1.1:1.0:0.5, as well as trace of galacturonic acid . Above all, the main monosaccharides of the Dendrobium species were glucose and mannose. Now, most studies of polysaccharides are related to their structure, composition and bioactivity, but there are few studies that focus on the key enzymes and genes in the biosynthesis of polysaccharide. Therefore, it is important to study the molecular mechanism of polysaccharide synthesis in D. huoshanense.
The contents, structures and pharmacological effects of Dendrobium alkaloids have already been clarified. In 2000, Morita et al. isolated three new Dendrobium-type sesquiterpenoid alkaloids from the whole plants of Dendrobium Snowflake . Through the study of the anti-inflammatory mechanism of total alkaloids from D. nobile, we found that the alkaloid enriched extract from D. nobile attenuated lipopolysaccharide-induced hyperphosphorylation of tau protein in rat’s hippocampus and protected against lipopolysaccharide-induced apoptosis in rat brain . Bioactive components and pharmacological effects play important roles in the study of medicinal plants. Therefore, biosynthetic pathway analysis and key enzyme gene mining have also become the main aim of transcriptome research in medicinal plants. Unveiling the D. officinale transcriptome allowed discovery of alkaloids that belong to the terpenoid indole alkaloid class, and to identify five key enzymes involved in the construction of the backbone of terpene indole alkaloids . It is known that the mevalonate (MVA) pathway and the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway supply the prenyl diphosphates in plants .
So far, several D. officinale transcriptomes [17, 19, 20] and a few genomes [21, 22] have been sequenced. A lot of key enzyme-encoding genes involved in the synthesis and metabolic pathways have been identified. However, there are some differences in the genomes of D. huoshanense and D. officinale, and D. huoshanense that has not been subjected to genome sequencing, so there are few studies on D. huoshanense. Here, we measured the content of polysaccharides and alkaloids in the roots, stems and leaves of D. huoshanense. It was found that the contents of polysaccharides and alkaloids were significantly different in different tissues. In order to explore the causes of these differences, we collected the roots, stems and leaves of D. huoshanense for transcriptome sequencing. Through comparative analysis of different tissues, we obtained a large number of differentially expressed genes, including key genes and transcription factors involved in the synthesis of polysaccharides and alkaloids. Furthermore, we identified and analyzed the key metabolic pathways MVA and MEP involved in the synthesis of alkaloids from D. huoshanense. The expression levels of key genes on the MVA and MEP metabolic pathways were further verified by real-time quantitative PCR. These results provide sufficient data resources and new insights for further study of the molecular mechanisms of the synthesis of polysaccharides and alkaloids from D. huoshanense.
D. huoshanense plants were artificially cultivated and collected from the greenhouse of Anhui Tongjisheng Biotechnology Company, Lu’an, China. Seed germination and protocorm-like bodies growth were cultured on half-strength Murashige and Skoog (MS) medium  adding 6-BA 0.1 mg·L− 1, NAA 0.5 mg·L− 1 and 1% additives (30 g·L− 1 sucrose + 4 g·L− 1 agar + 20% potato) under a 12/12 h light–dark cycle (approx. 30 μmol m− 2·S− 1) at 25 ± 2 °C. After 6 months, the plants were transplanted into pots and placed in the greenhouse at a temperature of 25–27 °C with a light/dark cycle of 12/12 h and 60–70% relative humidity. Roots, stems and leaves were collected from 2-year-old D. huoshanense plants for RNA extraction and determination of polysaccharides and total alkaloids content in March, 2017. In the determination of polysaccharides and total alkaloids, we selected the roots, stems and leaves of five independent plants for determination, that is, five biological replicates. In the transcriptome sequencing, we did three biological replicates of samples from independent plants. All D. huoshanense samples were stored at − 80 °C in an ultra-low temperature freezer.
Determination of polysaccharide and total alkaloid contents
Leaf, stem and root samples were collected from 2-year-old D. huoshanense at maturation stage. The phenol-sulfuric acid method was applied to determine the polysaccharide contents in different tissues. The polysaccharide content was determined by using a glucose standard.
The content of total alkaloid was determined using a method that is described as follows. 0.5 g fresh stems were soaked with 5 ml ammonia water, and they were sealed for 30 min. Then 10 ml of chloroform was added and the weight was recorded. And the mixture was in Soxhlet extractor for 2 h with water bath at 75 °C. Then the mixture was cooled and the chloroform was added to the original weight, and then it was filtered. The lower filtrate (1 ml) was placed in a 100-ml volumetric flask and the volume was increased with chloroform. 10 ml sample solution was put in the 50-ml centrifuge tube and 5 ml potassium hydrogen phthalate buffer solution (pH 4.5) was added. Then 1 ml 0.04% bromocresol green solution was added and it was oscillated for 5 min. Then it was placed for 30 min and 5 ml supernviaould like liquid was added in the test tubes and 1 ml 0.01 M NaOH anhydrous ethanol solution was added. The liquid was fully mixed and the absorbance was measured at 620 nm. The total alkaloid content was calculated by using a Dendrobine standard. There were five biological replicates for each tissue in this experiment.
Total RNA extraction, cDNA library preparation and transcriptome sequencing
Total RNA was extracted using OmniPlant RNA Kit (Cwbio, China) according to the manufacturer’s protocol. The OD260/280 should range from 1.8–2.0 to ensure the purity of the RNA sample. RNA integrity was monitored by agarose gel electrophoresis (1%) and using an Agilent 2100 Bioanalyzer (Agilent Technologies, USA) with RIN number > 7.5. A total of 5 μg RNA was used as an input per sample. Following the manufacturer’s instructions, the samples for transcriptome analysis were prepared using the TruSeqTM RNA Sample Preparation Kit (Illumina, USA). In short, the protocol consisted of the following steps: mRNA was isolated from total RNA using oligo (dT) magnetic beads and cut into short fragments by adding fragmentation buffer. First-strand cDNA was synthesized using random hexamer-primers, taking these short fragments as templates, which were then used to synthesize second-strand cDNA. The products were purified and enriched by PCR to create the final cDNA libraries. The prepared libraries were sequenced on an Illumina Hiseq 2500 platform and 125 bp paired-end reads were generated.
De novo assembly and functional annotation
To get clean reads, low-quality reads and adapter sequences were removed using SeqPrep (https://github.com/jstjohn/SeqPrep) and Sickle (https://github.com/najoshi/sickle). All clean reads were assembled using Trinity software  based on the left.fq and right.fq, with the min_kmer_cov set to 2 and all other parameters set as their defaults. For function annotation, the longest transcript of each gene was defined as the ‘unigene’. Nucleotide sequences of all unigenes were searched against the current version of Nr (NCBI non-redundant protein sequences) , GO (Gene ontology) , KEGG (Kyoto encyclopedia of genes and genomes) databases , and Swiss-Prot (A manually annotated and reviewed protein sequence database) , Pfam 32.0 (Protein family)  and String 10.0 (Search tool for the retrieval of interacting genes)  using BLAST2GO of version 2.5 with a cut-off E-value of 10− 5 . In order to identify genes related to carbohydrate activity, all unigenes were searched against CAZy database using BLAST 2.7.1 with a cut-off E-value of 10− 5.
Identification of differentially expressed genes (DEGs)
Gene expression level of all samples was estimated by mapping clean reads to the Trinity transcripts assembly using RSEM version 1.2.15  with the bowtie2 parameter set at mismatch 0. Differential expression analysis of two samples was performed using the edgeR . The thresholds for significant differential expression were a false discovery rate (FDR) < 0.05 and a |log2(fold change)| of ≥1. The identified DEGs were used for GO and KEGG enrichment analyses, which were performed using the Goatools version 0.5.9 (https://github.com/tanghaibao/Goatools) and KOBAS version 2.0.12 with default settings, respectively [34, 35].
Quantitative real-time PCR (qRT-PCR)
Five genes involved in the MVA and MEP pathways synthesis were analyzed using qRT-PCR. The five key enzyme genes are FPS, HMGR, DXR, HDR and DXS. All gene IDs and the primer sequences are listed in the Additional file 1: Table S1. The remaining total RNA from transcriptome sequencing was used for qRT-PCR analysis. The specific primers for DEGs were designed by Oligo 7 software. The qRT-PCR was performed on an ABI 7500 Real-time PCR system (Applied Biosystems, USA) using SYBR Premix Ex Taq (Takara, Japan) according to the manufacturer’s protocol. The expressions of DXR, HDR, DXS, HMGR and FPS were normalized against Actin and were calculated by the 2−ΔΔCT method . Pearson correlation analysis was used to calculate the consistency of RNA-seq and qRT-PCR data.
Determination of polysaccharides and total alkaloids in three different tissues
Polysaccharide and total alkaloids were determined in three different tissues; leaf, stem and root. The results showed that polysaccharide was mainly concentrated in stems, while total alkaloid was at high concentration in leaves. The highest content of polysaccharide was 32.88% for stems and the highest content of total alkaloid was 0.034% for leaves (Fig. 1). However, the content of polysaccharide and total alkaloids is the lowest in the roots (Fig. 1). These results indicate that the polysaccharides and alkaloids of D. huoshanense have significant differences in different tissues.
Transcriptome sequencing and de novo assembly
In this study, nine cDNA libraries (root cDNA libraries: Dh_R1, Dh_R2, Dh_R3; stem cDNA libraries: Dh_S1, Dh_S2, Dh_S3; leaf cDNA libraries: Dh_L1, Dh_L2, Dh_L3) were constructed. The sequencing raw data were deposited in the NCBI Sequence Read Archive (SRA) database under the accession number SRP122499. A total of 736,904,076 reads constituted the raw data file. After trimming the adapter and low-quality sequences, we obtained 716,634,006 clean reads with a total of 105.58 Gb nucleotides. The base average error rate was 0.01%, and the average Q20 and Q30 values were 97.81 and 93.77%, respectively. The average GC content was 47.15% (Table 1).
A total of 595,635 transcripts generated; among them, the shortest transcript was 201 bp and the longest transcript was 16,740 bp. The average length was 764.73 bp, and the N50 was 1241 bp. In total, 478,361 unigenes were obtained in the range of 201–16,740 bp with an N50 length of 902 bp (Table 1). Of these, 271,521 unigenes were 200–400 bp, 132,811 unigenes were 400–1000 bp, 53,004 unigenes were 1–2 kb and the remaining 21,025 unigenes were > 2 kb (Table 2).
Functional annotation and classification
Of the total 478,361 unigenes, 196,739, 101,132, 108,417, 69,529 and 37,775 had a hit in the NR, Swiss-Prot, KEGG, Pfam and String database, respectively (Table 3). A total of 14,743 unigenes were annotated to five major databases in common (Additional file 2: Figure S3). A total of 91,252 annotated unigenes were grouped into 63 functional groups by using BLAST2GO . Among these groups, 25 groups were involved in ‘biological processes’ (BP), 20 groups in ‘cellular components’ (CC) and 18 groups in ‘molecular functions’ (MF). Metabolic processes, cellular processes and catalytic activities are important life activities in plants. In D. huoshanense, most of the important active ingredients belong to metabolites such as polysaccharides and alkaloids, and their production depends on the catalysis of key enzymes. Therefore, the three top GO terms were metabolic process, cellular process and catalytic activity (Fig. 2). A total of 27,764 unigenes were annotated to 25 groups in the COG database, the three top terms were J (Translation, ribosomal structure and biogenesis), O (Posttranslational modification, protein turnover, chaperones) and R (General function prediction only) (Fig. 3). In total, 108,417 unigenes annotated to 33 pathways in the KEGG database (Fig. 4). All of the unigenes were divided into five branches according to the KEGG metabolic pathway: Cellular Processes (A); Environmental Information Processing (B), Genetic Information Processing (C), Metabolism (D), and Organismal Systems (E) (Fig. 4). For KEGG, 28,139 unigenes were annotated as the “global and overview maps” pathway, 15,786 unigenes were annotated as the “translation” pathway and 10,830 unigenes were annotated as the “carbohydrate metabolism” pathway (Fig. 4).
Identification of differentially expressed genes (DEGs), GO and KEGG enrichment analysis
We used edgeR for examining differential expression of replicated count data. For these DEGs, if FDR < 0.05 and log2Fold change ≥1, the DEG was considered as up-regulated but if FDR < 0.05 and log2Fold change ≤ − 1, it was considered as down-regulated. There were 34,964 DEGs identified between Dh_L and Dh_S, including 33,763 up-regulated and 1201 down-regulated DEGs (Additional file 3). A total of 34,125 DEGs were identified between Dh_R and Dh_L, including 2206 up-regulated and 31,919 down-regulated DEGs (Additional file 4). There were 29,037 DEGs between Dh_R and Dh_S, 15,352 of which were up-regulated and 13,685 of which were down-regulated (Additional file 5) (Fig. 5a-c). Using a Venn diagram, we compared the three data sets from different comparison groups (Dh_L vs. Dh_S, Dh_R vs. Dh_L and Dh_R vs. Dh_S). In all three comparison groups, 3915 DEGs were identified in common (Fig. 5d). In detail, 11,641 DEGs were identified in both “Dh_L vs. Dh_S” and “Dh_R vs. Dh_L” comparisons; 11,855 DEGs were identified in both “Dh_L vs. Dh_S” and “Dh_R vs. Dh_S” comparisons; while 12,752 DEGs were identified in both “Dh_R vs. Dh_S” and “Dh_R vs. Dh_L” comparisons.
To get a better understanding of the DEGs, we performed GO and KEGG enrichment analysis on all DEGs. The GO enrichment of DEGs was conducted by the hypergeometric Fisher exact test, in which p-value was calculated and adjusted as corrected p-value. GO terms with corrected p-value < 0.05 were considered to be significantly enriched. For pathway enrichment analysis, all DEGs were assigned to terms in KEGG database and searched for significantly enriched KEGG terms with the same analytic approach. The GO enrichment is shown in Additional file 2: Figure S1 and Additional files 6, 7 and 8. In the KEGG enrichment analysis, the top enriched pathway was flavonoid biosynthesis with 21 DEGs and 37 background unigenes in the Dh_L vs. Dh_S comparison. In the Dh_R vs. Dh_L comparison, the top pathway was diterpenoid biosynthesis with 12 DEGs and 22 background unigenes. Finally, in the Dh_R vs. Dh_S comparison, the top pathway was brassinosteroid biosynthesis with 9 DEGs and 18 background unigenes. The top 20 KEGG enrichment pathways are shown in Fig. 6.
Identifying D. huoshanense carbohydrate-active related genes and DEGs related to polysaccharides
The 1407 carbohydrate-active related unigenes were identified and divided into five gene families (Fig. 7). These unigenes include 447 glycosyltransferase genes (GTs), 818 glycoside hydrolases (GHs), 60 carbohydrate esterases (CEs), 62 carbohydrate-binding modules (CBMs), and 20 polysaccharide lyases (PLs) (Fig. 7). Carbohydrate esterase accounts for the largest share of 58.14%. However, polysaccharide lyases accounts for the smallest, only 1.42%.
The glycosyltransferases (GTs) are an important and functionally diverse family of enzymes, the main role of GTs in plants is in the biosynthesis of polysaccharides in the plant cell wall . A total of 315 DEGs related to polysaccharide have been identified among these three comparisons. These polysaccharide-related DEGs include 197 GTs, 51 glucosyltransferases (GTFs), 12 fucosyltransferases (FucTs), 34 mannosyltransferases (MNTs), 8 xylosyltransferases (XTs) and 13 galactosyltransferases (GATs) DEGs (Fig. 8 and Additional file 2: Figure S2). In the Dh_L vs. Dh_S comparison, 115 DEGs related to polysaccharides were identified, including 64 up-regulated and 3 down-regulated DEGs of GTs; 22 up-regulated and 1 down-regulated DEGs of GTFs; 2 up-regulated DEGs of FucTs; 19 up-regulated DEGs of MNTs; 3 up-regulated DEGs of XTs and 4 up-regulated DEGs of GATs. In the Dh_R vs. Dh_L comparison, 99 DEGs were identified, including 6 up-regulated and 59 down-regulated of GTs; 2 up-regulated and 14 down-regulated DEGs of GTFs; 2 down-regulated DEGs of FucTs; 1 up-regulated and 6 down-regulated DEGs of MNTs; 1 up-regulated and 3 down-regulated DEGs of XTs; 1 up-regulated and 4 down-regulated DEGs of GATs. In the Dh_R vs. Dh_S comparison, 101 DEGs were identified, including 34 up-regulated and 34 down-regulated DEGs of GTs; 5 up-regulated and 7 down-regulated DEGs of GTFs; 3 up-regulated and 5 down-regulated DEGs of FucTs; 7 up-regulated and 1 down-regulated DEGs of MNTs; 1 up-regulated of XTs and 4 down-regulated DEGs of GATs.
DEGs related to the biosynthesis of alkaloids in D. huoshanense
At present, some studies have shown that the alkaloids in the genus Dendrobium are mostly sesquiterpenoid alkaloids or terpenoid indole alkaloids. Both are derived from the terpenoids pathway (shikimate, MVA or MEP pathways) [17, 20]. DEGs associated with alkaloid biosynthesis in D. huoshanense are shown in Additional file 1: Table S1. In total, 53 unigenes associated with six enzymes were targeted to the shikimate pathway, including 3-deoxy-D-arabinoheptulosonate-DHS7-phosphate (DHS), 3-dehydroquinate synthase (DHQS), 3-dehydroquinate acid dehydratase (DHD), shikimate dehydrogenase (SKDH), 5-enolpyruvylshikimate-3-phosphate synthase (EPSP) and fanesyl diphosphase synthase (FPS). Forty-six unigenes associated with 12 enzymes were located to MEP and MVA pathway. What’s more, we identified two key enzyme-encoding genes involved in strictosidine, including tryptophan decarboxylase (TDC) and β-subunit of tryptophan synthase (TSB) (Additional file 9). Based on FPKM of unigenes, the average expression level of the unigenes associated with each enzyme was also determined. Most of encoding enzymes were predominantly expressed in stems, except for TDC and DXS. These two encoding enzymes were highly expressed in roots than stems (Fig. 9).
The downstream synthetic pathway of D. huoshanense alkaloid is uncertain due to lack of corresponding compound support. It is speculated that the downstream synthesis pathway is mainly the modification of strictosidine, of which CYP450 may catalyze a subset of the monooxygenase reactions and hydroxylation reactions [38, 39]. In our study, 229 unigenes were identified as putative P450 superfamily members, which are listed in Additional file 1: Table S2. The majority are CYP71 family members (7.8%), followed by CYP3A family members (6.1%) and CYP4 family members (4.8%). According to clan classification, the CYP71 subfamily is likely to be involved in secondary metabolism . Alkaloids are nitrogenous amino acid derivatives whose synthesis requires an aminotransferase enzyme. There were 47 putative unigenes annotated as associated with four independent transaminases (Additional file 1: Table S3).
Transcription factors involved in D. huoshanense transcriptome dataset
Transcription factors (TFs) play an important role in regulating the activity of polysaccharide biosynthesis and other secondary metabolism pathways. By comparison with the TFs from the iTAK database (http://itak.feilab.net/cgi-bin/itak/index.cgi), we identified a total of 2579 expressed TFs from our transcriptome. Each of these TFs belongs to 67 known TF families (Table 4). The most abundant TF family is the zinc finger C2H2 TF including 489 unigenes. The C2H2 TF family is one of the largest family of transcription factors in plants and regulates many biological processes such as plant morphogenesis, transcriptional activation and stress . In addition, some TFs, such as C3H, bHLH, bZIP, MYB-related and WRKY, have been previously demonstrated to play a role in D. officinale.
The expression analysis of key enzyme genes by qRT-PCR
To validate changes in gene expression patterns, we identified and examined five key enzyme-encoding genes associated with alkaloid biosynthesis including DXS (c382607_g2, c385678_g1, c383947_g1, c382607_g1, c406425_g1, c385029_g2), DXR (c415809_g2, c415809_g1), HDR (c421755_g2, c353149_g1), HMGR (c406740_g1, c418084_g7, c301772_g1, c368929_g1, c282476_g1), FPS (c343028_g1, c402550_g1, c387734_g1, c210849_g1, c387496_g1, c410441_g3, c379694_g1) using qRT-PCR at the transcriptional level. Primers and sequences are shown in Additional file 1: Table S4. The expression patterns of five key enzyme-encoding genes showed that the highest expression levels of HDR and DXR were in the stem, while the root material had the lowest expression level. In general, the expression of FPS was mostly high in the stem, moderate in leaf material, and relatively low in root material. We found that the transcript level of HMGR and DXS in root was higher than that in the leaf. However, the levels of HMGR and DXS expressed in stem were partly higher than those in root. The expression patterns of the 22 DEGs were consistent with the transcriptome data (r = 0.69078, p-value = 7.70629E-8). These results indicate that our transcriptomic analysis was highly reproducible and reliable (Fig. 10f).
D. huoshanense is a traditional Chinese herbal medicine, and its polysaccharides and alkaloids are main medicinal components. At present, a large number of putative polysaccharides and alkaloids biosynthetic genes have been identified in D. officinale [17, 19]. In D. huoshanense, there are only a few studies on gene cloning, components structure and isolation of polysaccharides and alkaloids. This report is to identify putative genes and metabolic pathways involved in polysaccharides and alkaloid biosynthesis using high-throughput transcriptome sequencing in D. huoshanense. In our study, nine libraries of transcriptome data were obtained, which assembled a total of 478,361 unigenes. Using the BLAST algorithm (E-value<1E-5), they were searched against the Nr, KEGG, Swiss-Prot, Pfam and String databases. Many genes involved polysaccharides and alkaloids biosynthesis have been identified and characterized [17, 19]. Our unigene numbers are almost twofold greater than the unigene numbers from a previous study that focused on four different tissues of D. officinale .
Polysaccharides are the main medicinal ingredients of D. huoshanense, and the polysaccharides of D. huoshanense are composed of glucose, galactose, mannose, arabinose and rhamnose . Glycosyltransferase (GT) is a class of enzymes that catalyze the transfer of sugar from an active donor to a specific receptor molecule, which is ubiquitous in organisms and forms a superfamily of genes. Glycosyltransferase is involved in development, signal transduction, defense and other biological processes. There are a large number of GT families in Arabidopsis and rice (Oryza sativa), with more than 450 genes [43,44,45]. In D. huoshanense, a large number of carbohydrate-related DEGs were identified which includes glycosyltransferase, glucosyltransferase, mannosyltransferase, fucosyltransferase, xylosyltransferase and galactosylatransferase. These genes may play important roles in the synthesis of polysaccharides from D. huoshanense. The results showed that most of these genes were expressed more highly in stems than in leaves and roots, which is consistent with previous work . The differences of D. huoshanense and D. officinale are not only in carbohydrate content but also in the number of genes related to carbohydrate synthesis. In a previous study, 1081 genes related to carbohydrate synthesis were identified in D. officinale, of which 430 GTs, 405 GHs, 150 CEs, 77 CBMs and 19 PLs . Compared with D. officinale, the number of genes related to carbohydrate synthesis of D. huoshanense was more than that of D. officinale, especially in GTs, GHs and PLs. Most of all, the number of GHs in D. huoshanense is twice as great as in D. officinale. D. huoshanense glycoside hydrolysis capacity may be stronger than D. officinale. These results indicate that differences in the number and expression levels of carbohydrate-related genes may influence the content of carbohydrates in different Dendrobium plants.
Based on previous studies, sucrose plays a central role in the growth and differentiation of plants . Sucrose hydrolysis can be derived from a lot of monosaccharides, and these monosaccharides participate in the synthesis of polysaccharides D. officinale . In higher plants, sucrose metabolism is mainly catalyzed by two enzymes: sucrose phosphate synthase (SPS) and sucrose synthase (Susy). In general, the synthesis of sucrose is thought to be catalyzed by SPS, while sucrose decomposition is mainly catalyzed by Susy . In the current market, the quality of Dendrobium depends on the content of soluble polysaccharides. One study has shown that soluble polysaccharides and sucrose metabolic enzymes are closely linked, sucrose metabolism is an important part of polysaccharide synthesis . Liang Yan analyzed the genome sequence of D. officinale, the essential process of polysaccharide synthesis was obtained, and the key enzymes in sucrose synthesis were identified: 10 SPS genes and 15 Susy genes. Compared with other species, these genes were in significantly higher expression, suggesting that the expansions of these genes might related to the richness of polysaccharides . We identified 13 SPS genes and 18 Susy genes in D. huoshanense, for the next step to provide the basis for verifying whether these genes are involved in polysaccharides synthesis in D. huoshanense.
Plant alkaloids are other active ingredients of Dendrobium. In the published D. officinale transcriptome data, the alkaloid synthesis pathway of D. officinale begins with shikimate pathway and is related to the mevalonate (MVA) pathway or the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway, but the complete synthesis pathway is not clear . In D. huoshanense transcriptome, searching through the KEGG database, 6 key enzyme genes involved in the shikimate pathway and 12 key enzyme genes involved in the MVA or MEP pathway were identified. In the terpenoid pathway, the terpene compound precursors mainly have the mevalonate (MVA) pathway and the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway; the end products of both pathways are isopentenyl pyrophosphate (IPP). We selected five key enzyme genes (DXS, DXR, HDR, HMGR and FPS) involved in the MVA and MEP pathways for quantitative validation.
The production processes of the shikimate pathway in plants are not only essential components of protein synthesis, but are also precursors for a wide range of secondary metabolites [50, 51]. EPSP is a shikimate pathway key enzyme involved in the formation of enolpyruvylshikimate 3-phosphate . In D. huoshanense, the stem-specific expressed EPSP encoding genes may increase the metabolic rate of shikimate pathway to produce more tryptamines, which are precursors for strictosidine biosynthesis. But the difference is that in D. officinale, EPSP encoding genes are the leaf-specific expressed genes, suggest EPSP in leaf-specific accumulation of tryptamine in D. officinale. The biosynthesis of secologanin is catalyzed by a series of enzymes associated with MVA and MEP pathways . In the MEP pathway, DXS is the first key enzyme , and DXR is the secondary acting and rate limiting enzyme, which catalyzes the formation of straight chain pentose sugars with a branched isovaleric precursor . In D. huoshanense, the genes encoding DXS and DXR were mainly expressed in stems. Studies have shown that the excessive expression of DXS in Spike lavender results in a significant increase in terpenoid content ; the expression of DXR is enhanced in the hairy roots of Salvia miltiorrhiza, which could effectively control and improve the synthesis of terpenoids . We suggest that there is a role of DXS and DXR in stem-specific accumulation of alkaloid in D. huoshanense. But in D. officinale, they have differences in tissue-specific accumulation of alkaloids, and they are mainly expressed in leaves not in stems . In the MVA pathway, HMGR is the key enzyme in the biosynthesis of plant terpenoids, which is capable of catalyzing 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) to form mevalonate . The overexpression of HMGR1 in ginseng increases the steroids and triterpenes in plants, indicating the effect of promoting ginsenoside biosynthesis . However, the levels of HMGR expressed in stem were partly higher than those in root. This is probably because the biosynthesis of D. huoshanense was completed after a multistep biochemical reaction, the catalytic role of HMGR in the upstream, followed by a variety of enzyme synergies that would affect the synthesis of alkaloids. In D. officinale, the levels of HMGR expressed in root were partly higher than those in stem. All in all, regarding the synthetic encoding genes expressed in D. huoshanense alkaloid, most of genes were expressed highly in stem, but in D. officinale, the result was exactly the opposite. They suggested that the D. officinale alkaloid was in leaf-specific accumulation.
Following the generation of strictosidine, the alkaloid biosynthesis pathway remains unclear. It is speculated that the downstream synthesis pathway is mainly the modification of strictosidine, which may participate in the oxidation reaction. At present, CYP450 involved in alkaloid biosynthesis remain uncharacterized. Cytochrome P450 (CYP450s), a superfamily of monooxygenase, plays critical roles in biosynthesis of plant secondary metabolites such as triterpenes, alkaloids, and sterols . However, no cytochrome P450 involved in downstream in the alkaloid biosynthetic pathway has been cloned or identified in Dendrobium plants. In our study, 229 unigenes were identified as putative P450 superfamily members in D. huoshanense. This result is consistent with D. officinale . In detail, the largest subfamily (subfamily 71) consists of 18 P450 genes, the second largest subfamily (subfamily 3A) consist of 14 P450 genes and the third subfamily (subfamily 4) consists of 11 P450 genes. CYP71 subfamily was reported to be involved in bioactive secondary metabolism in plants . The CYP3A subfamily of enzymes is the most important to drug metabolism in humans because these enzymes metabolize the majority of commercially available drugs . It is seen from the side that D. huoshanense plays a role in drug metabolism in humans, but these are not identified in D. officinale. At present, the research on CYP4 subfamily is mainly focused on animals, which mainly metabolize endogenous substances, and some members may also play a role in the metabolism of exogenous substances [61,62,63]. CYP450 is a major metabolic enzyme for drugs and other endogenous and exogenous substances, which need to be further identified in Dendrobium plants.
Through the functional analysis of D. huoshanense, the candidate genes related to the synthesis of a large number of polysaccharides and alkaloids were identified, and the sequence information was provided for the gene cloning, gene structure analysis and functional verification of these genes. The results of the study will fill the gaps in functional genomics of D. huoshanense, and provide the basis for functional identification, breeding of D. huoshanense and optimization of germplasm resources.
Differential expression genes
3-dehydroquinate acid dehydratase
Environmental Information Processing
False discovery rate
Fragments per kilobase million
Famesyl diphosphase synthase
Genetic Information Processing
Kyoto encyclopedia of genes and genomes
NCBI non-redundant protein sequences
Search tool for the retrieval of interacting genes
β-subunit of tryptophan synthase
Fernando SS, Ormerod P. An annotated checklist of the orchids of Sri Lanka. Rheedea. 2008;18(1):1–28.
Chang CC, Ku AF, Tseng YY, Yang WB, Fang JM, Wong CH. 6,8-Di-C-glycosyl flavonoids from Dendrobium huoshanense. J Nat Prod. 2010;73(2):229–32.
Ng TB, Liu J, Wong JH, Ye X, Wing Sze SC, Tong Y, Zhang KY. Review of research on Dendrobium, a prized folk medicine. Appl Microbiol Biotechnol. 2012;93(5):1795–803.
Mudiam MKR, Jin Q, Jiao C, Sun S, Song C, Cai Y, Lin Y, Fan H, Zhu Y. Metabolic analysis of medicinal Dendrobium officinale and Dendrobium huoshanense during different growth years. PLoS One. 2016;11(1):e0146607.
Hsieh YS, Chien C, Liao SK, Liao SF, Hung WT, Yang WB, Lin CC, Cheng TJ, Chang CC, Fang JM. Structure and bioactivity of the polysaccharides in medicinal plant Dendrobium huoshanense. Bioorg Med Chem. 2008;16(11):6054.
Luo JP, Deng YY, Zha XQ. Mechanism of polysaccharides from Dendrobium huoshanense. On Streptozotocin-induced diabetic cataract. Pharm Biol. 2008;46(4):243–9.
Fan Y, He XJ, Zhou S, Luo A, He T, Chun Z. Composition analysis and antioxidant activity of polysaccharide from Dendrobium denneanum. Int J Biol Macromol. 2009;45(2):169–73.
Luo AX, He XJ, Zhou SD, Fan YJ, Luo AS, Chun Z. Purification, composition analysis and antioxidant activity of the polysaccharides from Dendrobium nobile Lindl. Carbohydr Polym. 2010;79(4):1014–9.
Wang JH, Luo JP, Zha XQ, Feng BJ. Comparison of antitumor activities of different polysaccharide fractions from the stems of Dendrobium nobile Lindl. Carbohydr Polym. 2010;79(1):114–8.
Liu XF, Zhu J, Ge SY, Xia LJ, Yang HY, Qian YT, Ren FZ. Orally administered Dendrobium officinale and its polysaccharides enhance immune functions in BALB/c mice. Nat Prod Commun. 2011;6(6):867.
Xia L, Liu X, Guo H, Zhang H, Zhu J, Ren F. Partial characterization and immunomodulatory activity of polysaccharides from the stem of Dendrobium officinale ( Tiepishihu ) in vitro. J Funct Foods. 2012;4(1):294–301.
Zha XQ, Luo JP, Luo SZ, Jiang ST. Structure identification of a new immunostimulating polysaccharide from the stems of Dendrobium huoshanense. Carbohydr Polym. 2007;69(1):86–93.
Qian XP, Zha XQ, Xiao JJ, Zhang HL, Pan LH, Luo JP. Sulfated modification can enhance antiglycation abilities of polysaccharides from Dendrobium huoshanense. Carbohydr Polym. 2014;101(1):982–9.
Pan LH, Feng BJ, Wang JH, Zha XQ, Luo JP. Structural characterization and anti-glycation activity in vitro of a water-soluble polysaccharide from Dendrobium huoshanense. J Food Biochem. 2013;37(3):313–21.
Morita H, Fujiwara M, Yoshida N, Kobayashi JI. New Picrotoxinin-type and Dendrobine-type Sesquiterpenoids from Dendrobium snowflake ‘red star’. Tetrahedron. 2000;56(32):5801–5.
Shu Y, Qihai G, Qin W, Fei L, Yuanfu L, Jingshan S. Alkaloids enriched extract from Dendrobium nobile Lindl. Attenuates tau protein hyperphosphorylation and apoptosis induced by lipopolysaccharide in rat brain. Phytomedicine. 2014;21(5):712.
Guo X, Li Y, Li C, Luo H, Wang L, Qian J, Luo X, Xiang L, Song J, Sun C. Analysis of the Dendrobium officinale transcriptome reveals putative alkaloid biosynthetic genes and genetic markers. Gene. 2013;527(1):131.
Lichtenthaler HK. The 1-deoxy-d-xylulose-5-phosphate pathway of isoprenoid biosynthesis in plants. Annu Rev Plant Physiol Plant Mol Biol. 1999;50(50):47.
Zhang J, He C, Wu K, Teixeira da Silva JA, Zeng S, Zhang X, Yu Z, Xia H, Duan J. Transcriptome analysis of Dendrobium officinale and its application to the identification of genes associated with polysaccharide synthesis. Front Plant Sci. 2016;7:5.
Shen C, Guo H, Chen H, Shi Y, Meng Y, Lu J, Feng S, Wang H. Identification and analysis of genes associated with the synthesis of bioactive constituents in Dendrobium officinale using RNA-Seq. Sci Rep. 2017;7(1):187.
Yan L, Wang X, Liu H, Tian Y, Lian J, Yang R, Hao S, Wang X, Yang S, Li Q, et al. The genome of Dendrobium officinale illuminates the biology of the important traditional Chinese orchid herb. Mol Plant. 2015;8(6):922–34.
Zhang GQ, Xu Q, Bian C, Tsai WC, Yeh CM, Liu KW, Yoshida K, Zhang LS, Chang SB, Chen F, et al. The Dendrobium catenatum Lindl. Genome sequence provides insights into polysaccharide synthase, floral development and adaptive evolution. Sci Rep. 2016;6:19029.
Murashige TS, Skoog FA. A revised medium for rapid growth and bioassays with tobaco tissue cultures. Physiol Plant. 1962;15(3):473–97.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644.
Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35(Database issue):D61–5.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;27(1):29–34.
Boeckmann B, Bairoch A, Apweiler R, Blatter MC, Estreicher A, Gasteiger E, Martin MJ, Michoud K, O’Donovan C, Phan I. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003;31(1):365–70.
Finn RD, Penelope C, Eberhardt RY, Eddy SR, Jaina M, Mitchell AL, Potter SC, Marco P, Matloob Q, Amaia SV. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 2016;44(Database issue):D279–85.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huertacepas J, Simonovic M, Roth A, Santos A, Tsafou KP. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43(Database issue):D447.
Götz S, Garcíagómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008;36(10):3420–35.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12(1):323.
Robinson MD, Mccarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139.
Tang H, Flick P, Sato K, Ramirez F, Klopfenstein D, Mungall C, Yunes J, Pedersen B. GOATOOLS: tools for gene ontology. Zenodo. 2015.
Mao X, Cai T, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21(19):3787.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
Lao J, Oikawa A, Bromley J, McInerney P, Suttangkakul A, Smith-Moritz A, Plahar H, Chiu T, González Fernández-Niño S, Ebert B, et al. The plant glycosyltransferase clone collection for functional genomics. Plant J. 2014;79(3):517–29.
Morant M, Bak S, Møller BL, Werck-Reichhart D. Plant cytochromes P450: tools for pharmacology, plant protection and phytoremediation. Curr Opin Biotechnol. 2003;14(2):151.
Coon MJ. Cytochrome P450: nature's most versatile biological catalyst. Annu Rev Pharmacol Toxicol. 2005;45(45):1.
Nelson DR, Schuler MA, Paquette SM, Werck-Reichhart D, Bak S. Comparative genomics of Rice and Arabidopsis. Analysis of 727 cytochrome P450 genes and pseudogenes from a monocot and a dicot. Plant Physiol. 2004;135(2):756–72.
Laity JH, Lee BM, Wright PE. Zinc finger proteins: new insights into structural and functional diversity. Curr Opin Struct Biol. 2001;11(1):39.
Meng Y, Yu D, Xue J, Lu J, Feng S, Shen C, Wang H. A transcriptome-wide, organ-specific regulatory map of Dendrobium officinale, an important traditional Chinese orchid herb. Sci Rep. 2016;6:18864.
Lee C, Ye ZH. Focus issue on plant cell walls: the Arabidopsis Family GT43 glycosyltransferases form two functionally nonredundant groups essential for the elongation of Glucuronoxylan backbone. Plant Physiol. 2010;153(2):526–41.
Zhang BL, Zhao TM, Yu WG, Kuang BQ, Yao Y, Liu TL, Chen XY, Zhang WH, Wu AM. Functional conservation of the glycosyltransferase gene GT47A in the monocot rice. J Plant Res. 2014;127(3):423–32.
Xu J, Guan J, Chen XJ, Zhao J, Li SP. Comparison of polysaccharides from different Dendrobium using saccharide mapping. J Pharm Biomed Anal. 2011;55(5):977–83.
Jang JC, Sheen J. Sugar sensing in higher plants. Plant Cell. 1994;6(11):1665–79.
Yang J, Meng H, Yang S, Zhang W, Zha YH, Wen GS. Correlation between soluble polysaccharide and sucrose metabolic enzymes in Dendrobium officinale. J West Chin For Sci. 2012;41:62–7.
Huber SC, Huber JL. Role and regulation of sucrose-phosphate synthase in higher plants. Annu Rev Plant Physiol Plant Mol Biol. 1996;47(2):431.
Wang B, Pan L, Luo J, Zha X. Effect of salicylic acid on cell growth and polysaccharide production in suspension cultures of protocorm-like bodies from Dendrobium huoshanense. Chin J Biotechnol. 2009;25(7):1062–8.
Tzin V, Malitsky S, Zvi MMB, Bedair M, Sumner L, Aharoni A, Galili G. Expression of a bacterial feedback-insensitive 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase of the shikimate pathway in Arabidopsis elucidates potential metabolic bottlenecks between primary and secondary metabolism. New Phytol. 2012;194(2):430–9.
Takayuki T, Mutsumi W, Rainer H, Fernie AR. Shikimate and phenylalanine biosynthesis in the green lineage. Front Plant Sci. 2013;4(1):62.
Klee HJ, Muskopf YM, Gasser CS. Cloning of an Arabidopsis thaliana gene encoding 5-enolpyruvylshikimate-3-phosphate synthase: sequence analysis and manipulation to obtain glyphosate-tolerant plants. Mol Gen Genet. 1987;210(3):437–42.
Estévez J, Cantero A, Reindl A, Reichler S, León P. 1-deoxy-D-xylulose-5-phosphate synthase, a limiting enzyme for plastidic isoprenoid biosynthesis in plants. J Biol Chem. 2001;276(25):22901–9.
Ramak P, Kazempour Osaloo S, Ebrahimzadeh H, Sharifi M, Behmanesh M. Inhibition of the mevalonate pathway enhances carvacrol biosynthesis and DXR gene expression in shoot cultures of Satureja khuzistanica Jamzad. J Plant Physiol. 2013;170(13):1187–93.
Muñoz-Bertomeu J, Arrillaga I, Ros R, Segura J. Up-regulation of 1-deoxy-D-xylulose-5-phosphate synthase enhances production of essential oils in transgenic spike lavender. Plant Physiol. 2006;142(3):890–900.
Yan X, Zhang L, Wang J, Liao P, Zhang Y, Zhang R, Kai G. Molecular characterization and expression of 1-deoxy-D-xylulose 5-osphate reductoisomerase (DXR) gene from Salvia miltiorrhiza. Acta Physiol Plant. 2009;31(5):1015–22.
Chappell J, Wolf F, Proulx J, Cuellar R, Saunders C. Is the reaction catalyzed by 3-Hydroxy-3-Methylglutaryl coenzyme a reductase a rate-limiting step for isoprenoid biosynthesis in plants? Plant Physiol. 1995;109(4):1337–43.
Kim YJ, Lee OR, Ji YO, Jang MG, Yang DC. Functional analysis of HMGR encoding genes in triterpene saponin-producing Panax ginseng Meyer. Plant Physiol. 2014;165(1):373-87.
Niu Y, Luo H, Huang L, Chen S, He S. Advances in the study of CYP450 involved in Ginsenoside biosynthesis. World Sci Technol. 2012.
Wrighton SA, Schuetz EG, Thummel KE, Shen DD, Korzekwa KR, Watkins PB. The human CYP3A subfamily: practical consideration. Drug Metab Rev. 2000;32(3–4):339.
Claudianos C, Ranson H, Johnson RM, Biswas S, Schuler MA, Berenbaum MR, Feyereisen R, Oakeshott JG. A deficit of detoxification enzymes: pesticide sensitivity and environmental response in the honeybee. Insect Mol Biol. 2006;15(5):615–36.
Maïbèchecoisne M. cDNA cloning of biotransformation enzymes belonging to the cytochrome P450 family in the antennae of the noctuid moth Mamestra brassicae. Insect Mol Biol. 2002;11(3):273–81.
Jørgensen A, Rasmussen LJ, Andersen O. Characterisation of two novel CYP4 genes from the marine polychaete Nereis virens and their involvement in pyrene hydroxylase activity. Biochem Biophys Res Commun. 2005;336(3):890.
We would thank Prof. Donald L. DeAngelis, from the Wetland and Aquatic Research Center, U.S. Geological Survey, for his valuable comments and suggestions that have greatly improved the quality of this manuscript.
This project was supported by National Special Fund for Forestry Scientific Research in the Public Interest (Grant No. 201504406), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), Jiangsu Agriculture Science and Technology Innovation Fund (Grant No. CX (17)1004), the National traditional Chinese medicine standardization project (ZYYS-2014 ).
Availability of data and materials
All sequence data have been deposited at the National Center for Biotechnology Information short read archive (SRA) under the accession number SRP122499.
Ethics approval and consent to participate
Dendrobium huoshanense used in this study were cultivated by Maoyun Yu’s greenhouse. No permits were required for the collection of the samples.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. DEGs associated with alkaloid biosynthesis in D. huoshanense. Table S2. Putative cytochrome P450s involved in D. huoshanense transcriptome. Table S3. The information of 47 putative unigenes associated with four independent transaminases. Table S4. Genes IDs and primers used in the quantitative real-time PCR (qRT-PCR) experiments. (DOCX 34 kb)
Figure S1. Functional gene ontology classification of unigenes. Figure S2. The heatmap of important DEGs associated with glycosyltransferase. Figure S3 Venn diagram of all unigenes with annotations against five public databases. (DOCX 1469 kb)
Differential expression genes between Dh_L and Dh_S. (XLS 37782 kb)
Differential expression genes between Dh_R and Dh_L. (XLS 37941 kb)
Differential expression genes between Dh_R and Dh_S. (XLS 38011 kb)
Gene Ontology (GO) enrichment of D. huoshanense transcriptome between Dh_L and Dh_S. (XLS 4088 kb)
Gene Ontology (GO) enrichment of D. huoshanense transcriptome between Dh_R and Dh_L. (XLS 4064 kb)
Gene Ontology (GO) enrichment of D. huoshanense transcriptome between Dh_R and Dh_S. (XLS 2487 kb)
Detailed information about unigenes associated with putative alkaloid biosynthesis in D. huoshanense. (XLSX 60 kb)
About this article
Cite this article
Yuan, Y., Yu, M., Jia, Z. et al. Analysis of Dendrobium huoshanense transcriptome unveils putative genes associated with active ingredients synthesis. BMC Genomics 19, 978 (2018). https://doi.org/10.1186/s12864-018-5305-6
- Dendrobium huoshanense
- Polysaccharides synthesis
- Alkaloids synthesis