MS-GF+ makes progress towards a universal database search tool for proteomics

PgxSAVy: A tool for comprehensive evaluation of variant peptide quality in proteogenomics - catching the (un)usual suspects

Variant peptides resulting from single nucleotide polymorphisms (SNPs) can lead to aberrant protein functions and have translational potential for disease diagnosis and personalized therapy. Variant peptides detected by proteogenomics are fraught with high number of false positives, but there is no uniform and comprehensive approach to assess variant quality across analysis pipelines. Despite class-specific FDR along with ad-hoc filters, the problem is far from solved. These protocols are typically manual and tedious, and thus not uniform across labs. We demonstrate that variant peptide rescoring, integrated with intensity, variant event information and search result features, allows better discrimination of correct variant peptides. Implemented into PgxSAVy - a tool for quality control of variant peptides, this method can tackle the high rate of false positives. PgxSAVy provides a rigorous framework for quality control and annotations of variant peptides on the basis of (i) variant quality, (ii) isobaric masses, and (iii) disease annotation. PgxSAVy demonstrated high accuracy by identifying true variants with 98.43% accuracy on simulated data. Large-scale proteogenomic reanalysis of ∼2.8 million spectra (PXD004010 and PXD001468) resulted in 12,705 variant peptide spectrum matches (PSMs), of which PgxSAVy evaluated 3028 (23.8%), 1409 (11.1%) and 8268 (65.1%) as confident, semi-confident and doubtful respectively. PgxSAVy also annotates the variants based on their pathogenicity and provides support for assisted manual validation. The analysis of proteins carrying variants can provide fine granularity in discovering important pathways. PgxSAVy will advance personalized medicine by providing a comprehensive framework for quality control and prioritization of proteogenomics variants. PgxSAVy is freely available at https://pgxsavy.igib.res.in/ as a webserver and https://github.com/anuragraj/PgxSAVy as a stand-alone tool.

MOSCA 2.0: A bioinformatics framework for metagenomics, metatranscriptomics and metaproteomics data analysis and visualization

The analysis of meta-omics data requires the utilization of several bioinformatics tools and proficiency in informatics. The integration of multiple meta-omics data is even more challenging, and the outputs of existing bioinformatics solutions are not always easy to interpret. Here, we present a meta-omics bioinformatics pipeline, Meta-Omics Software for Community Analysis (MOSCA), which aims to overcome these limitations. MOSCA was initially developed for analysing metagenomics (MG) and metatranscriptomics (MT) data. Now, it also performs MG and metaproteomics (MP) integrated analysis, and MG/MT analysis was upgraded with an additional iterative binning step, metabolic pathways mapping, and several improvements regarding functional annotation and data visualization. MOSCA handles raw sequencing data and mass spectra and performs pre-processing, assembly, annotation, binning and differential gene/protein expression analysis. MOSCA shows taxonomic and functional analysis in large tables, performs metabolic pathways mapping, generates Krona plots and shows gene/protein expression results in heatmaps, improving omics data visualization. MOSCA is easily run from a single command while also providing a web interface (MOSGUITO). Relevant features include an extensive set of customization options, allowing tailored analyses to suit specific research objectives, and the ability to restart the pipeline from intermediary checkpoints using alternative configurations. Two case studies showcased MOSCA results, giving a complete view of the anaerobic microbial communities from anaerobic digesters and insights on the role of specific microorganisms. MOSCA represents a pivotal advancement in meta-omics research, offering an intuitive, comprehensive, and versatile solution for researchers seeking to unravel the intricate tapestry of microbial communities.

pyRBDome: a comprehensive computational platform for enhancing RNA-binding proteome data

High-throughput proteomics approaches have revolutionised the identification of RNA-binding proteins (RBPome) and RNA-binding sequences (RBDome) across organisms. Yet, the extent of noise, including false positives, associated with these methodologies, is difficult to quantify as experimental approaches for validating the results are generally low throughput. To address this, we introduce pyRBDome, a pipeline for enhancing RNA-binding proteome data in silico. It aligns the experimental results with RNA-binding site (RBS) predictions from distinct machine-learning tools and integrates high-resolution structural data when available. Its statistical evaluation of RBDome data enables quick identification of likely genuine RNA-binders in experimental datasets. Furthermore, by leveraging the pyRBDome results, we have enhanced the sensitivity and specificity of RBS detection through training new ensemble machine-learning models. pyRBDome analysis of a human RBDome dataset, compared with known structural data, revealed that although UV-cross-linked amino acids were more likely to contain predicted RBSs, they infrequently bind RNA in high-resolution structures. This discrepancy underscores the limitations of structural data as benchmarks, positioning pyRBDome as a valuable alternative for increasing confidence in RBDome datasets.

Comparative proteomic analyses of potato leaves from field-grown plants grown under extremely long days

There are limited molecular data and few biomarkers available for studies of field-grown plants, especially for plants grown during extremely long days. In this study we present quantitative proteomics data from 3 years of field trials on potato, conducted in northern and southern Sweden and analyze over 3000 proteins per year of the study and complement the proteomic analysis with metabolomic and transcriptomic analyses. Small but consistent differences linked to the longer days (an average of four more hours of light per day) in northern Sweden (20 h light/day) compared to southern Sweden can be observed, with a high correlation between the mRNA determined by RNA-seq and protein abundances. The majority of the proteins with differential abundances between northern and southern Sweden could be divided into three groups: metabolic enzymes (especially GABA metabolism), proteins involved in redox metabolism, and hydrolytic enzymes. The observed differences in metabolic enzyme abundances corresponded well with untargeted metabolite data determined by GC and LC mass-spectrometry. We also analyzed differences in protein abundance between potato varieties that performed relatively well in northern Sweden in terms of yield with those that performed relatively less well. This comparison indicates that the proteins with higher abundance in the high-yield quotient group are more anabolic in their character, whereas the proteins with lower abundance are more catabolic. Our results create a base of information about potato "field-omics" for improved understanding of physiological and molecular processes in field-grown plants, and our data indicate that the potato plant is not generally stressed by extremely long days.

C-TERMINAL DOMAIN PHOSPHATASE-LIKE 1 promotes flowering with TAF15b by repressing the floral repressor gene FLOWERING LOCUS C

Arabidopsis TATA-BINDING PROTEIN-ASSOCIATED FACTOR15b (TAF15b) is a plant-specific component of the transcription factor IID complex. TAF15b is involved in the autonomous pathway for flowering and represses the transcription of FLOWERING LOCUS C (FLC), a major floral repressor in Arabidopsis. While components of the autonomous flowering pathway have been extensively studied, scant attention has been directed toward elucidating the direct transcriptional regulators responsible for repressing FLC transcription. Here, we demonstrate that C-TERMINAL DOMAIN PHOSPHATASE-LIKE 1 (CPL1) is a physical and functional partner of TAF15b, playing a role in FLC repression. CPL1 is a protein phosphatase that dephosphorylates the C-terminal domain of RNA polymerase II (Pol II). Through the immunoprecipitation and mass spectrometry technique, we identified CPL1 as an interacting partner of TAF15b. Similar to taf15b, the cpl1 mutant showed a late-flowering phenotype caused by an increase in FLC levels. Additionally, the increase in cpl1 was correlated with the enrichment of phosphorylated Pol II in the FLC chromatin, as expected. We also discovered that CPL1 and TAF15b share additional common target genes through transcriptome analysis. These results suggest that TAF15b and CPL1 cooperatively repress transcription through the dephosphorylation of Pol II, especially at the FLC locus.

metaExpertPro: A Computational Workflow for Metaproteomics Spectral Library Construction and Data-Independent Acquisition Mass Spectrometry Data Analysis

