Abstract
Two contrasting cold response rubber tree clones, the cold-resistant '93-114' and cold-sensitive 'Reken501', were subject to a global transcriptome response assessing via high-throughput RNA-seq technique and comprehensive bioinformatics analysis using the referenced rubber tree genome with the purpose of exploring the potential molecular cues underlying the tolerance of rubber trees to cold stress. As a result, a total of 1919 genes had significantly higher expression, while 2929 genes had significantly lower expression in '93-114' than in 'Reken501' without cold stress. Upon cold stress, the numbers of genes with significantly higher expression decreased to 1501 at 1 h treatment and to 1285 at 24 h treatment in '93-114' than that of 'Reken501', conversely, the numbers of genes with significantly lower expression increased to 7567 at 1 h treatment and to 5482 at 24 h treatment. Functional annotation of the differentially expressed genes between '93-114' and 'Reken501' suggests that down-regulation of auxin and ethylene signaling and activation of heat shock module and ROS scavengers is a primary strategy for H. brasiliensis to cope with cold stress. Our identified vital differentially expressed genes may be beneficial for elucidation of the molecular mechanisms underlying cold tolerance and for genetic improvement of H. brasiliensis clones.