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.

Proteogenomic analysis of chemo-refractory high-grade serous ovarian cancer.

To improve the understanding of chemo-refractory high-grade serous ovarian cancers (HGSOCs), we characterized the proteogenomic landscape of 242 (refractory and sensitive) HGSOCs, representing one discovery and two validation cohorts across two biospecimen types (formalin-fixed paraffin-embedded and frozen). We identified a 64-protein signature that predicts with high specificity a subset of HGSOCs refractory to initial platinum-based therapy and is validated in two independent patient cohorts. We detected significant association between lack of Ch17 loss of heterozygosity (LOH) and chemo-refractoriness. Based on pathway protein expression, we identified 5 clusters of HGSOC, which validated across two independent patient cohorts and patient-derived xenograft (PDX) models. These clusters may represent different mechanisms of refractoriness and implicate putative therapeutic vulnerabilities.

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.

Clinical Proteomics for Solid Organ Tissues.

The evaluation of biopsied solid organ tissue has long relied on visual examination using a microscope. Immunohistochemistry is critical in this process, labeling and detecting cell lineage markers and therapeutic targets. However, while the practice of immunohistochemistry has reshaped diagnostic pathology and facilitated improvements in cancer treatment, it has also been subject to pervasive challenges with respect to standardization and reproducibility. Efforts are ongoing to improve immunohistochemistry, but for some applications, the benefit of such initiatives could be impeded by its reliance on monospecific antibody-protein reagents and limited multiplexing capacity. This perspective surveys the relevant challenges facing traditional immunohistochemistry and describes how mass spectrometry, particularly liquid chromatography-tandem mass spectrometry, could help alleviate problems. In particular, targeted mass spectrometry assays could facilitate measurements of individual proteins or analyte panels, using internal standards for more robust quantification and improved interlaboratory reproducibility. Meanwhile, untargeted mass spectrometry, showcased to date clinically in the form of amyloid typing, is inherently multiplexed, facilitating the detection and crude quantification of 100s to 1000s of proteins in a single analysis. Further, data-independent acquisition has yet to be applied in clinical practice, but offers particular strengths that could appeal to clinical users. Finally, we discuss the guidance that is needed to facilitate broader utilization in clinical environments and achieve standardization.

Monitoring Both Extended and Tryptic Forms of Stable Isotope-Labeled Standard Peptides Provides an Internal Quality Control of Proteolytic Digestion in Targeted Mass Spectrometry-Based Assays.

Targeted mass spectrometry (MS)-based proteomic assays, such as multiplexed multiple reaction monitoring (MRM)-MS assays, enable sensitive and specific quantification of proteotypic peptides as stoichiometric surrogates for proteins. Efforts are underway to expand the use of MRM-MS assays in clinical environments, which requires a reliable strategy to monitor proteolytic digestion efficiency within individual samples. Towards this goal, extended stable isotope-labeled standard (SIS) peptides (hE), which incorporate native proteolytic cleavage sites, can be spiked into protein lysates prior to proteolytic (trypsin) digestion, and release of the tryptic SIS peptide (hT) can be monitored. However, hT measurements alone cannot monitor the extent of digestion and may be confounded by matrix effects specific to individual patient samples; therefore, they are not sufficient to monitor sample-to-sample digestion variability. We hypothesized that measuring undigested hE, along with its paired hT, would improve detection of digestion issues compared to only measuring hT. We tested the ratio of the SIS pair measurements, or hE/hT, as a quality control (QC) metric of trypsin digestion for two MRM assays: a direct-MRM (398 targets) and an immuno-MRM (126 targets requiring immunoaffinity peptide enrichment) assay, with extended SIS peptides observable for 54% (216) and 62% (78) of the targets, respectively. We evaluated the quantitative bias for each target in a series of experiments that adversely affected proteolytic digestion (e.g., variable digestion times, pH, and temperature). We identified a subset of SIS pairs (36 for the direct-MRM, 7 for the immuno-MRM assay) for which the hE/hT ratio reliably detected inefficient digestion that resulted in decreased assay sensitivity and unreliable endogenous quantification. The hE/hT ratio was more responsive to a decrease in digestion efficiency than a metric based on hT measurements alone. For clinical-grade MRM-MS assays, this study describes a ready-to-use QC panel and also provides a road map for designing custom QC panels.

