Transcriptome profiling of rubber tree (Hevea brasiliensis) discovers candidate regulators of the cold stress response

  • Gong, Xiao-Xiao (Hainan Key Laboratory for Sustainable Utilization of Tropical Bioresources, Tropical Agriculture and Forestry Institute, Hainan University) ;
  • Yan, Bing-Yu (Hainan Key Laboratory for Sustainable Utilization of Tropical Bioresources, Tropical Agriculture and Forestry Institute, Hainan University) ;
  • Hu, Jin (Hainan Key Laboratory for Sustainable Utilization of Tropical Bioresources, Tropical Agriculture and Forestry Institute, Hainan University) ;
  • Yang, Cui-Ping (Hainan Key Laboratory for Sustainable Utilization of Tropical Bioresources, Tropical Agriculture and Forestry Institute, Hainan University) ;
  • Li, Yi-Jian (Service Center of Science and Technology, Rubber Research Institute, Chinese Academy of Tropical Agricultural Sciences) ;
  • Liu, Jin-Ping (Hainan Key Laboratory for Sustainable Utilization of Tropical Bioresources, Tropical Agriculture and Forestry Institute, Hainan University) ;
  • Liao, Wen-Bin (Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences)
  • Received : 2017.10.27
  • Accepted : 2018.02.28
  • Published : 2018.11.30


Tropical plant rubber tree (Hevea brasiliensis) is the sole source of commercial natural rubber and low-temperature stress is the most important limiting factor for its cultivation. To characterize the gene expression profiles of H. brasiliensis under the cold stress and discover the key cold stress-induced genes. Three cDNA libraries, CT (control), LT2 (cold treatment at $4^{\circ}C$ for 2 h) and LT24 (cold treatment at $4^{\circ}C$ for 24 h) were constructed for RNA sequencing (RNA-Seq) and gene expression profiling. Quantitative real time PCR (qRT-PCR) was conducted to validate the RNA-Seq and gene differentially expression results. A total of 1457 and 2328 differentially expressed genes (DEGs) in LT2 and LT24 compared with CT were respectively detected. Most significantly enriched KEGG pathways included flavonoid biosynthesis, phenylpropanoid biosynthesis, plant hormone signal transduction, cutin, suberine and wax biosynthesis, Pentose and glucuronate interconversions, phenylalanine metabolism and starch and sucrose metabolism. A total of 239 transcription factors (TFs) were differentially expressed following 2 h or/and 24 h of cold treatment. Cold-response transcription factor families included ARR-B, B3, BES1, bHLH, C2H, CO-like, Dof, ERF, FAR1, G2-like, GRAS, GRF, HD-ZIP, HSF, LBD, MIKC-MADS, M-type MADS, MYB, MYB-related, NAC, RAV, SRS, TALE, TCP, Trihelix, WOX, WRKY, YABBY and ZF-HD. The genome-wide transcriptional response of rubber tree to the cold treatments were determined and a large number of DEGs were characterized including 239 transcription factors, providing important clues for further elucidation of the mechanisms of cold stress responses in rubber tree.


