Abstract

Background 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. Methods The 6-plex tandem mass tag (TMT) labeling reagents coupled with TiO 2 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. Results 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