Evaluating False Transfer Rates from the Match-between-Runs Algorithm with a Two-Proteome Model

Enhanced feature matching in single-cell proteomics characterizes IFN-γ response and co-existence of cell states

Proteome analysis by data-independent acquisition (DIA) has become a powerful approach to obtain deep proteome coverage, and has gained recent traction for label-free analysis of single cells. However, optimal experimental design for DIA-based single-cell proteomics has not been fully explored, and performance metrics of subsequent data analysis tools remain to be evaluated. Therefore, we here formalize and comprehensively evaluate a DIA data analysis strategy that exploits the co-analysis of low-input samples with a so-called matching enhancer (ME) of higher input, to increase sensitivity, proteome coverage, and data completeness. We assess the matching specificity of DIA-ME by a two-proteome model, and demonstrate that false discovery and false transfer are maintained at low levels when using DIA-NN software, while preserving quantification accuracy. We apply DIA-ME to investigate the proteome response of U-2 OS cells to interferon gamma (IFN-γ) in single cells, and recapitulate the time-resolved induction of IFN-γ response proteins as observed in bulk material. Moreover, we uncover co- and anti-correlating patterns of protein expression within the same cell, indicating mutually exclusive protein modules and the co-existence of different cell states. Collectively our data show that DIA-ME is a powerful, scalable, and easy-to-implement strategy for single-cell proteomics.

Interrogating data-independent acquisition LC-MS/MS for affinity proteomics

Data-Independent Acquisition (DIA) LC-MS/MS is an attractive partner for co-immunoprecipitation (co-IP) and affinity proteomics in general. Reducing the variability of quantitation by DIA could increase the statistical contrast for detecting specific interactors versus what has been achieved in Data-Dependent Acquisition (DDA). By interrogating affinity proteomes featuring both DDA and DIA experiments, we sought to evaluate the spectral libraries, the missingness of protein quantity tables, and the CV of protein quantities in six studies representing three different instrument manufacturers. We examined four contemporary bioinformatics workflows for DIA: FragPipe, DIA-NN, Spectronaut, and MaxQuant. We determined that (1) identifying spectral libraries directly from DIA experiments works well enough that separate DDA experiments do not produce larger spectral libraries when given equivalent instrument time; (2) experiments involving mock pull-downs or IgG controls may feature such indistinct signals that contemporary software will struggle to quantify them; (3) measured CV values were well controlled by Spectronaut and DIA-NN (and FragPipe, which implements DIA-NN for the quantitation step); and (4) when FragPipe builds spectral libraries and quantifies proteins from DIA experiments rather than performing both operations in DDA experiments, the DIA route results in a larger number of proteins quantified without missing values as well as lower CV for measured protein quantities.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42485-024-00166-4.

The impact of cooking and burial on proteins: a characterisation of experimental foodcrusts and ceramics

Foodcrusts have received relatively little attention in the burgeoning field of proteomic analysis of ancient cuisine. We remain ignorant of how cooking and burial impact protein survival, and crucially, the extent to which the extractome reflects the composition of input ingredients. Therefore, through experimental analogues, we explore the extent of protein survival in unburied and buried foodcrusts and ceramics using 'typical' Mesolithic ingredients (red deer, Atlantic salmon and sweet chestnut). We then explore a number of physicochemical properties theorised to aid protein preservation. The results reveal that proteins were much more likely to be detected in foodcrusts than ceramics using the methodology employed, that input ingredient strongly influences protein preservation, and that degradation is not universal nor linear between proteins, indicating that multiple protein physicochemical properties are at play. While certain properties such as hydrophobicity apparently aid protein preservation, none single-handedly explain why particular proteins/peptides survive in buried foodcrusts: this complex interplay requires further investigation. The findings demonstrate that proteins indicative of the input ingredient can be identifiable in foodcrust, but that the full proteome is unlikely to preserve. While this shows promise for the survival of proteins in archaeological foodcrust, further research is needed to accurately interpret foodcrust extractomes.

Effect of secondary and tertiary wastewater treatment methods on opioids and the subsequent environmental impact of effluent and biosolids

Opioids are widely distributed, potent prescription analgesics that are known to be diverted for illicit use. Their prevalence of use is reflected by high concentrations of parent compounds and/or metabolites found in samples collected from wastewater treatment plants. Given that treatment byproducts enter the environment through several routes, the consequences of insufficient removal by treatment methods include unwanted environmental exposure and potential to disrupt ecosystems. Activated sludge treatment has been widely investigated for a large suite of prescription opioids but the same cannot be said for UV and chlorination. Additionally, the biosolid cycle of opioids has been overlooked previously. This study aimed to determine the extent to which secondary and tertiary wastewater treatment methods remove opioids from influent, and the associated environmental exposure for those persistent, as well as the fate of opioids in biosolids. Membrane bioreactor treatment proved effective for natural and semi-synthetic opioids while the effect of UV treatment was negligible. Chlorination was the most effective treatment method resulting in effluent with concentrations below theoretical predicted no-effect concentration. Biosolids are not subjected to any additional biological or chemical treatment after membrane bioreactor treatment and the levels detected in biosolid used as fertiliser had several opioids at potentially hazardous concentrations, indicated by a QSAR theoretical model. This data indicates a potential issue regarding the treatment process of biosolids and reliance on chlorination for effluent treatment that should be investigated in other treatment plants.

Comprehensive Overview of Bottom-Up Proteomics Using Mass Spectrometry

Proteomics is the large scale study of protein structure and function from biological systems through protein identification and quantification. "Shotgun proteomics" or "bottom-up proteomics" is the prevailing strategy, in which proteins are hydrolyzed into peptides that are analyzed by mass spectrometry. Proteomics studies can be applied to diverse studies ranging from simple protein identification to studies of proteoforms, protein-protein interactions, protein structural alterations, absolute and relative protein quantification, post-translational modifications, and protein stability. To enable this range of different experiments, there are diverse strategies for proteome analysis. The nuances of how proteomic workflows differ may be challenging to understand for new practitioners. Here, we provide a comprehensive overview of different proteomics methods. We cover from biochemistry basics and protein extraction to biological interpretation and orthogonal validation. We expect this Review will serve as a handbook for researchers who are new to the field of bottom-up proteomics.

Automated single-cell proteomics providing sufficient proteome depth to study complex biology beyond cell type classifications

The recent technological and computational advances in mass spectrometry-based single-cell proteomics have pushed the boundaries of sensitivity and throughput. However, reproducible quantification of thousands of proteins within a single cell remains challenging. To address some of those limitations, we present a dedicated sample preparation chip, the proteoCHIP EVO 96 that directly interfaces with the Evosep One. This, in combination with the Bruker timsTOF demonstrates double the identifications without manual sample handling and the newest generation timsTOF Ultra identifies up to 4000 with an average of 3500 protein groups per single HEK-293T without a carrier or match-between runs. Our workflow spans 4 orders of magnitude, identifies over 50 E3 ubiquitin-protein ligases, and profiles key regulatory proteins upon small molecule stimulation. This study demonstrates that the proteoCHIP EVO 96-based sample preparation with the timsTOF Ultra provides sufficient proteome depth to study complex biology beyond cell-type classifications.

SEMQuant: Extending Sipros-Ensemble with Match-Between-Runs for Comprehensive Quantitative Metaproteomics

