Cloning and Characterization of the 2- C -Methyl- D -erythritol 4-Phosphate (MEP) Pathway Genes of a Natural-Rubber Producing Plant, Hevea brasiliensis

Integrated physiology, transcriptome and proteome analyses highlight the potential roles of multiple hormone-mediated signaling pathways involved in tapping panel dryness in rubber tree

Currently, one of the most serious threats to rubber tree is the tapping panel dryness (TPD) that greatly restricts natural rubber production. Over-tapping or excessive ethephon stimulation is regarded as the main cause of TPD occurrence. Although extensive studies have been carried out, the molecular mechanism underlying TPD remains puzzled. An attempt was made to compare the levels of endogenous hormones and the profiles of transcriptome and proteome between healthy and TPD trees. Results showed that most of endogenous hormones such as jasmonic acid (JA), 1-aminocyclopropanecarboxylic acid (ACC), indole-3-acetic acid (IAA), trans-zeatin (tZ) and salicylic acid (SA) in the barks were significantly altered in TPD-affected rubber trees. Accordingly, multiple hormone-mediated signaling pathways were changed. In total, 731 differentially expressed genes (DEGs) and 671 differentially expressed proteins (DEPs) were identified, of which 80 DEGs were identified as putative transcription factors (TFs). Further analysis revealed that 12 DEGs and five DEPs regulated plant hormone synthesis, and that 16 DEGs and six DEPs were involved in plant hormone signal transduction pathway. Nine DEGs and four DEPs participated in rubber biosynthesis and most DEGs and all the four DEPs were repressed in TPD trees. All these results highlight the potential roles of endogenous hormones, signaling pathways mediated by these hormones and rubber biosynthesis pathway in the defense response of rubber trees to TPD. The present study extends our understanding of the nature and mechanism underlying TPD and provides some candidate genes and proteins related to TPD for further research in the future.

Molecular cloning and heterologous expression analysis of 1-Deoxy-D-Xylulose-5-Phosphate Synthase gene in Centella asiatica L

Many genes involved in triterpenoid saponins in plants control isoprenoid flux and constitute the precursor pool, which is channeled into various downstream pathways leading to the synthesis of triterpenoid saponins in C. asiatica. Full-length 1-Deoxy-D-Xylulose-5-Phosphate-Synthase (CaDXS) gene was isolated for the study from the previously annotated Centella asiatica leaves transcriptomic data. The CaDXS gene sequence was submitted to the NCBI databases with GenBank accession number MZ997832. The full-length CaDXS gene contained a 2244 base pair open reading frame that encoded a 747 amino acid polypeptide. The predicted molecular weight (MW) and theoretical pI of DXS are 76.28 kDa and 6.86, respectively. Multiple amino acid sequence alignment of amino acids and phylogenetic studies suggest that CaDXS shares high similarities with DXS from other plants DXS belonging to different families. A phylogenetic tree was constructed using Molecular Evolutionary Genetic Analysis (MEGA) version 10.1.6. Structural analysis provided fundamental information about the three-dimensional features and physicochemical parameters of the CaDXS protein. Quantitative expression analysis showed that CaDXS transcripts were maximally expressed in leaf, followed by petiole, roots, and node tissues. CaDXS was cloned into the expression vector pET28a, expressed heterologously in DH5α bacteria, confirmed by sequencing, and subsequently characterized by protein expression and functional complementation. The study focused on understanding the protein structure, biological significance, regulatory mechanism, functional analysis, and gene characterization of the centellosides biosynthetic pathway gene DXS for the first time in the plant. It would provide new information about the metabolic pathway and its relative contribution to isoprenoid biosynthesis.

Transcriptome analysis reveals the molecular mechanisms of rubber biosynthesis and laticifer differentiation during rubber seed germination

The molecular mechanisms underlying the initiation of natural rubber synthesis and laticifer differentiation have not been fully elucidated. In this study, we conducted a time-series transcriptome analysis of five rubber tree tissues at four stages of seed germination. A total of 161,199 DEGs were identified between the two groups, including most 16,673 DEGs (A3 vs B3 and A3 vs C3) and lest 1,210 DEGs (C2 vs D2). We found that the maturation of the seed is accompanied by the formation of laticifer cells in cotyledon. Meanwhile, the analysis of hormones related genes expression may provide effective clues for us to promote the differentiation of laticifer cells in seeds by hormones in the future. In this study, hormone-related gene enrichment analyses revealed that IAA, GA, and CTK were activated in laticifer containing tissues. Similarly, GO and GEGG analysis showed that hormone pathways, especially the auxin pathway, are enriched. Gene expression clustering was analyzed using the short time-series expression miner (STEM), and the analysis revealed four distinct trends in the gene expression profiles. Moreover, we enriched transcription factor (TF) enrichment in cotyledon and embryonic axis tissues, and the MYB type exhibited the most significant difference. Furthermore, our findings revealed that genes related to rubber synthesis exhibited tissue-specific expression patterns during seed germination. Notably, key genes associated with rubber biosynthesis, specifically small rubber particle protein (SRPP) and cis-prenyltransferase (CPT), exhibited significant changes in expression in cotyledon and embryonic axis tissues, suggesting synchronous rubber synthesis with seed germination. Our staining results reveled that laticifer cells were exits in the cotyledon before seed imbibition stage. In conclusion, these results lay the foundation for exploring the molecular mechanisms underlying laticifer differentiation and rubber synthesis during seed germination, deepening our understanding of the initiation stages of rubber biosynthesis and laticifer differentiation.

Downregulation of HbFPS1 affects rubber biosynthesis of Hevea brasiliensis suffering from tapping panel dryness

Tapping panel dryness (TPD) is a century-old problem that has plagued the natural rubber production of Hevea brasiliensis. TPD may result from self-protective mechanisms of H. brasiliensis in response to stresses such as excessive hormone stimulation and mechanical wounding (bark tapping). It has been hypothesized that TPD impairs rubber biosynthesis; however, the underlying mechanisms remain poorly understood. In the present study, we firstly verified that TPD-affected rubber trees exhibited lower rubber biosynthesis activity and greater rubber molecular weight compared to healthy rubber trees. We then demonstrated that HbFPS1, a key gene of rubber biosynthesis, and its expression products were downregulated in the latex of TPD-affected rubber trees, as revealed by transcriptome sequencing and iTRAQ-based proteome analysis. We further discovered that the farnesyl diphosphate synthase HbFPS1 could be recruited to small rubber particles by HbSRPP1 through protein-protein interactions to catalyze farnesyl diphosphate (FPP) synthesis and facilitate rubber biosynthesis initiation. FPP content in the latex of TPD-affected rubber trees was significantly decreased with the downregulation of HbFPS1, ultimately resulting in abnormal development of rubber particles, decreased rubber biosynthesis activity, and increased rubber molecular weight. Upstream regulator assays indicated that a novel regulator, MYB2-like, may be an important regulator of downregulation of HbFPS1 in the latex of TPD-affected rubber trees. Our findings not only provide new directions for studying the molecular events involved in rubber biosynthesis and TPD syndrome and contribute to rubber management strategies, but also broaden our knowledge of plant isoprenoid metabolism and its regulatory networks.

Functional characterization of a farnesyl diphosphate synthase from Dendrobium nobile Lindl

Dendrobium nobile Lindl. has been used as a traditional Chinese medicine for a long time, in which the most important compound is dendrobine functioning in a variety of pharmacological activities. Farnesyl diphosphate synthase (FPPS) is one of the key enzymes in the biosynthetic pathway of dendrobine. In this work, we found the expression profiles of DnFPPS were correlated with the contents of dendrobine under the methyl jasmonate (MeJA) treatments at different time. Then, the cloning and functional identification of a novel FPPS from D. nobile. The full length of DnFPPS is 1231 bp with an open reading frame of 1047 bp encoding 348 amino acids. The sequence similarity analysis demonstrated that DnFPPS was in the high homology with Dendrobium huoshanense and Dendrobium catenatum and contained four conserved domains. Phylogenetic analysis showed that DnFPPS was the close to the DhFPPS. Then, DnFPPS was induced to express in Escherichia coli, purified, and identified by SDS-PAGE electrophoresis. Gas chromatography-mass spectrometry analysis indicated that DnFPPS could catalyze dimethylallyl pyrophosphate and isopentenyl pyrophosphate to produce farnesyl diphosphate. Taken together, a novel DnFPPS was cloned and functionally identified, which supplied a candidate gene for the biosynthetic pathway of dendrobine.

