MSqRob Takes the Missing Hurdle: Uniting Intensity- and Count-Based Proteomics

How mass spectrometry can be exploited to study AMPK

AMP-activated protein kinase (AMPK) is a key regulator of metabolism and a recognised target for the treatment of metabolic diseases such as Type 2 diabetes (T2D). Here, we review how mass spectrometry (MS) can be used to study short-term control by AMPK via protein phosphorylation and long-term control due to changes in protein expression. We discuss how MS can quantify AMPK subunit levels in tissues from different species. We propose hydrogen-deuterium exchange (HDX)-MS to investigate molecular mechanisms of AMPK activation and thermoproteomic profiling (TPP) to assess off-target effects of pharmacological AMPK activators/inhibitors. Lastly, because large MS data sets are generated, we consider different approaches that can be used for their interpretation.

An MSstats workflow for detecting differentially abundant proteins in large-scale data-independent acquisition mass spectrometry experiments with FragPipe processing

Technological advances in mass spectrometry and proteomics have made it possible to perform larger-scale and more-complex experiments. The volume and complexity of the resulting data create major challenges for downstream analysis. In particular, next-generation data-independent acquisition (DIA) experiments enable wider proteome coverage than more traditional targeted approaches but require computational workflows that can manage much larger datasets and identify peptide sequences from complex and overlapping spectral features. Data-processing tools such as FragPipe, DIA-NN and Spectronaut have undergone substantial improvements to process spectral features in a reasonable time. Statistical analysis tools are needed to draw meaningful comparisons between experimental samples, but these tools were also originally designed with smaller datasets in mind. This protocol describes an updated version of MSstats that has been adapted to be compatible with large-scale DIA experiments. A very large DIA experiment, processed with FragPipe, is used as an example to demonstrate different MSstats workflows. The choice of workflow depends on the user's computational resources. For datasets that are too large to fit into a standard computer's memory, we demonstrate the use of MSstatsBig, a companion R package to MSstats. The protocol also highlights key decisions that have a major effect on both the results and the processing time of the analysis. The MSstats processing can be expected to take 1-3 h depending on the usage of MSstatsBig. The protocol can be run in the point-and-click graphical user interface MSstatsShiny or implemented with minimal coding expertise in R.

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.

Collagen IV assembly is influenced by fluid flow in kidney cell-derived matrices

Kidney podocytes and endothelial cells assemble a complex and dynamic basement membrane that is essential for kidney filtration. Whilst many components of this specialised matrix are known, the influence of fluid flow on its assembly and organisation remains poorly understood. Using the coculture of podocytes and glomerular endothelial cells in a low-shear stress, high-flow bioreactor, we investigated the effect of laminar fluid flow on the composition and assembly of cell-derived matrix. With immunofluorescence and matrix image analysis we found flow-mediated remodelling of collagen IV. Using proteomic analysis of the cell-derived matrix we identified changes in both abundance and composition of matrix proteins under flow, including the collagen-modifying enzyme, prolyl 4-hydroxylase (P4HA1). To track collagen IV assembly, we used CRISPR-Cas9 to knock in the luminescent marker HiBiT to the endogenous COL4A2 gene in podocytes. With this system, we found that collagen IV was secreted and accumulated consistently under both static and flow conditions. However knockdown of P4HA1 in podocytes led to a reduction in the secretion of collagen IV and this was more pronounced under flow. Together, this work demonstrates the effect of fluid flow on the composition, modification, and organisation of kidney cell-derived matrix and provides an in vitro system for investigating flow-induced matrix alteration in the context of kidney development and disease.

Kidney resident macrophages have distinct subsets and multifunctional roles

BACKGROUND

The kidney contains distinct glomerular and tubulointerstitial compartments with diverse cell types and extracellular matrix components. The role of immune cells in glomerular environment is crucial for dampening inflammation and maintaining homeostasis. Macrophages are innate immune cells that are influenced by their tissue microenvironment. However, the multifunctional role of kidney macrophages remains unclear.

METHODS

Flow and imaging cytometry were used to determine the relative expression of CD81 and CX3CR1 (C-X3-C motif chemokine receptor 1) in kidney macrophages. Monocyte replenishment was assessed in Cx3cr1CreER X R26-yfp-reporter and shielded chimeric mice. Bulk RNA-sequencing and mass spectrometry-based proteomics were performed on isolated kidney macrophages from wild type and Col4a5-/- (Alport) mice. RNAscope was used to visualize transcripts and macrophage purity in bulk RNA assessed by CIBERSORTx analyses.