Analysis of large-scale data-independent acquisition mass spectrometry metaproteomics data remains a computational challenge. Here, we present a computational pipeline called metaExpertPro for metaproteomics data analysis. This pipeline encompasses spectral library generation using data-dependent acquisition MS, protein identification and quantification using data-independent acquisition mass spectrometry, functional and taxonomic annotation, as well as quantitative matrix generation for both microbiota and hosts. By integrating FragPipe and DIA-NN, metaExpertPro offers compatibility with both Orbitrap and timsTOF MS instruments. To evaluate the depth and accuracy of identification and quantification, we conducted extensive assessments using human fecal samples and benchmark tests. Performance tests conducted on human fecal samples indicated that metaExpertPro quantified an average of 45,000 peptides in a 60-min diaPASEF injection. Notably, metaExpertPro outperformed three existing software tools by characterizing a higher number of peptides and proteins. Importantly, metaExpertPro maintained a low factual false discovery rate of approximately 5% for protein groups across four benchmark tests. Applying a filter of five peptides per genus, metaExpertPro achieved relatively high accuracy (F-score = 0.67-0.90) in genus diversity and showed a high correlation (rSpearman = 0.73-0.82) between the measured and true genus relative abundance in benchmark tests. Additionally, the quantitative results at the protein, taxonomy, and function levels exhibited high reproducibility and consistency across the commonly adopted public human gut microbial protein databases IGC and UHGP. In a metaproteomic analysis of dyslipidemia patients, metaExpertPro revealed characteristic alterations in microbial functions and potential interactions between the microbiota and the host.

Molecular Profiling Defines Three Subtypes of Synovial Sarcoma

Synovial Sarcomas (SS) are characterized by the presence of the SS18::SSX fusion gene, which protein product induce chromatin changes through remodeling of the BAF complex. To elucidate the genomic events that drive phenotypic diversity in SS, we performed RNA and targeted DNA sequencing on 91 tumors from 55 patients. Our results were verified by proteomic analysis, public gene expression cohorts and single-cell RNA sequencing. Transcriptome profiling identified three distinct SS subtypes resembling the known histological subtypes: SS subtype I and was characterized by hyperproliferation, evasion of immune detection and a poor prognosis. SS subtype II and was dominated by a vascular-stromal component and had a significantly better outcome. SS Subtype III was characterized by biphasic differentiation, increased genomic complexity and immune suppression mediated by checkpoint inhibition, and poor prognosis despite good responses to neoadjuvant therapy. Chromosomal abnormalities were an independent significant risk factor for metastasis. KRT8 was identified as a key component for epithelial differentiation in biphasic tumors, potentially controlled by OVOL1 regulation. Our findings explain the histological grounds for SS classification and indicate that a significantly larger proportion of patients have high risk tumors (corresponding to SS subtype I) than previously believed.

Machine Learning Strategies to Tackle Data Challenges in Mass Spectrometry-Based Proteomics

In computational proteomics, machine learning (ML) has emerged as a vital tool for enhancing data analysis. Despite significant advancements, the diversity of ML model architectures and the complexity of proteomics data present substantial challenges in the effective development and evaluation of these tools. Here, we highlight the necessity for high-quality, comprehensive data sets to train ML models and advocate for the standardization of data to support robust model development. We emphasize the instrumental role of key data sets like ProteomeTools and MassIVE-KB in advancing ML applications in proteomics and discuss the implications of data set size on model performance, highlighting that larger data sets typically yield more accurate models. To address data scarcity, we explore algorithmic strategies such as self-supervised pretraining and multitask learning. Ultimately, we hope that this discussion can serve as a call to action for the proteomics community to collaborate on data standardization and collection efforts, which are crucial for the sustainable advancement and refinement of ML methodologies in the field.

Regulation of bacterial stringent response by an evolutionarily conserved ribosomal protein L11 methylation

Lysine and arginine methylation is an important regulator of enzyme activity and transcription in eukaryotes. However, little is known about this covalent modification in bacteria. In this work, we investigated the role of methylation in bacteria. By reanalyzing a large phyloproteomics data set from 48 bacterial strains representing six phyla, we found that almost a quarter of the bacterial proteome is methylated. Many of these methylated proteins are conserved across diverse bacterial lineages, including those involved in central carbon metabolism and translation. Among the proteins with the most conserved methylation sites is ribosomal protein L11 (bL11). bL11 methylation has been a mystery for five decades, as the deletion of its methyltransferase PrmA causes no cell growth defects. Comparative proteomics analysis combined with inorganic polyphosphate and guanosine tetra/pentaphosphate assays of the ΔprmA mutant in Escherichia coli revealed that bL11 methylation is important for stringent response signaling. In the stationary phase, we found that the ΔprmA mutant has impaired guanosine tetra/pentaphosphate production. This leads to a reduction in inorganic polyphosphate levels, accumulation of RNA and ribosomal proteins, and an abnormal polysome profile. Overall, our investigation demonstrates that the evolutionarily conserved bL11 methylation is important for stringent response signaling and ribosomal activity regulation and turnover.

IMPORTANCE

Protein methylation in bacteria was first identified over 60 years ago. Since then, its functional role has been identified for only a few proteins. To better understand the functional role of methylation in bacteria, we analyzed a large phyloproteomics data set encompassing 48 diverse bacteria. Our analysis revealed that ribosomal proteins are often methylated at conserved residues, suggesting that methylation of these sites may have a functional role in translation. Further analysis revealed that methylation of ribosomal protein L11 is important for stringent response signaling and ribosomal homeostasis.

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.

Exploring the dynamic landscape of immunopeptidomics: Unravelling posttranslational modifications and navigating bioinformatics terrain

Immunopeptidomics is becoming an increasingly important field of study. The capability to identify immunopeptides with pivotal roles in the human immune system is essential to shift the current curative medicine towards personalized medicine. Throughout the years, the field has matured, giving insight into the current pitfalls. Nowadays, it is commonly accepted that generalizing shotgun proteomics workflows is malpractice because immunopeptidomics faces numerous challenges. While many of these difficulties have been addressed, the road towards the ideal workflow remains complicated. Although the presence of Posttranslational modifications (PTMs) in the immunopeptidome has been demonstrated, their identification remains highly challenging despite their significance for immunotherapies. The large number of unpredictable modifications in the immunopeptidome plays a pivotal role in the functionality and these challenges. This review provides a comprehensive overview of the current advancements in immunopeptidomics. We delve into the challenges associated with identifying PTMs within the immunopeptidome, aiming to address the current state of the field.

Coupling Microdroplet-Based Sample Preparation, Multiplexed Isobaric Labeling, and Nanoflow Peptide Fractionation for Deep Proteome Profiling of the Tissue Microenvironment

There is increasing interest in developing in-depth proteomic approaches for mapping tissue heterogeneity in a cell-type-specific manner to better understand and predict the function of complex biological systems such as human organs. Existing spatially resolved proteomics technologies cannot provide deep proteome coverage due to limited sensitivity and poor sample recovery. Herein, we seamlessly combined laser capture microdissection with a low-volume sample processing technology that includes a microfluidic device named microPOTS (microdroplet processing in one pot for trace samples), multiplexed isobaric labeling, and a nanoflow peptide fractionation approach. The integrated workflow allowed us to maximize proteome coverage of laser-isolated tissue samples containing nanogram levels of proteins. We demonstrated that the deep spatial proteomics platform can quantify more than 5000 unique proteins from a small-sized human pancreatic tissue pixel (∼60,000 μm2) and differentiate unique protein abundance patterns in pancreas. Furthermore, the use of the microPOTS chip eliminated the requirement for advanced microfabrication capabilities and specialized nanoliter liquid handling equipment, making it more accessible to proteomic laboratories.

MS2Rescore 3.0 Is a Modular, Flexible, and User-Friendly Platform to Boost Peptide Identifications, as Showcased with MS Amanda 3.0

