Optimization of Statistical Methods Impact on Quantitative Proteomics Data

Optimizing differential expression analysis for proteomics data via high-performing rules and ensemble inference

Identification of differentially expressed proteins in a proteomics workflow typically encompasses five key steps: raw data quantification, expression matrix construction, matrix normalization, missing value imputation (MVI), and differential expression analysis. The plethora of options in each step makes it challenging to identify optimal workflows that maximize the identification of differentially expressed proteins. To identify optimal workflows and their common properties, we conduct an extensive study involving 34,576 combinatoric experiments on 24 gold standard spike-in datasets. Applying frequent pattern mining techniques to top-ranked workflows, we uncover high-performing rules that demonstrate optimality has conserved properties. Via machine learning, we confirm optimal workflows are indeed predictable, with average cross-validation F1 scores and Matthew's correlation coefficients surpassing 0.84. We introduce an ensemble inference to integrate results from individual top-performing workflows for expanding differential proteome coverage and resolve inconsistencies. Ensemble inference provides gains in pAUC (up to 4.61%) and G-mean (up to 11.14%) and facilitates effective aggregation of information across varied quantification approaches such as topN, directLFQ, MaxLFQ intensities, and spectral counts. However, further development and evaluation are needed to establish acceptable frameworks for conducting ensemble inference on multiple proteomics workflows.

Exploring novel protein biomarkers for early-stage diagnosis and prognosis of T-acute lymphoblastic leukemia (T-ALL)

BACKGROUND

Efficient classification of T-acute lymphoblastic leukemia (T-ALL) involves considering various factors, such as age, white blood cell count, and chromosomal alterations. However, studying protein markers are crucial to improving T-ALL patients' diagnosis and treatment. A study analyzing the expression of proteomes was conducted to identify promising early-stage biomarkers for T-ALL patients METHODS: Label-free liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to analyze the blood proteins of both patients and healthy individuals to identify new biomarkers for T-ALL. The findings were validated by RT-PCR, ELISA and computational analysis RESULTS: The study identified 1467 proteins in the blood, of which nine were upregulated and 35 were downregulated by more than 2-fold. T-ALL patients showed a significant increase in specific disease-related proteins, such as eleven-nineteen lysine-rich leukemia protein, triggering receptor expressed on myeloid cells 1, cisplatin resistance-associated-overexpressed protein, X-ray radiation resistance-associated protein 1, tumor necrosis factor receptor superfamily member 10D, protein S100-A8, and copine-4, by more than 3-fold CONCLUSION: The findings of this study provide a valuable protein map of leukemic cells and identify potential biomarkers for leukemic aggressiveness. However, further studies using larger T-ALL patient samples must confirm these preliminary results.

CloudProteoAnalyzer: scalable processing of big data from proteomics using cloud computing

Summary

Shotgun proteomics is widely used in many system biology studies to determine the global protein expression profiles of tissues, cultures, and microbiomes. Many non-distributed computer algorithms have been developed for users to process proteomics data on their local computers. However, the amount of data acquired in a typical proteomics study has grown rapidly in recent years, owing to the increasing throughput of mass spectrometry and the expanding scale of study designs. This presents a big data challenge for researchers to process proteomics data in a timely manner. To overcome this challenge, we developed a cloud-based parallel computing application to offer end-to-end proteomics data analysis software as a service (SaaS). A web interface was provided to users to upload mass spectrometry-based proteomics data, configure parameters, submit jobs, and monitor job status. The data processing was distributed across multiple nodes in a supercomputer to achieve scalability for large datasets. Our study demonstrated SaaS for proteomics as a viable solution for the community to scale up the data processing using cloud computing.

Availability and implementation

This application is available online at https://sipros.oscer.ou.edu/ or https://sipros.unt.edu for free use. The source code is available at https://github.com/Biocomputing-Research-Group/CloudProteoAnalyzer under the GPL version 3.0 license.

Reporting and reproducibility: Proteomics of fish models in environmental toxicology and ecotoxicology

Environmental toxicology and ecotoxicology research efforts are employing proteomics with fish models as New Approach Methodologies, along with in silico, in vitro and other omics techniques to elucidate hazards of toxicants and toxins. We performed a critical review of toxicology studies with fish models using proteomics and reported fundamental parameters across experimental design, sample preparation, mass spectrometry, and bioinformatics of fish, which represent alternative vertebrate models in environmental toxicology, and routinely studied animals in ecotoxicology. We observed inconsistencies in reporting and methodologies among experimental designs, sample preparations, data acquisitions and bioinformatics, which can affect reproducibility of experimental results. We identified a distinct need to develop reporting guidelines for proteomics use in environmental toxicology and ecotoxicology, increased QA/QC throughout studies, and method optimization with an emphasis on reducing inconsistencies among studies. Several recommendations are offered as logical steps to advance development and application of this emerging research area to understand chemical hazards to public health and the environment.

Normalization methods in mass spectrometry-based analytical proteomics: A case study based on renal cell carcinoma datasets

Normalization is a crucial step in proteomics data analysis as it enables data adjustment and enhances comparability between datasets by minimizing multiple sources of variability, such as sampling, sample handling, storage, treatment, and mass spectrometry measurements. In this study, we investigated different normalization methods, including Z-score normalization, median divide normalization, and quantile normalization, to evaluate their performance using a case study based on renal cell carcinoma datasets. Our results demonstrate that when comparing datasets by pairs, both the Z-score and quantile normalization methods consistently provide better results in terms of the number of proteins identified and quantified as well as in identifying statistically significant up or down-regulated proteins. However, when three or more datasets are compared at the same time the differences are found to be negligible.

Benefits of FAIMS to Improve the Proteome Coverage of Deteriorated and/or Cross-Linked TMT 10-Plex FFPE Tissue and Plasma-Derived Exosomes Samples

The proteome characterization of complex, deteriorated, or cross-linked protein mixtures as paired clinical FFPE or exosome samples isolated from low plasma volumes (250 µL) might be a challenge. In this work, we aimed at investigating the benefits of FAIMS technology coupled to the Orbitrap Exploris 480 mass spectrometer for the TMT quantitative proteomics analyses of these complex samples in comparison to the analysis of protein extracts from cells, frozen tissue, and exosomes isolated from large volume plasma samples (3 mL). TMT experiments were performed using a two-hour gradient LC-MS/MS with or without FAIMS and two compensation voltages (CV = -45 and CV = -60). In the TMT experiments of cells, frozen tissue, or exosomes isolated from large plasma volumes (3 mL) with FAIMS, a limited increase in the number of identified and quantified proteins accompanied by a decrease in the number of peptides identified and quantified was observed. However, we demonstrated here a noticeable improvement (>100%) in the number of peptide and protein identifications and quantifications for the plasma exosomes isolated from low plasma volumes (250 µL) and FFPE tissue samples in TMT experiments with FAIMS in comparison to the LC-MS/MS analysis without FAIMS. Our results highlight the potential of mass spectrometry analyses with FAIMS to increase the depth into the proteome of complex samples derived from deteriorated, cross-linked samples and/or those where the material was scarce, such as FFPE and plasma-derived exosomes from low plasma volumes (250 µL), which might aid in the characterization of their proteome and proteoforms and in the identification of dysregulated proteins that could be used as biomarkers.

