Overexpression of an isopentenyl diphosphate isomerase gene to enhance trans-polyisoprene production in Eucommia ulmoides Oliver

EuHDZ25 positively affects rubber biosynthesis by targeting EuFPS1 in Eucommia leaves

Eucommia ulmoides is a temperate gum source plant that produces trans-polyisoprene (TPI), also known as Eucommia rubber. The structural configuration and function of TPI offer a new material with important potential for industrial development. In this study, we detected the TPI content in the leaves of diploid and triploid E. ulmoides plants. The average TPI content in the leaves of triploid E. ulmoides was significantly higher than that of diploid. Transcriptome data and weighted gene co-expression network analyses identified a significant positive correlation between the EuFPS1 gene and TPI content. Overexpression of EuFPS1 increased the density of rubber particles and TPI content, indicating its crucial role in TPI biosynthesis. In addition, the expression of EuHDZ25 in E. ulmoides was significantly positively correlated with EuFPS1 expression. Yeast one-hybrid and dual-luciferase assays demonstrated that EuHDZ25 mainly promotes TPI biosynthesis through positive regulation of EuFPS1 expression. The significantly up-regulated expression of EuHDZ25 and its consequent upregulation of EuFPS1 during the biosynthesis of TPI may partially explain the increased TPI content of triploids. This study provides an important theoretical foundation for further exploring the molecular mechanism of secondary metabolites content variation in polyploids and can help to promote the development and utilization of rubber resources.

Positive regulation of the Eucommia rubber biosynthesis-related gene EuFPS1 by EuWRKY30 in Eucommia ulmoides

Eucommia rubber is a secondary metabolite from Eucommia ulmoides that has attracted much attention because of its unique properties and enormous potential for application. However, the transcriptional mechanism regulating its biosynthesis has not yet been determined. Farnesyl pyrophosphate synthase is a key enzyme in the Eucommia rubber biosynthesis. In this study, the promoter of EuFPS1 was used as bait, EuWRKY30 was screened from the cDNA library of EuFPS1 via a yeast one-hybrid system. EuWRKY30 belongs to the WRKY IIa subfamily and contains a WRKY domain and a C2H2 zinc finger motif, and the expressed protein is located in the nucleus. EuWRKY30 and EuFPS1 exhibited similar tissue expression patterns, and yeast one-hybrid and dual-luciferase experiments confirmed that EuWRKY30 directly binds to the W-box element in the EuFPS1 promoter and activates its expression. Moreover, the overexpression of EuWRKY30 significantly upregulated the expression level of EuFPS1, further increasing the density of the rubber particles and Eucommia rubber content. The results of this study indicated that EuWRKY30 positively regulates EuFPS1, which plays a critical role in the synthesis of Eucommia rubber, provided a basis for further analysis of the underlying transcriptional regulatory mechanisms.

Full-length transcriptome profiling of Acanthopanax gracilistylus provides new insight into the kaurenoic acid biosynthesis pathway

Acanthopanax gracilistylus is a deciduous plant in the family Araliaceae, which is commonly used in Chinese herbal medicine, as the root bark has functions of nourishing the liver and kidneys, removing dampness and expelling wind, and strengthening the bones and tendons. Kaurenoic acid (KA) is the main effective substance in the root bark of A. gracilistylus with strong anti-inflammatory effects. To elucidate the KA biosynthesis pathway, second-generation (DNA nanoball) and third-generation (Pacific Biosciences) sequencing were performed to analyze the transcriptomes of the A. gracilistylus leaves, roots, and stems. Among the total 505,880 isoforms, 408,954 were annotated by seven major databases. Sixty isoforms with complete open reading frames encoding 11 key enzymes involved in the KA biosynthesis pathway were identified. Correlation analysis between isoform expression and KA content identified a total of eight key genes. Six key enzyme genes involved in KA biosynthesis were validated by real-time quantitative polymerase chain reaction. Based on the sequence analysis, the spatial structure of ent-kaurene oxidase was modeled, which plays roles in the three continuous oxidations steps of KA biosynthesis. This study greatly enriches the transcriptome data of A. gracilistylus and facilitates further analysis of the function and regulation mechanism of key enzymes in the KA biosynthesis pathway.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01436-7.

Integrated metabolome and transcriptome analysis identifies candidate genes involved in triterpenoid saponin biosynthesis in leaves of Centella asiatica (L.) Urban