RESULTS

In wild type mice we identified three distinct kidney macrophage subsets using CD81 and CX3CR1 and these subsets showed dependence on monocyte replenishment. In addition to their immune function, bulk RNA-sequencing of macrophages showed enrichment of biological processes associated with extracellular matrix. Proteomics identified collagen IV and laminins in kidney macrophages from wild type mice whilst other extracellular matrix proteins including cathepsins, ANXA2 and LAMP2 were enriched in Col4a5-/- (Alport) mice. A subset of kidney macrophages co-expressed matrix and macrophage transcripts.

CONCLUSIONS

We identified CD81 and CX3CR1 positive kidney macrophage subsets with distinct dependence for monocyte replenishment. Multiomic analysis demonstrated that these cells have diverse functions that underscore the importance of macrophages in kidney health and disease.

msqrob2PTM: Differential Abundance and Differential Usage Analysis of MS-Based Proteomics Data at the Posttranslational Modification and Peptidoform Level

In the era of open-modification search engines, more posttranslational modifications than ever can be detected by LC-MS/MS-based proteomics. This development can switch proteomics research into a higher gear, as PTMs are key in many cellular pathways important in cell proliferation, migration, metastasis, and aging. However, despite these advances in modification identification, statistical methods for PTM-level quantification and differential analysis have yet to catch up. This absence can partly be explained by statistical challenges inherent to the data, such as the confounding of PTM intensities with its parent protein abundance. Therefore, we have developed msqrob2PTM, a new workflow in the msqrob2 universe capable of differential abundance analysis at the PTM and at the peptidoform level. The latter is important for validating PTMs found as significantly differential. Indeed, as our method can deal with multiple PTMs per peptidoform, there is a possibility that significant PTMs stem from one significant peptidoform carrying another PTM, hinting that it might be the other PTM driving the perceived differential abundance. Our workflows can flag both differential peptidoform abundance (DPA) and differential peptidoform usage (DPU). This enables a distinction between direct assessment of differential abundance of peptidoforms (DPA) and differences in the relative usage of peptidoforms corrected for corresponding protein abundances (DPU). For DPA, we directly model the log2-transformed peptidoform intensities, while for DPU, we correct for parent protein abundance by an intermediate normalization step which calculates the log2-ratio of the peptidoform intensities to their summarized parent protein intensities. We demonstrated the utility and performance of msqrob2PTM by applying it to datasets with known ground truth, as well as to biological PTM-rich datasets. Our results show that msqrob2PTM is on par with, or surpassing the performance of, the current state-of-the-art methods. Moreover, msqrob2PTM is currently unique in providing output at the peptidoform level.

Bioinformatics Pipeline for Processing Single-Cell Data

Single-cell proteomics can offer valuable insights into dynamic cellular interactions, but identifying proteins at this level is challenging due to their low abundance. In this chapter, we present a state-of-the-art bioinformatics pipeline for single-cell proteomics that combines the search engine Sage (via SearchGUI), identification rescoring with MS2Rescore, quantification through FlashLFQ, and differential expression analysis using MSqRob2. MS2Rescore leverages LC-MS/MS behavior predictors, such as MS2PIP and DeepLC, to recalibrate scores with Percolator or mokapot. Combining these tools into a unified pipeline, this approach improves the detection of low-abundance peptides, resulting in increased identifications while maintaining stringent FDR thresholds.

Standardized Workflow for Mass-Spectrometry-Based Single-Cell Proteomics Data Processing and Analysis Using the scp Package

Mass-spectrometry (MS)-based single-cell proteomics (SCP) explores cellular heterogeneity by focusing on the functional effectors of the cells-proteins. However, extracting meaningful biological information from MS data is far from trivial, especially with single cells. Currently, data analysis workflows are substantially different from one research team to another. Moreover, it is difficult to evaluate pipelines as ground truths are missing. Our team has developed the R/Bioconductor package called scp to provide a standardized framework for SCP data analysis. It relies on the widely used QFeatures and SingleCellExperiment data structures. In addition, we used a design containing cell lines mixed in known proportions to generate controlled variability for data analysis benchmarking. In this chapter, we provide a flexible data analysis protocol for SCP data using the scp package together with comprehensive explanations at each step of the processing. Our main steps are quality control on the feature and cell level, aggregation of the raw data into peptides and proteins, normalization, and batch correction. We validate our workflow using our ground truth data set. We illustrate how to use this modular, standardized framework and highlight some crucial steps.