Metaproteomics, utilizing high-throughput LC-MS, offers a profound understanding of microbial communities. Quantitative metaproteomics further enriches this understanding by measuring relative protein abundance and revealing dynamic changes under different conditions. However, the challenge of missing peptide quantification persists in metaproteomics analysis, particularly in data-dependent acquisition mode, where high-intensity precursors for MS2 scans are selected. To tackle this issue, the match-between-runs (MBR) technique is used to transfer peptides between LC-MS runs. Inspired by the benefits of MBR and the need for streamlined metaproteomics data analysis, we developed SEMQuant, an end-to-end software integrating Sipros-Ensemble's robust peptide identifications with IonQuant's MBR function. The experiments show that SEMQuant consistently obtains the highest or second highest number of quantified proteins with notable precision and accuracy. This demonstrates SEMQuant's effectiveness in conducting comprehensive and accurate quantitative metaproteomics analyses across diverse datasets and highlights its potential to propel advancements in microbial community studies. SEMQuant is freely available under the GNU GPL license at https://github.com/Biocomputing-Research-Group/SEMQuant.

Optimizing Spectronaut Search Parameters to Improve Data Quality with Minimal Proteome Coverage Reductions in DIA Analyses of Heterogeneous Samples

Data-independent acquisition has seen breakthroughs that enable comprehensive proteome profiling using short gradients. As the proteome coverage continues to increase, the quality of the data generated becomes much more relevant. Using Spectronaut, we show that the default search parameters can be easily optimized to minimize the occurrence of false positives across different samples. Using an immunological infection model system to demonstrate the impact of adjusting search settings, we analyzed Mus musculus macrophages and compared their proteome to macrophages spiked withCandida albicans. This experimental system enabled the identification of "false positives" as Candida albicans peptides and proteins should not be present in the Mus musculus-only samples. We show that adjusting the search parameters reduced "false positive" identifications by 89% at the peptide and protein level, thereby considerably increasing the quality of the data. We also show that these optimized parameters incurred a moderate cost, only reducing the overall number of "true positive" identifications across each biological replicate by <6.7% at both the peptide and protein level. We believe the value of our updated search parameters extends beyond a two-organism analysis and would be of great value to any DIA experiment analyzing heterogeneous populations of cell types or tissues.

Increasing protein identifications in bottom-up proteomics of T. castaneum - Exploiting synergies of protein biochemistry and bioinformatics

Depending on the respective research question, LC-MS/MS based bottom-up proteomics poses challenges from the initial biological sample all the way to data evaluation. The focus of this study was to investigate the influence of sample preparation techniques and data analysis parameters on protein identification in Tribolium castaneum by applying free software proteomics platform Max Quant. Multidimensional protein extraction strategies in combination with electrophoretic or chromatographic off-line protein pre-fractionation were applied to enhance the spectrum of isolated proteins from T. castaneum and reduce the effect of co-elution and ion suppression effects during nano-LC-MS/MS measurements of peptides. For comprehensive data analysis, MaxQuant was used for protein identification and R for data evaluation. A wide range of parameters were evaluated to gain reproducible, reliable, and significant protein identifications. A simple phosphate buffer, pH 8, containing protease and phosphatase inhibitor cocktail and application of gentle extraction conditions were used as a first extraction step for T.castaneum proteins. Furthermore, a two-dimensional extraction procedure in combination with electrophoretic pre-fractionation of extracted proteins and subsequent in-gel digest resulted in almost 100% increase of identified proteins when compared to chromatographic fractionation as well as one-pot-analysis. The additionally identified proteins could be assigned to new molecular functions or cell compartments, emphasizing the positive effect of extended sample preparation in bottom-up proteomics. Besides the number of peptides during post-processing, MaxQuant's Match between Runs exhibited a crucial effect on the number of identified proteins. A maximum relative standard deviation of 2% must be considered for the data analysis. Our work with Tribolium castaneum larvae demonstrates that sometimes - depending on matrix and research question - more complex and time-consuming sample preparation can be advantageous for isolation and identification of additional proteins in bottom-up proteomics.

High-Resolution Proteomics Unravel a Native Functional Complex of Cav1.3, SK3, and Hyperpolarization-Activated Cyclic Nucleotide-Gated Channels in Midbrain Dopaminergic Neurons

Pacemaking activity in substantia nigra dopaminergic neurons is generated by the coordinated activity of a variety of distinct somatodendritic voltage- and calcium-gated ion channels. We investigated whether these functional interactions could arise from a common localization in macromolecular complexes where physical proximity would allow for efficient interaction and co-regulations. For that purpose, we immunopurified six ion channel proteins involved in substantia nigra neuron autonomous firing to identify their molecular interactions. The ion channels chosen as bait were Cav1.2, Cav1.3, HCN2, HCN4, Kv4.3, and SK3 channel proteins, and the methods chosen to determine interactions were co-immunoprecipitation analyzed through immunoblot and mass spectrometry as well as proximity ligation assay. A macromolecular complex composed of Cav1.3, HCN, and SK3 channels was unraveled. In addition, novel potential interactions between SK3 channels and sclerosis tuberous complex (Tsc) proteins, inhibitors of mTOR, and between HCN4 channels and the pro-degenerative protein Sarm1 were uncovered. In order to demonstrate the presence of these molecular interactions in situ, we used proximity ligation assay (PLA) imaging on midbrain slices containing the substantia nigra, and we could ascertain the presence of these protein complexes specifically in substantia nigra dopaminergic neurons. Based on the complementary functional role of the ion channels in the macromolecular complex identified, these results suggest that such tight interactions could partly underly the robustness of pacemaking in dopaminergic neurons.

Proteomic analysis of the urothelial cancer landscape

Urothelial bladder cancer (UC) has a wide tumor biological spectrum with challenging prognostic stratification and relevant therapy-associated morbidity. Most molecular classifications relate only indirectly to the therapeutically relevant protein level. We improve the pre-analytics of clinical samples for proteome analyses and characterize a cohort of 434 samples with 242 tumors and 192 paired normal mucosae covering the full range of UC. We evaluate sample-wise tumor specificity and rank biomarkers by target relevance. We identify robust proteomic subtypes with prognostic information independent from histopathological groups. In silico drug prediction suggests efficacy of several compounds hitherto not in clinical use. Both in silico and in vitro data indicate predictive value of the proteomic clusters for these drugs. We underline that proteomics is relevant for personalized oncology and provide abundance and tumor specificity data for a large part of the UC proteome ( www.cancerproteins.org ).

A Workflow for Improved Analysis of Cross-Linking Mass Spectrometry Data Integrating Parallel Accumulation-Serial Fragmentation with MeroX and Skyline

Cross-linking mass spectrometry (XL-MS) has evolved into a pivotal technique for probing protein interactions. This study describes the implementation of Parallel Accumulation-Serial Fragmentation (PASEF) on timsTOF instruments, enhancing the detection and analysis of protein interactions by XL-MS. Addressing the challenges in XL-MS, such as the interpretation of complex spectra, low abundant cross-linked peptides, and a data acquisition bias, our current study integrates a peptide-centric approach for the analysis of XL-MS data and presents the foundation for integrating data-independent acquisition (DIA) in XL-MS with a vendor-neutral and open-source platform. A novel workflow is described for processing data-dependent acquisition (DDA) of PASEF-derived information. For this, software by Bruker Daltonics is used, enabling the conversion of these data into a format that is compatible with MeroX and Skyline software tools. Our approach significantly improves the identification of cross-linked products from complex mixtures, allowing the XL-MS community to overcome current analytical limitations.

