Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection.

BACKGROUND: Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors. METHODS: We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers. RESULTS: We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations in TP53, PLXNB3, Zinc finger genes, and fewer mutations in SDC2, CDH1, PIK3CA, SVIL, and PTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA-protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, including KIAA1522, DCK, FOXO3, MYO9B, ARID1A, EPRS, ZC3HAV1, and RBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors. CONCLUSIONS: This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

Mapping the proteogenomic landscape enables prediction of drug response in acute myeloid leukemia.

Acute myeloid leukemia is a poor-prognosis cancer commonly stratified by genetic aberrations, but these mutations are often heterogeneous and fail to consistently predict therapeutic response. Here, we combine transcriptomic, proteomic, and phosphoproteomic datasets with ex vivo drug sensitivity data to help understand the underlying pathophysiology of AML beyond mutations. We measure the proteome and phosphoproteome of 210 patients and combine them with genomic and transcriptomic measurements to identify four proteogenomic subtypes that complement existing genetic subtypes. We build a predictor to classify samples into subtypes and map them to a "landscape" that identifies specific drug response patterns. We then build a drug response prediction model to identify drugs that target distinct subtypes and validate our findings on cell lines representing various stages of quizartinib resistance. Our results show how multiomics data together with drug sensitivity data can inform therapy stratification and drug combinations in AML.

Mass Spectrometry-Based Proteogenomics: New Therapeutic Opportunities for Precision Medicine.

Proteogenomics refers to the integration of comprehensive genomic, transcriptomic, and proteomic measurements from the same samples with the goal of fully understanding the regulatory processes converting genotypes to phenotypes, often with an emphasis on gaining a deeper understanding of disease processes. Although specific genetic mutations have long been known to drive the development of multiple cancers, gene mutations alone do not always predict prognosis or response to targeted therapy. The benefit of proteogenomics research is that information obtained from proteins and their corresponding pathways provides insight into therapeutic targets that can complement genomic information by providing an additional dimension regarding the underlying mechanisms and pathophysiology of tumors. This review describes the novel insights into tumor biology and drug resistance derived from proteogenomic analysis while highlighting the clinical potential of proteogenomic observations and advances in technique and analysis tools.

Validation of a Proteomic Signature of Lung Cancer Risk from Bronchial Specimens of Risk-Stratified Individuals.

A major challenge in lung cancer prevention and cure hinges on identifying the at-risk population that ultimately develops lung cancer. Previously, we reported proteomic alterations in the cytologically normal bronchial epithelial cells collected from the bronchial brushings of individuals at risk for lung cancer. The purpose of this study is to validate, in an independent cohort, a selected list of 55 candidate proteins associated with risk for lung cancer with sensitive targeted proteomics using selected reaction monitoring (SRM). Bronchial brushings collected from individuals at low and high risk for developing lung cancer as well as patients with lung cancer, from both a subset of the original cohort (batch 1: n = 10 per group) and an independent cohort of 149 individuals (batch 2: low risk (n = 32), high risk (n = 34), and lung cancer (n = 83)), were analyzed using multiplexed SRM assays. ALDH3A1 and AKR1B10 were found to be consistently overexpressed in the high-risk group in both batch 1 and batch 2 brushing specimens as well as in the biopsies of batch 1. Validation of highly discriminatory proteins and metabolic enzymes by SRM in a larger independent cohort supported their use to identify patients at high risk for developing lung cancer.

Integrative multi-omic cancer profiling reveals DNA methylation patterns associated with therapeutic vulnerability and cell-of-origin.

DNA methylation plays a critical role in establishing and maintaining cellular identity. However, it is frequently dysregulated during tumor development and is closely intertwined with other genetic alterations. Here, we leveraged multi-omic profiling of 687 tumors and matched non-involved adjacent tissues from the kidney, brain, pancreas, lung, head and neck, and endometrium to identify aberrant methylation associated with RNA and protein abundance changes and build a Pan-Cancer catalog. We uncovered lineage-specific epigenetic drivers including hypomethylated FGFR2 in endometrial cancer. We showed that hypermethylated STAT5A is associated with pervasive regulon downregulation and immune cell depletion, suggesting that epigenetic regulation of STAT5A expression constitutes a molecular switch for immunosuppression in squamous tumors. We further demonstrated that methylation subtype-enrichment information can explain cell-of-origin, intra-tumor heterogeneity, and tumor phenotypes. Overall, we identified cis-acting DNA methylation events that drive transcriptional and translational changes, shedding light on the tumor's epigenetic landscape and the role of its cell-of-origin.