Full-Length Transcriptomic Sequencing and Temporal Transcriptome Expression Profiling Analyses Offer Insights into Terpenoid Biosynthesis in Artemisia argyi

Artemisiae argyi Folium is a traditional herbal medicine used for moxibustion heat therapy in China. The volatile oils in A.argyi leaves are closely related to its medicinal value. Records suggest that the levels of these terpenoids components within the leaves vary as a function of harvest time, with June being the optimal time for A. argyi harvesting, owing to the high levels of active ingredients during this month. However, the molecular mechanisms governing terpenoid biosynthesis and the time-dependent changes in this activity remain unclear. In this study, GC–MS analysis revealed that volatile oil levels varied across four different harvest months (April, May, June, and July) in A. argyi leaves, and the primarily terpenoids components (including both monoterpenes and sesquiterpenes) reached peak levels in early June. Through single-molecule real-time (SMRT) sequencing, corrected by Illumina RNA-sequencing (RNA-Seq), 44 full-length transcripts potentially involved in terpenoid biosynthesis were identified in this study. Differentially expressed genes (DEGs) exhibiting time-dependent expression patterns were divided into 12 coexpression clusters. Integrated chemical and transcriptomic analyses revealed distinct time-specific transcriptomic patterns associated with terpenoid biosynthesis. Subsequent hierarchical clustering and correlation analyses ultimately identified six transcripts that were closely linked to the production of these two types of terpenoid within A. argyi leaves, revealing that the structural diversity of terpenoid is related to the generation of the diverse terpene skeletons by prenyltransferase (TPS) family of enzymes. These findings can guide further studies of the molecular mechanisms underlying the quality of A. argyi leaves, aiding in the selection of optimal timing for harvests of A. argyi.

Integrative transcriptomic and metabolomic analyses unveil tanshinone biosynthesis in Salvia miltiorrhiza root under N starvation stress

Salvia miltiorrhiza is a model plant for Chinese herbal medicine with significant pharmacologic effects due to its tanshinone components. Our previous study indicated that nitrogen starvation stress increased its tanshinone content. However, the molecular mechanism of this low nitrogen-induced tanshinone biosynthesis is still unclear. Thus, this study aimed to elucidate the molecular mechanism of tanshinone biosynthesis in S. miltiorrhiza under different N conditions [N-free (N0), low-N (Nl), and full-N (Nf, as control) conditions] by using transcriptome and metabolome analyses. Our results showed 3,437 and 2,274 differentially expressed unigenes between N0 and Nf as well as Nl and Nf root samples, respectively. N starvation (N0 and Nl) promoted the expression of the genes involved in the MVA and MEP pathway of tanshinone and terpenoid backbone biosynthesis. Gene ontology and KEGG analyses revealed that terpenoid backbone biosynthesis, hormone signal transduction, and phenylpropanoid biosynthesis were promoted under N starvation conditions, whereas starch and sucrose metabolisms, nitrogen and phosphorus metabolisms, as well as membrane development were inhibited. Furthermore, metabolome analysis showed that metabolite compounds and biosynthesis of secondary metabolites were upregulated. This study provided a novel insight into the molecular mechanisms of tanshinone production in S. miltiorrhiza in response to nitrogen stress.

RNASeq analysis of drought-stressed guayule reveals the role of gene transcription for modulating rubber, resin, and carbohydrate synthesis

The drought-adapted shrub guayule (Parthenium argentatum) produces rubber, a natural product of major commercial importance, and two co-products with potential industrial use: terpene resin and the carbohydrate fructan. The rubber content of guayule plants subjected to water stress is higher compared to that of well-irrigated plants, a fact consistently reported in guayule field evaluations. To better understand how drought influences rubber biosynthesis at the molecular level, a comprehensive transcriptome database was built from drought-stressed guayule stem tissues using de novo RNA-seq and genome-guided assembly, followed by annotation and expression analysis. Despite having higher rubber content, most rubber biosynthesis related genes were down-regulated in drought-stressed guayule, compared to well-irrigated plants, suggesting post-transcriptional effects may regulate drought-induced rubber accumulation. On the other hand, terpene resin biosynthesis genes were unevenly affected by water stress, implying unique environmental influences over transcriptional control of different terpene compounds or classes. Finally, drought induced expression of fructan catabolism genes in guayule and significantly suppressed these fructan biosynthesis genes. It appears then, that in guayule cultivation, irrigation levels might be calibrated in such a regime to enable tunable accumulation of rubber, resin and fructan.

Transcriptome analysis of Pará rubber tree (H. brasiliensis) seedlings under ethylene stimulation

BACKGROUND

Natural rubber (cis-1,4-polyioprene, NR) is an indispensable industrial raw material obtained from the Pará rubber tree (H. brasiliensis). Natural rubber cannot be replaced by synthetic rubber compounds because of the superior resilience, elasticity, abrasion resistance, efficient heat dispersion, and impact resistance of NR. In NR production, latex is harvested by periodical tapping of the trunk bark. Ethylene enhances and prolongs latex flow and latex regeneration. Ethephon, which is an ethylene-releasing compound, applied to the trunk before tapping usually results in a 1.5- to 2-fold increase in latex yield. However, intense mechanical damage to bark tissues by excessive tapping and/or over-stimulation with ethephon induces severe oxidative stress in laticifer cells, which often causes tapping panel dryness (TPD) syndrome. To enhance NR production without causing TPD, an improved understanding of the molecular mechanism of the ethylene response in the Pará rubber tree is required. Therefore, we investigated gene expression in response to ethephon treatment using Pará rubber tree seedlings as a model system.

RESULTS

After ethephon treatment, 3270 genes showed significant differences in expression compared with the mock treatment. Genes associated with carotenoids, flavonoids, and abscisic acid biosynthesis were significantly upregulated by ethephon treatment, which might contribute to an increase in latex flow. Genes associated with secondary cell wall formation were downregulated, which might be because of the reduced sugar supply. Given that sucrose is an important molecule for NR production, a trade-off may arise between NR production and cell wall formation for plant growth and for wound healing at the tapping panel.

CONCLUSIONS

Dynamic changes in gene expression occur specifically in response to ethephon treatment. Certain genes identified may potentially contribute to latex production or TPD suppression. These data provide valuable information to understand the mechanism of ethylene stimulation, and will contribute to improved management practices and/or molecular breeding to attain higher yields of latex from Pará rubber trees.

Transcriptional Analysis of Metabolic Pathways and Regulatory Mechanisms of Essential Oil Biosynthesis in the Leaves of Cinnamomum camphora (L.) Presl

The roots, bark, and leaves of Cinnamomum camphora are rich in essential oils, which mainly comprised monoterpenes and sesquiterpenes. Although the essential oils obtained from C. camphora have been widely used in pharmaceutical, medicinal, perfume, and food industries, the molecular mechanisms underlying terpenoid biosynthesis are poorly understood. To address this lack of knowledge, we performed transcriptome analysis to investigate the key regulatory genes involved in terpenoid biosynthesis in C. camphora. High-oil-yield trees of linalool type and low-oil-yield trees were used to assemble a de novo transcriptome of C. camphora. A total of 121,285 unigenes were assembled, and the total length, average length, N50, and GC content of unigenes were 87,869,987, 724, 1,063, and 41.1%, respectively. Comparison of the transcriptome profiles of linalool-type C. camphora with trees of low oil yield resulted in a total of 3,689 differentially expressed unigenes, among them 31 candidate genes had annotations associated with metabolism of terpenoids and polyketides, including four in the monoterpenoid biosynthesis pathway and three in the terpenoid backbone biosynthesis pathway. Collectively, this genome-wide transcriptome provides a valuable tool for future identification of genes related to essential oil biosynthesis. Additionally, the identification of a cohort of genes in the biosynthetic pathways of terpenoids provides a theoretical basis for metabolic engineering of essential oils in C. camphora.

