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

Construction of the first high-density SNP genetic map and identification of QTLs for the natural rubber content in Taraxacum kok-saghyz Rodin

BACKGROUND

Taraxacum kok-saghyz Rodin (TKS) is a promising commercial alternative natural rubber (NR) yielding plant. Cultivating TKS with a high NR content is an important breeding target, and developing molecular markers related to NR content can effectively accelerate the breeding process of TKS.

RESULTS

To construct a high-density SNP genetic map and uncover genomic regions related to the NR content in TKS, an F₁ mapping population of TKS was constructed by crossing two parents (l66 and X51) with significant differences in NR contents. The NR content of the F₁ plants ranged from 0.30 to 15.14% and was distributed normally with a coefficient of variation of 47.61%, indicating quantitative trait inheritance. Then, employing whole-genome resequencing (WGR), a TKS genetic linkage map of 12,680 bin markers comprising 322,439 SNPs was generated. Based on the genetic map and NR content of the F₁ population, six quantitative trait loci (QTLs) for NR content with LOD > 4.0 were identified on LG01/Chr01 and LG06/Chr06. Of them, the 2.17 Mb genomic region between qHRC-C6-1 and qHRC-C6-2 on ChrA06, with 65.62% PVE in total, was the major QTL region. In addition, the six QTLs have significant additive genetic effects on NR content and could be used to develop markers for marker-assisted selection (MAS) in TKS with a high NR content.

CONCLUSION

This work constructed the first high-density TKS genetic map and identified the QTLs and genomic regions controlling the NR content, which provides useful information for fine mapping, map-based cloning, and MAS in TKS.

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

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

Identification of cis conformation natural rubber and proteins in Ficus altissima Blume latex

Ficus altissima Blume, a horticultural plant in tropical and warm subtropical regions, can produce natural rubber with comparable molecular weight to the Hevea brasiliensis rubber. The F. altissima latex has an acidic pH (about 4.89). The rubber particle size distribution is a unimodal profile, and the peak frequency is at a size of 4.5 μm. The natural rubber of F. altissima was determined to be a cis conformation via ¹³C NMR. The Mp (molecular weight of the peak maxima) of the deproteinized F. altissima rubber was 9.34 × 10⁵ Da. LC-MS was used to identify the proteins of rubber particles and serum. The most abundant protein of the creamy rubber particle layer is an acid phosphatase, while the most abundant proteins of serum were an (R)-mandelonitrilelyase and a polygalacturonase inhibitor. Pharmaceutical proteins (ficins) or enzymes related to the biosynthesis of natural medicines (a cannabidiolic acid synthase and two lupeol synthase) were identified in F. altissima latex. The data of this study may be helpful for research on the functions of latex in latex-borne plants and the biosynthesis mechanism of natural rubber.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Hevea brasiliensis latex proteomics: a review of analytical methods and the way forward

Natural rubber or latex from the Hevea brasiliensis is an important commodity in various economic sectors in today's modern society. Proteins have been detected in latex since the early twentieth century, and they are known to regulate various biological pathways within the H. brasiliensis trees such as the natural rubber biosynthesis, defence against pathogens, wound healing, and stress tolerance. However, the exact mechanisms of the pathways are still not clear. Proteomic analyses on latex have found various proteins and revealed how they fit into the mechanisms of the biological pathways. In the past three decades, there has been rapid latex protein identification due to the improvement of latex protein extraction methods, as well as the emergence of two-dimensional gel electrophoresis (2-DE) and mass spectrometry (MS). In this manuscript, we reviewed the methods of latex protein extraction that keeps on improving over the past three decades as well as the results of numerous latex protein identification and quantitation.

Comparison of Morphological Characteristics and Determination of Different Patterns for Rubber Particles in Dandelion and Different Rubber Grass Varieties

