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.)

Comparative proteome and transcriptome analyses suggest the regulation of starch and sucrose metabolism and rubber biosynthesis pathways in the recovery of tapping panel dryness in rubber tree

BACKGROUND

Tapping panel dryness (TPD) in rubber tree has become the most severe restricting factor of natural rubber production. To date, there is no effective measures to prevent and control TPD. Previous studies primarily focused on analyzing the molecular mechanism underlying TPD occurrence. However, there is no research on the molecular mechanism of TPD recovery.

RESULTS

In this study, the TPD trees were recovered by treatment with TPD rehabilitation nutrient agents that could promote the recovery of latex flow on the tapping panel of TPD trees. The genes and proteins involved in TPD recovery were first identified by employing integrated transcriptomics and proteomics analyses. In total, 2029 differentially expressed genes (DEGs) and 951 differentially expressed proteins (DEPs) were detected in the bark of recovery trees compared to that of TPD trees. Among them, 19 DEPs and 11 DEGs were found to be involved in the starch and sucrose metabolism pathway, suggesting their important roles in regulating the syntheses of sucrose and D-glucose, which were the key precursors of natural rubber biosynthesis. Furthermore, 16 DEPs and 15 DEGs were identified in the rubber biosynthesis pathway. Interestingly, almost all the DEPs and DEGs related to rubber biosynthesis exhibited significantly up-regulated expressions in the recovery trees, indicating that latex biosynthesis were probably markedly enhanced during TPD recovery.

CONCLUSIONS

These results provide new insights into the molecular mechanisms underlying TPD recovery, as well as excellent supplements to the mechanisms of TPD occurrence, which will contribute to the development of more effective agents for the prevention and treatment of TPD in the future.

Genome-Wide Identification and Characterization of the HMGR Gene Family in Taraxacum kok-saghyz Provide Insights into Its Regulation in Response to Ethylene and Methyl Jsamonate Treatments

HMGR (3-hydroxy-3-methylglutaryl-CoA reductase) plays a crucial role as the first rate-limiting enzyme in the mevalonate (MVA) pathway, which is the upstream pathway of natural rubber biosynthesis. In this study, we carried out whole-genome identification of Taraxacum kok-saghyz (TKS), a novel rubber-producing alternative plant, and obtained six members of the TkHMGR genes. Bioinformatic analyses were performed including gene structure, protein properties, chromosomal localization, evolutionary relationships, and cis-acting element analyses. The results showed that HMGR genes were highly conserved during evolution with a complete HMG-CoA reductase conserved domain and were closely related to Asteraceae plants during the evolutionary process. The α-helix is the most prominent feature of the secondary structure of the TkHMGR proteins. Collinearity analyses demonstrated that a whole-genome duplication (WGD) event and tandem duplication event play a key role in the expansion of this family and TkHMGR1 and TkHMGR6 have more homologous gene between other species. Cis-acting element analysis revealed that the TkHMGR gene family had a higher number of MYB-related, light-responsive, hormone-responsive elements. In addition, we investigated the expression patterns of family members induced by ethylene (ETH) and methyl jasmonate (MeJA), and their expression levels at different stages of T. kok-saghyz root development. Finally, subcellular localization results showed that six TkHMGR members were all located in the endoplasmic reticulum. In conclusion, the results of our study lay a certain theoretical basis for the subsequent improvement of rubber yield, molecular breeding of rubber-producing plants, and genetic improvement of T. kok-saghyz.

Bark transcriptome analyses reveals molecular mechanisms involved in tapping panel dryness occurrence and development in rubber tree (Hevea brasiliensis)

Tapping panel dryness (TPD) has become the mostimportant limiting factor for increasing natural rubber yield, whereas illuminating the molecular mechanisms underlying TPD is the prerequisite for solving the problem of TPD. However, molecular mechanisms underlying TPD are largely unknown. In this study, healthy and different stages of TPD-affected rubber trees were utilized to analyze TPD for the first time. We found that the changing tendencies of key latex physiological parameters were closely related to TPD occurrence and development. To reveal the molecular mechanisms underlying TPD, we sequenced and compared bark transcriptomes among healthy rubber tree, and TPD-affected ones at initial and advanced stages. In total, 8607 genes were identified as TPD-related genes in contrast to healthy rubber tree. According to gene expression profiles, the five samples were divided into three groups including healthy rubber tree, and TPD-affected rubber tree in the initial and advanced stages, which was consistent with the stages of TPD occurrence and development. Interestingly, only asmall proportionof the TPD-related genes were constantly down- or up-regulated with TPD occurrence and development. The TPD-related genes in KEGG pathways significantly enriched were closely associated with protein metabolism, cell division and differentiation, PCD, stress responses, terpene biosynthesis, and various metabolism processes. Moreover, overexpression of HbAPX2 identified as a TPD-related gene enhanced oxidative stress tolerance in S. cerevisiae. The typical symptoms of TPD, partial or complete dry zone (no latex flow) on tapping panel, might attribute to lower IPP available for rubber biosynthesis, and downregulation of the genes in post-IPP steps of rubber biosynthesis and the genes involved in latex flow. Our results not only provide new insights into molecular mechanisms underlying TPD occurrence and development but also contribute to developing effective measures to control TPD in rubber trees.

