Trans-Proteomic Pipeline: Robust Mass Spectrometry-Based Proteomics Data Analysis Suite

iSanXoT: A standalone application for the integrative analysis of mass spectrometry-based quantitative proteomics data

Many bioinformatics tools are available for the quantitative analysis of proteomics experiments. Most of these tools use a dedicated statistical model to derive absolute quantitative protein values from mass spectrometry (MS) data. Here, we present iSanXoT, a standalone application that processes relative abundances between MS signals and then integrates them sequentially to upper levels using the previously published Generic Integration Algorithm (GIA). iSanXoT offers unique capabilities that complement conventional quantitative software applications, including statistical weighting and independent modeling of error distributions in each integration, aggregation of technical or biological replicates, quantification of posttranslational modifications, and analysis of coordinated protein behavior. iSanXoT is a standalone, user-friendly application that accepts output from popular proteomics pipelines and enables unrestricted creation of quantification workflows and fully customizable reports that can be reused across projects or shared among users. Numerous publications attest the successful application of diverse integrative workflows constructed using the GIA for the analysis of high-throughput quantitative proteomics experiments. iSanXoT has been tested with the main operating systems. Download links for the corresponding distributions are available at https://github.com/CNIC-Proteomics/iSanXoT/releases.

Comparative analysis of the longissimus muscle proteome of European wild boar and domestic pig in response to thermal processing

This study tries to fill the knowledge gap regarding differences in the expression of proteins in the meat of European wild boar (Sus scrofa scrofa) and domestic pig (Sus scrofa domestica), considering the impact of thermally induced degradation. We assessed relative protein changes between cooked longissimus thoracis et lumborum (LTL) muscle proteomes by using mass spectrometry, chemometric, label-free proteomic, and bioinformatic tools. Among 30 differentially abundant proteins identified MyHC-2a, ATPs-α, CK-S, ADP/ATPt1, IDH2, and MyBP-C1 were upregulated (x > 1) whereas NEB, γ-ENO and EPSF were downregulated (x < 1) in wild boar. ShinyGO and KEGG database pathway analyses revealed that these proteins are mainly involved in processes related to muscle contraction and various pathways of glucose metabolism and energy production. Protein expression changes could have been caused by the different muscle activity of wild animals in response to prolonged movement associated with foraging for food in the natural environment.

Metaproteomic analysis from cervical biopsies and cytologies identifies proteinaceous biomarkers representing both human and microbial species

The detection of HPV infection and microbial colonization in cervical lesions is currently done through PCR-based viral or bacterial DNA amplification. Our objective was to develop a methodology to expand the metaproteomic landscape of cervical disease and determine if protein biomarkers from both human and microbes could be detected in distinct cervical samples. This would lead to the development of multi-species proteomics, which includes protein-based lateral flow diagnostics that can define patterns of microbes and/or human proteins relevant to disease status. In this study, we collected both non-frozen tissue biopsy and exfoliative non-fixed cytology samples to assess the consistency of detecting human proteomic signatures between the cytology and biopsy samples. Our results show that proteomics using biopsies or cytologies can detect both human and microbial organisms. Across patients, Lumican and Galectin-1 were most highly expressed human proteins in the tissue biopsy, whilst IL-36 and IL-1RA were most highly expressed human proteins in the cytology. We also used mass spectrometry to assess microbial proteomes known to reside based on prior 16S rRNA gene signatures. Lactobacillus spp. was the most highly expressed proteome in patient samples and specific abundant Lactobacillus proteins were identified. These methodological approaches can be used in future metaproteomic clinical studies to interrogate the vaginal human and microbiome structure and metabolic diversity in cytologies or biopsies from the same patients who have pre-invasive cervical intraepithelial neoplasia, invasive cervical cancer, as well as in healthy controls to assess how human and pathogenic proteins may correlate with disease presence and severity.

Fish species authentication in commercial fish products using mass spectrometry and spectral library matching approach

Seafood fraud has become a global issue, threatening food security and safety. Adulteration, substitution, dilution, and incorrect labeling of seafood products are fraudulent practices that violate consumer safety. In this context, developing sensitive, robust, and high-throughput molecular tools for food and feed authentication is becoming crucial for regulatory purposes. Analytical approaches such as proteomics mass spectrometry have shown promise in detecting incorrectly labeled products. For the application of these tools, genome information is crucial, but currently, for many marine species of commercial importance, such information is unavailable. However, when combining proteomic analysis with spectral library matching, commercially important fish species were successfully identified, differentiated, and quantified in pure muscle samples and mixtures, even when genome information was scarce. This study further tested the previously developed spectral library matching approach to differentiate between 29 fish species from the North Sea and examined samples including individual fish, laboratory-prepared mixtures and commercial products. For authenticating libraries generated from 29 fish species, fresh muscle samples from the fish samples were matched against the reference spectral libraries. Species of the fresh fish samples were correctly authenticated using the spectral library approach. The same result was obtained when evaluating the laboratory-prepared mixtures. Furthermore, processed commercial products containing mixtures of two or three fish species were matched against these reference spectral libraries to test the accuracy and robustness of this method for authentication of fish species. The results indicated that the method is suitable for the authentication of fish species from highly processed samples such as fish cakes and burgers. The study shows that current and future challenges in food and feed authentication can efficiently be tackled by reference spectral libraries method when prospecting new resources in the Arctic.

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.