Supported by : National Natural Science Foundation of China


  1. Ahmed NU, Park JI, Jung HJ, Hur Y, Nou IS (2015) Anthocyanin biosynthesis for cold and freezing stress tolerance and desirable color in Brassica rapa. Funct Integr Genom 15:383-394
  2. An D, Yang J, Zhang P (2012) Transcriptome profiling of low temperature-treated Cassava apical shoots showed dynamic responses of tropical plant to cold stress. BMC Genom 13:64
  3. Azuma A, Yakushiji H, Koshita Y, Kobayashi S (2012) Flavonoid biosynthesis-related genes in grape skin are differentially regulated by temperature and light conditions. Planta 236:1067-1080
  4. Bahieldin A, Atef A, Edris S, Gadalla NO, Ali HM, Hassan SM, Al-Kordy MA, Ramadan AM, Makki RM, Al-Hajar AS et al (2016) Ethylene responsive transcription factor ERF109 retards PCD and improves salt tolerance in plant. BMC Plant Biol 16:216
  5. Brown RS (2005) Zinc finger proteins: getting a grip on RNA. Curr Opin Struct Biol 15:94-98
  6. Calzadilla PI, Maiale SJ, Ruiz OA, Escaray FJ (2016) Transcriptome response mediated by cold stress in Lotus japonicus. Front Plant Sci 7:374
  7. Chalker-Scott L (1999) Environmental significance of anthocyanins in plant stress responses. Photochem Photobiol 70:1-9
  8. Chen B-Q, Cao J-H, Wang J-K, Wu Z-X, Tao Z-L, Chen J-M, Yang C-A, Xie G-S (2012) Estimation of rubber stand age in typhoon and chilling injury afflicted area with Landsat TM data: a case study in Hainan Island, China. For Ecol Manag 274:222-230
  9. Chinnusamy V, Zhu J, Zhu JK (2007) Cold stress regulation of gene expression in plants. Trends Plant Sci 12:444-451
  10. Chinnusamy V, Zhu JK, Sunkar R (2010) Gene regulation during cold stress acclimation in plants. Methods Mol Biol 639:39-55
  11. Christie PJ, Alfenito MR, Walbot V (1994) Impact of low-temperature stress on general phenylpropanoid and anthocyanin pathways: enhancement of transcript abundance and anthocyanin pigmentation in maize seedlings. Planta 194:541-549
  12. Croft D, O'Kelly G, Wu G, Haw R, Gillespie M, Matthews L, Caudy M, Garapati P, Gopinath G, Jassal B et al (2011) Reactome: a database of reactions, pathways and biological processes. Nucleic Acids Res 39:D691-D697
  13. Ding J, Wang L, Ruan C (2017) Comparative transcriptome analysis of lipid biosynthesis in seeds and non-seed tissues of sea buckthorn. Genes Genom 39:1021-1033
  14. Du Z, Zhou X, Ling Y, Zhang Z, Su Z (2010) agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res 38:W64-W70
  15. Eremina M, Rozhon W, Poppenberger B (2016) Hormonal control of cold stress responses in plants. Cell Mol Life Sci 73:797-810
  16. Fernandez-Calvo P, Chini A, Fernandez-Barbero G, Chico JM, Gimenez-Ibanez S, Geerinck J, Eeckhout D, Schweizer F, Godoy M, Franco-Zorrilla JM et al (2011) The Arabidopsis bHLH transcription factors MYC3 and MYC4 are targets of JAZ repressors and act additively with MYC2 in the activation of jasmonate responses. Plant Cell 23:701-715
  17. Figueiredo DD, Barros PM, Cordeiro AM, Serra TS, Lourenco T, Chander S, Oliveira MM, Saibo NJM (2012) Seven zinc-finger transcription factors are novel regulators of the stress responsive gene OsDREB1B. J Exp Bot 63:3643-3656
  18. Fowler S, Thomashow MF (2002) Arabidopsis transcriptome profiling indicates that multiple regulatory pathways are activated during cold acclimation in addition to the CBF cold response pathway. Plant Cell 14:1675-1690
  19. Gilmour SJ, Sebolt AM, Salazar MP, Everard JD, Thomashow MF (2000) Overexpression of the Arabidopsis CBF3 transcriptional activator mimics multiple biochemical changes associated with cold acclimation. Plant Physiol 24:1854-1865
  20. Gilmour SJ, Fowler SG, Thomashow MF (2004) Arabidopsis transcriptional activators CBF1, CBF2, and CBF3 have matching functional activities. Plant Mol Biol 54:767-781
  21. Gronover CS, Wahler D, Prufer D (2011) Natural rubber biosynthesis and physic-chemical studies on plant derived latex. In: Elnashar M (ed) Biotechnology of biopolymers. InTech, Croatia, pp 75-88
  22. Hannah MA, Heyer AG, Hincha DK (2005) A global survey of gene regulation during cold acclimation in Arabidopsis thaliana. PLoS Genet 1:e26
  23. He J, Yang H-J, Lin M-X (1982) Studies on cellular injuries in rubber trees under chilling stress. Acta Sci Nat Univ Amoi 21:84-91
  24. Hu Y, Jiang Y, Han X, Wang H, Pan J, Yu D (2017) Jasmonate regulates leaf senescence and tolerance to cold stress: crosstalk with other phytohormones. J Exp Bot 68:1361-1369
  25. Hua W, Kong W, Cao X-Y, Chen C, Liu Q, Li X, Wang Z (2017) Transcriptome analysis of Dioscorea zingiberensis identifies genes involved in diosgenin biosynthesis. Genes Genom 39:509-520
  26. Huang ZD, Pan YQ (1992) Rubber cultivation under climatic stresses in China. In: Sethuraj MR, Mathew NM (eds) Natural rubber: biology, cultivation and technology. Elsevier, Dordrecht, pp 220-238
  27. Jaradat MR, Feurtado JA, Huang D, Lu Y, Cutler AJ (2013) Multiple roles of the transcription factor AtMYBR1/AtMYB44 in ABA signaling, stress responses, and leaf senescence. BMC Plant Biol 13:192
  28. Jin J, Zhang H, Zhang J, Liu P, Chen X, Li Z, Xu Y, Lu P, Cao P (2017) Integrated transcriptomics and metabolomics analysis to characterize cold stress responses in Nicotiana tabacum. BMC Genom 18:496
  29. Kan L-Y, Xie G-S, Tao Z-L, Yang L-F, Cui Z-F (2009) Analysis on rubber tree cold injury in 2007/2008 winter in Hainan. Chin Agric Sci Bull 25:251-257
  30. Kanehisa M, Goto S (2000) KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 28:27-30
  31. Kaplan F, Kopka J, Sung DY, Zhao W, Popp M, Porat R, Guy CL (2007) Transcript and metabolite profiling during cold acclimation of Arabidopsis reveals an intricate relationship of cold-regulated gene expression with modifications in metabolite content. Plant J 50:967-981
  32. Karp PD, Ouzounis CA, Moore-Kochlacs C, Goldovsky L, Kaipa P, Ahren D, Tsoka S, Darzentas N, Kunin V, Lopez-Bigas N (2005) Expansion of the BioCyc collection of pathway/genome databases to 160 genomes. Nucleic Acids Res 33:6083-6089
  33. Khaledian Y, Maali-Amiri R, Talei A (2015) Phenylpropanoid and antioxidant changes in chickpea plants during cold stress. Russ J Plant Physiol 62:772-778
  34. Knight MR, Knight H (2012) Low-temperature perception leading to gene expression and cold tolerance in higher plants. New Phytol 195:737-751
  35. Lee BH, Henderson DA, Zhu JK (2005) The Arabidopsis coldresponsive transcriptome and its regulation by ICE1. Plant Cell 17:3155-3175
  36. Li J, Yang S, Gai J (2017) Transcriptome comparative analysis between the cytoplasmic male sterile line and fertile line in soybean (Glycine max (L.) Merr.). Genes Genom 39:1117-1127
  37. Lin M-X, Yang H-J (1994) Physiological responses of Hevea brasiliensis during the chilling injury. Chin J Trop Crops 15:7-11
  38. Lindemose S, O'Shea C, Jensen MK, Skriver K (2013) Structure, function and networks of transcription factors involved in abiotic stress responses. Int J Mol Sci 14:5842-5878
  39. Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, Shinozaki K (1998) Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain, separate two cellular signal transduction pathways in drought- and low temperature responsive gene expression, respectively, in Arabidopsis. Plant Cell 10:1391-1406
  40. Liu J-P, Xia Z-Q, Tian X-Y, Li Y-J (2015) Transcriptome sequencing and analysis of rubber tree (Hevea brasiliensis Muell.) to discover putative genes associated with tapping panel dryness (TPD). BMC Genom 16:398
  41. Liu J-P, Zhuang Y-F, Guo X-L, Li Y-J (2016) Molecular mechanism of ethylene stimulation of latex yield in rubber tree (Hevea brasiliensis) revealed by de novo sequencing and transcriptome analysis. BMC Genom 17:257
  42. Luo P, He J-J, Yao Y-L, Dai X-H, Cheng R-X, Wang L-F (2012) Differential responses of two rubber tree clones to chilling stress. Afr J Biotech 11:13466-13471
  43. Mai J, Herbette S, Vandame M, Kositsup B, Kasemsap P, Cavaloc E, Julien J, Ameglio T, Roeckel-Drevet P (2009) Effect of chilling on photosynthesis and antioxidant enzymes in Hevea brasiliensis Muell. Arg Trees 23:863-874
  44. Matthews L, Gopinath G, Gillespie M, Caudy M, Croft D, de Bono B, Garapati P, Hemish J, Hermjakob H, Jassal B et al (2009) Reactome knowledgebase of human biological pathways and processes. Nucleic Acids Res 37:D619-D622
  45. Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K (2012) AP2/ERF family transcription factors in plant abiotic stress responses. Biochim Biophys Acta 1819:86-96
  46. Niu Y, Figueroa P, Browse J (2011) Characterization of JAZ-interacting bHLH transcription factors that regulate jasmonate responses in Arabidopsis. J Exp Bot 62:2143-2154
  47. Park S, Lee CM, Doherty CJ, Gilmour SJ, Kim Y, Thomashow MF (2015) Regulation of the Arabidopsis CBF regulon by a complex low-temperature regulatory network. Plant J 82:193-207
  48. Persak H, Pitzschke A (2014) Dominant repression by Arabidopsis transcription factor MYB44 causes oxidative damage and hypersensitivity to abiotic stress. Int J Mol Sci 15:2517-2537
  49. Petrussa E, Braidot E, Zancani M, Peresson C, Bertolini A, Patui S, Vianello A (2013) Plant flavonoids-biosynthesis, transport and involvement in stress responses. Int J Mol Sci 14:14950-14973
  50. Priyadarshan PM (2011) Biology of hevea rubber. CAB International, Wallingford, pp 1-163
  51. Provart NJ, Gil P, Chen W, Han B, Chang HS, Wang X, Zhu T (2003) Gene expression phenotypes of Arabidopsis associated with sensitivity to low temperatures. Plant Physiol 132:893-906
  52. Pushparajah E (1983) Problems and potentials for establishing Hevea under difficult environmental conditions. Planter 59:242-251
  53. Ramamoorthy R, Jiang SY, Kumar N, Venkatesh PN, Ramachandran S (2008) A comprehensive transcriptional profiling of the WRKY gene family in rice under various abiotic and phytohormone treatments. Plant Cell Physiol 49:865-879
  54. Roberts A, Pimentel H, Trapnell C, Pachter L (2011) Identification of novel transcripts in annotated genomes using RNA-SEq. Bioinformatics 27:2325-2329
  55. Roy S (2016) Function of MYB domain transcription factors in abiotic stress and epigenetic control of stress response in plant genome. Plant Signal Behav 11:e1117723
  56. Schaefer CF, Anthony K, Krupa S, Buchoff J, Day M, Hannay T, Buetow KH (2009) PID: the pathway interaction database. Nucleic Acids Res 37:D674-D679
  57. Schulz E, Tohge T, Zuther E, Fernie AR, Hincha DK (2016) Flavonoids are determinants of freezing tolerance and cold acclimation in Arabidopsis thaliana. Sci Rep 6:34027
  58. Seo J, Sohn H, Noh K, Jung C, An J, Donovan C, Somers D, Kim D, Jeong S-C, Kim C-G et al (2012) Expression of the Arabidopsis AtMYB44 gene confers drought/salt-stress tolerance in transgenic soybean. Mol Breed 29:601-608
  59. Sharma M, Laxmi A (2016) Jasmonates: emerging players in controlling temperature stress tolerance. Front Plant Sci 6:1129
  60. Solecka D (1997) Role of phenylpropanoid compounds in plant responses to different stress factors. Acta Physiol Plant 19:257-268
  61. Song S, Qi T, Fan M, Zhang X, Gao H, Huang H, Wu D, Guo H, Xie D (2013) The bHLH subgroup IIId factors negatively regulate jasmonate-mediated plant defense and development. PLoS Genet 9:e1003653
  62. Stockinger EJ, Gilmour SJ, Thomashow MF (1997) Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcription activator that binds to the C repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit. Proc Natl Acad Sci USA 94:1035-1040
  63. Strickler SR, Bombarely A, Mueller LA (2012) Designing a transcriptome next generation sequencing project for a nonmodel plant species. Am J Bot 99:257-266
  64. Tang C, Yang M, Fang Y, Luo Y, Gao S, Xiao X, An Z, Zhou B, Zhang B, Tan X et al (2016) The rubber tree genome reveals new insights into rubber production and species adaptation. Nat Plants 2:16073
  65. Theocharis A, Clement C, Barka EA (2012) Physiological and molecular changes in plants grown at low temperatures. Planta 235:1091-1105
  66. Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, Diemer K, Muruganujan A, Narechania A (2003) PANTHER: a library of protein families and subfamilies indexed by function. Genome Res 13:2129-2141
  67. Thomashow MF (1999) Plant cold acclimation freezing tolerance genes and regulatory mechanisms. Annu Rev Plant Physiol 50:571-599
  68. Tian Y-H, Yuan H-F, Xie J, Deng J-W, Dao X-S, Zheng Y-L (2016) Effect of diurnal irradiance on night-chilling tolerance of six rubber cultivars. Photosynthetica 54:374-380
  69. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L (2012a) Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc 7:562-578
  70. Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L (2012b) Differential analysis of gene regulation at transcript resolution with RNA-sEq. Nat Biotechnol 31:46-53
  71. Van Buskirk HA, Thomashow MF (2006) Arabidopsis transcription factors regulating cold acclimation. Physiol Plant 126:72-80
  72. Vannini C, Locatelli F, Bracale M, Magnani E, Marsoni M, Osnato M, Mattana M, Baldoni E, Coraggio I (2004) Overexpression of the rice Osmyb4 gene increases chilling and freezing tolerance of Arabidopsis thaliana plants. Plant J 37:115-127
  73. Varshney RK, Nayak SN, May GD, Jackson SA (2009) Next-generation sequencing technologies and their implications for crop genetics and breeding. Trends Biotechnol 27:522-530
  74. Vogel JT, Zarka DG, Van Buskirk HA, Fowler SG, Thomashow MF (2005) Roles of the CBF2 and ZAT12 transcription factors in configuring the low temperature transcriptome of Arabidopsis. Plant J 41:195-211
  75. Wang Y-R, Liu H-X, Guo Z-Y (1978) Effects of chilling temperature of plant metabolism of Hevea brasiliensis. Acta Bot Sin 20:44-53
  76. Wang X-J, Li W-G, Gao X-S, Wu C-T, Zhang X-F (2012) Physiological characteristics of Hevea brasiliensis in response to low temperature stress and its regulation mechanisms. Plant Physiol J 48:318-324
  77. Wang L, Zhu W, Fang L, Sun X, Su L, Liang Z, Wang N, Londo JP, Li S, Xin H (2014) Genome-wide identification of WRKY family genes and their response to cold stress in Vitis vinifera. BMC Plant Biol 14:103
  78. Wang H, Wang H, Shao H, Tang X (2016) Recent advances in utilizing transcription factors to improve plant abiotic stress tolerance by transgenic technology. Front Plant Sci 7:67
  79. Wang P, Ma L, Li Y, Wang S, Li L, Yang R (2017a) Transcriptome analysis reveals sunflower cytochrome P450 CYP93A1 responses to high salinity treatment at the seedling stage. Genes Genom 39:581-591
  80. Wang X, Zhao M, Wu W, Korir NK, Qian Y, Wang Z (2017b) Comparative transcriptome analysis of berry-sizing effects of gibberellin (GA3) on seedless Vitis vinifera L. Genes Genom 39:493-507
  81. Wei Y, Shi H, Xia Z, Tie W, Ding Z, Yan Y, Wang W, Hu W, Li K (2016) Genome-wide identification and expression analysis of the WRKY gene family in cassava. Front Plant Sci 7:25
  82. Wu ZG, Jiang W, Chen SL, Mantri N, Tao ZM, Jiang CX (2016) Insights from the cold transcriptome and metabolome of Dendrobium officinale: global reprogramming of metabolic and gene regulation networks during cold acclimation. Front Plant Sci 7:1653
  83. Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li CY, Wei L (2011) KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res 39:W316-322
  84. Yang QS, Gao J, He WD, Dou TX, Ding LJ, Wu JH, Li CY, Peng XX, Zhang S, Yi GJ (2015) Comparative transcriptomics analysis reveals difference of key gene expression between banana and plantain in response to cold stress. BMC Genom 16:446
  85. Zeng C, Chen Z, Xia J, Zhang K, Chen X, Zhou Y, Bo W, Song S, Deng D, Guo X et al (2014) Chilling acclimation provides immunity to stress by altering regulatory networks and inducing genes with protective functions in cassava. BMC Plant Biol 14:207
  86. Zhan X, Yang L, Wang D, Zhu JK, Lang Z (2016) De novo assembly and analysis of the transcriptome of Ocimum americanum var. pilosum under cold stress. BMC Genom 17:209
  87. Zhang Y, Yu H, Yang X, Li Q, Ling J, Wang H, Gu X, Huang S, Jiang W (2016) CsWRKY46, a WRKY transcription factor from cucumber, confers cold resistance in transgenic-plant by regulating a set of cold-stress responsive genes in an ABA-dependent manner. Plant Physiol Biochem 108:478-487
  88. Zhao C, Lang Z, Zhu JK (2015) Cold responsive gene transcription becomes more complex. Trends Plant Sci 20:466-468
  89. Zhou Q-Y, Tian A-G, Zou H-F, Xie ZM, Lei G, Huang J, Wang CM, Wang HW, Zhang JS, Chen SY (2008) Soybean WRKY type transcription factor genes, GmWRKY13, GmWRKY21, and GmWRKY54, confer differential tolerance to abiotic stresses in transgenic Arabidopsis plants. Plant Biotech J 6:486-503