Russian dandelion Taraxacum kok-saghyz (TKS) is one promising alternative crop for natural rubber production. However, it is easily confused with other dandelions. In this study, we performed a systematical comparison of the morphological characteristics for different TKS varieties and common dandelion Taraxacum officinale (TO). Our results demonstrated that several obvious differences in morphology can be found between TKS and TO. TO leaf is a pinnate shape, its margin is heavily jagged and its base is cuneate, but TKS leaf is more cuneate and its leaf margin is nearly smooth and round. There are obvious differences for the outer bracts of TO and TKS flower buds. TKS bracts are oblanceolate, apex obtuse, margin smooth and sinuate, and its outer layer of flower buds and faceplate involucre sepal is buckled inward to form a certain angle. TKS is self-incompatible, and its seeds are spindle-shaped achene and show upright plumpness. A large amount of laticifer cells and rubber particles can be detected from many TKS tissues, and dry roots of TKS contain high contents of natural rubber. Laticifer cells and rubber particles can only be examined in the vein, stem, and roots of TKS. Our statical results also revealed that the numbers of laticifer cells and rubber particles have a positive relationship with the rubber content in TKS roots. These morphological features can help us to easily distinguish TKS from common dandelion and approximately estimate the rubber content in the roots of different TKS varieties for TKS breeding in future.

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.

Alteration of metabolic profiles in Lemna paucicostata culture and enhanced production of GABA and ferulic acid by ethephon treatment

Lemna species have been used in the food, feed, and pharmaceutical industries, as they are inexpensive sources of proteins, starches, and fatty acids. In this study, we treated L. paucicostata with different concentrations (0.05, 0.1, 0.2, 0.5, or 1 mM) of ethephon. The total dry weight decreased in all ethephon-treated groups compared to the control group. We also investigated the alteration of metabolic profiles induced by ethephon treatment by using gas chromatography-mass spectrometry. This analysis identified 48 metabolites, and the relative levels of most of alcohols, amino acids, fatty acids, and phenols increased by the ethephon treatment, whereas levels of organic acids and sugars decreased. Among these, the highest production of γ-aminobutyric acid (GABA, 5.041 ± 1.373 mg/L) and ferulic acid (0.640 ± 0.071 mg/L) was observed in the 0.5 mM and the 0.2 mM ethephon treatment groups, respectively. These results could be useful for large-scale culture of L. paucicostata with enhanced GABA and ferulic acid content for utilization in the food, feed, cosmetic, and pharmaceutical industries.

Complementary transcriptome and proteome profiling in the mature seeds of Camellia oleifera from Hainan Island

Camellia oleifera Abel. (C. oleifera), as an important woody tree species producing edible oils in China, has attracted enormous attention due to its abundant unsaturated fatty acids and their associated benefits to human health. To reveal novel insights into the characters during the maturation period of this plant as well as the molecular basis of fatty acid biosynthesis and degradation, we conducted a conjoint analysis of the transcriptome and proteome of C. oleifera seeds from Hainan Island. Using RNA sequencing (RNA-seq) technology and shotgun proteomic method, 59,391 transcripts and 40,500 unigenes were obtained by TIGR Gene Indices Clustering Tools (TGICL), while 1691 protein species were identified from Mass Spectrometry (MS). Subsequently, all genes and proteins were employed in euKaryotic Orthologous Groups (KOG) classification, Gene Ontology (GO) annotation, and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis to investigate their essential functions. The results indicated that the most abundant pathways were biological metabolic processes. There were 946 unigenes associated with lipid metabolism at the transcriptome level, with 116 proteins at the proteome level; among these, 38 specific proteins were involved in protein-protein interactions, with the majority being related to fatty acid catabolic process. The expression levels of 21 candidate unigenes encoding target proteins were further detected by quantitative real-time polymerase chain reaction (qRT-PCR). Finally, Gas Chromatography Mass Spectrometry (GC-MS) was carried out to determine the fatty acid composition of C. oleifera oil. These findings not only deepened our understanding about the molecular mechanisms of fatty acid metabolism but also offered new evidence concerning the roles of relevant proteins in oil-bearing crops. Furthermore, the lipid-associated proteins recognized in this research might be helpful in providing a reference for the synthetic regulation of C. oleifera oil quality by genetic engineering techniques, thus resulting in potential application in agriculture.

