Transcriptomic investigation unveils the role of energy metabolism under low phosphorus and salt combined stress in soybean (Glycine max)

Soybean (Glycine max) is economically significant, but the mechanisms underlying its adaptation to simultaneous low phosphorus and salt stresses are unclear. We employed the Shennong 94-1-8 soybean germplasm to conduct a comprehensive analysis, integrating both physiochemical and transcriptomic approaches, to unravel the response mechanisms of soybean when subjected to simultaneous low phosphorus and salt stresses. Remarkably, the combined stress exhibited the most pronounced impact on the soybean root system, which led to a substantial reduction in total soluble sugar (TSS) and total soluble protein (TSP) within the plants under this treatment. A total of 20,953 differentially expressed genes were identified through pairwise comparisons. Heatmap analysis of genes related to energy metabolism pathways demonstrated a significant down-regulation in expression under salt and low phosphorus + salt treatments, while low phosphorus treatment did not exhibit similar expression trends. Furthermore, the weighted gene co-expression network analysis (WGCNA) indicated that the blue module had a strong positive correlation with TSS and TSP. Notably, 2,3-bisphosphoglycerate-dependent phosphoglycerate mutase 1, FCS-Like Zinc finger 8, auxin response factor 18 isoform X2, and NADP-dependent malic enzyme emerged as hub genes associated with energy metabolism. In summary, our findings indicate that soybean roots are more adversely affected by salt and combined stress than by low phosphorus alone due to reduced activity in energy metabolism-related pathways and hub genes. These results offer novel insights into the adaptive mechanisms of soybeans when facing the combined stress of low phosphorus and salinity.