lesSDRF is more: maximizing the value of proteomics data through streamlined metadata annotation

Public proteomics data often lack essential metadata, limiting its potential. To address this, we present lesSDRF, a tool to simplify the process of metadata annotation, thereby ensuring that data leave a lasting, impactful legacy well beyond its initial publication.

A Bioconductor workflow for processing, evaluating, and interpreting expression proteomics data

Background: Expression proteomics involves the global evaluation of protein abundances within a system. In turn, differential expression analysis can be used to investigate changes in protein abundance upon perturbation to such a system. Methods: Here, we provide a workflow for the processing, analysis and interpretation of quantitative mass spectrometry-based expression proteomics data. This workflow utilizes open-source R software packages from the Bioconductor project and guides users end-to-end and step-by-step through every stage of the analyses. As a use-case we generated expression proteomics data from HEK293 cells with and without a treatment. Of note, the experiment included cellular proteins labelled using tandem mass tag (TMT) technology and secreted proteins quantified using label-free quantitation (LFQ). Results: The workflow explains the software infrastructure before focusing on data import, pre-processing and quality control. This is done individually for TMT and LFQ datasets. The application of statistical differential expression analysis is demonstrated, followed by interpretation via gene ontology enrichment analysis. Conclusions: A comprehensive workflow for the processing, analysis and interpretation of expression proteomics is presented. The workflow is a valuable resource for the proteomics community and specifically beginners who are at least familiar with R who wish to understand and make data-driven decisions with regards to their analyses.

Revisiting the Thorny Issue of Missing Values in Single-Cell Proteomics

Missing values are a notable challenge when analyzing mass spectrometry-based proteomics data. While the field is still actively debating the best practices, the challenge increased with the emergence of mass spectrometry-based single-cell proteomics and the dramatic increase in missing values. A popular approach to deal with missing values is to perform imputation. Imputation has several drawbacks for which alternatives exist, but currently, imputation is still a practical solution widely adopted in single-cell proteomics data analysis. This perspective discusses the advantages and drawbacks of imputation. We also highlight 5 main challenges linked to missing value management in single-cell proteomics. Future developments should aim to solve these challenges, whether it is through imputation or data modeling. The perspective concludes with recommendations for reporting missing values, for reporting methods that deal with missing values, and for proper encoding of missing values.

Acute Deletion of the Glucocorticoid Receptor in Hepatocytes Disrupts Postprandial Lipid Metabolism in Male Mice

Hepatic lipid metabolism is highly dynamic, and disruption of several circadian transcriptional regulators results in hepatic steatosis. This includes genetic disruption of the glucocorticoid receptor (GR) as the liver develops. To address the functional role of GR in the adult liver, we used an acute hepatocyte-specific GR knockout model to study temporal hepatic lipid metabolism governed by GR at several preprandial and postprandial circadian timepoints. Lipidomics analysis revealed significant temporal lipid metabolism, where GR disruption results in impaired regulation of specific triglycerides, nonesterified fatty acids, and sphingolipids. This correlates with increased number and size of lipid droplets and mildly reduced mitochondrial respiration, most noticeably in the postprandial phase. Proteomics and transcriptomics analyses suggest that dysregulated lipid metabolism originates from pronounced induced expression of enzymes involved in fatty acid synthesis, β-oxidation, and sphingolipid metabolism. Integration of GR cistromic data suggests that induced gene expression is a result of regulatory actions secondary to direct GR effects on gene transcription.

Secretory Proteomic Responses of Endometrial Epithelial Cells to Trophoblast-Derived Extracellular Vesicles

