Proteogenomic characterization of skull-base chordoma.

Skull-base chordoma is a rare, aggressive bone cancer with a high recurrence rate. Despite advances in genomic studies, its molecular characteristics and effective therapies remain unknown. Here, we conduct integrative genomics, transcriptomics, proteomics, and phosphoproteomics analyses of 187 skull-base chordoma tumors. In our study, chromosome instability is identified as a prognostic predictor and potential therapeutic target. Multi-omics data reveals downstream effects of chromosome instability, with RPRD1B as a putative target for radiotherapy-resistant patients. Chromosome 1q gain, associated with chromosome instability and upregulated mitochondrial functions, lead to poorer clinical outcomes. Immune subtyping identify an immune cold subtype linked to chromosome 9p/10q loss and immune evasion. Proteomics-based classification reveals subtypes (P-II and P-III) with high chromosome instability and immune cold features, with P-II tumors showing increased invasiveness. These findings, confirmed in 17 paired samples, provide insights into the biology and treatment of skull-base chordoma.

Proteomic characterization of esophageal squamous cell carcinoma response to immunotherapy reveals potential therapeutic strategy and predictive biomarkers.

Immunotherapy is the first-line therapy for esophageal squamous cell carcinoma (ESCC), yet many patients do not respond due to drug resistance and the lack of reliable predictive markers. We collected 73 ESCC patients (including discovery cohort and validation cohort) without immune thrombocytopenia and undergoing anti-PD1 immunotherapy. Proteomic and phosphoproteomic analysis of 73 ESCC treatment-naive samples by mass spectrometry-based label-free quantification were applied to explore the potential resistant and sensitive mechanisms, and identify predictive markers of ESCC immunotherapy. Comparative analysis found the pathways related to immune and mitochondrial functions were associated with ESCC immunotherapy sensitivity; while platelet activation bioprocess showed negative correlation with CD8+ T cells and related to ESCC immunotherapy non-sensitivity. Finally, we identified 10 ESCC immunotherapy predictive biomarkers with high accuracy (≥ 0.90) to predict the immunotherapeutic response, which was validated in the independent cohort.

Longitudinal plasma proteome profiling reveals the diversity of biomarkers for diagnosis and cetuximab therapy response of colorectal cancer.

Cetuximab therapy is the major treatment for colorectal cancer (CRC), but drug resistance limits its effectiveness. Here, we perform longitudinal and deep proteomic profiling of 641 plasma samples originated from 147 CRC patients (CRCs) undergoing cetuximab therapy with multi-course treatment, and 90 healthy controls (HCs). COL12A1, THBS2, S100A8, and S100A9 are screened as potential proteins to distinguish CRCs from HCs both in plasma and tissue validation cohorts. We identify the potential biomarkers (RRAS2, MMP8, FBLN1, RPTOR, and IMPDH2) for the initial response prediction. In a longitudinal setting, we identify two clusters with distinct fluctuations and construct the model with high accuracy to predict the longitudinal response, further validated in the independent cohort. This study reveals the heterogeneity of different biomarkers for tumor diagnosis, the initial and longitudinal response prediction respectively in the first course and multi-course cetuximab treatment, may ultimately be useful in monitoring and intervention strategies for CRC.

Proteomic characterization of the colorectal cancer response to chemoradiation and targeted therapies reveals potential therapeutic strategies.

Chemoradiation and targeted therapies are the major treatments for colorectal cancer (CRC); however, molecular properties associated with therapy resistance are incompletely characterized. Here, we profile the proteome of 254 tumor tissues from patients with CRC undergoing chemotherapy, chemoradiation, or chemotherapy combined with targeted therapy. Proteome-based classification reveals four subtypes featured with distinct biological and therapeutic characteristics. The integrative analysis of CRC cell lines and clinical samples indicates that immune regulation is significantly associated with drug sensitivity. HSF1 can increase DNA damage repair and cell cycle, thus inducing resistance to radiation, while high expression of HDAC6 is negatively associated with response of cetuximab. Furthermore, we develop prognostic models with high accuracy to predict the therapeutic response, further validated by parallel reaction monitoring (PRM) assay in an independent validation cohort. This study provides a rich resource for investigating the mechanisms and indicators of chemoradiation and targeted therapy in CRC.

Integrative proteomics and metabolomics study reveal enhanced immune responses by COVID-19 vaccine booster shot against Omicron SARS-CoV-2 infection.

