Integrated proteomic and phosphoproteomic analysis reveals the MAPK cascade as a key regulator of ethylene-induced latex production in Hevea brasiliensis

Abstract

Natural rubber is a crucial industrial raw material globally, with wide applications in industries such as automotive and healthcare. Application of ethylene promotes latex production in the Hevea rubber tree (Hevea brasiliensis). However, its potential protein-regulation mechanism remains unclear. To elucidate the role of ethylene in latex production, we employed high-resolution proteomics and phosphoproteomics to determine the global changes in protein abundance and phosphorylation during ethylene-stimulated latex production. Using latex samples collected from ethephon-treated rubber tree, we identified more than 3,700 quantifiable proteins and 2,000 phosphorylated proteins with over 6,000 phosphorylation sites. Proteins involved in the mitogen-activated protein kinase (MAPK) cascade, metabolic pathways, and vesicular trafficking, like rubber elongation factors (REFs), small rubber particle proteins (SRPPs), and 14–3-3 proteins, exhibit pronounced robust phosphorylation dynamics upon ethephon application. Interestingly, the serine/threonine–proline (S/T-P) motif, which could be recognized by MAPKs, was highly enriched in REFs, 14–3-3, and ubiquitin-associated DENN (UDENN) domain-containing proteins, suggesting a potential role of MAPKs in ethylene-induced rubber biosynthesis. Consistently, the MAPKs show altered activation in ethephon-treated Hevea rubber trees, compared to water-treated controls. Additionally, activation of UDENN domain proteins implicated Ras-related Rab GTPases in membrane trafficking and rubber particle formation. These results provide comprehensive information on global changes in protein abundance and phosphorylation upon ethephon application in latex production, and may provide valuable clues for understanding the molecular basis of ethephon-regulated natural rubber biosynthesis in Hevea rubber trees.

Keywords

Proteomic / Phosphoproteomic / Latex / Ethylene / Rubber tree

References

[1] Amalou Z, Bangratz J, Chrestin H. Ethrel (ethylene releaser)-induced increases in the adenylate pool and transtonoplast deltapH within hevea latex cells. Plant Physiol, 1992, 98: 1270-1276

[2] Amerik AY, Martirosyan YT, Gachok IV. Regulation of natural rubber biosynthesis by proteins associated with rubber particles. Russ J Bioorg Chem, 2018, 44: 140-149

[3] Audley BG, Archer BL, Carruthers IB. Metabolism of ethephon (2-chloroethylphosphonic acid) and related compounds in Hevea brasiliensis. Arch Environ Contam Toxicol, 1976, 4: 183-200

[4] Bakshi A, Shemansky JM, Chang C, Binder BM. History of research on the plant hormone ethylene. J Plant Growth Regul, 2015, 34: 809-827

[5] Benschop JJ, Mohammed S, O'Flaherty M, Heck AJ, Slijper M, Menke FL. Quantitative phosphoproteomics of early elicitor signaling in Arabidopsis. Mol Cell Proteomics, 2007, 6: 1198-1214

[6] Berthelot K, Lecomte S, Estevez Y, Coulary-Salin B, Peruch F. Homologous Hevea brasiliensis REF (Hevb1) and SRPP (Hevb3) present different auto-assembling. Biochim Biophys Acta, 2014, 1844: 473-485

[7] Buakong W, Suwanmanee P, Vasinayanuwatana T, Ruttajorn K, Assawatreeratanakul K, Leake JE. Characterization and expression analysis of HbC3H66: implications for transcriptional regulation in rubber biosynthesis and abiotic stress responses in Hevea brasiliensis. Trends Sci, 2023, 20 6835–6835

[8] Camoni L, Visconti S, Aducci P, Marra M. 14–3-3 proteins in plant hormone signaling: doing several things at once. Front Plant Sci, 2018, 9: 297

[9] Cantalapiedra CP, Hernández-Plaza A, Letunic I, Bork P, Huerta-Cepas J. eggNOG-mapper v2: Functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol Biol Evol, 2021, 38: 5825-5829