Regulatory Potential of bHLH-Type Transcription Factors on the Road to Rubber Biosynthesis in Hevea brasiliensis

Natural rubber is the main component of latex obtained from laticifer cells of Hevea brasiliensis. For improving rubber yield, it is essential to understand the genetic molecular mechanisms responsible for laticifer differentiation and rubber biosynthesis. Jasmonate enhances both secondary laticifer differentiation and rubber biosynthesis. Here, we carried out time-course RNA-seq analysis in suspension-cultured cells treated with methyljasmonic acid (MeJA) to characterize the gene expression profile. Gene Ontology (GO) analysis showed that the term "cell differentiation" was enriched in upregulated genes at 24 hours after treatment, but inversely, the term was enriched in downregulated genes at 5 days, indicating that MeJA could induce cell differentiation at an early stage of the response. Jasmonate signaling is activated by MYC2, a basic helix-loop-helix (bHLH)-type transcription factor (TF). The aim of this work was to find any links between transcriptomic changes after MeJA application and regulation by TFs. Using an in vitro binding assay, we traced candidate genes throughout the whole genome that were targeted by four bHLH TFs: Hb_MYC2-1, Hb_MYC2-2, Hb_bHLH1, and Hb_bHLH2. The latter two are highly expressed in laticifer cells. Their physical binding sites were found in the promoter regions of a variety of other TF genes, which are differentially expressed upon MeJA exposure, and rubber biogenesis-related genes including SRPP1 and REF3. These studies suggest the possibilities that Hb_MYC2-1 and Hb_MYC2-2 regulate cell differentiation and that Hb_bHLH1 and Hb_bHLH2 promote rubber biosynthesis. We expect that our findings will help to increase natural rubber yield through genetic control in the future.

Molecular Mechanisms of Natural Rubber Biosynthesis

Natural rubber (NR), principally comprising cis-1,4-polyisoprene, is an industrially important natural hydrocarbon polymer because of its unique physical properties, which render it suitable for manufacturing items such as tires. Presently, industrial NR production depends solely on latex obtained from the Pará rubber tree, Hevea brasiliensis. In latex, NR is enclosed in rubber particles, which are specialized organelles comprising a hydrophobic NR core surrounded by a lipid monolayer and membrane-bound proteins. The similarity of the basic carbon skeleton structure between NR and dolichols and polyprenols, which are found in most organisms, suggests that the NR biosynthetic pathway is related to the polyisoprenoid biosynthetic pathway and that rubber transferase, which is the key enzyme in NR biosynthesis, belongs to the cis-prenyltransferase family. Here, we review recent progress in the elucidation of molecular mechanisms underlying NR biosynthesis through the identification of the enzymes that are responsible for the formation of the NR backbone structure.

Molecular cloning, characterization, and heterologous expression of an acetyl-CoA acetyl transferase gene from Sanghuangporus baumii

Acetyl-CoA C-acetyltransferase synthase gene (AACT) cDNA, DNA and promoter were cloned from Sanghuangporus baumii. The gene ORF (1260 bp) encoded 419 amino acids. The AACT DNA includes five exons (1-84 bp, 140-513 bp, 570-1027 bp, 1090-1282 bp, 1344-1494 bp) and four introns (85-139 bp, 514-569 bp, 1028-1089 bp, 1283-1343 bp). The molecular weight of AACT protein is 43.40 kDa, it is hydrophilic with a theoretical isoelectric point of 8.96. Furthermore, The region of the transcription start site is 1997-2047 bp of AACT promoter, and it contained promoter elements (TATA Boxs, CAAT Boxs, CAAT-box, ABRE, G-Boxs, Sp1, MSA-like, LTR). AACT recombinant protein (43.40 KDa + Tag protein 22.68 KDa) was subjected in SDS-PAGE. AACT the transcription levels of in different development stages were investigated. The expression of AACT in primordia (2.4-fold) and 15 d mycelia (2.3- fold) were significantly higher than 9 d mycelia (contral). The expression level of the AACT downstream genes and triterpenoids content were determined at different developmental stages. Triterpenoid content reached its peak on day 15(7.21 mg/g).

Single-nucleotide polymorphism markers within MVA and MEP pathways among Hevea brasiliensis clones through transcriptomic analysis

In this study, RNA sequencing of several Hevea brasiliensis clones grown in Malaysia with different annual rubber production yields and disease resistance was performed on the Illumina platform. A total of 29,862,548 reads were generated, resulting in 101,269 assembled transcripts that were used as the reference transcripts. A similarity search against the non-redundant (nr) protein databases presented 83,771 (83%) positive BLASTx hits. The transcriptome was annotated using gene ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the Pfam database. A search for putative molecular markers was performed to identify single-nucleotide polymorphisms (SNPs). Overall, 3,210,629 SNPs were detected and a total of 1314 SNPs associated with the genes involved in MVA and MEP pathways were identified. A total of 176 SNP primer pairs were designed from sequences that were related to the MVA and MEP pathways. The transcriptome of RRIM 3001 and RRIM 712 were subjected to pairwise comparison and the results revealed that there were 1262 significantly differentially expressed genes unique to RRIM 3001, 1499 significantly differentially expressed genes unique to RRIM 712 and several genes related to the MVA and MEP pathways such as AACT, HMGS, PMK, MVD, DXS and HDS were included. The results will facilitate the characterization of H. brasiliensis transcriptomes and the development of a new set of molecular markers in the form of SNPs from transcriptome assembly for the genotype identification of various rubber varieties with superior traits in Malaysia.

Natural rubber biosynthesis in plants, the rubber transferase complex, and metabolic engineering progress and prospects

Natural rubber (NR) is a nonfungible and valuable biopolymer, used to manufacture ~50 000 rubber products, including tires and medical gloves. Current production of NR is derived entirely from the para rubber tree (Hevea brasiliensis). The increasing demand for NR, coupled with limitations and vulnerability of H. brasiliensis production systems, has induced increasing interest among scientists and companies in potential alternative NR crops. Genetic/metabolic pathway engineering approaches, to generate NR-enriched genotypes of alternative NR plants, are of great importance. However, although our knowledge of rubber biochemistry has significantly advanced, our current understanding of NR biosynthesis, the biosynthetic machinery and the molecular mechanisms involved remains incomplete. Two spatially separated metabolic pathways provide precursors for NR biosynthesis in plants and their genes and enzymes/complexes are quite well understood. In contrast, understanding of the proteins and genes involved in the final step(s)-the synthesis of the high molecular weight rubber polymer itself-is only now beginning to emerge. In this review, we provide a critical evaluation of recent research developments in NR biosynthesis, in vitro reconstitution, and the genetic and metabolic pathway engineering advances intended to improve NR content in plants, including H. brasiliensis, two other prospective alternative rubber crops, namely the rubber dandelion and guayule, and model species, such as lettuce. We describe a new model of the rubber transferase complex, which integrates these developments. In addition, we highlight the current challenges in NR biosynthesis research and future perspectives on metabolic pathway engineering of NR to speed alternative rubber crop commercial development.

Proteomic Landscape of the Mature Roots in a Rubber-Producing Grass Taraxacum Kok-saghyz