Functional genomic analysis of genes important for Candida albicans fitness in diverse environmental conditions

Fungal pathogens such as Candida albicans pose a significant threat to human health with limited treatment options available. One strategy to expand the therapeutic target space is to identify genes important for pathogen growth in host-relevant environments. Here, we leverage a pooled functional genomic screening strategy to identify genes important for fitness of C. albicans in diverse conditions. We identify an essential gene with no known Saccharomyces cerevisiae homolog, C1_09670C, and demonstrate that it encodes subunit 3 of replication factor A (Rfa3). Furthermore, we apply computational analyses to identify functionally coherent gene clusters and predict gene function. Through this approach, we predict the cell-cycle-associated function of C3_06880W, a previously uncharacterized gene required for fitness specifically at elevated temperatures, and follow-up assays confirm that C3_06880W encodes Iml3, a component of the C. albicans kinetochore with roles in virulence in vivo. Overall, this work reveals insights into the vulnerabilities of C. albicans.

High-Performance Workflow for Identifying Site-Specific Crosslinks Originating from a Genetically Incorporated, Photoreactive Amino Acid

In conventional crosslinking mass spectrometry, proteins are crosslinked using a highly selective, bifunctional chemical reagent, which limits crosslinks to residues that are accessible and reactive to the reagent. Genetically incorporating a photoreactive amino acid offers two key advantages: any site can be targeted, including those that are inaccessible to conventional crosslinking reagents, and photoreactive amino acids can potentially react with a broad range of interaction partners. However, broad reactivity imposes additional challenges for crosslink identification. In this study, we incorporate benzoylphenylalanine (BPA), a photoreactive amino acid, at selected sites in an intrinsically disordered region of the human protein HSPB5. We report and characterize a workflow for identifying and visualizing residue-level interactions originating from BPA. We routinely identify 30 to 300 crosslinked peptide spectral matches with this workflow, which is up to ten times more than existing tools for residue-level BPA crosslink identification. Most identified crosslinks are assigned to a precision of one or two residues, which is supported by a high degree of overlap between replicate analyses. Based on these results, we anticipate that this workflow will support the more general use of genetically incorporated, photoreactive amino acids for characterizing the structures of proteins that have resisted high-resolution characterization.

Meta-Analysis of Rice Phosphoproteomics Data to Understand Variation in Cell Signaling Across the Rice Pan-Genome

Phosphorylation is the most studied post-translational modification, and has multiple biological functions. In this study, we have reanalyzed publicly available mass spectrometry proteomics data sets enriched for phosphopeptides from Asian rice (Oryza sativa). In total we identified 15,565 phosphosites on serine, threonine, and tyrosine residues on rice proteins. We identified sequence motifs for phosphosites, and link motifs to enrichment of different biological processes, indicating different downstream regulation likely caused by different kinase groups. We cross-referenced phosphosites against the rice 3,000 genomes, to identify single amino acid variations (SAAVs) within or proximal to phosphosites that could cause loss of a site in a given rice variety and clustered the data to identify groups of sites with similar patterns across rice family groups. The data has been loaded into UniProt Knowledge-Base─enabling researchers to visualize sites alongside other data on rice proteins, e.g., structural models from AlphaFold2, PeptideAtlas, and the PRIDE database─enabling visualization of source evidence, including scores and supporting mass spectra.

A eukaryotic-like ubiquitination system in bacterial antiviral defence

Ubiquitination pathways have crucial roles in protein homeostasis, signalling and innate immunity1-3. In these pathways, an enzymatic cascade of E1, E2 and E3 proteins conjugates ubiquitin or a ubiquitin-like protein (Ubl) to target-protein lysine residues4. Bacteria encode ancient relatives of E1 and Ubl proteins involved in sulfur metabolism5,6, but these proteins do not mediate Ubl-target conjugation, leaving open the question of whether bacteria can perform ubiquitination-like protein conjugation. Here we demonstrate that a bacterial operon associated with phage defence islands encodes a complete ubiquitination pathway. Two structures of a bacterial E1-E2-Ubl complex reveal striking architectural parallels with canonical eukaryotic ubiquitination machinery. The bacterial E1 possesses an amino-terminal inactive adenylation domain and a carboxy-terminal active adenylation domain with a mobile α-helical insertion containing the catalytic cysteine (CYS domain). One structure reveals a pre-reaction state with the bacterial Ubl C terminus positioned for adenylation, and a second structure mimics an E1-to-E2 transthioesterification state with the E1 CYS domain adjacent to the bound E2. We show that a deubiquitinase in the same pathway preprocesses the bacterial Ubl, exposing its C-terminal glycine for adenylation. Finally, we show that the bacterial E1 and E2 collaborate to conjugate Ubl to target-protein lysine residues. Together, these data reveal that bacteria possess bona fide ubiquitination systems with strong mechanistic and architectural parallels to canonical eukaryotic ubiquitination pathways, suggesting that these pathways arose first in bacteria.