[10] Cao YH, Chen QW, Xu X, Fernie A, Li J, Zhang Y. Insights from natural rubber biosynthesis evolution for pathway engineering. Trends Plant Sci, 2025, 30: 897-912

[11] Chao JQ, Chen Y, Wu S, Tian WM. Comparative transcriptome analysis of latex from rubber tree clone CATAS8-79 and PR107 reveals new cues for the regulation of latex regeneration and duration of latex flow. BMC Plant Biol, 2015, 15: 104

[12] Chen HM, Xiong F, Wangler AM, Bischoff T, Wang K, Miao Yet al.. Phosphorylation-dependent activation of the bHLH transcription factor ICE1/SCRM promotes polarization of the Arabidopsis zygote. New Phytol, 2025, 245: 1029-1039

[13] Cheng A, Grant CE, Noble WS, Bailey TL. MoMo: Discovery of statistically significant post-translational modification motifs. Bioinformatics, 2019, 35: 2774-2782

[14] Cherian S, Ryu SB, Cornish K. Natural rubber biosynthesis in plants, the rubber transferase complex, and metabolic engineering progress and prospects. Plant Biotechnol J, 2019, 17: 2041-2061

[15] Chevalier D, Morris ER, Walker JC. 14–3-3 and FHA domains mediate phosphoprotein interactions. Annu Rev Plant Biol, 2009, 60: 67-91

[16] Chrispeels MJ, Herman EM. Endoplasmic reticulum-derived compartments function in storage and as mediators of vacuolar remodeling via a new type of organelle, precursor protease vesicles. Plant Physiol, 2000, 123: 1227-1234

[17] Dai LJ, Kang GJ, Nie ZY, Li Y, Zeng R. Comparative proteomic analysis of latex from Hevea brasiliensis treated with Ethrel and methyl jasmonate using iTRAQ-coupled two-dimensional LC-MS/MS. J Proteomics, 2016, 132: 167-175

[18] Coupé M, Chrestin H (1989) Physico-chemical and biochemical mechanisms of hormonal (ethylene) stimulation. In D‘Auzac J, Jacob JL, Chrestin H, Physiology of rubber tree latex: the laticiferous cell and latex, a model of cytoplasm. CRC Press, Boca Raton. pp 295–318. https://doi.org/10.1201/9781351075695

[19] Day EK, Sosale NG, Lazzara MJ. Cell signaling regulation by protein phosphorylation: a multivariate, heterogeneous, and context-dependent process. Curr Opin Biotechnol, 2016, 40: 185-192

[20] Denison F, Paul A, Zupanska A, Ferl R. 14–3-3 proteins in plant physiology. Semin Cell Dev Biol, 2011, 22: 720-727

[21] Ding H, Wu S, Xia Z, Huang X, Shi M. Aggregation-related proteins binding to rubber particles with different diameter in rubber tree. Chin J Trop Crops, 2022, 43: 2413-2421

[22] Duan CF, Rio M, Leclercq J, Bonnot F, Oliver G, Montoro P. Gene expression pattern in response to wounding, methyl jasmonate and ethylene in the bark of Hevea brasiliensis. Tree Physiol, 2010, 30: 1349-1359

[23] Duan CF, Nie Z, Zeng R (2006) Establishment of 2-DE system and primary analyses on the membrane proteins of rubber particles in Hevea brasiliensis by MALDI-TOF. Chin J Trop Crops 27:22–29. https://api.semanticscholar.org/CorpusID:87436544

[24] Fan Y, Qi J, Xiao X, Li H, Lan J, Tang CRet al.. Transcript and protein profiling provides insights into the molecular mechanisms of harvesting-induced latex production in rubber tree. Front Genet, 2022, 13: 756270

