Using Peptide-Level Proteomics Data for Detecting Differentially Expressed Proteins

ASPSCR1-TFE3 reprograms transcription by organizing enhancer loops around hexameric VCP/p97

The t(X,17) chromosomal translocation, generating the ASPSCR1::TFE3 fusion oncoprotein, is the singular genetic driver of alveolar soft part sarcoma (ASPS) and some Xp11-rearranged renal cell carcinomas (RCCs), frustrating efforts to identify therapeutic targets for these rare cancers. Here, proteomic analysis identifies VCP/p97, an AAA+ ATPase with known segregase function, as strongly enriched in co-immunoprecipitated nuclear complexes with ASPSCR1::TFE3. We demonstrate that VCP is a likely obligate co-factor of ASPSCR1::TFE3, one of the only such fusion oncoprotein co-factors identified in cancer biology. Specifically, VCP co-distributes with ASPSCR1::TFE3 across chromatin in association with enhancers genome-wide. VCP presence, its hexameric assembly, and its enzymatic function orchestrate the oncogenic transcriptional signature of ASPSCR1::TFE3, by facilitating assembly of higher-order chromatin conformation structures demonstrated by HiChIP. Finally, ASPSCR1::TFE3 and VCP demonstrate co-dependence for cancer cell proliferation and tumorigenesis in vitro and in ASPS and RCC mouse models, underscoring VCP's potential as a novel therapeutic target.

Uncovering hidden cancer self-dependencies through analysis of shRNA-level dependency scores

Large-scale short hairpin RNA (shRNA) screens on well-characterized human cancer cell lines have been widely used to identify novel cancer dependencies. However, the off-target effects of shRNA reagents pose a significant challenge in the analysis of these screens. To mitigate these off-target effects, various approaches have been proposed that aggregate different shRNA viability scores targeting a gene into a single gene-level viability score. Most computational methods for discovering cancer dependencies rely on these gene-level scores. In this paper, we propose a computational method, named NBDep, to find cancer self-dependencies by directly analyzing shRNA-level dependency scores instead of gene-level scores. The NBDep algorithm begins by removing known batch effects of the shRNAs and selecting a subset of concordant shRNAs for each gene. It then uses negative binomial random effects models to statistically assess the dependency between genetic alterations and the viabilities of cell lines by incorporating all shRNA dependency scores of each gene into the model. We applied NBDep to the shRNA dependency scores available at Project DRIVE, which covers 26 different types of cancer. The proposed method identified more well-known and putative cancer genes compared to alternative gene-level approaches in pan-cancer and cancer-specific analyses. Additionally, we demonstrated that NBDep controls type-I error and outperforms statistical tests based on gene-level scores in simulation studies.

ASPSCR1-TFE3 reprograms transcription by organizing enhancer loops around hexameric VCP/p97

The t(X,17) chromosomal translocation, generating the ASPSCR1-TFE3 fusion oncoprotein, is the singular genetic driver of alveolar soft part sarcoma (ASPS) and some Xp11-rearranged renal cell carcinomas (RCC), frustrating efforts to identify therapeutic targets for these rare cancers. Proteomic analysis showed that VCP/p97, an AAA+ ATPase with known segregase function, was strongly enriched in co-immunoprecipitated nuclear complexes with ASPSCR1-TFE3. We demonstrate that VCP is a likely obligate co-factor of ASPSCR1-TFE3, one of the only such fusion oncoprotein co-factors identified in cancer biology. Specifically, VCP co-distributed with ASPSCR1-TFE3 across chromatin in association with enhancers genome-wide. VCP presence, its hexameric assembly, and its enzymatic function orchestrated the oncogenic transcriptional signature of ASPSCR1-TFE3, by facilitating assembly of higher-order chromatin conformation structures as demonstrated by HiChIP. Finally, ASPSCR1-TFE3 and VCP demonstrated co-dependence for cancer cell proliferation and tumorigenesis in vitro and in ASPS and RCC mouse models, underscoring VCP's potential as a novel therapeutic target.

A MYCN-independent mechanism mediating secretome reprogramming and metastasis in MYCN-amplified neuroblastoma

MYCN amplification (MNA) is a defining feature of high-risk neuroblastoma (NB) and predicts poor prognosis. However, whether genes within or in close proximity to the MYCN amplicon also contribute to MNA+ NB remains poorly understood. Here, we identify that GREB1, a transcription factor encoding gene neighboring the MYCN locus, is frequently coexpressed with MYCN and promotes cell survival in MNA+ NB. GREB1 controls gene expression independently of MYCN, among which we uncover myosin 1B (MYO1B) as being highly expressed in MNA+ NB and, using a chick chorioallantoic membrane (CAM) model, as a crucial regulator of invasion and metastasis. Global secretome and proteome profiling further delineates MYO1B in regulating secretome reprogramming in MNA+ NB cells, and the cytokine MIF as an important pro-invasive and pro-metastatic mediator of MYO1B activity. Together, we have identified a putative GREB1-MYO1B-MIF axis as an unconventional mechanism promoting aggressive behavior in MNA+ NB and independently of MYCN.