The rubber grass Taraxacum kok-saghyz (TKS) contains large amounts of natural rubber (cis-1,4-polyisoprene) in its enlarged roots and it is an alternative crop source of natural rubber. Natural rubber biosynthesis (NRB) and storage in the mature roots of TKS is a cascade process involving many genes, proteins and their cofactors. The TKS genome has just been annotated and many NRB-related genes have been determined. However, there is limited knowledge about the protein regulation mechanism for NRB in TKS roots. We identified 371 protein species from the mature roots of TKS by combining two-dimensional gel electrophoresis (2-DE) and mass spectrometry (MS). Meanwhile, a large-scale shotgun analysis of proteins in TKS roots at the enlargement stage was performed, and 3545 individual proteins were determined. Subsequently, all identified proteins from 2-DE gel and shotgun MS in TKS roots were subject to gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses and most proteins were involved in carbon metabolic process with catalytic activity in membrane-bounded organelles, followed by proteins with binding ability, transportation and phenylpropanoid biosynthesis activities. Fifty-eight NRB-related proteins, including eight small rubber particle protein (SRPP) and two rubber elongation factor(REF) members, were identified from the TKS roots, and these proteins were involved in both mevalonate acid (MVA) and methylerythritol phosphate (MEP) pathways. To our best knowledge, it is the first high-resolution draft proteome map of the mature TKS roots. Our proteomics of TKS roots revealed both MVA and MEP pathways are important for NRB, and SRPP might be more important than REF for NRB in TKS roots. These findings would not only deepen our understanding of the TKS root proteome, but also provide new evidence on the roles of these NRB-related proteins in the mature TKS roots.

Transcriptomic profiling reveals MEP pathway contributing to ginsenoside biosynthesis in Panax ginseng

Background

Panax ginseng C. A. Mey is one of famous medicinal herb plant species. Its major bioactive compounds are various ginsenosides in roots and rhizomes. It is commonly accepted that ginsenosides are synthesized from terpene precursors, IPP and DMAPP, through the cytoplasmic mevalonate (MVA) pathway. Another plastic 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway was proved also contributing to ginsenoside generation in the roots of P. ginseng by using specific chemical inhibitors recently. But their gene expression characteristics are still under reveal in P. ginseng. With the development of the high-throughput next generation sequencing (NGS) technologies, we have opportunities to discover more about the complex ginsenoside biosynthesis pathways in P. ginseng.

Results

We carried out deep RNA sequencing and comprehensive analyses on the ginseng root samples of 1–5 years old and five different tissues of 5 years old ginseng plants. The de novo assembly totally generated 48,165 unigenes, including 380 genes related to ginsenoside biosynthesis and all the genes encoding the enzymes of the MEP pathway and the MVA pathway. We further illustrated the gene expression profiles related to ginsenoside biosynthesis among 1–5 year-old roots and different tissues of 5 year-old ginseng plants. Particularly for the first time, we revealed that the gene transcript abundances of the MEP pathway were similar to those of the MVA pathway in ginseng roots but higher in ginseng leaves. The IspD was predicated to be the rate-limiting enzyme in the MEP pathway through both co-expression network and gene expression profile analyses.

Conclusions

At the transcriptional level, the MEP pathway has similar contribution to ginsenoside biosynthesis in ginseng roots, but much higher in ginseng leaves, compared with the MVA pathway. The IspD might be the key enzyme for ginsenoside generation through the MEP pathway. These results provide new information for further synthetic biology study on ginsenoside metabolic regulation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5718-x) contains supplementary material, which is available to authorized users.

Biosynthesis of Natural Rubber: Current State and Perspectives

Natural rubber is a kind of indispensable biopolymers with great use and strategic importance in human society. However, its production relies almost exclusively on rubber-producing plants Hevea brasiliensis, which have high requirements for growth conditions, and the mechanism of natural rubber biosynthesis remains largely unknown. In the past two decades, details of the rubber chain polymerization and proteins involved in natural rubber biosynthesis have been investigated intensively. Meanwhile, omics and other advanced biotechnologies bring new insight into rubber production and development of new rubber-producing plants. This review summarizes the achievements of the past two decades in understanding the biosynthesis of natural rubber, especially the massive information obtained from the omics analyses. Possibilities of natural rubber biosynthesis in vitro or in genetically engineered microorganisms are also discussed.

Over-expression of 3-hydroxy-3- methylglutaryl-coenzyme A reductase 1 (hmgr1) gene under super-promoter for enhanced latex biosynthesis in rubber tree (Hevea brasiliensis Muell. Arg.)

Natural rubber (cis-1, 4-polyisoprene) is being produced from bark laticifer cells of Hevea brasiliensis and the popular high latex yielding Indian rubber clones are easily prone to onset of tapping panel dryness syndrome (TPD) which is considered as a physiological syndrome affecting latex production either partially or completely. This report describes an efficient protocol for development of transgenic rubber plants by over-expression of 3-hydroxy 3-methylglutaryl Co-enzyme A reductase 1 (hmgr1) gene which is considered as rate limiting factor for latex biosynthesis via Agrobacterium-mediated transformation. The pBIB plasmid vector containing hmgr1 gene cloned under the control of a super-promoter was used for genetic transformation using embryogenic callus. Putatively transgenic cell lines were obtained on selection medium and produced plantlets with 44% regeneration efficiency. Transgene integration was confirmed by PCR amplification of 1.8 kb hmgr1 and 0.6 kb hpt genes from all putatively transformed callus lines as well as transgenic plants. Southern blot analysis showed the stable integration and presence of transgene in the transgenic plants. Over expression of hmgr1 transgene was determined by Northern blot hybridization, semi-quantitative PCR and real-time PCR (qRT-PCR) analysis. Accumulation of hmgr1 mRNA transcripts was more abundant in transgenic plants than control. Increased level of photosynthetic pigments, protein contents and HMGR enzyme activity was also noticed in transgenic plants over control. Interestingly, the latex yield was significantly enhanced in all transgenic plants compared to the control. The qRT-PCR results exhibit that the hmgr1 mRNA transcript levels was 160-fold more abundance in transgenic plants over untransformed control. These results altogether suggest that there is a positive correlation between latex yield and accumulation of mRNA transcripts level as well as HMGR enzyme activity in transgenic rubber plants. It is presumed that there is a possibility for enhanced level of latex biosynthesis in transgenic plants as the level of mRNA transcripts and HMGR enzyme activity is directly correlated with latex yield in rubber tree. Further, the present results clearly suggest that the quantification of HMGR enzyme activity in young seedlings will be highly beneficial for early selection of high latex yielding plants in rubber breeding programs.

Molecular cloning and expression analysis of two key genes, HDS and HDR, in the MEP pathway in Pyropia haitanensis

The 1-hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate synthase (HDS) gene and the 1-hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate reductase (HDR) gene are two important genes in the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway. In this study, we reported the isolation and characterization of full-length HDS (MF101802) and HDR (MF159558) from Pyropia haitanensis. Characteristics of 3-D structures of the PhHDS and PhHDR proteins were analysed respectively. The results showed that the full-length cDNA of PhHDS, which is 1801 bp long, contained a 1455 bp open reading frame (ORF) encoding a putative 484 amino acid residue protein with a predicted molecular mass of 51.60 kDa. Meanwhile, the full-length cDNA of PhHDR was 1668 bp and contained a 1434 bp ORF encoding a putative 477 amino acid 2 residue protein with a predicted molecular mass of 51.49 kDa. The expression levels of the two genes were higher in conchocelis than that in leafy thallus. Additionally, the expression levels could be influenced by light, temperature and salinity and induced by methyl jasmonate (MJ) and salicylic acid (SA). This study contributed to our in-depth understanding of the roles of PhHDS and PhHDR in terpenoid biosynthesis in Pyropia haitanensis and the regulation of the two genes by external environments.

Isoprenoid biosynthesis in dandelion latex is enhanced by the overexpression of three key enzymes involved in the mevalonate pathway

BACKGROUND

Latex from the dandelion species Taraxacum brevicorniculatum contains many high-value isoprenoid end products, e.g. triterpenes and polyisoprenes such as natural rubber. The isopentenyl pyrophosphate units required as precursors for these isoprenoids are provided by the mevalonate (MVA) pathway. The key enzyme in this pathway is 3-hydroxy-methyl-glutaryl-CoA reductase (HMGR) and its activity has been thoroughly characterized in many plant species including dandelion. However, two enzymes acting upstream of HMGR have not been characterized in dandelion latex: ATP citrate lyase (ACL), which provides the acetyl-CoA utilized in the MVA pathway, and acetoacetyl-CoA thiolase (AACT), which catalyzes the first step in the pathway to produce acetoacetyl-CoA. Here we isolated ACL and AACT genes from T. brevicorniculatum latex and characterized their expression profiles. We also overexpressed the well-characterized HMGR, ACL and AACT genes from Arabidopsis thaliana in T. brevicorniculatum to determine their impact on isoprenoid end products in the latex.