[25] Fang Y, Xiao X, Lin J, Lin Q, Wang J, Liu Ket al.. Pan-genome and phylogenomic analyses highlight Hevea species delineation and rubber trait evolution. Nat Commun, 2024, 15: 7232

[26] Fang FY, Yuan BX, He LX, He MM, Wang XC. Roles of rubber elongation factor and small rubber particle protein in rubber particles. Plant Mol Biol, 2025, 115: 71

[27] Farooq A, Zhou MM. Structure and regulation of MAPK phosphatases. Cell Signal, 2004, 16: 769-779

[28] Gong L, Liu X, Jiao L, Yang X, Lemoff A, Liu X. CK2-mediated phosphorylation of SUZ12 promotes PRC2 function by stabilizing enzyme active site. Nat Commun, 2022, 13: 6781

[29] Harker M, Holmberg N, Clayton JC, Gibbard CL, Wallace AD, Rawlins Set al.. Enhancement of seed phytosterol levels by expression of an N-terminal truncated Hevea brasiliensis (rubber tree) 3-hydroxy-3-methylglutaryl-CoA reductase. Plant Biotechnol J, 2003, 1: 113-121

[30] He H, Wang J, Meng Z, Dijkwel PP, Du P, Shi Set al.. Genome-wide analysis of the SRPP/REF gene family in Taraxacum kok-saghyz provides insights into its expression patterns in response to ethylene and methyl jasmonate treatments. Int J Mol Sci, 2024, 25: 6864

[31] Heintz D, Gallien S, Compagnon V, Berna A, Suzuki M, Yoshida Set al.. Phosphoproteome exploration reveals a reformatting of cellular processes in response to low sterol biosynthetic capacity in Arabidopsis. J Proteome Res, 2012, 11: 1228-1239

[32] Hillebrand A, Post JJ, Wurbs D, Wahler D, Lenders M, Krzyzanek Vet al.. Down-regulation of small rubber particle protein expression affects integrity of rubber particles and rubber content in Taraxacum brevicorniculatum. PLoS ONE, 2018, 7: e41874

[33] Huang XJ, Sun MM, Du XY, Quan L, Chao J, Deng X, Liao H (2025) Genome-wide identification of the regulatory network of mitogen-activated protein kinase signaling cascades gene families in Hevea brasiliensis. BMC Plant Biol 25:581. https://doi.org/10.1186/s12870-025-06615-6

[34] Jiang Z, Yao L, Zhu X, Hao G, Ding Y, Zhao Het al.. Ethylene signaling modulates air humidity responses in plants. Plant J, 2024, 117: 653-668

[35] Kitata R, Yang J, Chen Y. Advances in data-independent acquisition mass spectrometry towards comprehensive digital proteome landscape. Mass Spectrom Rev, 2023, 42: 2324-2348

[36] Laibach N, Schmidl S, Müller B, Bergmann M, Prüfer D, Schulze Gronover C. Small rubber particle proteins from Taraxacum brevicorniculatum promote stress tolerance and influence the size and distribution of lipid droplets and artificial poly(cis-1,4-isoprene) bodies. Plant J, 2018, 93: 1045-1061

[37] Lestari R, Rio M, Martin F, Leclercq J, Woraathasin N, Roques Set al.. Overexpression of Hevea brasiliensis ethylene response factor HbERF-IXc5 enhances growth and tolerance to abiotic stress and affects laticifer differentiation. Plant Biotechnol J, 2018, 16: 322-336

[38] Li JN, Silva-Sanchez C, Zhang T, Chen S, Li H. Phosphoproteomics technologies and applications in plant biology research. Front Plant Sci, 2015, 6: 430

[39] Li PP, Ma H, Xiao N, Zhang Y, Xu TY, Xia T. Overexpression of the ZmSUS1 gene alters the content and composition of endosperm starch in maize (Zea mays L.). Planta, 2023, 257: 97

[40] Lin L, Zheng J, Yu Q, Chen W, Xing J, Chen Cet al.. High throughput and accurate serum proteome profiling by integrated sample preparation technology and single-run data independent mass spectrometry analysis. J Proteomics, 2018, 174: 9-16