pepDESC: A method for the detection of differentially expressed proteins for mass spectrometry-based single-cell proteomics using peptide-level information

Single-cell proteomics as an emerging field has exhibited potential in revealing cellular heterogeneity at the functional level. However, accurate interpretation of single-cell proteomics data suffers from challenges such as measurement noise, internal heterogeneity, and the limited sample size of label-free quantitative mass spectrometry. Herein, the author describes pepDESC, a method for detecting differentially expressed proteins using peptide-level information designed for label-free quantitative mass spectrometry-based single-cell proteomics. While in this study, the author focuses on the heterogeneity among limited number of samples, pepDESC is also applicable to regular-size proteomics data. By balancing proteome coverage and quantification accuracy using peptide quantification, pepDESC is demonstrated to be effective in real-world single-cell and spike -in benchmark datasets. By applying pepDESC to a published single-mouse macrophages data, the author found a large fraction of differentially expressed proteins among three types of cells, which remarkably revealed distinct dynamics of different cellular functions responding to lipopolysaccharide stimulation.

Challenges and opportunities for single cell computational proteomics

Single-cell proteomics is growing rapidly and has made several technological advancements. As most research has been focused on improving instrumentation and sample preparation methods, very little attention has been given to algorithms responsible for identifying and quantifying proteins. Given the inherent difference between bulk data and single-cell data, it's necessary to realize that current algorithms being employed on single-cell data were designed for bulk data, and have underlying assumptions that may not hold true for single-cell data. In order to develop and optimize algorithms for single-cell data, we need to characterize the differences between single-cell data and bulk data, and assess how current algorithms perform on single-cell data. Here, we present a review of algorithms responsible for identifying and quantifying peptides and proteins. We will give a review of how each type of algorithm works, assumptions it relies on, how it performs on single-cell data, and possible optimizations and solutions that could be used to address the differences in single-cell data.

The proteome of clear cell ovarian carcinoma

Clear cell ovarian carcinoma (CCOC) is the second most common subtype of epithelial ovarian carcinoma. Late-stage CCOC is not responsive to gold-standard chemotherapy and results in suboptimal outcomes for patients. In-depth molecular insight is urgently needed to stratify the disease and drive therapeutic development. We conducted global proteomics for 192 cases of CCOC and compared these with other epithelial ovarian carcinoma subtypes. Our results showed distinct proteomic differences in CCOC compared with other epithelial ovarian cancer subtypes including alterations in lipid and purine metabolism pathways. Furthermore, we report potential clinically significant proteomic subgroups within CCOC, suggesting the biologic plausibility of stratified treatment for this cancer. Taken together, our results provide a comprehensive understanding of the CCOC proteomic landscape to facilitate future understanding and research of this disease. © 2022 The Pathological Society of Great Britain and Ireland.

Accounting for multiple imputation-induced variability for differential analysis in mass spectrometry-based label-free quantitative proteomics

Imputing missing values is common practice in label-free quantitative proteomics. Imputation aims at replacing a missing value with a user-defined one. However, the imputation itself may not be optimally considered downstream of the imputation process, as imputed datasets are often considered as if they had always been complete. Hence, the uncertainty due to the imputation is not adequately taken into account. We provide a rigorous multiple imputation strategy, leading to a less biased estimation of the parameters’ variability thanks to Rubin’s rules. The imputation-based peptide’s intensities’ variance estimator is then moderated using Bayesian hierarchical models. This estimator is finally included in moderated t-test statistics to provide differential analyses results. This workflow can be used both at peptide and protein-level in quantification datasets. Indeed, an aggregation step is included for protein-level results based on peptide-level quantification data. Our methodology, named mi4p, was compared to the state-of-the-art limma workflow implemented in the DAPARR package, both on simulated and real datasets. We observed a trade-off between sensitivity and specificity, while the overall performance of mi4p outperforms DAPAR in terms of F-Score.

Statistical inference methods commonly used in quantitative proteomics are based on the measurement of peptide intensities. They allow the deduction of protein abundances provided that sufficient peptides per protein are available. However, they do not satisfactorily consider peptides or proteins whose intensities are missing under certain conditions, even though they are particularly interesting from a biological or medical point of view, since they may explain a difference between the groups being compared. Some state-of-the-art statistical proteomics data processing software proposes to impute these missing values, while others simply remove proteins with too many missing peptides. The statistical treatment is not entirely satisfactory when imputation methods are used, notably multiple imputation techniques. Indeed, even if these statistical tools are relevant in this context, the data sets once imputed are considered as having always been complete in the subsequent analyses: the uncertainty caused by the imputation is not taken into account. These analyses generally conclude with a study of the differences in protein abundances between the different conditions, either using Student’s or Welch’s test for the most rudimentary approaches or using the t-tempered testing techniques based on empirical Bayesian approaches. Thus, we propose a new methodology that starts by imputing missing values at the peptide level and estimating the uncertainty associated with this imputation and naturally extends by incorporating this uncertainty into the current moderated variance estimation techniques.