Since its outbreak in late 2021, the Omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been widely reported to be able to evade neutralizing antibodies, becoming more transmissible while causing milder symptoms than previous SARS-CoV-2 strains. Understanding the underlying molecular changes of Omicron SARS-CoV-2 infection and corresponding host responses are important to the control of Omicron COVID-19 pandemic. In this study, we report an integrative proteomics and metabolomics investigation of serum samples from 80 COVID-19 patients infected with Omicron SARS-CoV-2, as well as 160 control serum samples from 80 healthy individuals and 80 patients who had flu-like symptoms but were negative for SARS-CoV-2 infection. The multiomics results indicated that Omicron SARS-CoV-2 infection caused significant changes to host serum proteome and metabolome comparing to the healthy controls and patients who had flu-like symptoms without COVID-19. Protein and metabolite changes also pointed to liver dysfunctions and potential damage to other host organs by Omicron SARS-CoV-2 infection. The Omicron COVID-19 patients could be roughly divided into two subgroups based on their proteome differences. Interestingly, the subgroup who mostly had received full vaccination with booster shot had fewer coughing symptom, changed sphingomyelin lipid metabolism, and stronger immune responses including higher numbers of lymphocytes, monocytes, neutrophils, and upregulated proteins related to CD4+ T cells, CD8+ effector memory T cells (Tem), and conventional dendritic cells, revealing beneficial effects of full COVID-19 vaccination against Omicron SARS-CoV-2 infection through molecular changes.

Quartet protein reference materials and datasets for multi-platform assessment of label-free proteomics.

BACKGROUND: Quantitative proteomics is an indispensable tool in life science research. However, there is a lack of reference materials for evaluating the reproducibility of label-free liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based measurements among different instruments and laboratories. RESULTS: Here, we develop the Quartet standard as a proteome reference material with built-in truths, and distribute the same aliquots to 15 laboratories with nine conventional LC-MS/MS platforms across six cities in China. Relative abundance of over 12,000 proteins on 816 mass spectrometry files are obtained and compared for reproducibility among the instruments and laboratories to ultimately generate proteomics benchmark datasets. There is a wide dynamic range of proteomes spanning about 7 orders of magnitude, and the injection order has marked effects on quantitative instead of qualitative characteristics. CONCLUSION: Overall, the Quartet offers valuable standard materials and data resources for improving the quality control of proteomic analyses as well as the reproducibility and reliability of research findings.

Proteogenomics of different urothelial bladder cancer stages reveals distinct molecular features for papillary cancer and carcinoma in situ.

The progression of urothelial bladder cancer (UC) is a complicated multi-step process. We perform a comprehensive multi-omics analysis of 448 samples from 190 UC patients, covering the whole spectrum of disease stages and grades. Proteogenomic integration analysis indicates the mutations of HRAS regulated mTOR signaling to form urothelial papilloma rather than papillary urothelial cancer (PUC). DNA damage is a key signaling pathway in the progression of carcinoma in situ (CIS) and related to APOBEC signature. Glucolipid metabolism increase and lower immune cell infiltration are associated with PUC compared to CIS. Proteomic analysis distinguishes the origins of invasive tumors (PUC-derived and CIS-derived), related to distinct clinical prognosis and molecular features. Additionally, loss of RBPMS, associated with CIS-derived tumors, is validated to increase the activity of AP-1 and promote metastasis. This study reveals the characteristics of two distinct branches (PUC and CIS) of UC progression and may eventually benefit clinical practice.

Plasma proteomic profiling discovers molecular features associated with upper tract urothelial carcinoma.

Upper tract urothelial carcinoma (UTUC) is often diagnosed late and exhibits poor prognosis. Limited data are available on potential non-invasive biomarkers for disease monitoring. Here, we investigate the proteomic profile of plasma in 362 UTUC patients and 239 healthy controls. We present an integrated tissue-plasma proteomic approach to infer the signature proteins for identifying patients with muscle-invasive UTUC. We discover a protein panel that reflects lymph node metastasis, which is of interest in identifying UTUC patients with high risk and poor prognosis. We also identify a ten-protein classifier and establish a progression clock predicting progression-free survival of UTUC patients. Finally, we further validate the signature proteins by parallel reaction monitoring assay in an independent cohort. Collectively, this study portrays the plasma proteomic landscape of a UTUC cohort and provides a valuable resource for further biological and diagnostic research in UTUC.

Modulation of cellular metabolism by protein crotonylation regulates pancreatic cancer progression.

