Covariation of Peptide Abundances Accurately Reflects Protein Concentration Differences

Evaluating different methods for kidney recellularization

This study explores a potential solution to the shortage of kidneys for transplantation in end-stage renal disease (ESRD). Currently, kidney transplantation stands as the optimal option, yet the scarcity of organs persists. Employing tissue engineering, researchers sought to assess the feasibility of generating kidneys for transplantation. Pig kidneys were utilized since they possess higher similarities to human kidneys. Cells were removed via decellularization, which maintains the organ's microarchitecture. Subsequently, pig kidney cells and human red blood cells were perfused into the vacant kidney structure to reconstitute it. The methodologies employed showed promising results, suggesting a viable approach to increase the recellularization rate in whole pig kidneys. This proof-of-concept establishes a groundwork for potentially extending this technology to human kidneys, tackling the organ shortage, thus positively enhancing outcomes for ESRD patients by increasing the availability of transplantable organs.

On the excessive use of coefficient of variation as a metric of quantitation quality in proteomics

The coefficient of variation (CV) is often used in proteomics as a proxy to characterize the performance of a quantitation method and/or the related software. In this note, we question the excessive reliance on this metric in quantitative proteomics that may result in erroneous conclusions. We support this note using a ground-truth Human-Yeast-E. coli dataset demonstrating in a number of cases that erroneous data processing methods may lead to a low CV which has nothing to do with these methods' performances in quantitation.

Ultra-Fast Mass Spectrometry in Plant Biochemistry: Response of Winter Wheat Proteomics to Pre-Sowing Treatment with Iron Compounds

In recent years, ultrafast liquid chromatography/mass spectrometry methods have been extensively developed for the use in proteome profiling in biochemical studies. These methods are intended for express monitoring of cell response to biotic stimuli and elucidation of correlation of molecular changes with biological processes and phenotypical changes. New technologies, including the use of nanomaterials, are actively introduced to increase agricultural production. However, this requires complex approbation of new fertilizers and investigation of mechanisms underlying the biotic effects on the germination, growth, and development of plants. The aim of this work was to adapt the method of ultrafast chromatography/mass spectrometry for rapid quantitative profiling of molecular changes in 7-day-old wheat seedlings in response to pre-sowing seed treatment with iron compounds. The used method allows to analyze up to 200 samples per day; its practical value lies in the possibility of express proteomic diagnostics of the biotic action of new treatments, including those intended for agricultural needs. Changes in the regulation of photosynthesis, biosynthesis of chlorophyll and porphyrin- and tetrapyrrole-containing compounds, glycolysis (in shoot tissues), and polysaccharide metabolism (in root tissues) were shown after seed treatment with suspensions containing film-forming polymers (PEG 400, Na-CMC, Na2-EDTA), iron (II, III) nanoparticles, or iron (II) sulfate. Observations at the protein levels were consistent with the results of morphometry, superoxide dismutase activity assay, and microelement analysis of 3-day-old germinated seeds and shoots and roots of 7-day-old seedlings. A characteristic molecular signature involving proteins participating in the regulation of photosynthesis and glycolytic process was suggested as a potential marker of the biotic effects of seed treatment with iron compounds, which will be confirmed in further studies.

IQMMA: Efficient MS1 Intensity Extraction Pipeline Using Multiple Feature Detection Algorithms for DDA Proteomics

One of the key steps in data dependent acquisition (DDA) proteomics is detection of peptide isotopic clusters, also called "features", in MS1 spectra and matching them to MS/MS-based peptide identifications. A number of peptide feature detection tools became available in recent years, each relying on its own matching algorithm. Here, we provide an integrated solution, the intensity-based Quantitative Mix and Match Approach (IQMMA), which integrates a number of untargeted peptide feature detection algorithms and returns the most probable intensity values for the MS/MS-based identifications. IQMMA was tested using available proteomic data acquired for both well-characterized (ground truth) and real-world biological samples, including a mix of Yeast and E. coli digests spiked at different concentrations into the Human K562 digest used as a background, and a set of glioblastoma cell lines. Three open-source feature detection algorithms were integrated: Dinosaur, biosaur2, and OpenMS FeatureFinder. None of them was found optimal when applied individually to all the data sets employed in this work; however, their combined use in IQMMA improved efficiency of subsequent protein quantitation. The software implementing IQMMA is freely available at https://github.com/PostoenkoVI/IQMMA under Apache 2.0 license.

Accurate Label-Free Quantification by directLFQ to Compare Unlimited Numbers of Proteomes

Recent advances in mass spectrometry-based proteomics enable the acquisition of increasingly large datasets within relatively short times, which exposes bottlenecks in the bioinformatics pipeline. Although peptide identification is already scalable, most label-free quantification (LFQ) algorithms scale quadratic or cubic with the sample numbers, which may even preclude the analysis of large-scale data. Here we introduce directLFQ, a ratio-based approach for sample normalization and the calculation of protein intensities. It estimates quantities via aligning samples and ion traces by shifting them on top of each other in logarithmic space. Importantly, directLFQ scales linearly with the number of samples, allowing analyses of large studies to finish in minutes instead of days or months. We quantify 10,000 proteomes in 10 min and 100,000 proteomes in less than 2 h, a 1000-fold faster than some implementations of the popular LFQ algorithm MaxLFQ. In-depth characterization of directLFQ reveals excellent normalization properties and benchmark results, comparing favorably to MaxLFQ for both data-dependent acquisition and data-independent acquisition. In addition, directLFQ provides normalized peptide intensity estimates for peptide-level comparisons. It is an important part of an overall quantitative proteomic pipeline that also needs to include high sensitive statistical analysis leading to proteoform resolution. Available as an open-source Python package and a graphical user interface with a one-click installer, it can be used in the AlphaPept ecosystem as well as downstream of most common computational proteomics pipelines.

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