Proteomic analysis of archival breast cancer clinical specimens identifies biological subtypes with distinct survival outcomes

Despite advances in genomic classification of breast cancer, current clinical tests and treatment decisions are commonly based on protein level information. Formalin-fixed paraffin-embedded (FFPE) tissue specimens with extended clinical outcomes are widely available. Here, we perform comprehensive proteomic profiling of 300 FFPE breast cancer surgical specimens, 75 of each PAM50 subtype, from patients diagnosed in 2008-2013 (n = 178) and 1986-1992 (n = 122) with linked clinical outcomes. These two cohorts are analyzed separately, and we quantify 4214 proteins across all 300 samples. Within the aggressive PAM50-classified basal-like cases, proteomic profiling reveals two groups with one having characteristic immune hot expression features and highly favorable survival. Her2-Enriched cases separate into heterogeneous groups differing by extracellular matrix, lipid metabolism, and immune-response features. Within 88 triple-negative breast cancers, four proteomic clusters display features of basal-immune hot, basal-immune cold, mesenchymal, and luminal with disparate survival outcomes. Our proteomic analysis characterizes the heterogeneity of breast cancer in a clinically-applicable manner, identifies potential biomarkers and therapeutic targets, and provides a resource for clinical breast cancer classification.

Protein level information enables the identification of potential biomarkers and therapeutic targets for breast cancer. Here, the authors perform proteomic analysis of 2 cohorts of breast cancer surgical specimens and identify distinct subtypes, immune features and survival outcomes.

De novo and cell line models of human mammary cell transformation reveal an essential role for Yb-1 in multiple stages of human breast cancer

Breast cancer heterogeneity has made it challenging to identify mechanisms critical to the initial stages of their genesis in vivo. Here, we sought to interrogate the role of YB-1 in newly arising human breast cancers as well as in established cell lines. In a first series of experiments, we found that short-hairpin RNA-mediated knockdown of YB-1 in MDA-MB-231 cells blocked both their local tumour-forming and lung-colonising activity in immunodeficient mice. Conversely, upregulated expression of YB-1 enhanced the poor in vivo tumorigenicity of T47D cells. We then found that YB-1 knockdown also inhibits the initial generation in mice of invasive ductal carcinomas and ductal carcinomas in situ from freshly isolated human mammary cells transduced, respectively, with KRAS^G12D or myristoylated-AKT1. Interestingly, increased expression of HIF1α and G3BP1, two YB-1 translational targets and elements of a stress-adaptive programme, mirrored the levels of YB-1 in both transformed primary and established MDA-MB-231 breast cancer cells.

IceR improves proteome coverage and data completeness in global and single-cell proteomics

Label-free proteomics by data-dependent acquisition enables the unbiased quantification of thousands of proteins, however it notoriously suffers from high rates of missing values, thus prohibiting consistent protein quantification across large sample cohorts. To solve this, we here present IceR (Ion current extraction Re-quantification), an efficient and user-friendly quantification workflow that combines high identification rates of data-dependent acquisition with low missing value rates similar to data-independent acquisition. Specifically, IceR uses ion current information for a hybrid peptide identification propagation approach with superior quantification precision, accuracy, reliability and data completeness compared to other quantitative workflows. Applied to plasma and single-cell proteomics data, IceR enhanced the number of reliably quantified proteins, improved discriminability between single-cell populations, and allowed reconstruction of a developmental trajectory. IceR will be useful to improve performance of large scale global as well as low-input proteomics applications, facilitated by its availability as an easy-to-use R-package.

Benchmarking Quantitative Performance in Label-Free Proteomics