RESULTS

The spatial and temporal expression profiles of T. brevicorniculatum ACL and AACT revealed their pivotal role in the synthesis of precursors necessary for isoprenoid biosynthesis in latex. The overexpression of A. thaliana ACL and AACT and HMGR in T. brevicorniculatum latex resulted in the accumulation of all three enzymes, increased the corresponding enzymatic activities and ultimately increased sterol levels by ~5-fold and pentacyclic triterpene and cis-1,4-isoprene levels by ~2-fold. Remarkably high levels of the triterpene precursor squalene were also detected in the triple-transgenic lines (up to 32 mg/g root dry weight) leading to the formation of numerous lipid droplets which were observed in root cross-sections.

CONCLUSIONS

We could show the effective expression of up to three transgenes in T. brevicorniculatum latex which led to increased enzymatic activity and resulted in high level squalene accumulation in the dandelion roots up to an industrially relevant amount. Our data provide insight into the regulation of the MVA pathway in dandelion latex and can be used as a basis for metabolic engineering to enhance the production of isoprenoid end products in this specialized tissue.

The rubber tree genome reveals new insights into rubber production and species adaptation

The Para rubber tree (Hevea brasiliensis) is an economically important tropical tree species that produces natural rubber, an essential industrial raw material. Here we present a high-quality genome assembly of this species (1.37 Gb, scaffold N50 = 1.28 Mb) that covers 93.8% of the genome (1.47 Gb) and harbours 43,792 predicted protein-coding genes. A striking expansion of the REF/SRPP (rubber elongation factor/small rubber particle protein) gene family and its divergence into several laticifer-specific isoforms seem crucial for rubber biosynthesis. The REF/SRPP family has isoforms with sizes similar to or larger than SRPP1 (204 amino acids) in 17 other plants examined, but no isoforms with similar sizes to REF1 (138 amino acids), the predominant molecular variant. A pivotal point in Hevea evolution was the emergence of REF1, which is located on the surface of large rubber particles that account for 93% of rubber in the latex (despite constituting only 6% of total rubber particles, large and small). The stringent control of ethylene synthesis under active ethylene signalling and response in laticifers resolves a longstanding mystery of ethylene stimulation in rubber production. Our study, which includes the re-sequencing of five other Hevea cultivars and extensive RNA-seq data, provides a valuable resource for functional genomics and tools for breeding elite Hevea cultivars.

Transcriptome sequencing and analysis of rubber tree (Hevea brasiliensis Muell.) to discover putative genes associated with tapping panel dryness (TPD)

Background

Tapping panel dryness (TPD) involves in the partial or complete cessation of latex flow thus seriously affect latex production in rubber tree (Hevea brasiliensis). Numerous studies have been conducted to define its origin and nature, but the molecular nature and mechanism of TPD occurrence remains unknown. This study is committed to de novo sequencing and comparative analysis of the transcriptomes of healthy (H) and TPD-affected (T) rubber trees to identify the genes and pathways related to the TPD.

Results

Total raw reads of 34,632,012 and 35,913,020 bp were obtained from H and T library, respectively using Illumina Hiseq 2000 sequencing technology. De novo assemblies yielded 141,456 and 169,285 contigs, and 96,070 and 112,243 unigenes from H and T library, respectively.

Among 73597 genes, 22577 genes were identified as differential expressed genes between H and T library via comparative transcript profiling. A majority of genes involved in natural rubber biosynthesis and jasmonate synthesis with most potential relevance in TPD occurrence were found to be differentially expressed.

Conclusions

In TPD-affected trees, the expression of most genes related to the latex biosynthesis and jasmonate synthesis was severely inhibited and is probably the direct cause of the TPD. These new de novo transcriptome data sets provide a significant resource for the discovery of genes related to TPD and improve our understanding of the occurrence and maintainace of TPD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1562-9) contains supplementary material, which is available to authorized users.

UPLC/Q-TOF MS-based metabolomics and qRT-PCR in enzyme gene screening with key role in triterpenoid saponin biosynthesis of Polygala tenuifolia

BACKGROUND

The dried root of Polygala tenuifolia, named Radix Polygalae, is a well-known traditional Chinese medicine. Triterpenoid saponins are some of the most important components of Radix Polygalae extracts and are widely studied because of their valuable pharmacological properties. However, the relationship between gene expression and triterpenoid saponin biosynthesis in P. tenuifolia is unclear.

METHODOLOGY/FINDINGS

In this study, ultra-performance liquid chromatography (UPLC) coupled with quadrupole time-of-flight mass spectrometry (Q-TOF MS)-based metabolomic analysis was performed to identify and quantify the different chemical constituents of the roots, stems, leaves, and seeds of P. tenuifolia. A total of 22 marker compounds (VIP>1) were explored, and significant differences in all 7 triterpenoid saponins among the different tissues were found. We also observed an efficient reference gene GAPDH for different tissues in this plant and determined the expression level of some genes in the triterpenoid saponin biosynthetic pathway. Results showed that MVA pathway has more important functions in the triterpenoid saponin biosynthesis of P. tenuifolia. The expression levels of squalene synthase (SQS), squalene monooxygenase (SQE), and beta-amyrin synthase (β-AS) were highly correlated with the peak area intensity of triterpenoid saponins compared with data from UPLC/Q-TOF MS-based metabolomic analysis.

CONCLUSIONS/SIGNIFICANCE

This finding suggested that a combination of UPLC/Q-TOF MS-based metabolomics and gene expression analysis can effectively elucidate the mechanism of triterpenoid saponin biosynthesis and can provide useful information on gene discovery. These findings can serve as a reference for using the overexpression of genes encoding for SQS, SQE, and/or β-AS to increase the triterpenoid saponin production of P. tenuifolia.

De novo assembly and transcriptome analysis of the rubber tree (Hevea brasiliensis) and SNP markers development for rubber biosynthesis pathways

Hevea brasiliensis (Willd. Ex Adr. Juss.) Muell.-Arg. is the primary source of natural rubber that is native to the Amazon rainforest. The singular properties of natural rubber make it superior to and competitive with synthetic rubber for use in several applications. Here, we performed RNA sequencing (RNA-seq) of H. brasiliensis bark on the Illumina GAIIx platform, which generated 179,326,804 raw reads on the Illumina GAIIx platform. A total of 50,384 contigs that were over 400 bp in size were obtained and subjected to further analyses. A similarity search against the non-redundant (nr) protein database returned 32,018 (63%) positive BLASTx hits. The transcriptome analysis was annotated using the clusters of orthologous groups (COG), gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Pfam databases. A search for putative molecular marker was performed to identify simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs). In total, 17,927 SSRs and 404,114 SNPs were detected. Finally, we selected sequences that were identified as belonging to the mevalonate (MVA) and 2-C-methyl-D-erythritol 4-phosphate (MEP) pathways, which are involved in rubber biosynthesis, to validate the SNP markers. A total of 78 SNPs were validated in 36 genotypes of H. brasiliensis. This new dataset represents a powerful information source for rubber tree bark genes and will be an important tool for the development of microsatellites and SNP markers for use in future genetic analyses such as genetic linkage mapping, quantitative trait loci identification, investigations of linkage disequilibrium and marker-assisted selection.

Leaf-, panel- and latex-expressed sequenced tags from the rubber tree (Hevea brasiliensis) under cold-stressed and suboptimal growing conditions: the development of gene-targeted functional markers for stress response

