Impact of an intragenic retrotransposon on the structural integrity and evolution of a major isoprenoid biosynthesis pathway gene in Hevea brasiliensis

Development of Novel Markers for Yield in Hevea brasiliensis Muell. Arg. Based on Candidate Genes from Biosynthetic Pathways Associated with Latex Production

Scarcity of functional genetic markers associated with candidate genes (CGs) is a serious constraint for marker-assisted selection in the natural rubber producing tree, Hevea brasiliensis. In order to develop markers associated with rubber yield, five CGs involved in latex biosynthesis were characterized from 16 popular Hevea varieties. Novel SNPs and indels were identified and developed into markers using simple genotyping techniques like allele-specific PCR, CAPS, etc. A progeny population was genotyped using these markers to validate them, to understand their segregation pattern and to map them to a genetic linkage map. Parent-specific maps were constructed using pseudo-test cross strategy with the help of additional markers. The sequence structure information generated will be useful for future studies on gene mapping, functional relevance of coding SNPs and evolution of rubber biosynthesis genes in Hevea. Concurrently, the markers developed may serve as powerful tools for yield-based selection and for genetic diversity and pedigree studies in Hevea. Above all, the marker assays designed for genotyping could be economically carried out in any laboratory having basic molecular biology infrastructure and expertise.

Metabolite Profiling and Transcriptome Analysis Explains Difference in Accumulation of Bioactive Constituents in Licorice (Glycyrrhiza uralensis) Under Salt Stress

Salinity stress significantly affects the contents of bioactive constituents in licorice Glycyrrhiza uralensis. To elucidate the molecular mechanism underlying the difference in the accumulation of these constituents under sodium chloride (NaCl, salt) stress, licorice seedlings were treated with NaCl and then subjected to an integrated transcriptomic and metabolite profiling analysis. The transcriptomic analysis results identified 3,664 differentially expressed genes (DEGs) including transcription factor family MYB and basic helix-loop-helix (bHLH). Most DEGs were involved in flavonoid and terpenoid biosynthesis pathways. In addition, 121 compounds including a triterpenoid and five classes of flavonoids (isoflavone, flavone, flavanone, isoflavan, and chalcone) were identified, and their relative levels were compared between the stressed and control groups using data from the ultrafast liquid chromatography (UFLC)-triple quadrupole-time of flight-tandem mass spectrometry (TOF-MS/MS) analysis. Putative biosynthesis networks of the flavonoids and triterpenoids were created and combined with structural DEGs such as phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase [4CL], cinnamate 4-hydroxylase [C4H], chalcone synthase [CHS], chalcone-flavanone isomerase [CHI], and flavonoid-3',5' hydroxylase (F3',5'H) for flavonoids, and CYP88D6 and CYP72A154 for glycyrrhizin biosynthesis. Notably, significant upregulation of UDP-glycosyltransferase genes (UGT) in salt-stressed licorice indicated that postmodification of glycosyltransferase may participate in downstream biosynthesis of flavonoid glycosides and triterpenoid saponins. Accordingly, the expression trend of the DEGs is positively correlated with the accumulation of glycosides. Our study findings indicate that key DEGs and crucial UGT genes co-regulate flavonoid and saponin biosynthesis in licorice under salt stress.

Structural and Functional Annotation of Transposable Elements Revealed a Potential Regulation of Genes Involved in Rubber Biosynthesis by TE-Derived siRNA Interference in Hevea brasiliensis

The natural rubber biosynthetic pathway is well described in Hevea, although the final stages of rubber elongation are still poorly understood. Small Rubber Particle Proteins and Rubber Elongation Factors (SRPPs and REFs) are proteins with major function in rubber particle formation and stabilization. Their corresponding genes are clustered on a scaffold1222 of the reference genomic sequence of the Hevea brasiliensis genome. Apart from gene expression by transcriptomic analyses, to date, no deep analyses have been carried out for the genomic environment of SRPPs and REFs loci. By integrative analyses on transposable element annotation, small RNAs production and gene expression, we analysed their role in the control of the transcription of rubber biosynthetic genes. The first in-depth annotation of TEs (Transposable Elements) and their capacity to produce TE-derived siRNAs (small interfering RNAs) is presented, only possible in the Hevea brasiliensis clone PB 260 for which all data are available. We observed that 11% of genes are located near TEs and their presence may interfere in their transcription at both genetic and epigenetic level. We hypothesized that the genomic environment of rubber biosynthesis genes has been shaped by TE and TE-derived siRNAs with possible transcriptional interference on their gene expression. We discussed possible functionalization of TEs as enhancers and as donors of alternative transcription start sites in promoter sequences, possibly through the modelling of genetic and epigenetic landscapes.

Comparative Proteomic Analysis Reveals the Molecular Mechanisms Underlying the Accumulation Difference of Bioactive Constituents in Glycyrrhiza uralensis Fisch under Salt Stress

