Global, site-resolved analysis of ubiquitylation occupancy and turnover rate reveals systems properties.

Ubiquitylation regulates most proteins and biological processes in a eukaryotic cell. However, the site-specific occupancy (stoichiometry) and turnover rate of ubiquitylation have not been quantified. Here we present an integrated picture of the global ubiquitylation site occupancy and half-life. Ubiquitylation site occupancy spans over four orders of magnitude, but the median ubiquitylation site occupancy is three orders of magnitude lower than that of phosphorylation. The occupancy, turnover rate, and regulation of sites by proteasome inhibitors are strongly interrelated, and these attributes distinguish sites involved in proteasomal degradation and cellular signaling. Sites in structured protein regions exhibit longer half-lives and stronger upregulation by proteasome inhibitors than sites in unstructured regions. Importantly, we discovered a surveillance mechanism that rapidly and site-indiscriminately deubiquitylates all ubiquitin-specific E1 and E2 enzymes, protecting them against accumulation of bystander ubiquitylation. The work provides a systems-scale, quantitative view of ubiquitylation properties and reveals general principles of ubiquitylation-dependent governance.

Pan-cancer proteogenomics characterization of tumor immunity.

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan-cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype-specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

Pan-cancer proteogenomics connects oncogenic drivers to functional states.

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation.

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Workflow enabling deepscale immunopeptidome, proteome, ubiquitylome, phosphoproteome, and acetylome analyses of sample-limited tissues.

Serial multi-omic analysis of proteome, phosphoproteome, and acetylome provides insights into changes in protein expression, cell signaling, cross-talk and epigenetic pathways involved in disease pathology and treatment. However, ubiquitylome and HLA peptidome data collection used to understand protein degradation and antigen presentation have not together been serialized, and instead require separate samples for parallel processing using distinct protocols. Here we present MONTE, a highly sensitive multi-omic native tissue enrichment workflow, that enables serial, deep-scale analysis of HLA-I and HLA-II immunopeptidome, ubiquitylome, proteome, phosphoproteome, and acetylome from the same tissue sample. We demonstrate that the depth of coverage and quantitative precision of each 'ome is not compromised by serialization, and the addition of HLA immunopeptidomics enables the identification of peptides derived from cancer/testis antigens and patient specific neoantigens. We evaluate the technical feasibility of the MONTE workflow using a small cohort of patient lung adenocarcinoma tumors.

Proteogenomic Markers of Chemotherapy Resistance and Response in Triple-Negative Breast Cancer.

Microscaled proteogenomics was deployed to probe the molecular basis for differential response to neoadjuvant carboplatin and docetaxel combination chemotherapy for triple-negative breast cancer (TNBC). Proteomic analyses of pretreatment patient biopsies uniquely revealed metabolic pathways, including oxidative phosphorylation, adipogenesis, and fatty acid metabolism, that were associated with resistance. Both proteomics and transcriptomics revealed that sensitivity was marked by elevation of DNA repair, E2F targets, G2-M checkpoint, interferon-gamma signaling, and immune-checkpoint components. Proteogenomic analyses of somatic copy-number aberrations identified a resistance-associated 19q13.31-33 deletion where LIG1, POLD1, and XRCC1 are located. In orthogonal datasets, LIG1 (DNA ligase I) gene deletion and/or low mRNA expression levels were associated with lack of pathologic complete response, higher chromosomal instability index (CIN), and poor prognosis in TNBC, as well as carboplatin-selective resistance in TNBC preclinical models. Hemizygous loss of LIG1 was also associated with higher CIN and poor prognosis in other cancer types, demonstrating broader clinical implications. SIGNIFICANCE: Proteogenomic analysis of triple-negative breast tumors revealed a complex landscape of chemotherapy response associations, including a 19q13.31-33 somatic deletion encoding genes serving lagging-strand DNA synthesis (LIG1, POLD1, and XRCC1), that correlate with lack of pathologic response, carboplatin-selective resistance, and, in pan-cancer studies, poor prognosis and CIN. This article is highlighted in the In This Issue feature, p. 2483.

Cancer proteogenomics: current impact and future prospects.