Transcriptome analysis of Taraxacum kok-saghyz reveals the role of exogenous methyl jasmonate in regulating rubber biosynthesis and drought tolerance

Taraxacum kok-saghyz has been identified as one of the most promising alternative rubber crops, with laticifer cells that produce high-quality rubber. To uncover the underlying molecular mechanisms regulating natural rubber biosynthesis under MeJA induction, a reference transcriptome was constructed from nine samples of T. kok-saghyz. MeJA treatment was applied for 0 h (control), 6 h, and 24 h. A total of 7452 differentially expressed genes (DEGs) were identified in response to MeJA stress, relative to the control. Functional enrichment showed that these DEGs were primarily related to hormone signaling, defensive responses, and secondary metabolism. Combined analysis of the DEGs induced by MeJA and high-expression genes in laticifer cells further identified seven DEGs related to natural rubber biosynthesis that were upregulated in latex tissue, suggesting that these candidate genes could prove valuable in studying the mechanism of MeJA-mediated natural rubber biosynthesis. In addition, 415 MeJA-responsive DEGs were from several transcription factor families associated with drought resistance. This study helps to elucidate the mechanism of natural rubber biosynthesis in T. kok-saghyz in response to MeJA stress and identifies key candidate MeJA-induced DEGs in laticifer tissue, as well as a candidate drought-response target gene, whose knowledge will promote the breeding of T. kok-saghyz in the aspect of rubber yields and quality, and drought tolerance.

Transcriptomic analysis reveals the significant effects of fertilization on the biosynthesis of sesquiterpenes in Phoebe bournei

Phoebe bournei is a potential medicinal plant. Its essential oils (Eos) are mainly composed of sesquiterpenes that has potential activities of anti-bacteria and anti-tumors. In this study, we evaluated the effects of compost and compound fertilizer on the total amount and main components of Eos in P. bournei, we also studied the molecular mechanism undergoing this process by deep sequencing the genes involved in the biosynthesis of sesquiterpenes. Fertilization enhanced the total amount of main components in Eos from both leaves and twigs. Bicyclogermacrene, the primary sesquiterpene in the leaf EO, was significantly increased under compost treatment, while bicyclogermacrene and δ-cadinene (the second most abundant sesquiterpene) were decreased under compound fertilizer treatment. The two fertilizers had no significant effect on the abundance of the primary (+) - δ-cadinene in the twig EO, but had a positive effect on the second most abundant sesquiterpene copaene. Significant differences were observed in the number of differentially expressed genes (DEGs) with the leaves showing greater number of DEGs as compared to the twigs after compost treatment. Terpenoid backbone biosynthesis (TBB) is a key pathway of sesquiterpenes synthesis. The expression of genes regulating several important enzymes in TBB was altered after fertilization. After the compost treatment, the expression of the leaf DXS gene (ACQ66107.1), being closely related to the sesquiterpene biosynthesis in P. bournei leaves, was decreased. Compost and compound fertilizer altered the expression of the two important branch-point enzymes (FPPS and GGPPS) genes (ART33314.1 and ATT59265.1), which contributed to the changes of the total amount and components of P. bournei sesquiterpenes. This study provides a new insight into the future use of P. bournei for Eos.

Multilevel interactions between native and ectopic isoprenoid pathways affect global metabolism in rice