Centella asiatica (L.) Urban is a well-known medicinal plant which has multiple pharmacological properties. Notably, the leaves of C. asiatica contain large amounts of triterpenoid saponins. However, there have only been a few studies systematically elucidating the metabolic dynamics and transcriptional differences regarding triterpenoid saponin biosynthesis during the leaf development stages of C. asiatica. Here, we performed a comprehensive analysis of the metabolome and transcriptome to reveal the dynamic patterns of triterpenoid saponin accumulation and identified the key candidate genes associated with their biosynthesis in C. asiatica leaves. In this study, we found that the key precursors in the synthesis of terpenoids, including DMAPP, IPP and β-amyrin, as well as 22 triterpenes and eight triterpenoid saponins were considered as differentially accumulated metabolites. The concentrations of DMAPP, IPP and β-amyrin showed significant increases during the entire stage of leaf development. The levels of 12 triterpenes decreased only during the later stages of leaf development, but five triterpenoid saponins rapidly accumulated at the early stages, and later decreased to a constant level. Furthermore, 48 genes involved in the MVA, MEP and 2, 3-oxidosqualene biosynthetic pathways were selected following gene annotation. Then, 17 CYP450s and 26 UGTs, which are respectively responsible for backbone modifications, were used for phylogenetic-tree construction and time-specific expression analysis. From these data, by integrating metabolomics and transcriptomics analyses, we identified CaHDR1 and CaIDI2 as the candidate genes associated with DMAPP and IPP synthesis, respectively, and CaβAS1 as the one regulating β-amyrin synthesis. Two genes from the CYP716 family were confirmed as CaCYP716A83 and CaCYP716C11. We also selected two UGT73 families as candidate genes, associated with glycosylation of the terpenoid backbone at C-3 in C. asiatica. These findings will pave the way for further research on the molecular mechanisms associated with triterpenoid saponin biosynthesis in C. asiatica.

Comparative analysis of latex transcriptomes reveals the potential mechanisms underlying rubber molecular weight variations between the Hevea brasiliensis clones RRIM600 and Reyan7-33-97

BACKGROUND

The processabilities and mechanical properties of natural rubber depend greatly on its molecular weight (MW) and molecular weight distribution (MWD). However, the mechanisms underlying the regulation of molecular weight during rubber biosynthesis remain unclear.

RESULTS

In the present study, we determined the MW and particle size of latex from 1-year-old virgin trees and 30-year-old regularly tapped trees of the Hevea clones Reyan7-33-97 and RRIM600. The results showed that both the MW and the particle size of latex varied between these two clones and increased with tree age. Latex from RRIM600 trees had a smaller average particle size than that from Reyan7-33-97 trees of the same age. In 1-year-old trees, the Reyan7-33-97 latex displayed a slightly higher MW than that of RRIM600, whereas in 30-year-old trees, the RRIM600 latex had a significantly higher MW than the Reyan7-33-97 latex. Comparative analysis of the transcriptome profiles indicated that the average rubber particle size is negatively correlated with the expression levels of rubber particle associated proteins, and that the high-MW traits of latex are closely correlated with the enhanced expression of isopentenyl pyrophosphate (IPP) monomer-generating pathway genes and downstream allylic diphosphate (APP) initiator-consuming non-rubber pathways. By bioinformatics analysis, we further identified a group of transcription factors that potentially regulate the biosynthesis of IPP.

CONCLUSIONS

Altogether, our results revealed the potential regulatory mechanisms involving gene expression variations in IPP-generating pathways and the non-rubber isoprenoid pathways, which affect the ratios and contents of IPP and APP initiators, resulting in significant rubber MW variations among same-aged trees of the Hevea clones Reyan7-33-97 and RRIM600. Our findings provide a better understanding of rubber biosynthesis and lay the foundation for genetic improvement of rubber quality in H. brasiliensis.

Natural Polymer Eucommia Ulmoides Rubber: A Novel Material

As the second natural rubber resource, Eucommia ulmoides rubber (EUR) from Eucommia ulmoides Oliver is mainly composed of trans-1,4-polyisoprene, which is the isomer of natural rubber cis-1,4-polyisoprene from Hevea brasiliensis. In the past few years, the great potential application of EUR has received increasing attention, and there is a growing awareness that the natural polymer EUR could become an emerging research topic in field of the novel materials due to its unique and excellent duality of both rubber and plastic. To gain insight into its further development, in this review, the extraction, structure, physicochemical properties, and modification of EUR are discussed in detail. More emphasis on the potential applications in the fields of the environment, agriculture, engineering, and biomedical engineering is summarized. Finally, some insights into the challenges and perspectives of EUR are also suggested.