Rescoring of peptide-spectrum matches (PSMs) has emerged as a standard procedure for the analysis of tandem mass spectrometry data. This emphasizes the need for software maintenance and continuous improvement for such algorithms. We introduce MS2Rescore 3.0, a versatile, modular, and user-friendly platform designed to increase peptide identifications. Researchers can install MS2Rescore across various platforms with minimal effort and benefit from a graphical user interface, a modular Python API, and extensive documentation. To showcase this new version, we connected MS2Rescore 3.0 with MS Amanda 3.0, a new release of the well-established search engine, addressing previous limitations on automatic rescoring. Among new features, MS Amanda now contains additional output columns that can be used for rescoring. The full potential of rescoring is best revealed when applied on challenging data sets. We therefore evaluated the performance of these two tools on publicly available single-cell data sets, where the number of PSMs was substantially increased, thereby demonstrating that MS2Rescore offers a powerful solution to boost peptide identifications. MS2Rescore's modular design and user-friendly interface make data-driven rescoring easily accessible, even for inexperienced users. We therefore expect the MS2Rescore to be a valuable tool for the wider proteomics community. MS2Rescore is available at https://github.com/compomics/ms2rescore.

Mapping microhabitats of lignocellulose decomposition by a microbial consortium

The leaf-cutter ant fungal garden ecosystem is a naturally evolved model system for efficient plant biomass degradation. Degradation processes mediated by the symbiotic fungus Leucoagaricus gongylophorus are difficult to characterize due to dynamic metabolisms and spatial complexity of the system. Herein, we performed microscale imaging across 12-µm-thick adjacent sections of Atta cephalotes fungal gardens and applied a metabolome-informed proteome imaging approach to map lignin degradation. This approach combines two spatial multiomics mass spectrometry modalities that enabled us to visualize colocalized metabolites and proteins across and through the fungal garden. Spatially profiled metabolites revealed an accumulation of lignin-related products, outlining morphologically unique lignin microhabitats. Metaproteomic analyses of these microhabitats revealed carbohydrate-degrading enzymes, indicating a prominent fungal role in lignocellulose decomposition. Integration of metabolome-informed proteome imaging data provides a comprehensive view of underlying biological pathways to inform our understanding of metabolic fungal pathways in plant matter degradation within the micrometer-scale environment.

Alternative proteoforms and proteoform-dependent assemblies in humans and plants

The variability of proteins at the sequence level creates an enormous potential for proteome complexity. Exploring the depths and limits of this complexity is an ongoing goal in biology. Here, we systematically survey human and plant high-throughput bottom-up native proteomics data for protein truncation variants, where substantial regions of the full-length protein are missing from an observed protein product. In humans, Arabidopsis, and the green alga Chlamydomonas, approximately one percent of observed proteins show a short form, which we can assign by comparison to RNA isoforms as either likely deriving from transcript-directed processes or limited proteolysis. While some detected protein fragments align with known splice forms and protein cleavage events, multiple examples are previously undescribed, such as our observation of fibrocystin proteolysis and nuclear translocation in a green alga. We find that truncations occur almost entirely between structured protein domains, even when short forms are derived from transcript variants. Intriguingly, multiple endogenous protein truncations of phase-separating translational proteins resemble cleaved proteoforms produced by enteroviruses during infection. Some truncated proteins are also observed in both humans and plants, suggesting that they date to the last eukaryotic common ancestor. Finally, we describe novel proteoform-specific protein complexes, where the loss of a domain may accompany complex formation.

Stress-induced nuclear translocation of ONAC023 improves drought and heat tolerance through multiple processes in rice

Drought and heat are major abiotic stresses frequently coinciding to threaten rice production. Despite hundreds of stress-related genes being identified, only a few have been confirmed to confer resistance to multiple stresses in crops. Here we report ONAC023, a hub stress regulator that integrates the regulations of both drought and heat tolerance in rice. ONAC023 positively regulates drought and heat tolerance at both seedling and reproductive stages. Notably, the functioning of ONAC023 is obliterated without stress treatment and can be triggered by drought and heat stresses at two layers. The expression of ONAC023 is induced in response to stress stimuli. We show that overexpressed ONAC23 is translocated to the nucleus under stress and evidence from protoplasts suggests that the dephosphorylation of the remorin protein OSREM1.5 can promote this translocation. Under drought or heat stress, the nuclear ONAC023 can target and promote the expression of diverse genes, such as OsPIP2;7, PGL3, OsFKBP20-1b, and OsSF3B1, which are involved in various processes including water transport, reactive oxygen species homeostasis, and alternative splicing. These results manifest that ONAC023 is fine-tuned to positively regulate drought and heat tolerance through the integration of multiple stress-responsive processes. Our findings provide not only an underlying connection between drought and heat responses, but also a promising candidate for engineering multi-stress-resilient rice.

The Open Pediatric Cancer Project

Background

In 2019, the Open Pediatric Brain Tumor Atlas (OpenPBTA) was created as a global, collaborative open-science initiative to genomically characterize 1,074 pediatric brain tumors and 22 patient-derived cell lines. Here, we extend the OpenPBTA to create the Open Pediatric Cancer (OpenPedCan) Project, a harmonized open-source multi-omic dataset from 6,112 pediatric cancer patients with 7,096 tumor events across more than 100 histologies. Combined with RNA-Seq from the Genotype-Tissue Expression (GTEx) and The Cancer Genome Atlas (TCGA), OpenPedCan contains nearly 48,000 total biospecimens (24,002 tumor and 23,893 normal specimens).

Findings

We utilized Gabriella Miller Kids First (GMKF) workflows to harmonize WGS, WXS, RNA-seq, and Targeted Sequencing datasets to include somatic SNVs, InDels, CNVs, SVs, RNA expression, fusions, and splice variants. We integrated summarized CPTAC whole cell proteomics and phospho-proteomics data, miRNA-Seq data, and have developed a methylation array harmonization workflow to include m-values, beta-vales, and copy number calls. OpenPedCan contains reproducible, dockerized workflows in GitHub, CAVATICA, and Amazon Web Services (AWS) to deliver harmonized and processed data from over 60 scalable modules which can be leveraged both locally and on AWS. The processed data are released in a versioned manner and accessible through CAVATICA or AWS S3 download (from GitHub), and queryable through PedcBioPortal and the NCI's pediatric Molecular Targets Platform. Notably, we have expanded PBTA molecular subtyping to include methylation information to align with the WHO 2021 Central Nervous System Tumor classifications, allowing us to create research- grade integrated diagnoses for these tumors.

Conclusions

OpenPedCan data and its reproducible analysis module framework are openly available and can be utilized and/or adapted by researchers to accelerate discovery, validation, and clinical translation.

Multi-scale signaling and tumor evolution in high-grade gliomas

Although genomic anomalies in glioblastoma (GBM) have been well studied for over a decade, its 5-year survival rate remains lower than 5%. We seek to expand the molecular landscape of high-grade glioma, composed of IDH-wildtype GBM and IDH-mutant grade 4 astrocytoma, by integrating proteomic, metabolomic, lipidomic, and post-translational modifications (PTMs) with genomic and transcriptomic measurements to uncover multi-scale regulatory interactions governing tumor development and evolution. Applying 14 proteogenomic and metabolomic platforms to 228 tumors (212 GBM and 16 grade 4 IDH-mutant astrocytoma), including 28 at recurrence, plus 18 normal brain samples and 14 brain metastases as comparators, reveals heterogeneous upstream alterations converging on common downstream events at the proteomic and metabolomic levels and changes in protein-protein interactions and glycosylation site occupancy at recurrence. Recurrent genetic alterations and phosphorylation events on PTPN11 map to important regulatory domains in three dimensions, suggesting a central role for PTPN11 signaling across high-grade gliomas.

APIR: Aggregating Universal Proteomics Database Search Algorithms for Peptide Identification with FDR Control