MS-PyCloud: A Cloud Computing-Based Pipeline for Proteomic and Glycoproteomic Data Analyses

Rapid development and wide adoption of mass spectrometry-based glycoproteomic technologies have empowered scientists to study proteins and protein glycosylation in complex samples on a large scale. This progress has also created unprecedented challenges for individual laboratories to store, manage, and analyze proteomic and glycoproteomic data, both in the cost for proprietary software and high-performance computing and in the long processing time that discourages on-the-fly changes of data processing settings required in explorative and discovery analysis. We developed an open-source, cloud computing-based pipeline, MS-PyCloud, with graphical user interface (GUI), for proteomic and glycoproteomic data analysis. The major components of this pipeline include data file integrity validation, MS/MS database search for spectral assignments to peptide sequences, false discovery rate estimation, protein inference, quantitation of global protein levels, and specific glycan-modified glycopeptides as well as other modification-specific peptides such as phosphorylation, acetylation, and ubiquitination. To ensure the transparency and reproducibility of data analysis, MS-PyCloud includes open-source software tools with comprehensive testing and versioning for spectrum assignments. Leveraging public cloud computing infrastructure via Amazon Web Services (AWS), MS-PyCloud scales seamlessly based on analysis demand to achieve fast and efficient performance. Application of the pipeline to the analysis of large-scale LC-MS/MS data sets demonstrated the effectiveness and high performance of MS-PyCloud. The software can be downloaded at https://github.com/huizhanglab-jhu/ms-pycloud.

A phosphorylation-controlled switch confers cell cycle-dependent protein relocalization

Tools for acute manipulation of protein localization enable elucidation of spatiotemporally defined functions, but their reliance on exogenous triggers can interfere with cell physiology. This limitation is particularly apparent for studying mitosis, whose highly choreographed events are sensitive to perturbations. Here we exploit the serendipitous discovery of a phosphorylation-controlled, cell cycle-dependent localization change of the adaptor protein PLEKHA5 to develop a system for mitosis-specific protein recruitment to the plasma membrane that requires no exogenous stimulus. Mitosis-enabled Anchor-away/Recruiter System (MARS) comprises an engineered, 15-kDa module derived from PLEKHA5 capable of recruiting functional protein cargoes to the plasma membrane during mitosis, either through direct fusion or via GFP-GFP nanobody interaction. Applications of MARS include both knock sideways to rapidly extract proteins from their native localizations during mitosis and conditional recruitment of lipid-metabolizing enzymes for mitosis-selective editing of plasma membrane lipid content, without the need for exogenous triggers or perturbative synchronization methods.

CLIPPER 2.0: Peptide-Level Annotation and Data Analysis for Positional Proteomics

Positional proteomics methodologies have transformed protease research, and have brought mass spectrometry (MS)-based degradomics studies to the forefront of protease characterization and system-wide interrogation of protease signaling. Considerable advancements in both sensitivity and throughput of liquid chromatography (LC)-MS/MS instrumentation enable the generation of enormous positional proteomics datasets of natural and protein termini and neo-termini of cleaved protease substrates. However, concomitant progress has not been observed to the same extent in data analysis and post-processing steps, arguably constituting the largest bottleneck in positional proteomics workflows. Here, we present a computational tool, CLIPPER 2.0, that builds on prior algorithms developed for MS-based protein termini analysis, facilitating peptide-level annotation and data analysis. CLIPPER 2.0 can be used with several sample preparation workflows and proteomics search algorithms and enables fast and automated database information retrieval, statistical and network analysis, as well as visualization of terminomic datasets. We demonstrate the applicability of our tool by analyzing GluC and MMP9 cleavages in HeLa lysates. CLIPPER 2.0 is available at https://github.com/UadKLab/CLIPPER-2.0.

LiF-MS+, a revised technique for mapping peptide-protein interactions

Short linear motifs are sequences of amino acids present in unstructured polypeptide regions that function as ligands for specific sites on folded protein domains. These interactions, which often occur with low to modest affinity, modulate dynamic biological processes such as signal transduction and membrane trafficking. We recently described Ligand Footprinting-Mass Spectrometry (LiF-MS), a technique that rapidly and precisely maps sites at which short peptide ligands bind their biologically relevant recognition sites on folded protein domains. This approach marks the binding location of a peptide ligand on a structured protein using a cleavable crosslinker appended to the ligand that leaves behind a stable chemical modification following cleavage. This modification serves as a mass tag detectable by mass spectrometry, pinpointing sites of peptide ligand binding. Here we present LiF-MS+, an improved version of the footprinting technique that replaces the butanol mass tag with 1-butylpyrrolidine, which is positively charged at neutral pH and thus aids in ionization of the crosslinked peptide for analysis by mass spectrometry. We show ligand-mediated butylpyrrolidine footprinting effectively maps the well characterized binding interaction of the p38α mitogen-activated protein kinase (MAPK) with a MKK6 D-motif short linear motif peptide ligand, uncovering additional binding site information not observed in our original experiment. LiF-MS+ is thus a straightforward improvement of our previously published methodology for mapping the binding of short linear motifs to folded protein domains.