Genome-wide association study reveals the genes associated with the leaf inclusion contents in Chinese medical tree Eucommia ulmoides

Eucommia ulmoides is an economic tree that can biosynthesize secondary metabolites with pharmacological functions. Genetic basis of biosynthesis of these compounds is almost unknown. Therefore, genomic-wide association study was performed to exploit the genetic loci maybe involved in biosynthetic pathways of 5 leaf inclusions (aucubin, chlorogenic acid, gutta-percha, polyphenols, total flavonoids). It was shown that contents of the 5 leaf metabolites have a wide variation following normal distribution. A total of 2 013 102 single nucleotide polymorphism (SNP) markers were identified in a population containing 62 individual clones. Through genome-wide association study analysis, many SNP loci were identified perhaps associated with phenotypes of the leaf inclusions. Higher transcriptional levels of the candidate genes denoted by significant SNPs in leaves suggested they may be involved in biosynthesis of the leaf inclusions. These genetic loci provide with invaluable information for further studies on the gene functions in biosynthesis of the leaf inclusions and selective breeding of the plus trees.

Molecular Sex Identification in the Hardy Rubber Tree (Eucommia ulmoides Oliver) via ddRAD Markers

Eucommia ulmoides, also known as the industrially and medicinally important hardy rubber tree, is the sole species of Eucommiaceae. Nevertheless, its dioecious property hinders sex recognition by traditional morphological observation at very early developmental stages, thus inhibiting breeding and economic cropping. In this study, double-digest restriction site-associated DNA sequencing (ddRAD-seq) was applied to screen sex-linked molecular markers for sex identification and investigation of the sex determination system in 20 male and female E. ulmoides individual plants, respectively. In consequence, five candidate male-specific loci but no female-specific loci were predicated among the 183,752 male and 147,122 female catalogue loci by bioinformatics analysis. Subsequent PCR (polymerase chain reaction) amplification and Sanger sequencing examinations were performed on another 24 individuals, 12 for each sex, from a separate population. One ideal sex-linked locus, MSL4, was identified among the five putative male-specific loci that were found using ddRAD data. MSL4 is 479 bp in length and highly conserved in all the male individuals, suggesting its feature of being stable and repeatable. Our results also indicated that the sex of E. ulmoides is likely determined genetically. In short, this study provides a consistent and reproducible ddRAD marker (MSL4) that is able to discriminate male from female seedlings in E. ulmoides, which will be valuable for rapid breeding practice and better commercial production of this economically important tree.

Integration of Transcriptomes, Small RNAs, and Degradome Sequencing to Identify Putative miRNAs and Their Targets Related to Eu-Rubber Biosynthesis in Eucommia ulmoides

Eucommia ulmoides has attracted much attention as a valuable natural rubber (Eu-rubber) production tree. As a strategic material, Eu-rubber plays a vital role in general and defence industries. However, the study of Eu-rubber biosynthesis at a molecular level is scarce, and the regulatory network between microRNAs (miRNAs) and messenger RNAs (mRNAs) in Eu-rubber biosynthesis has not been assessed. In this study, we comprehensively analyzed the transcriptomes, small RNAs (sRNAs) and degradome to reveal the regulatory network of Eu-rubber biosynthesis in E. ulmoides. A total of 82,065 unigenes and 221 miRNAs were identified using high-throughput sequencing; 20,815 targets were predicted using psRNATarget software. Of these targets, 779 miRNA-target pairs were identified via degradome sequencing. Thirty-one miRNAs were differentially expressed; 22 targets of 34 miRNAs were annotated in the terpenoid backbone biosynthesis pathway (ko00900) based on the Kyoto Encyclopedia of Genes and Genomes (KEGG). These miRNAs were putatively related to Eu-rubber biosynthesis. A regulatory network was constructed according to the expression profiles of miRNAs and their targets. These results provide a comprehensive analysis of transcriptomics, sRNAs and degradome to reveal the Eu-rubber accumulation, and provide new insights into genetic engineering techniques which may improve the content of Eu-rubber in E. ulmoides.

