Integrated Analysis of Metabolome and Transcriptome Reveals Insights for Low Phosphorus Tolerance in Wheat Seedling

Abstract

Low phosphorus (LP) stress leads to a significant reduction in wheat yield, primarily in the reduction of biomass, the number of tillers and spike grains, the delay in heading and flowering, and the inhibition of starch synthesis and grouting. However, the differences in regulatory pathway responses to low phosphorus stress among different wheat genotypes are still largely unknown. In this study, metabolome and transcriptome analyses of G28 (LP-tolerant) and L143 (LP-sensitive) wheat varieties after 72 h of normal phosphorus (CK) and LP stress were performed. A total of 181 and 163 differentially accumulated metabolites (DAMs) were detected for G28CK vs. G28LP and L143CK vs. L143LP, respectively. Notably, the expression of pilocarpine (C07474) in G28CK vs. G28LP was significantly downregulated 4.77-fold, while the expression of neochlorogenic acid (C17147) in L143CK vs. L143LP was significantly upregulated 2.34-fold. A total of 4023 differentially expressed genes (DEGs) were acquired between G28 and L143, of which 1120 DEGs were considered as the core DEGs of LP tolerance of wheat after LP treatment. The integration of metabolomics and transcriptomic data further revealed that the LP tolerance of wheat was closely related to 15 metabolites and 18 key genes in the sugar and amino acid metabolism pathway. The oxidative phosphorylation pathway was enriched to four ATPases, two cytochrome c reductase genes, and fumaric acid under LP treatment. Moreover, PHT1;1, TFs (ARFA, WRKY40, MYB4, MYB85), and IAA20 genes were related to the Pi starvation stress of wheat roots. Therefore, the differences in LP tolerance of different wheat varieties were related to energy metabolism, amino acid metabolism, phytohormones, and PHT proteins, and precisely regulated by the levels of various molecular pathways to adapt to Pi starvation stress. Taken together, this study may help to reveal the complex regulatory process of wheat adaptation to Pi starvation and provide new genetic clues for further study on improving plant Pi utilization efficiency.