Real-Time Search-Assisted Acquisition on a Tribrid Mass Spectrometer Improves Coverage in Multiplexed Single-Cell Proteomics

In the young field of single-cell proteomics (scMS), there is a great need for improved global proteome characterization, both in terms of proteins quantified per cell and quantitative performance thereof. The recently introduced real-time search (RTS) on the Orbitrap Eclipse Tribrid mass spectrometer in combination with SPS-MS3 acquisition has been shown to be beneficial for the measurement of samples that are multiplexed using isobaric tags. Multiplexed scMS requires high ion injection times and high-resolution spectra to quantify the single-cell signal; however, the carrier channel facilitates peptide identification and thus offers the opportunity for fast on-the-fly precursor filtering before committing to the time-intensive quantification scan. Here, we compared classical MS2 acquisition against RTS-SPS-MS3, both using the Orbitrap Eclipse Tribrid MS with the FAIMS Pro ion mobility interface and present a new acquisition strategy termed RETICLE (RTS enhanced quant of single cell spectra) that makes use of fast real-time searched linear ion trap scans to preselect MS1 peptide precursors for quantitative MS2 Orbitrap acquisition. We show that classical MS2 acquisition is outperformed by both RTS-SPS-MS3 through increased quantitative accuracy at similar proteome coverage, and RETICLE through higher proteome coverage, with the latter enabling the quantification of over 1000 proteins per cell at an MS2 injection time of 750 ms using a 2 h gradient.     

A Novel Spectral Annotation Strategy Streamlines Reporting of Mono-ADP-ribosylated Peptides Derived from Mouse Liver and Spleen in Response to IFN-γ

Mass-spectrometry-enabled ADP-ribosylation workflows are developing rapidly, providing researchers a variety of ADP-ribosylome enrichment strategies and mass spectrometric acquisition options. Despite the growth spurt in upstream technologies, systematic ADP-ribosyl (ADPr) peptide mass spectral annotation methods are lacking. HCD-dependent ADP-ribosylome studies are common, but the resulting MS2 spectra are complex, owing to a mixture of b/y-ions and the m/p-ion peaks representing one or more dissociation events of the ADPr moiety (m-ion) and peptide (p-ion). In particular, p-ions that dissociate further into one or more fragment ions can dominate HCD spectra but are not recognized by standard spectral annotation workflows. As a result, annotation strategies that are solely reliant upon the b/y-ions result in lower spectral scores that in turn reduce the number of reportable ADPr peptides. To improve the confidence of spectral assignments, we implemented an ADPr peptide annotation and scoring strategy. All MS2 spectra are scored for the ADPr m-ions, but once spectra are assigned as an ADPr peptide, they are further annotated and scored for the p-ions. We implemented this novel workflow to ADPr peptides enriched from the liver and spleen isolated from mice post 4 h exposure to systemic IFN-γ. HCD collision energy experiments were first performed on the Orbitrap Fusion Lumos and the Q Exactive, with notable ADPr peptide dissociation properties verified with CID (Lumos). The m-ion and p-ion series score distributions revealed that ADPr peptide dissociation properties vary markedly between instruments and within instrument collision energy settings, with consequences on ADPr peptide reporting and amino acid localization. Consequentially, we increased the number of reportable ADPr peptides by 25% (liver) and 17% (spleen) by validation and the inclusion of lower confidence ADPr peptide spectra. This systematic annotation strategy will streamline future reporting of ADPr peptides that have been sequenced using any HCD/CID-based method.     

Improving Statistical Certainty of Glycosylation Similarity between Influenza A Virus Variants Using Data-Independent Acquisition Mass Spectrometry