[41] Liu J. Molecular mechanism underlying ethylene stimulation of latex production in rubber tree (Hevea brasiliensis). Trees, 2016, 30: 1913-1921

[42] Liu JP, Zhuang Y, Guo X, Li Y. Molecular mechanism of ethylene stimulation of latex yield in rubber tree (Hevea brasiliensis) revealed by de novo sequencing and transcriptome analysis. BMC Genomics, 2016, 17: 257

[43] Liu J, Shi C, Shi CC, Li W, Zhang QJ, Zhang Yet al.. The chromosome-based rubber tree genome provides new insights into spurge genome evolution and rubber biosynthesis. Mol Plant, 2020, 13: 336-350

[44] Liu H, Yuan K, Hu Y, Wang S, He Q, Feng Cet al.. Construction and analysis of the tapping panel dryness-related lncRNA/circRNA-miRNA-mRNA ceRNA network in latex ofHevea brasiliensis. Plant Physiol Biochem, 2023, 205: 108156

[45] Marat AL, Dokainish H, McPherson PS. Denn domain proteins: regulators of Rab GTPases. J Biol Chem, 2011, 286: 13791-13800

[46] Mehrpouyan S, Menon U, Tetlow IJ, Emes MJ. Protein phosphorylation regulates maize endosperm starch synthase IIa activity and protein-protein interactions. Plant J, 2021, 105: 1098-1112

[47] Men X, Wang F, Chen GQ, Zhang HB, Xian M. Biosynthesis of natural rubber: current state and perspectives. Int J Mol Sci, 2018, 20: 50

[48] Müller-Dott S, Jaehnig EJ, Munchic KP, Jiang W, Yaron-Barir TM, Savage SRet al.. Comprehensive evaluation of phosphoproteomic-based kinase activity inference. Nat Commun, 2025, 16: 4771

[49] Nawamawat K, Sakdapipanich JT, Ho CC, Ma Y, Song J, Vancso JG. Surface nanostructure of Hevea brasiliensis natural rubber latex particles. Colloids Surf A Physicochem Eng Asp, 2011, 390: 157-166

[50] Okuda S, Yoshizawa AC, Kobayashi D, Takahashi Y, Watanabe Y, Moriya Yet al.. jPOST environment accelerates the reuse and reanalysis of public proteome mass spectrometry data. Nucleic Acids Res, 2025, 53: D462-D467

[51] Owens DM, Keyse SM. Differential regulation of MAP kinase signalling by dual-specificity protein phosphatases. Oncogene, 2007, 26: 3203-3213

[52] Pellegrin F, Nandris D, Chrestin H, Duran-Vila N. Rubber tree (Hevea brasiliensis) bark necrosis sy-ndrome I: still no evidence of a biotic causal agent. Plant Dis, 2004, 88: 1046

[53] Priya P, Venkatachalam P, Thulaseedharan A. Differential expression pattern of rubber elongation factor (REF) mRNA transcripts from high and low yielding clones of rubber tree (Hevea brasiliensis Muell. Arg.). Plant Cell Rep, 2007, 26: 1833-1838

[54] Qin YX, Shi F, Tang CR. Molecular characterization and expression analysis of cDNAs encoding four Rab and two Arf GTPases in the latex of Hevea brasiliensis. Plant Physiol Biochem, 2011, 49: 729-737

[55] Qin YX, Yang JH, Wang ZY, Fang Y, Hu B, Xiao Xet al.. The role of citrate synthase HbCS4 in latex regeneration of Hevea brasiliensis (para rubber tree). Ind Crops Prod, 2023, 206: 117637

[56] Richter S, Voss U, Jürgens G. Post-golgi traffic in plants. Traffic, 2009, 10: 819-828