Varroa destructor parasitism and Deformed wing virus infection in honey bees are linked to peroxisome-induced pathways

The ectoparasitic mite Varroa destructor transmits and triggers viral infections that have deleterious effects on honey bee colonies worldwide. We performed a manipulative experiment in which worker bees collected at emergence were exposed to Varroa for 72 h, and their proteomes were compared with those of untreated control bees. Label-free quantitative proteomics identified 77 differentially expressed A. mellifera proteins (DEPs). In addition, viral proteins were identified by orthogonal analysis, and most importantly, Deformed wing virus (DWV) was found at high levels/intensity in Varroa-exposed bees. Pathway enrichment analysis suggested that the main pathways affected included peroxisomal metabolism, cyto-/exoskeleton reorganization, and cuticular proteins. Detailed examination of individual DEPs revealed that additional changes in DEPs were associated with peroxisomal function. In addition, the proteome data support the importance of TGF-β signaling in Varroa-DWV interaction and the involvement of the mTORC1 and Hippo pathways. These results suggest that the effect of DWV on bees associated with Varroa feeding results in aberrant autophagy. In particular, autophagy is selectively modulated by peroxisomes, to which the observed proteome changes strongly corresponded. This study complements previous research with different study designs and suggests the importance of the peroxisome, which plays a key role in viral infections.

Computational approaches to identify sites of phosphorylation

Due to their oftentimes ambiguous nature, phosphopeptide positional isomers can present challenges in bottom-up mass spectrometry-based workflows as search engine scores alone are often not enough to confidently distinguish them. Additional scoring algorithms can remedy this by providing confidence metrics in addition to these search results, reducing ambiguity. Here we describe challenges to interpreting phosphoproteomics data and review several different approaches to determine sites of phosphorylation for both data-dependent and data-independent acquisition-based workflows. Finally, we discuss open questions regarding neutral losses, gas-phase rearrangement, and false localization rate estimation experienced by both types of acquisition workflows and best practices for managing ambiguity in phosphosite determination.

A critical evaluation of ultrasensitive single-cell proteomics strategies

Success of mass spectrometry characterization of the proteome of single cells allows us to gain a greater understanding than afforded by transcriptomics alone but requires clear understanding of the tradeoffs between analytical throughput and precision. Recent advances in mass spectrometry acquisition techniques, including updated instrumentation and sample preparation, have improved the quality of peptide signals obtained from single cell data. However, much of the proteome remains uncharacterized, and higher throughput techniques often come at the expense of reduced sensitivity and coverage, which diminish the ability to measure proteoform heterogeneity, including splice variants and post-translational modifications, in single cell data analysis. Here, we assess the growing body of ultrasensitive single-cell approaches and their tradeoffs as researchers try to balance throughput and precision in their experiments.

Untargeted metabolomic, and proteomic analysis identifies metabolic biomarkers and pathway alterations in individuals with 22q11.2 deletion syndrome

INTRODUCTION

The chromosome 22q11.2 deletion syndrome (22q11.2DS) is characterized by a well-defined microdeletion and is associated with a wide range of brain-related phenotypes including schizophrenia spectrum disorders (SCZ), autism spectrum disorders (ASD), anxiety disorders and attention deficit disorders (ADHD). The typically deleted region in 22q11.2DS contains multiple genes which haploinsufficiency has the potential of altering the protein and the metabolic profiles.

OBJECTIVES

Alteration in metabolic processes and downstream protein pathways during the early brain development may help to explain the increased prevalence of the observed neurodevelopmental phenotypes in 22q11.2DS. However, relatively little is known about the correlation of dysregulated protein/metabolite expression and neurobehavioral impairments in individuals who developed them over time.

METHODS

In this study, we performed untargeted metabolic and proteomic analysis in plasma samples derived from 30 subjects including 16 participants with 22q11.2DS and 14 healthy controls (TD) enrolled in a longitudinal study, aiming to identify a metabolic and protein signature informing about the underlying mechanisms involved in disease development and progression. The metabolic and proteomic profiles were also compared between the participants with 22q11.2DS with and without various comorbidities, such as medical involvement, psychiatric conditions, and autism spectrum disorder (ASD) to detect potential changes among multiple specimens, collected overtime, with the aim to understand the basic underlying mechanisms involved in disease development and progression.

RESULTS

We observed a large number of statistically significant differences in metabolites between the two groups. Among them, the levels of taurine and arachidonic acid were significantly lower in 22q11.2DS compared to the TD group. In addition, we identified 16 proteins that showed significant changes in expression levels (adjusted P < 0.05) in 22q11.2DS as compared to TD, including those involved in 70 pathways such as gene expression, the PI3K-Akt signaling pathway and the complement system. Within participants with 22q11.2DS, no significant changes in those with and without medical or psychiatric conditions were observed.

CONCLUSION

To our knowledge, this is the first report on plasma metabolic and proteomic profiling and on the identification of unique biomarkers in 22q11.2DS. These findings may suggest the potential role of the identified metabolites and proteins as biomarkers for the onset of comorbid conditions in 22q11.2DS. Ultimately, the altered protein pathways in 22q11.2DS may provide insights of the biological mechanisms underlying the neurodevelopmental phenotype and may provide missing molecular outcome measures in future clinical trials to assess early-diagnosis treatment and the efficacy of response to targeted treatment.

Automated single-cell proteomics providing sufficient proteome depth to study complex biology beyond cell type classifications

Mass spectrometry (MS)-based single-cell proteomics (SCP) has gained massive attention as a viable complement to other single cell approaches. The rapid technological and computational advances in the field have pushed the boundaries of sensitivity and throughput. However, reproducible quantification of thousands of proteins within a single cell at reasonable proteome depth to characterize biological phenomena remains a challenge. To address some of those limitations we present a combination of fully automated single cell sample preparation utilizing a dedicated chip within the picolitre dispensing robot, the cellenONE. The proteoCHIP EVO 96 can be directly interfaced with the Evosep One chromatographic system for in-line desalting and highly reproducible separation with a throughput of 80 samples per day. This, in combination with the Bruker timsTOF MS instruments, demonstrates double the identifications without manual sample handling. Moreover, relative to standard high-performance liquid chromatography, the Evosep One separation provides further 2-fold improvement in protein identifications. The implementation of the newest generation timsTOF Ultra with our proteoCHIP EVO 96-based sample preparation workflow reproducibly identifies up to 4,000 proteins per single HEK-293T without a carrier or match-between runs. Our current SCP depth spans over 4 orders of magnitude and identifies over 50 biologically relevant ubiquitin ligases. We complement our highly reproducible single-cell proteomics workflow to profile hundreds of lipopolysaccharide (LPS)-perturbed THP-1 cells and identified key regulatory proteins involved in interleukin and interferon signaling. This study demonstrates that the proteoCHIP EVO 96-based SCP sample preparation with the timsTOF Ultra provides sufficient proteome depth to study complex biology beyond cell-type classifications.

Label-Free Quantitation of Endogenous Peptides