Deep learning integrates histopathology and proteogenomics at a pan-cancer level.

We introduce a pioneering approach that integrates pathology imaging with transcriptomics and proteomics to identify predictive histology features associated with critical clinical outcomes in cancer. We utilize 2,755 H&E-stained histopathological slides from 657 patients across 6 cancer types from CPTAC. Our models effectively recapitulate distinctions readily made by human pathologists: tumor vs. normal (AUROC = 0.995) and tissue-of-origin (AUROC = 0.979). We further investigate predictive power on tasks not normally performed from H&E alone, including TP53 prediction and pathologic stage. Importantly, we describe predictive morphologies not previously utilized in a clinical setting. The incorporation of transcriptomics and proteomics identifies pathway-level signatures and cellular processes driving predictive histology features. Model generalizability and interpretability is confirmed using TCGA. We propose a classification system for these tasks, and suggest potential clinical applications for this integrated human and machine learning approach. A publicly available web-based platform implements these models.

Proteome Mapping of the Human Pancreatic Islet Microenvironment Reveals Endocrine-Exocrine Signaling Sphere of Influence.

The need for a clinically accessible method with the ability to match protein activity within heterogeneous tissues is currently unmet by existing technologies. Our proteomics sample preparation platform, named microPOTS (Microdroplet Processing in One pot for Trace Samples), can be used to measure relative protein abundance in micron-scale samples alongside the spatial location of each measurement, thereby tying biologically interesting proteins and pathways to distinct regions. However, given the smaller pixel/voxel number and amount of tissue measured, standard mass spectrometric analysis pipelines have proven inadequate. Here we describe how existing computational approaches can be adapted to focus on the specific biological questions asked in spatial proteomics experiments. We apply this approach to present an unbiased characterization of the human islet microenvironment comprising the entire complex array of cell types involved while maintaining spatial information and the degree of the islet's sphere of influence. We identify specific functional activity unique to the pancreatic islet cells and demonstrate how far their signature can be detected in the adjacent tissue. Our results show that we can distinguish pancreatic islet cells from the neighboring exocrine tissue environment, recapitulate known biological functions of islet cells, and identify a spatial gradient in the expression of RNA processing proteins within the islet microenvironment.

Use of Longitudinal Serum Analysis and Machine Learning to Develop a Classifier for Cancer Early Detection.

Early detection of solid tumors through a simple screening process, such as the proteomic analysis of biofluids, has the potential to significantly alter the management and outcomes of cancers. The application of advanced targeted proteomics measurements and data analysis strategies to uniformly collected serum or plasma samples would enable longitudinal studies of cancer risk, progression, and response to therapy that have the potential to significantly reduce cancer burden in general. In this article, we describe a generalizable workflow combining robust, multiplexed targeted proteomics measurements applied to longitudinal samples from the Department of Defense Serum Repository with a Random Forest machine learning method for developing and initially evaluating the performance of candidate biomarker panels for early detection of cancers. The effectiveness of this approach was demonstrated in a cohort of 175 head and neck squamous cell carcinoma patients. The outlined protocols include methods for sample preparation, instrument analysis, and data analysis and interpretation using this workflow.

A streamlined tandem tip-based workflow for sensitive nanoscale phosphoproteomics.

Effective phosphoproteome of nanoscale sample analysis remains a daunting task, primarily due to significant sample loss associated with non-specific surface adsorption during enrichment of low stoichiometric phosphopeptide. We develop a tandem tip phosphoproteomics sample preparation method that is capable of sample cleanup and enrichment without additional sample transfer, and its integration with our recently developed SOP (Surfactant-assisted One-Pot sample preparation) and iBASIL (improved Boosting to Amplify Signal with Isobaric Labeling) approaches provides a streamlined workflow enabling sensitive, high-throughput nanoscale phosphoproteome measurements. This approach significantly reduces both sample loss and processing time, allowing the identification of >3000 (>9500) phosphopeptides from 1 (10) µg of cell lysate using the label-free method without a spectral library. It also enables precise quantification of ~600 phosphopeptides from 100 sorted cells (single-cell level input for the enriched phosphopeptides) and ~700 phosphopeptides from human spleen tissue voxels with a spatial resolution of 200 µm (equivalent to ~100 cells) in a high-throughput manner. The new workflow opens avenues for phosphoproteome profiling of mass-limited samples at the low nanogram level.

Histopathologic and proteogenomic heterogeneity reveals features of clear cell renal cell carcinoma aggressiveness.