The Differential Expression of Mevalonate Pathway Genes in the Gut of the Bark Beetle Dendroctonus rhizophagus (Curculionidae: Scolytinae) Is Unrelated to the de Novo Synthesis of Terpenoid Pheromones

Bark beetles commonly produce de novo terpenoid pheromones using precursors synthesized through the mevalonate pathway. This process is regulated by Juvenile Hormone III (JH III). In this work, the expression levels of mevalonate pathway genes were quantified after phloem feeding-to induce the endogenous synthesis of JH III-and after the topical application of a JH III solution. The mevalonate pathway genes from D. rhizophagus were cloned, molecularly characterized, and their expression levels were quantified. Also, the terpenoid compounds produced in the gut were identified and quantified by Gas Chromatography Mass Spectrometry (GC-MS). The feeding treatment produced an evident upregulation, mainly in acetoacetyl-CoA thiolase (AACT), 3-hydroxy-3-methylglutaryl-CoA synthase (HMGS), 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), phosphomevalonate kinase (PMK), and isopentenyl diphosphate isomerase (IPPI) genes, and males reached higher expression levels compared to females. In contrast, the JH III treatment did not present a clear pattern of upregulation in any sex or time. Notably, the genes responsible for the synthesis of frontalin and ipsdienol precursors (geranyl diphosphate synthase/farnesyl diphosphate synthase (GPPS/FPPS) and geranylgeranyl diphosphate synthase (GGPPS)) were not clearly upregulated, nor were these compounds further identified. Furthermore, trans-verbenol and myrtenol were the most abundant compounds in the gut, which are derived from an α-pinene transformation rather than de novo synthesis. Hence, the expression of mevalonate pathway genes in D. rhizophagus gut is not directed to the production of terpenoid pheromones, regardless of their frequent occurrence in the genus Dendroctonus.

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.

Identification of the Sex-Biased Gene Expression and Putative Sex-Associated Genes in Eucommia ulmoides Oliver Using Comparative Transcriptome Analyses

Eucommia ulmoides is a model representative of the dioecious plants with sex differentiation at initiation. Nevertheless, the genetic mechanisms of sexual dimorphism and sex determination in E. ulmoides remain poorly understood. In this study de novo transcriptome sequencing on Illumina platform generated >45 billion high-quality bases from fresh leaves of six male and female individuals of E. ulmoides. A total of 148,595 unigenes with an average length of 801 base-pairs (bp) were assembled. Through comparative transcriptome analyses, 116 differentially expressed genes (DEGs) between the males and the females were detected, including 73 male-biased genes and 43 female-biased genes. Of these DEGs, three female-biased genes were annotated to be related with the sexually dimorphic gutta content in E. ulmoides. One male-biased DEG was identified as putative MADS box gene APETALA3, a B class floral organ identity gene in the flowering plants. SNPs calling analyses further confirmed that the APETALA3-like gene was probably involved in the sex determination in E. ulmoides. Four other male-biased DEGs were potential sex-associated genes as well with segregated SNPs in accord with sex type. In addition, the SNPs density was 1.02 per kilobase (kb) in the expressed genes of E. ulmoides, implying a relatively high genetic diversity.

Role of isopentenyl-diphosphate isomerase in heterologous cyanobacterial (Synechocystis) isoprene production

Heterologous production of isoprene (C₅H₈) hydrocarbons in cyanobacteria, emanating from sunlight, CO₂, and water, is now attracting increasing attention. The concept entails application of an isoprene synthase transgene from terrestrial plants, heterologously expressed in cyanobacteria, aiming to reprogram carbon flux in the terpenoid biosynthetic pathway toward formation and spontaneous release of this volatile chemical from the cell and liquid culture. However, flux manipulations and carbon-partitioning reactions between isoprene (the product) and native terpenoid biosynthesis for cellular needs are not yet optimized for isoprene yield. The primary reactant for isoprene biosynthesis is dimethylallyl diphosphate (DMAPP), whereas both DMAPP and its isopentenyl diphosphate (IPP) isomer are needed for cellular terpenoid biosynthesis. The present work addressed the function of an isopentenyl diphosphate (IPP) isomerase in cyanobacteria and its role in carbon partitioning between IPP and DMAPP, both of which serve, in variable ratios, as reactants for the synthesis of different cellular terpenoids. The work was approached upon the heterologous expression in Synechocystis of the "isopentenyl diphosphate isomerase" gene (FNI) from Streptococcus pneumoniae, using isoprene production as a "reporter process" for substrate partitioning between DMAPP and IPP. It is shown that transgenic expression of the FNI gene in Synechocystis resulted in a 250 % increase in the "reporter isoprene" rate and yield, suggesting that the FNI isomerase shifted the endogenous DMAPP-IPP steady-state pool size toward DMAPP, thereby enhancing rates and yield of isoprene production. The work provides insight into the significance and functional role of the IPP isomerase in these photosynthetic microorganisms.