Previous benchmarking studies have demonstrated the importance of instrument acquisition methodology and statistical analysis on quantitative performance in label-free proteomics. However, the effects of these parameters in combination with replicate number and false discovery rate (FDR) corrections are not known. Using a benchmarking standard, we systematically evaluated the combined impact of acquisition methodology, replicate number, statistical approach, and FDR corrections. These analyses reveal a complex interaction between these parameters that greatly impacts the quantitative fidelity of protein- and peptide-level quantification. At a high replicate number (n = 8), both data-dependent acquisition (DDA) and data-independent acquisition (DIA) methodologies yield accurate protein quantification across statistical approaches. However, at a low replicate number (n = 4), only DIA in combination with linear models for microarrays (LIMMA) and reproducibility-optimized test statistic (ROTS) produced a high level of quantitative fidelity. Quantitative accuracy at low replicates is also greatly impacted by FDR corrections, with Benjamini-Hochberg and Storey corrections yielding variable true positive rates for DDA workflows. For peptide quantification, replicate number and acquisition methodology are even more critical. A higher number of replicates in combination with DIA and LIMMA produce high quantitative fidelity, while DDA performs poorly regardless of replicate number or statistical approach. These results underscore the importance of pairing instrument acquisition methodology with the appropriate replicate number and statistical approach for optimal quantification performance.

Re-expression of SMARCA4/BRG1 in small cell carcinoma of ovary, hypercalcemic type (SCCOHT) promotes an epithelial-like gene signature through an AP-1-dependent mechanism

Small cell carcinoma of the ovary, hypercalcemic type (SCCOHT) is a rare and aggressive form of ovarian cancer. SCCOHT tumors have inactivating mutations in SMARCA4 (BRG1), one of the two mutually exclusive ATPases of the SWI/SNF chromatin remodeling complex. To address the role that BRG1 loss plays in SCCOHT tumorigenesis, we performed integrative multi-omic analyses in SCCOHT cell lines +/- BRG1 reexpression. BRG1 reexpression induced a gene and protein signature similar to an epithelial cell and gained chromatin accessibility sites correlated with other epithelial originating TCGA tumors. Gained chromatin accessibility and BRG1 recruited sites were strongly enriched for transcription-factor-binding motifs of AP-1 family members. Furthermore, AP-1 motifs were enriched at the promoters of highly upregulated epithelial genes. Using a dominant-negative AP-1 cell line, we found that both AP-1 DNA-binding activity and BRG1 reexpression are necessary for the gene and protein expression of epithelial genes. Our study demonstrates that BRG1 reexpression drives an epithelial-like gene and protein signature in SCCOHT cells that depends upon by AP-1 activity.

Detecting differential protein abundance by combining peptide level P-values

The majority of methods for detecting differentially abundant proteins between samples in label-free LC-MS bottom-up proteomics experiments rely on statistically testing inferred protein abundances derived from peptide ionization intensities or averaging peptide level statistics. Here, we statistically test peptide ionization intensities directly and combine the resulting dependent P-values using the Empirical Brown's Method (EBM), avoiding error introduced through the estimation of protein abundances or summarizing test statistics. We show that on a spike-in proteomics dataset, a peptide level approach using EBM outperforms differential abundance detection using a protein level approach and several analysis workflows, including MSstats. Additionally, we demonstrate the effectiveness of this approach by detecting enriched proteins from an activity-based protein profiling dataset.

Whole-proteome analysis of mesonephric-derived cancers describes new potential biomarkers

Mesonephric carcinomas (MEs) and female adnexal tumors of probable Wolffian origin (FATWO) are derived from embryologic remnants of Wolffian/mesonephric ducts. Mesonephric-like carcinomas (MLCs) show identical morphology to ME of the cervix but occur in the uterus and ovary without convincing mesonephric remnants. ME, MLC, and FATWO are challenging to diagnose due to their morphologic similarities to Müllerian/paramesonephric tumors, contributing to a lack of evidence-based and tumor-specific treatments. We performed whole-proteomic analysis on 9 ME/MLC and 56 endometrial carcinomas (ECs) to identify potential diagnostic biomarkers. Although there were no convincing differences between ME and MLC, 543 proteins showed increased expression in ME/MLC relative to EC. From these proteins, euchromatic histone lysine methyltransferase 2 (EHMT2), glutathione S-transferase Mu 3 (GSTM3), eukaryotic translation elongation factor 1 alpha 2 (EEF1A2), and glycogen synthase kinase 3 beta were identified as putative biomarkers. Immunohistochemistry was performed on these candidates and GATA3 in 14 ME/MLC, 8 FATWO, 155 EC, and normal tissues. Of the candidates, only GATA3 and EHMT2 were highly expressed in mesonephric remnants and mesonephric-derived male tissues. GATA3 had the highest sensitivity and specificity for ME/MLC versus EC (93% and 99%) but was absent in FATWO. EHMT2 was 100% sensitive for ME/MLC & FATWO but was not specific (65%). Similarly, EEF1A2 was reasonably sensitive to ME/MLC (92%) and FATWO (88%) but was the least specific (38%). GSTM3 performed intermediately (sensitivity for ME/MLC and FATWO: 83% and 38%, respectively; specificity 67%). Although GATA3 remained the best diagnostic biomarker for ME/MLC, we have identified EHMT2, EEF1A2, and GSTM3 as proteins of interest in these cancers. FATWO's cell of origin is uncertain and remains an area for future research.