Advances in proteomics: characterization of the innate immune system after birth and during inflammation

Proteomics is the characterization of the protein composition, the proteome, of a biological sample. It involves the large-scale identification and quantification of proteins, peptides, and post-translational modifications. This review focuses on recent developments in mass spectrometry-based proteomics and provides an overview of available methods for sample preparation to study the innate immune system. Recent advancements in the proteomics workflows, including sample preparation, have significantly improved the sensitivity and proteome coverage of biological samples including the technically difficult blood plasma. Proteomics is often applied in immunology and has been used to characterize the levels of innate immune system components after perturbations such as birth or during chronic inflammatory diseases like rheumatoid arthritis (RA) and inflammatory bowel disease (IBD). In cancers, the tumor microenvironment may generate chronic inflammation and release cytokines to the circulation. In these situations, the innate immune system undergoes profound and long-lasting changes, the large-scale characterization of which may increase our biological understanding and help identify components with translational potential for guiding diagnosis and treatment decisions. With the ongoing technical development, proteomics will likely continue to provide increasing insights into complex biological processes and their implications for health and disease. Integrating proteomics with other omics data and utilizing multi-omics approaches have been demonstrated to give additional valuable insights into biological systems.

Proteomics Identifies Circulating TIMP-1 as a Prognostic Biomarker for Diffuse Large B-Cell Lymphoma

Diffuse large B-cell lymphoma (DLBCL) is a heterogeneous disease, although disease stratification using in-depth plasma proteomics has not been performed to date. By measuring more than 1000 proteins in the plasma of 147 DLBCL patients using data-independent acquisition mass spectrometry and antibody array, DLBCL patients were classified into four proteomic subtypes (PS-I-IV). Patients with the PS-IV subtype and worst prognosis had increased levels of proteins involved in inflammation, including a high expression of metalloproteinase inhibitor-1 (TIMP-1) that was associated with poor survival across two validation cohorts (n = 180). Notably, the combination of TIMP-1 with the international prognostic index (IPI) identified 64.00% to 88.24% of relapsed and 65.00% to 80.49% of deceased patients in the discovery and two validation cohorts, which represents a 24.00% to 41.67% and 20.00% to 31.70% improvement compared to the IPI score alone, respectively. Taken together, we demonstrate that DLBCL heterogeneity is reflected in the plasma proteome and that TIMP-1, together with the IPI, could improve the prognostic stratification of patients.

Proteomics in Childhood Acute Lymphoblastic Leukemia: Challenges and Opportunities

Acute lymphoblastic leukemia (ALL) is the most common cancer in children and one of the success stories in cancer therapeutics. Risk-directed therapy based on clinical, biologic and genetic features has played a significant role in this accomplishment. Despite the observed improvement in survival rates, leukemia remains one of the leading causes of cancer-related deaths. Implementation of next-generation genomic and transcriptomic sequencing tools has illustrated the genomic landscape of ALL. However, the underlying dynamic changes at protein level still remain a challenge. Proteomics is a cutting-edge technology aimed at deciphering the mechanisms, pathways, and the degree to which the proteome impacts leukemia subtypes. Advances in mass spectrometry enable high-throughput collection of global proteomic profiles, representing an opportunity to unveil new biological markers and druggable targets. The purpose of this narrative review article is to provide a comprehensive overview of studies that have utilized applications of proteomics in an attempt to gain insight into the pathogenesis and identification of biomarkers in childhood ALL.

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.

prolfqua: A Comprehensive R-Package for Proteomics Differential Expression Analysis

Mass spectrometry is widely used for quantitative proteomics studies, relative protein quantification, and differential expression analysis of proteins. There is a large variety of quantification software and analysis tools. Nevertheless, there is a need for a modular, easy-to-use application programming interface in R that transparently supports a variety of well principled statistical procedures to make applying them to proteomics data, comparing and understanding their differences easy. The prolfqua package integrates essential steps of the mass spectrometry-based differential expression analysis workflow: quality control, data normalization, protein aggregation, statistical modeling, hypothesis testing, and sample size estimation. The package makes integrating new data formats easy. It can be used to model simple experimental designs with a single explanatory variable and complex experiments with multiple factors and hypothesis testing. The implemented methods allow sensitive and specific differential expression analysis. Furthermore, the package implements benchmark functionality that can help to compare data acquisition, data preprocessing, or data modeling methods using a gold standard data set. The application programmer interface of prolfqua strives to be clear, predictable, discoverable, and consistent to make proteomics data analysis application development easy and exciting. Finally, the prolfqua R-package is available on GitHub https://github.com/fgcz/prolfqua, distributed under the MIT license. It runs on all platforms supported by the R free software environment for statistical computing and graphics.

Label-free proteome quantification and evaluation

The label-free quantification (LFQ) has emerged as an exceptional technique in proteomics owing to its broad proteome coverage, great dynamic ranges and enhanced analytical reproducibility. Due to the extreme difficulty lying in an in-depth quantification, the LFQ chains incorporating a variety of transformation, pretreatment and imputation methods are required and constructed. However, it remains challenging to determine the well-performing chain, owing to its strong dependence on the studied data and the diverse possibility of integrated chains. In this study, an R package EVALFQ was therefore constructed to enable a performance evaluation on >3000 LFQ chains. This package is unique in (a) automatically evaluating the performance using multiple criteria, (b) exploring the quantification accuracy based on spiking proteins and (c) discovering the well-performing chains by comprehensive assessment. All in all, because of its superiority in assessing from multiple perspectives and scanning among over 3000 chains, this package is expected to attract broad interests from the fields of proteomic quantification. The package is available at https://github.com/idrblab/EVALFQ.

Benchmarking tools for detecting longitudinal differential expression in proteomics data allows establishing a robust reproducibility optimization regression approach