Licorice (Glycyrrhiza uralensis Fisch) possesses a substantial share of the global markets for its unique sweet flavor and diverse pharmacological compounds. Cultivated licorice is widely distributed in northwest regions of China, covered with land with a broad range of salinities. A preliminary study indicated that suitable salt stress significantly increased the content of bioactive constituents in licorice. However, the molecular mechanisms underlying the influence of salinity on the accumulation of these constituents remain unclear, which hinders quality breeding of cultivated licorice. In our study, flavonoid-related structural genes were obtained, and most of them, such as phenylalanine ammonia-lyases, cinnamate 4-hydroxylases, 4-coumarate: CoA ligases, chalcone synthases, chalcone-flavanone isomerase, and flavonol synthase, showed high levels after salt treatment. In the biosynthesis of glycyrrhizin, three key enzymes (bAS, CYP88D6, and CYP72A154) were identified as differentially expressed proteins and remarkably upregulated in the salt-stressed group. Combining these results with the contents of 14 bioactive constituents, we also found that the expression patterns of those structural proteins were logically consistent with changes in bioactive constituent profiles. Thus, we believe that suitable salt stress increased the accumulation of bioactive constituents in licorice by upregulating proteins involved in the related biosynthesis pathways. This work provided valuable proteomic information for unraveling the molecular mechanism of flavonoid and glycyrrhizin metabolism and offered fundamental resources for quality breeding in licorice.

Characterization and Function of 3-Hydroxy-3-Methylglutaryl-CoA Reductase in Populus trichocarpa: Overexpression of PtHMGR Enhances Terpenoids in Transgenic Poplar

In the mevalonic acid (MVA) pathway, 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) is considered the first rate-limiting enzyme in isoprenoid biosynthesis. In this study, we cloned a full-length cDNA from Populus trichocarpa with an open reading frame of 1,734 bp. The deduced PtHMGR sequence contained two HMG-CoA motifs and two NADPH motifs, which exhibited homology with HMGR proteins from other species. Subsequently, truncated PtHMGR was expressed in Escherichia coli BL21 (DE3) cells, and enzyme activity analysis revealed that the truncated PtHMGR protein could catalyze the reaction of HMG-CoA and NADPH to form MVA. Relative expression analysis suggests that PtHMGR expression varies among tissues and that PtHMGR responds significantly to abscisic acid (ABA), NaCl, PEG6000, hydrogen peroxide (H2O2), and cold stresses. We used polymerase chain reaction (PCR) analysis to select transgenic Nanlin 895 poplars (Populus× euramericana cv.) and quantitative reverse-transcription PCR (qRT-PCR) to show that PtHMGR expression levels were 3- to 10-fold higher in transgenic lines than in wild-type (WT) poplars. qRT-PCR was also used to determine transcript levels of methylerythritol phosphate (MEP)-, MVA-, and downstream-related genes, indicating that overexpression of PtHMGR not only affects expression levels of MVA-related genes, but also those of MEP-related genes. We also measured the content of terpenoids including ABA, gibberellic acid (GA), carotenes, and lycopene. PtHMGR overexpression significantly increased ABA, GA, carotene, and lycopene content, indicating that PtHMGR participates in the regulation of terpenoid compound synthesis.

Proteomic Landscape Has Revealed Small Rubber Particles Are Crucial Rubber Biosynthetic Machines for Ethylene-Stimulation in Natural Rubber Production

Rubber particles are a specific organelle for natural rubber biosynthesis (NRB) and storage. Ethylene can significantly improve rubber latex production by increasing the generation of small rubber particles (SRPs), regulating protein accumulation, and activating many enzyme activities. We conducted a quantitative proteomics study of different SRPs upon ethylene stimulation by differential in-gel electrophoresis (DIGE) and using isobaric tags for relative and absolute quantification (iTRAQ) methods. In DIGE, 79 differentially accumulated proteins (DAPs) were determined as ethylene responsive proteins. Our results show that the abundance of many NRB-related proteins has been sharply induced upon ethylene stimulation. Among them, 23 proteins were identified as rubber elongation factor (REF) and small rubber particle protein (SRPP) family members, including 16 REF and 7 SRPP isoforms. Then, 138 unique phosphorylated peptides, containing 129 phosphorylated amino acids from the 64 REF/SRPP family members, were identified, and most serine and threonine were phosphorylated. Furthermore, we identified 226 DAPs from more than 2000 SRP proteins by iTRAQ. Integrative analysis revealed that almost all NRB-related proteins can be detected in SRPs, and many proteins are positively responsive to ethylene stimulation. These results indicate that ethylene may stimulate latex production by regulating the accumulation of some key proteins. The phosphorylation modification of REF and SRPP isoforms might be crucial for NRB, and SRP may act as a complex natural rubber biosynthetic machine.

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.

