Integration of quantitative phosphoproteomics and transcriptomics revealed phosphorylation-mediated molecular events as useful tools for a potential patient stratification and personalized treatment of human nonfunctional pituitary adenomas.

Abstract

Invasiveness is a very challenging clinical problem in nonfunctional pituitary adenomas (NFPAs), and currently, there are no effective invasiveness-related molecular biomarkers. The post-neurosurgery treatment is much different as for invasive and noninvasive NFPAs. The aim of this study was to integrate phosphoproteomics and transcriptomics data to reveal phosphorylation-mediated molecular events for invasive characteristics of NFPAs to achieve a potential tool for patient stratification, and prognostic/predictive assessment to discriminate invasive from noninvasive NFPAs for personalized attitude. The 6-plex tandem mass tag (TMT) labeling reagents coupled with TiO2 enrichment of phosphopeptides and liquid chromatography-tandem mass spectrometry (LC-MS/MS) were used to identify and quantify each phosphoprotein and phosphosite in NFPAs and controls. Differentially expressed genes (DEGs) between invasive NFPA and control tissues were obtained from the Gene Expression Omnibus (GEO) database. The overlapping analysis was performed between phosphoprotiens and invasive DEGs. Gene Ontology (GO) enrichment, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and protein-protein interaction (PPI) analyses were used to analyze these overlapped molecules. In total, 1035 phosphoproteins with 2982 phosphorylation sites were identified in NFPAs vs. controls, and 2751 DEGs were identified in invasive NFPAs vs. controls. Overlapping analysis of these phosphoproteins and DEGs exposed 130 overlapped molecules (phosphoproteins; invasive DEGs). GO enrichment and KEGG pathway analyses of 130 overlapped molecules revealed multiple biological processes and signaling pathway network alterations, including cell-cell adhesion, platelet activation, GTPase signaling pathway, protein kinase signaling, calcium signaling pathway, estrogen signaling pathway, glucagon signaling pathway, cGMP-PKG signaling pathway, GnRH signaling pathway, inflammatory mediator regulation of TRP channels, vascular smooth muscle contraction, and Fc gamma R-mediated phagocytosis, which were obviously associated with tumor invasive characteristics. For 130 overlapped molecules, PPI network-based molecular complex detection (MCODE) identified 10 hub molecules, namely SLC2A4, TSC2, AKT1, SCG3, ALB, APOL1, ACACA, SPARCL1, CHGB, and IGFBP5. These hub molecules are involved in multiple signaling pathways and represent potential predictive/prognostic markers in NFPA patients as well as they represent potential therapeutic targets. This study provided the first large-scale phosphoprotein profiling and phosphorylation-related signaling pathway network alterations in human NFPA tissues. Further, overlapping analysis of phosphoproteins and invasive DEGs revealed the phosphorylation-mediated signaling pathway network changes in invasive NFPAs. These findings are the precious resource for in-depth insight into the molecular mechanisms of NFPAs, as well as for the discovery of effective phosphoprotein biomarkers and therapeutic targets for invasive NFPAs.