Yeast Ribonucleotide Reductase Is a Direct Target of the Proteasome and Provides Hyper Resistance to the Carcinogen 4-NQO

Various external and internal factors damaging DNA constantly disrupt the stability of the genome. Cells use numerous dedicated DNA repair systems to detect damage and restore genomic integrity in a timely manner. Ribonucleotide reductase (RNR) is a key enzyme providing dNTPs for DNA repair. Molecular mechanisms of indirect regulation of yeast RNR activity are well understood, whereas little is known about its direct regulation. The study was aimed at elucidation of the proteasome-dependent mechanism of direct regulation of RNR subunits in Saccharomyces cerevisiae. Proteome analysis followed by Western blot, RT-PCR, and yeast plating analysis showed that upregulation of RNR by proteasome deregulation is associated with yeast hyper resistance to 4-nitroquinoline-1-oxide (4-NQO), a UV-mimetic DNA-damaging drug used in animal models to study oncogenesis. Inhibition of RNR or deletion of RNR regulatory proteins reverses the phenotype of yeast hyper resistance to 4-NQO. We have shown for the first time that the yeast Rnr1 subunit is a substrate of the proteasome, which suggests a common mechanism of RNR regulation in yeast and mammals.

Proteomics-based scoring of cellular response to stimuli for improved characterization of signaling pathway activity

Omics technologies focus on uncovering the complex nature of molecular mechanisms in cells and organisms, including biomarkers and drug targets discovery. Aiming at these tasks, we see that information extracted from omics data is still underused. In particular, characteristics of differentially regulated molecules can be combined in a single score to quantify the signaling pathway activity. Such a metric can be useful for comprehensive data interpretation to follow: (1) developing molecular responses in time; (2) potency of a drug on a certain cell culture; (3) ranking the signaling pathway activity in stimulated cells; and (4) integration of the omics data and assay-based measurements of cell viability, cytotoxicity, and proliferation. With recent advances in ultrafast mass spectrometry for quantitative proteomics allowing data collection in a few minutes, proteomics score for cellular response to stimuli can become a fast, accurate, and informative complement to bioassays. Here, we utilized an interquartile-based selection of differentially regulated features and a variety of schemes for quantifying cellular responses to come up with the quantitative metric for total cellular response and pathway activity. Validation was performed using antiproliferative and virus assays and label-free proteomics data collected for cancer cells subjected to drug stimulation.

On the Feasibility of Using an Ultra-Fast DirectMS1 Method of Proteome-Wide Analysis for Searching Drug Targets in Chemical Proteomics

Protein quantitation in tissue cells or physiological fluids based on liquid chromatography/mass spectrometry is one of the key sources of information on the mechanisms of cell functioning during chemotherapeutic treatment. Information on significant changes in protein expression upon treatment can be obtained by chemical proteomics and requires analysis of the cellular proteomes, as well as development of experimental and bioinformatic methods for identification of the drug targets. Low throughput of whole proteome analysis based on liquid chromatography and tandem mass spectrometry is one of the main factors limiting the scale of these studies. The method of direct mass spectrometric identification of proteins, DirectMS1, is one of the approaches developed in recent years allowing ultrafast proteome-wide analyses employing minute-scale gradients for separation of proteolytic mixtures. Aim of this work was evaluation of both possibilities and limitations of the method for identification of drug targets at the level of whole proteome and for revealing cellular processes activated by the treatment. Particularly, the available literature data on chemical proteomics obtained earlier for a large set of onco-pharmaceuticals using multiplex quantitative proteome profiling were analyzed. The results obtained were further compared with the proteome-wide data acquired by the DirectMS1 method using ultrashort separation gradients to evaluate efficiency of the method in identifying known drug targets. Using ovarian cancer cell line A2780 as an example, a whole-proteome comparison of two cell lysis techniques was performed, including the freeze-thaw lysis commonly employed in chemical proteomics and the one based on ultrasonication for cell disruption, which is the widely accepted as a standard in proteomic studies. Also, the proteome-wide profiling was performed using ultrafast DirectMS1 method for A2780 cell line treated with lonidamine, followed by gene ontology analyses to evaluate capabilities of the method in revealing regulation of proteins in the cellular processes associated with drug treatment.

DirectMS1Quant: Ultrafast Quantitative Proteomics with MS/MS-Free Mass Spectrometry

Recently, we presented the DirectMS1 method of ultrafast proteome-wide analysis based on minute-long LC gradients and MS1-only mass spectra acquisition. Currently, the method provides the depth of human cell proteome coverage of 2500 proteins at a 1% false discovery rate (FDR) when using 5 min LC gradients and 7.3 min runtime in total. While the standard MS/MS approaches provide 4000-5000 protein identifications within a couple of hours of instrumentation time, we advocate here that the higher number of identified proteins does not always translate into better quantitation quality of the proteome analysis. To further elaborate on this issue, we performed a one-on-one comparison of quantitation results obtained using DirectMS1 with three popular MS/MS-based quantitation methods: label-free (LFQ) and tandem mass tag quantitation (TMT), both based on data-dependent acquisition (DDA) and data-independent acquisition (DIA). For comparison, we performed a series of proteome-wide analyses of well-characterized (ground truth) and biologically relevant samples, including a mix of UPS1 proteins spiked at different concentrations into an Echerichia coli digest used as a background and a set of glioblastoma cell lines. MS1-only data was analyzed using a novel quantitation workflow called DirectMS1Quant developed in this work. The results obtained in this study demonstrated comparable quantitation efficiency of 5 min DirectMS1 with both TMT and DIA methods, yet the latter two utilized a 10-20-fold longer instrumentation time.