Genome-wide identification and expression analysis of the phosphatase 2A family in rubber tree (Hevea brasiliensis)

The protein phosphatase 2As (PP2As) play a key role in manipulating protein phosphorylation. Although a number of proteins in the latex of laticifers are phosphorylated during latex regeneration in rubber tree, information about the PP2A family is limited. In the present study, 36 members of the HbPP2A family were genome-wide identified. They were clustered into five subgroups: the subgroup HbPP2AA (4), HbPP2AB' (14), HbPP2AB'' (6), HbPP2AB55 (4), and HbPP2AC (8). The members within the same subgroup shared highly conserved gene structures and protein motifs. Most of HbPP2As possessed ethylene- and wounding-responsive cis-acting elements. The transcripts of 29 genes could be detected in latex by using published high-throughput sequencing data. Of the 29 genes, seventeen genes were significantly down-regulated while HbPP2AA1-1 and HbPP2AB55α/Bα-1were up-regulated by tapping. Of the 17 genes, 14 genes were further significantly down-regulated by ethrel application. The down-regulated expression of a large number of HbPP2As may attribute to the enhanced phosphorylation of the proteins in latex from the tapped trees and the trees treated with ethrel application.

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

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

Proteomic Landscape 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.

Comparative Proteomic Analysis of the Hevea brasiliensis Latex under Ethylene and Calcium Stimulation

BACKGROUND

Calcium ions usually act as a second messenger in the signal transmission process and a major element required by plants. In Hevea, calcium ion could alleviate the negative effects of long-term ethylene application to a certain extent. However, the molecular mechanisms remain unclear.

METHODS

Two-dimensional electrophoresis was used to determine the pattern of protein changes in latex after treatments with calcium and/or ethylene. Quantitative real-time polymerase chain reaction and Western blotting were used to determine the expression levels of some proteins and genes. STRING software was used to determine the protein-protein interaction network of the identified proteins.

RESULTS

Comparative proteomics identified 145 differentially expressed proteins, which represented 103 unique proteins. The abundance change patterns of some proteins involved in signal transduction, rubber particle aggregation, and natural rubber biosynthesis were altered upon calcium stimulation. Quantitative real-time polymerase chain reaction analysis of 29 proteins showed that gene expression did not always maintain the same trend as protein expression. The increased enzyme activities of superoxide dismutase, ascorbate peroxidase, and glutathione reductase suggested that calcium can enhance the antistress ability of plants by increasing the activity of their antioxidant enzyme systems.

CONCLUSION

These results supplement the rubber latex proteome, and provide evidence for investigating the molecular mechanisms by which calcium alleviates the negative effects of ethylene stimulation.

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.

An improved protein extraction method applied to cotton leaves is compatible with 2-DE and LC-MS

Background

Two-dimensional electrophoresis (2-DE) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) are widely used in plant proteomics research. However, these two techniques cannot be simultaneously satisfied by traditional protein extraction methods when investigate cotton leaf proteome.

Results

Here, we evaluated the efficiency of three different protein extraction methods for 2-DE and LC-MS/MS analyses of total proteins obtained from cotton leaves. The protein yield of the borax/PVPP/phenol (BPP) method (0.14%) was significantly lower than the yields of the trichloroacetic acid/acetone (TCA) precipitation method (1.42%) and optimized TCA combined with BPP (TCA-B) method (0.47%). The BPP method was failed to get a clear 2-DE electrophoretogram. Fifty pairs of protein spots were randomly selected from the 2-DE gels of TCA- and TCA-B-extracted proteins for identification by MALDI TOF/TOF, and the results of 42 pairs were consistent. High-throughput proteomic analysis showed that 6339, 9282 and 9697 unique proteins were identified from the total cotton leaf proteins extracted by the TCA, BPP and TCA-B methods, respectively. Gene Ontology (GO) analysis revealed that the proteins specifically identified by TCA method were primarily distributed in the plasma membrane, while BPP and TCA-B methods specific proteins distributed in the cytosol, indicating the sub-cellular preference of different protein extraction methods. Further, ATP-dependent zinc metalloprotease FTSH 8 could be observed in the 2-DE gels of TCA and TCA-B methods, and could only be detected in the LC-MS/MS results of the BPP and TCA-B methods, showing that TCA-B method might be the optimized choice for both 2-DE and LC-MS/MS.