Distinguishing CPT gene family members and vetting the sequence structure of a putative rubber synthesizing variant in Hevea brasiliensis

cis-Prenyltransferases (cis-PTs) constitute a large family of enzymes conserved during evolution and present in all domains of life. cis-PTs catalyze the cis-1,4-polymerization of isoprene units to generate isoprenoids with carbon skeletons varying from C10 (neryl pyrophosphate) to C > 10,000 (natural rubber). Though the previously reported CPTs in Hevea are designated based on sequence variations, their classification was done mostly by phylogenetic analysis using a mixture of partial as well as full length sequences often excluding the UTRs. In this context an attempt was made to reclassify the CPTs strictly based on their sequence similarity and distinguish the members putatively associated with rubber biosynthesis from the others. Extensive computational analysis was carried out on CPT sequences obtained from public resources and whole genome assemblies of Hevea. Based on the results from BLAST analysis, multiple sequence alignments of protein, nucleotide and untranslated regions, open reading frame analysis, gene prediction analysis and sequence length variations, we conclude that there exists mainly three CPTs namely RubCPT1, RubCPT2 and RubCPT3 putatively associated with rubber biosynthesis in Hevea brasiliensis. The rest were categorised as variants of dehydrodolichyl diphosphate synthase (DHDDS) involved in the synthesis of dolichols having short chain isoprenoids. Analysis of the sequence structure of the most highly expressed RubCPT1 in latex revealed the allele richness and diversity of this important variant prevailing in the popular rubber clones. Haplotypes consisting of SNPs with high degree of heterozygosity were also identified. Segregation and linkage disequilibrium analysis confirmed that recombination is the major contributor towards the generation of allelic diversity rather than point mutations. Alternatively, gene expression analysis indicated the possibility of association between specific haplotypes and RubCPT1 expression in Hevea clones which may have downstream impact up to the level of rubber production. The conclusions from this study may pave way for the identification and better understanding of CPTs directly involved with natural rubber biosynthesis in Hevea and the SNP data generated may aid in the development of molecular markers putatively associated with yield in rubber.

The relationship between latex metabolism gene expression with rubber yield and related traits in Hevea brasiliensis

Background

Expression patterns of many laticifer-specific gens are closely correlative with rubber yield of Hevea brasiliensis (para rubber tree). To unveil the mechanisms underlying the rubber yield, transcript levels of nine major latex metabolism-related genes, i.e., HMG-CoA synthase (HMGS), HMG-CoA reductase (HMGR), diphosphomevalonate decarboxylase (PMD), farnesyl diphosphate synthase (FPS), cis-prenyltransferase (CPT), rubber elongation factor (REF), small rubber particle protein (SRPP), dihydroxyacid dehydratase (DHAD) and actin depolymerizing factor (ADF), were dertermined, and the relationship between rubber yield with their expression levels was analysed.

Results

Except HbHMGR1, HbPMD and HbDHAD, most of these genes were predominantly expressed in latex, and bark tapping markedly elevated the transcript abundance of the analyzed genes, with the 7th tapping producing the greatest expression levels. Both ethephon (ETH) and methyl jasmonate (MeJA) stimulation greatly induced the expression levels of the examined genes, at least at one time point, except HbDHAD, which was unresponsive to MeJA. The genes’ expression levels, as well as the rubber yields and two yield characteristics differed significantly among the different genotypes examined. Additionally, the latex and dry rubber yields increased gradually but the dry rubber content did not. Rubber yields and/or yield characteristics were significantly positively correlated with HbCPT, HbFPS, HbHMGS, HbHMGR1 and HbDHAD expression levels, negatively correlated with that of HbREF, but not significantly correlated with HbPMD, HbSRPP and HbADF expression levels. In addition, during rubber production, significantly positive correlations existed between the expression level of HbPMD and the levels of HbREF and HbHMGR1, between HbSRPP and the levels of HbHMGS and HbHMGR1, and between HbADF and HbFPS.

Conclusions

The up-regulation of these genes might be related to the latex production of rubber trees under the stress of bark tapping and latex metabolism. The various correlations among the genes implied that there are differences in their synergic interactions. Thus, these nine genes might be related to rubber yield and yield-related traits in H. brasiliensis, and this work increases our understanding of their complex functions and how they are expressed in both high-and medium-yield rubber tree varieties and low-yield wild rubber tree germplasm.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5242-4) contains supplementary material, which is available to authorized users.

Subcellular proteome profiles of different latex fractions revealed washed solutions from rubber particles contain crucial enzymes for natural rubber biosynthesis