Amino acid sequences of immunodominant domains of hemagglutinin (HA) on the surface of influenza A virus (IAV) evolve rapidly, producing viral variants. HA mediates receptor recognition, binding and cell entry, and serves as the target for IAV vaccines. Glycosylation, a post-translational modification that places large branched polysaccharide molecules on proteins, can modulate the function of HA and shield antigenic regions allowing for viral evasion from immune responses. Our previous work showed that subtle changes in the HA protein sequence can have a measurable change in glycosylation. Thus, being able to quantitatively measure glycosylation changes in variants is critical for understanding how HA function may change throughout viral evolution. Moreover, understanding quantitatively how the choice of viral expression systems affects glycosylation can help in the process of vaccine design and manufacture. Although IAV vaccines are most commonly expressed in chicken eggs, cell-based vaccines have many advantages, and the adoption of more cell-based vaccines would be an important step in mitigating seasonal influenza and protecting against future pandemics. Here, we have investigated the use of data-independent acquisition (DIA) mass spectrometry for quantitative glycoproteomics. We found that DIA improved the sensitivity of glycopeptide detection for four variants of A/Switzerland/9715293/2013 (H3N2): WT and mutant, each expressed in embryonated chicken eggs and Madin–Darby canine kidney cells. We used the Tanimoto similarity metric to quantify changes in glycosylation between WT and mutant and between egg-expressed and cell-expressed virus. Our DIA site-specific glycosylation similarity comparison of WT and mutant expressed in eggs confirmed our previous analysis while achieving greater depth of coverage. We found that sequence variations and changing viral expression systems affected distinct glycosylation sites of HA. Our methods can be applied to track glycosylation changes in circulating IAV variants to bolster genomic surveillance already being done, for a more complete understanding of IAV evolution.     

Multiattribute Glycan Identification and FDR Control for Glycoproteomics

Rapidly improving methods for glycoproteomics have enabled increasingly large-scale analyses of complex glycopeptide samples, but annotating the resulting mass spectrometry data with high confidence remains a major bottleneck. We recently introduced a fast and sensitive glycoproteomics search method in our MSFragger search engine, which reports glycopeptides as a combination of a peptide sequence and the mass of the attached glycan. In samples with complex glycosylation patterns, converting this mass to a specific glycan composition is not straightforward; however, as many glycans have similar or identical masses. Here, we have developed a new method for determining the glycan composition of N-linked glycopeptides fragmented by collisional or hybrid activation that uses multiple sources of information from the spectrum, including observed glycan B-type (oxonium) and Y-type ions and mass and precursor monoisotopic selection errors to discriminate between possible glycan candidates. Combined with false discovery rate estimation for the glycan assignment, we show that this method is capable of specifically and sensitively identifying glycans in complex glycopeptide analyses and effectively controls the rate of false glycan assignments. The new method has been incorporated into the PTM-Shepherd modification analysis tool to work directly with the MSFragger glyco search in the FragPipe graphical user interface, providing a complete computational pipeline for annotation of N-glycopeptide spectra with false discovery rate control of both peptide and glycan components that is both sensitive and robust against false identifications.     

Recent advancements in mass spectrometry–based tools to investigate newly synthesized proteins

Tight regulation of protein translation drives the proteome to undergo changes under influence of extracellular or intracellular signals. Despite mass spectrometry–based proteomics being an excellent method to study differences in protein abundance in complex proteomes, analyzing minute or rapid changes in protein synthesis and abundance remains challenging. Therefore, several dedicated techniques to directly detect and quantify newly synthesized proteins have been developed, notably puromycin-based, bio-orthogonal noncanonical amino acid tagging–based, and stable isotope labeling by amino acids in cell culture–based methods, combined with mass spectrometry. These techniques have enabled the investigation of perturbations, stress, or stimuli on protein synthesis. Improvements of these methods are still necessary to overcome various remaining limitations. Recent improvements include enhanced enrichment approaches and combinations with various stable isotope labeling techniques, which allow for more accurate analysis and comparison between conditions on shorter timeframes and in more challenging systems. Here, we aim to review the current state in this field.     

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.     

Open Modification Searching of SARS-CoV-2–Human Protein Interaction Data Reveals Novel Viral Modification Sites