Genomic analyses in cancer have been enormously impactful, leading to the identification of driver mutations and development of targeted therapies. But the functions of the vast majority of somatic mutations and copy number variants in tumours remain unknown, and the causes of resistance to targeted therapies and methods to overcome them are poorly defined. Recent improvements in mass spectrometry-based proteomics now enable direct examination of the consequences of genomic aberrations, providing deep and quantitative characterization of tumour tissues. Integration of proteins and their post-translational modifications with genomic, epigenomic and transcriptomic data constitutes the new field of proteogenomics, and is already leading to new biological and diagnostic knowledge with the potential to improve our understanding of malignant transformation and therapeutic outcomes. In this Review we describe recent developments in proteogenomics and key findings from the proteogenomic analysis of a wide range of cancers. Considerations relevant to the selection and use of samples for proteogenomics and the current technologies used to generate, analyse and integrate proteomic with genomic data are described. Applications of proteogenomics in translational studies and immuno-oncology are rapidly emerging, and the prospect for their full integration into therapeutic trials and clinical care seems bright.

A highly multiplexed quantitative phosphosite assay for biology and preclinical studies.

Reliable methods to quantify dynamic signaling changes across diverse pathways are needed to better understand the effects of disease and drug treatment in cells and tissues but are presently lacking. Here, we present SigPath, a targeted mass spectrometry (MS) assay that measures 284 phosphosites in 200 phosphoproteins of biological interest. SigPath probes a broad swath of signaling biology with high throughput and quantitative precision. We applied the assay to investigate changes in phospho-signaling in drug-treated cancer cell lines, breast cancer preclinical models, and human medulloblastoma tumors. In addition to validating previous findings, SigPath detected and quantified a large number of differentially regulated phosphosites newly associated with disease models and human tumors at baseline or with drug perturbation. Our results highlight the potential of SigPath to monitor phosphoproteomic signaling events and to nominate mechanistic hypotheses regarding oncogenesis, response, and resistance to therapy.

Genomic Profiling of Lung Adenocarcinoma in Never-Smokers.

PURPOSE: Approximately 10%-40% of patients with lung cancer report no history of tobacco smoking (never-smokers). We analyzed whole-exome and RNA-sequencing data of 160 tumor and normal lung adenocarcinoma (LUAD) samples from never-smokers to identify clinically actionable alterations and gain insight into the environmental and hereditary risk factors for LUAD among never-smokers. METHODS: We performed whole-exome and RNA-sequencing of 88 and 69 never-smoker LUADs. We analyzed these data in conjunction with data from 76 never-smoker and 299 smoker LUAD samples sequenced by The Cancer Genome Atlas and Clinical Proteomic Tumor Analysis Consortium. RESULTS: We observed a high prevalence of clinically actionable driver alterations in never-smoker LUADs compared with smoker LUADs (78%-92% v 49.5%; P < .0001). Although a subset of never-smoker samples demonstrated germline alterations in DNA repair genes, the frequency of samples showing germline variants in cancer predisposing genes was comparable between smokers and never-smokers (6.4% v 6.9%; P = .82). A subset of never-smoker samples (5.9%) showed mutation signatures that were suggestive of passive exposure to cigarette smoke. Finally, analysis of RNA-sequencing data showed distinct immune transcriptional subtypes of never-smoker LUADs that varied in their expression of clinically relevant immune checkpoint molecules and immune cell composition. CONCLUSION: In this comprehensive genomic and transcriptome analysis of never-smoker LUADs, we observed a potential role for germline variants in DNA repair genes and passive exposure to cigarette smoke in the pathogenesis of a subset of never-smoker LUADs. Our findings also show that clinically actionable driver alterations are highly prevalent in never-smoker LUADs, highlighting the need for obtaining biopsies with adequate cellularity for clinical genomic testing in these patients.

Automating UbiFast for High-throughput and Multiplexed Ubiquitin Enrichment.

