AbstractPituitary neuroendocrine tumor (PitNET) is one of the most common intracranial tumors. Due to its extensive tumor heterogeneity and the lack of high-quality tissues for biomarker discovery, the causative molecular mechanisms are far from being fully defined. Therefore, more studies are needed to improve the current clinicopathological classification system, and advanced treatment strategies such as targeted therapy and immunotherapy are yet to be explored. Here, we performed the largest integrative genomics, transcriptomics, proteomics, and phosphoproteomics analysis reported to date for a cohort of 200 PitNET patients. Genomics data indicate that GNAS copy number gain can serve as a reliable diagnostic marker for hyperproliferation of the PIT1 lineage. Proteomics-based classification of PitNETs identified 7 clusters, among which, tumors overexpressing epithelial-mesenchymal transition (EMT) markers clustered into a more invasive subgroup. Further analysis identified potential therapeutic targets, including CDK6, TWIST1, EGFR, and VEGFR2, for different clusters. Immune subtyping to explore the potential for application of immunotherapy in PitNET identified an association between alterations in the JAK1-STAT1-PDL1 axis and immune exhaustion, and between changes in the JAK3-STAT6-FOS/JUN axis and immune infiltration. These identified molecular markers and alternations in various clusters/subtypes were further confirmed in an independent cohort of 750 PitNET patients. This proteogenomic analysis across traditional histological boundaries improves our current understanding of PitNET pathophysiology and suggests novel therapeutic targets and strategies.