The outbreak of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of the coronavirus 2019 disease, has led to an ongoing global pandemic since 2019. Mass spectrometry can be used to understand the molecular mechanisms of viral infection by SARS-CoV-2, for example, by determining virus–host protein–protein interactions through which SARS-CoV-2 hijacks its human hosts during infection, and to study the role of post-translational modifications. We have reanalyzed public affinity purification–mass spectrometry data using open modification searching to investigate the presence of post-translational modifications in the context of the SARS-CoV-2 virus–host protein–protein interaction network. Based on an over twofold increase in identified spectra, our detected protein interactions show a high overlap with independent mass spectrometry-based SARS-CoV-2 studies and virus–host interactions for alternative viruses, as well as previously unknown protein interactions. In addition, we identified several novel modification sites on SARS-CoV-2 proteins that we investigated in relation to their interactions with host proteins. A detailed analysis of relevant modifications, including phosphorylation, ubiquitination, and S-nitrosylation, provides important hypotheses about the functional role of these modifications during viral infection by SARS-CoV-2.     

DIA-Based Proteomics Identifies IDH2 as a Targetable Regulator of Acquired Drug Resistance in Chronic Myeloid Leukemia

Drug resistance is a critical obstacle to effective treatment in patients with chronic myeloid leukemia. To understand the underlying resistance mechanisms in response to imatinib mesylate (IMA) and adriamycin (ADR), the parental K562 cells were treated with low doses of IMA or ADR for 2 months to generate derivative cells with mild, intermediate, and severe resistance to the drugs as defined by their increasing resistance index. PulseDIA-based (DIA [data-independent acquisition]) quantitative proteomics was then employed to reveal the proteome changes in these resistant cells. In total, 7082 proteins from 98,232 peptides were identified and quantified from the dataset using four DIA software tools including OpenSWATH, Spectronaut, DIA-NN, and EncyclopeDIA. Sirtuin signaling pathway was found to be significantly enriched in both ADR-resistant and IMA-resistant K562 cells. In particular, isocitrate dehydrogenase (NADP(+)) 2 was identified as a potential drug target correlated with the drug resistance phenotype, and its inhibition by the antagonist AGI-6780 reversed the acquired resistance in K562 cells to either ADR or IMA. Together, our study has implicated isocitrate dehydrogenase (NADP(+)) 2 as a potential target that can be therapeutically leveraged to alleviate the drug resistance in K562 cells when treated with IMA and ADR.     

ImShot: An Open-Source Software for Probabilistic Identification of Proteins In Situ and Visualization of Proteomics Data

Imaging mass spectrometry (IMS) has developed into a powerful tool allowing label-free detection of numerous biomolecules in situ. In contrast to shotgun proteomics, proteins/peptides can be detected directly from biological tissues and correlated to its morphology leading to a gain of crucial clinical information. However, direct identification of the detected molecules is currently challenging for MALDI–IMS, thereby compelling researchers to use complementary techniques and resource intensive experimental setups. Despite these strategies, sufficient information could not be extracted because of lack of an optimum data combination strategy/software. Here, we introduce a new open-source software ImShot that aims at identifying peptides obtained in MALDI–IMS. This is achieved by combining information from IMS and shotgun proteomics (LC–MS) measurements of serial sections of the same tissue. The software takes advantage of a two-group comparison to determine the search space of IMS masses after deisotoping the corresponding spectra. Ambiguity in annotations of IMS peptides is eliminated by introduction of a novel scoring system that identifies the most likely parent protein of a detected peptide in the corresponding IMS dataset. Thanks to its modular structure, the software can also handle LC–MS data separately and display interactive enrichment plots and enriched Gene Ontology terms or cellular pathways. The software has been built as a desktop application with a conveniently designed graphic user interface to provide users with a seamless experience in data analysis. ImShot can run on all the three major desktop operating systems and is freely available under Massachusetts Institute of Technology license.     

Sensitive UHPLC-MS/MS quantification method for 4β- and 4α-hydroxycholesterol in plasma for accurate CYP3A phenotyping