Mass Spectrometry-Based Chemical Proteomics for Drug Target Discoveries

Chemical proteomics, emerging rapidly in recent years, has become a main approach to identifying interactions between the small molecules and proteins in the cells on a proteome scale and mapping the signaling and/or metabolic pathways activated and regulated by these interactions. The methods of chemical proteomics allow not only identifying proteins targeted by drugs, characterizing their toxicity and discovering possible off-target proteins, but also elucidation of the fundamental mechanisms of cell functioning under conditions of drug exposure or due to the changes in physiological state of the organism itself. Solving these problems is essential for both basic research in biology and clinical practice, including approaches to early diagnosis of various forms of serious diseases or prediction of the effectiveness of therapeutic treatment. At the same time, recent developments in high-resolution mass spectrometry have provided the technology for searching the drug targets across the whole cell proteomes. This review provides a concise description of the main objectives and problems of mass spectrometry-based chemical proteomics, the methods and approaches to their solution, and examples of implementation of these methods in biomedical research.

VIQoR: a web service for visually supervised protein inference and protein quantification

MOTIVATION

In quantitative bottom-up mass spectrometry (MS)-based proteomics, the reliable estimation of protein concentration changes from peptide quantifications between different biological samples is essential. This estimation is not a single task but comprises the two processes of protein inference and protein abundance summarization. Furthermore, due to the high complexity of proteomics data and associated uncertainty about the performance of these processes, there is a demand for comprehensive visualization methods able to integrate protein with peptide quantitative data including their post-translational modifications. Hence, there is a lack of a suitable tool that provides post-identification quantitative analysis of proteins with simultaneous interactive visualization.

RESULTS

In this article, we present VIQoR, a user-friendly web service that accepts peptide quantitative data of both labeled and label-free experiments and accomplishes the crucial components protein inference and summarization and interactive visualization modules, including the novel VIQoR plot. We implemented two different parsimonious algorithms to solve the protein inference problem, while protein summarization is facilitated by a well-established factor analysis algorithm called fast-FARMS followed by a weighted average summarization function that minimizes the effect of missing values. In addition, summarization is optimized by the so-called Global Correlation Indicator (GCI). We test the tool on three publicly available ground truth datasets and demonstrate the ability of the protein inference algorithms to handle shared peptides. We furthermore show that GCI increases the accuracy of the quantitative analysis in datasets with replicated design.

AVAILABILITY AND IMPLEMENTATION

VIQoR is accessible at: http://computproteomics.bmb.sdu.dk/Apps/VIQoR/. The source code is available at: https://bitbucket.org/veitveit/viqor/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Multiomic Profiling Identified EGF Receptor Signaling as a Potential Inhibitor of Type I Interferon Response in Models of Oncolytic Therapy by Vesicular Stomatitis Virus

Cancer cell lines responded differentially to type I interferon treatment in models of oncolytic therapy using vesicular stomatitis virus (VSV). Two opposite cases were considered in this study, glioblastoma DBTRG-05MG and osteosarcoma HOS cell lines exhibiting resistance and sensitivity to VSV after the treatment, respectively. Type I interferon responses were compared for these cell lines by integrative analysis of the transcriptome, proteome, and RNA editome to identify molecular factors determining differential effects observed. Adenosine-to-inosine RNA editing was equally induced in both cell lines. However, transcriptome analysis showed that the number of differentially expressed genes was much higher in DBTRG-05MG with a specific enrichment in inflammatory proteins. Further, it was found that two genes, EGFR and HER2, were overexpressed in HOS cells compared with DBTRG-05MG, supporting recent reports that EGF receptor signaling attenuates interferon responses via HER2 co-receptor activity. Accordingly, combined treatment of cells with EGF receptor inhibitors such as gefitinib and type I interferon increases the resistance of sensitive cell lines to VSV. Moreover, sensitive cell lines had increased levels of HER2 protein compared with non-sensitive DBTRG-05MG. Presumably, the level of this protein expression in tumor cells might be a predictive biomarker of their resistance to oncolytic viral therapy.

Putting Humpty Dumpty Back Together Again: What Does Protein Quantification Mean in Bottom-Up Proteomics?

Bottom-up proteomics provides peptide measurements and has been invaluable for moving proteomics into large-scale analyses. Commonly, a single quantitative value is reported for each protein-coding gene by aggregating peptide quantities into protein groups following protein inference or parsimony. However, given the complexity of both RNA splicing and post-translational protein modification, it is overly simplistic to assume that all peptides that map to a singular protein-coding gene will demonstrate the same quantitative response. By assuming that all peptides from a protein-coding sequence are representative of the same protein, we may miss the discovery of important biological differences. To capture the contributions of existing proteoforms, we need to reconsider the practice of aggregating protein values to a single quantity per protein-coding gene.

Data-Independent Acquisition Protease-Multiplexing Enables Increased Proteome Sequence Coverage Across Multiple Fragmentation Modes

The use of multiple proteases has been shown to increase protein sequence coverage in proteomics experiments; however, due to the additional analysis time required, it has not been widely adopted in routine data-dependent acquisition (DDA) proteomic workflows. Alternatively, data-independent acquisition (DIA) has the potential to analyze multiplexed samples from different protease digests, but has been primarily optimized for fragmenting tryptic peptides. Here we evaluate a DIA multiplexing approach that combines three proteolytic digests (Trypsin, AspN, and GluC) into a single sample. We first optimize data acquisition conditions for each protease individually with both the canonical DIA fragmentation mode (beam type CID), as well as resonance excitation CID, to determine optimal consensus conditions across proteases. Next, we demonstrate that application of these conditions to a protease-multiplexed sample of human peptides results in similar protein identifications and quantitative performance as compared to trypsin alone, but enables up to a 63% increase in peptide detections, and a 45% increase in nonredundant amino acid detections. Nontryptic peptides enabled noncanonical protein isoform determination and resulted in 100% sequence coverage for numerous proteins, suggesting the utility of this approach in applications where sequence coverage is critical, such as protein isoform analysis.