SpectiCal: m/z Calibration of MS2 Peptide Spectra Using Known Low Mass Ions

Most tandem mass spectrometry fragmentation spectra have small calibration errors that can lead to suboptimal interpretation and annotation. We developed SpectiCal, a software tool that can read mzML files from data-dependent acquisition proteomics experiments in parallel, compute m/z calibrations for each file prior to identification analysis based on known low-mass ions, and produce information about frequently observed peaks and their explanations. Using calibration coefficients, the data can be corrected to generate new calibrated mzML files. SpectiCal was tested using five public data sets, creating a table of commonly observed low-mass ions and their identifications. Information about the calibration and individual peaks is written in PDF and TSV files. This includes information for each peak, such as the number of runs in which it appears, the percentage of spectra in which it appears, and a plot of the aggregated region surrounding each peak. SpectiCal can be used to compute MS run calibrations, examine MS runs for artifacts that might hinder downstream analysis, and generate tables of detected low-mass ions for further analysis. SpectiCal is freely available at https://github.com/PlantProteomes/SpectiCal.

Diverging co-translational protein complex assembly pathways are governed by interface energy distribution

Protein-protein interactions are at the heart of all cellular processes, with the ribosome emerging as a platform, orchestrating the nascent-chain interplay dynamics. Here, to study the characteristics governing co-translational protein folding and complex assembly, we combine selective ribosome profiling, imaging, and N-terminomics with all-atoms molecular dynamics. Focusing on conserved N-terminal acetyltransferases (NATs), we uncover diverging co-translational assembly pathways, where highly homologous subunits serve opposite functions. We find that only a few residues serve as "hotspots," initiating co-translational assembly interactions upon exposure at the ribosome exit tunnel. These hotspots are characterized by high binding energy, anchoring the entire interface assembly. Alpha-helices harboring hotspots are highly thermolabile, folding and unfolding during simulations, depending on their partner subunit to avoid misfolding. In vivo hotspot mutations disrupted co-translational complexation, leading to aggregation. Accordingly, conservation analysis reveals that missense NATs variants, causing neurodevelopmental and neurodegenerative diseases, disrupt putative hotspot clusters. Expanding our study to include phosphofructokinase, anthranilate synthase, and nucleoporin subcomplex, we employ AlphaFold-Multimer to model the complexes' complete structures. Computing MD-derived interface energy profiles, we find similar trends. Here, we propose a model based on the distribution of interface energy as a strong predictor of co-translational assembly.

Global Release of Translational Repression Across Plasmodium's Host-to-Vector Transmission Event

Malaria parasites must be able to respond quickly to changes in their environment, including during their transmission between mammalian hosts and mosquito vectors. Therefore, before transmission, female gametocytes proactively produce and translationally repress mRNAs that encode essential proteins that the zygote requires to establish a new infection. This essential regulatory control requires the orthologues of DDX6 (DOZI), LSM14a (CITH), and ALBA proteins to form a translationally repressive complex in female gametocytes that associates with many of the affected mRNAs. However, while the release of translational repression of individual mRNAs has been documented, the details of the global release of translational repression have not. Moreover, the changes in spatial arrangement and composition of the DOZI/CITH/ALBA complex that contribute to translational control are also not known. Therefore, we have conducted the first quantitative, comparative transcriptomics and DIA-MS proteomics of Plasmodium parasites across the host-to-vector transmission event to document the global release of translational repression. Using female gametocytes and zygotes of P. yoelii, we found that nearly 200 transcripts are released for translation soon after fertilization, including those with essential functions for the zygote. However, we also observed that some transcripts remain repressed beyond this point. In addition, we have used TurboID-based proximity proteomics to interrogate the spatial and compositional changes in the DOZI/CITH/ALBA complex across this transmission event. Consistent with recent models of translational control, proteins that associate with either the 5' or 3' end of mRNAs are in close proximity to one another during translational repression in female gametocytes and then dissociate upon release of repression in zygotes. This observation is cross-validated for several protein colocalizations in female gametocytes via ultrastructure expansion microscopy and structured illumination microscopy. Moreover, DOZI exchanges its interaction from NOT1-G in female gametocytes to the canonical NOT1 in zygotes, providing a model for a trigger for the release of mRNAs from DOZI. Finally, unenriched phosphoproteomics revealed the modification of key translational control proteins in the zygote. Together, these data provide a model for the essential translational control mechanisms used by malaria parasites to promote their efficient transmission from their mammalian host to their mosquito vector.

Polyphosphate kinase regulates LPS structure and polymyxin resistance during starvation in E. coli