Synchronized crosstalk between the embryo and endometrium during the periconception period is integral to pregnancy establishment. Increasing evidence suggests that the exchange of extracellular vesicles (EVs) of both embryonic and endometrial origin is a critical component of embryo-maternal communication during peri-implantation. Here, we investigated whether embryonic signals in the form of EVs can modulate the endometrial epithelial cell secretome. Receptive endometrial analog RL95-2 cells were supplemented with trophoblast analog JAr cell-derived EVs, and the secretory protein changes occurring in the RL95-2 cells were analyzed using mass spectrometry. EVs of non-trophoblastic origin (HEK 293 cells) were used as the control EV source to supplement endometrial cells. Trophoblast cell-derived EVs enriched endometrial epithelial cell secretions with proteins that support embryo development, attachment, or implantation, whereas control EVs were unable to induce the same effect. The present study suggests that embryonic signals in the form of EVs may prime receptive endometrial epithelial cells to enrich their secretory proteome with critical proteomic molecules with functional importance for periconception milieu formation.

Tidyproteomics: an open-source R package and data object for quantitative proteomics post analysis and visualization

BACKGROUND

The analysis of mass spectrometry-based quantitative proteomics data can be challenging given the variety of established analysis platforms, the differences in reporting formats, and a general lack of approachable standardized post-processing analyses such as sample group statistics, quantitative variation and even data filtering. We developed tidyproteomics to facilitate basic analysis, improve data interoperability and potentially ease the integration of new processing algorithms, mainly through the use of a simplified data-object.

RESULTS

The R package tidyproteomics was developed as both a framework for standardizing quantitative proteomics data and a platform for analysis workflows, containing discrete functions that can be connected end-to-end, thus making it easier to define complex analyses by breaking them into small stepwise units. Additionally, as with any analysis workflow, choices made during analysis can have large impacts on the results and as such, tidyproteomics allows researchers to string each function together in any order, select from a variety of options and in some cases develop and incorporate custom algorithms.

CONCLUSIONS

Tidyproteomics aims to simplify data exploration from multiple platforms, provide control over individual functions and analysis order, and serve as a tool to assemble complex repeatable processing workflows in a logical flow. Datasets in tidyproteomics are easy to work with, have a structure that allows for biological annotations to be added, and come with a framework for developing additional analysis tools. The consistent data structure and accessible analysis and plotting tools also offers a way for researchers to save time on mundane data manipulation tasks.

LFQ-Based Peptide and Protein Intensity Differential Expression Analysis

Testing for significant differences in quantities at the protein level is a common goal of many LFQ-based mass spectrometry proteomics experiments. Starting from a table of protein and/or peptide quantities from a given proteomics quantification software, many tools and R packages exist to perform the final tasks of imputation, summarization, normalization, and statistical testing. To evaluate the effects of packages and settings in their substeps on the final list of significant proteins, we studied several packages on three public data sets with known expected protein fold changes. We found that the results between packages and even across different parameters of the same package can vary significantly. In addition to usability aspects and feature/compatibility lists of different packages, this paper highlights sensitivity and specificity trade-offs that come with specific packages and settings.

MSstats Version 4.0: Statistical Analyses of Quantitative Mass Spectrometry-Based Proteomic Experiments with Chromatography-Based Quantification at Scale

The MSstats R-Bioconductor family of packages is widely used for statistical analyses of quantitative bottom-up mass spectrometry-based proteomic experiments to detect differentially abundant proteins. It is applicable to a variety of experimental designs and data acquisition strategies and is compatible with many data processing tools used to identify and quantify spectral features. In the face of ever-increasing complexities of experiments and data processing strategies, the core package of the family, with the same name MSstats, has undergone a series of substantial updates. Its new version MSstats v4.0 improves the usability, versatility, and accuracy of statistical methodology, and the usage of computational resources. New converters integrate the output of upstream processing tools directly with MSstats, requiring less manual work by the user. The package's statistical models have been updated to a more robust workflow. Finally, MSstats' code has been substantially refactored to improve memory use and computation speed. Here we detail these updates, highlighting methodological differences between the new and old versions. An empirical comparison of MSstats v4.0 to its previous implementations, as well as to the packages MSqRob and DEqMS, on controlled mixtures and biological experiments demonstrated a stronger performance and better usability of MSstats v4.0 as compared to existing methods.

Mass spectrometry-based retina proteomics

A subfield of neuroproteomics, retina proteomics has experienced a transformative growth since its inception due to methodological advances in enabling chemical, biochemical, and molecular biology techniques. This review focuses on mass spectrometry's contributions to facilitate mammalian and avian retina proteomics to catalog and quantify retinal protein expressions, determine their posttranslational modifications, as well as its applications to study the proteome of the retina in the context of biology, health and diseases, and therapy developments.