Advances in mass spectrometry (MS) have enabled high-throughput analysis of proteomes in biological systems. The state-of-the-art MS data analysis relies on database search algorithms to quantify proteins by identifying peptide-spectrum matches (PSMs), which convert mass spectra to peptide sequences. Different database search algorithms use distinct search strategies and thus may identify unique PSMs. However, no existing approaches can aggregate all user-specified database search algorithms with a guaranteed increase in the number of identified peptides and a control on the false discovery rate (FDR). To fill in this gap, we proposed a statistical framework, Aggregation of Peptide Identification Results (APIR), that is universally compatible with all database search algorithms. Notably, under an FDR threshold, APIR is guaranteed to identify at least as many, if not more, peptides as individual database search algorithms do. Evaluation of APIR on a complex proteomics standard dataset showed that APIR outpowers individual database search algorithms and empirically controls the FDR. Real data studies showed that APIR can identify disease-related proteins and post-translational modifications missed by some individual database search algorithms. The APIR framework is easily extendable to aggregating discoveries made by multiple algorithms in other high-throughput biomedical data analysis, e.g., differential gene expression analysis on RNA sequencing data. The APIR R package is available at https://github.com/yiling0210/APIR.

Evaluation of a Reference-Free Collision Cross Section Calibration Strategy for Proteomics Using SLIM-Based High-Resolution Ion Mobility Spectrometry-Mass Spectrometry

Ion mobility spectrometry (IMS) is a gas-phase analytical technique that separates ions with different sizes and shapes and is compatible with mass spectrometry (MS) to provide an additional separation dimension. The rapid nature of the IMS separation combined with the high sensitivity of MS-based detection and the ability to derive structural information on analytes in the form of the property collision cross section (CCS) makes IMS particularly well-suited for characterizing complex samples in -omics applications. In such applications, the quality of CCS from IMS measurements is critical to confident annotation of the detected components in the complex -omics samples. However, most IMS instrumentation in mainstream use requires calibration to calculate CCS from measured arrival times, with the most notable exception being drift tube IMS measurements using multifield methods. The strategy for calibrating CCS values, particularly selection of appropriate calibrants, has important implications for CCS accuracy, reproducibility, and transferability between laboratories. The conventional approach to CCS calibration involves explicitly defining calibrants ahead of data acquisition and crucially relies upon availability of reference CCS values. In this work, we present a novel reference-free approach to CCS calibration which leverages trends among putatively identified features and computational CCS prediction to conduct calibrations post-data acquisition and without relying on explicitly defined calibrants. We demonstrated the utility of this reference-free CCS calibration strategy for proteomics application using high-resolution structures for lossless ion manipulations (SLIM)-based IMS-MS. We first validated the accuracy of CCS values using a set of synthetic peptides and then demonstrated using a complex peptide sample from cell lysate.

Rescoring Peptide Spectrum Matches: Boosting Proteomics Performance by Integrating Peptide Property Predictors Into Peptide Identification

Rescoring of peptide spectrum matches originating from database search engines enabled by peptide property predictors is exceeding the performance of peptide identification from traditional database search engines. In contrast to the peptide spectrum match scores calculated by traditional database search engines, rescoring peptide spectrum matches generates scores based on comparing observed and predicted peptide properties, such as fragment ion intensities and retention times. These newly generated scores enable a more efficient discrimination between correct and incorrect peptide spectrum matches. This approach was shown to lead to substantial improvements in the number of confidently identified peptides, facilitating the analysis of challenging datasets in various fields such as immunopeptidomics, metaproteomics, proteogenomics, and single-cell proteomics. In this review, we summarize the key elements leading up to the recent introduction of multiple data-driven rescoring pipelines. We provide an overview of relevant post-processing rescoring tools, introduce prominent data-driven rescoring pipelines for various applications, and highlight limitations, opportunities, and future perspectives of this approach and its impact on mass spectrometry-based proteomics.

The quantification of zebrafish ocular-associated proteins provides hints for sex-biased visual impairments and perception

Biochemical differences between sexes can also be seen in non-sexual organs and may affect organ functions and susceptibility to diseases. It has been shown that there are sex-biased visual perceptions and impairments. Abundance differences of eye proteins could provide explanations for some of these. Exploration of the ocular proteome was performed to find sex-based protein abundance differences in zebrafish Danio rerio. A label-free protein quantification workflow using high-resolution mass spectrometry was employed to find proteins with significant differences between the sexes. In total, 3740 unique master proteins were identified and quantified, and 49 proteins showed significant abundance differences between the eyes of male and female zebrafish. Those proteins belong to lipoproteins, immune system, blood coagulation, antioxidants, iron and heme-binding proteins, ion channels, pumps and exchangers, neuronal and photoreceptor proteins, and the cytoskeleton. An extensive literature review provided clues for the possible links between the sex-biased level of proteins and visual perception and impairments. In conclusion, sexual dimorphism at the protein level was discovered for the first time in the eye of zebrafish and should be accounted for in ophthalmological studies. Data are available via ProteomeXchange with identifier PXD033338.

Proteomes of plasmodium knowlesi early and late ring-stage parasites and infected host erythrocytes

The emerging malaria parasite Plasmodium knowlesi threatens the goal of worldwide malaria elimination due to its zoonotic spread in Southeast Asia. After brief ex-vivo culture we used 2D LC/MS/MS to examine the early and late ring stages of infected Macaca mulatta red blood cells harboring P. knowlesi. The M. mulatta clathrin heavy chain and T-cell and macrophage inhibitor ERMAP were overexpressed in the early ring stage; glutaredoxin 3 was overexpressed in the late ring stage; GO term differential enrichments included response to oxidative stress and the cortical cytoskeleton in the early ring stage. P. knowlesi clathrin heavy chain and 60S acidic ribosomal protein P2 were overexpressed in the late ring stage; GO term differential enrichments included vacuoles in the early ring stage, ribosomes and translation in the late ring stage, and Golgi- and COPI-coated vesicles, proteasomes, nucleosomes, vacuoles, ion-, peptide-, protein-, nucleocytoplasmic- and RNA-transport, antioxidant activity and glycolysis in both stages. SIGNIFICANCE: Due to its zoonotic spread, cases of the emerging human pathogen Plasmodium knowlesi in southeast Asia, and particularly in Malaysia, threaten regional and worldwide goals for malaria elimination. Infection by this parasite can be fatal to humans, and can be associated with significant morbidity. Due to zoonotic transmission from large macaque reservoirs that are untreatable by drugs, and outdoor biting mosquito vectors that negate use of preventive measures such as bed nets, its containment remains a challenge. Its biology remains incompletely understood. Thus we examine the expressed proteome of the early and late ex-vivo cultured ring stages, the first intraerythrocyte developmental stages after infection of host rhesus macaque erythrocytes. We used GO term enrichment strategies and differential protein expression to compare early and late ring stages. The early ring stage is characterized by the enrichment of P. knowlesi vacuoles, and overexpression of the M. mulatta clathrin heavy chain, important for clathrin-coated pits and vesicles, and clathrin-mediated endocytosis. The M. mulatta protein ERMAP was also overexpressed in the early ring stage, suggesting a potential role in early ring stage inhibition of T-cells and macrophages responding to P. knowlesi infection of reticulocytes. This could allow expansion of the host P. knowlesi cellular niche, allowing parasite adaptation to invasion of a wider age range of RBCs than the preferred young RBCs or reticulocytes, resulting in proliferation and increased pathogenesis in infected humans. Other GO terms differentially enriched in the early ring stage include the M. mulatta cortical cytoskeleton and response to oxidative stress. The late ring stage is characterized by overexpression of the P. knowlesi clathrin heavy chain. Combined with late ring stage GO term enrichment of Golgi-associated and coated vesicles, and enrichment of COPI-coated vesicles in both stages, this suggests the importance to P. knowlesi biology of clathrin-mediated endocytosis. P. knowlesi ribosomes and translation were also differentially enriched in the late ring stage. With expression of a variety of heat shock proteins, these results suggest production of folded parasite proteins is increasing by the late ring stage. M. mulatta endocytosis was differentially enriched in the late ring stage, as were clathrin-coated vesicles and endocytic vesicles. This suggests that M. mulatta clathrin-based endocytosis, perhaps in infected reticulocytes rather than mature RBC, may be an important process in the late ring stage. Additional ring stage biology from enriched GO terms includes M. mulatta proteasomes, protein folding and the chaperonin-containing T complex, actin and cortical actin cytoskeletons. P knowlesi biology also includes proteasomes, as well as nucleosomes, antioxidant activity, a variety of transport processes, glycolysis, vacuoles and protein folding. Mature RBCs have lost internal organelles, suggesting infection here may involve immature reticulocytes still retaining organelles. P. knowlesi parasite proteasomes and translational machinery may be ring stage drug targets for known selective inhibitors of these processes in other Plasmodium species. To our knowledge this is the first examination of more than one timepoint within the ring stage. Our results expand knowledge of both host and parasite proteins, pathways and organelles underlying P. knowlesi ring stage biology.