Multi-Omics Analysis of Glioblastoma Cells' Sensitivity to Oncolytic Viruses

Simple Summary

This study aims to uncover the contribution of interferon-dependent antiviral mechanisms preserved in tumor cells to the resistance of glioblastoma multiforme cells to oncolytic viruses. To characterize the functionality of interferon signaling, we used omics profiling and titration-based measurements of cell sensitivity to a panel of viruses of diverse oncolytic potential. This study shows why patient-derived glioblastoma cultures can acquire increased resistance to oncolytic viruses in the presence of interferons and suggests an approach to ranking glioblastoma cells by the acquired resistance. Our findings are important for monitoring the oncolytic potential of viruses to overcome IFN-induced resistance of tumor cells and contribute to successful therapy.

Abstract

Oncolytic viruses have gained momentum in the last decades as a promising tool for cancer treatment. Despite the progress, only a fraction of patients show a positive response to viral therapy. One of the key variable factors contributing to therapy outcomes is interferon-dependent antiviral mechanisms in tumor cells. Here, we evaluated this factor using patient-derived glioblastoma multiforme (GBM) cultures. Cell response to the type I interferons’ (IFNs) stimulation was characterized at mRNA and protein levels. Omics analysis revealed that GBM cells overexpress interferon-stimulated genes (ISGs) and upregulate their proteins, similar to the normal cells. A conserved molecular pattern unambiguously differentiates between the preserved and defective responses. Comparing ISGs’ portraits with titration-based measurements of cell sensitivity to a panel of viruses, the “strength” of IFN-induced resistance acquired by GBM cells was ranked. The study demonstrates that suppressing a single ISG and encoding an essential antiviral protein, does not necessarily increase sensitivity to viruses. Conversely, silencing IFIT3 and PLSCR1 genes in tumor cells can negatively affect the internalization of vesicular stomatitis and Newcastle disease viruses. We present evidence of a complex relationship between the interferon response genes and other factors affecting the sensitivity of tumor cells to viruses.

Ultra-fast proteomics with Scanning SWATH

Accurate quantification of the proteome remains challenging for large sample series and longitudinal experiments. We report a data-independent acquisition method, Scanning SWATH, that accelerates mass spectrometric (MS) duty cycles, yielding quantitative proteomes in combination with short gradients and high-flow (800 µl min-1) chromatography. Exploiting a continuous movement of the precursor isolation window to assign precursor masses to tandem mass spectrometry (MS/MS) fragment traces, Scanning SWATH increases precursor identifications by ~70% compared to conventional data-independent acquisition (DIA) methods on 0.5-5-min chromatographic gradients. We demonstrate the application of ultra-fast proteomics in drug mode-of-action screening and plasma proteomics. Scanning SWATH proteomes capture the mode of action of fungistatic azoles and statins. Moreover, we confirm 43 and identify 11 new plasma proteome biomarkers of COVID-19 severity, advancing patient classification and biomarker discovery. Thus, our results demonstrate a substantial acceleration and increased depth in fast proteomic experiments that facilitate proteomic drug screens and clinical studies.

Systematic detection of functional proteoform groups from bottom-up proteomic datasets

To a large extent functional diversity in cells is achieved by the expansion of molecular complexity beyond that of the coding genome. Various processes create multiple distinct but related proteins per coding gene - so-called proteoforms - that expand the functional capacity of a cell. Evaluating proteoforms from classical bottom-up proteomics datasets, where peptides instead of intact proteoforms are measured, has remained difficult. Here we present COPF, a tool for COrrelation-based functional ProteoForm assessment in bottom-up proteomics data. It leverages the concept of peptide correlation analysis to systematically assign peptides to co-varying proteoform groups. We show applications of COPF to protein complex co-fractionation data as well as to more typical protein abundance vs. sample data matrices, demonstrating the systematic detection of assembly- and tissue-specific proteoform groups, respectively, in either dataset. We envision that the presented approach lays the foundation for a systematic assessment of proteoforms and their functional implications directly from bottom-up proteomic datasets.

Proteomics of Cellular Response to Stress: Taking Control of False Positive Results

One of the main goals of quantitative proteomics is molecular profiling of cellular response to stress at the protein level. To perform this profiling, statistical analysis of experimental data involves multiple testing of a hypothesis about the equality of protein concentrations between the cells under normal and stress conditions. This analysis is then associated with the multiple testing problem dealing with the increased chance of obtaining false positive results. A number of solutions to this problem are known, yet, they may lead to the loss of potentially important biological information when applied with commonly accepted thresholds of statistical significance. Using the proteomic data obtained earlier for the yeast samples containing proteins at known concentrations and the biological models of early and late cellular responses to stress, we analyzed dependences of distributions of false positive and false negative rates on the protein fold changes and thresholds of statistical significance. Based on the analysis of the density of data points in the volcano plots, Benjamini-Hochberg method, and gene ontology analysis, visual approach for optimization of the statistical threshold and selection of the differentially regulated proteins has been suggested, which could be useful for researchers working in the field of quantitative proteomics.

A Tutorial for Variance-Sensitive Clustering and the Quantitative Analysis of Protein Complexes