Combining time-resolved transcriptomics and proteomics data for Adverse Outcome Pathway refinement in ecotoxicology

Conventional Environmental Risk Assessment (ERA) of pesticide pollution is based on soil concentrations and apical endpoints, such as the reproduction of test organisms, but has traditionally disregarded information along the organismal response cascade leading to an adverse outcome. The Adverse Outcome Pathway (AOP) framework includes response information at any level of biological organization, providing opportunities to use intermediate responses as a predictive read-out for adverse outcomes instead. Transcriptomic and proteomic data can provide thousands of data points on the response to toxic exposure. Combining multiple omics data types is necessary for a comprehensive overview of the response cascade and, therefore, AOP development. However, it is unclear if transcript and protein responses are synchronized in time or time lagged. To understand if analysis of multi-omics data obtained at the same timepoint reveal one synchronized response cascade, we studied time-resolved shifts in gene transcript and protein abundance in the springtail Folsomia candida, a soil ecotoxicological model, after exposure to the neonicotinoid insecticide imidacloprid. We analyzed transcriptome and proteome data every 12 h up to 72 h after onset of exposure. The most pronounced shift in both transcript and protein abundances was observed after 48 h exposure. Moreover, cross-correlation analyses indicate that most genes displayed the highest correlation between transcript and protein abundances without a time-lag. This demonstrates that a combined analysis of transcriptomic and proteomic data from the same time-point can be used for AOP improvement. This data will promote the development of biomarkers for the presence of neonicotinoid insecticides or chemicals with a similar mechanism of action in soils.

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.

The Current State of Single-Cell Proteomics Data Analysis

Sound data analysis is essential to retrieve meaningful biological information from single-cell proteomics experiments. This analysis is carried out by computational methods that are assembled into workflows, and their implementations influence the conclusions that can be drawn from the data. In this work, we explore and compare the computational workflows that have been used over the last four years and identify a profound lack of consensus on how to analyze single-cell proteomics data. We highlight the need for benchmarking of computational workflows and standardization of computational tools and data, as well as carefully designed experiments. Finally, we cover the current standardization efforts that aim to fill the gap, list the remaining missing pieces, and conclude with lessons learned from the replication of published single-cell proteomics analyses. © 2023 Wiley Periodicals LLC.

Proteomic Comparison of Three Wild-Type Pseudorabies Virus Strains and the Attenuated Bartha Strain Reveals Reduced Incorporation of Several Tegument Proteins in Bartha Virions

Pseudorabies virus (PRV) is a member of the alphaherpesvirus subfamily and the causative agent of Aujeszky's disease in pigs. Driven by the large economic losses associated with PRV infection, several vaccines and vaccine programs have been developed. To this day, the attenuated Bartha strain, generated by serial passaging, represents the golden standard for PRV vaccination. However, a proteomic comparison of the Bartha virion to wild-type (WT) PRV virions is lacking. Here, we present a comprehensive mass spectrometry-based proteome comparison of the attenuated Bartha strain and three commonly used WT PRV strains: Becker, Kaplan, and NIA3. We report the detection of 40 structural and 14 presumed nonstructural proteins through a combination of data-dependent and data-independent acquisition. Interstrain comparisons revealed that packaging of the capsid and most envelope proteins is largely comparable in-between all four strains, except for the envelope protein pUL56, which is less abundant in Bartha virions. However, distinct differences were noted for several tegument proteins. Most strikingly, we noted a severely reduced incorporation of the tegument proteins IE180, VP11/12, pUS3, VP22, pUL41, pUS1, and pUL40 in Bartha virions. Moreover, and likely as a consequence, we also observed that Bartha virions are on average smaller and more icosahedral compared to WT virions. Finally, we detected at least 28 host proteins that were previously described in PRV virions and noticed considerable strain-specific differences with regard to host proteins, arguing that the potential role of packaged host proteins in PRV replication and spread should be further explored. IMPORTANCE The pseudorabies virus (PRV) vaccine strain Bartha-an attenuated strain created by serial passaging-represents an exceptional success story in alphaherpesvirus vaccination. Here, we used mass spectrometry to analyze the Bartha virion composition in comparison to three established WT PRV strains. Many viral tegument proteins that are considered nonessential for viral morphogenesis were drastically less abundant in Bartha virions compared to WT virions. Interestingly, many of the proteins that are less incorporated in Bartha participate in immune evasion strategies of alphaherpesviruses. In addition, we observed a reduced size and more icosahedral morphology of the Bartha virions compared to WT PRV. Given that the Bartha vaccine strain elicits potent immune responses, our findings here suggest that differences in protein packaging may contribute to its immunogenicity. Further exploration of these observations could aid the development of efficacious vaccines against other alphaherpesvirus vaccines such as HSV-1/2 or EHV-1.

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

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

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