An algorithm for decoy-free false discovery rate estimation in XL-MS/MS proteomics

MOTIVATION

Cross-linking tandem mass spectrometry (XL-MS/MS) is an established analytical platform used to determine distance constraints between residues within a protein or from physically interacting proteins, thus improving our understanding of protein structure and function. To aid biological discovery with XL-MS/MS, it is essential that pairs of chemically linked peptides be accurately identified, a process that requires: (i) database search, that creates a ranked list of candidate peptide pairs for each experimental spectrum and (ii) false discovery rate (FDR) estimation, that determines the probability of a false match in a group of top-ranked peptide pairs with scores above a given threshold. Currently, the only available FDR estimation mechanism in XL-MS/MS is the target-decoy approach (TDA). However, despite its simplicity, TDA has both theoretical and practical limitations that impact the estimation accuracy and increase run time over potential decoy-free approaches (DFAs).

RESULTS

We introduce a novel decoy-free framework for FDR estimation in XL-MS/MS. Our approach relies on multi-sample mixtures of skew normal distributions, where the latent components correspond to the scores of correct peptide pairs (both peptides identified correctly), partially incorrect peptide pairs (one peptide identified correctly, the other incorrectly), and incorrect peptide pairs (both peptides identified incorrectly). To learn these components, we exploit the score distributions of first- and second-ranked peptide-spectrum matches for each experimental spectrum and subsequently estimate FDR using a novel expectation-maximization algorithm with constraints. We evaluate the method on ten datasets and provide evidence that the proposed DFA is theoretically sound and a viable alternative to TDA owing to its good performance in terms of accuracy, variance of estimation, and run time.

AVAILABILITY AND IMPLEMENTATION

https://github.com/shawn-peng/xlms.

SpecEncoder: deep metric learning for accurate peptide identification in proteomics

MOTIVATION

Tandem mass spectrometry (MS/MS) is a crucial technology for large-scale proteomic analysis. The protein database search or the spectral library search are commonly used for peptide identification from MS/MS spectra, which, however, may face challenges due to experimental variations between replicated spectra and similar fragmentation patterns among distinct peptides. To address this challenge, we present SpecEncoder, a deep metric learning approach to address these challenges by transforming MS/MS spectra into robust and sensitive embedding vectors in a latent space. The SpecEncoder model can also embed predicted MS/MS spectra of peptides, enabling a hybrid search approach that combines spectral library and protein database searches for peptide identification.

RESULTS

We evaluated SpecEncoder on three large human proteomics datasets, and the results showed a consistent improvement in peptide identification. For spectral library search, SpecEncoder identifies 1%-2% more unique peptides (and PSMs) than SpectraST. For protein database search, it identifies 6%-15% more unique peptides than MSGF+ enhanced by Percolator, Furthermore, SpecEncoder identified 6%-12% additional unique peptides when utilizing a combined library of experimental and predicted spectra. SpecEncoder can also identify more peptides when compared to deep-learning enhanced methods (MSFragger boosted by MSBooster). These results demonstrate SpecEncoder's potential to enhance peptide identification for proteomic data analyses.

AVAILABILITY AND IMPLEMENTATION

The source code and scripts for SpecEncoder and peptide identification are available on GitHub at https://github.com/lkytal/SpecEncoder. Contact: hatang@iu.edu.

Functional surface expression of immunoglobulin cleavage systems in a candidate Mycoplasma vaccine chassis

The Mycoplasma Immunoglobulin Binding/Protease (MIB-MIP) system is a candidate 'virulence factor present in multiple pathogenic species of the Mollicutes, including the fast-growing species Mycoplasma feriruminatoris. The MIB-MIP system cleaves the heavy chain of host immunoglobulins, hence affecting antigen-antibody interactions and potentially facilitating immune evasion. In this work, using -omics technologies and 5'RACE, we show that the four copies of the M. feriruminatoris MIB-MIP system have different expression levels and are transcribed as operons controlled by four different promoters. Individual MIB-MIP gene pairs of M. feriruminatoris and other Mollicutes were introduced in an engineered M. feriruminatoris strain devoid of MIB-MIP genes and were tested for their functionality using newly developed oriC-based plasmids. The two proteins are functionally expressed at the surface of M. feriruminatoris, which confirms the possibility to display large membrane-associated proteins in this bacterium. However, functional expression of heterologous MIB-MIP systems introduced in this engineered strain from phylogenetically distant porcine Mollicutes like Mesomycoplasma hyorhinis or Mesomycoplasma hyopneumoniae could not be achieved. Finally, since M. feriruminatoris is a candidate for biomedical applications such as drug delivery, we confirmed its safety in vivo in domestic goats, which are the closest livestock relatives to its native host the Alpine ibex.

A novel clinical metaproteomics workflow enables bioinformatic analysis of host-microbe dynamics in disease

Clinical metaproteomics has the potential to offer insights into the host-microbiome interactions underlying diseases. However, the field faces challenges in characterizing microbial proteins found in clinical samples, usually present at low abundance relative to the host proteins. As a solution, we have developed an integrated workflow coupling mass spectrometry-based analysis with customized bioinformatic identification, quantification, and prioritization of microbial proteins, enabling targeted assay development to investigate host-microbe dynamics in disease. The bioinformatics tools are implemented in the Galaxy ecosystem, offering the development and dissemination of complex bioinformatic workflows. The modular workflow integrates MetaNovo (to generate a reduced protein database), SearchGUI/PeptideShaker and MaxQuant [to generate peptide-spectral matches (PSMs) and quantification], PepQuery2 (to verify the quality of PSMs), Unipept (for taxonomic and functional annotation), and MSstatsTMT (for statistical analysis). We have utilized this workflow in diverse clinical samples, from the characterization of nasopharyngeal swab samples to bronchoalveolar lavage fluid. Here, we demonstrate its effectiveness via analysis of residual fluid from cervical swabs. The complete workflow, including training data and documentation, is available via the Galaxy Training Network, empowering non-expert researchers to utilize these powerful tools in their clinical studies.

IMPORTANCE

Clinical metaproteomics has immense potential to offer functional insights into the microbiome and its contributions to human disease. However, there are numerous challenges in the metaproteomic analysis of clinical samples, including handling of very large protein sequence databases for sensitive and accurate peptide and protein identification from mass spectrometry data, as well as taxonomic and functional annotation of quantified peptides and proteins to enable interpretation of results. To address these challenges, we have developed a novel clinical metaproteomics workflow that provides customized bioinformatic identification, verification, quantification, and taxonomic and functional annotation. This bioinformatic workflow is implemented in the Galaxy ecosystem and has been used to characterize diverse clinical sample types, such as nasopharyngeal swabs and bronchoalveolar lavage fluid. Here, we demonstrate its effectiveness and availability for use by the research community via analysis of residual fluid from cervical swabs.

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.

Untargeted, tandem mass spectrometry metaproteome of Columbia River sediments

Rivers are critical ecosystems that impact global biogeochemical cycles. Nonetheless, a mechanistic understanding of river microbial metabolisms and their influences on geochemistry is lacking. Here, we announce metaproteomes of river sediments that are paired with metagenomes and metabolites, enabling an understanding of the microbial underpinnings of river respiration.

Assessment of false discovery rate control in tandem mass spectrometry analysis using entrapment

