Global Transcriptome and Weighted Gene Co-Expression Network Analyses of Cold Stress Responses in Chinese Cabbage

TL;DR

This study identifies genes and pathways in Chinese cabbage that respond to cold stress, which could help improve cold tolerance in crops.

Contribution

The study provides novel insights into cold-stress-responsive genes and co-expression modules in Chinese cabbage.

Findings

1832 differentially expressed genes were identified across cold stress conditions.

Key pathways affected by cold stress include sulfur metabolism and photosynthesis.

Four co-expression modules were linked to vernalization-related processes.

Abstract

Background/Objectives: Chinese cabbage (Brassica rapa ssp. Pekinensis, AA) growth and development is highly sensitive to cold temperatures. Prolonged low-temperature exposure during early growth stages can induce premature bolting, which reduces market quality and yield. Methods: Here, using comparative leaf RNA-seq transcriptome analysis of plants grown at 6, 9, 12, and 15 °C, we explored key genes and metabolic pathways regulating Chinese cabbage cold response. Results: RNA-seq transcriptome analysis identified a total of 1832 differentially expressed genes (DEGs) in the three comparison groups, with 5452, 1861, and 752 DEGs specifically expressed in the A6_vs_A15, A9_vs_A15, and A12_vs_A15 groups, respectively. KEGG enrichment analysis of DEGs showed that sulfur metabolism, secondary metabolites biosynthesis and photosynthesis pathways were mostly affected by cold stress. K-means clustering revealed distinct expression profiles among the DEGs enriched in cold stress response-associated clusters. Subsequently, DEGs were divided into 18 modules by WGCNA, whereupon co-expression genes that clustered into similar modules exhibited diverse expression and were annotated to various GO terms at different temperatures. Module-trait association analysis revealed M1, M2, M3, and M6 modules as key clusters potentially linked to vernalization-related processes. These modules harbored candidate hub genes encoding transcription factors (including MYB, bZIP, and WRKY), protein kinases, and cold-stress-responsive genes. Additionally, phenotypic analysis showed that 12 °C to 15 °C supported optimal growth, whereas <9 °C temperature inhibited growth. Physiological measurements showed increased antioxidant enzyme activity and proline accumulation at 6 °C. Conclusions: Overall, our study provides a set of candidate cold-stress-responsive genes and co-expression modules that may support cold stress tolerance breeding in Chinese cabbage.

Figures