Conclusion

Our data provided an improved TCA-B method for protein extraction that is compatible with 2-DE and LC-MS/MS for cotton leaves and similar plant tissues which is rich in polysaccharides and polyphenols.

Electronic supplementary material

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

Biosynthesis of Natural Rubber: Current State and Perspectives

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

Transcriptome analysis in Hevea brasiliensis latex revealed changes in hormone signalling pathways during ethephon stimulation and consequent Tapping Panel Dryness

Tapping Panel Dryness (TPD) affects latex production in Hevea brasiliensis. This physiological syndrome involves the agglutination of rubber particles, which leads to partial or complete cessation of latex flow. Latex harvesting consists in tapping soft bark. Ethephon can be applied to stimulate latex flow and its regeneration in laticifers. Several studies have reported transcriptome changes in bark tissues. This study is the first report on deep RNA sequencing of latex to compare the effect of ethephon stimulation and TPD severity. Trees were carefully selected for paired-end sequencing using an Illumina HiSeq 2000. In all, 43 to 60 million reads were sequenced for each treatment in three biological replicates (slight TPD trees without ethephon stimulation, and slight and severe TPD trees with ethephon treatment). Differentially expressed genes were identified and annotated, giving 8,111 and 728 in response to ethephon in slight TPD trees and in ethephon-induced severe TPD trees, respectively. A biological network of responses to ethephon and TPD highlighted the major influence of metabolic processes and the response to stimulus, especially wounding and jasmonate depression in TPD-affected trees induced by ethephon stimulation.

Physiological and Proteomic Analyses of Molybdenum- and Ethylene-Responsive Mechanisms in Rubber Latex

Molybdenum (Mo) is an essential micronutrient in many plants. In the rubber tree Hevea brasiliensis, Mo application can reduce the shrinkage of the tapping line, decrease tapping panel dryness, and finally increase rubber latex yield. After combined Mo with ethylene (Eth), these effects become more obvious. However, the molecular mechanism remains unclear. Here, we compared the changed patterns of physiological parameters and protein accumulation in rubber latex after treated with Mo and/or Eth. Our results demonstrated that both Eth and Mo can improve the contents of thiol, sucrose, and dry yield in rubber latex. However, lutoid bursting is significantly inhibited by Mo. Comparative proteomics identified 169 differentially expressed proteins, including 114 unique proteins, which are mainly involved in posttranslational modification, carbohydrate metabolism, and energy production. The abundances of several proteins involved in rubber particle aggregation are decreased upon Mo stimulation, while many enzymes related to natural rubber biosynthesis are increased. Comparison of the accumulation patterns of 25 proteins revealed that a large portion of proteins have different changed patterns with their gene expression levels. Activity assays of six enzymes revealed that Mo stimulation can increase latex yield by improving the activity of some Mo-responsive enzymes. These results not only deepen our understanding of the rubber latex proteome but also provide new insights into the molecular mechanism of Mo-stimulated rubber latex yield.

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.

Transcriptome analysis of Hevea brasiliensis in response to exogenous methyl jasmonate provides novel insights into regulation of jasmonate-elicited rubber biosynthesis