A pressing statistical challenge in the field of mass spectrometry proteomics is how to assess whether a given software tool provides accurate error control. Each software tool for searching such data uses its own internally implemented methodology for reporting and controlling the error. Many of these software tools are closed source, with incompletely documented methodology, and the strategies for validating the error are inconsistent across tools. In this work, we identify three different methods for validating false discovery rate (FDR) control in use in the field, one of which is invalid, one of which can only provide a lower bound rather than an upper bound, and one of which is valid but under-powered. The result is that the field has a very poor understanding of how well we are doing with respect to FDR control, particularly for the analysis of data-independent acquisition (DIA) data. We therefore propose a new, more powerful method for evaluating FDR control in this setting, and we then employ that method, along with an existing lower bounding technique, to characterize a variety of popular search tools. We find that the search tools for analysis of data-dependent acquisition (DDA) data generally seem to control the FDR at the peptide level, whereas none of the DIA search tools consistently controls the FDR at the peptide level across all the datasets we investigated. Furthermore, this problem becomes much worse when the latter tools are evaluated at the protein level. These results may have significant implications for various downstream analyses, since proper FDR control has the potential to reduce noise in discovery lists and thereby boost statistical power.

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.

Multi-Omic Analysis of Esophageal Adenocarcinoma Uncovers Candidate Therapeutic Targets and Cancer-Selective Posttranscriptional Regulation

Efforts to address the poor prognosis associated with esophageal adenocarcinoma (EAC) have been hampered by a lack of biomarkers to identify early disease and therapeutic targets. Despite extensive efforts to understand the somatic mutations associated with EAC over the past decade, a gap remains in understanding how the atlas of genomic aberrations in this cancer impacts the proteome and which somatic variants are of importance for the disease phenotype. We performed a quantitative proteomic analysis of 23 EACs and matched adjacent normal esophageal and gastric tissues. We explored the correlation of transcript and protein abundance using tissue-matched RNA-seq and proteomic data from seven patients and further integrated these data with a cohort of EAC RNA-seq data (n = 264 patients), EAC whole-genome sequencing (n = 454 patients), and external published datasets. We quantified protein expression from 5879 genes in EAC and patient-matched normal tissues. Several biomarker candidates with EAC-selective expression were identified, including the transmembrane protein GPA33. We further verified the EAC-enriched expression of GPA33 in an external cohort of 115 patients and confirm this as an attractive diagnostic and therapeutic target. To further extend the insights gained from our proteomic data, an integrated analysis of protein and RNA expression in EAC and normal tissues revealed several genes with poorly correlated protein and RNA abundance, suggesting posttranscriptional regulation of protein expression. These outlier genes, including SLC25A30, TAOK2, and AGMAT, only rarely demonstrated somatic mutation, suggesting post-transcriptional drivers for this EAC-specific phenotype. AGMAT was demonstrated to be overexpressed at the protein level in EAC compared to adjacent normal tissues with an EAC-selective, post-transcriptional mechanism of regulation of protein abundance proposed. Integrated analysis of proteome, transcriptome, and genome in EAC has revealed several genes with tumor-selective, posttranscriptional regulation of protein expression, which may be an exploitable vulnerability.

Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection

BACKGROUND

Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors.

METHODS

We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers.

RESULTS

We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations in TP53, PLXNB3, Zinc finger genes, and fewer mutations in SDC2, CDH1, PIK3CA, SVIL, and PTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA-protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, including KIAA1522, DCK, FOXO3, MYO9B, ARID1A, EPRS, ZC3HAV1, and RBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors.

CONCLUSIONS

This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

DrugMap: A quantitative pan-cancer analysis of cysteine ligandability

Cysteine-focused chemical proteomic platforms have accelerated the clinical development of covalent inhibitors for a wide range of targets in cancer. However, how different oncogenic contexts influence cysteine targeting remains unknown. To address this question, we have developed "DrugMap," an atlas of cysteine ligandability compiled across 416 cancer cell lines. We unexpectedly find that cysteine ligandability varies across cancer cell lines, and we attribute this to differences in cellular redox states, protein conformational changes, and genetic mutations. Leveraging these findings, we identify actionable cysteines in NF-κB1 and SOX10 and develop corresponding covalent ligands that block the activity of these transcription factors. We demonstrate that the NF-κB1 probe blocks DNA binding, whereas the SOX10 ligand increases SOX10-SOX10 interactions and disrupts melanoma transcriptional signaling. Our findings reveal heterogeneity in cysteine ligandability across cancers, pinpoint cell-intrinsic features driving cysteine targeting, and illustrate the use of covalent probes to disrupt oncogenic transcription-factor activity.

Metabolomics and Proteomics in Prostate Cancer Research: Overview, Analytical Techniques, Data Analysis, and Recent Clinical Applications

Prostate cancer (PCa) is a significant global contributor to mortality, predominantly affecting males aged 65 and above. The field of omics has recently gained traction due to its capacity to provide profound insights into the biochemical mechanisms underlying conditions like prostate cancer. This involves the identification and quantification of low-molecular-weight metabolites and proteins acting as crucial biochemical signals for early detection, therapy assessment, and target identification. A spectrum of analytical methods is employed to discern and measure these molecules, revealing their altered biological pathways within diseased contexts. Metabolomics and proteomics generate refined data subjected to detailed statistical analysis through sophisticated software, yielding substantive insights. This review aims to underscore the major contributions of multi-omics to PCa research, covering its core principles, its role in tumor biology characterization, biomarker discovery, prognostic studies, various analytical technologies such as mass spectrometry and Nuclear Magnetic Resonance, data processing, and recent clinical applications made possible by an integrative "omics" approach. This approach seeks to address the challenges associated with current PCa treatments. Hence, our research endeavors to demonstrate the valuable applications of these potent tools in investigations, offering significant potential for understanding the complex biochemical environment of prostate cancer and advancing tailored therapeutic approaches for further development.

Time-resolved profiling of RNA binding proteins throughout the mRNA life cycle

mRNAs continually change their protein partners throughout their lifetimes, yet our understanding of mRNA-protein complex (mRNP) remodeling is limited by a lack of temporal data. Here, we present time-resolved mRNA interactome data by performing pulse metabolic labeling with photoactivatable ribonucleoside in human cells, UVA crosslinking, poly(A)+ RNA isolation, and mass spectrometry. This longitudinal approach allowed the quantification of over 700 RNA binding proteins (RBPs) across ten time points. Overall, the sequential order of mRNA binding aligns well with known functions, subcellular locations, and molecular interactions. However, we also observed RBPs with unexpected dynamics: the transcription-export (TREX) complex recruited posttranscriptionally after nuclear export factor 1 (NXF1) binding, challenging the current view of transcription-coupled mRNA export, and stress granule proteins prevalent in aged mRNPs, indicating roles in late stages of the mRNA life cycle. To systematically identify mRBPs with unknown functions, we employed machine learning to compare mRNA binding dynamics with Gene Ontology (GO) annotations. Our data can be explored at chronology.rna.snu.ac.kr.

E-cadherin interacts with EGFR resulting in hyper-activation of ERK in multiple models of breast cancer

The loss of intercellular adhesion molecule E-cadherin is a hallmark of the epithelial-mesenchymal transition (EMT), during which tumor cells transition into an invasive phenotype. Accordingly, E-cadherin has long been considered a tumor suppressor gene; however, E-cadherin expression is paradoxically correlated with breast cancer survival rates. Using novel multi-compartment organoids and multiple in vivo models, we show that E-cadherin promotes a hyper-proliferative phenotype in breast cancer cells via interaction with the transmembrane receptor EGFR. The E-cad and EGFR interaction results in activation of the MEK/ERK signaling pathway, leading to a significant increase in proliferation via activation of transcription factors, including c-Fos. Pharmacological inhibition of MEK activity in E-cadherin positive breast cancer significantly decreases both tumor growth and macro-metastasis in vivo. This work provides evidence for a novel role of E-cadherin in breast tumor progression and identifies a new target to treat hyper-proliferative E-cadherin-positive breast tumors, thus providing the foundation to utilize E-cadherin as a biomarker for specific therapeutic success.