Hevea brasiliensis is a native species of the Amazon Basin of South America and the primary source of natural rubber worldwide. Due to the occurrence of South American Leaf Blight disease in this area, rubber plantations have been extended to suboptimal regions. Rubber tree breeding is time-consuming and expensive, but molecular markers can serve as a tool for early evaluation, thus reducing time and costs. In this work, we constructed six different cDNA libraries with the aim of developing gene-targeted molecular markers for the rubber tree. A total of 8,263 reads were assembled, generating 5,025 unigenes that were analyzed; 912 expressed sequence tags (ESTs) represented new transcripts, and two sequences were highly up-regulated by cold stress. These unigenes were scanned for microsatellite (SSR) regions and single nucleotide polymorphisms (SNPs). In total, 169 novel EST-SSR markers were developed; 138 loci were polymorphic in the rubber tree, and 98 % presented transferability to six other Hevea species. Locus duplication was observed in H. brasiliensis and other species. Additionally, 43 SNP markers in 13 sequences that showed similarity to proteins involved in stress response, latex biosynthesis and developmental processes were characterized. cDNA libraries are a rich source of SSR and SNP markers and enable the identification of new transcripts. The new markers developed here will be a valuable resource for linkage mapping, QTL identification and other studies in the rubber tree and can also be used to evaluate the genetic variability of other Hevea species, which are valuable assets in rubber tree breeding.

Draft genome sequence of the rubber tree Hevea brasiliensis

Background

Hevea brasiliensis, a member of the Euphorbiaceae family, is the major commercial source of natural rubber (NR). NR is a latex polymer with high elasticity, flexibility, and resilience that has played a critical role in the world economy since 1876.

Results

Here, we report the draft genome sequence of H. brasiliensis. The assembly spans ~1.1 Gb of the estimated 2.15 Gb haploid genome. Overall, ~78% of the genome was identified as repetitive DNA. Gene prediction shows 68,955 gene models, of which 12.7% are unique to Hevea. Most of the key genes associated with rubber biosynthesis, rubberwood formation, disease resistance, and allergenicity have been identified.

Conclusions

The knowledge gained from this genome sequence will aid in the future development of high-yielding clones to keep up with the ever increasing need for natural rubber.

Distinct expression patterns of two Ginkgo biloba 1-hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate reductase/isopentenyl diphospahte synthase (HDR/IDS) promoters in Arabidopsis model

1-Hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate reductase (HDR) or isopentenyl diphosphate synthase (IDS) is an enzyme at the final step of the MEP pathway. The multi-copy nature of IDS gene in a gymnosperm Ginkgo biloba is known. To evaluate the function of each isogene, the roles of the promoters were examined in Arabidopsis model. Among the promoters of GbIDS series, about 1.3 kb of GbIDS1pro and 1.5 kb of GbIDS2pro were cloned and fused with GUS. The GbIDS1pro::GUS was introduced into Arabidopsis to show GUS expression in most organs except for roots, petals, and stamina, whereas the GbIDS2pro::GUS was expressed only in the young leaves, internodes where the flower and shoot branched, and notably in primary root junction. This pattern of GUS expression correlated with high transcript level of GbIDS2 compared to that of GbIDS1 in Ginkgo roots. Methyl jasmonate (MeJA) treatment resulted in down-regulated GbIDS1pro activity in Arabidopsis leaves and upregulated GbIDS2pro activity in roots. The same pattern of gene regulation in roots was also seen upon treatments of gibberellins, abscisic acid, and indole butyric acid.

Construction and analysis of EST libraries of the trans-polyisoprene producing plant, Eucommia ulmoides Oliver

Eucommia ulmoides Oliver is one of a few woody plants capable of producing abundant quantities of trans-polyisoprene rubber in their leaves, barks, and seed coats. One cDNA library each was constructed from its outer stem tissue and inner stem tissue. They comprised a total of 27,752 expressed sequence tags (ESTs) representing 10,520 unigenes made up of 4,302 contigs and 6,218 singletons. Homologues of genes coding for rubber particle membrane proteins that participate in the synthesis of high-molecular poly-isoprene in latex were isolated, as well as those encoding known major latex proteins (MLPs). MLPs extensively shared ESTs, indicating their abundant expression during trans-polyisoprene rubber biosynthesis. The six mevalonate pathway genes which are implicated in the synthesis of isopentenyl diphosphate (IPP), a starting material of poly-isoprene biosynthesis, were isolated, and their role in IPP biosynthesis was confirmed by functional complementation of suitable yeast mutants. Genes encoding five full-length trans-isoprenyl diphosphate synthases were also isolated, and two among those synthesized farnesyl diphosphate from IPP and dimethylallyl diphosphate, an assumed intermediate of rubber biosynthesis. This study should provide a valuable resource for further studies of rubber synthesis in E. ulmoides.

Map-based cloning of zb7 encoding an IPP and DMAPP synthase in the MEP pathway of maize

IspH is a key enzyme in the last step of the methyl-D-erythritol-4-phosphate (MEP) pathway. Loss of function of IspH can often result in complete yellow or albino phenotype in many plants. Here, we report the characterization of a recessive mutant of maize, zebra7 (zb7), showing transverse green/yellow striped leaves in young plants. The yellow bands of the mutant have decreased levels of chlorophylls and carotenoids with delayed chloroplast development. Low temperature suppressed mutant phenotype, while alternate light/dark cycle or high temperature enlarged the yellow section. Map-based cloning demonstrated that zb7 encodes the IspH protein with a mis-sense mutation in a conserved region. Transgenic silencing of Zb7 in maize resulted in complete albino plantlets that are aborted in a few weeks, confirming that Zb7 is important in the early stages of maize chloroplast development. Zb7 is constitutively expressed and its expression subject to a 16-h light/8-h dark cycle regulation. Our results suggest that the less effective or unstable IspH in zb7 mutant, together with its diurnal expression, are mechanistically accounted for the zebra phenotype. The increased IspH mRNA in the leaves of zb7 at the late development stage may explain the restoration of mutant phenotype in mature stages.

Transcriptome and gene expression analysis in cold-acclimated guayule (Parthenium argentatum) rubber-producing tissue

Natural rubber biosynthesis in guayule (Parthenium argentatum Gray) is associated with moderately cold night temperatures. To begin to dissect the molecular events triggered by cold temperatures that govern rubber synthesis induction in guayule, the transcriptome of bark tissue, where rubber is produced, was investigated. A total of 11,748 quality expressed sequence tags (ESTs) were obtained. The vast majority of ESTs encoded proteins that are similar to stress-related proteins, whereas those encoding rubber biosynthesis-related proteins comprised just over one percent of the ESTs. Sequence information derived from the ESTs was used to design primers for quantitative analysis of the expression of genes that encode selected enzymes and proteins with potential impact on rubber biosynthesis in field-grown guayule plants, including 3-hydroxy-3-methylglutaryl-CoA synthase, 3-hydroxy-3-methylglutaryl-CoA reductase, farnesyl pyrophosphate synthase, squalene synthase, small rubber particle protein, allene oxide synthase, and cis-prenyl transferase. Gene expression was studied for field-grown plants during the normal course of seasonal variation in temperature (monthly average maximum 41.7 °C to minimum 0 °C, from November 2005 through March 2007) and rubber transferase enzymatic activity was also evaluated. Levels of gene expression did not correlate with air temperatures nor with rubber transferase activity. Interestingly, a sudden increase in night temperature 10 days before harvest took place in advance of the highest CPT gene expression level.

De novo assembly and characterization of bark transcriptome using Illumina sequencing and development of EST-SSR markers in rubber tree (Hevea brasiliensis Muell. Arg.)

Background

In rubber tree, bark is one of important agricultural and biological organs. However, the molecular mechanism involved in the bark formation and development in rubber tree remains largely unknown, which is at least partially due to lack of bark transcriptomic and genomic information. Therefore, it is necessary to carried out high-throughput transcriptome sequencing of rubber tree bark to generate enormous transcript sequences for the functional characterization and molecular marker development.

Results