Fanconi anemia-isogenic head and neck cancer cell line pairs: A basic and translational science resource.

Fanconi anemia (FA) is a heritable malformation, bone marrow failure and cancer predisposition syndrome that confers an exceptionally high risk of squamous carcinomas. These carcinomas originate in epithelia lining the mouth, proximal esophagus, vulva and anus: their origins are not understood, and no effective ways have been identified to prevent or delay their appearance. Many FA-associated carcinomas are also therapeutically challenging: they may be multi-focal and stage-advanced at diagnosis, and most individuals with FA cannot tolerate standard-of-care systemic therapies such as DNA cross-linking drugs or ionizing radiation due to constitutional DNA damage hypersensitivity. We developed the Fanconi Anemia Cancer Cell Line Resource (FA-CCLR) to foster new work on the origins, treatment and prevention of FA-associated carcinomas. The FA-CCLR consists of Fanconi-isogenic head and neck squamous cell carcinoma (HNSCC) cell line pairs generated from five individuals with FA-associated HNSCC, and five individuals with sporadic HNSCC. Sporadic, isogenic HNSCC cell line pairs were generated in parallel with FA patient-derived isogenic cell line pairs to provide comparable experimental material to use to identify cell and molecular phenotypes driven by germline or somatic loss of Fanconi pathway function, and the subset of these FA-dependent phenotypes that can be modified, complemented or suppressed. All 10 FANC-isogenic cell line pairs are available to academic, non-profit and industry investigators via the "Fanconi Anemia Research Materials" Resource and Repository at Oregon Health & Sciences University, Portland OR.

Hormone Receptor-Positive Breast Cancer Sensitive to Pembrolizumab: Evidence of the Pathogenicity of the MLH1 Variant 1835del3.

A woman with estrogen/progesterone receptor-positive, ERBB2-negative metastatic breast cancer developed progressive disease despite treatment with multiple hormonal and chemotherapeutic modalities. She carried a germline variant of MLH1 (1835del3), also known as c.1835_1837del and v612del, the pathogenicity of which has not been conclusively determined. MLH1 staining was not seen on immunohistochemical staining of her tumor tissue. The patient experienced a >5-year dramatic response to 4 doses of pembrolizumab. Family studies revealed multiple other relatives with the MLH1 1835del3 variant, as well as multiple relatives with colon cancer. The one relative with colon cancer who underwent genetic testing demonstrated the same variant. Laboratory studies revealed that the patient's tumor showed loss of heterozygosity (LOH) in the MLH1 region, high levels of microsatellite instability, and a high tumor mutational burden. LOH in the MLH1 region, along with the remarkable clinical response to pembrolizumab treatment and the presence of the same MLH1 variant in affected relatives, supports the hypothesis that the MLH1 1835del3 variant is pathogenic. Given the patient's family history, this likely represents an uncommon presentation of Lynch syndrome. Physicians should be alert to evaluate patients for targetable genetic variants even in unlikely clinical situations such as the one described here.

DAGBagM: learning directed acyclic graphs of mixed variables with an application to identify protein biomarkers for treatment response in ovarian cancer.