Robust methods for deep-scale enrichment and site-specific identification of ubiquitylation sites are necessary for characterizing the myriad roles of protein ubiquitylation. To this end we previously developed UbiFast, a sensitive method for highly multiplexed ubiquitylation profiling where K-ϵ-GG peptides are enriched with anti-K-ε-GG antibody and labeled on-antibody with isobaric labeling reagents for sample multiplexing. Here, we present robotic automation of the UbiFast method using a magnetic bead-conjugated K-ε-GG antibody (mK-ε-GG) and a magnetic particle processor. We report the identification of ∼20,000 ubiquitylation sites from a TMT10-plex with 500 µg input per sample processed in ∼2 h. Automation of the UbiFast method greatly increased the number of identified and quantified ubiquitylation sites, improved reproducibility, and significantly reduced processing time. The automated method also significantly reduced variability across process replicates compared with the manual method. The workflow enables processing of up to 96 samples in a single day making it suitable to study ubiquitylation in large sample sets. Here we demonstrate the applicability of the method to profile small amounts of tissue using breast cancer patient-derived xenograft (PDX) tissue samples.

A proteogenomic portrait of lung squamous cell carcinoma.

Lung squamous cell carcinoma (LSCC) remains a leading cause of cancer death with few therapeutic options. We characterized the proteogenomic landscape of LSCC, providing a deeper exposition of LSCC biology with potential therapeutic implications. We identify NSD3 as an alternative driver in FGFR1-amplified tumors and low-p63 tumors overexpressing the therapeutic target survivin. SOX2 is considered undruggable, but our analyses provide rationale for exploring chromatin modifiers such as LSD1 and EZH2 to target SOX2-overexpressing tumors. Our data support complex regulation of metabolic pathways by crosstalk between post-translational modifications including ubiquitylation. Numerous immune-related proteogenomic observations suggest directions for further investigation. Proteogenomic dissection of CDKN2A mutations argue for more nuanced assessment of RB1 protein expression and phosphorylation before declaring CDK4/6 inhibition unsuccessful. Finally, triangulation between LSCC, LUAD, and HNSCC identified both unique and common therapeutic vulnerabilities. These observations and proteogenomics data resources may guide research into the biology and treatment of LSCC.

SIK2 orchestrates actin-dependent host response upon Salmonella infection.

Salmonella is an intracellular pathogen of a substantial global health concern. In order to identify key players involved in Salmonella infection, we performed a global host phosphoproteome analysis subsequent to bacterial infection. Thereby, we identified the kinase SIK2 as a central component of the host defense machinery upon Salmonella infection. SIK2 depletion favors the escape of bacteria from the Salmonella-containing vacuole (SCV) and impairs Xenophagy, resulting in a hyperproliferative phenotype. Mechanistically, SIK2 associates with actin filaments under basal conditions; however, during bacterial infection, SIK2 is recruited to the SCV together with the elements of the actin polymerization machinery (Arp2/3 complex and Formins). Notably, SIK2 depletion results in a severe pathological cellular actin nucleation and polymerization defect upon Salmonella infection. We propose that SIK2 controls the formation of a protective SCV actin shield shortly after invasion and orchestrates the actin cytoskeleton architecture in its entirety to control an acute Salmonella infection after bacterial invasion.

Enhancers are activated by p300/CBP activity-dependent PIC assembly, RNAPII recruitment, and pause release.

The metazoan-specific acetyltransferase p300/CBP is involved in activating signal-induced, enhancer-mediated transcription of cell-type-specific genes. However, the global kinetics and mechanisms of p300/CBP activity-dependent transcription activation remain poorly understood. We performed genome-wide, time-resolved analyses to show that enhancers and super-enhancers are dynamically activated through p300/CBP-catalyzed acetylation, deactivated by the opposing deacetylase activity, and kinetic acetylation directly contributes to maintaining cell identity at very rapid (minutes) timescales. The acetyltransferase activity is dispensable for the recruitment of p300/CBP and transcription factors but essential for promoting the recruitment of TFIID and RNAPII at virtually all enhancers and enhancer-regulated genes. This identifies pre-initiation complex assembly as a dynamically controlled step in the transcription cycle and reveals p300/CBP-catalyzed acetylation as the signal that specifically promotes transcription initiation at enhancer-regulated genes. We propose that p300/CBP activity uses a "recruit-and-release" mechanism to simultaneously promote RNAPII recruitment and pause release and thereby enables kinetic activation of enhancer-mediated transcription.