Quantitative proteomics has matured into an established tool and longitudinal proteomics experiments have begun to emerge. However, no effective, simple-to-use differential expression method for longitudinal proteomics data has been released. Typically, such data is noisy, contains missing values, and has only few time points and biological replicates. To address this need, we provide a comprehensive evaluation of several existing differential expression methods for high-throughput longitudinal omics data and introduce a Robust longitudinal Differential Expression (RolDE) approach. The methods are evaluated using over 3000 semi-simulated spike-in proteomics datasets and three large experimental datasets. In the comparisons, RolDE performs overall best; it is most tolerant to missing values, displays good reproducibility and is the top method in ranking the results in a biologically meaningful way. Furthermore, RolDE is suitable for different types of data with typically unknown patterns in longitudinal expression and can be applied by non-experienced users.

Comparative tissue proteomics reveals unique action mechanisms of vaccine adjuvants

Radiofrequency adjuvant (RFA) was recently developed to boost influenza vaccination without the safety concerns of chemical adjuvants due to their physical nature. Yet, the action mechanisms of RFA remain largely unknown. Omics techniques offer new opportunities to identify molecular mechanisms of RFA. This study utilized comparative tissue proteomics to explore molecular mechanisms of the physical RFA. Comparison of RFA and chemical adjuvant (Alum, AddaVax, MPL, MPL/Alum)-induced tissue proteome changes identified 14 exclusively induced proteins by RFA, among which heat shock protein (HSP) 70 was selected for further analysis due to its known immune-modulating functions. RFA showed much weakened ability to boost ovalbumin and pandemic influenza vaccination in HSP70 knockout than wild-type mice, hinting crucial roles of HSP70 in RFA effects. This study supports comparative tissue proteomics to be an effective tool to study molecular mechanisms of vaccine adjuvants.

Normics: Proteomic Normalization by Variance and Data-Inherent Correlation Structure

Several algorithms for the normalization of proteomic data are currently available, each based on a priori assumptions. Among these is the extent to which differential expression (DE) can be present in the dataset. This factor is usually unknown in explorative biomarker screens. Simultaneously, the increasing depth of proteomic analyses often requires the selection of subsets with a high probability of being DE to obtain meaningful results in downstream bioinformatical analyses. Based on the relationship of technical variation and (true) biological DE of an unknown share of proteins, we propose the “Normics” algorithm: Proteins are ranked based on their expression level–corrected variance and the mean correlation with all other proteins. The latter serves as a novel indicator of the non-DE likelihood of a protein in a given dataset. Subsequent normalization is based on a subset of non-DE proteins only. No a priori information such as batch, clinical, or replicate group is necessary. Simulation data demonstrated robust and superior performance across a wide range of stochastically chosen parameters. Five publicly available spike-in and biologically variant datasets were reliably and quantitively accurately normalized by Normics with improved performance compared to standard variance stabilization as well as median, quantile, and LOESS normalizations. In complex biological datasets Normics correctly determined proteins as being DE that had been cross-validated by an independent transcriptome analysis of the same samples. In both complex datasets Normics identified the most DE proteins. We demonstrate that combining variance analysis and data-inherent correlation structure to identify non-DE proteins improves data normalization. Standard normalization algorithms can be consolidated against high shares of (one-sided) biological regulation. The statistical power of downstream analyses can be increased by focusing on Normics-selected subsets of high DE likelihood.

•

Normics is a tool for the normalization of proteomic data based on existing algorithms.

•

Specifically addresses data with high shares of differential expression.

•

Combines variance and data-inherent correlation structure.

•

Provides a ranking of differential expression likelihood.

•

Enables normalization based on the most stable proteins.

Benchmarking of analysis strategies for data-independent acquisition proteomics using a large-scale dataset comprising inter-patient heterogeneity

Numerous software tools exist for data-independent acquisition (DIA) analysis of clinical samples, necessitating their comprehensive benchmarking. We present a benchmark dataset comprising real-world inter-patient heterogeneity, which we use for in-depth benchmarking of DIA data analysis workflows for clinical settings. Combining spectral libraries, DIA software, sparsity reduction, normalization, and statistical tests results in 1428 distinct data analysis workflows, which we evaluate based on their ability to correctly identify differentially abundant proteins. From our dataset, we derive bootstrap datasets of varying sample sizes and use the whole range of bootstrap datasets to robustly evaluate each workflow. We find that all DIA software suites benefit from using a gas-phase fractionated spectral library, irrespective of the library refinement used. Gas-phase fractionation-based libraries perform best against two out of three reference protein lists. Among all investigated statistical tests non-parametric permutation-based statistical tests consistently perform best.

Data independent acquisition (DIA) has been gaining momentum in clinical proteomics. Here, the authors create a benchmark dataset comprising inter-patient heterogeneity to compare popular DIA data analysis workflows for identifying differentially abundant proteins.

Benchmarking differential expression, imputation and quantification methods for proteomics data

Data analysis is a critical part of quantitative proteomics studies in interpreting biological questions. Numerous computational tools for protein quantification, imputation and differential expression (DE) analysis were generated in the past decade and the search for optimal tools is still going on. Moreover, due to the rapid development of RNA sequencing (RNA-seq) technology, a vast number of DE analysis methods were created for that purpose. The applicability of these newly developed RNA-seq-oriented tools to proteomics data remains in doubt. In order to benchmark these analysis methods, a proteomics dataset consisting of proteins derived from humans, yeast and drosophila, in defined ratios, was generated in this study. Based on this dataset, DE analysis tools, including microarray- and RNA-seq-based ones, imputation algorithms and protein quantification methods were compared and benchmarked. Furthermore, applying these approaches to two public datasets showed that RNA-seq-based DE tools achieved higher accuracy (ACC) in identifying DEPs. This study provides useful guidelines for analyzing quantitative proteomics datasets. All the methods used in this study were integrated into the Perseus software, version 2.0.3.0, which is available at https://www.maxquant.org/perseus.

StatsPro: Systematic integration and evaluation of statistical approaches for detecting differential expression in label-free quantitative proteomics

