Hevea brasiliensis, a tropical tree species from the Amazon rainforest, is the main source of natural rubber worldwide. Due to the high pressure of fungal diseases in hot, humid regions, rubber plantations have been moved to escape areas, which are dryer and have lower temperatures during the winter. Here, we combined gene expression data of a primary (GT1) and a secondary (RRIM600) young rubber tree clones, which present different cold tolerance strategies, to analyze rubber tree gene expression regulation during 24 h of cold exposure (10 degrees Celsius). Together with traditional differential expression approaches, a RNA sequencing (RNA-seq) gene coexpression network (GCN) comprising 27,220 genes was established in which the genes were grouped into 832 clusters. In the GCN, most of the rubber tree molecular responses to cold stress were grouped in 26 clusters, which were divided into three GCN modules: a downregulated group comprising 12 clusters and two upregulated groups comprising eleven and three clusters. Considering the three modules identified, the general Hevea response to short-term cold exposure involved downregulation of gibberellin (GA) signaling, complex regulation of jasmonic acid (JA) stress responses and programmed cell death (PCD) and upregulation of ethylene responsive genes. The hub genes of the cold-responsive modules were subsequently identified and analyzed. As a result of the GCN strategy applied in this study, we could not only access individual DEGs related to the Hevea cold response, but also provide insights into a deeper cascade of associated mechanisms involved in the response to cold stress in young rubber trees. Our results may represent the genetic stress responses of the species, developed during its evolution, since the varieties chosen for this work are genotypes that were selected during the early years of rubber tree domestication. The understanding of H. brasiliensis cold response mechanisms can greatly improve the breeding strategies for this crop, which has a narrow genetic base, is impacted by climate change and is the only source for large-scale rubber production.
HbSnRK2.6 Functions in ABA-Regulated Cold Stress Response by Promoting HbICE2 Transcriptional Activity in Hevea brasiliensis
