Quantitative Phosphoproteomic Analysis Provides Insights into the Sodium Bicarbonate Responsiveness of Glycine max

Sodium bicarbonate stress caused by NaHCO₃ is one of the most severe abiotic stresses affecting agricultural production worldwide. However, little attention has been given to the molecular mechanisms underlying plant responses to sodium bicarbonate stress. To understand phosphorylation events in signaling pathways triggered by sodium bicarbonate stress, TMT-labeling-based quantitative phosphoproteomic analyses were performed on soybean leaf and root tissues under 50 mM NaHCO₃ treatment. In the present study, a total of 7856 phosphopeptides were identified from cultivated soybeans (Glycine max L. Merr.), representing 3468 phosphoprotein groups, in which 2427 phosphoprotein groups were newly identified. These phosphoprotein groups contained 6326 unique high-probability phosphosites (UHPs), of which 77.2% were newly identified, increasing the current soybean phosphosite database size by 43.4%. Among the phosphopeptides found in this study, we determined 67 phosphopeptides (representing 63 phosphoprotein groups) from leaf tissue and 554 phosphopeptides (representing 487 phosphoprotein groups) from root tissue that showed significant changes in phosphorylation levels under sodium bicarbonate stress (fold change >1.2 or <0.83, respectively; p < 0.05). Localization prediction showed that most phosphoproteins localized in the nucleus for both leaf and root tissues. GO and KEGG enrichment analyses showed quite different enriched functional terms between leaf and root tissues, and more pathways were enriched in the root tissue than in the leaf tissue. Moreover, a total of 53 different protein kinases and 7 protein phosphatases were identified from the differentially expressed phosphoproteins (DEPs). A protein kinase/phosphatase interactor analysis showed that the interacting proteins were mainly involved in/with transporters/membrane trafficking, transcriptional level regulation, protein level regulation, signaling/stress response, and miscellaneous functions. The results presented in this study reveal insights into the function of post-translational modification in plant responses to sodium bicarbonate stress.