The phytohomorne methyl jasmonate (MeJA) is known to trigger extensive reprogramming of gene expression leading to transcriptional activation of many secondary metabolic pathways. However, natural rubber is a commercially important secondary metabolite and little is known about the genetic and genomic basis of jasmonate-elicited rubber biosynthesis in rubber tree (Hevea brasiliensis). RNA sequencing (RNA-seq) of H. brasiliensis bark treated with 1 g lanolin paste containing 0.02% w/w MeJA for 24 h (M2) and 0.04% w/w MeJA for 24 h (M4) was performed. A total of 2950 and 2850 differentially expressed genes in M2 and M4 compared with control (C) were respectively detected. Key genes involved in 2-C-methyl-D-erythritol 4-phosphate, rubber biosynthesis, glycolysis and carbon fixation (Calvin cycle) pathway were found to be up-regulated by MeJA treatment. Particularly, the expression of 3-hydroxy-3-metylglutaryl coenzyme A reductase in MVA pathway was down-regulated by MeJA treatment, but the expression of farnesyl diphosphate synthase (FPS) and cis-prenyltransferase (CPT, or rubber transferase) in rubber biosynthesis pathway were up-regulated by MeJA treatment. Up-regulation of critical genes in JA biosynthesis in response to MeJA treatment exhibited the self-activation of JA biosynthesis. In addition, up-regulated genes of great regulatory importance in cross-talk between JA and other hormone signaling, and of transcriptional regulation were identified. The increased expression levels of FPS and CPT in rubber biosynthesis pathway possibly resulted in an increased latex production in rubber tree treated with MeJA. The present results provide insights into the mechanism by which MeJA activates the rubber biosynthesis and the transcriptome data can also serve as the foundation for future research into the molecular basis for MeJA regulation of other cellular processes.

A Multi-Omics Analysis of Glycine max Leaves Reveals Alteration in Flavonoid and Isoflavonoid Metabolism Upon Ethylene and Abscisic Acid Treatment

Phytohormones are central to plant growth and development. Despite the advancement in our knowledge of hormone signaling, downstream targets, and their interactions upon hormones action remain largely fragmented, especially at the protein and metabolite levels. With an aim to get new insight into the effects of two hormones, ethylene (ET) and abscisic acid (ABA), this study utilizes an integrated proteomics and metabolomics approach to investigate their individual and combined (ABA+ET) signaling in soybean leaves. Targeting low-abundance proteins, our previously established protamine sulfate precipitation method was applied, followed by label-free quantification of identified proteins. A total of 4129 unique protein groups including 1083 differentially modulated in one (individual) or other (combined) treatments were discerned. Functional annotation of the identified proteins showed an increased abundance of proteins related to the flavonoid and isoflavonoid biosynthesis and MAPK signaling pathway in response to ET treatment. HPLC analysis showed an accumulation of isoflavones (genistin, daidzein, and genistein) upon ET treatment, in agreement with the proteomics results. A metabolome analysis assigned 79 metabolites and further confirmed the accumulation of flavonoids and isoflavonoids in response to ET. A potential cross-talk between ET and MAPK signaling, leading to the accumulation of flavonoids and isoflavonoids in soybean leaves is suggested.

Calcium-Dependent Protein Kinase Family Genes Involved in Ethylene-Induced Natural Rubber Production in Different Hevea brasiliensis Cultivars