Liquid chromatography-mass spectrometry (LC-MS)-based peptidomics methods allow for the detection and identification of many peptides in a complex biological mixture in an untargeted manner. Quantitative peptidomics approaches allow for comparisons of peptide abundance between different samples, allowing one to draw conclusions about peptide differences as a function of experimental treatment or physiology. While stable isotope labeling is a powerful approach for quantitative proteomics and peptidomics, advances in mass spectrometry instrumentation and analysis tools have allowed label-free methods to gain popularity in recent years. In a general label-free quantitative peptidomics experiment, peak intensity information for each peptide is compared across multiple LC-MS runs. Here, we outline a general approach for label-free quantitative peptidomics experiments, including steps for sample preparation, LC-MS data acquisition, data processing, and statistical analysis. Special attention is paid to address run-to-run variability, which can lead to several major problems in label-free experiments. Overall, our method provides researchers with a framework for the development of their own quantitative peptidomics workflows applicable to quantitation of peptides from a wide variety of different biological sources.

DeepRTAlign: toward accurate retention time alignment for large cohort mass spectrometry data analysis

Retention time (RT) alignment is a crucial step in liquid chromatography-mass spectrometry (LC-MS)-based proteomic and metabolomic experiments, especially for large cohort studies. The most popular alignment tools are based on warping function method and direct matching method. However, existing tools can hardly handle monotonic and non-monotonic RT shifts simultaneously. Here, we develop a deep learning-based RT alignment tool, DeepRTAlign, for large cohort LC-MS data analysis. DeepRTAlign has been demonstrated to have improved performances by benchmarking it against current state-of-the-art approaches on multiple real-world and simulated proteomic and metabolomic datasets. The results also show that DeepRTAlign can improve identification sensitivity without compromising quantitative accuracy. Furthermore, using the MS features aligned by DeepRTAlign, we trained and validated a robust classifier to predict the early recurrence of hepatocellular carcinoma. DeepRTAlign provides an advanced solution to RT alignment in large cohort LC-MS studies, which is currently a major bottleneck in proteomics and metabolomics research.

Comprehensive Overview of Bottom-Up Proteomics using Mass Spectrometry

Proteomics is the large scale study of protein structure and function from biological systems through protein identification and quantification. "Shotgun proteomics" or "bottom-up proteomics" is the prevailing strategy, in which proteins are hydrolyzed into peptides that are analyzed by mass spectrometry. Proteomics studies can be applied to diverse studies ranging from simple protein identification to studies of proteoforms, protein-protein interactions, protein structural alterations, absolute and relative protein quantification, post-translational modifications, and protein stability. To enable this range of different experiments, there are diverse strategies for proteome analysis. The nuances of how proteomic workflows differ may be challenging to understand for new practitioners. Here, we provide a comprehensive overview of different proteomics methods to aid the novice and experienced researcher. We cover from biochemistry basics and protein extraction to biological interpretation and orthogonal validation. We expect this work to serve as a basic resource for new practitioners in the field of shotgun or bottom-up proteomics.

Flexible Quality Control for Protein Turnover Rates Using d2ome

Bioinformatics tools are used to estimate in vivo protein turnover rates from the LC-MS data of heavy water labeled samples in high throughput. The quantification includes peak detection and integration in the LC-MS domain of complex input data of the mammalian proteome, which requires the integration of results from different experiments. The existing software tools for the estimation of turnover rate use predefined, built-in, stringent filtering criteria to select well-fitted peptides and determine turnover rates for proteins. The flexible control of filtering and quality measures will help to reduce the effects of fluctuations and interferences to the signals from target peptides while retaining an adequate number of peptides. This work describes an approach for flexible error control and filtering measures implemented in the computational tool d2ome for automating protein turnover rates. The error control measures (based on spectral properties and signal features) reduced the standard deviation and tightened the confidence intervals of the estimated turnover rates.

SPPUSM: An MS/MS spectra merging strategy for improved low-input and single-cell proteome identification

Single and rare cell analysis provides unique insights into the investigation of biological processes and disease progress by resolving the cellular heterogeneity that is masked by bulk measurements. Although many efforts have been made, the techniques used to measure the proteome in trace amounts of samples or in single cells still lag behind those for DNA and RNA due to the inherent non-amplifiable nature of proteins and the sensitivity limitation of current mass spectrometry. Here, we report an MS/MS spectra merging strategy termed SPPUSM (same precursor-produced unidentified spectra merging) for improved low-input and single-cell proteome data analysis. In this method, all the unidentified MS/MS spectra from multiple test files are first extracted. Then, the corresponding MS/MS spectra produced by the same precursor ion from different files are matched according to their precursor mass and retention time (RT) and are merged into one new spectrum. The newly merged spectra with more fragment ions are next searched against the database to increase the MS/MS spectra identification and proteome coverage. Further improvement can be achieved by increasing the number of test files and spectra to be merged. Up to 18.2% improvement in protein identification was achieved for 1 ng HeLa peptides by SPPUSM. Reliability evaluation by the "entrapment database" strategy using merged spectra from human and E. coli revealed a marginal error rate for the proposed method. For application in single cell proteome (SCP) study, identification enhancement of 28%-61% was achieved for proteins for different SCP data. Furthermore, a lower abundance was found for the SPPUSM-identified peptides, indicating its potential for more sensitive low sample input and SCP studies.

Proteomics analysis of circulating small extracellular vesicles: Focus on the contribution of EVs to tumor metabolism

The term small extracellular vesicle (sEV) is a comprehensive term that includes any type of cell-derived, membrane-delimited particle that has a diameter < 200 nm, and which includes exosomes and smaller microvesicles. sEVs transfer bioactive molecules between cells and are crucial for cellular homeostasis and particularly during tumor development, where sEVs provide important contributions to the formation of the premetastic niche and to their altered metabolism. sEVs are thus legitimate targets for intervention and have also gained increasing interest as an easily accessible source of biomarkers because they can be rapidly isolated from serum/plasma and their molecular cargo provides information on their cell-of origin. To target sEVs that are specific for a given cell/disease it is essential to identify EV surface proteins that are characteristic of that cell/disease. Mass-spectrometry based proteomics is widely used for the identification and quantification of sEV proteins. The methods used for isolating the sEVs, preparing the sEV sample for proteomics analysis, and mass spectrometry analysis, can have a strong influence on the results and requires careful consideration. This review provides an overview of the approaches used for sEV proteomics and discusses the inherent compromises regarding EV purity versus depth of coverage. Additionally, it discusses the practical applications of the methods to unravel the involvement of sEVs in regulating the metabolism of pancreatic ductal adenocarcinoma (PDAC). The metabolic reprogramming in PDAC includes enhanced glycolysis, elevated glutamine metabolism, alterations in lipid metabolism, mitochondrial dysfunction and hypoxia, all of which are crucial in promoting tumor cell growth. A thorough understanding of these metabolic adaptations is imperative for the development of targeted therapies to exploit PDAC's vulnerabilities.

A Complete Workflow for High Throughput Human Single Skeletal Muscle Fiber Proteomics