Quantitative label-free mass spectrometry (MS) is an increasingly powerful technology for profiling thousands of proteins from complex biological samples. One of the primary goals of analyses performed on such proteomics data is to detect differentially expressed proteins (DEPs) under different experimental conditions. Many statistical methods have been developed and assessed for DEP detection in various proteomics studies. However, it remains a challenge for many proteomics scientists to choose an appropriate statistical procedure. Therefore, in this study, we organized 12 common testing algorithms and 6 P-value combination methods and further provided Cohen's d effect size for every protein and three evaluation criteria to help proteomics scientists investigate their influence on DEP detection in a systematic manner. To promote the widespread use of these methods, we developed a user-friendly web tool, StatsPro, and presented two case studies involving label-free quantitative proteomics data obtained using data-dependent acquisition and data-independent acquisition to illustrate its practicability. This tool is freely available in our GitHub repository (https://github.com/YanglabWCH/StatsPro/). SIGNIFICANCE: One of the primary goals of analyses performed on liquid chromatography-mass spectrometry (LC-MS) based proteomics data is to detect differentially expressed proteins (DEPs) under different experimental conditions. Despite of many research efforts have been proposed to detect DEPs, to date, there is a scarcity of efficient, systematic, and easy-to-handle tools that are tailored for proteomics scientists to choose an appropriate statistical procedure. Herein, we present a new tool, StatsPro, to enable implementation and evaluation of different statistical methods for proteomics scientists. This tool has two significant advances compared to existing software: a) It integrates up to 18 common statistical approaches (12 statistical tests and 6 P-value combination strategies) and performs Cohen's d effect size systematically for users, moreover, it provides a web-based interface and can be quite conveniently operated by users, even those with less profound computational background. b) It supports three performance evaluation criteria (e.g. number of DEPs, correlation coefficient between P-values and effect sizes, Area under the ROC curve) for users to review the final statistical results, which may guide the method selection for DEPs detection.

Differential ATAC-seq and ChIP-seq peak detection using ROTS

Changes in cellular chromatin states fine-tune transcriptional output and ultimately lead to phenotypic changes. Here we propose a novel application of our reproducibility-optimized test statistics (ROTS) to detect differential chromatin states (ATAC-seq) or differential chromatin modification states (ChIP-seq) between conditions. We compare the performance of ROTS to existing and widely used methods for ATAC-seq and ChIP-seq data using both synthetic and real datasets. Our results show that ROTS outperformed other commonly used methods when analyzing ATAC-seq data. ROTS also displayed the most accurate detection of small differences when modeling with synthetic data. We observed that two-step methods that require the use of a separate peak caller often more accurately called enrichment borders, whereas one-step methods without a separate peak calling step were more versatile in calling sub-peaks. The top ranked differential regions detected by the methods had marked correlation with transcriptional differences of the closest genes. Overall, our study provides evidence that ROTS is a useful addition to the available differential peak detection methods to study chromatin and performs especially well when applied to study differential chromatin states in ATAC-seq data.

ANPELA: analysis and performance assessment of the label-free quantification workflow for metaproteomic studies