Proteome dynamics at broken replication forks reveal a distinct ATM-directed repair response suppressing DNA double-strand break ubiquitination.

Cells have evolved an elaborate DNA repair network to ensure complete and accurate DNA replication. Defects in these repair machineries can fuel genome instability and drive carcinogenesis while creating vulnerabilities that may be exploited in therapy. Here, we use nascent chromatin capture (NCC) proteomics to characterize the repair of replication-associated DNA double-strand breaks (DSBs) triggered by topoisomerase 1 (TOP1) inhibitors. We reveal profound changes in the fork proteome, including the chromatin environment and nuclear membrane interactions, and identify three classes of repair factors according to their enrichment at broken and/or stalled forks. ATM inhibition dramatically rewired the broken fork proteome, revealing that ataxia telangiectasia mutated (ATM) signalling stimulates DNA end resection, recruits PLK1, and concomitantly suppresses the canonical DSB ubiquitination response by preventing accumulation of RNF168 and BRCA1-A. This work and collection of replication fork proteomes provide a new framework to understand how cells orchestrate homologous recombination repair of replication-associated DSBs.

Proteogenomic Landscape of Breast Cancer Tumorigenesis and Targeted Therapy.

The integration of mass spectrometry-based proteomics with next-generation DNA and RNA sequencing profiles tumors more comprehensively. Here this "proteogenomics" approach was applied to 122 treatment-naive primary breast cancers accrued to preserve post-translational modifications, including protein phosphorylation and acetylation. Proteogenomics challenged standard breast cancer diagnoses, provided detailed analysis of the ERBB2 amplicon, defined tumor subsets that could benefit from immune checkpoint therapy, and allowed more accurate assessment of Rb status for prediction of CDK4/6 inhibitor responsiveness. Phosphoproteomics profiles uncovered novel associations between tumor suppressor loss and targetable kinases. Acetylproteome analysis highlighted acetylation on key nuclear proteins involved in the DNA damage response and revealed cross-talk between cytoplasmic and mitochondrial acetylation and metabolism. Our results underscore the potential of proteogenomics for clinical investigation of breast cancer through more accurate annotation of targetable pathways and biological features of this remarkably heterogeneous malignancy.

Proteogenomic Characterization Reveals Therapeutic Vulnerabilities in Lung Adenocarcinoma.

To explore the biology of lung adenocarcinoma (LUAD) and identify new therapeutic opportunities, we performed comprehensive proteogenomic characterization of 110 tumors and 101 matched normal adjacent tissues (NATs) incorporating genomics, epigenomics, deep-scale proteomics, phosphoproteomics, and acetylproteomics. Multi-omics clustering revealed four subgroups defined by key driver mutations, country, and gender. Proteomic and phosphoproteomic data illuminated biology downstream of copy number aberrations, somatic mutations, and fusions and identified therapeutic vulnerabilities associated with driver events involving KRAS, EGFR, and ALK. Immune subtyping revealed a complex landscape, reinforced the association of STK11 with immune-cold behavior, and underscored a potential immunosuppressive role of neutrophil degranulation. Smoking-associated LUADs showed correlation with other environmental exposure signatures and a field effect in NATs. Matched NATs allowed identification of differentially expressed proteins with potential diagnostic and therapeutic utility. This proteogenomics dataset represents a unique public resource for researchers and clinicians seeking to better understand and treat lung adenocarcinomas.

Proteogenomic Characterization of Endometrial Carcinoma.

We undertook a comprehensive proteogenomic characterization of 95 prospectively collected endometrial carcinomas, comprising 83 endometrioid and 12 serous tumors. This analysis revealed possible new consequences of perturbations to the p53 and Wnt/ß-catenin pathways, identified a potential role for circRNAs in the epithelial-mesenchymal transition, and provided new information about proteomic markers of clinical and genomic tumor subgroups, including relationships to known druggable pathways. An extensive genome-wide acetylation survey yielded insights into regulatory mechanisms linking Wnt signaling and histone acetylation. We also characterized aspects of the tumor immune landscape, including immunogenic alterations, neoantigens, common cancer/testis antigens, and the immune microenvironment, all of which can inform immunotherapy decisions. Collectively, our multi-omic analyses provide a valuable resource for researchers and clinicians, identify new molecular associations of potential mechanistic significance in the development of endometrial cancers, and suggest novel approaches for identifying potential therapeutic targets.