Protein lysine crotonylation has been recently identified as a vital posttranslational modification in cellular processes, particularly through the modification of histones. We show that lysine crotonylation is an important modification of the cytoplastic and mitochondria proteins. Enzymes in glycolysis, the tricarboxylic acid (TCA) cycle, fatty acid metabolism, glutamine metabolism, glutathione metabolism, the urea cycle, one-carbon metabolism, and mitochondrial fusion/fission dynamics are found to be extensively crotonylated in pancreatic cancer cells. This modulation is mainly controlled by a pair of crotonylation writers and erasers including CBP/p300, HDAC1, and HDAC3. The dynamic crotonylation of metabolic enzymes is involved in metabolism regulation, which is linked with tumor progression. Interestingly, the activation of MTHFD1 by decrotonylation at Lys354 and Lys553 promotes the development of pancreatic cancer by increasing resistance to ferroptosis. Our study suggests that crotonylation represents a metabolic regulatory mechanism in pancreatic cancer progression.

Proteomic characteristics reveal the signatures and the risks of T1 colorectal cancer metastasis to lymph nodes.

The presence of lymph node metastasis (LNM) affects treatment strategy decisions in T1NxM0 colorectal cancer (CRC), but the currently used clinicopathological-based risk stratification cannot predict LNM accurately. In this study, we detected proteins in formalin-fixed paraffin-embedded (FFPE) tumor samples from 143 LNM-negative and 78 LNM-positive patients with T1 CRC and revealed changes in molecular and biological pathways by label-free liquid chromatography tandem mass spectrometry (LC-MS/MS) and established classifiers for predicting LNM in T1 CRC. An effective 55-proteins prediction model was built by machine learning and validated in a training cohort (N=132) and two validation cohorts (VC1, N=42; VC2, N=47), achieved an impressive AUC of 1.00 in the training cohort, 0.96 in VC1 and 0.93 in VC2, respectively. We further built a simplified classifier with nine proteins, and achieved an AUC of 0.824. The simplified classifier was performed excellently in two external validation cohorts. The expression patterns of 13 proteins were confirmed by immunohistochemistry, and the IHC score of five proteins was used to build an IHC predict model with an AUC of 0.825. RHOT2 silence significantly enhanced migration and invasion of colon cancer cells. Our study explored the mechanism of metastasis in T1 CRC and can be used to facilitate the individualized prediction of LNM in patients with T1 CRC, which may provide a guidance for clinical practice in T1 CRC.

Most patients with early-stage colorectal cancer can be treated with a minimally invasive procedure. Surgeons use a flexible tool to remove precancerous or cancerous cells, cutting the risk of death from colorectal cancer in half. But a small number of early-stage colorectal cancer patients are at risk of their cancer spreading to the lymph nodes. These patients need more extensive surgery. Clinicians use risk stratification tools to decide which patients need more extensive surgery. Unfortunately, the existing risk stratification tools are not very accurate. The current approach, which analyzes colon tissue for cancerous changes, classifies 70% to 80% of early-stage colorectal cancer patients as high risk for cancer spread. But only about 8% to 16% of patients in the high risk group have lymph node metastasis. As a result, many patients undergo unnecessary, invasive surgery. Zhuang, Zhuang, Chen, Qin, et al. developed a more accurate way to predict which patients are at risk of lymph node metastasis using proteins. In the experiments, the team analyzed the proteins in tumor samples from 143 patients with early colorectal cancer who did not have lymph node metastases and 78 patients with metastases. Zhuang et al. then used machine learning to build a prediction tool that used 55 proteins to identify patients at risk of metastases. The new approach was more accurate than existing tools and simplified versions with only nine or five proteins also performed better than existing tools. This work provides preliminary evidence that protein-based models using as few as five proteins can more accurately identify which patients are at risk of metastasis. These models may reduce the number of patients who undergo unnecessary invasive surgery. The experiments also identified potential targets for therapies to prevent or treat lymph metastases. For example, they showed that low levels of the RHOT2 protein predict metastasis.

Integrative proteogenomic characterization of early esophageal cancer.

Esophageal squamous cell carcinoma (ESCC) is malignant while the carcinogenesis is still unclear. Here, we perform a comprehensive multi-omics analysis of 786 trace-tumor-samples from 154 ESCC patients, covering 9 histopathological stages and 3 phases. Proteogenomics elucidates cancer-driving waves in ESCC progression, and reveals the molecular characterization of alcohol drinking habit associated signatures. We discover chromosome 3q gain functions in the transmit from nontumor to intraepithelial neoplasia phases, and find TP53 mutation enhances DNA replication in intraepithelial neoplasia phase. The mutations of AKAP9 and MCAF1 upregulate glycolysis and Wnt signaling, respectively, in advanced-stage ESCC phase. Six major tracks related to different clinical features during ESCC progression are identified, which is validated by an independent cohort with another 256 samples. Hyperphosphorylated phosphoglycerate kinase 1 (PGK1, S203) is considered as a drug target in ESCC progression. This study provides insight into the understanding of ESCC molecular mechanism and the development of therapeutic targets.