BACKGROUND: Applying directed acyclic graph (DAG) models to proteogenomic data has been shown effective for detecting causal biomarkers of complex diseases. However, there remain unsolved challenges in DAG learning to jointly model binary clinical outcome variables and continuous biomarker measurements. RESULTS: In this paper, we propose a new tool, DAGBagM, to learn DAGs with both continuous and binary nodes. By using appropriate models, DAGBagM allows for either continuous or binary nodes to be parent or child nodes. It employs a bootstrap aggregating strategy to reduce false positives in edge inference. At the same time, the aggregation procedure provides a flexible framework to robustly incorporate prior information on edges. CONCLUSIONS: Through extensive simulation experiments, we demonstrate that DAGBagM has superior performance compared to alternative strategies for modeling mixed types of nodes. In addition, DAGBagM is computationally more efficient than two competing methods. When applying DAGBagM to proteogenomic datasets from ovarian cancer studies, we identify potential protein biomarkers for platinum refractory/resistant response in ovarian cancer. DAGBagM is made available as a github repository at https://github.com/jie108/dagbagM .

Internal Standard Triggered-Parallel Reaction Monitoring Mass Spectrometry Enables Multiplexed Quantification of Candidate Biomarkers in Plasma.

Despite advances in proteomic technologies, clinical translation of plasma biomarkers remains low, partly due to a major bottleneck between the discovery of candidate biomarkers and costly clinical validation studies. Due to a dearth of multiplexable assays, generally only a few candidate biomarkers are tested, and the validation success rate is accordingly low. Previously, mass spectrometry-based approaches have been used to fill this gap but feature poor quantitative performance and were generally limited to hundreds of proteins. Here, we demonstrate the capability of an internal standard triggered-parallel reaction monitoring (IS-PRM) assay to greatly expand the numbers of candidates that can be tested with improved quantitative performance. The assay couples immunodepletion and fractionation with IS-PRM and was developed and implemented in human plasma to quantify 5176 peptides representing 1314 breast cancer biomarker candidates. Characterization of the IS-PRM assay demonstrated the precision (median % CV of 7.7%), linearity (median R2 > 0.999 over 4 orders of magnitude), and sensitivity (median LLOQ < 1 fmol, approximately) to enable rank-ordering of candidate biomarkers for validation studies. Using three plasma pools from breast cancer patients and three control pools, 893 proteins were quantified, of which 162 candidate biomarkers were verified in at least one of the cancer pools and 22 were verified in all three cancer pools. The assay greatly expands capabilities for quantification of large numbers of proteins and is well suited for prioritization of viable candidate biomarkers.

Proteogenomics connects somatic mutations to signalling in breast cancer.

Somatic mutations have been extensively characterized in breast cancer, but the effects of these genetic alterations on the proteomic landscape remain poorly understood. Here we describe quantitative mass-spectrometry-based proteomic and phosphoproteomic analyses of 105 genomically annotated breast cancers, of which 77 provided high-quality data. Integrated analyses provided insights into the somatic cancer genome including the consequences of chromosomal loss, such as the 5q deletion characteristic of basal-like breast cancer. Interrogation of the 5q trans-effects against the Library of Integrated Network-based Cellular Signatures, connected loss of CETN3 and SKP1 to elevated expression of epidermal growth factor receptor (EGFR), and SKP1 loss also to increased SRC tyrosine kinase. Global proteomic data confirmed a stromal-enriched group of proteins in addition to basal and luminal clusters, and pathway analysis of the phosphoproteome identified a G-protein-coupled receptor cluster that was not readily identified at the mRNA level. In addition to ERBB2, other amplicon-associated highly phosphorylated kinases were identified, including CDK12, PAK1, PTK2, RIPK2 and TLK2. We demonstrate that proteogenomic analysis of breast cancer elucidates the functional consequences of somatic mutations, narrows candidate nominations for driver genes within large deletions and amplified regions, and identifies therapeutic targets.

Quantification of Human Epidermal Growth Factor Receptor 2 by Immunopeptide Enrichment and Targeted Mass Spectrometry in Formalin-Fixed Paraffin-Embedded and Frozen Breast Cancer Tissues.