Novel protein pathways in development and progression of pulmonary sarcoidosis

Pulmonary involvement occurs in up to 95% of sarcoidosis cases. In this pilot study, we examine lung compartment-specific protein expression to identify pathways linked to development and progression of pulmonary sarcoidosis. We characterized bronchoalveolar lavage (BAL) cells and fluid (BALF) proteins in recently diagnosed sarcoidosis cases. We identified 4,306 proteins in BAL cells, of which 272 proteins were differentially expressed in sarcoidosis compared to controls. These proteins map to novel pathways such as integrin-linked kinase and IL-8 signaling and previously implicated pathways in sarcoidosis, including phagosome maturation, clathrin-mediated endocytic signaling and redox balance. In the BALF, the differentially expressed proteins map to several pathways identified in the BAL cells. The differentially expressed BALF proteins also map to aryl hydrocarbon signaling, communication between innate and adaptive immune response, integrin, PTEN and phospholipase C signaling, serotonin and tryptophan metabolism, autophagy, and B cell receptor signaling. Additional pathways that were different between progressive and non-progressive sarcoidosis in the BALF included CD28 signaling and PFKFB4 signaling. Our studies demonstrate the power of contemporary proteomics to reveal novel mechanisms operational in sarcoidosis. Application of our workflows in well-phenotyped large cohorts maybe beneficial to identify biomarkers for diagnosis and prognosis and therapeutically tenable molecular mechanisms.

Precursor Intensity-Based Label-Free Quantification Software Tools for Proteomic and Multi-Omic Analysis within the Galaxy Platform

For mass spectrometry-based peptide and protein quantification, label-free quantification (LFQ) based on precursor mass peak (MS1) intensities is considered reliable due to its dynamic range, reproducibility, and accuracy. LFQ enables peptide-level quantitation, which is useful in proteomics (analyzing peptides carrying post-translational modifications) and multi-omics studies such as metaproteomics (analyzing taxon-specific microbial peptides) and proteogenomics (analyzing non-canonical sequences). Bioinformatics workflows accessible via the Galaxy platform have proven useful for analysis of such complex multi-omic studies. However, workflows within the Galaxy platform have lacked well-tested LFQ tools. In this study, we have evaluated moFF and FlashLFQ, two open-source LFQ tools, and implemented them within the Galaxy platform to offer access and use via established workflows. Through rigorous testing and communication with the tool developers, we have optimized the performance of each tool. Software features evaluated include: (a) match-between-runs (MBR); (b) using multiple file-formats as input for improved quantification; (c) use of containers and/or conda packages; (d) parameters needed for analyzing large datasets; and (e) optimization and validation of software performance. This work establishes a process for software implementation, optimization, and validation, and offers access to two robust software tools for LFQ-based analysis within the Galaxy platform.

A Bayesian Null Interval Hypothesis Test Controls False Discovery Rates and Improves Sensitivity in Label-Free Quantitative Proteomics

Statistical significance tests are a common feature in quantitative proteomics workflows. The Student's t-test is widely used to compute the statistical significance of a protein's change between two groups of samples. However, the t-test's null hypothesis asserts that the difference in means between two groups is exactly zero, often marking small but uninteresting fold-changes as statistically significant. Compensations to address this issue are widely used in quantitative proteomics, but we suggest that a replacement of the t-test with a Bayesian approach offers a better path forward. In this article, we describe a Bayesian hypothesis test in which the null hypothesis is an interval rather than a single point at zero; the width of the interval is estimated from population statistics. The improved sensitivity of the method substantially increases the number of truly changing proteins detected in two benchmark data sets (ProteomeXchange identifiers PXD005590 and PXD016470). The method has been implemented within FlashLFQ, an open-source software program that quantifies bottom-up proteomics search results obtained from any search tool. FlashLFQ is rapid, sensitive, and accurate and is available both as an easy-to-use graphical user interface (Windows) and as a command-line tool (Windows/Linux/OSX).

MS1 ion current-based quantitative proteomics: A promising solution for reliable analysis of large biological cohorts

The rapidly-advancing field of pharmaceutical and clinical research calls for systematic, molecular-level characterization of complex biological systems. To this end, quantitative proteomics represents a powerful tool but an optimal solution for reliable large-cohort proteomics analysis, as frequently involved in pharmaceutical/clinical investigations, is urgently needed. Large-cohort analysis remains challenging owing to the deteriorating quantitative quality and snowballing missing data and false-positive discovery of altered proteins when sample size increases. MS1 ion current-based methods, which have become an important class of label-free quantification techniques during the past decade, show considerable potential to achieve reproducible protein measurements in large cohorts with high quantitative accuracy/precision. Nonetheless, in order to fully unleash this potential, several critical prerequisites should be met. Here we provide an overview of the rationale of MS1-based strategies and then important considerations for experimental and data processing techniques, with the emphasis on (i) efficient and reproducible sample preparation and LC separation; (ii) sensitive, selective and high-resolution MS detection; iii)accurate chromatographic alignment; (iv) sensitive and selective generation of quantitative features; and (v) optimal post-feature-generation data quality control. Prominent technical developments in these aspects are discussed. Finally, we reviewed applications of MS1-based strategy in disease mechanism studies, biomarker discovery, and pharmaceutical investigations.