Comprehensive proteogenomic characterization of early duodenal cancer reveals the carcinogenesis tracks of different subtypes.

The subtypes of duodenal cancer (DC) are complicated and the carcinogenesis process is not well characterized. We present comprehensive characterization of 438 samples from 156 DC patients, covering 2 major and 5 rare subtypes. Proteogenomics reveals LYN amplification at the chromosome 8q gain functioned in the transmit from intraepithelial neoplasia phase to infiltration tumor phase via MAPK signaling, and illustrates the DST mutation improves mTOR signaling in the duodenal adenocarcinoma stage. Proteome-based analysis elucidates stage-specific molecular characterizations and carcinogenesis tracks, and defines the cancer-driving waves of the adenocarcinoma and Brunner's gland subtypes. The drug-targetable alanyl-tRNA synthetase (AARS1) in the high tumor mutation burden/immune infiltration is significantly enhanced in DC progression, and catalyzes the lysine-alanylation of poly-ADP-ribose polymerases (PARP1), which decreases the apoptosis of cancer cells, eventually promoting cell proliferation and tumorigenesis. We assess the proteogenomic landscape of early DC, and provide insights into the molecular features corresponding therapeutic targets.

Proteogenomic characterization of cholangiocarcinoma.

BACKGROUND AND AIMS: Cholangiocarcinoma (CCA) is a highly heterogeneous cancer with limited understanding and few effective therapeutic approaches. We aimed at providing a proteogenomic CCA characterization to inform biological processes and treatment vulnerabilities. APPROACH AND RESULTS: Integrative genomic analysis with functional validation uncovered biological perturbations downstream of driver events including DPCR1 , RBM47 mutations, SH3BGRL2 copy number alterations, and FGFR2 fusions in CCA. Proteomic clustering identified three subtypes with distinct clinical outcomes, molecular features, and potential therapeutics. Phosphoproteomics characterized targetable kinases in CCA, suggesting strategies for effective treatment with CDK and MAPK inhibitors. Patients with CCA with HBV infection showed increased antigen processing and presentation (APC) and T cell infiltration, conferring a favorable prognosis compared with those without HBV infection. The characterization of extrahepatic CCA recommended the feasible application of vascular endothelial-derived growth factor inhibitors. Multiomics profiling presented distinctive molecular characteristics of the large bile duct and the small bile duct of intrahepatic CCA. The immune landscape further revealed diverse tumor immune microenvironments, suggesting immune subtypes C1 and C5 might benefit from immune checkpoint therapy. TCN1 was identified as a potential CCA prognostic biomarker, promoting cell growth by enhancing vitamin B12 metabolism. CONCLUSIONS: We characterized the proteogenomic landscape of 217 CCAs with 197 paired normal adjacent tissues and identified their subtypes and potential therapeutic targets. The multiomics analyses with other databases and some functional validations have indicated strategies regarding the clinical, biological, and therapeutic approaches to the management of CCA.

GwAAP: A genome-wide amino acid coding-decoding quantitative proteomics system.

Mass spectrometry-based proteomic technology has greatly improved and has been widely applied in various biological science fields. However, proteome-wide accurate quantification of proteins in signaling pathways remains challenging. Here, we report a genome-wide amino acid coding-decoding quantitative proteomic (GwAAP) system to facilitate precise proteome quantification. For each protein, a unique code peptide was assigned and incorporated into the N-terminus of the targeted protein and used for identification and quantification. As a proof of principle, we systematically tagged 40 yeast proteins with codes and employed mass spectrometry to decode. We successfully recovered all 40 code peptides with a large and consistent quantitative dynamic range (CV slope <10%, R2 > 0.8). We further verified the alteration of the glucose and galactose metabolism pathways in yeast under different carbon source conditions. The GwAAP system could potentially provide a strategy to achieve absolute quantification of the entire yeast proteome without bias.

Integrated proteomic and transcriptomic landscape of macrophages in mouse tissues.