Polyphosphates (polyP) are chains of inorganic phosphates that can reach over 1,000 residues in length. In Escherichia coli, polyP is produced by the polyP kinase (PPK) and is thought to play a protective role during the response to cellular stress. However, the molecular pathways impacted by PPK activity and polyP accumulation remain poorly characterized. In this work, we used label-free mass spectrometry to study the response of bacteria that cannot produce polyP (Δppk) during starvation to identify novel pathways regulated by PPK. In response to starvation, we found 92 proteins significantly differentially expressed between wild-type and Δppk mutant cells. Wild-type cells were enriched for proteins related to amino acid biosynthesis and transport, while Δppk mutants were enriched for proteins related to translation and ribosome biogenesis, suggesting that without PPK, cells remain inappropriately primed for growth even in the absence of the required building blocks. From our data set, we were particularly interested in Arn and EptA proteins, which were down-regulated in Δppk mutants compared to wild-type controls, because they play a role in lipid A modifications linked to polymyxin resistance. Using western blotting, we confirm differential expression of these and related proteins in K-12 strains and a uropathogenic isolate, and provide evidence that this mis-regulation in Δppk cells stems from a failure to induce the BasRS two-component system during starvation. We also show that Δppk mutants unable to up-regulate Arn and EptA expression lack the respective L-Ara4N and pEtN modifications on lipid A. In line with this observation, loss of ppk restores polymyxin sensitivity in resistant strains carrying a constitutively active basR allele. Overall, we show a new role for PPK in lipid A modification during starvation and provide a rationale for targeting PPK to sensitize bacteria towards polymyxin treatment. We further anticipate that our proteomics work will provide an important resource for researchers interested in the diverse pathways impacted by PPK.

The T2T-CHM13 reference assembly uncovers essential WASH1 and GPRIN2 paralogues

Summary

The recently published T2T-CHM13 reference assembly completed the annotation of the final 8% of the human genome. It introduced 1956 genes, close to 100 of which are predicted to be coding because they have a protein coding parent gene. Here, we confirm the coding status and functional relevance of two of these genes, paralogues of WASHC1 and GPRIN2. We find that LOC124908094, one of four novel subtelomeric WASH1 genes uncovered in the new assembly, produces the WASH1 protein that forms part of the vital actin-regulatory WASH complex. Its coding status is supported by abundant proteomics, conservation, and cDNA evidence. It was previously assumed that gene WASHC1 produced the functional WASH1 protein, but new evidence shows that WASHC1 is a human-derived duplication and likely to be one of 12 WASH1 pseudogenes in the human gene set. We also find that the T2T-CHM13 assembly has added a functionally important copy of GPRIN2 to the human gene set. We demonstrate that uniquely mapping peptides from proteomics databases support the novel LOC124900631 rather than the GRCh38 assembly GPRIN2 gene. These new additions to the set of human coding genes underlines the importance of the new T2T-CHM13 assembly.

Availability and implementation

None.

The 2023 Report on the Proteome from the HUPO Human Proteome Project

Since 2010, the Human Proteome Project (HPP), the flagship initiative of the Human Proteome Organization (HUPO), has pursued two goals: (1) to credibly identify the protein parts list and (2) to make proteomics an integral part of multiomics studies of human health and disease. The HPP relies on international collaboration, data sharing, standardized reanalysis of MS data sets by PeptideAtlas and MassIVE-KB using HPP Guidelines for quality assurance, integration and curation of MS and non-MS protein data by neXtProt, plus extensive use of antibody profiling carried out by the Human Protein Atlas. According to the neXtProt release 2023-04-18, protein expression has now been credibly detected (PE1) for 18,397 of the 19,778 neXtProt predicted proteins coded in the human genome (93%). Of these PE1 proteins, 17,453 were detected with mass spectrometry (MS) in accordance with HPP Guidelines and 944 by a variety of non-MS methods. The number of neXtProt PE2, PE3, and PE4 missing proteins now stands at 1381. Achieving the unambiguous identification of 93% of predicted proteins encoded from across all chromosomes represents remarkable experimental progress on the Human Proteome parts list. Meanwhile, there are several categories of predicted proteins that have proved resistant to detection regardless of protein-based methods used. Additionally there are some PE1-4 proteins that probably should be reclassified to PE5, specifically 21 LINC entries and ∼30 HERV entries; these are being addressed in the present year. Applying proteomics in a wide array of biological and clinical studies ensures integration with other omics platforms as reported by the Biology and Disease-driven HPP teams and the antibody and pathology resource pillars. Current progress has positioned the HPP to transition to its Grand Challenge Project focused on determining the primary function(s) of every protein itself and in networks and pathways within the context of human health and disease.

Acquisition and Analysis of DIA-Based Proteomic Data: A Comprehensive Survey in 2023

Data-independent acquisition (DIA) mass spectrometry (MS) has emerged as a powerful technology for high-throughput, accurate, and reproducible quantitative proteomics. This review provides a comprehensive overview of recent advances in both the experimental and computational methods for DIA proteomics, from data acquisition schemes to analysis strategies and software tools. DIA acquisition schemes are categorized based on the design of precursor isolation windows, highlighting wide-window, overlapping-window, narrow-window, scanning quadrupole-based, and parallel accumulation-serial fragmentation-enhanced DIA methods. For DIA data analysis, major strategies are classified into spectrum reconstruction, sequence-based search, library-based search, de novo sequencing, and sequencing-independent approaches. A wide array of software tools implementing these strategies are reviewed, with details on their overall workflows and scoring approaches at different steps. The generation and optimization of spectral libraries, which are critical resources for DIA analysis, are also discussed. Publicly available benchmark datasets covering global proteomics and phosphoproteomics are summarized to facilitate performance evaluation of various software tools and analysis workflows. Continued advances and synergistic developments of versatile components in DIA workflows are expected to further enhance the power of DIA-based proteomics.