4β-Hydroxycholesterol (4β-OHC) is formed by Cytochrome P450 (CYP)3A and has drawn attention as an endogenous phenotyping probe for CYP3A activity. However, 4β-OHC is also increased by cholesterol autooxidation occurring in vitro due to dysregulated storage and in vivo by oxidative stress or inflammation, independent of CYP3A activity. 4α-hydroxycholesterol (4α-OHC), a stereoisomer of 4β-OHC, is also formed via autooxidation of cholesterol, not by CYP3A, and thus may have clinical potential in reflecting the state of cholesterol autooxidation. In this study, we establish a sensitive method for simultaneous quantification of 4β-OHC and 4α-OHC in human plasma using ultra-high performance liquid chromatography coupled to tandem mass spectrometry. Plasma samples were prepared by saponification, two-step liquid-liquid extraction, and derivatization using picolinic acid. Intense [M+H]+ signals for 4β-OHC and 4α-OHC di-picolinyl esters were monitored using electrospray ionization. The assay fulfilled the requirements of the US Food and Drug Administration guidance for bioanalytical method validation, with a lower limit of quantification of 0.5 ng/ml for both 4β-OHC and 4α-OHC. Apparent recovery rates from human plasma ranged from 88.2% to 101.5% for 4β-OHC, and 91.8% to 114.9% for 4α-OHC. Additionally, matrix effects varied between 86.2% and 117.6% for 4β-OHC and between 89.5% and 116.9% for 4α-OHC. Plasma 4β-OHC and 4α-OHC concentrations in healthy volunteers, stage 3–5 chronic kidney disease (CKD) patients, and stage 5D CKD patients as measured by the validated assay were within the calibration ranges in all samples. We propose this novel quantification method may contribute to accurate evaluation of in vivo CYP3A activity.     

Immunopeptidomic Analyses of Colorectal Cancers With and Without Microsatellite Instability

Colorectal cancer is the second leading cause of cancer death worldwide, and the incidence of this disease is expected to increase as global socioeconomic changes occur. Immune checkpoint inhibition therapy is effective in treating a minority of colorectal cancer tumors; however, microsatellite stable tumors do not respond well to this treatment. Emerging cancer immunotherapeutic strategies aim to activate a cytotoxic T cell response against tumor-specific antigens, presented exclusively at the cell surface of cancer cells. These antigens are rare and are most effectively identified with a mass spectrometry–based approach, which allows the direct sampling and sequencing of these peptides. Although the few tumor-specific antigens identified to date are derived from coding regions of the genome, recent findings indicate that a large proportion of tumor-specific antigens originate from allegedly noncoding regions. Here, we employed a novel proteogenomic approach to identify tumor antigens in a collection of colorectal cancer–derived cell lines and biopsy samples consisting of matched tumor and normal adjacent tissue. The generation of personalized cancer databases paired with mass spectrometry analyses permitted the identification of more than 30,000 unique MHC I–associated peptides. We identified 19 tumor-specific antigens in both microsatellite stable and unstable tumors, over two-thirds of which were derived from noncoding regions. Many of these peptides were derived from source genes known to be involved in colorectal cancer progression, suggesting that antigens from these genes could have therapeutic potential in a wide range of tumors. These findings could benefit the development of T cell–based vaccines, in which T cells are primed against these antigens to target and eradicate tumors. Such a vaccine could be used in tandem with existing immune checkpoint inhibition therapies, to bridge the gap in treatment efficacy across subtypes of colorectal cancer with varying prognoses. Data are available via ProteomeXchange with identifier PXD028309.     

Proteomic Discovery of Plasma Protein Biomarkers and Development of Models Predicting Prognosis of High-Grade Serous Ovarian Carcinoma