Clear cell renal cell carcinomas (ccRCCs) represent ∼75% of RCC cases and account for most RCC-associated deaths. Inter- and intratumoral heterogeneity (ITH) results in varying prognosis and treatment outcomes. To obtain the most comprehensive profile of ccRCC, we perform integrative histopathologic, proteogenomic, and metabolomic analyses on 305 ccRCC tumor segments and 166 paired adjacent normal tissues from 213 cases. Combining histologic and molecular profiles reveals ITH in 90% of ccRCCs, with 50% demonstrating immune signature heterogeneity. High tumor grade, along with BAP1 mutation, genome instability, increased hypermethylation, and a specific protein glycosylation signature define a high-risk disease subset, where UCHL1 expression displays prognostic value. Single-nuclei RNA sequencing of the adverse sarcomatoid and rhabdoid phenotypes uncover gene signatures and potential insights into tumor evolution. In vitro cell line studies confirm the potential of inhibiting identified phosphoproteome targets. This study molecularly stratifies aggressive histopathologic subtypes that may inform more effective treatment strategies.

A prognostic risk model for glioma patients by systematic evaluation of genomic variations.

The overall survival rate of gliomas has not significantly improved despite new effective treatments, mainly due to tumor heterogeneity and drug delivery. Here, we perform an integrated clinic-genomic analysis of 1, 477 glioma patients from a Chinese cohort and a TCGA cohort and propose a potential prognostic model for gliomas. We identify that SBS11 and SBS23 mutational signatures are associated with glioma recurrence and indicate worse prognosis only in low-grade type of gliomas and IDH-Mut subtype. We also identify 42 genomic features associated with distinct clinical outcome and successfully used ten of these to develop a prognostic risk model of gliomas. The high-risk glioma patients with shortened survival were characterized by high level of frequent copy number alterations including PTEN, CDKN2A/B deletion, EGFR amplification, less IDH1 or CIC gene mutations, high infiltration levels of immunosuppressive cells and activation of G2M checkpoint and Oxidative phosphorylation oncogenic pathway.

New Views of Old Proteins: Clarifying the Enigmatic Proteome.

All human diseases involve proteins, yet our current tools to characterize and quantify them are limited. To better elucidate proteins across space, time, and molecular composition, we provide a >10 years of projection for technologies to meet the challenges that protein biology presents. With a broad perspective, we discuss grand opportunities to transition the science of proteomics into a more propulsive enterprise. Extrapolating recent trends, we describe a next generation of approaches to define, quantify, and visualize the multiple dimensions of the proteome, thereby transforming our understanding and interactions with human disease in the coming decade.

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.

Evaluation of Differential Peptide Loading on Tandem Mass Tag-Based Proteomic and Phosphoproteomic Data Quality.

Global and phosphoproteome profiling has demonstrated great utility for the analysis of clinical specimens. One barrier to the broad clinical application of proteomic profiling is the large amount of biological material required, particularly for phosphoproteomics─currently on the order of 25 mg wet tissue weight. For hematopoietic cancers such as acute myeloid leukemia (AML), the sample requirement is ≥10 million peripheral blood mononuclear cells (PBMCs). Across large study cohorts, this requirement will exceed what is obtainable for many individual patients/time points. For this reason, we were interested in the impact of differential peptide loading across multiplex channels on proteomic data quality. To achieve this, we tested a range of channel loading amounts (approximately the material obtainable from 5E5, 1E6, 2.5E6, 5E6, and 1E7 AML patient cells) to assess proteome coverage, quantification precision, and peptide/phosphopeptide detection in experiments utilizing isobaric tandem mass tag (TMT) labeling. As expected, fewer missing values were observed in TMT channels with higher peptide loading amounts compared to lower loadings. Moreover, channels with a lower loading have greater quantitative variability than channels with higher loadings. A statistical analysis showed that decreased loading amounts result in an increase in the type I error rate. We then examined the impact of differential loading on the detection of known differences between distinct AML cell lines. Similar patterns of increased data missingness and higher quantitative variability were observed as loading was decreased resulting in fewer statistical differences; however, we found good agreement in features identified as differential, demonstrating the value of this approach.

CD63-mediated cloaking of VEGF in small extracellular vesicles contributes to anti-VEGF therapy resistance.

Despite wide use of anti-vascular endothelial growth factor (VEGF) therapy for many solid cancers, most individuals become resistant to this therapy, leading to disease progression. Therefore, new biomarkers and strategies for blocking adaptive resistance of cancer to anti-VEGF therapy are needed. As described here, we demonstrate that cancer-derived small extracellular vesicles package increasing quantities of VEGF and other factors in response to anti-VEGF therapy. The packaging process of VEGF into small extracellular vesicles (EVs) is mediated by the tetraspanin CD63. Furthermore, small EV-VEGF (eVEGF) is not accessible to anti-VEGF antibodies and can trigger intracrine VEGF signaling in endothelial cells. eVEGF promotes angiogenesis and enhances tumor growth despite bevacizumab treatment. These data demonstrate a mechanism where VEGF is partitioned into small EVs and promotes tumor angiogenesis and progression. These findings have clinical implications for biomarkers and therapeutic strategies for ovarian cancer.