WOMBAT-P: Benchmarking Label-Free Proteomics Data Analysis Workflows

The inherent diversity of approaches in proteomics research has led to a wide range of software solutions for data analysis. These software solutions encompass multiple tools, each employing different algorithms for various tasks such as peptide-spectrum matching, protein inference, quantification, statistical analysis, and visualization. To enable an unbiased comparison of commonly used bottom-up label-free proteomics workflows, we introduce WOMBAT-P, a versatile platform designed for automated benchmarking and comparison. WOMBAT-P simplifies the processing of public data by utilizing the sample and data relationship format for proteomics (SDRF-Proteomics) as input. This feature streamlines the analysis of annotated local or public ProteomeXchange data sets, promoting efficient comparisons among diverse outputs. Through an evaluation using experimental ground truth data and a realistic biological data set, we uncover significant disparities and a limited overlap in the quantified proteins. WOMBAT-P not only enables rapid execution and seamless comparison of workflows but also provides valuable insights into the capabilities of different software solutions. These benchmarking metrics are a valuable resource for researchers in selecting the most suitable workflow for their specific data sets. The modular architecture of WOMBAT-P promotes extensibility and customization. The software is available at https://github.com/wombat-p/WOMBAT-Pipelines.

A multivariate data analysis approach for the investigation of in vitro derived metabolites of ACP-105 in comparison with human in vivo metabolites

Selective androgen receptor modulators (SARMs) such as ACP-105 are prohibited in sports due to their anabolic properties. ACP-105 has in previous equine studies shown to undergo extensive metabolism, which makes its metabolite profile important to investigate in humans, since the metabolism is unknown in this species. The aims of the study were to systematically optimize in vitro microsome incubations for improved metabolite yield and to utilize a multivariate data analysis (MVDA) approach to aid the metabolite discovery. Microsomes together with S9 fractions were used at optimal conditions, both with and without phase II additives. Furthermore, the relevance of the in vitro derived metabolites was evaluated as analytical targets in doping control by comparison with results from a human post-administration urine sample collected after a single dose of 100 µg ACP-105. All samples were analyzed with liquid chromatography - Orbitrap mass spectrometry. The use of the systematical optimization and MVDA greatly simplified the search and a total of 18 in vitro metabolites were tentatively identified. The yield of the two main monohydroxylated isomers increased by 24 and 10 times, respectively. In the human urine sample, a total of seven metabolites of ACP-105, formed by a combination of hydroxylations and glucuronic acid conjugations, were tentatively identified. The main metabolites were two monohydroxylated forms that are suggested as analytical targets for human doping control after hydrolysis. All the in vivo metabolites could be detected with the MVDA approach on the in vitro models, demonstrating its usefulness for prediction of the in vivo metabolite profile.

A meta-analysis of rice phosphoproteomics data to understand variation in cell signalling across the rice pan-genome

Phosphorylation is the most studied post-translational modification, and has multiple biological functions. In this study, we have re-analysed publicly available mass spectrometry proteomics datasets enriched for phosphopeptides from Asian rice (Oryza sativa). In total we identified 15,522 phosphosites on serine, threonine and tyrosine residues on rice proteins. We identified sequence motifs for phosphosites, and link motifs to enrichment of different biological processes, indicating different downstream regulation likely caused by different kinase groups. We cross-referenced phosphosites against the rice 3,000 genomes, to identify single amino acid variations (SAAVs) within or proximal to phosphosites that could cause loss of a site in a given rice variety. The data was clustered to identify groups of sites with similar patterns across rice family groups, for example those highly conserved in Japonica, but mostly absent in Aus type rice varieties - known to have different responses to drought. These resources can assist rice researchers to discover alleles with significantly different functional effects across rice varieties. The data has been loaded into UniProt Knowledge-Base - enabling researchers to visualise sites alongside other data on rice proteins e.g. structural models from AlphaFold2, PeptideAtlas and the PRIDE database - enabling visualisation of source evidence, including scores and supporting mass spectra.

A comprehensive spectral assay library to quantify the Halobacterium salinarum NRC-1 proteome by DIA/SWATH-MS

Data-Independent Acquisition (DIA) is a mass spectrometry-based method to reliably identify and reproducibly quantify large fractions of a target proteome. The peptide-centric data analysis strategy employed in DIA requires a priori generated spectral assay libraries. Such assay libraries allow to extract quantitative data in a targeted approach and have been generated for human, mouse, zebrafish, E. coli and few other organisms. However, a spectral assay library for the extreme halophilic archaeon Halobacterium salinarum NRC-1, a model organism that contributed to several notable discoveries, is not publicly available yet. Here, we report a comprehensive spectral assay library to measure 2,563 of 2,646 annotated H. salinarum NRC-1 proteins. We demonstrate the utility of this library by measuring global protein abundances over time under standard growth conditions. The H. salinarum NRC-1 library includes 21,074 distinct peptides representing 97% of the predicted proteome and provides a new, valuable resource to confidently measure and quantify any protein of this archaeon. Data and spectral assay libraries are available via ProteomeXchange (PXD042770, PXD042774) and SWATHAtlas (SAL00312-SAL00319).