In this study, more than 30 million sequencing reads were generated using Illumina paired-end sequencing technology. In total, 22,756 unigenes with an average length of 485 bp were obtained with de novo assembly. The similarity search indicated that 16,520 and 12,558 unigenes showed significant similarities to known proteins from NCBI non-redundant and Swissprot protein databases, respectively. Among these annotated unigenes, 6,867 and 5,559 unigenes were separately assigned to Gene Ontology (GO) and Clusters of Orthologous Group (COG). When 22,756 unigenes searched against the Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) database, 12,097 unigenes were assigned to 5 main categories including 123 KEGG pathways. Among the main KEGG categories, metabolism was the biggest category (9,043, 74.75%), suggesting the active metabolic processes in rubber tree bark. In addition, a total of 39,257 EST-SSRs were identified from 22,756 unigenes, and the characterizations of EST-SSRs were further analyzed in rubber tree. 110 potential marker sites were randomly selected to validate the assembly quality and develop EST-SSR markers. Among 13 Hevea germplasms, PCR success rate and polymorphism rate of 110 markers were separately 96.36% and 55.45% in this study.

Conclusion

By assembling and analyzing de novo transcriptome sequencing data, we reported the comprehensive functional characterization of rubber tree bark. This research generated a substantial fraction of rubber tree transcriptome sequences, which were very useful resources for gene annotation and discovery, molecular markers development, genome assembly and annotation, and microarrays development in rubber tree. The EST-SSR markers identified and developed in this study will facilitate marker-assisted selection breeding in rubber tree. Moreover, this study also supported that transcriptome analysis based on Illumina paired-end sequencing is a powerful tool for transcriptome characterization and molecular marker development in non-model species, especially those with large and complex genomes.

Genome-wide identification and characterization of novel genes involved in terpenoid biosynthesis in Salvia miltiorrhiza

Terpenoids are the largest class of plant secondary metabolites and have attracted widespread interest. Salvia miltiorrhiza, belonging to the largest and most widely distributed genus in the mint family, is a model medicinal plant with great economic and medicinal value. Diterpenoid tanshinones are the major lipophilic bioactive components in S. miltiorrhiza. Systematic analysis of genes involved in terpenoid biosynthesis has not been reported to date. Searching the recently available working draft of the S. miltiorrhiza genome, 40 terpenoid biosynthesis-related genes were identified, of which 27 are novel. These genes are members of 19 families, which encode all of the enzymes involved in the biosynthesis of the universal isoprene precursor isopentenyl diphosphate and its isomer dimethylallyl diphosphate, and two enzymes associated with the biosynthesis of labdane-related diterpenoids. Through a systematic analysis, it was found that 20 of the 40 genes could be involved in tanshinone biosynthesis. Using a comprehensive approach, the intron/exon structures and expression patterns of all identified genes and their responses to methyl jasmonate treatment were analysed. The conserved domains and phylogenetic relationships among the deduced S. miltiorrhiza proteins and their homologues isolated from other plant species were revealed. It was discovered that some of the key enzymes, such as 1-deoxy-D-xylulose 5-phosphate synthase, 4-hydroxy-3-methylbut-2-enyl diphosphate reductase, hydroxymethylglutaryl-CoA reductase, and geranylgeranyl diphosphate synthase, are encoded by multiple gene members with different expression patterns and subcellular localizations, and both homomeric and heteromeric geranyl diphosphate synthases exist in S. miltiorrhiza. The results suggest the complexity of terpenoid biosynthesis and the existence of metabolic channels for diverse terpenoids in S. miltiorrhiza and provide useful information for improving tanshinone production through genetic engineering.

Laticifer-specific cis-prenyltransferase silencing affects the rubber, triterpene, and inulin content of Taraxacum brevicorniculatum

Certain Taraxacum species, such as Taraxacum koksaghyz and Taraxacum brevicorniculatum, produce large amounts of high-quality natural rubber in their latex, the milky cytoplasm of specialized cells known as laticifers. This high-molecular mass biopolymer consists mainly of poly(cis-1,4-isoprene) and is deposited in rubber particles by particle-bound enzymes that carry out the stereospecific condensation of isopentenyl diphosphate units. The polymer configuration suggests that the chain-elongating enzyme (rubber transferase; EC 2.5.1.20) is a cis-prenyltransferase (CPT). Here, we present a comprehensive analysis of transgenic T. brevicorniculatum plants in which the expression of three recently isolated CPTs known to be associated with rubber particles (TbCPT1 to -3) was heavily depleted by laticifer-specific RNA interference (RNAi). Analysis of the CPT-RNAi plants by nuclear magnetic resonance, size-exclusion chromatography, and gas chromatography-mass spectrometry indicated a significant reduction in rubber biosynthesis and a corresponding 50% increase in the levels of triterpenes and the main storage carbohydrate, inulin. Transmission electron microscopy revealed that the laticifers in CPT-RNAi plants contained fewer and smaller rubber particles than wild-type laticifers. We also observed lower activity of hydroxymethylglutaryl-coenzyme A reductase, the key enzyme in the mevalonate pathway, reflecting homeostatic control of the isopentenyl diphosphate pool. To our knowledge, this is the first in planta demonstration of latex-specific CPT activity in rubber biosynthesis.

Metabolic routes affecting rubber biosynthesis in Hevea brasiliensis latex

The cytosolic mevalonate (MVA) pathway in Hevea brasiliensis latex is the conventionally accepted pathway which provides isopentenyl diphosphate (IPP) for cis-polyisoprene (rubber) biosynthesis. However, the plastidic 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway may be an alternative source of IPP since its more recent discovery in plants. Quantitative RT-PCR (qRT-PCR) expression profiles of genes from both pathways in latex showed that subcellular compartmentalization of IPP for cis-polyisoprene synthesis is related to the degree of plastidic carotenoid synthesis. From this, the occurrence of two schemes of IPP partitioning and utilization within one species is proposed whereby the supply of IPP for cis-polyisoprene from the MEP pathway is related to carotenoid production in latex. Subsequently, a set of latex unique gene transcripts was sequenced and assembled and they were then mapped to IPP-requiring pathways. Up to eight such pathways, including cis-polyisoprene biosynthesis, were identified. Our findings on pre- and post-IPP metabolic routes form an important aspect of a pathway knowledge-driven approach to enhancing cis-polyisoprene biosynthesis in transgenic rubber trees.

Comprehensive transcriptome analysis reveals novel genes involved in cardiac glycoside biosynthesis and mlncRNAs associated with secondary metabolism and stress response in Digitalis purpurea

Abstract

Conclusions

Through comprehensive transcriptome analysis, we not only identified 29 novel gene families potentially involved in the biosynthesis of cardiac glycosides but also characterized a large number of mlncRNAs. Our results suggest the importance of mlncRNAs in secondary metabolism and stress response in D. purpurea.

Two distinct pathways for essential metabolic precursors for isoprenoid biosynthesis

Isoprenoids are a diverse group of molecules found in all organisms, where they perform such important biological functions as hormone signaling (e.g., steroids) in mammals, antioxidation (e.g., carotenoids) in plants, electron transport (e.g., ubiquinone), and cell wall biosynthesis intermediates in bacteria. All isoprenoids are synthesized by the consecutive condensation of the five-carbon monomer isopentenyl diphosphate (IPP) to its isomer, dimethylallyl diphosphate (DMAPP). The biosynthetic pathway for the formation of IPP from acetyl-CoA (i.e., the mevalonate pathway) had been established mainly in mice and the budding yeast Saccharomyces cerevisiae. Curiously, most prokaryotic microorganisms lack homologs of the genes in the mevalonate pathway, even though IPP and DMAPP are essential for isoprenoid biosynthesis in bacteria. This observation provided an impetus to search for an alternative pathway to synthesize IPP and DMAPP, ultimately leading to the discovery of the mevalonate-independent 2-C-methyl-D-erythritol 4-phosphate pathway. This review article focuses on our significant contributions to a comprehensive understanding of the biosynthesis of IPP and DMAPP.

A raison d'être for two distinct pathways in the early steps of plant isoprenoid biosynthesis?