A Data Analysis Protocol for Quantitative Data-Independent Acquisition Proteomics

Data-independent acquisition (DIA) mode of mass spectrometry, such as the SWATH-MS technology, enables accurate and consistent measurement of proteins, which is crucial for comparative proteomics studies. However, there is lack of free and easy to implement data analysis protocols that can handle the different data processing steps from raw spectrum files to peptide intensity matrix and its downstream analysis. Here, we provide a data analysis protocol, named diatools, covering all these steps from spectral library building to differential expression analysis of DIA proteomics data. The data analysis tools used in this protocol are open source and the protocol is distributed at Docker Hub as a complete software environment that supports Linux, Windows, and macOS operating systems.

Robust determination of differential abundance in shotgun proteomics using nonparametric statistics

Label-free shotgun mass spectrometry enables the detection of significant changes in protein abundance between different conditions. Due to often limited cohort sizes or replication, large ratios of potential protein markers to number of samples, as well as multiple null measurements pose important technical challenges to conventional parametric models. From a statistical perspective, a scenario similar to that of unlabeled proteomics is encountered in genomics when looking for differentially expressed genes. Still, the difficulty of detecting a large fraction of the true positives without a high false discovery rate is arguably greater in proteomics due to even smaller sample sizes and peptide-to-peptide variability in detectability. These constraints argue for nonparametric (or distribution-free) tests on normalized peptide values, thus minimizing the number of free parameters, as well as for measuring significance with permutation testing. We propose such a procedure with a class-based statistic, no parametric assumptions, and no parameters to select other than a nominal false discovery rate. Our method was tested on a new dataset which is available via ProteomeXchange with identifier PXD006447. The dataset was prepared using a standard proteolytic digest of a human protein mixture at 1.5-fold to 3-fold protein concentration changes and diluted into a constant background of yeast proteins. We demonstrate its superiority relative to other approaches in terms of the realized sensitivity and realized false discovery rates determined by ground truth, and recommend it for detecting differentially abundant proteins from MS data.

Histone Deacetylase Inhibitors Synergize with Catalytic Inhibitors of EZH2 to Exhibit Antitumor Activity in Small Cell Carcinoma of the Ovary, Hypercalcemic Type

Small cell carcinoma of the ovary, hypercalcemic type (SCCOHT) is a rare but extremely lethal malignancy that mainly impacts young women. SCCOHT is characterized by a diploid genome with loss of SMARCA4 and lack of SMARCA2 expression, two mutually exclusive ATPases of the SWI/SNF chromatin-remodeling complex. We and others have identified the histone methyltransferase EZH2 as a promising therapeutic target for SCCOHT, suggesting that SCCOHT cells depend on the alternation of epigenetic pathways for survival. In this study, we found that SCCOHT cells were more sensitive to pan-HDAC inhibitors compared with other ovarian cancer lines or immortalized cell lines tested. Pan-HDAC inhibitors, such as quisinostat, reversed the expression of a group of proteins that were deregulated in SCCOHT cells due to SMARCA4 loss, leading to growth arrest, apoptosis, and differentiation in vitro and suppressed tumor growth of xenografted tumors of SCCOHT cells. Moreover, combined treatment of HDAC inhibitors and EZH2 inhibitors at sublethal doses synergistically induced histone H3K27 acetylation and target gene expression, leading to rapid induction of apoptosis and growth suppression of SCCOHT cells and xenografted tumors. Therefore, our preclinical study highlighted the therapeutic potential of combined treatment of HDAC inhibitors with EZH2 catalytic inhibitors to treat SCCOHT.

Selective aggregation of the splicing factor Hsh155 suppresses splicing upon genotoxic stress

Upon genotoxic stress, dynamic relocalization events control DNA repair as well as alterations of the transcriptome and proteome, enabling stress recovery. How these events may influence one another is only partly known. Beginning with a cytological screen of genome stability proteins, we find that the splicing factor Hsh155 disassembles from its partners and localizes to both intranuclear and cytoplasmic protein quality control (PQC) aggregates under alkylation stress. Aggregate sequestration of Hsh155 occurs at nuclear and then cytoplasmic sites in a manner that is regulated by molecular chaperones and requires TORC1 activity signaling through the Sfp1 transcription factor. This dynamic behavior is associated with intron retention in ribosomal protein gene transcripts, a decrease in splicing efficiency, and more rapid recovery from stress. Collectively, our analyses suggest a model in which some proteins evicted from chromatin and undergoing transcriptional remodeling during stress are targeted to PQC sites to influence gene expression changes and facilitate stress recovery.