Spectroscape enables real-time query and visualization of a spectral archive in proteomics

In proteomics, spectral archives organize the enormous amounts of publicly available peptide tandem mass spectra by similarity, offering opportunities for error correction and novel discoveries. Here we adapt an indexing algorithm developed by Facebook for organizing online multimedia resources to tandem mass spectra and achieve practically instantaneous retrieval and clustering of approximate nearest neighbors in a large spectral archive. An interactive web-based graphical user interface enables the user to view a query spectrum in its clustered neighborhood, which facilitates contextual validation of peptide identifications and exploration of the dark proteome.

Bacterial antiviral defense pathways encode eukaryotic-like ubiquitination systems

Ubiquitination and related pathways play crucial roles in protein homeostasis, signaling, and innate immunity1-3. In these pathways, an enzymatic cascade of E1, E2, and E3 proteins conjugates ubiquitin or a ubiquitin-like protein (Ubl) to target-protein lysine residues4. Bacteria encode ancient relatives of E1 and Ubl proteins involved in sulfur metabolism5,6 but these proteins do not mediate Ubl-target conjugation, leaving open the question of whether bacteria can perform ubiquitination-like protein conjugation. Here, we demonstrate that a bacterial antiviral immune system encodes a complete ubiquitination pathway. Two structures of a bacterial E1:E2:Ubl complex reveal striking architectural parallels with canonical eukaryotic ubiquitination machinery. The bacterial E1 encodes an N-terminal inactive adenylation domain (IAD) and a C-terminal active adenylation domain (AAD) with a mobile α-helical insertion containing the catalytic cysteine (CYS domain). One structure reveals a pre-reaction state with the bacterial Ubl C-terminus positioned for adenylation, and the E1 CYS domain poised nearby for thioester formation. A second structure mimics an E1-to-E2 transthioesterification state, with the E1 CYS domain rotated outward and its catalytic cysteine adjacent to the bound E2. We show that a deubiquitinase (DUB) in the same pathway pre-processes the bacterial Ubl, exposing its C-terminal glycine for adenylation. Finally, we show that the bacterial E1 and E2 collaborate to conjugate Ubl to target-protein lysine residues. Together, these data reveal that bacteria possess bona fide ubiquitination systems with strong mechanistic and architectural parallels to canonical eukaryotic ubiquitination pathways, suggesting that these pathways arose first in bacteria.

Database Creator for Mass Analysis of Peptides and Proteins, DC-MAPP: A Standalone Tool for Simplifying Manual Analysis of Mass Spectral Data to Identify Peptide/Protein Sequences

Proteomic studies typically involve the use of different types of software for annotating experimental tandem mass spectrometric data (MS/MS) and thereby simplifying the process of peptide and protein identification. For such annotations, these softwares calculate the m/z values of the peptide/protein precursor and fragment ions, for which a database of protein sequences must be provided as an input file. The calculated m/z values are stored as another database, which the user usually cannot view. Database Creator for Mass Analysis of Peptides and Proteins (DC-MAPP) is a novel standalone software that can create custom databases for "viewing" the calculated m/z values of precursor and fragment ions, prior to the database search. It contains three modules. Peptide/Protein sequences as per user's choice can be entered as input to the first module for creating a custom database. In the second module, m/z values must be queried-in, which are searched within the custom database to identify protein/peptide sequences. The third module is suited for peptide mass fingerprinting, which can be used to analyze both ESI and MALDI mass spectral data. The feature of "viewing" the custom database can be helpful not only for better understanding the search engine processes, but also for designing multiple reaction monitoring (MRM) methods. Post-translational modifications and protein isoforms can also be analyzed. Since, DC-MAPP relies on the protein/peptide "sequences" for creating custom databases, it may not be applicable for the searches involving spectral libraries. Python language was used for implementation, and the graphical user interface was built with Page/Tcl, making this tool more user-friendly. It is freely available at https://vit.ac.in/DC-MAPP/.

XCP1 cleaves Pathogenesis-related protein 1 into CAPE9 for systemic immunity in Arabidopsis

Proteolytic activation of cytokines regulates immunity in diverse organisms. In animals, cysteine-dependent aspartate-specific proteases (caspases) play central roles in cytokine maturation. Although the proteolytic production of peptide cytokines is also essential for plant immunity, evidence for cysteine-dependent aspartate-specific proteases in regulating plant immunity is still limited. In this study, we found that the C-terminal proteolytic processing of a caspase-like substrate motif "CNYD" within Pathogenesis-related protein 1 (PR1) generates an immunomodulatory cytokine (CAPE9) in Arabidopsis. Salicylic acid enhances CNYD-targeted protease activity and the proteolytic release of CAPE9 from PR1 in Arabidopsis. This process involves a protease exhibiting caspase-like enzyme activity, identified as Xylem cysteine peptidase 1 (XCP1). XCP1 exhibits a calcium-modulated pH-activity profile and a comparable activity to human caspases. XCP1 is required to induce systemic immunity triggered by pathogen-associated molecular patterns. This work reveals XCP1 as a key protease for plant immunity, which produces the cytokine CAPE9 from the canonical salicylic acid signaling marker PR1 to activate systemic immunity.