Genome-Wide Identification of MicroRNAs and Their Targets in the Leaves and Fruits of Eucommia ulmoides Using High-Throughput Sequencing

MicroRNAs (miRNAs), a group of endogenous small non-coding RNAs, play important roles in plant growth, development, and stress response processes. Eucommia ulmoides Oliver (hardy rubber tree) is one of the few woody plants capable of producing trans-1, 4-polyisoprene (TPI), also known as Eu-rubber, which has been utilized as an industrial raw material and is extensively cultivated in China. However, the mechanism of TPI biosynthesis has not been identified in E. ulmoides. To characterize small RNAs and their targets with potential biological roles involved in the TPI biosynthesis in E. ulmoides, in the present study, eight small RNA libraries were constructed and sequenced from young and mature leaves and fruits of E. ulmoides. Further analysis identified 34 conserved miRNAs belonging to 20 families (two unclassified families), and 115 novel miRNAs seemed to be specific to E. ulmoides. Among these miRNAs, fourteen conserved miRNAs and 49 novel miRNAs were significantly differentially expressed and identified as Eu-rubber accumulation related miRNAs. Based on the E. ulmoides genomic data, 202 and 306 potential target genes were predicted for 33 conserved and 92 novel miRNAs, respectively; the predicted targets are mostly transcription factors and functional genes, which were enriched in metabolic pathways and biosynthesis of secondary metabolites. Noticeably, based on the expression patterns of miRNAs and their target genes in combination with the Eu-rubber accumulation, the negative correlation of expression of six miRNAs (Eu-miR14, Eu-miR91, miR162a, miR166a, miR172c, and miR396a) and their predicted targets serving as potential regulators in Eu-rubber accumulation. This study is the first to detect conserved and novel miRNAs and their potential targets in E. ulmoides and identify several candidate genes potentially controlling rubber accumulation, and thus provide molecular evidence for understanding the roles of miRNAs in regulating the TPI biosynthesis in E. ulmoides.

Large-Scale Evolutionary Analysis of Genes and Supergene Clusters from Terpenoid Modular Pathways Provides Insights into Metabolic Diversification in Flowering Plants

An important component of plant evolution is the plethora of pathways producing more than 200,000 biochemically diverse specialized metabolites with pharmacological, nutritional and ecological significance. To unravel dynamics underlying metabolic diversification, it is critical to determine lineage-specific gene family expansion in a phylogenomics framework. However, robust functional annotation is often only available for core enzymes catalyzing committed reaction steps within few model systems. In a genome informatics approach, we extracted information from early-draft gene-space assemblies and non-redundant transcriptomes to identify protein families involved in isoprenoid biosynthesis. Isoprenoids comprise terpenoids with various roles in plant-environment interaction, such as pollinator attraction or pathogen defense. Combining lines of evidence provided by synteny, sequence homology and Hidden-Markov-Modelling, we screened 17 genomes including 12 major crops and found evidence for 1,904 proteins associated with terpenoid biosynthesis. Our terpenoid genes set contains evidence for 840 core terpene-synthases and 338 triterpene-specific synthases. We further identified 190 prenyltransferases, 39 isopentenyl-diphosphate isomerases as well as 278 and 219 proteins involved in mevalonate and methylerithrol pathways, respectively. Assessing the impact of gene and genome duplication to lineage-specific terpenoid pathway expansion, we illustrated key events underlying terpenoid metabolic diversification within 250 million years of flowering plant radiation. By quantifying Angiosperm-wide versatility and phylogenetic relationships of pleiotropic gene families in terpenoid modular pathways, our analysis offers significant insight into evolutionary dynamics underlying diversification of plant secondary metabolism. Furthermore, our data provide a blueprint for future efforts to identify and more rapidly clone terpenoid biosynthetic genes from any plant species.