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

Abstract

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.

TcellSubC: An Atlas of the Subcellular Proteome of Human T Cells

“…To identify potential stimulation-induced translocations we analyzed the MS data using paired analysis in Limma (24) and the DeqMS R-package (25). The list of overlapping proteins from the 3 sets (3 donors, n = 7,122 proteins), with full quantitation in all channels ( n = 6,572 proteins was imported into and used as a matrix.…”

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

“… 58 DeqMs, designed for TMT, is based on log2 normalization and LIMMA. 59 Isobar is an R package designed for TMT and iTRAQ. 60 Extensive benchmarking is still necessary to obtain consensus pipelines used by the proteomics community.…”

Identification of a Proteomic Signature of Senescence in Primary Human Mammary Epithelial Cells

“…The optimized number of neighbors was determined to be n = 10. Protein abundance log2 ratios and statistical significance were calculated using DEqMS in R software (Zhu et al 2020). Briefly, peptide sequences were aggregated into protein log2 ratios by the median sweeping method: raw intensity values were log2 transformed, the median of log2 intensity was subtracted for each PSM, and then for each protein, the relative log2 ratio was calculated as the median of log2 ratio of the PSMs assigned to that protein.…”