Data clustering facilitates the identification of biologically relevant molecular features in quantitative proteomics experiments with thousands of measurements over multiple conditions. It finds groups of proteins or peptides with similar quantitative behavior across multiple experimental conditions. This co-regulatory behavior suggests that the proteins of such a group share their functional behavior and thus often can be mapped to the same biological processes and molecular subnetworks.While usual clustering approaches dismiss the variance of the measured proteins, VSClust combines statistical testing with pattern recognition into a common algorithm. Here, we show how to use the VSClust web service on a large proteomics data set and present further tools to assess the quantitative behavior of protein complexes.

Peptide Correlation Analysis (PeCorA) Reveals Differential Proteoform Regulation

Shotgun proteomics techniques infer the presence and quantity of proteins using peptide proxies produced by cleavage of the proteome with a protease. Most protein quantitation strategies assume that multiple peptides derived from a protein will behave quantitatively similar across treatment groups, but this assumption may be false due to (1) heterogeneous proteoforms and (2) technical artifacts. Here we describe a strategy called peptide correlation analysis (PeCorA) that detects quantitative disagreements between peptides mapped to the same protein. PeCorA fits linear models to assess whether a peptide's change across treatment groups differs from all other peptides assigned to the same protein. PeCorA revealed that ∼15% of proteins in a mouse microglia stress data set contain at least one discordant peptide. Inspection of the discordant peptides shows the utility of PeCorA for the direct and indirect detection of regulated post-translational modifications (PTMs) and also for the discovery of poorly quantified peptides. The exclusion of poorly quantified peptides before protein quantity summarization decreased false-positives in a benchmark data set. Finally, PeCorA suggests that the inactive isoform of prothrombin, a coagulation cascade protease, is more abundant in plasma from COVID-19 patients relative to non-COVID-19 controls. PeCorA is freely available as an R package that works with arbitrary tables of quantified peptides.

What Room for Two-Dimensional Gel-Based Proteomics in a Shotgun Proteomics World?

Two-dimensional gel electrophoresis was instrumental in the birth of proteomics in the late 1980s. However, it is now often considered as an outdated technique for proteomics—a thing of the past. Although this opinion may be true for some biological questions, e.g., when analysis depth is of critical importance, for many others, two-dimensional gel electrophoresis-based proteomics still has a lot to offer. This is because of its robustness, its ability to separate proteoforms, and its easy interface with many powerful biochemistry techniques (including western blotting). This paper reviews where and why two-dimensional gel electrophoresis-based proteomics can still be profitably used. It emerges that, rather than being a thing of the past, two-dimensional gel electrophoresis-based proteomics is still highly valuable for many studies. Thus, its use cannot be dismissed on simple fashion arguments and, as usual, in science, the tree is to be judged by the fruit.

Focus on the spectra that matter by clustering of quantification data in shotgun proteomics

In shotgun proteomics, the analysis of label-free quantification experiments is typically limited by the identification rate and the noise level in the quantitative data. This generally causes a low sensitivity in differential expression analysis. Here, we propose a quantification-first approach for peptides that reverses the classical identification-first workflow, thereby preventing valuable information from being discarded in the identification stage. Specifically, we introduce a method, Quandenser, that applies unsupervised clustering on both MS1 and MS2 level to summarize all analytes of interest without assigning identities. This reduces search time due to the data reduction. We can now employ open modification and de novo searches to identify analytes of interest that would have gone unnoticed in traditional pipelines. Quandenser+Triqler outperforms the state-of-the-art method MaxQuant+Perseus, consistently reporting more differentially abundant proteins for all tested datasets. Software is available for all major operating systems at https://github.com/statisticalbiotechnology/quandenser, under Apache 2.0 license.

Selection of Features with Consistent Profiles Improves Relative Protein Quantification in Mass Spectrometry Experiments

In bottom-up mass spectrometry-based proteomics, relative protein quantification is often achieved with data-dependent acquisition (DDA), data-independent acquisition (DIA), or selected reaction monitoring (SRM). These workflows quantify proteins by summarizing the abundances of all the spectral features of the protein (e.g. precursor ions, transitions or fragments) in a single value per protein per run. When abundances of some features are inconsistent with the overall protein profile (for technological reasons such as interferences, or for biological reasons such as post-translational modifications), the protein-level summaries and the downstream conclusions are undermined. We propose a statistical approach that automatically detects spectral features with such inconsistent patterns. The detected features can be separately investigated, and if necessary, removed from the data set. We evaluated the proposed approach on a series of benchmark-controlled mixtures and biological investigations with DDA, DIA and SRM data acquisitions. The results demonstrated that it could facilitate and complement manual curation of the data. Moreover, it can improve the estimation accuracy, sensitivity and specificity of detecting differentially abundant proteins, and reproducibility of conclusions across different data processing tools. The approach is implemented as an option in the open-source R-based software MSstats.

BayesENproteomics: Bayesian Elastic Nets for Quantification of Peptidoforms in Complex Samples

Multivariate regression modelling provides a statistically powerful means of quantifying the effects of a given treatment while compensating for sources of variation and noise, such as variability between human donors and the behavior of different peptides during mass spectrometry. However, methods to quantify endogenous post-translational modifications (PTMs) are typically reliant on summary statistical methods that fail to consider sources of variability such as changes in the levels of the parent protein. Here, we compare three multivariate regression methods, including a novel Bayesian elastic net algorithm (BayesENproteomics) that enables assessment of relative protein abundances while also quantifying identified PTMs for each protein. We tested the ability of these methods to accurately quantify expression of proteins in a mixed-species benchmark experiment and to quantify synthetic PTMs induced by stable isotope labelling. Finally, we extended our regression pipeline to calculate fold changes at the pathway level, providing a complement to commonly used enrichment analysis. Our results show that BayesENproteomics can quantify changes to protein levels across a broad dynamic range while also accurately quantifying PTM and pathway-level fold changes.

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