Macrophages are involved in tissue homeostasis and are critical for innate immune responses, yet distinct macrophage populations in different tissues exhibit diverse gene expression patterns and biological processes. While tissue-specific macrophage epigenomic and transcriptomic profiles have been reported, proteomes of different macrophage populations remain poorly characterized. Here we use mass spectrometry and bulk RNA sequencing to assess the proteomic and transcriptomic patterns, respectively, of 10 primary macrophage populations from seven mouse tissues, bone marrow-derived macrophages and the cell line RAW264.7. The results show distinct proteomic landscape and protein copy numbers between tissue-resident and recruited macrophages. Construction of a hierarchical regulatory network finds cell-type-specific transcription factors of macrophages serving as hubs for denoting tissue and functional identity of individual macrophage subsets. Finally, Il18 is validated to be essential in distinguishing molecular signatures and cellular function features between tissue-resident and recruited macrophages in the lung and liver. In summary, these deposited datasets and our open proteome server ( http://macrophage.mouseprotein.cn ) integrating all information will provide a valuable resource for future functional and mechanistic studies of mouse macrophages.

Proteogenomic insights into the biology and treatment of pancreatic ductal adenocarcinoma.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is a devastating disease with poor prognosis. Proteogenomic characterization and integrative proteomic analysis provide a functional context to annotate genomic abnormalities with prognostic value. METHODS: We performed an integrated multi-omics analysis, including whole-exome sequencing, RNA-seq, proteomic, and phosphoproteomic analysis of 217 PDAC tumors with paired non-tumor adjacent tissues. In vivo functional experiments were performed to further illustrate the biological events related to PDAC tumorigenesis and progression. RESULTS: A comprehensive proteogenomic landscape revealed that TP53 mutations upregulated the CDK4-mediated cell proliferation process and led to poor prognosis in younger patients. Integrative multi-omics analysis illustrated the proteomic and phosphoproteomic alteration led by genomic alterations such as KRAS mutations and ADAM9 amplification of PDAC tumorigenesis. Proteogenomic analysis combined with in vivo experiments revealed that the higher amplification frequency of ADAM9 (8p11.22) could drive PDAC metastasis, though downregulating adhesion junction and upregulating WNT signaling pathway. Proteome-based stratification of PDAC revealed three subtypes (S-I, S-II, and S-III) related to different clinical and molecular features. Immune clustering defined a metabolic tumor subset that harbored FH amplicons led to better prognosis. Functional experiments revealed the role of FH in altering tumor glycolysis and in impacting PDAC tumor microenvironments. Experiments utilizing both in vivo and in vitro assay proved that loss of HOGA1 promoted the tumor growth via activating LARP7-CDK1 pathway. CONCLUSIONS: This proteogenomic dataset provided a valuable resource for researchers and clinicians seeking for better understanding and treatment of PDAC.

Integrated proteogenomic characterization across major histological types of pituitary neuroendocrine tumors.

Pituitary 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.

Proteomic characterization of gastric cancer response to chemotherapy and targeted therapy reveals new therapeutic strategies.

Chemotherapy and targeted therapy are the major treatments for gastric cancer (GC), but drug resistance limits its effectiveness. Here, we profile the proteome of 206 tumor tissues from patients with GC undergoing either chemotherapy or anti-HER2-based therapy. Proteome-based classification reveals four subtypes (G-I-G-IV) related to different clinical and molecular features. MSI-sig high GC patients benefit from docetaxel combination treatment, accompanied by anticancer immune response. Further study reveals patients with high T cell receptor signaling respond to anti-HER2-based therapy; while activation of extracellular matrix/PI3K-AKT pathway impair anti-tumor effect of trastuzumab. We observe CTSE functions as a cell intrinsic enhancer of chemosensitivity of docetaxel, whereas TKTL1 functions as an attenuator. Finally, we develop prognostic models with high accuracy to predict therapeutic response, further validated in an independent validation cohort. This study provides a rich resource for investigating the mechanisms and indicators of chemotherapy and targeted therapy in GC.

Proteomic analysis reveals key differences between squamous cell carcinomas and adenocarcinomas across multiple tissues.

Squamous cell carcinoma (SCC) and adenocarcinoma (AC) are two main histological subtypes of solid cancer; however, SCCs are derived from different organs with similar morphologies, and it is challenging to distinguish the origin of metastatic SCCs. Here we report a deep proteomic analysis of 333 SCCs of 17 organs and 69 ACs of 7 organs. Proteomic comparison between SCCs and ACs identifies distinguishable pivotal pathways and molecules in those pathways play consistent adverse or opposite prognostic roles in ACs and SCCs. A comparison between common and rare SCCs highlights lipid metabolism may reinforce the malignancy of rare SCCs. Proteomic clusters reveal anatomical features, and kinase-transcription factor networks indicate differential SCC characteristics, while immune subtyping reveals diverse tumor microenvironments across and within diagnoses and identified potential druggable targets. Furthermore, tumor-specific proteins provide candidates with differentially diagnostic values. This proteomics architecture represents a public resource for researchers seeking a better understanding of SCCs and ACs.