[57] Sang T, Chen CW, Lin Z, Ma Y, Du Y, Lin PYet al.. DIA-based phosphoproteomics identifies early phosphorylation events in response to EGTA and mannitol in Arabidopsis. Mol Cell Proteomics, 2024, 23: 100804

[58] Sang T, Zhang Z, Liu GT, Wang PC. Navigating the landscape of plant proteomics. J Integr Plant Biol, 2025, 67: 740-761

[59] Schaller H, Grausem B, Benveniste P, Chye ML, Tan CTet al.. Expression of the Hevea brasiliensis (H.B.K) Mull. Arg. 3-hydroxy-3-methylglutaryl-coenzyme A reductase 1 in tobacco results in sterol overproduction. Plant Physiol, 1995, 109: 761-770

[60] Southorn WA. Thread-like reticulum in latex from Hevea brasiliensis and its relation to latex particles. Nature, 1961, 189: 1000-1001

[61] Southorn WA, Yip E. Latex flow studies. III. Electrostatic considerations in the colloidal stability of fresh hevea latex. J Rubber Res Inst Malaya, 1968, 20: 201-215

[62] Sugiyama N, Nakagami H, Mochida K, Daudi A, Tomita M, Shirasu Ket al.. Large-scale phosphorylation mapping reveals the extent of tyrosine phosphorylation in Arabidopsis. Mol Syst Biol, 2008, 4: 193

[63] Suzuki H, MacDonald J, Syed K, Salamov A, Hori C, Aerts Aet al.. Comparative genomics of the white-rot fungi, Phanerochaete carnosa and P. chrysosporium, to elucidate the genetic basis of the distinct wood types they colonize. BMC Genomics, 2012, 13: 444

[64] Tan JJ, Feng HQ, Deng Z. DIA-based phosphoproteomics profiling reveals rapid cold-responsive phosphorylation events in tomato fruit. Postharvest Biol Technol, 2025, 219: 113207

[65] Tang CR, Yang M, Fang YJ, Luo Y, Gao S, Xiao Xet al.. The rubber tree genome reveals new insights into rubber production and species adaptation. Nat Plants, 2016, 2: 1-10

[66] Tang DD, Chen MJ, Huang XH, Zhang GH, Zeng L, Zhang Get al.. Srplot: a free online platform for data visualization and graphing. PLoS ONE, 2023, 18: e0294236

[67] Tetlow IJ, Wait R, Lu Z, Akkasaeng R, Bowsher CG, Esposito Set al.. Protein phosphorylation in amyloplasts regulates starch branching enzyme activity and protein-protein interactions. Plant Cell, 2004, 16: 694-708

[68] Tsai C, Hsu C, Hung J, Wang Y, Choong W, Zeng Met al.. Sequential phosphoproteomic enrichment through complementary metal-directed immobilized metal ion affinity chromatography. Anal Chem, 2014, 86: 685-693

[69] Wadeesirisak K, Castano S, Vaysse L, Bonfils F, Peruch F, Rattanaporn Ket al.. Interactions of REF1 and SRPP1 rubber particle proteins from Hevea brasiliensis with synthetic phospholipids: effect of charge and size of lipid headgroup. Biochem Biophys Res Commun, 2023, 679: 205-214

[70] Wang XC, Shi MJ, Wang D, Chen Y, Cai F, Zhang Set al.. Comparative proteomics of primary and secondary lutoids reveals that chitinase and glucanase play a crucial combined role in rubber particle aggregation in Hevea brasiliensis. J Proteome Res, 2013, 12: 5146-5159

[71] Wang XC, Wang D, Sun Y, Yang Q, Chang L, Wang Let al.. Comprehensive proteomics analysis of laticifer latex reveals new insights into ethylene stimulation of natural rubber production. Sci Rep, 2015, 5: 13778

[72] Wang D, Xie QL, Sun Y, Tong Z, Chang L, Yu Let al.. Proteomic landscape has revealed small rubber particles are crucial rubber biosynthetic machines for ethylene-stimulation in natural rubber production. Int J Mol Sci, 2019, 20: 5082

