Transcriptomic and weighted gene co-expression network analysis of tropic and temperate maize inbred lines recovering from heat stress

Heat stress (HS) causes imbalance of cellular homeostasis, growth impairment and extensively yield loss in crop production. In the present study, the tropic maize inbred CIMBL55 showed more thermotolerance than the maize temperate inbred B73, with less leaf damage rate and ROS accumulation. Transcriptome profiling of CIMBL55 and B73 upon (exposing at 45 ℃ for 0, 1, and 6 h) and post (recovering at 28 ℃ for 1 and 6 h) HS were further assessed and a total of 20204 DEGs were identified. Functional annotation revealed that HS activated unfolded protein response in endoplasmic reticulum in both two inbreds. Moreover, in CIMBL55, far more primary and secondary metabolism pathways were transcriptional altered. Afterwards, weighted gene co-expression analysis grouped all expressed genes into eighteen co-expressed modules. Four HS responsive and four CIMBL55 recovery-related modules were subsequently identified. Highly connected genes (hub genes) in these modules were characterized as transcription factors, heat shock proteins, Ca2+ signaling related genes and various enzymes. Moreover, one hub gene, ZmHsftf13 was verified to positively regulate thermotolerance by heterologous expressing in Arabidopsis and its Mu insertion mutant. The present research provides promising genes related to HS response in maize and is of great significance for breeding.