Comparative analysis of latex transcriptome reveals putative molecular mechanisms underlying super productivity of 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.

Transcript and Protein Profiling Provides Insights Into the Molecular Mechanisms of Harvesting-Induced Latex Production in Rubber Tree

Natural rubber, an important industrial raw material with wide applications, is harvested in the form of latex (cytoplasm of rubber-producing laticifers) from Hevea brasiliensis (para rubber tree) by the way of tapping. Conspicuous stimulation on latex production is observed for the first few tappings conducted on virgin (untapped before) or resting (tapped before but no tapping for a period) rubber trees. To understand the underlying mechanisms, an integrative analysis of the latex transcriptome and proteome was conducted on virgin or resting Hevea trees for the first five tappings. A total of 505 non-redundant differentially expressed (DE) transcript-derived fragments (TDFs) were identified by silver-staining cDNA-AFLP, with 217 exhibiting patterns of upregulated, 180 downregulated and 108 irregularly-regulated. Meanwhile, 117 two dimensional gel electrophoresis DE-protein spots were isolated and subjected to mass spectrometry analysis, with 89 and 57 being successfully identified by MALDI-TOF and MALDI-TOF/TOF, respectively. About 72.5% DE-TDFs and 76.1% DE-proteins were functionally annotated and categorized. Noteworthily, most of the DE-TDFs implicated in sugar transport and metabolism as well as rubber biosynthesis were upregulated by the tapping treatment. The importance of sugar metabolism in harvesting-induced latex production was reinforced by the identification of abundant relevant DE-protein spots. About 83.8% of the randomly selected DE-TDFs were validated for expression patterns by semi-quantitative RT-PCR, and an 89.7% consistency for the 29 latex regeneration-related DE-TDFs examined by quantitative RT-PCR analysis. In brief, our results reveal extensive physiological and molecular changes in Hevea laticifers incurred by the tapping treatment, and the vast number of DE genes and proteins identified here contribute to unraveling the gene regulatory network of tapping-stimulated latex production.

Identification and Characterization of Glycoproteins and Their Responsive Patterns upon Ethylene Stimulation in the Rubber Latex

Natural rubber is an important industrial material, which is obtained from the only commercially cultivated rubber tree, Hevea brasiliensis. In rubber latex production, ethylene has been extensively used as a stimulant. Recent research showed that post-translational modifications (PTMs) of latex proteins, such as phosphorylation, glycosylation and ubiquitination, are crucial in natural rubber biosynthesis. In this study, comparative proteomics was performed to identify the glycosylated proteins in rubber latex treated with ethylene for different days. Combined with Pro-Q Glycoprotein gel staining and mass spectrometry techniques, we provided the first visual profiling of glycoproteomics of rubber latex and finally identified 144 glycosylated protein species, including 65 differentially accumulated proteins (DAPs) after treating with ethylene for three and/or five days. Gene Ontology (GO) functional annotation showed that these ethylene-responsive glycoproteins are mainly involved in cell parts, membrane components and metabolism. Pathway analysis demonstrated that these glycosylated rubber latex proteins are mainly involved in carbohydrate metabolism, energy metabolism, degradation function and cellular processes in rubber latex metabolism. Protein-protein interaction analysis revealed that these DAPs are mainly centered on acetyl-CoA acetyltransferase and hydroxymethylglutaryl-CoA synthase (HMGS) in the mevalonate pathway for natural rubber biosynthesis. In our glycoproteomics, three protein isoforms of HMGS2 were identified from rubber latex, and only one HMGS2 isoform was sharply increased in rubber latex by ethylene treatment for five days. Furthermore, the HbHMGS2 gene was over-expressed in a model rubber-producing grass Taraxacum Kok-saghyz and rubber content in the roots of transgenic rubber grass was significantly increased over that in the wild type plant, indicating HMGS2 is the key component for natural rubber production.

Identification and evaluation of suitable reference genes for gene expression analysis in rubber tree leaf

Gene expression profiles are increasingly applied to investigate molecular mechanism for which, normalization with suitable reference genes is critical. Previously we have reported several suitable reference genes for laticifer samples from rubber tree, however, little is known in leaf. The main objective of this current study was to identify some stable expression reference genes at various developmental stages of leaf, as well as during abiotic (high and low temperature extremes) and biotic stresses (pathogen stress). Gene expression profilings identified the ubiquitin-proteasome system as excellent potential as reference genes for rubber tree leaf. Among a total of 30 tested genes investigated, 24 new candidate (including 11 genes involved in the ubiquitin-proteasome system), 4 previously identified and 2 specific genes, were further evaluated using quantitative real-time PCR. Our results indicated that the new candidate genes had better expression stability comparing with others. For instance, an ubiquitin conjugating enzyme (RG0099) and three ubiquitin-protein ligases (RG0928, RG2190 and RG0118) expressed stably in all samples, and were confirmed to be suitable reference genes for rubber tree leaf under four different conditions. Finally, we suggest that using more than one reference gene may be appropriate in gene expression studies when employing different software to normalize gene expression data. Our findings have significant implications for the reliability of data obtained from genomics studies in rubber tree and perhaps in other species.

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.