EPIFANY: A Method for Efficient High-Confidence Protein Inference

Accurate protein inference in the presence of shared peptides is still one of the key problems in bottom-up proteomics. Most protein inference tools employing simple heuristic inference strategies are efficient but exhibit reduced accuracy. More advanced probabilistic methods often exhibit better inference quality but tend to be too slow for large data sets. Here, we present a novel protein inference method, EPIFANY, combining a loopy belief propagation algorithm with convolution trees for efficient processing of Bayesian networks. We demonstrate that EPIFANY combines the reliable protein inference of Bayesian methods with significantly shorter runtimes. On the 2016 iPRG protein inference benchmark data, EPIFANY is the only tested method that finds all true-positive proteins at a 5% protein false discovery rate (FDR) without strict prefiltering on the peptide-spectrum match (PSM) level, yielding an increase in identification performance (+10% in the number of true positives and +14% in partial AUC) compared to previous approaches. Even very large data sets with hundreds of thousands of spectra (which are intractable with other Bayesian and some non-Bayesian tools) can be processed with EPIFANY within minutes. The increased inference quality including shared peptides results in better protein inference results and thus increased robustness of the biological hypotheses generated. EPIFANY is available as open-source software for all major platforms at https://OpenMS.de/epifany.

ComplexBrowser: A Tool for Identification and Quantification of Protein Complexes in Large-scale Proteomics Datasets

We have developed ComplexBrowser, an open source, online platform for supervised analysis of quantitative proteomic data (label free and isobaric mass tag based) that focuses on protein complexes. The software uses manually curated information from CORUM and Complex Portal databases to identify protein complex components. For the first time, we provide a Complex Fold Change (CFC) factor that identifies up- and downregulated complexes based on the level of complex subunits coregulation. The software provides interactive visualization of protein complexes' composition and expression for exploratory analysis and incorporates a quality control step that includes normalization and statistical analysis based on the limma package. ComplexBrowser was tested on two published studies identifying changes in protein expression within either human adenocarcinoma tissue or activated mouse T-cells. The analysis revealed 1519 and 332 protein complexes, of which 233 and 41 were found coordinately regulated in the respective studies. The adopted approach provided evidence for a shift to glucose-based metabolism and high proliferation in adenocarcinoma tissues, and the identification of chromatin remodeling complexes involved in mouse T-cell activation. The results correlate with the original interpretation of the experiments and provide novel biological details about the protein complexes affected. ComplexBrowser is, to our knowledge, the first tool to automate quantitative protein complex analysis for high-throughput studies, providing insights into protein complex regulation within minutes of analysis.

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.

Barocycler-Based Concurrent Multiomics Method To Assess Molecular Changes Associated with Atherosclerosis Using Small Amounts of Arterial Tissue from a Single Mouse

Atherosclerosis is a complex, multifactorial disease characterized by the buildup of plaque in the arterial wall. Apolipoprotein E gene deficient (Apoe^-/-) mice serve as a commonly used tool to elucidate the pathophysiology of atherosclerosis because of their propensity to spontaneously develop arterial lesions. To date, however, an integrated omics assessment of atherosclerotic lesions in individual Apoe^-/- mice has been challenging because of the small amount of diseased and nondiseased tissue available. To address this current limitation, we developed a multiomics method (Multi-ABLE) based on the proteomic method called accelerated Barocycler lysis and extraction (ABLE) to assess the depth of information that can be obtained from arterial tissue derived from a single mouse by splitting ABLE to allow for a combined proteomics-metabolomics-lipidomics analysis (Multi-ABLE). The new method includes tissue lysis via pressure cycling technology (PCT) in a Barocycler, followed by proteomic analysis of half the sample by nanoLC-MS and sequential extraction of lipids (organic extract) and metabolites (aqueous extract) combined with HILIC and reversed phase chromatography and time-of-flight mass spectrometry on the other half. Proteomic analysis identified 845 proteins, 93 of which were significantly altered in lesion-containing arteries. Lipidomic and metabolomic analyses detected 851 lipid and 362 metabolite features, which included 215 and 65 identified lipids and metabolites, respectively. The Multi-ABLE method is the first to apply a concurrent multiomics pipeline to cardiovascular disease using small (<5 mg) tissue samples, and it is applicable to other diseases where limited size samples are available at specific points during disease progression.

Bayesian Confidence Intervals for Multiplexed Proteomics Integrate Ion-statistics with Peptide Quantification Concordance

Multiplexed proteomics has emerged as a powerful tool to measure relative protein expression levels across multiple conditions. The relative protein abundances are inferred by comparing the signals generated by isobaric tags, which encode the samples' origins. Intuitively, the trust associated with a protein measurement depends on the similarity of ratios from the protein's peptides and the signal-strength of these measurements. However, typically the average peptide ratio is reported as the estimate of relative protein abundance, which is only the most likely ratio with a very naive model. Moreover, there is no sense on the confidence in these measurements. Here, we present a mathematically rigorous approach that integrates peptide signal strengths and peptide-measurement agreement into an estimation of the true protein ratio and the associated confidence (BACIQ). The main advantages of BACIQ are: (1) It removes the need to threshold reported peptide signal based on an arbitrary cut-off, thereby reporting more measurements from a given experiment; (2) Confidence can be assigned without replicates; (3) For repeated experiments BACIQ provides confidence intervals for the union, not the intersection, of quantified proteins; (4) For repeated experiments, BACIQ confidence intervals are more predictive than confidence intervals based on protein measurement agreement. To demonstrate the power of BACIQ we reanalyzed previously published data on subcellular protein movement on treatment with an Exportin-1 inhibiting drug. We detect ∼2× more highly significant movers, down to subcellular localization changes of ∼1%. Thus, our method drastically increases the value obtainable from quantitative proteomics experiments, helping researchers to interpret their data and prioritize resources. To make our approach easily accessible we distribute it via a Python/Stan package.