Ovarian cancer is one of the most lethal female cancers. For accurate prognosis prediction, this study aimed to investigate novel, blood-based prognostic biomarkers for high-grade serous ovarian carcinoma (HGSOC) using mass spectrometry–based proteomics methods. We conducted label-free liquid chromatography–tandem mass spectrometry using frozen plasma samples obtained from patients with newly diagnosed HGSOC (n = 20). Based on progression-free survival (PFS), the samples were divided into two groups: good (PFS ≥18 months) and poor prognosis groups (PFS <18 months). Proteomic profiles were compared between the two groups. Referring to proteomics data that we previously obtained using frozen cancer tissues from chemotherapy-naïve patients with HGSOC, overlapping protein biomarkers were selected as candidate biomarkers. Biomarkers were validated using an independent set of HGSOC plasma samples (n = 202) via enzyme-linked immunosorbent assay (ELISA). To construct models predicting the 18-month PFS rate, we performed stepwise selection based on the area under the receiver operating characteristic curve (AUC) with 5-fold cross-validation. Analysis of differentially expressed proteins in plasma samples revealed that 35 and 61 proteins were upregulated in the good and poor prognosis groups, respectively. Through hierarchical clustering and bioinformatic analyses, GSN, VCAN, SND1, SIGLEC14, CD163, and PRMT1 were selected as candidate biomarkers and were subjected to ELISA. In multivariate analysis, plasma GSN was identified as an independent poor prognostic biomarker for PFS (adjusted hazard ratio, 1.556; 95% confidence interval, 1.073–2.256; p = 0.020). By combining clinical factors and ELISA results, we constructed several models to predict the 18-month PFS rate. A model consisting of four predictors (FIGO stage, residual tumor after surgery, and plasma levels of GSN and VCAN) showed the best predictive performance (mean validated AUC, 0.779). The newly developed model was converted to a nomogram for clinical use. Our study results provided insights into protein biomarkers, which might offer clues for developing therapeutic targets.     

Proteome Remodeling of the Eye Lens at 50 Years Identified With Data-Independent Acquisition

The eye lens is responsible for focusing and transmitting light to the retina. The lens does this in the absence of organelles, yet maintains transparency for at least 5 decades before onset of age-related nuclear cataract (ARNC). It is hypothesized that oxidative stress contributes significantly to ARNC formation. It is in addition hypothesized that transparency is maintained by a microcirculation system that delivers antioxidants to the lens nucleus and exports small molecule waste. Common data-dependent acquisition methods are hindered by dynamic range of lens protein expression and provide limited context to age-related changes in the lens. In this study, we utilized data-independent acquisition mass spectrometry to analyze the urea-insoluble membrane protein fractions of 16 human lenses subdivided into three spatially distinct lens regions to characterize age-related changes, particularly concerning the lens microcirculation system and oxidative stress response. In this pilot cohort, we measured 4788 distinct protein groups, 46,681 peptides, and 7592 deamidated sequences, more than in any previous human lens data-dependent acquisition approach. Principally, we demonstrate that a significant proteome remodeling event occurs at approximately 50 years of age, resulting in metabolic preference for anaerobic glycolysis established with organelle degradation, decreased abundance of protein networks involved in calcium-dependent cell–cell contacts while retaining networks related to oxidative stress response. Furthermore, we identified multiple antioxidant transporter proteins not previously detected in the human lens and describe their spatiotemporal and age-related abundance changes. Finally, we demonstrate that aquaporin-5, among other proteins, is modified with age by post-translational modifications including deamidation and truncation. We suggest that the continued accumulation of each of these age-related outcomes in proteome remodeling contribute to decreased fiber cell permeability and result in ARNC formation.     

基于数据非依赖采集的以肽为中心分析算法和软件的研究进展