Proteomic Assessment of C57BL/6 Hippocampi after Non-Selective Pharmacological Inhibition of Nitric Oxide Synthase Activity: Implications of Seizure-like Neuronal Hyperexcitability Followed by Tauopathy

Nitric oxide (NO) is a small gaseous signaling molecule responsible for maintaining homeostasis in a myriad of tissues and molecular pathways in neurology and the cardiovasculature. In recent years, there has been increasing interest in the potential interaction between arterial stiffness (AS), an independent cardiovascular risk factor, and neurodegenerative syndromes given increasingly epidemiological study reports. For this reason, we previously investigated the mechanistic convergence between AS and neurodegeneration via the progressive non-selective inhibition of all nitric oxide synthase (NOS) isoforms with N(G)-nitro-L-arginine methyl ester (L-NAME) in C57BL/6 mice. Our previous results showed progressively increased AS in vivo and impaired visuospatial learning and memory in L-NAME-treated C57BL/6 mice. In the current study, we sought to further investigate the progressive molecular signatures in hippocampal tissue via LC–MS/MS proteomic analysis. Our data implicate mitochondrial dysfunction due to progressive L-NAME treatment. Two weeks of L-NAME treatment implicates altered G-protein-coupled-receptor signaling in the nerve synapse and associated presence of seizures and altered emotional behavior. Furthermore, molecular signatures implicate the cerebral presence of seizure-related hyperexcitability after short-term (8 weeks) treatment followed by ribosomal dysfunction and tauopathy after long-term (16 weeks) treatment.

A combined test for feature selection on sparse metaproteomics data-an alternative to missing value imputation

One of the difficulties encountered in the statistical analysis of metaproteomics data is the high proportion of missing values, which are usually treated by imputation. Nevertheless, imputation methods are based on restrictive assumptions regarding missingness mechanisms, namely "at random" or "not at random". To circumvent these limitations in the context of feature selection in a multi-class comparison, we propose a univariate selection method that combines a test of association between missingness and classes, and a test for difference of observed intensities between classes. This approach implicitly handles both missingness mechanisms. We performed a quantitative and qualitative comparison of our procedure with imputation-based feature selection methods on two experimental data sets, as well as simulated data with various scenarios regarding the missingness mechanisms and the nature of the difference of expression (differential intensity or differential presence). Whereas we observed similar performances in terms of prediction on the experimental data set, the feature ranking and selection from various imputation-based methods were strongly divergent. We showed that the combined test reaches a compromise by correlating reasonably with other methods, and remains efficient in all simulated scenarios unlike imputation-based feature selection methods.

Challenges and Opportunities for Bayesian Statistics in Proteomics

Proteomics is a data-rich science with complex experimental designs and an intricate measurement process. To obtain insights from the large data sets produced, statistical methods, including machine learning, are routinely applied. For a quantity of interest, many of these approaches only produce a point estimate, such as a mean, leaving little room for more nuanced interpretations. By contrast, Bayesian statistics allows quantification of uncertainty through the use of probability distributions. These probability distributions enable scientists to ask complex questions of their proteomics data. Bayesian statistics also offers a modular framework for data analysis by making dependencies between data and parameters explicit. Hence, specifying complex hierarchies of parameter dependencies is straightforward in the Bayesian framework. This allows us to use a statistical methodology which equals, rather than neglects, the sophistication of experimental design and instrumentation present in proteomics. Here, we review Bayesian methods applied to proteomics, demonstrating their potential power, alongside the challenges posed by adopting this new statistical framework. To illustrate our review, we give a walk-through of the development of a Bayesian model for dynamic organic orthogonal phase-separation (OOPS) data.

PBLMM: Peptide-based linear mixed models for differential expression analysis of shotgun proteomics data