Target-decoy false discovery rate estimation using Crema

Assigning statistical confidence estimates to discoveries produced by a tandem mass spectrometry proteomics experiment is critical to enabling principled interpretation of the results and assessing the cost/benefit ratio of experimental follow-up. The most common technique for computing such estimates is to use target-decoy competition (TDC), in which observed spectra are searched against a database of real (target) peptides and a database of shuffled or reversed (decoy) peptides. TDC procedures for estimating the false discovery rate (FDR) at a given score threshold have been developed for application at the level of spectra, peptides, or proteins. Although these techniques are relatively straightforward to implement, it is common in the literature to skip over the implementation details or even to make mistakes in how the TDC procedures are applied in practice. Here we present Crema, an open-source Python tool that implements several TDC methods of spectrum-, peptide- and protein-level FDR estimation. Crema is compatible with a variety of existing database search tools and provides a straightforward way to obtain robust FDR estimates.

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.

Machine learning-based peptide-spectrum match rescoring opens up the immunopeptidome

Immunopeptidomics is a key technology in the discovery of targets for immunotherapy and vaccine development. However, identifying immunopeptides remains challenging due to their non-tryptic nature, which results in distinct spectral characteristics. Moreover, the absence of strict digestion rules leads to extensive search spaces, further amplified by the incorporation of somatic mutations, pathogen genomes, unannotated open reading frames, and post-translational modifications. This inflation in search space leads to an increase in random high-scoring matches, resulting in fewer identifications at a given false discovery rate. Peptide-spectrum match rescoring has emerged as a machine learning-based solution to address challenges in mass spectrometry-based immunopeptidomics data analysis. It involves post-processing unfiltered spectrum annotations to better distinguish between correct and incorrect peptide-spectrum matches. Recently, features based on predicted peptidoform properties, including fragment ion intensities, retention time, and collisional cross section, have been used to improve the accuracy and sensitivity of immunopeptide identification. In this review, we describe the diverse bioinformatics pipelines that are currently available for peptide-spectrum match rescoring and discuss how they can be used for the analysis of immunopeptidomics data. Finally, we provide insights into current and future machine learning solutions to boost immunopeptide identification.

Challenges in the discovery of tumor-specific alternative splicing-derived cell-surface antigens in glioma

Despite advancements in cancer immunotherapy, solid tumors remain formidable challenges. In glioma, profound inter- and intra-tumoral heterogeneity of antigen landscape hampers therapeutic development. Therefore, it is critical to consider alternative sources to expand the repertoire of targetable (neo-)antigens and improve therapeutic outcomes. Accumulating evidence suggests that tumor-specific alternative splicing (AS) could be an untapped reservoir of antigens. In this study, we investigated tumor-specific AS events in glioma, focusing on those predicted to generate major histocompatibility complex (MHC)-presentation-independent, cell-surface antigens that could be targeted by antibodies and chimeric antigen receptor-T cells. We systematically analyzed bulk RNA-sequencing datasets comparing 429 tumor samples (from The Cancer Genome Atlas) and 9166 normal tissue samples (from the Genotype-Tissue Expression project), and identified 13 AS events in 7 genes predicted to be expressed in more than 10% of the patients, including PTPRZ1 and BCAN, which were corroborated by an external RNA-sequencing dataset. Subsequently, we validated our predictions and elucidated the complexity of the isoforms using full-length transcript amplicon sequencing on patient-derived glioblastoma cells. However, analyses of the RNA-sequencing datasets of spatially mapped and longitudinally collected clinical tumor samples unveiled remarkable spatiotemporal heterogeneity of the candidate AS events. Furthermore, proteomics analysis did not reveal any peptide spectra matching the putative antigens. Our investigation illustrated the diverse characteristics of the tumor-specific AS events and the challenges of antigen exploration due to their notable spatiotemporal heterogeneity and elusive nature at the protein levels. Redirecting future efforts toward intracellular, MHC-presented antigens could offer a more viable avenue.

A Comprehensive Understanding of Post-Translational Modification of Sox2 via Acetylation and O-GlcNAcylation in Colorectal Cancer

Aberrant expression of the pluripotency-associated transcription factor Sox2 is associated with poor prognosis in colorectal cancer (CRC). We investigated the regulatory roles of major post-translational modifications in Sox2 using two CRC cell lines, SW480 and SW620, derived from the same patient but with low and high Sox2 expression, respectively. Acetylation of K75 in the Sox2 nuclear export signal was relatively increased in SW480 cells and promotes Sox2 nucleocytoplasmic shuttling and proteasomal degradation of Sox2. LC-MS-based proteomics analysis identified HDAC4 and p300 as binding partners involved in the acetylation-mediated control of Sox2 expression in the nucleus. Sox2 K75 acetylation is mediated by the acetyltransferase activity of CBP/p300 and ACSS3. In SW620 cells, HDAC4 deacetylates K75 and is regulated by miR29a. O-GlcNAcylation on S246, in addition to K75 acetylation, also regulates Sox2 stability. These findings provide insights into the regulation of Sox2 through multiple post-translational modifications and pathways in CRC.

Exact p-value calculation for XCorr scoring of high-resolution MS/MS data

Exact p-value (XPV)-based methods for dot product-like score functions-such as the XCorr score implemented in Tide, SEQUEST, Comet or shared peak count-based scoring in MSGF+ and ASPV-provide a fairly good calibration for peptide-spectrum-match (PSM) scoring in database searching-based MS/MS spectrum data identification. Unfortunately, standard XPV methods, in practice, cannot handle high-resolution fragmentation data produced by state-of-the-art mass spectrometers because having smaller bins increases the number of fragment matches that are assigned to incorrect bins and scored improperly. In this article, we present an extension of the XPV method, called the high-resolution exact p-value (HR-XPV) method, which can be used to calibrate PSM scores of high-resolution MS/MS spectra obtained with dot product-like scoring such as the XCorr. The HR-XPV carries remainder masses throughout the fragmentation, allowing them to greatly increase the number of fragments that are properly assigned to the correct bin and, thus, taking advantage of high-resolution data. Using four mass spectrometry data sets, our experimental results demonstrate that HR-XPV produces well-calibrated scores, which in turn results in more trusted spectrum annotations at any false discovery rate level.

Proteome changes in larval zebrafish (Danio rerio) and fathead minnow (Pimephales promelas) exposed to (±) anatoxin-a

