Comparative Transcriptomic Analyses Reveal Key Pathways in Response to Cold Stress at the Germination Stage of Quinoa (<i>Chenopodium quinoa</i> Willd.) Seeds

Quinoa (<i>Chenopodium quinoa</i> Willd.) has been widely grown as a cash crop. However, the molecular mechanism by which it responds to cold stress at the seed germination stage is still largely unknown. In this study, we performed a comparative transcriptomic analysis between the cold-tolerant cultivar XCq and cold-sensitive cultivar QCq in response to cold stress. A total number of 4552 and 4845 differentially expressed genes (<i>DEGs</i>) were identified in XCq and QCq upon the treatment of cold stress, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis demonstrated that the mitogen-activated protein kinase (MAPK) signaling pathway was identified only among the up-regulated DEGs in XCq.The expression of DEGs, which encoding transcription factors, such as AP2/ERF, bHLH, bZIP, MYB, ICEs, and CORs related to cold response, were higher in XCq than in QCq in response to cold stress. Weighted gene co-expression network analysis (WGCNA) showed that DEGs clustered in the co-expression modules positively correlated with the factors of quinoa variety and temperature were significantly enriched in the oxidative phosphorylation metabolic pathway. Further physiochemical analyses showed that the activities of superoxide dismutase and peroxidase as well as the contents of soluble protein and sugar, were significantly higher in XCq than in QCq. In summary, MAPK signaling and oxidative metabolism were the key pathways in quinoa upon cold stress. Our findings revealed that the enhanced activities of antioxidant enzymes alleviate the lipid peroxidation of membranes and promote the accumulation of osmotic adjustment substances, thereby enabling seeds to better resist oxidative damage under cold stress.