Enhanced differential expression statistics for data-independent acquisition proteomics

We describe a new reproducibility-optimization method ROPECA for statistical analysis of proteomics data with a specific focus on the emerging data-independent acquisition (DIA) mass spectrometry technology. ROPECA optimizes the reproducibility of statistical testing on peptide-level and aggregates the peptide-level changes to determine differential protein-level expression. Using a ‘gold standard’ spike-in data and a hybrid proteome benchmark data we show the competitive performance of ROPECA over conventional protein-based analysis as well as state-of-the-art peptide-based tools especially in DIA data with consistent peptide measurements. Furthermore, we also demonstrate the improved accuracy of our method in clinical studies using proteomics data from a longitudinal human twin study.

HIF prolyl hydroxylase PHD3 regulates translational machinery and glucose metabolism in clear cell renal cell carcinoma

BACKGROUND

A key feature of clear cell renal cell carcinoma (ccRCC) is the inactivation of the von Hippel-Lindau tumour suppressor protein (pVHL) that leads to the activation of hypoxia-inducible factor (HIF) pathway also in well-oxygenated conditions. Important regulator of HIF-α, prolyl hydroxylase PHD3, is expressed in high amounts in ccRCC. Although several functions and downstream targets for PHD3 in cancer have been suggested, the role of elevated PHD3 expression in ccRCC is not clear.

METHODS

To gain insight into the functions of high PHD3 expression in ccRCC, we used PHD3 knockdown by siRNA in 786-O cells under normoxic and hypoxic conditions and performed discovery mass spectrometry (LC-MS/MS) of the purified peptide samples. The LC-MS/MS results were analysed by label-free quantification of proteome data using a peptide-level expression-change averaging procedure and subsequent gene ontology enrichment analysis.

RESULTS

Our data reveals an intriguingly widespread effect of PHD3 knockdown with 91 significantly regulated proteins. Under hypoxia, the response to PHD3 silencing was wider than under normoxia illustrated by both the number of regulated proteins and by the range of protein expression levels. The main cellular functions regulated by PHD3 expression were glucose metabolism, protein translation and messenger RNA (mRNA) processing. PHD3 silencing led to downregulation of most glycolytic enzymes from glucose transport to lactate production supported by the reduction in extracellular acidification and lactate production and increase in cellular oxygen consumption rate. Moreover, upregulation of mRNA processing-related proteins and alteration in a number of ribosomal proteins was seen as a response to PHD3 silencing. Further studies on upstream effectors of the translational machinery revealed a possible role for PHD3 in regulation of mTOR pathway signalling.

CONCLUSIONS

Our findings suggest crucial involvement of PHD3 in the maintenance of key cellular functions including glycolysis and protein synthesis in ccRCC.

The histone methyltransferase EZH2 is a therapeutic target in small cell carcinoma of the ovary, hypercalcaemic type

Small cell carcinoma of the ovary, hypercalcaemic type (SCCOHT) is a rare but aggressive and untreatable malignancy affecting young women. We and others recently discovered that SMARCA4, a gene encoding the ATPase of the SWI/SNF chromatin-remodelling complex, is the only gene recurrently mutated in the majority of SCCOHT. The low somatic complexity of SCCOHT genomes and the prominent role of the SWI/SNF chromatin-remodelling complex in transcriptional control of genes suggest that SCCOHT cells may rely on epigenetic rewiring for oncogenic transformation. Herein, we report that approximately 80% (19/24) of SCCOHT tumour samples have strong expression of the histone methyltransferase EZH2 by immunohistochemistry, with the rest expressing variable amounts of EZH2. Re-expression of SMARCA4 suppressed the expression of EZH2 in SCCOHT cells. In comparison to other ovarian cell lines, SCCOHT cells displayed hypersensitivity to EZH2 shRNAs and two selective EZH2 inhibitors, GSK126 and EPZ-6438. EZH2 inhibitors induced cell cycle arrest, apoptosis, and cell differentiation in SCCOHT cells, along with the induction of genes involved in cell cycle regulation, apoptosis, and neuron-like differentiation. EZH2 inhibitors suppressed tumour growth and improved the survival of mice bearing SCCOHT xenografts. Therefore, our data suggest that loss of SMARCA4 creates a dependency on the catalytic activity of EZH2 in SCCOHT cells and that pharmacological inhibition of EZH2 is a promising therapeutic strategy for treating this disease. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Covariation of Peptide Abundances Accurately Reflects Protein Concentration Differences