Natural rubber latex production can be improved by ethylene stimulation in the rubber tree (Hevea brasiliensis). However, the expression levels of most functional proteins for natural rubber biosynthesis are not induced after ethylene application, indicating that post-translational modifications, especially protein phosphorylation, may play important roles in ethylene signaling in Hevea. Here, we performed a comprehensive investigation on evolution, ethylene-induced expression and protein-protein interaction of calcium-dependent protein kinases (CPKs), an important serine/threonine protein kinase family, in Hevea. Nine duplication events were determined in the 30 identified HbCPK genes. Expression profiling of HbCPKs in three rubber tree cultivars with low, medium and high ethylene sensitivity showed that HbCPK6, 17, 20, 22, 24, 28 and 30 are induced by ethylene in at least one cultivar. Evolution rate analysis suggested accelerated evolution rates in two paralogue pairs, HbCPK9/18 and HbCPK19/20. Analysis of proteomic data for rubber latex after ethylene treatment showed that seven HbCPK proteins could be detected, including six ethylene-induced ones. Protein-protein interaction analysis of the 493 different abundant proteins revealed that protein kinases, especially calcium-dependent protein kinases, possess most key nodes of the interaction network, indicating that protein kinase and protein phosphorylation play important roles in ethylene signaling in latex of Hevea. In summary, our data revealed the expression patterns of HbCPK family members and functional divergence of two HbCPK paralogue pairs, as well as the potential important roles of HbCPKs in ethylene-induced rubber production improvement in Hevea.

Expression Profiling of Mitogen-Activated Protein Kinase Genes Reveals Their Evolutionary and Functional Diversity in Different Rubber Tree (Hevea brasiliensis) Cultivars

Rubber tree (Heveabrasiliensis) is the only commercially cultivated plant for producing natural rubber, one of the most essential industrial raw materials. Knowledge of the evolutionary and functional characteristics of kinases in H. brasiliensis is limited because of the long growth period and lack of well annotated genome information. Here, we reported mitogen-activated protein kinases in H. brasiliensis (HbMPKs) by manually checking and correcting the rubber tree genome. Of the 20 identified HbMPKs, four members were validated by proteomic data. Protein motif and phylogenetic analyses classified these members into four known groups comprising Thr-Glu-Tyr (TEY) and Thr-Asp-Tyr (TDY) domains, respectively. Evolutionary and syntenic analyses suggested four duplication events: HbMPK3/HbMPK6, HbMPK8/HbMPK9/HbMPK15, HbMPK10/HbMPK12 and HbMPK11/HbMPK16/HbMPK19. Expression profiling of the identified HbMPKs in roots, stems, leaves and latex obtained from three cultivars with different latex yield ability revealed tissue- and variety-expression specificity of HbMPK paralogues. Gene expression patterns under osmotic, oxidative, salt and cold stresses, combined with cis-element distribution analyses, indicated different regulation patterns of HbMPK paralogues. Further, Ka/Ks and Tajima analyses suggested an accelerated evolutionary rate in paralogues HbMPK10/12. These results revealed HbMPKs have diverse functions in natural rubber biosynthesis, and highlighted the potential possibility of using MPKs to improve stress tolerance in future rubber tree breeding.

Comparative Proteomics of Rubber Latex Revealed Multiple Protein Species of REF/SRPP Family Respond Diversely to Ethylene Stimulation among Different Rubber Tree Clones

Rubber elongation factor (REF) and small rubber particle protein (SRPP) are two key factors for natural rubber biosynthesis. To further understand the roles of these proteins in rubber formation, six different genes for latex abundant REF or SRPP proteins, including REF_138,175,258 and SRPP_117,204,243, were characterized from Hevea brasiliensis Reyan (RY) 7-33-97. Sequence analysis showed that REFs have a variable and long N-terminal, whereas SRPPs have a variable and long C-terminal beyond the REF domain, and REF₂₅₈ has a β subunit of ATPase in its N-terminal. Through two-dimensional electrophoresis (2-DE), each REF/SRPP protein was separated into multiple protein spots on 2-DE gels, indicating they have multiple protein species. The abundance of REF/SRPP proteins was compared between ethylene and control treatments or among rubber tree clones with different levels of latex productivity by analyzing 2-DE gels. The total abundance of each REF/SRPP protein decreased or changed a little upon ethylene stimulation, whereas the abundance of multiple protein species of the same REF/SRPP changed diversely. Among the three rubber tree clones, the abundance of the protein species also differed significantly. Especially, two protein species of REF₁₇₅ or REF₂₅₈ were ethylene-responsive only in the high latex productivity clone RY 8-79 instead of in RY 7-33-97 and PR 107. Some individual protein species were positively related to ethylene stimulation and latex productivity. These results suggested that the specific protein species could be more important than others for rubber production and post-translational modifications might play important roles in rubber biosynthesis.

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 and functional characterization of HbOsmotin from Hevea brasiliensis