Proteomics analysis reveals a role for E. coli polyphosphate kinase in membrane structure and polymyxin resistance during starvation

Polyphosphates (polyP) are chains of inorganic phosphates that can reach over 1000 residues in length. In Escherichia coli, polyP is produced by the polyP kinase (PPK) and is thought to play a protective role during the response to cellular stress. However, the molecular pathways impacted by PPK activity and polyP accumulation remain poorly characterized. In this work we used label-free mass spectrometry to study the response of bacteria that cannot produce polyP (∆ppk) during starvation to identify novel pathways regulated by PPK. In response to starvation, we found 92 proteins significantly differentially expressed between wild-type and ∆ppk mutant cells. Wild-type cells were enriched for proteins related to amino acid biosynthesis and transport, while Δppk mutants were enriched for proteins related to translation and ribosome biogenesis, suggesting that without PPK, cells remain inappropriately primed for growth even in the absence of required building blocks. From our dataset, we were particularly interested in Arn and EptA proteins, which were downregulated in ∆ppk mutants compared to wild-type controls, because they play a role in lipid A modifications linked to polymyxin resistance. Using western blotting, we confirm differential expression of these and related proteins, and provide evidence that this mis-regulation in ∆ppk cells stems from a failure to induce the BasS/BasR two-component system during starvation. We also show that ∆ppk mutants unable to upregulate Arn and EptA expression lack the respective L-Ara4N and pEtN modifications on lipid A. In line with this observation, loss of ppk restores polymyxin sensitivity in resistant strains carrying a constitutively active basR allele. Overall, we show a new role for PPK in lipid A modification during starvation and provide a rationale for targeting PPK to sensitize bacteria towards polymyxin treatment. We further anticipate that our proteomics work will provide an important resource for researchers interested in the diverse pathways impacted by PPK.

Lost in the WASH. The functional human WASH complex 1 gene is on chromosome 20

The WASH1 gene produces a protein that forms part of the developmentally important WASH complex. The WASH complex activates the Arp2/3 complex to initiate branched actin networks at the surface of endosomes. As a curiosity, the human reference gene set includes nine WASH1 genes. How many of these are pseudogenes and how many are bona fide coding genes is not clear. Eight of the nine WASH1 genes reside in rearrangement and duplication-prone subtelomeric regions. Many of these subtelomeric regions had gaps in the GRCh38 human genome assembly, but the recently published T2T-CHM13 assembly from the Telomere to Telomere (T2T) Consortium has filled in the gaps. As a result, the T2T Consortium has added four new WASH1 paralogues in previously unannotated subtelomeric regions. Here we show that one of these four novel WASH1 genes, LOC124908094, is the gene most likely to produce the functional WASH1 protein. We also demonstrate that the other twelve WASH1 genes derived from a single WASH8P pseudogene on chromosome 12. These 12 genes include WASHC1, the gene currently annotated as the functional WASH1 gene. We propose LOC124908094 should be annotated as a coding gene and all functional information relating to the WASHC1 gene on chromosome 9 should be transferred to LOC124908094. The remaining WASH1 genes, including WASHC1. should be annotated as pseudogenes. This work confirms that the T2T assembly has added at least one functionally relevant coding gene to the human reference set. It remains to be seen whether other important coding genes are missing from the GRCh38 reference assembly.

Does the Ubiquitination Degradation Pathway Really Reach inside of the Chloroplast? A Re-Evaluation of Mass Spectrometry-Based Assignments of Ubiquitination

A recent paper in Science Advances by Sun et al. claims that intra-chloroplast proteins in the model plant Arabidopsis can be polyubiquitinated and then extracted into the cytosol for subsequent degradation by the proteasome. Most of this conclusion hinges on several sets of mass spectrometry (MS) data. If the proposed results and conclusion are true, this would be a major change in the proteolysis/proteostasis field, breaking the long-standing dogma that there are no polyubiquitination mechanisms within chloroplast organelles (nor in mitochondria). Given its importance, we reanalyzed their raw MS data using both open and closed sequence database searches and encountered many issues not only with the results but also discrepancies between stated methods (e.g., use of alkylating agent iodoacetamide (IAA)) and observed mass modifications. Although there is likely enrichment of ubiquitination signatures in a subset of the data (probably from ubiquitination in the cytosol), we show that runaway alkylation with IAA caused extensive artifactual modifications of N termini and lysines to the point that a large fraction of the desired ubiquitination signatures is indistinguishable from artifactual acetamide signatures, and thus, no intra-chloroplast polyubiquitination conclusions can be drawn from these data. We provide recommendations on how to avoid such perils in future work.