Facile One-Pot Nanoproteomics for Label-Free Proteome Profiling of 50-1000 Mammalian Cells.

Recent advances in sample preparation enable label-free mass spectrometry (MS)-based proteome profiling of small numbers of mammalian cells. However, specific devices are often required to downscale sample processing volume from the standard 50-200 µL to sub-µL for effective nanoproteomics, which greatly impedes the implementation of current nanoproteomics methods by the proteomics research community. Herein, we report a facile one-pot nanoproteomics method termed SOPs-MS (surfactant-assisted one-pot sample processing at the standard volume coupled with MS) for convenient robust proteome profiling of 50-1000 mammalian cells. Building upon our recent development of SOPs-MS for label-free single-cell proteomics at a low µL volume, we have systematically evaluated its processing volume at 10-200 µL using 100 human cells. The processing volume of 50 µL that is in the range of volume for standard proteomics sample preparation has been selected for easy sample handling with a benchtop micropipette. SOPs-MS allows for reliable label-free quantification of ∼1200-2700 protein groups from 50 to 1000 MCF10A cells. When applied to small subpopulations of mouse colon crypt cells, SOPs-MS has revealed protein signatures between distinct subpopulation cells with identification of ∼1500-2500 protein groups for each subpopulation. SOPs-MS may pave the way for routine deep proteome profiling of small numbers of cells and low-input samples.

The AML microenvironment catalyzes a stepwise evolution to gilteritinib resistance.

Our study details the stepwise evolution of gilteritinib resistance in FLT3-mutated acute myeloid leukemia (AML). Early resistance is mediated by the bone marrow microenvironment, which protects residual leukemia cells. Over time, leukemia cells evolve intrinsic mechanisms of resistance, or late resistance. We mechanistically define both early and late resistance by integrating whole-exome sequencing, CRISPR-Cas9, metabolomics, proteomics, and pharmacologic approaches. Early resistant cells undergo metabolic reprogramming, grow more slowly, and are dependent upon Aurora kinase B (AURKB). Late resistant cells are characterized by expansion of pre-existing NRAS mutant subclones and continued metabolic reprogramming. Our model closely mirrors the timing and mutations of AML patients treated with gilteritinib. Pharmacological inhibition of AURKB resensitizes both early resistant cell cultures and primary leukemia cells from gilteritinib-treated AML patients. These findings support a combinatorial strategy to target early resistant AML cells with AURKB inhibitors and gilteritinib before the expansion of pre-existing resistance mutations occurs.

Assessment of TMT Labeling Efficiency in Large-Scale Quantitative Proteomics: The Critical Effect of Sample pH.

Isobaric labeling via tandem mass tag (TMT) reagents enables sample multiplexing prior to LC-MS/MS, facilitating high-throughput large-scale quantitative proteomics. Consistent and efficient labeling reactions are essential to achieve robust quantification; therefore, embedded in our clinical proteomic protocol is a quality control (QC) sample that contains a small aliquot from each sample within a TMT set, referred to as "Mixing QC." This Mixing QC enables the detection of TMT labeling issues by LC-MS/MS before combining the full samples to allow for salvaging of poor TMT labeling reactions. While TMT labeling is a valuable tool, factors leading to poor reactions are not fully studied. We observed that relabeling does not necessarily rescue TMT reactions and that peptide samples sometimes remained acidic after resuspending in 50 mM HEPES buffer (pH 8.5), which coincided with low labeling efficiency (LE) and relatively low median reporter ion intensities (MRIIs). To obtain a more resilient TMT labeling procedure, we investigated LE, reporter ion missingness, the ratio of mean TMT set MRII to individual channel MRII, and the distribution of log 2 reporter ion ratios of Mixing QC samples. We discovered that sample pH is a critical factor in LE, and increasing the buffer concentration in poorly labeled samples before relabeling resulted in the successful rescue of TMT labeling reactions. Moreover, resuspending peptides in 500 mM HEPES buffer for TMT labeling resulted in consistently higher LE and lower missing data. By better controlling the sample pH for labeling and implementing multiple methods for assessing labeling quality before combining samples, we demonstrate that robust TMT labeling for large-scale quantitative studies is achievable.