Latex in the laticiferous cell network of Hevea brasiliensis tree is composed of cytoplasm that synthesizes natural rubber. Ethylene stimulation of the tree bark enhances latex production partly by prolonging the duration of latex flow during the tapping process. Here, we identified an osmotin-like cDNA sequence (HbOsmotin) from H. brasiliensis that belongs to the pathogenesis-related 5 (PR-5) gene family. The HbOsmotin protein is present in the lutoids of latex in H. brasiliensis, whereas in onion epidermal cells, this protein is predominantly distributed around the cell wall, suggesting that it may be secreted from the cytoplasm. We investigated the effects of exogenous ethylene on HbOsmotin transcription and protein accumulation in rubber latex, and further determined the protein function after osmotic stress in Arabidopsis. In regularly tapped trees, HbOsmotin expression was drastically inhibited in rubber latex after tapping, although the expression was subsequently recovered by ethylene stimulation. However, in virgin plants that had never been tapped, exogenous ethylene application slightly decreased HbOsmotin expression. HbOsmotin overexpression in Arabidopsis showed that HbOsmotin reduced the osmotic stress tolerance of the plant, which likely occurred by raising the water potential. These data indicated that HbOsmotin may contribute to osmotic regulation in laticiferous cells.

A protein extraction method for low protein concentration solutions compatible with the proteomic analysis of rubber particles

The extraction of high-purity proteins from the washing solution (WS) of rubber particles (also termed latex-producing organelles) from laticifer cells in rubber tree for proteomic analysis is challenging due to the low concentration of proteins in the WS. Recent studies have revealed that proteins in the WS might play crucial roles in natural rubber biosynthesis. To further examine the involvement of these proteins in natural rubber biosynthesis, we designed an efficiency method to extract high-purity WS proteins. We improved our current borax and phenol-based method by adding reextraction steps with phenol (REP) to improve the yield from low protein concentration samples. With this new method, we extracted WS proteins that were suitable for proteomics. Indeed, compared to the original borax and phenol-based method, the REP method improved both the quality and quantity of isolated proteins. By repeatedly extracting from low protein concentration solutions using the same small amount of phenol, the REP method yielded enough protein of sufficiently high-quality from starting samples containing less than 0.02 mg of proteins per milliliter. This method was successfully applied to extract the rubber particle proteins from the WS of natural rubber latex samples. The REP-extracted WS proteins were resolved by 2DE, and 28 proteins were positively identified by MS. This method has the potential to become widely used for the extraction of proteins from low protein concentration solutions for proteomic analysis.

Evaluation of Reference Genes for Quantitative Real-Time PCR Analysis of the Gene Expression in Laticifers on the Basis of Latex Flow in Rubber Tree (Hevea brasiliensis Muell. Arg.)

Latex exploitation-caused latex flow is effective in enhancing latex regeneration in laticifer cells of rubber tree. It should be suitable for screening appropriate reference gene for analysis of the expression of latex regeneration-related genes by quantitative real-time PCR (qRT-PCR). In the present study, the expression stability of 23 candidate reference genes was evaluated on the basis of latex flow by using geNorm and NormFinder algorithms. Ubiquitin-protein ligase 2a (UBC2a) and ubiquitin-protein ligase 2b (UBC2b) were the two most stable genes among the selected candidate references in rubber tree clones with differential duration of latex flow. The two genes were also high-ranked in previous reference gene screening across different tissues and experimental conditions. By contrast, the transcripts of latex regeneration-related genes fluctuated significantly during latex flow. The results suggest that screening reference gene during latex flow should be an efficient and effective clue for selection of reference genes in qRT-PCR.

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

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