Mass Spectrometry-based Absolute Quantification of 20S Proteasome Status for Controlled Ex-vivo Expansion of Human Adipose-derived Mesenchymal Stromal/Stem Cells

The proteasome controls a multitude of cellular processes through protein degradation and has been identified as a therapeutic target in oncology. However, our understanding of its function and the development of specific modulators are hampered by the lack of a straightforward method to determine the overall proteasome status in biological samples. Here, we present a method to determine the absolute quantity and stoichiometry of ubiquitous and tissue-specific human 20S proteasome subtypes based on a robust, absolute SILAC-based multiplexed LC-Selected Reaction Monitoring (SRM) quantitative mass spectrometry assay with high precision, accuracy, and sensitivity. The method was initially optimized and validated by comparison with a reference ELISA assay and by analyzing the dynamics of catalytic subunits in HeLa cells following IFNγ-treatment and in range of human tissues. It was then successfully applied to reveal IFNγ- and O2-dependent variations of proteasome status during primary culture of Adipose-derived-mesenchymal Stromal/Stem Cells (ADSCs). The results show the critical importance of controlling the culture conditions during cell expansion for future therapeutic use in humans. We hypothesize that a shift from the standard proteasome to the immunoproteasome could serve as a predictor of immunosuppressive and differentiation capacities of ADSCs and, consequently, that quality control should include proteasomal quantification in addition to examining other essential cell parameters. The method presented also provides a new powerful tool to conduct more individualized protocols in cancer or inflammatory diseases where selective inhibition of the immunoproteasome has been shown to reduce side effects.

Integrated Identification and Quantification Error Probabilities for Shotgun Proteomics

Protein quantification by label-free shotgun proteomics experiments is plagued by a multitude of error sources. Typical pipelines for identifying differential proteins use intermediate filters to control the error rate. However, they often ignore certain error sources and, moreover, regard filtered lists as completely correct in subsequent steps. These two indiscretions can easily lead to a loss of control of the false discovery rate (FDR). We propose a probabilistic graphical model, Triqler, that propagates error information through all steps, employing distributions in favor of point estimates, most notably for missing value imputation. The model outputs posterior probabilities for fold changes between treatment groups, highlighting uncertainty rather than hiding it. We analyzed 3 engineered data sets and achieved FDR control and high sensitivity, even for truly absent proteins. In a bladder cancer clinical data set we discovered 35 proteins at 5% FDR, whereas the original study discovered 1 and MaxQuant/Perseus 4 proteins at this threshold. Compellingly, these 35 proteins showed enrichment for functional annotation terms, whereas the top ranked proteins reported by MaxQuant/Perseus showed no enrichment. The model executes in minutes and is freely available at https://pypi.org/project/triqler/.

CoExpresso: assess the quantitative behavior of protein complexes in human cells

Background

Translational and post-translational control mechanisms in the cell result in widely observable differences between measured gene transcription and protein abundances. Herein, protein complexes are among the most tightly controlled entities by selective degradation of their individual proteins. They furthermore act as control hubs that regulate highly important processes in the cell and exhibit a high functional diversity due to their ability to change their composition and their structure. Better understanding and prediction of these functional states demands methods for the characterization of complex composition, behavior, and abundance across multiple cell states. Mass spectrometry provides an unbiased approach to directly determine protein abundances across different cell populations and thus to profile a comprehensive abundance map of proteins.

Results

We provide a tool to investigate the behavior of protein subunits in known complexes by comparing their abundance profiles across up to 140 cell types available in ProteomicsDB. Thorough assessment of different randomization methods and statistical scoring algorithms allows determining the significance of concurrent profiles within a complex, therefore providing insights into the conservation of their composition across human cell types as well as the identification of intrinsic structures in complex behavior to determine which proteins orchestrate complex function. This analysis can be extended to investigate common profiles within arbitrary protein groups. CoExpresso can be accessed through http://computproteomics.bmb.sdu.dk/Apps/CoExpresso.

Conclusions

With the CoExpresso web service, we offer a potent scoring scheme to assess proteins for their co-regulation and thereby offer insight into their potential for forming functional groups like protein complexes.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2573-8) contains supplementary material, which is available to authorized users.

Pyteomics 4.0: Five Years of Development of a Python Proteomics Framework

Many of the novel ideas that drive today's proteomic technologies are focused essentially on experimental or data-processing workflows. The latter are implemented and published in a number of ways, from custom scripts and programs, to projects built using general-purpose or specialized workflow engines; a large part of routine data processing is performed manually or with custom scripts that remain unpublished. Facilitating the development of reproducible data-processing workflows becomes essential for increasing the efficiency of proteomic research. To assist in overcoming the bioinformatics challenges in the daily practice of proteomic laboratories, 5 years ago we developed and announced Pyteomics, a freely available open-source library providing Python interfaces to proteomic data. We summarize the new functionality of Pyteomics developed during the time since its introduction.

Urinary afamin levels are associated with the progression of diabetic nephropathy

AIMS

In this study, we applied quantitative proteomic analysis to identify urinary proteins associated with diabetic nephropathy (DN).

METHODS

Two-dimensional image-converted analysis of liquid chromatography and mass spectrometry detected the proteins differentially excreted between normoalbuminuric and macroalbuminuric patients with type 2 diabetes mellitus (T2DM) (n = 6 each). Urinary levels of excreted proteins were measured by multiple reaction monitoring (MRM) analysis using an independent sample set (n = 77). Urinary afamin levels were measured by ELISA in T2DM and DN patients enrolled in this cohort study (n = 203).