Here, we present a peptide-based linear mixed models tool-PBLMM, a standalone desktop application for differential expression analysis of proteomics data. We also provide a Python package that allows streamlined data analysis workflows implementing the PBLMM algorithm. PBLMM is easy to use without scripting experience and calculates differential expression by peptide-based linear mixed regression models. We show that peptide-based models outperform classical methods of statistical inference of differentially expressed proteins. In addition, PBLMM exhibits superior statistical power in situations of low effect size and/or low sample size. Taken together our tool provides an easy-to-use, high-statistical-power method to infer differentially expressed proteins from proteomics data.

Comparative assessment and novel strategy on methods for imputing proteomics data

Missing values are a major issue in quantitative proteomics analysis. While many methods have been developed for imputing missing values in high-throughput proteomics data, a comparative assessment of imputation accuracy remains inconclusive, mainly because mechanisms contributing to true missing values are complex and existing evaluation methodologies are imperfect. Moreover, few studies have provided an outlook of future methodological development. We first re-evaluate the performance of eight representative methods targeting three typical missing mechanisms. These methods are compared on both simulated and masked missing values embedded within real proteomics datasets, and performance is evaluated using three quantitative measures. We then introduce fused regularization matrix factorization, a low-rank global matrix factorization framework, capable of integrating local similarity derived from additional data types. We also explore a biologically-inspired latent variable modeling strategy—convex analysis of mixtures—for missing value imputation and present preliminary experimental results. While some winners emerged from our comparative assessment, the evaluation is intrinsically imperfect because performance is evaluated indirectly on artificial missing or masked values not authentic missing values. Nevertheless, we show that our fused regularization matrix factorization provides a novel incorporation of external and local information, and the exploratory implementation of convex analysis of mixtures presents a biologically plausible new approach.

Replication of single-cell proteomics data reveals important computational challenges

INTRODUCTION

Mass spectrometry-based proteomics is actively embracing quantitative, single-cell level analyses. Indeed, recent advances in sample preparation and mass spectrometry (MS) have enabled the emergence of quantitative MS-based single-cell proteomics (SCP). While exciting and promising, SCP still has many rough edges. The current analysis workflows are custom and built from scratch. The field is therefore craving for standardized software that promotes principled and reproducible SCP data analyses.

AREAS COVERED

This special report is the first step toward the formalization and standardization of SCP data analysis. scp, the software that accompanies this work, successfully replicates one of the landmark SCP studies and is applicable to other experiments and designs. We created a repository containing the replicated workflow with comprehensive documentation in order to favor further dissemination and improvements of SCP data analyses.

EXPERT OPINION

Replicating SCP data analyses uncovers important challenges in SCP data analysis. We describe two such challenges in detail: batch correction and data missingness. We provide the current state-of-the-art and illustrate the associated limitations. We also highlight the intimate dependence that exists between batch effects and data missingness and offer avenues for dealing with these exciting challenges.

Comprehensive histone epigenetics: A mass spectrometry based screening assay to measure epigenetic toxicity

Evidence of the involvement of epigenetics in pathologies such as cancer, diabetes, and neurodegeneration has increased global interest in epigenetic modifications. For nearly thirty years, it has been known that cancer cells exhibit abnormal DNA methylation patterns. In contrast, the large-scale analysis of histone post-translational modifications (hPTMs) has lagged behind because classically, histone modification analysis has relied on site specific antibody-based techniques. Mass spectrometry (MS) is a technique that holds the promise to picture the histone code comprehensively in a single experiment. Therefore, we developed an MS-based method that is capable of tracking all possible hPTMs in an untargeted approach. In this way, trends in single and combinatorial hPTMs can be reported and enable prediction of the epigenetic toxicity of compounds. Moreover, this method is based on the use of human cells to provide preliminary data, thereby omitting the need to sacrifice laboratory animals. Improving the workflow and the user-friendliness in order to become a high throughput, easily applicable, toxicological screening assay is an ongoing effort. Still, this novel toxicoepigenetic assay and the data it generates holds great potential for, among others, pharmaceutical industry, food science, clinical diagnostics and, environmental toxicity screening. •There is a growing interest in epigenetic modifications, and more specifically in histone post-translational modifications (hPTMs).•We describe an MS-based workflow that is capable of tracking all possible hPTMs in an untargeted approach that makes use of human cells.•Improving the workflow and the user-friendliness in order to become a high throughput, easily applicable, toxicological screening assay is an ongoing effort.