Most implementations of mass spectrometry-based proteomics involve enzymatic digestion of proteins, expanding the analysis to multiple proteolytic peptides for each protein. Currently, there is no consensus of how to summarize peptides' abundances to protein concentrations, and such efforts are complicated by the fact that error control normally is applied to the identification process, and do not directly control errors linking peptide abundance measures to protein concentration. Peptides resulting from suboptimal digestion or being partially modified are not representative of the protein concentration. Without a mechanism to remove such unrepresentative peptides, their abundance adversely impacts the estimation of their protein's concentration. Here, we present a relative quantification approach, Diffacto, that applies factor analysis to extract the covariation of peptides' abundances. The method enables a weighted geometrical average summarization and automatic elimination of incoherent peptides. We demonstrate, based on a set of controlled label-free experiments using standard mixtures of proteins, that the covariation structure extracted by the factor analysis accurately reflects protein concentrations. In the 1% peptide-spectrum match-level FDR data set, as many as 11% of the peptides have abundance differences incoherent with the other peptides attributed to the same protein. If not controlled, such contradicting peptide abundance have a severe impact on protein quantifications. When adding the quantities of each protein's three most abundant peptides, we note as many as 14% of the proteins being estimated as having a negative correlation with their actual concentration differences between samples. Diffacto reduced the amount of such obviously incorrectly quantified proteins to 1.6%. Furthermore, by analyzing clinical data sets from two breast cancer studies, our method revealed the persistent proteomic signatures linked to three subtypes of breast cancer. We conclude that Diffacto can facilitate the interpretation and enhance the utility of most types of proteomics data.

Quantitative Profiling of Single Formalin Fixed Tumour Sections: proteomics for translational research

Although re-sequencing of gene panels and mRNA expression profiling are now firmly established in clinical laboratories, in-depth proteome analysis has remained a niche technology, better suited for studying model systems rather than challenging materials such as clinical trial samples. To address this limitation, we have developed a novel and optimized platform called SP3-Clinical Tissue Proteomics (SP3-CTP) for in-depth proteome profiling of practical quantities of tumour tissues, including formalin fixed and paraffin embedded (FFPE). Using single 10 μm scrolls of clinical tumour blocks, we performed in-depth quantitative analyses of individual sections from ovarian tumours covering the high-grade serous, clear cell, and endometrioid histotypes. This examination enabled the generation of a novel high-resolution proteome map of ovarian cancer histotypes from clinical tissues. Comparison of the obtained proteome data with large-scale genome and transcriptome analyses validated the observed proteome biology for previously validated hallmarks of this disease, and also identified novel protein features. A tissue microarray analysis validated cystathionine gamma-lyase (CTH) as a novel clear cell carcinoma feature with potential clinical relevance. In addition to providing a milestone in the understanding of ovarian cancer biology, these results show that in-depth proteomic analysis of clinically annotated FFPE materials can be effectively used as a biomarker discovery tool and perhaps ultimately as a diagnostic approach.

Oncogenic KRAS and BRAF Drive Metabolic Reprogramming in Colorectal Cancer

Metabolic reprogramming, in which altered utilization of glucose and glutamine supports rapid growth, is a hallmark of most cancers. Mutations in the oncogenes KRAS and BRAF drive metabolic reprogramming through enhanced glucose uptake, but the broader impact of these mutations on pathways of carbon metabolism is unknown. Global shotgun proteomic analysis of isogenic DLD-1 and RKO colon cancer cell lines expressing mutant and wild type KRAS or BRAF, respectively, failed to identify significant differences (at least 2-fold) in metabolic protein abundance. However, a multiplexed parallel reaction monitoring (PRM) strategy targeting 73 metabolic proteins identified significant protein abundance increases of 1.25-twofold in glycolysis, the nonoxidative pentose phosphate pathway, glutamine metabolism, and the phosphoserine biosynthetic pathway in cells with KRAS G13D mutations or BRAF V600E mutations. These alterations corresponded to mutant KRAS and BRAF-dependent increases in glucose uptake and lactate production. Metabolic reprogramming and glucose conversion to lactate in RKO cells were proportional to levels of BRAF V600E protein. In DLD-1 cells, these effects were independent of the ratio of KRAS G13D to KRAS wild type protein. A study of 8 KRAS wild type and 8 KRAS mutant human colon tumors confirmed the association of increased expression of glycolytic and glutamine metabolic proteins with KRAS mutant status. Metabolic reprogramming is driven largely by modest (<2-fold) alterations in protein expression, which are not readily detected by the global profiling methods most commonly employed in proteomic studies. The results indicate the superiority of more precise, multiplexed, pathway-targeted analyses to study functional proteome systems. Data are available through MassIVE Accession MSV000079486 at ftp://MSV000079486@massive.ucsd.edu.