Reactive oxygen species in Hevea brasiliensis latex and relevance to Tapping Panel Dryness

Environmental stress can lead to oxidative stress resulting from an increase in reactive oxygen species (ROS) and involves redox adjustments. Natural rubber is synthesized in laticifers, which is a non-photosynthetic tissue particularly prone to oxidative stress. This paper reviews the current state of knowledge on the ROS production and ROS-scavenging systems in laticifers. These regulations have been the subject of intense research into a physiological syndrome, called Tapping Panel Dryness (TPD), affecting latex production in Hevea brasiliensis. In order to prevent TPD occurrence, monitoring thiol content appeared to be a crucial factor of latex diagnosis. Thiols, ascorbate and γ-tocotrienol are the major antioxidants in latex. They are involved in membrane protection from ROS and likely have an effect on the quality of raw rubber. Some transcription factors might play a role in the redox regulatory network in Hevea, in particular ethylene response factors, which have been the most intensively studied given the role of ethylene on rubber production. Current challenges for rubber research and development with regard to redox systems will involve improving antioxidant capacity using natural genetic variability.

Identification of differentially expressed genes and signalling pathways in bark of Hevea brasiliensis seedlings associated with secondary laticifer differentiation using gene expression microarray

The natural rubber of Para rubber tree, Hevea brasiliensis, is the main crop involved in industrial rubber production due to its superior quality. The Hevea bark is commercially exploited to obtain latex, which is produced from the articulated secondary laticifer. The laticifer is well defined in the aspect of morphology; however, only some genes associated with its development have been reported. We successfully induced secondary laticifer in the jasmonic acid (JA)-treated and linolenic acid (LA)-treated Hevea bark but secondary laticifer is not observed in the ethephon (ET)-treated and untreated Hevea bark. In this study, we analysed 27,195 gene models using NimbleGen microarrays based on the Hevea draft genome. 491 filtered differentially expressed (FDE) transcripts that are common to both JA- and LA-treated bark samples but not ET-treated bark samples were identified. In the Eukaryotic Orthologous Group (KOG) analysis, 491 FDE transcripts belong to different functional categories that reflect the diverse processes and pathways involved in laticifer differentiation. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) and KOG analysis, the profile of the FDE transcripts suggest that JA- and LA-treated bark samples have a sufficient molecular basis for secondary laticifer differentiation, especially regarding secondary metabolites metabolism. FDE genes in this category are from the cytochrome (CYP) P450 family, ATP-binding cassette (ABC) transporter family, short-chain dehydrogenase/reductase (SDR) family, or cinnamyl alcohol dehydrogenase (CAD) family. The data includes many genes involved in cell division, cell wall synthesis, and cell differentiation. The most abundant transcript in FDE list was SDR65C, reflecting its importance in laticifer differentiation. Using the Basic Local Alignment Search Tool (BLAST) as part of annotation and functional prediction, several characterised as well as uncharacterized transcription factors and genes were found in the dataset. Hence, the further characterization of these genes is necessary to unveil their role in laticifer differentiation. This study provides a platform for the further characterization and identification of the key genes involved in secondary laticifer differentiation.

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.

The Hevea brasiliensis XIP aquaporin subfamily: genomic, structural and functional characterizations with relevance to intensive latex harvesting