BACKGROUND: Conventional HER2-targeting therapies improve outcomes for patients with HER2-positive breast cancer (BC), defined as tumors showing HER2 protein overexpression by immunohistochemistry and/or ERBB2 gene amplification determined by in situ hybridization (ISH). Emerging HER2-targeting compounds show benefit in some patients with neither HER2 protein overexpression nor ERBB2 gene amplification, creating a need for new assays to select HER2-low tumors for treatment with these compounds. We evaluated the analytical performance of a targeted mass spectrometry-based assay for quantifying HER2 protein in formalin-fixed paraffin-embedded (FFPE) and frozen BC biopsies. METHODS: We used immunoaffinity-enrichment coupled to multiple reaction monitoring-mass spectrometry (immuno-MRM-MS) to quantify HER2 protein (as peptide GLQSLPTHDPSPLQR) in 96 frozen and 119 FFPE BC biopsies. We characterized linearity, lower limit of quantification (LLOQ), and intra- and inter-day variation of the assay in frozen and FFPE tissue matrices. We determined concordance between HER2 immuno-MRM-MS and predicate immunohistochemistry and ISH assays and examined the benefit of multiplexing the assay to include proteins expressed in tumor subcompartments (e.g., stroma, adipose, lymphocytes, epithelium) to account for tissue heterogeneity. RESULTS: HER2 immuno-MRM-MS assay linearity was ≥103, assay coefficient of variation was 7.8% (FFPE) and 5.9% (frozen) for spiked-in analyte, and 7.7% (FFPE) and 7.9% (frozen) for endogenous measurements. Immuno-MRM-MS-based HER2 measurements strongly correlated with predicate assay HER2 determinations, and concordance was improved by normalizing to glyceraldehyde-3-phosphate dehydrogenase. HER2 was quantified above the LLOQ in all tumors. CONCLUSIONS: Immuno-MRM-MS can be used to quantify HER2 in FFPE and frozen BC biopsies, even at low HER2 expression levels.

Integrative Proteo-genomic Analysis to Construct CNA-protein Regulatory Map in Breast and Ovarian Tumors.

Recent development in high throughput proteomics and genomics profiling enable one to study regulations of genome alterations on protein activities in a systematic manner. In this article, we propose a new statistical method, ProMAP, to systematically characterize the regulatory relationships between proteins and DNA copy number alterations (CNA) in breast and ovarian tumors based on proteogenomic data from the CPTAC-TCGA studies. Because of the dynamic nature of mass spectrometry instruments, proteomics data from labeled mass spectrometry experiments usually have non-ignorable batch effects. Moreover, mass spectrometry based proteomic data often possesses high percentages of missing values and non-ignorable missing-data patterns. Thus, we use a linear mixed effects model to account for the batch structure and explicitly incorporate the abundance-dependent-missing-data mechanism of proteomic data in ProMAP. In addition, we employ a multivariate regression framework to characterize the multiple-to-multiple regulatory relationships between CNA and proteins. Further, we use proper statistical regularization to facilitate the detection of master genetic regulators, which affect the activities of many proteins and often play important roles in genetic regulatory networks. Improved performance of ProMAP over existing methods were illustrated through extensive simulation studies and real data examples. Applying ProMAP to the CPTAC-TCGA breast and ovarian cancer data sets, we identified many genome regions, including a few novel ones, whose CNA were associated with protein and or phosphoprotein abundances. For example, in breast tumors, a small region in 8p11.21 was recognized as the second biggest hub in the CNA-phosphoprotein regulatory map, and further investigation of the regulatory targets suggests the potential role of 8p11.21 CNA in perturbing oxygen binding and transport activities in tumor cells. This and other findings from our analyses help to characterize the impacts of CNAs on protein activity landscapes and cast light on the genetic regulation mechanisms underlying these tumors.

Proteome Profiling Outperforms Transcriptome Profiling for Coexpression Based Gene Function Prediction.