Isoprenoids are natural products derived from isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP). In plants, these precursors are synthesized via the cytosolic mevalonate (MVA) and plastidial methylerythritol phosphate (MEP) pathways. The regulation of these pathways must therefore be understood in detail to develop effective strategies for isoprenoid metabolic engineering. We hypothesized that the strict regulation of the native MVA pathway could be circumvented by expressing an ectopic plastidial MVA pathway that increases the accumulation of IPP and DMAPP in plastids. We therefore introduced genes encoding the plastid-targeted enzymes HMGS, tHMGR, MK, PMK and MVD and the nuclear-targeted transcription factor WR1 into rice and evaluated the impact of their endosperm-specific expression on (1) endogenous metabolism at the transcriptomic and metabolomic levels, (2) the synthesis of phytohormones, carbohydrates and fatty acids, and (3) the macroscopic phenotype including seed morphology. We found that the ectopic plastidial MVA pathway enhanced the expression of endogenous cytosolic MVA pathway genes while suppressing the native plastidial MEP pathway, increasing the production of certain sterols and tocopherols. Plants carrying the ectopic MVA pathway only survived if WR1 was also expressed to replenish the plastid acetyl-CoA pool. The transgenic plants produced higher levels of fatty acids, abscisic acid, gibberellins and lutein, reflecting crosstalk between phytohormones and secondary metabolism.

Tobacco rattle virus-induced gene silencing in Hevea brasiliensis

Virus-induced gene silencing (VIGS) is a powerful gene-silencing tool that has been intensively applied in plants. To data, the application of VIGS in rubber tree has not yet been reported. In this study, we described the efficient gene silencing in rubber tree by VIGS. The gene encoding Hevea brasiliensis phytoene desaturase (HbPDS) was identified in rubber tree genome. Small interfering RNAs from HbPDS and the silencing gene fragment were predicted and a length of 399 bp was selected to be tested. We showed that the tobacco rattle virus (TRV)-VIGS could induce effective HbPDS silencing in rubber tree. This study was the first to report VIGS in rubber tree. The present TRV-VIGS method could be used to perform reverse genetic approaches to identify unknown gene functions and might be further applied to produce gene silenced rubber tree plants, to advance functional gene of rubber tree.

Strategies for the production of biochemicals in bioenergy crops

Industrial crops are grown to produce goods for manufacturing. Rather than food and feed, they supply raw materials for making biofuels, pharmaceuticals, and specialty chemicals, as well as feedstocks for fabricating fiber, biopolymer, and construction materials. Therefore, such crops offer the potential to reduce our dependency on petrochemicals that currently serve as building blocks for manufacturing the majority of our industrial and consumer products. In this review, we are providing examples of metabolites synthesized in plants that can be used as bio-based platform chemicals for partial replacement of their petroleum-derived counterparts. Plant metabolic engineering approaches aiming at increasing the content of these metabolites in biomass are presented. In particular, we emphasize on recent advances in the manipulation of the shikimate and isoprenoid biosynthetic pathways, both of which being the source of multiple valuable compounds. Implementing and optimizing engineered metabolic pathways for accumulation of coproducts in bioenergy crops may represent a valuable option for enhancing the commercial value of biomass and attaining sustainable lignocellulosic biorefineries.

The 3-hydroxy-3-methylglutaryl-coenzyme A reductase 5 gene from Malus domestica enhances oxidative stress tolerance in Arabidopsis thaliana

3-Hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) is the first rate-limiting enzyme regulating the synthesis of terpenoids upstream of the mevalonate (MVA) pathway. In higher plants, members of the HMGR genes families play an important role in plant growth and development and in response to various environmental stresses. In the present study, a novel HMGR gene, designated MdHMGR5, was isolated from apple (Malus domestica L.) and characterized. Expression of MdHMGR5 enhanced the activity of HMGR enzyme in transgenic Arabidopsis thaliana L. plants. Under oxidative stress, transgenic A. thaliana plants over-expressing MdHMGR5 had a higher germination rate, a longer main root length, higher chlorophyll and proline content, and higher activities of antioxidant enzymes. On the other hand, malondialdehyde (MDA) content, relative conductivity and reactive oxygen species (ROS) production rate were significantly lower than in wild type plants. These results indicated that over-expression of MdHMGR5 enhanced plant tolerance to oxidative stress by scavenging ROS in transgenic plants. Over-expression of MdHMGR5 also affected the expression levels of genes in mevalonic acid and 2C-methyl-D-erythritol 4-phosphate (MVA and MEP) pathways of A. thaliana plants. These results indicate that over-expression of MdHMGR5 enhances tolerance to oxidative stress by maintaining photosynthesis and scavenging ROS in transgenic A. thaliana plants.

Isolation, characterization and in silico analysis of 3-Hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) gene from Andrographis paniculata (Burm. f) Nees