Rubber particle (RP) is a specific organelle for natural rubber biosynthesis (NRB) and storage in rubber tree Hevea brasiliensis. NRB is processed by RP membrane-localized proteins, which were traditionally purified by repeated washing. However, we noticed many proteins in the discarded washing solutions (WS) from RP. Here, we compared the proteome profiles of WS, C-serum (CS) and RP by 2-DE, and identified 233 abundant proteins from WS by mass spectrometry. Many spots on 2-DE gels were identified as different protein species. We further performed shotgun analysis of CS, WS and RP and identified 1837, 1799 and 1020 unique proteins, respectively. Together with 2-DE, we finally identified 1825 proteins from WS, 246 were WS-specific. These WS-specific proteins were annotated in Gene Ontology, indicating most abundant pathways are organic substance metabolic process, protein degradation, primary metabolic process, and energy metabolism. Protein-protein interaction analysis revealed these WS-specific proteins are mainly involved in ribosomal metabolism, proteasome system, vacuolar protein sorting and endocytosis. Label free and Western blotting revealed many WS-specific proteins and protein complexes are crucial for NRB initiation. These findings not only deepen our understanding of WS proteome, but also provide new evidences on the roles of RP membrane proteins in NRB.

SIGNIFICANCE

Natural rubber is stored in rubber particle from the rubber tree. Rubber particles were traditionally purified by repeated washing, but many proteins were identified from the washing solutions (WS). We obtained the first visualization proteome profiles with 1825 proteins from WS, including 246 WS-specific ones. These WS proteins contain almost all enzymes for polyisoprene initiation and may play important roles in rubber biosynthesis.

Molecular evolution and functional characterisation of haplotypes of an important rubber biosynthesis gene in Hevea brasiliensis

Hydroxy-methylglutaryl coenzyme-A synthase (HMGS) is a rate-limiting enzyme in the cytoplasmic isoprenoid biosynthesis pathway leading to natural rubber production in Hevea brasiliensis (rubber). Analysis of the structural variants of this gene is imperative to understand their functional significance in rubber biosynthesis so that they can be properly utilised for ongoing crop improvement programmes in Hevea. We report here allele richness and diversity of the HMGS gene in selected popular rubber clones. Haplotypes consisting of single nucleotide polymorphisms (SNPs) from the coding and non-coding regions with a high degree of heterozygosity were identified. Segregation and linkage disequilibrium analysis confirmed that recombination is the major contributor to the generation of allelic diversity, rather than point mutations. The evolutionarily conserved nature of some SNPs was identified by comparative DNA sequence analysis of HMGS orthologues from diverse taxa, demonstrating the molecular evolution of rubber biosynthesis genes in general. In silico three-dimensional structural studies highlighting the structural positioning of non-synonymous SNPs from different HMGS haplotypes revealed that the ligand-binding site on the enzyme remains impervious to the reported sequence variations. In contrast, gene expression results indicated the possibility of association between specific haplotypes and HMGS expression in Hevea clones, which may have a downstream impact up to the level of rubber production. Moreover, haplotype diversity of the HMGS gene and its putative association with gene expression can be the basis for further genetic association studies in rubber. Furthermore, the data also show the role of SNPs in the evolution of candidate genes coding for functional traits in plants.

Characterization and function of the 3-hydroxy-3-methylglutaryl-CoA reductase gene in Alisma orientale (Sam.) Juz. and its relationship with protostane triterpene production

Protostane triterpenes from Alisma orientale (Sam.) Juz. have exhibited distinct pharmacological properties that are currently in high demand. 3-Hydroxy-3-methylglutaryl-CoA reductase (HMGR) is considered the first rate-limiting enzyme in isoprenoid biosynthesis via the mevalonic acid (MVA) pathway. In this study, we cloned a full-length cDNA of A. orientale (Sam.) Juz. HMGR (AoHMGR; 2252 bp; GenBank accession no. KP342318) with an open reading frame (ORF) of 1809 bp. The deduced protein sequence contained four conserved motifs and exhibited homology with HMGR proteins from other plants. We next expressed the cloned gene in Escherichia coli BL21 (Rosetta) cells, collected the expressed products, and incubated those with 3-hydroxy-3-methylglutaryl-CoA (HMG-CoA) to determine enzymatic activity. GC/MS analysis revealed that the products were able to catalyze HMG-CoA and NADPH to form MVA. The purified protein was used to immunize New Zealand rabbits and prepare an antibody against AoHMGR. Western blot results demonstrated that the antibodies specifically recognized AoHMGR protein in A. orientale (Sam.) Juz. We then established a rapid test to detect AoHMGR protein in the plant, and found the tuber to be the most AoHMGR protein-abundant organ in A. orientale (Sam.) Juz. Furthermore, we detected the expression level of AoHMGR and contents of the main active component, Alisol B 23-acetate, at different growth phases of A. orientale (Sam.) Juz. A significant positive correlation was identified, indicating that AoHMGR represents a key enzyme in the synthetic pathway of protostane triterpenes.