Coexpression of mRNAs under multiple conditions is commonly used to infer cofunctionality of their gene products despite well-known limitations of this "guilt-by-association" (GBA) approach. Recent advancements in mass spectrometry-based proteomic technologies have enabled global expression profiling at the protein level; however, whether proteome profiling data can outperform transcriptome profiling data for coexpression based gene function prediction has not been systematically investigated. Here, we address this question by constructing and analyzing mRNA and protein coexpression networks for three cancer types with matched mRNA and protein profiling data from The Cancer Genome Atlas (TCGA) and the Clinical Proteomic Tumor Analysis Consortium (CPTAC). Our analyses revealed a marked difference in wiring between the mRNA and protein coexpression networks. Whereas protein coexpression was driven primarily by functional similarity between coexpressed genes, mRNA coexpression was driven by both cofunction and chromosomal colocalization of the genes. Functionally coherent mRNA modules were more likely to have their edges preserved in corresponding protein networks than functionally incoherent mRNA modules. Proteomic data strengthened the link between gene expression and function for at least 75% of Gene Ontology (GO) biological processes and 90% of KEGG pathways. A web application Gene2Net (http://cptac.gene2net.org) developed based on the three protein coexpression networks revealed novel gene-function relationships, such as linking ERBB2 (HER2) to lipid biosynthetic process in breast cancer, identifying PLG as a new gene involved in complement activation, and identifying AEBP1 as a new epithelial-mesenchymal transition (EMT) marker. Our results demonstrate that proteome profiling outperforms transcriptome profiling for coexpression based gene function prediction. Proteomics should be integrated if not preferred in gene function and human disease studies.

pRAD50: a novel and clinically applicable pharmacodynamic biomarker of both ATM and ATR inhibition identified using mass spectrometry and immunohistochemistry.

BACKGROUND: AZD0156 and AZD6738 are potent and selective inhibitors of ataxia-telangiectasia-kinase (ATM) and ataxia-telangiectasia-mutated and Rad3-related (ATR), respectively, important sensors/signallers of DNA damage. METHODS: We used multiplexed targeted-mass-spectrometry to select pRAD50(Ser635) as a pharmacodynamic biomarker for AZD0156-mediated ATM inhibition from a panel of 45 peptides, then developed and tested a clinically applicable immunohistochemistry assay for pRAD50(Ser635) detection in FFPE tissue. RESULTS: We found moderate pRAD50 baseline levels across cancer indications. pRAD50 was detectable in 100% gastric cancers (n = 23), 99% colorectal cancers (n = 102), 95% triple-negative-breast cancers (TNBC) (n = 40) and 87.5% glioblastoma-multiformes (n = 16). We demonstrated AZD0156 target inhibition in TNBC patient-derived xenograft models; where AZD0156 monotherapy or post olaparib treatment, resulted in a 34-72% reduction in pRAD50. Similar inhibition of pRAD50 (68%) was observed following ATM inhibitor treatment post irinotecan in a colorectal cancer xenograft model. ATR inhibition, using AZD6738, increased pRAD50 in the ATM-proficient models whilst in ATM-deficient models the opposite was observed, suggesting pRAD50 pharmacodynamics post ATR inhibition may be ATM-dependent and could be useful to determine ATM functionality in patients treated with ATR inhibitors. CONCLUSION: Together these data support clinical utilisation of pRAD50 as a biomarker of AZD0156 and AZD6738 pharmacology to elucidate clinical pharmacokinetic/pharmacodynamic relationships, thereby informing recommended Phase 2 dose/schedule.

Assessment of a method to characterize antibody selectivity and specificity for use in immunoprecipitation.

Antibodies are used in multiple cell biology applications, but there are no standardized methods to assess antibody quality-an absence that risks data integrity and reproducibility. We describe a mass spectrometry-based standard operating procedure for scoring immunoprecipitation antibody quality. We quantified the abundance of all the proteins in immunoprecipitates of 1,124 new recombinant antibodies for 152 chromatin-related human proteins by comparing normalized spectral abundance factors from the target antigen with those of all other proteins. We validated the performance of the standard operating procedure in blinded studies in five independent laboratories. Antibodies for which the target antigen or a member of its known protein complex was the most abundant protein were classified as 'IP gold standard'. This method generates quantitative outputs that can be stored and archived in public databases, and it represents a step toward a platform for community benchmarking of antibody quality.

Immobilized Metal Affinity Chromatography Coupled to Multiple Reaction Monitoring Enables Reproducible Quantification of Phospho-signaling.

A major goal in cell signaling research is the quantification of phosphorylation pharmacodynamics following perturbations. Traditional methods of studying cellular phospho-signaling measure one analyte at a time with poor standardization, rendering them inadequate for interrogating network biology and contributing to the irreproducibility of preclinical research. In this study, we test the feasibility of circumventing these issues by coupling immobilized metal affinity chromatography (IMAC)-based enrichment of phosphopeptides with targeted, multiple reaction monitoring (MRM) mass spectrometry to achieve precise, specific, standardized, multiplex quantification of phospho-signaling responses. A multiplex immobilized metal affinity chromatography- multiple reaction monitoring assay targeting phospho-analytes responsive to DNA damage was configured, analytically characterized, and deployed to generate phospho-pharmacodynamic curves from primary and immortalized human cells experiencing genotoxic stress. The multiplexed assays demonstrated linear ranges of ≥3 orders of magnitude, median lower limit of quantification of 0.64 fmol on column, median intra-assay variability of 9.3%, median inter-assay variability of 12.7%, and median total CV of 16.0%. The multiplex immobilized metal affinity chromatography- multiple reaction monitoring assay enabled robust quantification of 107 DNA damage-responsive phosphosites from human cells following DNA damage. The assays have been made publicly available as a resource to the community. The approach is generally applicable, enabling wide interrogation of signaling networks.

Quantification of ATP7B Protein in Dried Blood Spots by Peptide Immuno-SRM as a Potential Screen for Wilson's Disease.

Wilson's Disease (WD), a copper transport disorder caused by a genetic defect in the ATP7B gene, has been a long time strong candidate for newborn screening (NBS), since early interventions can give better results by preventing irreversible neurological disability or liver cirrhosis. Several previous pilot studies measuring ceruloplasmin (CP) in infants or children showed that this marker alone was insufficient to meet the universal screening for WD. WD results from mutations that cause absent or markedly diminished levels of ATP7B. Therefore, ATP7B could serve as a marker for the screening of WD, if the protein can be detected from dried blood spots (DBS). This study demonstrates that the immuno-SRM platform can quantify ATP7B in DBS in the picomolar range, and that the assay readily distinguishes affected cases from normal controls (p < 0.0001). The assay precision was <10% CV, and the protein was stable for a week in DBS at room temperature. These promising proof-of-concept data open up the possibility of screening WD in newborns and the potential for a multiplexed assay for screening a variety of congenital disorders using proteins as biomarkers in DBS.

Demonstrating the feasibility of large-scale development of standardized assays to quantify human proteins.

Multiple reaction monitoring (MRM) mass spectrometry has been successfully applied to monitor targeted proteins in biological specimens, raising the possibility that assays could be configured to measure all human proteins. We report the results of a pilot study designed to test the feasibility of a large-scale, international effort for MRM assay generation. We have configured, validated across three laboratories and made publicly available as a resource to the community 645 novel MRM assays representing 319 proteins expressed in human breast cancer. Assays were multiplexed in groups of >150 peptides and deployed to quantify endogenous analytes in a panel of breast cancer-related cell lines. The median assay precision was 5.4%, with high interlaboratory correlation (R(2) > 0.96). Peptide measurements in breast cancer cell lines were able to discriminate among molecular subtypes and identify genome-driven changes in the cancer proteome. These results establish the feasibility of a large-scale effort to develop an MRM assay resource.