[73] Wang C, Guo HB, He XW, Zhang S, Wang J, Guo XQet al.. Scaffold protein GhMORG1 enhances the resistance of cotton to Fusarium oxysporum by facilitating the MKK6-MPK4 cascade. Plant Biotechnol J, 2020, 18: 1421-1433

[74] Wang N, Wang XR, Qi Q, Iqbal A, Zhang H, Shi Jet al.. Analysis of the effects of mepiquat chloride priming on the seedling growth-promoting in cotton under salt stress by multi-omics. Ind Crops Prod, 2022, 186: 115296

[75] Wang YH, Kang GJ, Li Y, Qin HD, Zeng RZ (2018) Expression analysis of key genes REF, SRPP, HRT1 and HRT2 related to rubber biosynthesis of Hevea brasiliensis. Genom Appl Biol 37:3933–3943. https://doi.org/10.13417/j.gab.037.003933

[76] Wu SH, Zhang SX, Yang SG, Tian WM. Molecular cloning and expression of HbMKK4 from Hevea brasiliensis (para rubber tree). Acta Bot Boreali-Occident Sin, 2015, 35: 440-446

[77] Wu CT, Lan L, Li Y, Nie ZY, Zeng RZ. The relationship between latex metabolism gene expression with rubber yield and related traits in Hevea brasiliensis. BMC Genomics, 2018, 19: 897

[78] Wu BS, Gao L, Sun Y, Wu M, Wang D, Wei Jet al.. Comparative proteomic analysis of the Hevea brasiliensis latex under ethylene and calcium stimulation. Protein Pept Lett, 2019, 26: 834-847

[79] Xu TC, Li YM, Liu X, Yang X, Huang Z, Xing Jet al.. Rubber biosynthesis drives the biogenesis and development of rubber particles, the rubber-producing organelles. Plant Biotechnol J, 2025, 23: 2303-2316

[80] Yamashita S, Yamaguchi H, Waki T, Aoki Y, Mizuno M, Yanbe Fet al.. Identification and reconstitution of the rubber biosynthetic machinery on rubber particles from Hevea brasiliensis. Elife, 2016, 5: e19022

[81] Yang H, Dai LJ, Liu MY, Fan XK, Lu L, Guo Bet al.. Integrative analysis of transcriptome and metabolome reveals how ethylene increases natural rubber yield in Hevea brasiliensis. Front Plant Sci, 2024, 15: 1444693

[82] Yang SG, Du XY, Deng X, Tian W, Chao JQ (2024) Comparative phosphoproteome analyses provide new insight into ethephon stimulated latex yield in rubber tree. Trop Plants 3:e001. https://doi.org/10.48130/tp-0024-0002

[83] Zander M, Lewsey MG, Clark NM, Yin L, Bartlett A, Saldierna Guzmán JPet al.. Integrated multi-omics framework of the plant response to jasmonic acid. Nat Plants, 2020, 6: 290-302

[84] Zhang WJ, Zhou YW, Zhang Y, Su YH, Xu TD. Protein phosphorylation: a molecular switch in plant signaling. Cell Rep, 2023, 42: 112729

[85] Zhu JH, Zhang ZL. Ethylene stimulation of latex production in Hevea brasiliensis. Plant Signal Behav, 2009, 4: 1072-1074

[86] Zou Z, Zheng YJ, Xiao YH, Liu HY, Huang J, Zhao YG (2024) Molecular insights into PIP aquaporins in tigernut (Cyperus esculentus L.), a cyperaceae tuber plant. Trop Plants 3:e027. https://doi.org/10.48130/tp-0024-0030

Funding

the National Key Research and Development Program of China(2021YFA1300401)

the Key Technology Research and Development Program for Rubber-producing Plants of Guizhou Province(No. BQW2024-005)

the Special Fund for Talents of Guizhou University(LJ2023-03)

RIGHTS & PERMISSIONS

The Author(s)