Comparative transcriptomics and co-expression networks reveal cultivar-specific molecular signatures associated with reproductive-stage cold stress in rice.

key messageThe resistance of Huaidao5 results from the high constitutive expression of tolerance genes, while that of Huaidao9 is due to the cold-induced resistance in flag leaves and panicles. The regulation mechanism of rice seedlings' cold tolerance is relatively clear, and knowledge of its underlying mechanisms at the reproductive stage is limited. We performed differential expression and co-expression network analyses to transcriptomes from panicle and flag leaf tissues of a cold-tolerant cultivar (Huaidao5), and a sensitive cultivar (Huaidao9), under reproductive-stage cold stress. The results revealed that the expression levels of genes in stress-related pathways such as MAPK signaling pathway, diterpenoid biosynthesis, glutathione metabolism, plant-pathogen interaction and plant hormone signal transduction were constitutively highly expressed in Huaidao5, especially in panicles. Moreover, the Hudaidao5's panicle sample-specific (under cold) module contained some genes related to rice yield, such as GW5L, GGC2, SG1 and CTPS1. However, the resistance of Huaidao9 was derived from the induced resistance to cold in flag leaves and panicles. In the flag leaves, the responses included a series of stress response and signal transduction, while in the panicles nitrogen metabolism was severely affected, especially 66 endosperm-specific genes. Through integrating differential expression with co-expression networks, we predicted 161 candidate genes (79 cold-responsive genes common to both cultivars and 82 cold-tolerance genes associated with differences in cold tolerance between cultivars) potentially affecting cold response/tolerance, among which 85 (52.80%) were known to be cold-related genes. Moreover, 52 (65.82%) cold-responsive genes (e.g., TIFY11C, LSK1 and LPA) could be confirmed by previous transcriptome studies and 72 (87.80%) cold-tolerance genes (e.g., APX5, OsFbox17 and OsSTA109) were located within QTLs associated with cold tolerance. This study provides an efficient strategy for further discovery of mechanisms of cold tolerance in rice.