When compared to other organisms, plants are atypical with respect to isoprenoid biosynthesis: they utilize two distinct and separately compartmentalized pathways to build up isoprene units. The co-existence of these pathways in the cytosol and in plastids might permit the synthesis of many vital compounds, being essential for a sessile organism. While substrate exchange across membranes has been shown for a variety of plant species, lack of complementation of strong phenotypes, resulting from inactivation of either the cytosolic pathway (growth and development defects) or the plastidial pathway (pigment bleaching), seems to be surprising at first sight. Hundreds of isoprenoids have been analyzed to determine their biosynthetic origins. It can be concluded that in angiosperms, under standard growth conditions, C₂₀-phytyl moieties, C₃₀-triterpenes and C₄₀-carotenoids are made nearly exclusively within compartmentalized pathways, while mixed origins are widespread for other types of isoprenoid-derived molecules. It seems likely that this coexistence is essential for the interaction of plants with their environment. A major purpose of this review is to summarize such observations, especially within an ecological and functional context and with some emphasis on regulation. This latter aspect still requires more work and present conclusions are preliminary, although some general features seem to exist.

RNA-Seq analysis and de novo transcriptome assembly of Hevea brasiliensis

Hevea brasiliensis, being the only source of commercial natural rubber, is an extremely economically important crop. In an effort to facilitate biological, biochemical and molecular research in rubber biosynthesis, here we report the use of next-generation massively parallel sequencing technologies and de novo transcriptome assembly to gain a comprehensive overview of the H. brasiliensis transcriptome. The sequencing output generated more than 12 million reads with an average length of 90 nt. In total 48,768 unigenes (mean size = 436 bp, median size = 328 bp) were assembled through de novo transcriptome assembly. Out of 13,807 H. brasiliensis cDNA sequences deposited in Genbank of the National Center for Biotechnology Information (NCBI) (as of Feb 2011), 11,746 sequences (84.5%) could be matched with the assembled unigenes through nucleotide BLAST. The assembled sequences were annotated with gene descriptions, Gene Ontology (GO) and Clusters of Orthologous Group (COG) terms. In all, 37,432 unigenes were successfully annotated, of which 24,545 (65.5%) aligned to Ricinus communis proteins. Furthermore, the annotated uingenes were functionally classified according to the GO, COG and Kyoto Encyclopedia of Genes and Genomes databases. Our data provides the most comprehensive sequence resource available for the study of rubber trees as well as demonstrates effective use of Illumina sequencing and de novo transcriptome assembly in a species lacking genomic information.

Screening of valid reference genes for real-time RT-PCR data normalization in Hevea brasiliensis and expression validation of a sucrose transporter gene HbSUT3

Real-time RT-PCR (RT-qPCR) is a sensitive and precise method of quantifying gene expression, however, suitable reference genes are required. Here, a systematic reference gene screening was performed by RT-qPCR on 22 candidate genes in Hevea brasiliensis. Two ubiquitin-protein ligases (UBC2a and UBC4) were the most stable when all samples were analyzed together. A mitosis protein (YLS8) and a eukaryotic translation initiation factor (eIF1Aa) were the most stable in response to tapping. UBC2b and UBC1 were the most stable among different genotypes. UBC2b and a DEAD box RNA helicase (RH2b) were the most stable across individual trees. YLS8 and RH8 were most stably expressed in hormone-treated samples. Expression of the candidate reference genes varied significantly across different tissues, and at least four genes (RH2b, RH8, UBC2a and eIF2) were needed for expression normalization. In addition, examination of relative expression of a sucrose transporter HbSUT3 in different RNA samples demonstrated the importance of additional reference genes to ensure accurate quantitative expression analysis. Overall, our work serves as a guide for selection of reference genes in RT-qPCR gene expression studies in H. brasiliensis.

Ethylene stimulation of latex yield depends on the expression of a sucrose transporter (HbSUT1B) in rubber tree (Hevea brasiliensis)

Hevea brasiliensis is an important industrial crop for natural rubber production. Latex biosynthesis occurs in the cytoplasm of highly specialized latex cells and requires sucrose as the unique precursor. Ethylene stimulation of latex production results in high sugar flow from the surrounding cells of inner bark towards the latex cells. The aim of this work was to understand the role of seven sucrose transporters (HbSUTs) and one hexose transporter (HbHXT1) in this process. Two Hevea clones were used: PB217 and PB260, respectively described as high and low yielding clones. The expression pattern of these sugar transporters (HbSUTs and HbHXT1) was monitored under different physiological conditions and found to be maximal in latex cells. HbSUT1, one of the most abundant isoforms, displayed the greatest response to ethylene treatment. In clone PB217, ethylene treatment led to a higher accumulation of HbSUT1B in latex cells than in the inner bark tissues. Conversely, stronger expression of HbSUT1B was observed in inner bark tissues than in latex cells of PB260. A positive correlation with HbSUT1B transcript accumulation and increased latex production was further supported by its lower expression in latex cells of the virgin clone PB217.

The isogene 1-deoxy-D-xylulose 5-phosphate synthase 2 controls isoprenoid profiles, precursor pathway allocation, and density of tomato trichomes

Plant isoprenoids are formed from precursors synthesized by the mevalonate (MVA) pathway in the cytosol or by the methyl-D-erythritol 4-phosphate (MEP) pathway in plastids. Although some exchange of precursors occurs, cytosolic sesquiterpenes are assumed to derive mainly from MVA, while plastidial monoterpenes are produced preferentially from MEP precursors. Additional complexity arises in the first step of the MEP pathway, which is typically catalyzed by two divergent 1-deoxy-D-xylulose 5-phosphate synthase isoforms (DXS1, DXS2). In tomato (Solanum lycopersicum), the SlDXS1 gene is ubiquitously expressed with highest levels during fruit ripening, whereas SlDXS2 transcripts are abundant in only few tissues, including young leaves, petals, and isolated trichomes. Specific down-regulation of SlDXS2 expression was performed by RNA interference in transgenic plants to investigate feedback mechanisms. SlDXS2 down-regulation led to a decrease in the monoterpene β-phellandrene and an increase in two sesquiterpenes in trichomes. Moreover, incorporation of MVA-derived precursors into residual monoterpenes and into sesquiterpenes was elevated as determined by comparison of ¹³C to ¹²C natural isotope ratios. A compensatory up-regulation of SlDXS1 was not observed. Down-regulated lines also exhibited increased trichome density and showed less damage by leaf-feeding Spodoptera littoralis caterpillars. The results reveal novel, non-redundant roles of DXS2 in modulating isoprenoid metabolism and a pronounced plasticity in isoprenoid precursor allocation.

Characterization of 1-hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate synthase (HDS) gene from Ginkgo biloba

Diterpene trilactone ginkgolides, one of the major constituents of Ginkgo biloba extract, have shown interesting bioactivities including platelet-activating factor antagonistic activity. 1-Hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate synthase (HDS), converting 2-C-methyl-d-erythritol-2,4-cyclodiphosphate into 1-hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate, is the penultimate enzyme of the seven-step 2-C-methyl-d-erythritol 4-phosphate pathway that supplies building blocks for plant isoprenoids of plastid origin such as ginkgolides and carotenoids. Here, we report on the isolation and characterization of the full-length cDNA encoding HDS (GbHDS, GenBank accession number: DQ251630) from G. biloba. Full-length cDNA of GbHDS, 2,763 bp long, contained an ORF of 2,226 bp encoding a protein composed of 741 amino acids. The theoretical molecular weight and pI of the deduced mature GbHDS of 679 amino acid residues are 75.6 kDa and 5.5, respectively. From 2 weeks after initiation of the culture onward, transcription level of this gene in the ginkgo embryo roots increased to about two times higher than that in the leaves. GbHDS was predicted to possess chloroplast transit peptide of 62 amino acid residues, suggesting its putative localization in the plastids. The transient gene expression in Arabidopsis protoplasts confirmed that the transit peptide was capable of delivering the GbHDS protein from the cytosol into the chloroplasts. The isolation and characterization of GbHDS gene enabled us to further understand the role of GbHDS in the terpenoid biosynthesis in G. biloba.