X-Intrinsic Proteins (XIP) were recently identified in a narrow range of plants as a full clade within the aquaporins. These channels reportedly facilitate the transport of a wide range of hydrophobic solutes. The functional roles of XIP in planta remain poorly identified. In this study, we found three XIP genes (HbXIP1;1, HbXIP2;1 and HbXIP3;1) in the Hevea brasiliensis genome. Comprehensive bioinformatics, biochemical and structural analyses were used to acquire a better understanding of this AQP subfamily. Phylogenetic analysis revealed that HbXIPs clustered into two major groups, each distributed in a specific lineage of the order Malpighiales. Tissue-specific expression profiles showed that only HbXIP2;1 was expressed in all the vegetative tissues tested (leaves, stem, bark, xylem and latex), suggesting that HbXIP2;1 could take part in a wide range of cellular processes. This is particularly relevant to the rubber-producing laticiferous system, where this isoform was found to be up-regulated during tapping and ethylene treatments. Furthermore, the XIP transcriptional pattern is significantly correlated to latex production level. Structural comparison with SoPIP2;1 from Spinacia oleracea species provides new insights into the possible role of structural checkpoints by which HbXIP2;1 ensures glycerol transfer across the membrane. From these results, we discuss the physiological involvement of glycerol and HbXIP2;1 in water homeostasis and carbon stream of challenged laticifers. The characterization of HbXIP2;1 during rubber tree tapping lends new insights into molecular and physiological response processes of laticifer metabolism in the context of latex exploitation.

Molecular mechanism of ethylene stimulation of latex yield in rubber tree (Hevea brasiliensis) revealed by de novo sequencing and transcriptome analysis

Background

Rubber tree (Hevea brasiliensis) is an important industrial crop cultivated in tropical areas for natural rubber production. Treatment of the bark of rubber trees with ehephon (an ethylene releaser) has been a routine measure to increase latex yield, but the molecular mechanism behind the stimulation of rubber production by ethylene still remains a puzzle. Deciphering the enigma is of great importance for improvement of rubber tree for high yield.

Results

De novo sequencing and assembly of the bark transciptomes of Hevea brasiliensis induced with ethephon for 8 h (E8) and 24 h (E24) were performed. 51,965,770, 52,303,714 and 53,177,976 high-quality clean reads from E8, E24 and C (control) samples were assembled into 81,335, 80,048 and 80,800 unigenes respectively, with a total of 84,425 unigenes and an average length of 1,101 bp generated. 10,216 and 9,374 differentially expressed genes (DEGs) in E8 and E24 compared with C were respectively detected. The expression of several enzymes in crucial points of regulation in glycolysis were up-regulated and DEGs were not significantly enriched in isopentenyl diphosphate (IPP) biosynthesis pathway. In addition, up-regulated genes of great regulatory importance in carbon fixation (Calvin cycle) were identified.

Conclusions

The rapid acceleration of glycolytic pathway supplying precursors for the biosynthesis of IPP and natural rubber, instead of rubber biosynthesis per se, may be responsible for ethylene stimulation of latex yield in rubber tree. The elevated rate of flux throughout the Calvin cycle may account for some durability of ethylene-induced stimulation. Our finding lays the foundations for molecular diagnostic and genetic engineering for high-yielding improvement of rubber tree.

Electronic supplementary material

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

Profiling Ethylene-Responsive Genes Expressed in the Latex of the Mature Virgin Rubber Trees Using cDNA Microarray

Ethylene is commonly used as a latex stimulant of Hevea brasiliensis by application of ethephon (chloro-2-ethylphosphonic acid); however, the molecular mechanism by which ethylene increases latex production is not clear. To better understand the effects of ethylene stimulation on the laticiferous cells of rubber trees, a latex expressed sequence tag (EST)-based complementary DNA microarray containing 2,973 unique genes (probes) was first developed and used to analyze the gene expression changes in the latex of the mature virgin rubber trees after ethephon treatment at three different time-points: 8, 24 and 48 h. Transcript levels of 163 genes were significantly altered with fold-change values ≥ 2 or ≤ –2 (q-value < 0.05) in ethephon-treated rubber trees compared with control trees. Of the 163 genes, 92 were up-regulated and 71 down-regulated. The microarray results were further confirmed using real-time quantitative reverse transcript-PCR for 20 selected genes. The 163 ethylene-responsive genes were involved in several biological processes including organic substance metabolism, cellular metabolism, primary metabolism, biosynthetic process, cellular response to stimulus and stress. The presented data suggest that the laticifer water circulation, production and scavenging of reactive oxygen species, sugar metabolism, and assembly and depolymerization of the latex actin cytoskeleton might play important roles in ethylene-induced increase of latex production. The results may provide useful insights into understanding the molecular mechanism underlying the effect of ethylene on latex metabolism of H. brasiliensis.

Sequencing Crop Genomes: A Gateway to Improve Tropical Agriculture