Anatoxin-a and its analogues are potent neurotoxins produced by several genera of cyanobacteria. Due in part to its high toxicity and potential presence in drinking water, these toxins pose threats to public health, companion animals and the environment. It primarily exerts toxicity as a cholinergic agonist, with high affinity at neuromuscular junctions, but molecular mechanisms by which it elicits toxicological responses are not fully understood. To advance understanding of this cyanobacteria, proteomic characterization (DIA shotgun proteomics) of two common fish models (zebrafish and fathead minnow) was performed following  (±) anatoxin-a exposure. Specifically, proteome changes were identified and quantified in larval fish exposed for 96 h (0.01-3 mg/L (±) anatoxin-a and caffeine (a methodological positive control) with environmentally relevant treatment levels examined based on environmental exposure distributions of surface water data. Proteomic concentration - response relationships revealed 48 and 29 proteins with concentration - response relationships curves for zebrafish and fathead minnow, respectively. In contrast, the highest number of differentially expressed proteins (DEPs) varied between zebrafish (n = 145) and fathead minnow (n = 300), with only fatheads displaying DEPs at all treatment levels. For both species, genes associated with reproduction were significantly downregulated, with pathways analysis that broadly clustered genes into groups associated with DNA repair mechanisms. Importantly, significant differences in proteome response between the species was also observed, consistent with prior observations of differences in response using both behavioral assays and gene expression, adding further support to model specific differences in organismal sensitivity and/or response. When DEPs were read across from humans to zebrafish, disease ontology enrichment identified diseases associated with cognition and muscle weakness consistent with the prior literature. Our observations highlight limited knowledge of how (±) anatoxin-a, a commonly used synthetic racemate surrogate, elicits responses at a molecular level and advances its toxicological understanding.

Regulation of urea cycle by reversible high-stoichiometry lysine succinylation

The post-translational modification lysine succinylation is implicated in the regulation of various metabolic pathways. However, its biological relevance remains uncertain due to methodological difficulties in determining high-impact succinylation sites. Here, using stable isotope labelling and data-independent acquisition mass spectrometry, we quantified lysine succinylation stoichiometries in mouse livers. Despite the low overall stoichiometry of lysine succinylation, several high-stoichiometry sites were identified, especially upon deletion of the desuccinylase SIRT5. In particular, multiple high-stoichiometry lysine sites identified in argininosuccinate synthase (ASS1), a key enzyme in the urea cycle, are regulated by SIRT5. Mutation of the high-stoichiometry lysine in ASS1 to succinyl-mimetic glutamic acid significantly decreased its enzymatic activity. Metabolomics profiling confirms that SIRT5 deficiency decreases urea cycle activity in liver. Importantly, SIRT5 deficiency compromises ammonia tolerance, which can be reversed by the overexpression of wild-type, but not succinyl-mimetic, ASS1. Therefore, lysine succinylation is functionally important in ammonia metabolism.

CloudProteoAnalyzer: scalable processing of big data from proteomics using cloud computing

Summary

Shotgun proteomics is widely used in many system biology studies to determine the global protein expression profiles of tissues, cultures, and microbiomes. Many non-distributed computer algorithms have been developed for users to process proteomics data on their local computers. However, the amount of data acquired in a typical proteomics study has grown rapidly in recent years, owing to the increasing throughput of mass spectrometry and the expanding scale of study designs. This presents a big data challenge for researchers to process proteomics data in a timely manner. To overcome this challenge, we developed a cloud-based parallel computing application to offer end-to-end proteomics data analysis software as a service (SaaS). A web interface was provided to users to upload mass spectrometry-based proteomics data, configure parameters, submit jobs, and monitor job status. The data processing was distributed across multiple nodes in a supercomputer to achieve scalability for large datasets. Our study demonstrated SaaS for proteomics as a viable solution for the community to scale up the data processing using cloud computing.

Availability and implementation

This application is available online at https://sipros.oscer.ou.edu/ or https://sipros.unt.edu for free use. The source code is available at https://github.com/Biocomputing-Research-Group/CloudProteoAnalyzer under the GPL version 3.0 license.

Delineating functional and molecular impact of ex vivo sample handling in precision medicine

Consistent handling of samples is crucial for achieving reproducible molecular and functional testing results in translational research. Here, we used 229 acute myeloid leukemia (AML) patient samples to assess the impact of sample handling on high-throughput functional drug testing, mass spectrometry-based proteomics, and flow cytometry. Our data revealed novel and previously described changes in cell phenotype and drug response dependent on sample biobanking. Specifically, myeloid cells with a CD117 (c-KIT) positive phenotype decreased after biobanking, potentially distorting cell population representations and affecting drugs targeting these cells. Additionally, highly granular AML cell numbers decreased after freezing. Secondly, protein expression levels, as well as sensitivity to drugs targeting cell proliferation, metabolism, tyrosine kinases (e.g., JAK, KIT, FLT3), and BH3 mimetics were notably affected by biobanking. Moreover, drug response profiles of paired fresh and frozen samples showed that freezing samples can lead to systematic errors in drug sensitivity scores. While a high correlation between fresh and frozen for the entire drug library was observed, freezing cells had a considerable impact at an individual level, which could influence outcomes in translational studies. Our study highlights conditions where standardization is needed to improve reproducibility, and where validation of data generated from biobanked cohorts may be particularly important.

Identification of structural and regulatory cell-shape determinants in Haloferax volcanii

Archaea play indispensable roles in global biogeochemical cycles, yet many crucial cellular processes, including cell-shape determination, are poorly understood. Haloferax volcanii, a model haloarchaeon, forms rods and disks, depending on growth conditions. Here, we used a combination of iterative proteomics, genetics, and live-cell imaging to identify mutants that only form rods or disks. We compared the proteomes of the mutants with wild-type cells across growth phases, thereby distinguishing between protein abundance changes specific to cell shape and those related to growth phases. The results identified a diverse set of proteins, including predicted transporters, transducers, signaling components, and transcriptional regulators, as important for cell-shape determination. Through phenotypic characterization of deletion strains, we established that rod-determining factor A (RdfA) and disk-determining factor A (DdfA) are required for the formation of rods and disks, respectively. We also identified structural proteins, including an actin homolog that plays a role in disk-shape morphogenesis, which we named volactin. Using live-cell imaging, we determined volactin's cellular localization and showed its dynamic polymerization and depolymerization. Our results provide insights into archaeal cell-shape determination, with possible implications for understanding the evolution of cell morphology regulation across domains.

A Composite Biomarker Signature of Type 1 Diabetes Risk Identified via Augmentation of Parallel Multi-Omics Data from a Small Cohort

Background

Biomarkers of early pathogenesis of type 1 diabetes (T1D) are crucial to enable effective prevention measures in at-risk populations before significant damage occurs to their insulin producing beta-cell mass. We recently introduced the concept of integrated parallel multi-omics and employed a novel data augmentation approach which identified promising candidate biomarkers from a small cohort of high-risk T1D subjects. We now validate selected biomarkers to generate a potential composite signature of T1D risk.

Methods

Twelve candidate biomarkers, which were identified in the augmented data and selected based on their fold-change relative to healthy controls and cross-reference to proteomics data previously obtained in the expansive TEDDY and DAISY cohorts, were measured in the original samples by ELISA.

Results

All 12 biomarkers had established connections with lipid/lipoprotein metabolism, immune function, inflammation, and diabetes, but only 7 were found to be markedly changed in the high-risk subjects compared to the healthy controls: ApoC1 and PON1 were reduced while CETP, CD36, FGFR1, IGHM, PCSK9, SOD1, and VCAM1 were elevated.

Conclusions

Results further highlight the promise of our data augmentation approach in unmasking important patterns and pathologically significant features in parallel multi-omics datasets obtained from small sample cohorts to facilitate the identification of promising candidate T1D biomarkers for downstream validation. They also support the potential utility of a composite biomarker signature of T1D risk characterized by the changes in the above markers.

UniSpec: Deep Learning for Predicting the Full Range of Peptide Fragment Ion Series to Enhance the Proteomics Data Analysis Workflow

We present UniSpec, an attention-driven deep neural network designed to predict comprehensive collision-induced fragmentation spectra, thereby improving peptide identification in shotgun proteomics. Utilizing a training data set of 1.8 million unique high-quality tandem mass spectra (MS2) from 0.8 million unique peptide ions, UniSpec learned with a peptide fragmentation dictionary encompassing 7919 fragment peaks. Among these, 5712 are neutral loss peaks, with 2310 corresponding to modification-specific neutral losses. Remarkably, UniSpec can predict 73%-77% of fragment intensities based on our NIST reference library spectra, a significant leap from the 35%-45% coverage of only b and y ions. Comparative studies with Prosit elucidate that while both models are strong at predicting their respective fragment ion series, UniSpec particularly shines in generating more complex MS2 spectra with diverse ion annotations. The integration of UniSpec's predictions into shotgun proteomics data analysis boosts the identification rate of tryptic peptides by 48% at a 1% false discovery rate (FDR) and 60% at a more confident 0.1% FDR. Using UniSpec's predicted in-silico spectral library, the search results closely matched those from search engines and experimental spectral libraries used in peptide identification, highlighting its potential as a stand-alone identification tool. The source code and Python scripts are available on GitHub (https://github.com/usnistgov/UniSpec) and Zenodo (https://zenodo.org/records/10452792), and all data sets and analysis results generated in this work were deposited in Zenodo (https://zenodo.org/records/10052268).