Abstract
Cold tolerance is a complex trait that requires a critical perspective to understand its underpinning mechanism. To unravel the molecular framework underlying maize (Zea mays L.) cold stress tolerance, we conducted a comparative transcriptome profiling of 24 cold-tolerant and 22 cold-sensitive inbred lines affected by cold stress at the seedling stage. Using the RNA-seq method, we identified 2237 differentially expressed genes (DEGs), namely 1656 and 581 annotated and unannotated DEGs, respectively. Further analysis of the 1656 annotated DEGs mined out two critical sets of cold-responsive DEGs, namely 779 and 877 DEGs, which were significantly enhanced in the tolerant and sensitive lines, respectively. Functional analysis of the 1656 DEGs highlighted the enrichment of signaling, carotenoid, lipid metabolism, transcription factors (TFs), peroxisome, and amino acid metabolism. A total of 147 TFs belonging to 32 families, including MYB, ERF, NAC, WRKY, bHLH, MIKC MADS, and C(2)H(2), were strongly altered by cold stress. Moreover, the tolerant lines' 779 enhanced DEGs were predominantly associated with carotenoid, ABC transporter, glutathione, lipid metabolism, and amino acid metabolism. In comparison, the cold-sensitive lines' 877 enhanced DEGs were significantly enriched for MAPK signaling, peroxisome, ribosome, and carbon metabolism pathways. The biggest proportion of the unannotated DEGs was implicated in the roles of long non-coding RNAs (lncRNAs). Taken together, this study provides valuable insights that offer a deeper understanding of the molecular mechanisms underlying maize response to cold stress at the seedling stage, thus opening up possibilities for a breeding program of maize tolerance to cold stress.