Agricultural development in the tropics lags behind development in the temperate latitudes due to the lack of advanced technology, and various biotic and abiotic factors. To cope with the increasing demand for food and other plant-based products, improved crop varieties have to be developed. To breed improved varieties, a better understanding of crop genetics is necessary. With the advent of next-generation DNA sequencing technologies, many important crop genomes have been sequenced. Primary importance has been given to food crops, including cereals, tuber crops, vegetables, and fruits. The DNA sequence information is extremely valuable for identifying key genes controlling important agronomic traits and for identifying genetic variability among the cultivars. However, massive DNA re-sequencing and gene expression studies have to be performed to substantially improve our understanding of crop genetics. Application of the knowledge obtained from the genomes, transcriptomes, expression studies, and epigenetic studies would enable the development of improved varieties and may lead to a second green revolution. The applications of next generation DNA sequencing technologies in crop improvement, its limitations, future prospects, and the features of important crop genome projects are reviewed herein.

Comparative Transcriptome Analysis of Latex Reveals Molecular Mechanisms Underlying Increased Rubber Yield in Hevea brasiliensis Self-Rooting Juvenile Clones

Rubber tree (Hevea brasiliensis) self-rooting juvenile clones (JCs) are promising planting materials for rubber production. In a comparative trial between self-rooting JCs and donor clones (DCs), self-rooting JCs exhibited better performance in rubber yield. To study the molecular mechanism associated with higher rubber yield in self-rooting JCs, we sequenced and comparatively analyzed the latex of rubber tree self-rooting JCs and DCs at the transcriptome level. Total raw reads of 34,632,012 and 35,913,020 bp were obtained from the library of self-rooting JCs and DCs, respectively, by using Illumina HiSeq 2000 sequencing technology. De novo assemblies yielded 54689 unigenes from the library of self-rooting JCs and DCs. Among 54689 genes, 1716 genes were identified as differentially expressed between self-rooting JCs and DCs via comparative transcript profiling. Functional analysis showed that the genes related to the mass of categories were differentially enriched between the two clones. Several genes involved in carbohydrate metabolism, hormone metabolism and reactive oxygen species scavenging were up-regulated in self-rooting JCs, suggesting that the self-rooting JCs provide sufficient molecular basis for the increased rubber yielding, especially in the aspects of improved latex metabolisms and latex flow. Some genes encoding epigenetic modification enzymes were also differentially expressed between self-rooting JCs and DCs. Epigenetic modifications may lead to gene differential expression between self-rooting JCs and DCs. These data will provide new cues to understand the molecular mechanism underlying the improved rubber yield of H. brasiliensis self-rooting clones.

HbMADS4, a MADS-box Transcription Factor from Hevea brasiliensis, Negatively Regulates HbSRPP

In plants MADS-box transcription factors (TFs) play important roles in growth and development. However, no plant MADS-box gene has been identified to have a function related to secondary metabolites regulation. Here, a MADS-box TF gene, designated as HbMADS4, was isolated from Hevea brasiliensis by the yeast one-hybrid experiment to screen the latex cDNA library using the promoter of the gene encoding H. brasiliensis small rubber particle protein (HbSRPP) as bait. HbMADS4 was 984-bp containing 633-bp open reading frame encoding a deduced protein of 230 amino acid residues with a typical conserved MADS-box motif at the N terminus. HbMADS4 was preferentially expressed in the latex, but little expression was detected in the leaves, flowers, and roots. Its expression was inducible by methyl jasmonate and ethylene. Furthermore, transient over-expression and over-expression of HbMADS4 in transgenic tobacco plants significantly suppressed the activity of the HbSRP promoter. Altogether, it is proposed that HbMADS4 is a negative regulator of HbSRPP which participates in the biosynthesis of natural rubber.

Comprehensive Proteomics Analysis of Laticifer Latex Reveals New Insights into Ethylene Stimulation of Natural Rubber Production

Ethylene is a stimulant to increase natural rubber latex. After ethylene application, both fresh yield and dry matter of latex are substantially improved. Moreover, we found that ethylene improves the generation of small rubber particles. However, most genes involved in rubber biosynthesis are inhibited by exogenous ethylene. Therefore, we conducted a proteomics analysis of ethylene-stimulated rubber latex, and identified 287 abundant proteins as well as 143 ethylene responsive latex proteins (ERLPs) with mass spectrometry from the 2-DE and DIGE gels, respectively. In addition, more than 1,600 proteins, including 404 ERLPs, were identified by iTRAQ. Functional classification of ERLPs revealed that enzymes involved in post-translational modification, carbohydrate metabolism, hydrolase activity, and kinase activity were overrepresented. Some enzymes for rubber particle aggregation were inhibited to prolong latex flow, and thus finally improved latex production. Phosphoproteomics analysis identified 59 differential phosphoproteins; notably, specific isoforms of rubber elongation factor and small rubber particle protein that were phosphorylated mainly at serine residues. This post-translational modification and isoform-specific phosphorylation might be important for ethylene-stimulated latex production. These results not only deepen our understanding of the rubber latex proteome but also provide new insights into the use of ethylene to stimulate rubber latex production.