数据非依赖采集(data-independent acquisition,DIA)是一种高通量、无偏性的质谱数据采集方法,具有定量结果重现性好,对低丰度蛋白质友好的特点,是近年来进行大队列蛋白质组研究的首选方法之一。由于DIA产生的二级谱是混合谱,包含了多个肽段的碎片离子信息,使得蛋白质鉴定和定量更加困难。目前,DIA数据分析方法分为两大类,即以肽为中心和以谱图为中心。其中,以肽为中心的分析方法鉴定更灵敏,定量更准确,已成为DIA数据解析的主流方法。其分析流程包括构建谱图库、提取色谱峰群、特征打分和结果质控4个关键步骤。本文综述了以肽为中心的DIA数据分析流程,介绍了基于此流程的数据分析软件及相关比较评估工作,进一步总结了已有的算法改进工作,最后对未来发展方向进行了展望。     

Spatial single cell metabolomics: Current challenges and future developments

Single cell metabolomics is a rapidly advancing field of bio-analytical chemistry which aims to observe cellular biology with the greatest detail possible. Mass spectrometry imaging and selective cell sampling (e.g. using nanocapillaries) are two common approaches within the field. Recent achievements such as observation of cell–cell interactions, lipids determining cell states and rapid phenotypic identification demonstrate the efficacy of these approaches and the momentum of the field. However, single cell metabolomics can only continue with the same impetus if the universal challenges to the field are met, such as the lack of strategies for standardisation and quantification, and lack of specificity/sensitivity. Mass spectrometry imaging and selective cell sampling come with unique advantages and challenges which, in many cases are complementary to each other. We propose here that the challenges specific to each approach could be ameliorated with collaboration between the two communities driving these approaches.     

Early Prediction of COVID-19 Patient Survival by Targeted Plasma Multi-Omics and Machine Learning

The recent surge of coronavirus disease 2019 (COVID-19) hospitalizations severely challenges healthcare systems around the globe and has increased the demand for reliable tests predictive of disease severity and mortality. Using multiplexed targeted mass spectrometry assays on a robust triple quadrupole MS setup which is available in many clinical laboratories, we determined the precise concentrations of hundreds of proteins and metabolites in plasma from hospitalized COVID-19 patients. We observed a clear distinction between COVID-19 patients and controls and, strikingly, a significant difference between survivors and nonsurvivors. With increasing length of hospitalization, the survivors’ samples showed a trend toward normal concentrations, indicating a potential sensitive readout of treatment success. Building a machine learning multi-omic model that considers the concentrations of 10 proteins and five metabolites, we could predict patient survival with 92% accuracy (area under the receiver operating characteristic curve: 0.97) on the day of hospitalization. Hence, our standardized assays represent a unique opportunity for the early stratification of hospitalized COVID-19 patients.     

Proteomic Profiling of Intra-Islet Features Reveals Substructure-Specific Protein Signatures

Despite their diminutive size, islets of Langerhans play a large role in maintaining systemic energy balance in the body. New technologies have enabled us to go from studying the whole pancreas to isolated whole islets, to partial islet sections, and now to islet substructures isolated from within the islet. Using a microfluidic nanodroplet-based proteomics platform coupled with laser capture microdissection and field asymmetric waveform ion mobility spectrometry, we present an in-depth investigation of protein profiles specific to features within the islet. These features include the islet-acinar interface vascular tissue, inner islet vasculature, isolated endocrine cells, whole islet with vasculature, and acinar tissue from around the islet. Compared to interface vasculature, unique protein signatures observed in the inner vasculature indicate increased innervation and intra-islet neuron-like crosstalk. We also demonstrate the utility of these data for identifying localized structure-specific drug–target interactions using existing protein/drug binding databases.     

Quantitative interactome analysis with chemical cross-linking and mass spectrometry

Structural plasticity and dynamic protein–protein interactions are critical determinants of protein function within living systems. Quantitative chemical cross-linking with mass spectrometry (qXL-MS) is an emerging technology able to provide information on changes in protein conformations and interactions. Importantly, qXL-MS is applicable to complex biological systems, including living cells and tissues, thereby providing insights into proteins within their native environments. Here, we present an overview of recent technological developments and applications involving qXL-MS, including design and synthesis of isotope-labeled cross-linkers, development of new liquid chromatography–MS methodologies, and computational developments enabling interpretation of the data.