RESULTS

One-hundred-four proteins displayed significant alterations in excretion. Nine of these candidates were validated by MRM analysis. Among them, the levels of afamin, CD44 antigen, and lysosome-associated membrane glycoprotein 2, which have not previously been implicated in DN, were significantly associated with both the urinary albumin to creatinine ratio (ACR) and eGFR. We further measured afamin levels in urine collected from T2DM patients who did not yet have significant kidney disease (ACR < 300 mg/g or eGFR change rate ≤ 3.3%/year). The urinary afamin to creatinine ratio (Afa/Cre) was significantly higher in patients who progressed to a more severe DN stage or had early renal decline than in patients who did not.

CONCLUSIONS

Afa/Cre was significantly increased in T2DM patients who subsequently developed DN. Afa/Cre may be useful to predict patients with T2DM at high risk of nephropathy before the development of macroalbuminuria or reduced kidney function, although further validation studies in a larger population are needed.

Proteomic Analysis of Mouse Brain Subjected to Spaceflight

There is evidence that spaceflight poses acute and late risks to the central nervous system. To explore possible mechanisms, the proteomic changes following spaceflight in mouse brain were characterized. Space Shuttle Atlantis (STS-135) was launched from the Kennedy Space Center (KSC) on a 13-day mission. Within 3–5 h after landing, brain tissue was collected to evaluate protein expression profiles using quantitative proteomic analysis. Our results showed that there were 26 proteins that were significantly altered after spaceflight in the gray and/or white matter. While there was no overlap between the white and gray matter in terms of individual proteins, there was overlap in terms of function, synaptic plasticity, vesical activity, protein/organelle transport, and metabolism. Our data demonstrate that exposure to the spaceflight environment induces significant changes in protein expression related to neuronal structure and metabolic function. This might lead to a significant impact on brain structural and functional integrity that could affect the outcome of space missions.

Evaluation of Serum Glycoprotein Biomarker Candidates for Detection of Esophageal Adenocarcinoma and Surveillance of Barrett's Esophagus

Esophageal adenocarcinoma (EAC) is thought to develop from asymptomatic Barrett's esophagus (BE) with a low annual rate of conversion. Current endoscopy surveillance of BE patients is probably not cost-effective. Previously, we discovered serum glycoprotein biomarker candidates which could discriminate BE patients from EAC. Here, we aimed to validate candidate serum glycoprotein biomarkers in independent cohorts, and to develop a biomarker candidate panel for BE surveillance. Serum glycoprotein biomarker candidates were measured in 301 serum samples collected from Australia (4 states) and the United States (1 clinic) using previously established lectin magnetic bead array (LeMBA) coupled multiple reaction monitoring mass spectrometry (MRM-MS) tier 3 assay. The area under receiver operating characteristic curve (AUROC) was calculated as a measure of discrimination, and multivariate recursive partitioning was used to formulate a multi-marker panel for BE surveillance. Complement C9 (C9), gelsolin (GSN), serum paraoxonase/arylesterase 1 (PON1) and serum paraoxonase/lactonase 3 (PON3) were validated as diagnostic glycoprotein biomarkers in lectin pull-down samples for EAC across both cohorts. A panel of 10 serum glycoprotein biomarker candidates discriminated BE patients not requiring intervention (BE± low grade dysplasia) from those requiring intervention (BE with high grade dysplasia (BE-HGD) or EAC) with an AUROC value of 0.93. Tissue expression of C9 was found to be induced in BE, dysplastic BE and EAC. In longitudinal samples from subjects that have progressed toward EAC, levels of serum C9 were significantly (p < 0.05) increased with disease progression in EPHA (erythroagglutinin from Phaseolus vulgaris) and NPL (Narcissus pseudonarcissus lectin) pull-down samples. The results confirm alteration of complement pathway glycoproteins during BE-EAC pathogenesis. Further prospective clinical validation of the confirmed biomarker candidates in a large cohort is warranted, prior to development of a first-line BE surveillance blood test.

Proceedings of the EuBIC developer's meeting 2018

The inaugural European Bioinformatics Community (EuBIC) developer's meeting was held from January 9th to January 12th 2018 in Ghent, Belgium. While the meeting kicked off with an interactive keynote session featuring four internationally renowned experts in the field of computational proteomics, its primary focus were the hands-on hackathon sessions which featured six community-proposed projects revolving around three major topics: Here, we present an overview of the scientific program of the EuBIC developer's meeting and provide a starting point for follow-up on the covered projects.

Ultrafast Peptide Label-Free Quantification with FlashLFQ

The rapid and accurate quantification of peptides is a critical element of modern proteomics that has become increasingly challenging as proteomic data sets grow in size and complexity. We present here FlashLFQ, a computer program for high-speed label-free quantification of peptides following a search of bottom-up mass spectrometry data. FlashLFQ is approximately an order of magnitude faster than established label-free quantification methods. The increased speed makes it practical to base quantification upon all of the charge states for a given peptide rather than solely upon the charge state that was selected for MS2 fragmentation. This increases the number of quantified peptides, improves replicate-to-replicate reproducibility, and increases quantitative accuracy. We integrated FlashLFQ into the graphical user interface of the MetaMorpheus search software, allowing it to work together with the global post-translational modification discovery (G-PTM-D) engine to accurately quantify modified peptides. FlashLFQ is also available as a NuGet package, facilitating its integration into other software, and as a standalone command line software program for the quantification of search results from other programs (e.g., MaxQuant).