Integration of Metabolome and Transcriptome Reveals the Major Metabolic Pathways and Potential Biomarkers in Response to Freeze-Stress Regulation in Apple (Malus domestica)

Freezing stress is the main factor affecting the normal growth and distribution of plants. The safe overwintering of a perennial deciduous plant is a crucial link to ensuring its survival and yield. However, little is known about the molecular mechanism of its gene regulation metabolites as related to its freeze-tolerance. In order to enhance our comprehension of freeze-tolerance metabolites and gene expression in dormant apple trees, we examined the metabolic and transcriptomic differences between ‘Ralls’ and ‘Fuji’, two apple varieties with varying degrees of resistance to freezing. The results of the freezing treatment showed that ‘Ralls’ had stronger freeze-tolerance than ‘Fuji’. We identified 302, 334, and 267 up-regulated differentially accumulated metabolites (DAMs) and 408, 387, and 497 down-regulated DAMs between ‘Ralls’ and ‘Fuji’ under −10, −15, and −20 °C treatment, respectively. A total of 359 shared metabolites were obtained in the upward trend modules, of which 62 metabolites were associated with 89 pathways. The number of up-regulated genes accounted for 50.2%, 45.6%, and 43.2% of the total number of differentially expressed genes (DEGs), respectively, at −10, −15, and −20 °C. Through combined transcriptome and metabolome analysis, we identified 12 pathways that included 16 DAMs and 65 DEGs. Meanwhile, we found that 20 DEGs were identified in the phenylpropanoid biosynthesis pathway and its related pathways, involving the metabolism of p-Coumaroyl-CoA, 7, 4′-Dihydroxyflavone, and scolymoside. These discoveries advance our comprehension of the molecular mechanism underlying apple freeze-tolerance and provide genetic material for breeding apple cultivars with enhanced freeze-tolerance.