Whole-transcriptome survey of the putative ATP-binding cassette (ABC) transporter family genes in the latex-producing laticifers of Hevea brasiliensis

The ATP-binding cassette (ABC) proteins or transporters constitute a large protein family in plants and are involved in many different cellular functions and processes, including solute transportation, channel regulation and molecular switches, etc. Through transcriptome sequencing, a transcriptome-wide survey and expression analysis of the ABC protein genes were carried out using the laticiferous latex from Hevea brasiliensis (rubber tree). A total of 46 putative ABC family proteins were identified in the H. brasiliensis latex. These consisted of 12 ‘full-size’, 21 ‘half-size’ and 13 other putative ABC proteins, and all of them showed strong conservation with their Arabidopsis thaliana counterparts. This study indicated that all eight plant ABC protein paralog subfamilies were identified in the H. brasiliensis latex, of which ABCB, ABCG and ABCI were the most abundant. Real-time quantitative reverse transcription-polymerase chain reaction assays demonstrated that gene expression of several latex ABC proteins was regulated by ethylene, jasmonic acid or bark tapping (a wound stress) stimulation, and that HbABCB15, HbABCB19, HbABCD1 and HbABCG21 responded most significantly of all to the abiotic stresses. The identification and expression analysis of the latex ABC family proteins could facilitate further investigation into their physiological involvement in latex metabolism and rubber biosynthesis by H. brasiliensis.

Transcriptome sequencing and comparative analysis reveal long-term flowing mechanisms in Hevea brasiliensis latex

BACKGROUND

The rubber tree, Hevea brasiliensis, is a major commercial source of natural rubber. Increasing the rubber yield of rubber trees is a very serious problem since the demands for high quality rubber materials are great. Establishment of a tapping system is based on an estimate of tapping intensity from the rubber tree. Latex flowing time is one of the most critical factors that determine the rubber yield. Long-term flow is a type of phenomenon of the rubber tree latex with longer flowing time than normal latex flow, and is always caused by intensive tapping. Thus, transcriptome and expression profiling data for long-term flowing latex (LFL) are needed as an important resource to identify genes and to better understand the biological mechanisms of latex flow in rubber trees.

RESULTS

The transcripts were sequenced using the Illumina sequencing platform. After cleaning, quality checks and sequencing, 98,697 transcripts and 38,584 unigenes were assembled with the mean size of 1437.31bp and 923.86bp, respectively. In BLAST searches of our database against public databases, 65.17% (25,147) of the unigenes were annotated with gene descriptions, conserved protein domains, or gene ontology terms. Functional categorization further revealed 853 individual unigenes related to long-term flow. According to KEGG classification, the clusters for "cysteine and methionine metabolism", "energy", "oxidative phosphorylation", "terpenoid backbone biosynthesis", "plant hormone signal transduction" and "copper, potassium transporter" were significantly enriched metabolic pathways.

CONCLUSIONS

We conducted high-resolution transcriptome profiling related to LFL in H. brasiliensis. The research facilitates further studies on gene discovery and on the molecular mechanisms related to the estimation of tapping intensity and prolonging latex flowing time. We concluded that it was necessary to improve energy supplies for intensive tapping and the copper ion content of rubber tree latex could be considered as a standard to estimate tapping intensity.

Identification and characterization of the 14-3-3 gene family in Hevea brasiliensis

The 14-3-3 proteins are a family of conserved phospho-specific binding proteins involved in diverse physiological processes. Although the genome-wide analysis of this family has been carried out in certain plant species, little is known about 14-3-3 protein genes in rubber tree (Hevea brasiliensis). In this study, we identified 10 14-3-3 protein genes (designated as HbGF14a to HbGF14j) in the latest rubber tree genome. A phylogenetic tree was constructed and found to demonstrate that HbGF14s can be divided into two major groups. Tissue-specific expression profiles showed that 10 HbGF14 were expressed in at least one of the tissues, which suggested that HbGF14s participated in numerous cellular processes. The 10 HbGF14s responded to jasmonic acid (JA) and ethylene (ET) treatment, which suggested that these HbGF14s were involved in response to JA and ET signaling. The target of HbGF14c protein was related to small rubber particle protein, a major rubber particle protein that is involved in rubber biosynthesis. These findings suggested that 14-3-3 proteins may be involved in the regulation of natural rubber biosynthesis.