3-Hydroxy-3-methylglutaryl-coenzymeA reductase (HMGR), the first rate-limiting enzyme of Mevalonate (MVA) pathway was isolated from Andrographis paniculata (ApHMGR) and expressed in bacterial cells. Full length ApHMGR (1937 bp) was submitted to NCBI with accession number MG271748.1. The open reading frame (ORF) was flanked by a 31-bp 5'-UTR, 118-bp 3'-UTR and ApHMGR contained a 1787 bp ORF encoding protein of 595 amino acids. ApHMGR protein was approximately 64 kDa, with isoelectric point of 5.75. Isolated ApHMGR was cloned into pET102 vector and expressed in E. coli BL21 (DE 3) cells, and characterized by SDS-PAGE. HPLC analysis for andrographolide content in leaf, stem and root of A. paniculata revealed highest in leaf tissue. The expression patterns of ApHMGR in different plant tissues using qRT-PCR revealed high in root tissue correlating with HPLC data. Three dimensional (3D) structural model of ApHMGR displayed 90% of the amino acids in most favored regions of the Ramachandran plot with 93% overall quality factor. ApHMGR was highly conserved with plant specific N-terminal membrane domains and C-terminal catalytic regions. Phylogenetic analysis showed A. paniculata sharing common ancestor with Handroanthus impetiginosus. 3D model of ApHMGR was screened for the interaction with substrates NADPH, HMG CoA and inhibitor using Auto Dock Vina. In silico analysis revealed that full length ApHMGR had extensive similarities to other plant HMGRs. The present communication reports the isolation of full length HMGR from A. paniculata, its heterologous expression in bacterial cells and in silico structural and functional characterization providing valuable genomic information for future molecular interventions.

cDNA cloning, prokaryotic expression and functional analysis of 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR) in Pogostemon cablin

Pogostemon cablin is an important commercial source of patchouli oil, whose main active ingredient is patchouli alcohol. This sesquiterpene is a product of the mevalonate pathway, in which 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (HMGCR) is the rate-limiting enzyme. In this study, P. cablin HMGCR cDNA, comprising 2209 nucleotides encoding 425 amino acid residues was isolated, and bioinformatics analysis was used to analyze the protein sequence. Based on this analysis, a C-terminal truncated variant was engineered for recombinant expression in E. coli. The 38 kDa recombinant protein was identified by SDS-PAGE, and assayed for mevalonolactone production. According to the PcHMGCR1 gene sequence alignment with other species, the HMGCR protein had obvious resemblance with other plants HMG coenzyme A reductase and had homology with other species including plants, fungi, archaebacteria and animals. The prokaryotic expression vector was constructed by restriction enzyme digestion to be transformed into E. coli to express the recombinant protein, and 38 kDa recombinant protein was identified by the SDS-PAGE. Enzymatic activity was detected using GC-MS and, as a result, mevalonolactone was detected in the in vitro reaction mixture. Differential expression analysis showed that PcHMGCR1 expressed the highest amount in roots. The research results are of great significance for further research on the molecular biosynthesis mechanism of Patchouli alcohol in P. cablin.

Identification, expression, and phylogenetic analyses of terpenoid biosynthesis-related genes in secondary xylem of loblolly pine (Pinus taeda L.) based on transcriptome analyses

Loblolly pine (Pinus taeda L.) is one of the most important species for oleoresin (a mixture of terpenoids) in South China. The high oleoresin content of loblolly pine is associated with resistance to bark beetles and other economic benefits. In this study, we conducted transcriptome analyses of loblolly pine secondary xylem to gain insight into the genes involved in terpenoid biosynthesis. A total of 372 unigenes were identified as being critical for oleoresin production, including genes for ATP-binding cassette (ABC) transporters, the cytochrome P450 (CYP) protein family, and terpenoid backbone biosynthesis enzymes. Six key genes involved in terpenoid biosynthetic pathways were selected for multiple sequence alignment, conserved motif prediction, and phylogenetic and expression profile analyses. The protein sequences of all six genes exhibited a higher degree of sequence conservation, and upstream genes were relatively more conserved than downstream genes in terpenoid biosynthetic pathways. The N-terminal regions of these sequences were less conserved than the C-terminal ends, as the N-terminals were quite diverse in both length and composition. The phylogenetic analyses revealed that most genes originated from gene duplication after species divergence, and partial genes exhibited incomplete lineage sorting. In addition, the expression profile analyses showed that all six genes exhibited high expression levels during the high-oleoresin-yielding phase.