Skeletal muscle is a major regulatory tissue of whole-body metabolism and is composed of a diverse mixture of cell (fiber) types. Aging and several diseases differentially affect the various fiber types, and therefore, investigating the changes in the proteome in a fiber-type specific manner is essential. Recent breakthroughs in isolated single muscle fiber proteomics have started to reveal heterogeneity among fibers. However, existing procedures are slow and laborious, requiring 2 h of mass spectrometry time per single muscle fiber; 50 fibers would take approximately 4 days to analyze. Thus, to capture the high variability in fibers both within and between individuals requires advancements in high throughput single muscle fiber proteomics. Here we use a single cell proteomics method to enable quantification of single muscle fiber proteomes in 15 min total instrument time. As proof of concept, we present data from 53 isolated skeletal muscle fibers obtained from two healthy individuals analyzed in 13.25 h. Adapting single cell data analysis techniques to integrate the data, we can reliably separate type 1 and 2A fibers. Ninety-four proteins were statistically different between clusters indicating alteration of proteins involved in fatty acid oxidation, oxidative phosphorylation, and muscle structure and contractile function. Our results indicate that this method is significantly faster than prior single fiber methods in both data collection and sample preparation while maintaining sufficient proteome depth. We anticipate this assay will enable future studies of single muscle fibers across hundreds of individuals, which has not been possible previously due to limitations in throughput.

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

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

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

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

ProInfer: An interpretable protein inference tool leveraging on biological networks

In mass spectrometry (MS)-based proteomics, protein inference from identified peptides (protein fragments) is a critical step. We present ProInfer (Protein Inference), a novel protein assembly method that takes advantage of information in biological networks. ProInfer assists recovery of proteins supported only by ambiguous peptides (a peptide which maps to more than one candidate protein) and enhances the statistical confidence for proteins supported by both unique and ambiguous peptides. Consequently, ProInfer rescues weakly supported proteins thereby improving proteome coverage. Evaluated across THP1 cell line, lung cancer and RAW267.4 datasets, ProInfer always infers the most numbers of true positives, in comparison to mainstream protein inference tools Fido, EPIFANY and PIA. ProInfer is also adept at retrieving differentially expressed proteins, signifying its usefulness for functional analysis and phenotype profiling. Source codes of ProInfer are available at https://github.com/PennHui2016/ProInfer.

High-quality and robust protein quantification in large clinical/pharmaceutical cohorts with IonStar proteomics investigation