Label-free quantification (LFQ) with a specific and sequentially integrated workflow of acquisition technique, quantification tool and processing method has emerged as the popular technique employed in metaproteomic research to provide a comprehensive landscape of the adaptive response of microbes to external stimuli and their interactions with other organisms or host cells. The performance of a specific LFQ workflow is highly dependent on the studied data. Hence, it is essential to discover the most appropriate one for a specific data set. However, it is challenging to perform such discovery due to the large number of possible workflows and the multifaceted nature of the evaluation criteria. Herein, a web server ANPELA (https://idrblab.org/anpela/) was developed and validated as the first tool enabling performance assessment of whole LFQ workflow (collective assessment by five well-established criteria with distinct underlying theories), and it enabled the identification of the optimal LFQ workflow(s) by a comprehensive performance ranking. ANPELA not only automatically detects the diverse formats of data generated by all quantification tools but also provides the most complete set of processing methods among the available web servers and stand-alone tools. Systematic validation using metaproteomic benchmarks revealed ANPELA's capabilities in 1 discovering well-performing workflow(s), (2) enabling assessment from multiple perspectives and (3) validating LFQ accuracy using spiked proteins. ANPELA has a unique ability to evaluate the performance of whole LFQ workflow and enables the discovery of the optimal LFQs by the comprehensive performance ranking of all 560 workflows. Therefore, it has great potential for applications in metaproteomic and other studies requiring LFQ techniques, as many features are shared among proteomic studies.

DDASSQ: an open-source, multiple peptide sequencing strategy for label free quantification based on an OpenMS pipeline in the KNIME analytics platform

In this study we investigated the performance of a computational pipeline for protein identification and label free quantification (LFQ) of LC–MS/MS data sets from experimental animal tissue samples, as well as the impact of its specific peptide search combinatorial approach. The full pipeline workflow was composed of peptide search engine adapters based on different identification algorithms, in the frame of the open‐source OpenMS software running within the KNIME analytics platform. Two different in silico tryptic digestion, database‐search assisted approaches (X!Tandem and MS‐GF+), de novo peptide sequencing based on Novor and consensus library search (SpectraST), were tested for the processing of LC‐MS/MS raw data files obtained from proteomic LC‐MS experiments done on proteolytic extracts from mouse ex vivo liver samples. The results from proteomic LFQ were compared to those based on the application of the two software tools MaxQuant and Proteome Discoverer for protein inference and label‐free data analysis in shotgun proteomics. Data are available via ProteomeXchange with identifier PXD025097.

DIMA: Data-Driven Selection of an Imputation Algorithm

Imputation is a prominent strategy when dealing with missing values (MVs) in proteomics data analysis pipelines. However, it is difficult to assess the performance of different imputation methods and varies strongly depending on data characteristics. To overcome this issue, we present the concept of a data-driven selection of an imputation algorithm (DIMA). The performance and broad applicability of DIMA are demonstrated on 142 quantitative proteomics data sets from the PRoteomics IDEntifications (PRIDE) database and on simulated data consisting of 5-50% MVs with different proportions of missing not at random and missing completely at random values. DIMA reliably suggests a high-performing imputation algorithm, which is always among the three best algorithms and results in a root mean square error difference (ΔRMSE) ≤ 10% in 80% of the cases. DIMA implementation is available in MATLAB at github.com/kreutz-lab/OmicsData and in R at github.com/kreutz-lab/DIMAR.

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.

Differential production and secretion of potentially toxigenic extracellular proteins from hypervirulent Aeromonas hydrophila under biofilm and planktonic culture

Background

Hypervirulent Aeromonas hydrophila (vAh) is an emerging pathogen in freshwater aquaculture that results in the loss of over 3 million pounds of marketable channel catfish, Ictalurus punctatus, and channel catfish hybrids (I. punctatus, ♀ x blue catfish, I. furcatus, ♂) each year from freshwater catfish production systems in Alabama, U.S.A. vAh isolates are clonal in nature and are genetically unique from, and significantly more virulent than, traditional A. hydrophila isolates from fish. Even with the increased virulence, natural infections cannot be reproduced in aquaria challenges making it difficult to determine modes of infection and the pathophysiology behind the devastating mortalities that are commonly observed. Despite the intimate connection between environmental adaptation and plastic response, the role of environmental adaption on vAh pathogenicity and virulence has not been previously explored. In this study, secreted proteins of vAh cultured as free-living planktonic cells and within a biofilm were compared to elucidate the role of biofilm growth on virulence.

Results

Functional proteolytic assays found significantly increased degradative activity in biofilm secretomes; in contrast, planktonic secretomes had significantly increased hemolytic activity, suggesting higher toxigenic potential. Intramuscular injection challenges in a channel catfish model showed that in vitro degradative activity translated into in vivo tissue destruction. Identification of secreted proteins by HPLC-MS/MS revealed the presence of many putative virulence proteins under both growth conditions. Biofilm grown vAh produced higher levels of proteolytic enzymes and adhesins, whereas planktonically grown cells secreted higher levels of toxins, porins, and fimbrial proteins.

Conclusions

This study is the first comparison of the secreted proteomes of vAh when grown in two distinct ecological niches. These data on the adaptive physiological response of vAh based on growth condition increase our understanding of how environmental niche partitioning could affect vAh pathogenicity and virulence. Increased secretion of colonization factors and degradative enzymes during biofilm growth and residency may increase bacterial attachment and host invasiveness, while increased secretion of hemolysins, porins, and other potential toxins under planktonic growth (or after host invasion) could result in increased host mortality. The results of this research underscore the need to use culture methods that more closely mimic natural ecological habitat growth to improve our understanding of vAh pathogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-020-02065-2.

Proteomic analysis of radio-resistant breast cancer xenografts: Increased TGF-β signaling and metabolism

Our previous studies have shown that MCF-7 breast cancer cell line exposed to 6 Gy and allowed to recover for 7 days (D7-6G) developed radio-resistance. In this study, we have tested the ability of these cells to form tumors in severe combined immunodeficiency (SCID) mice and characterized these tumors by proteomic analyses. Untreated (MCF-C) and D7-6G cells (MCF-R) were injected s.c. in SCID mice and tumor growth monitored. On Day 18, the mice were killed and tumor tissues were fixed in formalin or RNA later. Expression of genes was assessed by reverse transcription-polymerase chain reaction and proteins by enzyme-linked immunosorbent assay/antibody labeling and flow cytometry. Label free proteomic analyses was carried out by liquid chromatography-mass spectrometry. Metabolic analysis was carried out using Seahorse analyzer. MCF-R cells had a shorter latency and formed larger tumors. These tumors were characterized by an increased expression of transforming growth factor β (TGF-β) isoforms; its downstream genes pSMAD3, Snail-1, Zeb-1, HMGA2; hybrid epithelial/mesenchymal phenotype; migration, enrichment of cancer stem cells and radioresistance following challenge dose of radiation. Proteomic analysis of MCF-7R tumors resulted in identification of a total of 649 differentially expressed proteins and pathway analyses using protein annotation through evolutionary relationship indicated enrichment of genes involved in metabolism. Data are available via ProteomeXchange with identifier PXD022506. Seahorse analyzer confirmed increased metabolism in these cells with increased oxidative phosphorylation as well as glycolysis. Increased uptake of 2-NBDG further confirmed increased glycolysis. In summary, we demonstrate that radioresistant breast cancer cells had an enrichment of TGF-β signaling and increased metabolism.

Utility of Proteomics in Emerging and Re-Emerging Infectious Diseases Caused by RNA Viruses

Emerging and re-emerging infectious diseases due to RNA viruses cause major negative consequences for the quality of life, public health, and overall economic development. Most of the RNA viruses causing illnesses in humans are of zoonotic origin. Zoonotic viruses can directly be transferred from animals to humans through adaptation, followed by human-to-human transmission, such as in human immunodeficiency virus (HIV), severe acute respiratory syndrome coronavirus (SARS-CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), and, more recently, SARS coronavirus 2 (SARS-CoV-2), or they can be transferred through insects or vectors, as in the case of Crimean-Congo hemorrhagic fever virus (CCHFV), Zika virus (ZIKV), and dengue virus (DENV). At the present, there are no vaccines or antiviral compounds against most of these viruses. Because proteins possess a vast array of functions in all known biological systems, proteomics-based strategies can provide important insights into the investigation of disease pathogenesis and the identification of promising antiviral drug targets during an epidemic or pandemic. Mass spectrometry technology has provided the capacity required for the precise identification and the sensitive and high-throughput analysis of proteins on a large scale and has contributed greatly to unravelling key protein-protein interactions, discovering signaling networks, and understanding disease mechanisms. In this Review, we present an account of quantitative proteomics and its application in some prominent recent examples of emerging and re-emerging RNA virus diseases like HIV-1, CCHFV, ZIKV, and DENV, with more detail with respect to coronaviruses (MERS-CoV and SARS-CoV) as well as the recent SARS-CoV-2 pandemic.

Eatomics: Shiny Exploration of Quantitative Proteomics Data

Quantitative proteomics data are becoming increasingly more available, and as a consequence are being analyzed and interpreted by a larger group of users. However, many of these users have less programming experience. Furthermore, experimental designs and setups are getting more complicated, especially when tissue biopsies are analyzed. Luckily, the proteomics community has already established some best practices on how to conduct quality control, differential abundance analysis and enrichment analysis. However, an easy-to-use application that wraps together all steps for the exploration and flexible analysis of quantitative proteomics data is not yet available. For Eatomics, we utilize the R Shiny framework to implement carefully chosen parts of established analysis workflows to (i) make them accessible in a user-friendly way, (ii) add a multitude of interactive exploration possibilities, and (iii) develop a unique experimental design setup module, which interactively translates a given research hypothesis into a differential abundance and enrichment analysis formula. In this, we aim to fulfill the needs of a growing group of inexperienced quantitative proteomics data analysts. Eatomics may be tested with demo data directly online via https://we.analyzegenomes.com/now/eatomics/ or with the user's own data by installation from the Github repository at https://github.com/Millchmaedchen/Eatomics.

MetaFS: Performance assessment of biomarker discovery in metaproteomics

Metaproteomics suffers from the issues of dimensionality and sparsity. Data reduction methods can maximally identify the relevant subset of significant differential features and reduce data redundancy. Feature selection (FS) methods were applied to obtain the significant differential subset. So far, a variety of feature selection methods have been developed for metaproteomic study. However, due to FS's performance depended heavily on the data characteristics of a given research, the well-suitable feature selection method must be carefully selected to obtain the reproducible differential proteins. Moreover, it is critical to evaluate the performance of each FS method according to comprehensive criteria, because the single criterion is not sufficient to reflect the overall performance of the FS method. Therefore, we developed an online tool named MetaFS, which provided 13 types of FS methods and conducted the comprehensive evaluation on the complex FS methods using four widely accepted and independent criteria. Furthermore, the function and reliability of MetaFS were systematically tested and validated via two case studies. In sum, MetaFS could be a distinguished tool for discovering the overall well-performed FS method for selecting the potential biomarkers in microbiome studies. The online tool is freely available at https://idrblab.org/metafs/.

DEqMS: A Method for Accurate Variance Estimation in Differential Protein Expression Analysis

Quantitative proteomics by mass spectrometry is widely used in biomarker research and basic biology research for investigation of phenotype level cellular events. Despite the wide application, the methodology for statistical analysis of differentially expressed proteins has not been unified. Various methods such as t test, linear model and mixed effect models are used to define changes in proteomics experiments. However, none of these methods consider the specific structure of MS-data. Choices between methods, often originally developed for other types of data, are based on compromises between features such as statistical power, general applicability and user friendliness. Furthermore, whether to include proteins identified with one peptide in statistical analysis of differential protein expression varies between studies. Here we present DEqMS, a robust statistical method developed specifically for differential protein expression analysis in mass spectrometry data. In all data sets investigated there is a clear dependence of variance on the number of PSMs or peptides used for protein quantification. DEqMS takes this feature into account when assessing differential protein expression. This allows for a more accurate data-dependent estimation of protein variance and inclusion of single peptide identifications without increasing false discoveries. The method was tested in several data sets including E. coli proteome spike-in data, using both label-free and TMT-labeled quantification. Compared with previous statistical methods used in quantitative proteomics, DEqMS showed consistently better accuracy in detecting altered protein levels compared with other statistical methods in both label-free and labeled quantitative proteomics data. DEqMS is available as an R package in Bioconductor.

SWATH-MS analysis of cerebrospinal fluid to generate a robust battery of biomarkers for Alzheimer's disease

Cerebrospinal fluid (CSF) Aβ42 and tau protein levels are established diagnostic biomarkers of Alzheimer's disease (AD). However, their inadequacy to represent clinical efficacy in drug trials indicates the need for new biomarkers. Sequential window acquisition of all theoretical fragment ion spectra (SWATH)-based mass spectrometry (MS) is an advanced proteomic tool for large-scale, high-quality quantification. In this study, SWATH-MS showed that VGF, chromogranin-A, secretogranin-1, and opioid-binding protein/cell adhesion molecule were significantly decreased in 42 AD patients compared to 39 controls, whereas 14-3-3ζ was increased (FDR < 0.05). In addition, 16 other proteins showed substantial changes (FDR < 0.2). The expressions of the top 21 analytes were closely interconnected, but were poorly correlated with CSF Aβ42, tTau, and pTau181 levels. Logistic regression analysis and data mining were used to establish the best algorithm for AD, which created novel biomarker panels with high diagnostic value (AUC = 0.889 and 0.924) and a strong correlation with clinical severity (all p < 0.001). Targeted proteomics was used to validate their usefulness in a different cohort (n = 36) that included patients with other brain disorders (all p < 0.05). This study provides a list of proteins (and combinations thereof) that could serve as new AD biomarkers.

A comparative proteomic study of plasma in Colombian childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Owing to the incorporation of risk-adapted therapy and the arrival of new directed agents, the cure rate and survival of patients with ALL have improved dramatically, get near to 90%. In Latin American countries, the mortality rates of ALL are high, for example in Colombia, during the last decade, ALL has been the most prevalent cancer among children between 0–14 years of age. In the face of this public health problem and coupled with the fact that the knowledge of the proteome of the child population is little, our investigation proposes the study of the plasma proteome of Colombian children diagnosed with B-cell ALL (B-ALL) to determine potential disease markers that could reflect processes altered by the presence of the disease or in response to it. A proteomic study by LC-MS/MS and quantification by label-free methods were performed in search of proteins differentially expressed between healthy children and those diagnosed with B-ALL. We quantified a total of 472 proteins in depleted blood plasma, and 25 of these proteins were differentially expressed (fold change >2, Bonferroni-adjusted P-values <0.05). Plasma Aggrecan core protein, alpha-2-HS-glycoprotein, coagulation factor XIII A chain and gelsolin protein were examined by ELISA assay and compared to shotgun proteomics results. Our data provide new information on the plasma proteome of Colombian children. Additionally, these proteins may also have certain potential as illness markers or as therapeutic targets in subsequent investigations.

Simultaneous Improvement in the Precision, Accuracy, and Robustness of Label-free Proteome Quantification by Optimizing Data Manipulation Chains

The label-free proteome quantification (LFQ) is multistep workflow collectively defined by quantification tools and subsequent data manipulation methods that has been extensively applied in current biomedical, agricultural, and environmental studies. Despite recent advances, in-depth and high-quality quantification remains extremely challenging and requires the optimization of LFQs by comparatively evaluating their performance. However, the evaluation results using different criteria (precision, accuracy, and robustness) vary greatly, and the huge number of potential LFQs becomes one of the bottlenecks in comprehensively optimizing proteome quantification. In this study, a novel strategy, enabling the discovery of the LFQs of simultaneously enhanced performance from thousands of workflows (integrating 18 quantification tools with 3,128 manipulation chains), was therefore proposed. First, the feasibility of achieving simultaneous improvement in the precision, accuracy, and robustness of LFQ was systematically assessed by collectively optimizing its multistep manipulation chains. Second, based on a variety of benchmark datasets acquired by various quantification measurements of different modes of acquisition, this novel strategy successfully identified a number of manipulation chains that simultaneously improved the performance across multiple criteria. Finally, to further enhance proteome quantification and discover the LFQs of optimal performance, an online tool (https://idrblab.org/anpela/) enabling collective performance assessment (from multiple perspectives) of the entire LFQ workflow was developed. This study confirmed the feasibility of achieving simultaneous improvement in precision, accuracy, and robustness. The novel strategy proposed and validated in this study together with the online tool might provide useful guidance for the research field requiring the mass-spectrometry-based LFQ technique.

Density distribution of gene expression profiles and evaluation of using maximal information coefficient to identify differentially expressed genes

The hypothesis of data probability density distributions has many effects on the design of a new statistical method. Based on the analysis of a group of real gene expression profiles, this study reveal that the primary density distributions of the real profiles are normal/log-normal and t distributions, accounting for 80% and 19% respectively. According to these distributions, we generated a series of simulation data to make a more comprehensive assessment for a novel statistical method, maximal information coefficient (MIC). The results show that MIC is not only in the top tier in the overall performance of identifying differentially expressed genes, but also exhibits a better adaptability and an excellent noise immunity in comparison with the existing methods.

Antagonistic Interactions between Benzo[a]pyrene and Fullerene (C60) in Toxicological Response of Marine Mussels

This study aimed to assess the ecotoxicological effects of the interaction of fullerene (C60) and benzo[a]pyrene (B[a]P) on the marine mussel, Mytilus galloprovincialis. The uptake of nC60, B[a]P and mixtures of nC60 and B[a]P into tissues was confirmed by Gas Chromatography-Mass Spectrometry (GC-MS), Liquid Chromatography-High Resolution Mass Spectrometry (LC-HRMS) and Inductively Coupled Plasma Mass Spectrometer (ICP-MS). Biomarkers of DNA damage as well as proteomics analysis were applied to unravel the interactive effect of B[a]P and C60. Antagonistic responses were observed at the genotoxic and proteomic level. Differentially expressed proteins (DEPs) were only identified in the B[a]P single exposure and the B[a]P mixture exposure groups containing 1 mg/L of C60, the majority of which were downregulated (~52%). No DEPs were identified at any of the concentrations of nC60 (p < 0.05, 1% FDR). Using DEPs identified at a threshold of (p < 0.05; B[a]P and B[a]P mixture with nC60), gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that these proteins were enriched with a broad spectrum of biological processes and pathways, including those broadly associated with protein processing, cellular processes and environmental information processing. Among those significantly enriched pathways, the ribosome was consistently the top enriched term irrespective of treatment or concentration and plays an important role as the site of biological protein synthesis and translation. Our results demonstrate the complex multi-modal response to environmental stressors in M. galloprovincialis.

A comprehensive evaluation of popular proteomics software workflows for label-free proteome quantification and imputation

Label-free mass spectrometry (MS) has developed into an important tool applied in various fields of biological and life sciences. Several software exist to process the raw MS data into quantified protein abundances, including open source and commercial solutions. Each software includes a set of unique algorithms for different tasks of the MS data processing workflow. While many of these algorithms have been compared separately, a thorough and systematic evaluation of their overall performance is missing. Moreover, systematic information is lacking about the amount of missing values produced by the different proteomics software and the capabilities of different data imputation methods to account for them.In this study, we evaluated the performance of five popular quantitative label-free proteomics software workflows using four different spike-in data sets. Our extensive testing included the number of proteins quantified and the number of missing values produced by each workflow, the accuracy of detecting differential expression and logarithmic fold change and the effect of different imputation and filtering methods on the differential expression results. We found that the Progenesis software performed consistently well in the differential expression analysis and produced few missing values. The missing values produced by the other software decreased their performance, but this difference could be mitigated using proper data filtering or imputation methods. Among the imputation methods, we found that the local least squares (lls) regression imputation consistently increased the performance of the software in the differential expression analysis, and a combination of both data filtering and local least squares imputation increased performance the most in the tested data sets.

A critical assessment of the feature selection methods used for biomarker discovery in current metaproteomics studies

Microbial community (MC) has great impact on mediating complex disease indications, biogeochemical cycling and agricultural productivities, which makes metaproteomics powerful technique for quantifying diverse and dynamic composition of proteins or peptides. The key role of biostatistical strategies in MC study is reported to be underestimated, especially the appropriate application of feature selection method (FSM) is largely ignored. Although extensive efforts have been devoted to assessing the performance of FSMs, previous studies focused only on their classification accuracy without considering their ability to correctly and comprehensively identify the spiked proteins. In this study, the performances of 14 FSMs were comprehensively assessed based on two key criteria (both sample classification and spiked protein discovery) using a variety of metaproteomics benchmarks. First, the classification accuracies of those 14 FSMs were evaluated. Then, their abilities in identifying the proteins of different spiked concentrations were assessed. Finally, seven FSMs (FC, LMEB, OPLS-DA, PLS-DA, SAM, SVM-RFE and T-Test) were identified as performing consistently superior or good under both criteria with the PLS-DA performing consistently superior. In summary, this study served as comprehensive analysis on the performances of current FSMs and could provide a valuable guideline for researchers in metaproteomics.

Maximal information coefficient applied to differentially expressed genes identification: A feasibility study

BACKGROUND:

The main obstacle encountered in microarray technology is how to mine the valuable information under the profiles and study the genes function.

OBJECTIVE:

Maximal information coefficient (MIC) is a novel, non-parametric statistic that has been successfully applied to genome-wide association studies and differentially gene and miRNA expression analysis. However, the data used in these applications are not gold standard but real data.

METHODS:

Therefore, this study attempts to test the feasibility of MIC for differentially expressed gene identification with simulation data.

RESULTS:

Our experiments indicate that, MIC perfermance is better than Limma always, which is almost the same level of SAM, ROTS or DESeq2. However, the count of AUC < 0.5 of MIC is significantly smaller than the three methods, and MIC does not exhibit an abnormal phenomenon in which the AUC increases as the noise increases.

CONCLUSIONS:

Compared to the existing methods, our experiments show that MIC is not only in the first tier in identifying differentially expressed genes and noise immunity, but also shows better robustness and stronger data/environment adaptability.

Analysis of significant protein abundance from multiple reaction-monitoring data

Background

Discovering reliable protein biomarkers is one of the most important issues in biomedical research. The ELISA is a traditional technique for accurate quantitation of well-known proteins. Recently, the multiple reaction-monitoring (MRM) mass spectrometry has been proposed for quantifying newly discovered protein and has become a popular alternative to ELISA. For the MRM data analysis, linear mixed modeling (LMM) has been used to analyze MRM data. MSstats is one of the most widely used tools for MRM data analysis that is based on the LMMs. However, LMMs often provide various significance results, depending on model specification. Sometimes it would be difficult to specify a correct LMM method for the analysis of MRM data. Here, we propose a new logistic regression-based method for Significance Analysis of Multiple Reaction Monitoring (LR-SAM).

Results

Through simulation studies, we demonstrate that LMM methods may not preserve type I error, thus yielding high false- positive errors, depending on how random effects are specified. Our simulation study also shows that the LR-SAM approach performs similarly well as LMM approaches, in most cases. However, LR-SAM performs better than the LMMs, particularly when the effects sizes of peptides from the same protein are heterogeneous. Our proposed method was applied to MRM data for identification of proteins associated with clinical responses of treatment of 115 hepatocellular carcinoma (HCC) patients with the tyrosine kinase inhibitor sorafenib. Of 124 candidate proteins, LMM approaches provided 6 results varying in significance, while LR-SAM, by contrast, yielded 18 significant results that were quite reproducibly consistent.

Conclusion

As exemplified by an application to HCC data set, LR-SAM more effectively identified proteins associated with clinical responses of treatment than LMM did.

A systematic evaluation of normalization methods in quantitative label-free proteomics

To date, mass spectrometry (MS) data remain inherently biased as a result of reasons ranging from sample handling to differences caused by the instrumentation. Normalization is the process that aims to account for the bias and make samples more comparable. The selection of a proper normalization method is a pivotal task for the reliability of the downstream analysis and results. Many normalization methods commonly used in proteomics have been adapted from the DNA microarray techniques. Previous studies comparing normalization methods in proteomics have focused mainly on intragroup variation. In this study, several popular and widely used normalization methods representing different strategies in normalization are evaluated using three spike-in and one experimental mouse label-free proteomic data sets. The normalization methods are evaluated in terms of their ability to reduce variation between technical replicates, their effect on differential expression analysis and their effect on the estimation of logarithmic fold changes. Additionally, we examined whether normalizing the whole data globally or in segments for the differential expression analysis has an effect on the performance of the normalization methods. We found that variance stabilization normalization (Vsn) reduced variation the most between technical replicates in all examined data sets. Vsn also performed consistently well in the differential expression analysis. Linear regression normalization and local regression normalization performed also systematically well. Finally, we discuss the choice of a normalization method and some qualities of a suitable normalization method in the light of the results of our evaluation.

Comparison of methods to detect differentially expressed genes between single-cell populations

We compared five statistical methods to detect differentially expressed genes between two distinct single-cell populations. Currently, it remains unclear whether differential expression methods developed originally for conventional bulk RNA-seq data can also be applied to single-cell RNA-seq data analysis. Our results in three diverse comparison settings showed marked differences between the different methods in terms of the number of detections as well as their sensitivity and specificity. They, however, did not reveal systematic benefits of the currently available single-cell-specific methods. Instead, our previously introduced reproducibility-optimization method showed good performance in all comparison settings without any single-cell-specific modifications.

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.

In-depth method assessments of differentially expressed protein detection for shotgun proteomics data with missing values

Considering as one of the major goals in quantitative proteomics, detection of the differentially expressed proteins (DEPs) plays an important role in biomarker selection and clinical diagnostics. There have been plenty of algorithms and tools focusing on DEP detection in proteomics research. However, due to the different application scopes of these methods, and various kinds of experiment designs, it is not very apparent about the best choice for large-scale proteomics data analyses. Moreover, given the fact that proteomics data usually contain high percentage of missing values (MVs), but few replicates, a systematic evaluation of the DEP detection methods combined with the MV imputation methods is essential and urgent. Here, we analyzed a total of four representative imputation methods and five DEP methods on different experimental and simulated datasets. The results showed that (i) MV imputation could not always improve the performances of DEP detection methods and the imputation effects differed in the missing value percentages; (ii) the DEP detection methods had different statistical powers affected by the percentage of MVs. Two statistical methods (i.e. the empirical Bayesian random censoring threshold model, and the significance analysis of microarray) performed better than the other evaluated methods in terms of accuracy and sensitivity.

ROTS: An R package for reproducibility-optimized statistical testing

Differential expression analysis is one of the most common types of analyses performed on various biological data (e.g. RNA-seq or mass spectrometry proteomics). It is the process that detects features, such as genes or proteins, showing statistically significant differences between the sample groups under comparison. A major challenge in the analysis is the choice of an appropriate test statistic, as different statistics have been shown to perform well in different datasets. To this end, the reproducibility-optimized test statistic (ROTS) adjusts a modified t-statistic according to the inherent properties of the data and provides a ranking of the features based on their statistical evidence for differential expression between two groups. ROTS has already been successfully applied in a range of different studies from transcriptomics to proteomics, showing competitive performance against other state-of-the-art methods. To promote its widespread use, we introduce here a Bioconductor R package for performing ROTS analysis conveniently on different types of omics data. To illustrate the benefits of ROTS in various applications, we present three case studies, involving proteomics and RNA-seq data from public repositories, including both bulk and single cell data. The package is freely available from Bioconductor (https://www.bioconductor.org/packages/ROTS).

Statistical Analysis of ATM-Dependent Signaling in Quantitative Mass Spectrometry Phosphoproteomics

Ataxia-telangiectasia mutated (ATM) is a serine/threonine protein kinase, which when perturbed is associated with modified protein signaling that ultimately leads to a range of neurological and DNA repair defects. Recent advances in phospho-proteomics coupled with high-resolution mass-spectrometry provide new opportunities to dissect signaling pathways that ATM utilize under a number of conditions. This chapter begins by providing a brief overview of ATM function, its various regulatory roles and then leads into a workflow focused on the use of the statistical programming language R, together with code, for the identification of ATM-dependent substrates in the cytoplasm. This chapter cannot cover statistical properties in depth nor the range of possible methods in great detail, but instead aims to equip researchers with a set of tools to perform analysis between two conditions through examples with R functions.

Thousand and one ways to quantify and compare protein abundances in label-free bottom-up proteomics

How to process and analyze MS data to quantify and statistically compare protein abundances in bottom-up proteomics has been an open debate for nearly fifteen years. Two main approaches are generally used: the first is based on spectral data generated during the process of identification (e.g. peptide counting, spectral counting), while the second makes use of extracted ion currents to quantify chromatographic peaks and infer protein abundances based on peptide quantification. These two approaches actually refer to multiple methods which have been developed during the last decade, but were submitted to deep evaluations only recently. In this paper, we compiled these different methods as exhaustively as possible. We also summarized the way they address the different problems raised by bottom-up protein quantification such as normalization, the presence of shared peptides, unequal peptide measurability and missing data. This article is part of a Special Issue entitled: Plant Proteomics--a bridge between fundamental processes and crop production, edited by Dr. Hans-Peter Mock.

Liver lipid metabolism is altered by increased circulating estrogen to androgen ratio in male mouse

Estrogens are suggested to lower the risk of developing metabolic syndrome in both sexes. In this study, we investigated how the increased circulating estrogen-to-androgen ratio (E/A) alters liver lipid metabolism in males. The cytochrome P450 aromatase (P450arom) is an enzyme converting androgens to estrogens. Male mice overexpressing human aromatase enzyme (AROM+ mice), and thus have high circulating E/A, were used as a model in this study. Proteomics and gene expression analyses indicated an increase in the peroxisomal β-oxidation in the liver of AROM+ mice as compared with their wild type littermates. Correspondingly, metabolomic analysis revealed a decrease in the amount of phosphatidylcholines with long-chain fatty acids in the plasma. With interest we noted that the expression of Cyp4a12a enzyme, which specifically metabolizes arachidonic acid (AA) to 20-hydroxy AA, was dramatically decreased in the AROM+ liver. As a consequence, increased amounts of phospholipids having AA as a fatty acid tail were detected in the plasma of the AROM+ mice. Overall, these observations demonstrate that high circulating E/A in males is linked to indicators of higher peroxisomal β-oxidation and lower AA metabolism in the liver. Furthermore, the plasma phospholipid profile reflects the changes in the liver lipid metabolism.