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

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 °C for 2 h) and LT24 (cold treatment at 4 °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.