Robust, reliable quantification of large sample cohorts is often essential for meaningful clinical or pharmaceutical proteomics investigations, but it is technically challenging. When analyzing very large numbers of samples, isotope labeling approaches may suffer from substantial batch effects, and even with label-free methods, it becomes evident that low-abundance proteins are not reliably measured owing to unsufficient reproducibility for quantification. The MS1-based quantitative proteomics pipeline IonStar was designed to address these challenges. IonStar is a label-free approach that takes advantage of the high sensitivity/selectivity attainable by ultrahigh-resolution (UHR)-MS1 acquisition (e.g., 120-240k full width at half maximum at m/z = 200) which is now widely available on ultrahigh-field Orbitrap instruments. By selectively and accurately procuring quantitative features of peptides within precisely defined, very narrow m/z windows corresponding to the UHR-MS1 resolution, the method minimizes co-eluted interferences and substantially enhances signal-to-noise ratio of low-abundance species by decreasing noise level. This feature results in high sensitivity, selectivity, accuracy and precision for quantification of low-abundance proteins, as well as fewer missing data and fewer false positives. This protocol also emphasizes the importance of well-controlled, robust experimental procedures to achieve high-quality quantification across a large cohort. It includes a surfactant cocktail-aided sample preparation procedure that achieves high/reproducible protein/peptide recoveries among many samples, and a trapping nano-liquid chromatography-mass spectrometry strategy for sensitive and reproducible acquisition of UHR-MS1 peptide signal robustly across a large cohort. Data processing and quality evaluation are illustrated using an example dataset ( http://proteomecentral.proteomexchange.org ), and example results from pharmaceutical project and one clinical project (patients with acute respiratory distress syndrome) are shown. The complete IonStar pipeline takes ~1-2 weeks for a sample cohort containing ~50-100 samples.

Complete Workflow for High Throughput Human Single Skeletal Muscle Fiber Proteomics

Skeletal muscle is a major regulatory tissue of whole-body metabolism and is composed of a diverse mixture of cell (fiber) types. Aging and several diseases differentially affect the various fiber types, and therefore, investigating the changes in the proteome in a fiber-type specific manner is essential. Recent breakthroughs in isolated single muscle fiber proteomics have started to reveal heterogeneity among fibers. However, existing procedures are slow and laborious requiring two hours of mass spectrometry time per single muscle fiber; 50 fibers would take approximately four days to analyze. Thus, to capture the high variability in fibers both within and between individuals requires advancements in high throughput single muscle fiber proteomics. Here we use a single cell proteomics method to enable quantification of single muscle fiber proteomes in 15 minutes total instrument time. As proof of concept, we present data from 53 isolated skeletal muscle fibers obtained from two healthy individuals analyzed in 13.25 hours. Adapting single cell data analysis techniques to integrate the data, we can reliably separate type 1 and 2A fibers. Sixty-five proteins were statistically different between clusters indicating alteration of proteins involved in fatty acid oxidation, muscle structure and regulation. Our results indicate that this method is significantly faster than prior single fiber methods in both data collection and sample preparation while maintaining sufficient proteome depth. We anticipate this assay will enable future studies of single muscle fibers across hundreds of individuals, which has not been possible previously due to limitations in throughput.

Retention Time Alignment for Protein Turnover Studies Using Heavy Water Metabolic Labeling

Retention time (RT) alignment has been important for robust protein identification and quantification in proteomics. In data-dependent acquisition mode, whereby the precursor ions are semistochastically chosen for fragmentation in MS/MS, the alignment is used in an approach termed matched between runs (MBR). MBR transfers peptides, which were fragmented and identified in one experiment, to a replicate experiment where they were not identified. Before the MBR transfer, the RTs of experiments are aligned to reduce the chance of erroneous transfers. Despite its widespread use in other areas of quantitative proteomics, RT alignment has not been applied in data analyses for protein turnover using an atom-based stable isotope-labeling agent such as metabolic labeling with deuterium oxide, D2O. Deuterium incorporation changes isotope profiles of intact peptides in full scans and their fragment ions in tandem mass spectra. It reduces the peptide identification rates in current database search engines. Therefore, the MBR becomes more important. Here, we report on an approach to incorporate RT alignment with peptide quantification in studies of proteome turnover using heavy water metabolic labeling and LC-MS. The RT alignment uses correlation-optimized time warping. The alignment, followed by the MBR, improves labeling time point coverage, especially for long labeling durations.

MS-DAP Platform for Downstream Data Analysis of Label-Free Proteomics Uncovers Optimal Workflows in Benchmark Data Sets and Increased Sensitivity in Analysis of Alzheimer's Biomarker Data

In the rapidly moving proteomics field, a diverse patchwork of data analysis pipelines and algorithms for data normalization and differential expression analysis is used by the community. We generated a mass spectrometry downstream analysis pipeline (MS-DAP) that integrates both popular and recently developed algorithms for normalization and statistical analyses. Additional algorithms can be easily added in the future as plugins. MS-DAP is open-source and facilitates transparent and reproducible proteome science by generating extensive data visualizations and quality reporting, provided as standardized PDF reports. Second, we performed a systematic evaluation of methods for normalization and statistical analysis on a large variety of data sets, including additional data generated in this study, which revealed key differences. Commonly used approaches for differential testing based on moderated t-statistics were consistently outperformed by more recent statistical models, all integrated in MS-DAP. Third, we introduced a novel normalization algorithm that rescues deficiencies observed in commonly used normalization methods. Finally, we used the MS-DAP platform to reanalyze a recently published large-scale proteomics data set of CSF from AD patients. This revealed increased sensitivity, resulting in additional significant target proteins which improved overlap with results reported in related studies and includes a large set of new potential AD biomarkers in addition to previously reported.

Software Tool for Visualization and Validation of Protein Turnover Rates Using Heavy Water Metabolic Labeling and LC-MS

Metabolic stable isotope labeling followed by liquid chromatography coupled with mass spectrometry (LC-MS) is a powerful tool for in vivo protein turnover studies of individual proteins on a large scale and with high throughput. Turnover rates of thousands of proteins from dozens of time course experiments are determined by data processing tools, which are essential components of the workflows for automated extraction of turnover rates. The development of sophisticated algorithms for estimating protein turnover has been emphasized. However, the visualization and annotation of the time series data are no less important. The visualization tools help to validate the quality of the model fits, their goodness-of-fit characteristics, mass spectral features of peptides, and consistency of peptide identifications, among others. Here, we describe a graphical user interface (GUI) to visualize the results from the protein turnover analysis tool, d2ome, which determines protein turnover rates from metabolic D₂O labeling followed by LC-MS. We emphasize the specific features of the time series data and their visualization in the GUI. The time series data visualized by the GUI can be saved in JPEG format for storage and further dissemination.

Conditional Fragment Ion Probabilities Improve Database Searching for Nonmonoisotopic Precursors

Stochastic, intensity-based precursor isolation can result in isotopically enriched fragment ions. This problem is exacerbated for large peptides and stable isotope labeling experiments using deuterium or ¹⁵N. For stable isotope labeling experiments, incomplete and ubiquitous labeling strategies result in the isolation of peptide ions composed of many distinct structural isomers. Unfortunately, existing proteomics search algorithms do not account for this variability in isotopic incorporation, and thus often yield poor peptide and protein identification rates. We sought to resolve this shortcoming by deriving the expected isotopic distributions of each fragment ion and incorporating them into the theoretical mass spectra used for peptide-spectrum-matching. We adapted the Comet search platform to integrate a modified spectral prediction algorithm we term Conditional fragment Ion Distribution Search (CIDS). Comet-CIDS uses a traditional database searching strategy, but for each candidate peptide we compute the isotopic distribution of each fragment to better match the observed m/z distributions. Evaluating previously generated D₂O and ¹⁵N labeled data sets, we found that Comet-CIDS identified more confident peptide spectral matches and higher protein sequence coverage compared to traditional theoretical spectra generation, with the magnitude of improvement largely determined by the amount of labeling in the sample.

Spectral Library-Based Single-Cell Proteomics Resolves Cellular Heterogeneity

Dissecting the proteome of cell types and states at single-cell resolution, while being highly challenging, has significant implications in basic science and biomedicine. Mass spectrometry (MS)-based single-cell proteomics represents an emerging technology for system-wide, unbiased profiling of proteins in single cells. However, significant challenges remain in analyzing an extremely small amount of proteins collected from a single cell, as a proteome-wide amplification of proteins is not currently feasible. Here, we report an integrated spectral library-based single-cell proteomics (SLB-SCP) platform that is ultrasensitive and well suited for a large-scale analysis. To overcome the low MS/MS signal intensity intrinsically associated with a single-cell analysis, this approach takes an alternative approach by extracting a breadth of information that specifically defines the physicochemical characteristics of a peptide from MS1 spectra, including monoisotopic mass, isotopic distribution, and retention time (hydrophobicity), and uses a spectral library for proteomic identification. This conceptually unique MS platform, coupled with the DIRECT sample preparation method, enabled identification of more than 2000 proteins in a single cell to distinguish different proteome landscapes associated with cellular types and heterogeneity. We characterized individual normal and cancerous pancreatic ductal cells (HPDE and PANC-1, respectively) and demonstrated the substantial difference in the proteomes between HPDE and PANC-1 at the single-cell level. A significant upregulation of multiple protein networks in cancer hallmarks was identified in the PANC-1 cells, functionally discriminating the PANC-1 cells from the HPDE cells. This integrated platform can be built on high-resolution MS and widely accepted proteomic software, making it possible for community-wide applications.

Comparative Assessment of Quantification Methods for Tumor Tissue Phosphoproteomics

With increasing sensitivity and accuracy in mass spectrometry, the tumor phosphoproteome is getting into reach. However, the selection of quantitation techniques best-suited to the biomedical question and diagnostic requirements remains a trial and error decision as no study has directly compared their performance for tumor tissue phosphoproteomics. We compared label-free quantification (LFQ), spike-in-SILAC (stable isotope labeling by amino acids in cell culture), and tandem mass tag (TMT) isobaric tandem mass tags technology for quantitative phosphosite profiling in tumor tissue. Compared to the classic SILAC method, spike-in-SILAC is not limited to cell culture analysis, making it suitable for quantitative analysis of tumor tissue samples. TMT offered the lowest accuracy and the highest precision and robustness toward different phosphosite abundances and matrices. Spike-in-SILAC offered the best compromise between these features but suffered from a low phosphosite coverage. LFQ offered the lowest precision but the highest number of identifications. Both spike-in-SILAC and LFQ presented susceptibility to matrix effects. Match between run (MBR)-based analysis enhanced the phosphosite coverage across technical replicates in LFQ and spike-in-SILAC but further reduced the precision and robustness of quantification. The choice of quantitative methodology is critical for both study design such as sample size in sample groups and quantified phosphosites and comparison of published cancer phosphoproteomes. Using ovarian cancer tissue as an example, our study builds a resource for the design and analysis of quantitative phosphoproteomic studies in cancer research and diagnostics.

SIMSI-Transfer: Software-assisted reduction of missing values in phosphoproteomic and proteomic isobaric labeling data using tandem mass spectrum clustering

Isobaric stable isotope labeling techniques such as tandem mass tags (TMTs) have become popular in proteomics because they enable the relative quantification of proteins with high precision from up to 18 samples in a single experiment. While missing values in peptide quantification are rare in a single TMT experiment, they rapidly increase when combining multiple TMT experiments. As the field moves toward analyzing ever higher numbers of samples, tools that reduce missing values also become more important for analyzing TMT datasets. To this end, we developed SIMSI-Transfer (Similarity-based Isobaric Mass Spectra 2 [MS2] Identification Transfer), a software tool that extends our previously developed software MaRaCluster (© Matthew The) by clustering similar tandem MS2 from multiple TMT experiments. SIMSI-Transfer is based on the assumption that similarity-clustered MS2 spectra represent the same peptide. Therefore, peptide identifications made by database searching in one TMT batch can be transferred to another TMT batch in which the same peptide was fragmented but not identified. To assess the validity of this approach, we tested SIMSI-Transfer on masked search engine identification results and recovered >80% of the masked identifications while controlling errors in the transfer procedure to below 1% false discovery rate. Applying SIMSI-Transfer to six published full proteome and phosphoproteome datasets from the Clinical Proteomic Tumor Analysis Consortium led to an increase of 26 to 45% of identified MS2 spectra with TMT quantifications. This significantly decreased the number of missing values across batches and, in turn, increased the number of peptides and proteins identified in all TMT batches by 43 to 56% and 13 to 16%, respectively.

•

Spectrum clustering enables peptide identification transfer between LC–MS/MS runs.

•

The SIMSI pipeline supports processing full proteome and phosphoproteome data.

•

SIMSI increases the number of quantifiable PSMs by 26 to 45%.

•

SIMSI reduces missing values in multibatch TMT labeling experiments by up to 21%.

Predicting missing proteomics values using machine learning: Filling the gap using transcriptomics and other biological features

Proteins are often considered the main biological element in charge of the different functions and structures of a cell. However, proteomics, the global study of all expressed proteins, often performed by mass spectrometry, is limited by its stochastic sampling and can only quantify a limited amount of protein per sample. Transcriptomics, which allows an exhaustive analysis of all expressed transcripts, is often used as a surrogate. However, the transcript level does not present a high level of correlation with the corresponding protein level, notably due to the existence of several post-transcriptional regulatory mechanisms. In this publication, we hypothesize that the missing protein values in proteomics could be predicted using machine learning regression methods, trained with many features extracted from transcriptomics, including known translational regulatory elements such as microRNAs and circular RNAs. After considering different machine learning algorithms applied on two different splitting strategies, we report that random forest can predict proteins in new samples out of transcriptomics data with good accuracy. The proposed pre-processing and model building scripts can be accessed on GitHub: https://github.com/jochotecoa/ml_proteomics.

Quantitative Accuracy and Precision in Multiplexed Single-Cell Proteomics

Single-cell proteomics workflows have considerably improved in sensitivity and reproducibility to characterize as-yet unknown biological phenomena. With the emergence of multiplexed single-cell proteomics, studies increasingly present single-cell measurements in conjunction with an abundant congruent carrier to improve the precursor selection and enhance identifications. While these extreme carrier spikes are often >100× more abundant than the investigated samples, the total ion current undoubtably increases but the quantitative accuracy possibly is affected. We here focus on narrowly titrated carrier spikes (i.e., <20×) and assess their elimination for a comparable sensitivity with superior accuracy. We find that subtle changes in the carrier ratio can severely impact the measurement variability and describe alternative multiplexing strategies to evaluate data quality. Lastly, we demonstrate elevated replicate overlap while preserving acquisition throughput at an improved quantitative accuracy with DIA-TMT and discuss optimized experimental designs for multiplexed proteomics of trace samples. This comprehensive benchmarking gives an overview of currently available techniques and guides the conceptualization of the optimal single-cell proteomics experiment.

Domont G, C. S. Nogueira F, Marko-Varga G, Sanchez A. Mapping the Melanoma Plasma Proteome (MPP) Using Single-Shot Proteomics Interfaced with the WiMT Database

Plasma analysis by mass spectrometry-based proteomics remains a challenge due to its large dynamic range of 10 orders in magnitude. We created a methodology for protein identification known as Wise MS Transfer (WiMT). Melanoma plasma samples from biobank archives were directly analyzed using simple sample preparation. WiMT is based on MS1 features between several MS runs together with custom protein databases for ID generation. This entails a multi-level dynamic protein database with different immunodepletion strategies by applying single-shot proteomics. The highest number of melanoma plasma proteins from undepleted and unfractionated plasma was reported, mapping >1200 proteins from >10,000 protein sequences with confirmed significance scoring. Of these, more than 660 proteins were annotated by WiMT from the resulting ~5800 protein sequences. We could verify 4000 proteins by MS1t analysis from HeLA extracts. The WiMT platform provided an output in which 12 previously well-known candidate markers were identified. We also identified low-abundant proteins with functions related to (i) cell signaling, (ii) immune system regulators, and (iii) proteins regulating folding, sorting, and degradation, as well as (iv) vesicular transport proteins. WiMT holds the potential for use in large-scale screening studies with simple sample preparation, and can lead to the discovery of novel proteins with key melanoma disease functions.

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

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

Evaluation of the Sensitivity of Proteomics Methods Using the Absolute Copy Number of Proteins in a Single Cell as a Metric

Proteomic technology has improved at a staggering pace in recent years, with even practitioners challenged to keep up with new methods and hardware. The most common metric used for method performance is the number of peptides and proteins identified. While this metric may be helpful for proteomics researchers shopping for new hardware, this is often not the most biologically relevant metric. Biologists often utilize proteomics in the search for protein regulators that are of a lower relative copy number in the cell. In this review, I re-evaluate untargeted proteomics data using a simple graphical representation of the absolute copy number of proteins present in a single cancer cell as a metric. By comparing single-shot proteomics data to the coverage of the most in-depth proteomic analysis of that cell line acquired to date, we can obtain a rapid metric of method performance. Using a simple copy number metric allows visualization of how proteomics has developed in both sensitivity and overall dynamic range when using both relatively long and short acquisition times. To enable reanalysis beyond what is presented here, two available web applications have been developed for single- and multi-experiment comparisons with reference protein copy number data for multiple cell lines and organisms.

Comprehensive proteomic analysis of murine terminal erythroid differentiation

Murine-based cellular models have provided and continue to provide many useful insights into the fundamental mechanisms of erythropoiesis, as well as insights into the pathophysiology of inherited and acquired red cell disorders. Although detailed information on many aspects of these cell models is available, comprehensive proteomic data are lacking. This is a critical knowledge gap, as proteins are effectors of most biologic processes. To address this critical unmet need, proteomes of the murine cell lines Friend erythroleukemia (MEL), GATA1 erythroid (G1ER), and embryonic stem cell-derived erythroid progenitor (MEDEP) and proteomes of cultured murine marrow-derived erythroblasts at different stages of terminal erythroid differentiation were analyzed. The proteomes of MEDEP cells and primary murine erythroid cells were most similar, whereas those of MEL and G1ER cells were more distantly related. We demonstrated that the overall cellular content of histones does not decrease during terminal differentiation, despite strong chromatin condensation. Comparison of murine and human proteomes throughout terminal erythroid differentiation revealed that many noted transcriptomic changes were significantly dampened at the proteome level, especially at the end of the terminal differentiation process. Analysis of the early events associated with induction of terminal differentiation in MEDEP cells revealed divergent alterations in associated transcriptomes and proteomes. These proteomic data are powerful and valuable tools for the study of fundamental mechanisms of normal and disordered erythropoiesis and will be of broad interest to a wide range of investigators for making the appropriate choice of various cell lines to study inherited and acquired diseases of the erythrocyte.

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

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

Pilot proteomic analysis of cerebrospinal fluid in Alzheimer's disease

Proteomic analysis of cerebrospinal fluid (CSF) holds great promise in understanding the progression of neurodegenerative diseases, including Alzheimer's disease (AD). As one of the primary reservoirs of neuronal biomolecules, CSF provides a window into the biochemical and cellular aspects of the neurological environment. CSF can be drawn from living participants allowing the potential alignment of clinical changes with these biochemical markers. Using cutting-edge mass spectrometry technologies, we perform a streamlined proteomic analysis of CSF. We quantify greater than 700 proteins across 10 pairs of age- and sex-matched participants in approximately one hour of analysis time each. Using the paired participant study structure, we identify a small group of biologically relevant proteins that show substantial changes in abundance between cognitive normal and AD participants, which were then analyzed at the peptide level using parallel reaction monitoring experiments. Our findings suggest the utility of fractionating a single sample and using matching to increase proteomic depth in cerebrospinal fluid, as well as the potential power of an expanded study.

IonQuant Enables Accurate and Sensitive Label-Free Quantification With FDR-Controlled Match-Between-Runs

Missing values weaken the power of label-free quantitative proteomic experiments to uncover true quantitative differences between biological samples or experimental conditions. Match-between-runs (MBR) has become a common approach to mitigate the missing value problem, where peptides identified by tandem mass spectra in one run are transferred to another by inference based on m/z, charge state, retention time, and ion mobility when applicable. Though tolerances are used to ensure such transferred identifications are reasonably located and meet certain quality thresholds, little work has been done to evaluate the statistical confidence of MBR. Here, we present a mixture model-based approach to estimate the false discovery rate (FDR) of peptide and protein identification transfer, which we implement in the label-free quantification tool IonQuant. Using several benchmarking datasets generated on both Orbitrap and timsTOF mass spectrometers, we demonstrate superior performance of IonQuant with FDR-controlled MBR compared with MaxQuant (19-38 times faster; 6-18% more proteins quantified and with comparable or better accuracy). We further illustrate the performance of IonQuant and highlight the need for FDR-controlled MBR, in two single-cell proteomics experiments, including one acquired with the help of high-field asymmetric ion mobility spectrometry separation. Fully integrated in the FragPipe computational environment, IonQuant with FDR-controlled MBR enables fast and accurate peptide and protein quantification in label-free proteomics experiments.

Feasibility of Phosphoproteomics on Leftover Samples After RNA Extraction With Guanidinium Thiocyanate

In daily practice, different types of biomolecules are usually extracted for large-scale "omics" analysis with tailored protocols. However, when sample material is limited, an all-in-one strategy is preferable. Although lysis of cells and tissues with urea is widely used for phosphoproteomic applications, DNA, RNA, and proteins can be simultaneously extracted from small samples using acid guanidinium thiocyanate-phenol-chloroform (AGPC). Use of AGPC for mass spectrometry-based phosphoproteomics was reported but has not yet been thoroughly evaluated against a classical phosphoproteomic protocol. Here we compared urea- with AGPC-based protein extraction, profiling phosphorylations in the DNA damage response pathway after ionizing irradiation of U2OS cells as proof of principle. On average we identified circa 9000 phosphosites per sample with both extraction methods. Moreover, we observed high similarity of phosphosite characteristics (e.g., 94% shared class 1 identifications) and deduced kinase activities (e.g., ATM, ATR, CHEK1/2, PRKDC). We furthermore extended our comparison to murine and human tissue samples yielding similar and highly correlated results for both extraction protocols. AGPC-based sample extraction can thus replace common cell lysates for phosphoproteomic workflows and may thus be an attractive way to obtain input material for multiple omics workflows, yielding several data types from a single sample.

The rise of single‐cell proteomics

Mass spectrometry‐based proteomics comprehensively defines proteome expression patterns in thousands of cells majorly contributing to our current understanding of many biological processes. More recently, single‐cell transcriptome and genome studies, however, have demonstrated overwhelming heterogeneity of tissues and cellular subpopulations. These studies have indicated different cellular functionality and identity, which are mainly driven by proteins and their posttranscriptional modifications. The rapidly emerging field of single‐cell proteomics aims at complementing transcriptome and genome data by generating comparative protein expression profiles from individual cells. Recent developments demonstrated tremendous improvements in sample preparation workflows and MS instrumentation, quantifying over 1000 proteins from a single cell. Efficient and reproducible sample processing in conjunction with sensitive MS acquisition strategies will allow to further increase the proteome coverage of tissues with single‐cell resolution. The required throughput and data reliability of such studies are still subject to further developments. Therefore, we herein discuss recent progress on specialized workflows and instrumentation next to advancements outside the field, which we expect to contribute to the development of comprehensive single‐cell proteomics.

Triqler for MaxQuant: Enhancing Results from MaxQuant by Bayesian Error Propagation and Integration

Error estimation for differential protein quantification by label-free shotgun proteomics is challenging due to the multitude of error sources, each contributing uncertainty to the final results. We have previously designed a Bayesian model, Triqler, to combine such error terms into one combined quantification error. Here we present an interface for Triqler that takes MaxQuant results as input, allowing quick reanalysis of already processed data. We demonstrate that Triqler outperforms the original processing for a large set of both engineered and clinical/biological relevant data sets. Triqler and its interface to MaxQuant are available as a Python module under an Apache 2.0 license from https://pypi.org/project/triqler/.

Global and Site-Specific Effect of Phosphorylation on Protein Turnover

To date, the effects of specific modification types and sites on protein lifetime have not been systematically illustrated. Here, we describe a proteomic method, DeltaSILAC, to quantitatively assess the impact of site-specific phosphorylation on the turnover of thousands of proteins in live cells. Based on the accurate and reproducible mass spectrometry-based method, a pulse labeling approach using stable isotope-labeled amino acids in cells (pSILAC), phosphoproteomics, and a unique peptide-level matching strategy, our DeltaSILAC profiling revealed a global, unexpected delaying effect of many phosphosites on protein turnover. We further found that phosphorylated sites accelerating protein turnover are functionally selected for cell fitness, enriched in Cyclin-dependent kinase substrates, and evolutionarily conserved, whereas the glutamic acids surrounding phosphosites significantly delay protein turnover. Our method represents a generalizable approach and provides a rich resource for prioritizing the effects of phosphorylation sites on protein lifetime in the context of cell signaling and disease biology.

Integrated analysis of the aging brain transcriptome and proteome in tauopathy

Background

Tau neurofibrillary tangle pathology characterizes Alzheimer’s disease and other neurodegenerative tauopathies. Brain gene expression profiles can reveal mechanisms; however, few studies have systematically examined both the transcriptome and proteome or differentiated Tau- versus age-dependent changes.

Methods

Paired, longitudinal RNA-sequencing and mass-spectrometry were performed in a Drosophila model of tauopathy, based on pan-neuronal expression of human wildtype Tau (Tau^WT) or a mutant form causing frontotemporal dementia (Tau^R406W). Tau-induced, differentially expressed transcripts and proteins were examined cross-sectionally or using linear regression and adjusting for age. Hierarchical clustering was performed to highlight network perturbations, and we examined overlaps with human brain gene expression profiles in tauopathy.

Results

Tau^WT induced 1514 and 213 differentially expressed transcripts and proteins, respectively. Tau^R406W had a substantially greater impact, causing changes in 5494 transcripts and 697 proteins. There was a ~ 70% overlap between age- and Tau-induced changes and our analyses reveal pervasive bi-directional interactions. Strikingly, 42% of Tau-induced transcripts were discordant in the proteome, showing opposite direction of change. Tau-responsive gene expression networks strongly implicate innate immune activation. Cross-species analyses pinpoint human brain gene perturbations specifically triggered by Tau pathology and/or aging, and further differentiate between disease amplifying and protective changes.

Conclusions

Our results comprise a powerful, cross-species functional genomics resource for tauopathy, revealing Tau-mediated disruption of gene expression, including dynamic, age-dependent interactions between the brain transcriptome and proteome.