Microscaled proteogenomic methods for precision oncology.

Cancer proteogenomics promises new insights into cancer biology and treatment efficacy by integrating genomics, transcriptomics and protein profiling including modifications by mass spectrometry (MS). A critical limitation is sample input requirements that exceed many sources of clinically important material. Here we report a proteogenomics approach for core biopsies using tissue-sparing specimen processing and microscaled proteomics. As a demonstration, we analyze core needle biopsies from ERBB2 positive breast cancers before and 48-72 h after initiating neoadjuvant trastuzumab-based chemotherapy. We show greater suppression of ERBB2 protein and both ERBB2 and mTOR target phosphosite levels in cases associated with pathological complete response, and identify potential causes of treatment resistance including the absence of ERBB2 amplification, insufficient ERBB2 activity for therapeutic sensitivity despite ERBB2 amplification, and candidate resistance mechanisms including androgen receptor signaling, mucin overexpression and an inactive immune microenvironment. The clinical utility and discovery potential of proteogenomics at biopsy-scale warrants further investigation.

Rapid and deep-scale ubiquitylation profiling for biology and translational research.

Protein ubiquitylation is involved in a plethora of cellular processes. While antibodies directed at ubiquitin remnants (K-ɛ-GG) have improved the ability to monitor ubiquitylation using mass spectrometry, methods for highly multiplexed measurement of ubiquitylation in tissues and primary cells using sub-milligram amounts of sample remains a challenge. Here, we present a highly sensitive, rapid and multiplexed protocol termed UbiFast for quantifying ~10,000 ubiquitylation sites from as little as 500 µg peptide per sample from cells or tissue in a TMT10plex in ca. 5 h. High-field Asymmetric Waveform Ion Mobility Spectrometry (FAIMS) is used to improve quantitative accuracy for posttranslational modification analysis. We use the approach to rediscover substrates of the E3 ligase targeting drug lenalidomide and to identify proteins modulated by ubiquitylation in models of basal and luminal human breast cancer. The sensitivity and speed of the UbiFast method makes it suitable for large-scale studies in primary tissue samples.

TMT Labeling for the Masses: A Robust and Cost-efficient, In-solution Labeling Approach.

Isobaric stable isotope labeling using, for example, tandem mass tags (TMTs) is increasingly being applied for large-scale proteomic studies. Experiments focusing on proteoform analysis in drug time course or perturbation studies or in large patient cohorts greatly benefit from the reproducible quantification of single peptides across samples. However, such studies often require labeling of hundreds of micrograms of peptides such that the cost for labeling reagents represents a major contribution to the overall cost of an experiment. Here, we describe and evaluate a robust and cost-effective protocol for TMT labeling that reduces the quantity of required labeling reagent by a factor of eight and achieves complete labeling. Under- and overlabeling of peptides derived from complex digests of tissues and cell lines were systematically evaluated using peptide quantities of between 12.5 and 800 µg and TMT-to-peptide ratios (wt/wt) ranging from 8:1 to 1:2 at different TMT and peptide concentrations. When reaction volumes were reduced to maintain TMT and peptide concentrations of at least 10 mm and 2 g/l, respectively, TMT-to-peptide ratios as low as 1:1 (wt/wt) resulted in labeling efficiencies of > 99% and excellent intra- and interlaboratory reproducibility. The utility of the optimized protocol was further demonstrated in a deep-scale proteome and phosphoproteome analysis of patient-derived xenograft tumor tissue benchmarked against the labeling procedure recommended by the TMT vendor. Finally, we discuss the impact of labeling reaction parameters for N-hydroxysuccinimide ester-based chemistry and provide guidance on adopting efficient labeling protocols for different peptide quantities.