PeptideShaker enables reanalysis of MS-derived proteomics data sets

Optimization of protein identifications through the use of different chromatographic approaches and bioinformatic pipelines

RATIONALE

Selection of proteomic workflows for a given project can be a daunting task. This research provides a guide outlining the impact on protein identification of different steps such as chromatographic separation, data acquisition strategies, and bioinformatic pipelines. The data presented here will help experts and nonexpert proteomic users to increase proteome coverage and peptide identification.

METHODS

HeLa protein digests were analyzed through different C18 chromatographic columns (15 and 50 cm in length), using top 12 data-dependent acquisition (DDA), top 20 DDA, and data-independent acquisition (DIA) with a nanospray source in positive mode in a Thermo Q Exactive instrument. The raw data were analyzed using different search engines, rescoring approaches, and multi-engine searches. The results were analyzed in the context of peptide and protein identifications, precursor properties, and computation requirements to understand the differences between methods.

RESULTS

Our results showed that higher column lengths and top N DDA approaches were able to significantly increase protein identifications. The use of multiple search engines yielded limited gains, whereas the use of rescoring methods clearly outperformed other strategies. Finally, DIA approaches, although successful at generating new identifications, had a limited performance influenced by the previous collection of DDA data, which could prohibitively increase instrument time. Nonetheless, the use of library-free methods showed promising results.

CONCLUSIONS

Our results highlight the impact of different experimental approaches on proteome coverage. Changes in chromatographic columns, data acquisition, or bioinformatic analysis can significantly increase the number of protein identifications (>400%). Thus, this research provides a reference upon which to build a successful proteomic workflow with different considerations at every step.

Discovering Novel Proteoforms Using Proteogenomic Workflows Within the Galaxy Bioinformatics Platform

Proteogenomics is a growing "multi-omics" research area that combines mass spectrometry-based proteomics and high-throughput nucleotide sequencing technologies. Proteogenomics has helped in genomic annotation for organisms whose complete genome sequences became available by using high-throughput DNA sequencing technologies. Apart from genome annotation, this multi-omics approach has also helped researchers confirm expression of variant proteins belonging to unique proteoforms that could have resulted from single-nucleotide polymorphism (SNP), insertion and deletions (Indels), splice isoforms, or other genome or transcriptome variations.A proteogenomic study depends on a multistep informatics workflow, requiring different software at each step. These integrated steps include creating an appropriate protein sequence database, matching spectral data against these sequences, and finally identifying peptide sequences corresponding to novel proteoforms followed by variant classification and functional analysis. The disparate software required for a proteogenomic study is difficult for most researchers to access and use, especially those lacking computational expertise. Furthermore, using them disjointedly can be error-prone as it requires setting up individual parameters for each software. Consequently, reproducibility suffers. Managing output files from each software is an additional challenge. One solution for these challenges in proteogenomics is the open-source Web-based computational platform Galaxy. Its capability to create and manage workflows comprised of disparate software while recording and saving all important parameters promotes both usability and reproducibility. Here, we describe a workflow that can perform proteogenomic analysis on a Galaxy-based platform. This Galaxy workflow facilitates matching of spectral data with a customized protein sequence database, identifying novel protein variants, assessing quality of results, and classifying variants along with visualization against the genome.

Small Molecule Assembly Agonist Alters the Dynamics of Hepatitis B Virus Core Protein Dimer and Capsid

Chronic hepatitis B virus (HBV) poses a significant public health burden worldwide, encouraging the search for curative antivirals. One approach is capsid assembly modulators (CAMs), which are assembly agonists. CAMs lead to empty and defective capsids, inhibiting the formation of new viruses, and can also lead to defects in the release of the viral genome, inhibiting new infections. In this study, we employed hydrogen-deuterium exchange mass spectrometry (HDX-MS) to assess the impact of one such CAM, HAP18, on HBV dimers, capsids composed of 120 (or 90) capsid protein dimers, and cross-linked capsids (xl-capsids). HDX analysis revealed hydrogen bonding networks within and between the dimers. HAP18 disrupted the hydrogen bonding network of dimers, demonstrating a previously unappreciated impact on the dimer structure. Conversely, HAP18 stabilized both unmodified and cross-linked capsids. Intriguingly, cross-linking the capsid, which was accomplished by forming disulfides between an engineered C-terminal cysteine, increased the overall rate of HDX. Moreover, HAP18 binding induced conformational changes beyond the binding sites. Our findings provide evidence for allosteric communication within and between capsid protein dimers. These results show that CAMs are capable of harnessing this allosteric network to modulate the dimer and capsid dynamics.

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.

quantms: a cloud-based pipeline for quantitative proteomics enables the reanalysis of public proteomics data

The volume of public proteomics data is rapidly increasing, causing a computational challenge for large-scale reanalysis. Here, we introduce quantms ( https://quant,ms.org/ ), an open-source cloud-based pipeline for massively parallel proteomics data analysis. We used quantms to reanalyze 83 public ProteomeXchange datasets, comprising 29,354 instrument files from 13,132 human samples, to quantify 16,599 proteins based on 1.03 million unique peptides. quantms is based on standard file formats improving the reproducibility, submission and dissemination of the data to ProteomeXchange.

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.

Giant polyketide synthase enzymes in the biosynthesis of giant marine polyether toxins

Prymnesium parvum are harmful haptophyte algae that cause massive environmental fish kills. Their polyketide polyether toxins, the prymnesins, are among the largest nonpolymeric compounds in nature and have biosynthetic origins that have remained enigmatic for more than 40 years. In this work, we report the "PKZILLAs," massive P. parvum polyketide synthase (PKS) genes that have evaded previous detection. PKZILLA-1 and -2 encode giant protein products of 4.7 and 3.2 megadaltons that have 140 and 99 enzyme domains. Their predicted polyene product matches the proposed pre-prymnesin precursor of the 90-carbon-backbone A-type prymnesins. We further characterize the variant PKZILLA-B1, which is responsible for the shorter B-type analog prymnesin-B1, from P. parvum RCC3426 and thus establish a general model of haptophyte polyether biosynthetic logic. This work expands expectations of genetic and enzymatic size limits in biology.

Exploring Ion Mobility Mass Spectrometry Data File Conversions to Leverage Existing Tools and Enable New Workflows

Ion mobility (IM) is often combined with LC-MS experiments to provide an additional dimension of separation for complex sample analysis. While highly complex samples are better characterized by the full dimensionality of LC-IM-MS experiments to uncover new information, downstream data analysis workflows are often not equipped to properly mine the additional IM dimension. For many samples the data acquisition benefits of including IM separations are all that is necessary to uncover sample information and the full dimensionality of the data is not required for data analysis. Postacquisition reduction and adaptation of the dimensions of LC-IM-MS and IM-MS experiments into an LC-MS format opens the possibility to use a plethora of existing software tools. In this work, we developed data file conversion tools to reduce the complexity of IM data analysis. Three data file transformations are introduced in the PNNL PreProcessor software: (1) mapping the IM axis to the LC axis for IM-MS data, (2) converting the drift time vs m/z space to CCS/z vs m/z space, and (3) transforming All Ions IM/MS mobility aligned fragmentation data to a standard LC-MS DDA data file format. These new data file conversions are demonstrated with corresponding lipidomics and proteomics workflows that leverage existing LC-MS data analysis software to highlight the benefits of the data transformations.

Soil Metaproteomics for Microbial Community Profiling: Methodologies and Challenges

Soil represents a complex and dynamic ecosystem, hosting a myriad of microorganisms that coexist and play vital roles in nutrient cycling and organic matter transformation. Among these microorganisms, bacteria and fungi are key members of the microbial community, profoundly influencing the fate of nitrogen, sulfur, and carbon in terrestrial environments. Understanding the intricacies of soil ecosystems and the biological processes orchestrated by microbial communities necessitates a deep dive into their composition and metabolic activities. The advent of next-generation sequencing and 'omics' techniques, such as metagenomics and metaproteomics, has revolutionized our understanding of microbial ecology and the functional dynamics of soil microbial communities. Metagenomics enables the identification of microbial community composition in soil, while metaproteomics sheds light on the current biological functions performed by these communities. However, metaproteomics presents several challenges, both technical and computational. Factors such as the presence of humic acids and variations in extraction methods can influence protein yield, while the absence of high-resolution mass spectrometry and comprehensive protein databases limits the depth of protein identification. Notwithstanding these limitations, metaproteomics remains a potent tool for unraveling the intricate biological processes and functions of soil microbial communities. In this review, we delve into the methodologies and challenges of metaproteomics in soil research, covering aspects such as protein extraction, identification, and bioinformatics analysis. Furthermore, we explore the applications of metaproteomics in soil bioremediation, highlighting its potential in addressing environmental challenges.

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.

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.

Critical steps in an environmental metaproteomics workflow

Environmental metaproteomics is a rapidly advancing field that provides insights into the structure, dynamics, and metabolic activity of microbial communities. As the field is still maturing, it lacks consistent workflows, making it challenging for non-expert researchers to navigate. This review aims to introduce the workflow of environmental metaproteomics. It outlines the standard practices for sample collection, processing, and analysis, and offers strategies to overcome the unique challenges presented by common environmental matrices such as soil, freshwater, marine environments, biofilms, sludge, and symbionts. The review also highlights the bottlenecks in data analysis that are specific to metaproteomics samples and provides suggestions for researchers to obtain high-quality datasets. It includes recent benchmarking studies and descriptions of software packages specifically built for metaproteomics analysis. The article is written without assuming the reader's familiarity with single-organism proteomic workflows, making it accessible to those new to proteomics or mass spectrometry in general. This primer for environmental metaproteomics aims to improve accessibility to this exciting technology and empower researchers to tackle challenging and ambitious research questions. While it is primarily a resource for those new to the field, it should also be useful for established researchers looking to streamline or troubleshoot their metaproteomics experiments.

An adaptive biomolecular condensation response is conserved across environmentally divergent species

Cells must sense and respond to sudden maladaptive environmental changes-stresses-to survive and thrive. Across eukaryotes, stresses such as heat shock trigger conserved responses: growth arrest, a specific transcriptional response, and biomolecular condensation of protein and mRNA into structures known as stress granules under severe stress. The composition, formation mechanism, adaptive significance, and even evolutionary conservation of these condensed structures remain enigmatic. Here we provide a remarkable view into stress-triggered condensation, its evolutionary conservation and tuning, and its integration into other well-studied aspects of the stress response. Using three morphologically near-identical budding yeast species adapted to different thermal environments and diverged by up to 100 million years, we show that proteome-scale biomolecular condensation is tuned to species-specific thermal niches, closely tracking corresponding growth and transcriptional responses. In each species, poly(A)-binding protein-a core marker of stress granules-condenses in isolation at species-specific temperatures, with conserved molecular features and conformational changes modulating condensation. From the ecological to the molecular scale, our results reveal previously unappreciated levels of evolutionary selection in the eukaryotic stress response, while establishing a rich, tractable system for further inquiry.

Venom diversity in Naja mossambica: Insights from proteomic and immunochemical analyses reveal intraspecific differences

BACKGROUND

Intraspecific variations in snake venom composition have been extensively documented, contributing to the diverse clinical effects observed in envenomed patients. Understanding these variations is essential for developing effective snakebite management strategies and targeted antivenom therapies. We aimed to comprehensively investigate venoms from three distinct populations of N. mossambica from Eswatini, Limpopo, and KwaZulu-Natal regions in Africa in terms of their protein composition and reactivity with three commercial antivenoms (SAIMR polyvalent, EchiTAb+ICP, and Antivipmyn Africa).

METHODOLOGY/PRINCIPAL FINDINGS

Naja mossambica venoms from Eswatini region exhibited the highest content of neurotoxic proteins, constituting 20.70% of all venom proteins, compared to Limpopo (13.91%) and KwaZulu-Natal (12.80%), and was characterized by the highest diversity of neurotoxic proteins, including neurotoxic 3FTxs, Kunitz-type inhibitors, vespryns, and mamba intestinal toxin 1. KwaZulu-Natal population exhibited considerably lower cytotoxic 3FTx, higher PLA2 content, and significant diversity in low-abundant proteins. Conversely, Limpopo venoms demonstrated the least diversity as demonstrated by electrophoretic and mass spectrometry analyses. Immunochemical assessments unveiled differences in venom-antivenom reactivity, particularly concerning low-abundance proteins. EchiTAb+ICP antivenom demonstrated superior reactivity in serial dilution ELISA assays compared to SAIMR polyvalent.

CONCLUSIONS/SIGNIFICANCE

Our findings reveal a substantial presence of neurotoxic proteins in N. mossambica venoms, challenging previous understandings of their composition. Additionally, the detection of numerous peptides aligning to uncharacterized proteins or proteins with unknown functions underscores a critical issue with existing venom protein databases, emphasizing the substantial gaps in our knowledge of snake venom protein components. This underscores the need for enhanced research in this domain. Moreover, our in vitro immunological assays suggest EchiTAb+ICP's potential as an alternative to SAIMR antivenom, requiring confirmation through prospective in vivo neutralization studies.

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.

Systematic profiling of ale yeast protein dynamics across fermentation and repitching

Studying the genetic and molecular characteristics of brewing yeast strains is crucial for understanding their domestication history and adaptations accumulated over time in fermentation environments, and for guiding optimizations to the brewing process itself. Saccharomyces cerevisiae (brewing yeast) is among the most profiled organisms on the planet, yet the temporal molecular changes that underlie industrial fermentation and beer brewing remain understudied. Here, we characterized the genomic makeup of a Saccharomyces cerevisiae ale yeast widely used in the production of Hefeweizen beers, and applied shotgun mass spectrometry to systematically measure the proteomic changes throughout 2 fermentation cycles which were separated by 14 rounds of serial repitching. The resulting brewing yeast proteomics resource includes 64,740 protein abundance measurements. We found that this strain possesses typical genetic characteristics of Saccharomyces cerevisiae ale strains and displayed progressive shifts in molecular processes during fermentation based on protein abundance changes. We observed protein abundance differences between early fermentation batches compared to those separated by 14 rounds of serial repitching. The observed abundance differences occurred mainly in proteins involved in the metabolism of ergosterol and isobutyraldehyde. Our systematic profiling serves as a starting point for deeper characterization of how the yeast proteome changes during commercial fermentations and additionally serves as a resource to guide fermentation protocols, strain handling, and engineering practices in commercial brewing and fermentation environments. Finally, we created a web interface (https://brewing-yeast-proteomics.ccbb.utexas.edu/) to serve as a valuable resource for yeast geneticists, brewers, and biochemists to provide insights into the global trends underlying commercial beer production.

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.

Deep Learning Prediction Boosts Phosphoproteomics-Based Discoveries Through Improved Phosphopeptide Identification

Shotgun phosphoproteomics enables high-throughput analysis of phosphopeptides in biological samples. One of the primary challenges associated with this technology is the relatively low rate of phosphopeptide identification during data analysis. This limitation hampers the full realization of the potential offered by shotgun phosphoproteomics. Here we present DeepRescore2, a computational workflow that leverages deep learning-based retention time and fragment ion intensity predictions to improve phosphopeptide identification and phosphosite localization. Using a state-of-the-art computational workflow as a benchmark, DeepRescore2 increases the number of correctly identified peptide-spectrum matches by 17% in a synthetic dataset and identifies 19% to 46% more phosphopeptides in biological datasets. In a liver cancer dataset, 30% of the significantly altered phosphosites between tumor and normal tissues and 60% of the prognosis-associated phosphosites identified from DeepRescore2-processed data could not be identified based on the state-of-the-art workflow. Notably, DeepRescore2-processed data uniquely identifies EGFR hyperactivation as a new target in poor-prognosis liver cancer, which is validated experimentally. Integration of deep learning prediction in DeepRescore2 improves phosphopeptide identification and facilitates biological discoveries.

Identification of protein components of the transformation system in the cell line of immortalized human keratinocytes HaCaT exposed to surfactants

Using the method of shotgun mass spectrometry, we have evaluated changes in the proteomic profile of HaCat cells in response to the treatment with sodium dodecyl sulfate (anionic surfactant) and Triton-X100 (non-ionic surfactant) in two concentrations (12.5 µg/ml and 25.0 µg/ml). The study revealed induction of orphan CYP2S1 (biotransformation phase I) in response to Triton-X100. We have identified proteins of II (glutathione-S-transferases, GSTs) and III (solute carrier proteins, SLCs) biotransformation phases, as well as antioxidant proteins (peroxiredoxins, PRDXs; catalase, CAT; thioredoxin, TXN). Thus, proteins of all three xenobiotic detoxification phases were detected. The presented results suggest a new prospect of using HaCaT keratinocytes as a model of human epidermis for studying the metabolism of drugs/toxicants in human skin in vitro.

Interaction of Whitefly Effector G4 with Tomato Proteins Impacts Whitefly Performance

The phloem-feeding insect Bemisia tabaci is an important pest, responsible for the transmission of several crop-threatening virus species. While feeding, the insect secretes a cocktail of effectors to modulate plant defense responses. Here, we present a set of proteins identified in an artificial diet on which B. tabaci was salivating. We subsequently studied whether these candidate effectors can play a role in plant immune suppression. Effector G4 was the most robust suppressor of an induced- reactive oxygen species (ROS) response in Nicotiana benthamiana. In addition, G4 was able to suppress ROS production in Solanum lycopersicum (tomato) and Capsicum annuum (pepper). G4 localized predominantly in the endoplasmic reticulum in N. benthamiana leaves and colocalized with two identified target proteins in tomato: REF-like stress related protein 1 (RSP1) and meloidogyne-induced giant cell protein DB141 (MIPDB141). Silencing of MIPDB141 in tomato reduced whitefly fecundity up to 40%, demonstrating that the protein is involved in susceptibility to B. tabaci. Together, our data demonstrate that effector G4 impairs tomato immunity to whiteflies by interfering with ROS production and via an interaction with tomato susceptibility protein MIPDB141. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

VEuPathDB: the eukaryotic pathogen, vector and host bioinformatics resource center in 2023

The Eukaryotic Pathogen, Vector and Host Informatics Resource (VEuPathDB, https://veupathdb.org) is a Bioinformatics Resource Center funded by the National Institutes of Health with additional funding from the Wellcome Trust. VEuPathDB supports >600 organisms that comprise invertebrate vectors, eukaryotic pathogens (protists and fungi) and relevant free-living or non-pathogenic species or hosts. Since 2004, VEuPathDB has analyzed omics data from the public domain using contemporary bioinformatic workflows, including orthology predictions via OrthoMCL, and integrated the analysis results with analysis tools, visualizations, and advanced search capabilities. The unique data mining platform coupled with >3000 pre-analyzed data sets facilitates the exploration of pertinent omics data in support of hypothesis driven research. Comparisons are easily made across data sets, data types and organisms. A Galaxy workspace offers the opportunity for the analysis of private large-scale datasets and for porting to VEuPathDB for comparisons with integrated data. The MapVEu tool provides a platform for exploration of spatially resolved data such as vector surveillance and insecticide resistance monitoring. To address the growing body of omics data and advances in laboratory techniques, VEuPathDB has added several new data types, searches and features, improved the Galaxy workspace environment, redesigned the MapVEu interface and updated the infrastructure to accommodate these changes.

Web of venom: exploration of big data resources in animal toxin research

Research on animal venoms and their components spans multiple disciplines, including biology, biochemistry, bioinformatics, pharmacology, medicine, and more. Manipulating and analyzing the diverse array of data required for venom research can be challenging, and relevant tools and resources are often dispersed across different online platforms, making them less accessible to nonexperts. In this article, we address the multifaceted needs of the scientific community involved in venom and toxin-related research by identifying and discussing web resources, databases, and tools commonly used in this field. We have compiled these resources into a comprehensive table available on the VenomZone website (https://venomzone.expasy.org/10897). Furthermore, we highlight the challenges currently faced by researchers in accessing and using these resources and emphasize the importance of community-driven interdisciplinary approaches. We conclude by underscoring the significance of enhancing standards, promoting interoperability, and encouraging data and method sharing within the venom research community.

Metaproteomics for Coinfections in the Upper Respiratory Tract: The Case of COVID-19

The upper respiratory tract (URT) is home to a diverse range of microbial species. Respiratory infections disturb the microbial flora in the URT, putting people at risk of secondary infections. The potential dangers and clinical effects of bacterial and fungal coinfections with SARS-CoV-2 support the need to investigate the microbiome of the URT using clinical samples. Mass spectrometry (MS)-based metaproteomics analysis of microbial proteins is a novel approach to comprehensively assess the clinical specimens with complex microbial makeup. The coronavirus that causes severe acute respiratory syndrome (SARS-CoV-2) is responsible for the COVID-19 pandemic resulting in a plethora of microbial coinfections impeding therapy, prognosis, and overall disease management. In this chapter, the corresponding workflows for MS-based shotgun proteomics and metaproteomic analysis are illustrated.

Bioinformatic Workflows for Metaproteomics

The strong influence of microbiomes on areas such as ecology and human health has become widely recognized in the past years. Accordingly, various techniques for the investigation of the composition and function of microbial community samples have been developed. Metaproteomics, the comprehensive analysis of the proteins from microbial communities, allows for the investigation of not only the taxonomy but also the functional and quantitative composition of microbiome samples. Due to the complexity of the investigated communities, methods developed for single organism proteomics cannot be readily applied to metaproteomic samples. For this purpose, methods specifically tailored to metaproteomics are required. In this work, a detailed overview of current bioinformatic solutions and protocols in metaproteomics is given. After an introduction to the proteomic database search, the metaproteomic post-processing steps are explained in detail. Ten specific bioinformatic software solutions are focused on, covering various steps including database-driven identification and quantification as well as taxonomic and functional assignment.

Biodiversity Analysis of Metaproteomics Samples with Unipept: A Comprehensive Tutorial

Metaproteomics has become a crucial omics technology for studying microbiomes. In this area, the Unipept ecosystem, accessible at https://unipept.ugent.be , has emerged as a valuable resource for analyzing metaproteomic data. It offers in-depth insights into both taxonomic distributions and functional characteristics of complex ecosystems. This tutorial explains essential concepts like Lowest Common Ancestor (LCA) determination and the handling of peptides with missed cleavages. It also provides a detailed, step-by-step guide on using the Unipept Web application and Unipept Desktop for thorough metaproteomics analyses. By integrating theoretical principles with practical methodologies, this tutorial empowers researchers with the essential knowledge and tools needed to fully utilize metaproteomics in their microbiome studies.

Forecasting the dynamics of a complex microbial community using integrated meta-omics

Predicting the behaviour of complex microbial communities is challenging. However, this is essential for complex biotechnological processes such as those in biological wastewater treatment plants (BWWTPs), which require sustainable operation. Here we summarize 14 months of longitudinal meta-omics data from a BWWTP anaerobic tank into 17 temporal signals, explaining 91.1% of the temporal variance, and link those signals to ecological events within the community. We forecast the signals over the subsequent five years and use 21 extra samples collected at defined time intervals for testing and validation. Our forecasts are correct for six signals and hint on phenomena such as predation cycles. Using all the 17 forecasts and the environmental variables, we predict gene abundance and expression, with a coefficient of determination ≥0.87 for the subsequent three years. Our study demonstrates the ability to forecast the dynamics of open microbial ecosystems using interactions between community cycles and environmental parameters.

Exploring the Alternative Proteome with OpenProt and Mass Spectrometry

Proteogenomics has revealed the translation of unannotated open reading frames (ORFs) present in mRNAs and in noncoding RNAs (ncRNAs). OpenProt annotates all ORFs with a minimum of 30 codons in the transcriptome of several species and displays many functional features associated with the corresponding proteins. Two types of proteins are annotated: reference or canonical proteins which are proteins already annotated in UniProt, RefSeq, or Ensembl and noncanonical proteins. Noncanonical proteins form two groups: predicted novel isoforms that display a significant level of homology with a reference protein and alternative proteins that are new proteins with no significant homology to known proteins. This chapter describes how to check whether a gene and/or transcript contains multiple open reading frames and how to use OpenProt databases for the detection of alternative proteins and novel isoforms by mass spectrometry-based proteomics.

The Application of Ejaculate-Based Shotgun Proteomics for Male Infertility Screening

Problems with the male reproductive system are of both medical and social significance. As a rule, spermatozoa and seminal plasma proteomes are investigated separately to assess sperm quality. The current study aimed to compare ejaculate proteomes with spermatozoa and seminal plasma protein profiles regarding the identification of proteins related to fertility scores. A total of 1779, 715, and 2163 proteins were identified in the ejaculate, seminal plasma, and spermatozoa, respectively. Among these datasets, 472 proteins were shared. GO enrichment analysis of the common proteins enabled us to distinguish biological processes such as single fertilization (GO:0007338), spermatid development (GO:0007286), and cell motility (GO:0048870). Among the abundant terms for GO cellular components, zona pellucida receptor complex, sperm fibrous sheath, and outer dense fiber were revealed. Overall, we identified 139 testis-specific proteins. For these proteins, PPI networks that are common in ejaculate, spermatozoa, and seminal plasma were related to the following GO biological processes: cilium movement (GO:0003341), microtubule-based movement (GO:0007018), and sperm motility (GO:0097722). For ejaculate and spermatozoa, they shared 15 common testis-specific proteins with spermatogenesis (GO:0007283) and male gamete generation (GO:0048232). Therefore, we speculated that ejaculate-based proteomics could yield new insights into the peculiar reproductive physiology and spermatozoa function of men and potentially serve as an explanation for male infertility screening.

Fe protein docking transduces conformational changes to MoFe nitrogenase active site in a nucleotide-dependent manner

The reduction of dinitrogen to ammonia catalyzed by nitrogenase involves a complex series of events, including ATP hydrolysis, electron transfer, and activation of metal clusters for N2 reduction. Early evidence shows that an essential part of the mechanism involves transducing information between the nitrogenase component proteins through conformational dynamics. Here, millisecond time-resolved hydrogen-deuterium exchange mass spectrometry was used to unravel peptide-level protein motion on the time scale of catalysis of Mo-dependent nitrogenase from Azotobacter vinelandii. Normal mode analysis calculations complemented this data, providing insights into the specific signal transduction pathways that relay information across protein interfaces at distances spanning 100 Å. Together, these results show that conformational changes induced by protein docking are rapidly transduced to the active site, suggesting a specific mechanism for activating the metal cofactor in the enzyme active site.

Sage: An Open-Source Tool for Fast Proteomics Searching and Quantification at Scale

The growing complexity and volume of proteomics data necessitate the development of efficient software tools for peptide identification and quantification from mass spectra. Given their central role in proteomics, it is imperative that these tools are auditable and extensible─requirements that are best fulfilled by open-source and permissively licensed software. This work presents Sage, a high-performance, open-source, and freely available proteomics pipeline. Scalable and cloud-ready, Sage matches the performance of state-of-the-art software tools while running an order of magnitude faster.

Differential Expression of Proteins Associated with Bipolar Disorder as Identified Using the PeptideShaker Software

The prevalence of bipolar disorder (BD) in modern society is growing rapidly, but due to the lack of paraclinical criteria, its differential diagnosis with other mental disorders is somewhat challenging. In this regard, the relevance of proteomic studies is increasing due to the development of methods for processing large data arrays; this contributes to the discovery of protein patterns of pathological processes and the creation of new methods of diagnosis and treatment. It seems promising to search for proteins involved in the pathogenesis of BD in an easily accessible material-blood serum. Sera from BD patients and healthy individuals were purified via affinity chromatography to isolate 14 major proteins and separated using 1D SDS-PAGE. After trypsinolysis, the proteins in the samples were identified via HPLC/mass spectrometry. Mass spectrometric data were processed using the OMSSA and X!Tandem search algorithms using the UniProtKB database, and the results were analyzed using PeptideShaker. Differences in proteomes were assessed via an unlabeled NSAF-based analysis using a two-tailed Bonferroni-adjusted t-test. When comparing the blood serum proteomes of BD patients and healthy individuals, 10 proteins showed significant differences in NSAF values. Of these, four proteins were predominantly present in BD patients with the maximum NSAF value: 14-3-3 protein zeta/delta; ectonucleoside triphosphate diphosphohydrolase 7; transforming growth factor-beta-induced protein ig-h3; and B-cell CLL/lymphoma 9 protein. Further exploration of the role of these proteins in BD is warranted; conducting such studies will help develop new paraclinical criteria and discover new targets for BD drug therapy.

Increasing taxonomic and functional characterization of host-microbiome interactions by DIA-PASEF metaproteomics

Introduction

Metaproteomics is a rapidly advancing field that offers unique insights into the taxonomic composition and the functional activity of microbial communities, and their effects on host physiology. Classically, data-dependent acquisition (DDA) mass spectrometry (MS) has been applied for peptide identification and quantification in metaproteomics. However, DDA-MS exhibits well-known limitations in terms of depth, sensitivity, and reproducibility. Consequently, methodological improvements are required to better characterize the protein landscape of microbiomes and their interactions with the host.

Methods

We present an optimized proteomic workflow that utilizes the information captured by Parallel Accumulation-Serial Fragmentation (PASEF) MS for comprehensive metaproteomic studies in complex fecal samples of mice.

Results and discussion

We show that implementing PASEF using a DDA acquisition scheme (DDA-PASEF) increased peptide quantification up to 5 times and reached higher accuracy and reproducibility compared to previously published classical DDA and data-independent acquisition (DIA) methods. Furthermore, we demonstrate that the combination of DIA, PASEF, and neuronal-network-based data analysis, was superior to DDA-PASEF in all mentioned parameters. Importantly, DIA-PASEF expanded the dynamic range towards low-abundant proteins and it doubled the quantification of proteins with unknown or uncharacterized functions. Compared to previous classical DDA metaproteomic studies, DIA-PASEF resulted in the quantification of up to 4 times more taxonomic units using 16 times less injected peptides and 4 times shorter chromatography gradients. Moreover, 131 additional functional pathways distributed across more and even uniquely identified taxa were profiled as revealed by a peptide-centric taxonomic-functional analysis. We tested our workflow on a validated preclinical mouse model of neuropathic pain to assess longitudinal changes in host-gut microbiome interactions associated with pain - an unexplored topic for metaproteomics. We uncovered the significant enrichment of two bacterial classes upon pain, and, in addition, the upregulation of metabolic activities previously linked to chronic pain as well as various hitherto unknown ones. Furthermore, our data revealed pain-associated dynamics of proteome complexes implicated in the crosstalk between the host immune system and the gut microbiome. In conclusion, the DIA-PASEF metaproteomic workflow presented here provides a stepping stone towards a deeper understanding of microbial ecosystems across the breadth of biomedical and biotechnological fields.

Similar proteome expression profiles of the aggregated lymphoid nodules area and Peyer's patches in Bactrian camel

BACKGROUND

The presence of Aggregated Lymphoid Nodules Area (ALNA) is a notable anatomical characteristic observed in the abomasum of Bactrian camels. This area is comprised of two separate regions, namely the Reticular Mucosal Folds Region (RMFR) and the Longitudinal Mucosal Folds Region (LMFR). The histological properties of ALNA exhibit significant similarities to those of Peyer's patches (PPs) found in the gastrointestinal system. The functional characteristics of ALNA were examined in relation to mucosal immunity in the gastrointestinal system.

RESULTS

We used iTRAQ-based proteomic analysis on twelve Bactrian camels to measure the amount of proteins expressed in ALNA. In the experiment, we sampled the RMFR and LMFR separately from the ALNA and compared their proteomic quantification results with samples from the PPs. A total of 1253 proteins were identified, among which 39 differentially expressed proteins (DEPs) were found between RMFR and PPs, 33 DEPs were found between LMFR and PPs, and 22 DEPs were found between LMFR and RMFR. The proteins FLNA, MYH11, and HSPB1 were chosen for validation using the enzyme-linked immunosorbent assay (ELISA), and the observed expression profiles were found to be in agreement with the results obtained from the iTRAQ study. The InnateDB database was utilized to get data pertaining to immune-associated proteins in ALNA. It was observed that a significant proportion, specifically 76.6%, of these proteins were found to be associated with the same orthogroups as human immune-related genes. These proteins are acknowledged to be associated with a diverse range of functions, encompassing the uptake, processing and presentation of antigens, activation of lymphocytes, the signaling pathways of T-cell and B-cell receptors, and the control of actin polymerization.

CONCLUSIONS

The experimental results suggest that there are parallels in the immune-related proteins found in ALNA and PPs. Although there are variations in the structures of LMFR and RMFR, the proteins produced in both structures exhibit a high degree of similarity and perform comparable functions in the context of mucosal immune responses.

Retention Time and Fragmentation Predictors Increase Confidence in Identification of Common Variant Peptides

Precision medicine focuses on adapting care to the individual profile of patients, for example, accounting for their unique genetic makeup. Being able to account for the effect of genetic variation on the proteome holds great promise toward this goal. However, identifying the protein products of genetic variation using mass spectrometry has proven very challenging. Here we show that the identification of variant peptides can be improved by the integration of retention time and fragmentation predictors into a unified proteogenomic pipeline. By combining these intrinsic peptide characteristics using the search-engine post-processor Percolator, we demonstrate improved discrimination power between correct and incorrect peptide-spectrum matches. Our results demonstrate that the drop in performance that is induced when expanding a protein sequence database can be compensated, hence enabling efficient identification of genetic variation products in proteomics data. We anticipate that this enhancement of proteogenomic pipelines can provide a more refined picture of the unique proteome of patients and thereby contribute to improving patient care.

Assessment and Comparison of Database Search Engines for Peptidomic Applications

Protein database search engines are an integral component of mass spectrometry-based peptidomic analyses. Given the unique computational challenges of peptidomics, many factors must be taken into consideration when optimizing search engine selection, as each platform has different algorithms by which tandem mass spectra are scored for subsequent peptide identifications. In this study, four different database search engines, PEAKS, MS-GF+, OMSSA, and X! Tandem, were compared with Aplysia californica and Rattus norvegicus peptidomics data sets, and various metrics were assessed such as the number of unique peptide and neuropeptide identifications, and peptide length distributions. Given the tested conditions, PEAKS was found to have the highest number of peptide and neuropeptide identifications out of the four search engines in both data sets. Furthermore, principal component analysis and multivariate logistic regression were employed to determine whether specific spectral features contribute to false C-terminal amidation assignments by each search engine. From this analysis, it was found that the primary features influencing incorrect peptide assignments were the precursor and fragment ion m/z errors. Finally, an assessment employing a mixed species protein database was performed to evaluate search engine precision and sensitivity when searched against an enlarged search space containing human proteins.

Systematic Profiling of Ale Yeast Protein Dynamics across Fermentation and Repitching

Studying the genetic and molecular characteristics of brewing yeast strains is crucial for understanding their domestication history and adaptations accumulated over time in fermentation environments, and for guiding optimizations to the brewing process itself. Saccharomyces cerevisiae (brewing yeast) is amongst the most profiled organisms on the planet, yet the temporal molecular changes that underlie industrial fermentation and beer brewing remain understudied. Here, we characterized the genomic makeup of a Saccharomyces cerevisiae ale yeast widely used in the production of Hefeweizen beers, and applied shotgun mass spectrometry to systematically measure the proteomic changes throughout two fermentation cycles which were separated by 14 rounds of serial repitching. The resulting brewing yeast proteomics resource includes 64,740 protein abundance measurements. We found that this strain possesses typical genetic characteristics of Saccharomyces cerevisiae ale strains and displayed progressive shifts in molecular processes during fermentation based on protein abundance changes. We observed protein abundance differences between early fermentation batches compared to those separated by 14 rounds of serial repitching. The observed abundance differences occurred mainly in proteins involved in the metabolism of ergosterol and isobutyraldehyde. Our systematic profiling serves as a starting point for deeper characterization of how the yeast proteome changes during commercial fermentations and additionally serves as a resource to guide fermentation protocols, strain handling, and engineering practices in commercial brewing and fermentation environments. Finally, we created a web interface (https://brewing-yeast-proteomics.ccbb.utexas.edu/) to serve as a valuable resource for yeast geneticists, brewers, and biochemists to provide insights into the global trends underlying commercial beer production.

ERFVII action and modulation through oxygen-sensing in Arabidopsis thaliana

Oxygen is a key signalling component of plant biology, and whilst an oxygen-sensing mechanism was previously described in Arabidopsis thaliana, key features of the associated PLANT CYSTEINE OXIDASE (PCO) N-degron pathway and Group VII ETHYLENE RESPONSE FACTOR (ERFVII) transcription factor substrates remain untested or unknown. We demonstrate that ERFVIIs show non-autonomous activation of root hypoxia tolerance and are essential for root development and survival under oxygen limiting conditions in soil. We determine the combined effects of ERFVIIs in controlling gene expression and define genetic and environmental components required for proteasome-dependent oxygen-regulated stability of ERFVIIs through the PCO N-degron pathway. Using a plant extract, unexpected amino-terminal cysteine sulphonic acid oxidation level of ERFVIIs was observed, suggesting a requirement for additional enzymatic activity within the pathway. Our results provide a holistic understanding of the properties, functions and readouts of this oxygen-sensing mechanism defined through its role in modulating ERFVII stability.

Interruption of p53-MDM2 Interaction by Nutlin-3a in Human Lymphoma Cell Models Initiates a Cell-Dependent Global Effect on Transcriptome and Proteome Level

In most lymphomas, p53 signaling pathway is inactivated by various mechanisms independent to p53 gene mutations or deletions. In many cases, p53 function is largely regulated by alterations in the protein abundance levels by the action of E3 ubiquitin-protein ligase MDM2, targeting p53 to proteasome-mediated degradation. In the present study, an integrating transcriptomics and proteomics analysis was employed to investigate the effect of p53 activation by a small-molecule MDM2-antagonist, nutlin-3a, on three lymphoma cell models following p53 activation. Our analysis revealed a system-wide nutlin-3a-associated effect in all examined lymphoma types, identifying in total of 4037 differentially affected proteins involved in a plethora of pathways, with significant heterogeneity among lymphomas. Our findings include known p53-targets and novel p53 activation effects, involving transcription, translation, or degradation of protein components of pathways, such as a decrease in key members of PI3K/mTOR pathway, heat-shock response, and glycolysis, and an increase in key members of oxidative phoshosphorylation, autophagy and mitochondrial translation. Combined inhibition of HSP90 or PI3K/mTOR pathway with nutlin-3a-mediated p53-activation enhanced the apoptotic effects suggesting a promising strategy against human lymphomas. Integrated omic profiling after p53 activation offered novel insights on the regulatory role specific proteins and pathways may have in lymphomagenesis.

Insight on physicochemical properties governing peptide MS1 response in HPLC-ESI-MS/MS: A deep learning approach

Accurate and absolute quantification of peptides in complex mixtures using quantitative mass spectrometry (MS)-based methods requires foreground knowledge and isotopically labeled standards, thereby increasing analytical expenses, time consumption, and labor, thus limiting the number of peptides that can be accurately quantified. This originates from differential ionization efficiency between peptides and thus, understanding the physicochemical properties that influence the ionization and response in MS analysis is essential for developing less restrictive label-free quantitative methods. Here, we used equimolar peptide pool repository data to develop a deep learning model capable of identifying amino acids influencing the MS1 response. By using an encoder-decoder with an attention mechanism and correlating attention weights with amino acid physicochemical properties, we obtain insight on properties governing the peptide-level MS1 response within the datasets. While the problem cannot be described by one single set of amino acids and properties, distinct patterns were reproducibly obtained. Properties are grouped in three main categories related to peptide hydrophobicity, charge, and structural propensities. Moreover, our model can predict MS1 intensity output under defined conditions based solely on peptide sequence input. Using a refined training dataset, the model predicted log-transformed peptide MS1 intensities with an average error of 9.7 ± 0.5% based on 5-fold cross validation, and outperformed random forest and ridge regression models on both log-transformed and real scale data. This work demonstrates how deep learning can facilitate identification of physicochemical properties influencing peptide MS1 responses, but also illustrates how sequence-based response prediction and label-free peptide-level quantification may impact future workflows within quantitative proteomics.

The salivary proteome in relation to oral mucositis in autologous hematopoietic stem cell transplantation recipients: a labelled and label-free proteomics approach

BACKGROUND

Oral mucositis is a frequently seen complication in the first weeks after hematopoietic stem cell transplantation recipients which can severely affects patients quality of life. In this study, a labelled and label-free proteomics approach were used to identify differences between the salivary proteomes of autologous hematopoietic stem cell transplantation (ASCT) recipients developing ulcerative oral mucositis (ULC-OM; WHO score ≥ 2) or not (NON-OM).

METHODS

In the TMT-labelled analysis we pooled saliva samples from 5 ULC-OM patients at each of 5 timepoints: baseline, 1, 2, 3 weeks and 3 months after ASCT and compared these with pooled samples from 5 NON-OM patients. For the label-free analysis we analyzed saliva samples from 9 ULC-OM and 10 NON-OM patients at 6 different timepoints (including 12 months after ASCT) with Data-Independent Acquisition (DIA). As spectral library, all samples were grouped (ULC-OM vs NON-OM) and analyzed with Data Dependent Analysis (DDA). PCA plots and a volcano plot were generated in RStudio and differently regulated proteins were analyzed using GO analysis with g:Profiler.

RESULTS

A different clustering of ULC-OM pools was found at baseline, weeks 2 and 3 after ASCT with TMT-labelled analysis. Using label-free analysis, week 1-3 samples clustered distinctly from the other timepoints. Unique and up-regulated proteins in the NON-OM group (DDA analysis) were involved in immune system-related processes, while those proteins in the ULC-OM group were intracellular proteins indicating cell lysis.

CONCLUSIONS

The salivary proteome in ASCT recipients has a tissue protective or tissue-damage signature, that corresponded with the absence or presence of ulcerative oral mucositis, respectively.

TRIAL REGISTRATION

The study is registered in the national trial register (NTR5760; automatically added to the International Clinical Trial Registry Platform).

A Galaxy of informatics resources for MS-based proteomics

INTRODUCTION

Continuous advances in mass spectrometry (MS) technologies have enabled deeper and more reproducible proteome characterization and a better understanding of biological systems when integrated with other 'omics data. Bioinformatic resources meeting the analysis requirements of increasingly complex MS-based proteomic data and associated multi-omic data are critically needed. These requirements included availability of software that would span diverse types of analyses, scalability for large-scale, compute-intensive applications, and mechanisms to ease adoption of the software.

AREAS COVERED

The Galaxy ecosystem meets these requirements by offering a multitude of open-source tools for MS-based proteomics analyses and applications, all in an adaptable, scalable, and accessible computing environment. A thriving global community maintains these software and associated training resources to empower researcher-driven analyses.

EXPERT OPINION

The community-supported Galaxy ecosystem remains a crucial contributor to basic biological and clinical studies using MS-based proteomics. In addition to the current status of Galaxy-based resources, we describe ongoing developments for meeting emerging challenges in MS-based proteomic informatics. We hope this review will catalyze increased use of Galaxy by researchers employing MS-based proteomics and inspire software developers to join the community and implement new tools, workflows, and associated training content that will add further value to this already rich ecosystem.

LFQ-Based Peptide and Protein Intensity Differential Expression Analysis

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

Small extracellular vesicles but not microvesicles from Opisthorchis viverrini promote cell proliferation in human cholangiocytes

Chronic infection with O. viverrini has been linked to the development of cholangiocarcinoma (CCA), which is a major public health burden in the Lower Mekong River Basin countries, including Thailand, Lao PDR, Vietnam and Cambodia. Despite its importance, the exact mechanisms by which O. viverrini promotes CCA are largely unknown. In this study, we characterized different extracellular vesicle populations released by O. viverrini (OvEVs) using proteomic and transcriptomic analyses and investigated their potential role in host-parasite interactions. While 120k OvEVs promoted cell proliferation in H69 cells at different concentrations, 15k OvEVs did not produce any effect compared to controls. The proteomic analysis of both populations showed differences in their composition that could contribute to this differential effect. Furthermore, the miRNAs present in 120k EVs were analysed and their potential interactions with human host genes was explored by computational target prediction. Different pathways involved in inflammation, immune response and apoptosis were identified as potentially targeted by the miRNAs present in this population of EVs. This is the first study showing specific roles for different EV populations in the pathogenesis of a parasitic helminth, and more importantly, an important advance towards deciphering the mechanisms used in establishment of opisthorchiasis and liver fluke infection-associated malignancy.

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

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

Comparison of SPEED, S-Trap, and In-Solution-Based Sample Preparation Methods for Mass Spectrometry in Kidney Tissue and Plasma

Mass spectrometry is a powerful technique for investigating renal pathologies and identifying biomarkers, and efficient protein extraction from kidney tissue is essential for bottom-up proteomic analyses. Detergent-based strategies aid cell lysis and protein solubilization but are poorly compatible with downstream protein digestion and liquid chromatography-coupled mass spectrometry, requiring additional purification and buffer-exchange steps. This study compares two well-established detergent-based methods for protein extraction (in-solution sodium deoxycholate (SDC); suspension trapping (S-Trap)) with the recently developed sample preparation by easy extraction and digestion (SPEED) method, which uses strong acid for denaturation. We compared the quantitative performance of each method using label-free mass spectrometry in both sheep kidney cortical tissue and plasma. In kidney tissue, SPEED quantified the most unique proteins (SPEED 1250; S-Trap 1202; SDC 1197). In plasma, S-Trap produced the most unique protein quantifications (S-Trap 150; SDC 148; SPEED 137). Protein quantifications were reproducible across biological replicates in both tissue (R2 = 0.85–0.90) and plasma (SPEED R2 = 0.84; SDC R2 = 0.76, S-Trap R2 = 0.65). Our data suggest SPEED as the optimal method for proteomic preparation in kidney tissue and S-Trap or SPEED as the optimal method for plasma, depending on whether a higher number of protein quantifications or greater reproducibility is desired.

Dynamic states of eIF6 and SDS variants modulate interactions with uL14 of the 60S ribosomal subunit

Assembly of ribosomal subunits into active ribosomal complexes is integral to protein synthesis. Release of eIF6 from the 60S ribosomal subunit primes 60S to associate with the 40S subunit and engage in translation. The dynamics of eIF6 interaction with the uL14 (RPL23) interface of 60S and its perturbation by somatic mutations acquired in Shwachman-Diamond Syndrome (SDS) is yet to be clearly understood. Here, by using a modified strategy to obtain high yields of recombinant human eIF6 we have uncovered the critical interface entailing eight key residues in the C-tail of uL14 that is essential for physical interactions between 60S and eIF6. Disruption of the complementary binding interface by conformational changes in eIF6 disease variants provide a mechanism for weakened interactions of variants with the 60S. Hydrogen-deuterium exchange mass spectrometry (HDX-MS) analyses uncovered dynamic configurational rearrangements in eIF6 induced by binding to uL14 and exposed an allosteric interface regulated by the C-tail of eIF6. Disrupting key residues in the eIF6-60S binding interface markedly limits proliferation of cancer cells, which highlights the significance of therapeutically targeting this interface. Establishing these key interfaces thus provide a therapeutic framework for targeting eIF6 in cancers and SDS.

Probing the E1o-E2o and E1a-E2o Interactions in Binary Subcomplexes of the Human 2-Oxoglutarate Dehydrogenase and 2-Oxoadipate Dehydrogenase Complexes by Chemical Cross-Linking Mass Spectrometry and Molecular Dynamics Simulation

The human 2-oxoglutarate dehydrogenase complex (hOGDHc) is a key enzyme in the tricarboxylic acid cycle and is one of the main regulators of mitochondrial metabolism through NADH and reactive oxygen species levels. Evidence was obtained for formation of a hybrid complex between the hOGDHc and its homologue the 2-oxoadipate dehydrogenase complex (hOADHc) in the L-lysine metabolic pathway, suggesting a crosstalk between the two distinct pathways. Findings raised fundamental questions about the assembly of hE1a (2-oxoadipate-dependent E1 component) and hE1o (2-oxoglutarate-dependent E1) to the common hE2o core component. Here we report chemical cross-linking mass spectrometry (CL-MS) and molecular dynamics (MD) simulation analyses to understand assembly in binary subcomplexes. The CL-MS studies revealed the most prominent loci for hE1o-hE2o and hE1a-hE2o interactions and suggested different binding modes. The MD simulation studies led to the following conclusions: (i) The N-terminal regions in E1s are shielded by, but do not interact directly with hE2o. (ii) The hE2o linker region exhibits the highest number of H-bonds with the N-terminus and α/β1 helix of hE1o, yet with the interdomain linker and α/β1 helix of hE1a. (iii) The C-termini are involved in dynamic interactions in complexes, suggesting the presence of at least two conformations in solution.

MetaProD: A Highly-Configurable Mass Spectrometry Analyzer for Multiplexed Proteomic and Metaproteomic Data

The microbiome has been shown to be important for human health because of its influence on disease and the immune response. Mass spectrometry is an important tool for evaluating protein expression and species composition in the microbiome but is technically challenging and time-consuming. Multiplexing has emerged as a way to make spectrometry workflows faster while improving results. Here, we present MetaProD (MetaProteomics in Django) as a highly configurable metaproteomic data analysis pipeline supporting label-free and multiplexed mass spectrometry. The pipeline is open-source, uses fully open-source tools, and is integrated with Django to offer a web-based interface for configuration and data access. Benchmarking of MetaProD using multiple metaproteomics data sets showed that MetaProD achieved fast and efficient identification of peptides and proteins. Application of MetaProD to a multiplexed cancer data set resulted in identification of more differentially expressed human proteins in cancer tissues versus healthy tissues as compared to previous studies; in addition, MetaProD identified bacterial proteins in those samples, some of which are differentially abundant.

MetaLab-MAG: A Metaproteomic Data Analysis Platform for Genome-Level Characterization of Microbiomes from the Metagenome-Assembled Genomes Database

The studies of microbial communities have drawn increased attention in various research fields such as agriculture, environment, and human health. Recently, metaproteomics has become a powerful tool to interpret the roles of the community members by investigating the expressed proteins of the microbes. However, analyzing the metaproteomic data sets at genome resolution is still challenging because of the lack of efficient bioinformatics tools. Here we develop MetaLab-MAG, a specially designed tool for the characterization of microbiomes from metagenome-assembled genomes databases. MetaLab-MAG was evaluated by analyzing various human gut microbiota data sets and performed comparably or better than searching the gene catalog protein database directly. MetaLab-MAG can quantify the genome-level microbiota compositions and supports both label-free and isobaric labeling-based quantification strategies. MetaLab-MAG removes the obstacles of metaproteomic data analysis and provides the researchers with in-depth and comprehensive information from the microbiomes.

[Antibody proteomics]

Antibodies represent an essential component of humoral immunity; therefore their study is important for molecular biology and medicine. The unique property of antibodies to specifically recognize and bind a certain molecular target (an antigen) determines their widespread application in treatment and diagnostics of diseases, as well as in laboratory and biotechnological practices. High specificity and affinity of antibodies is determined by the presence of primary structure variable regions, which are not encoded in the human genome and are unique for each antibody-producing B cell clone. Hence, there is little or no information about amino acid sequences of the variable regions in the databases. This differs identification of antibody primary structure from most of the proteomic studies because it requires either B cell genome sequencing or de novo amino acid sequencing of the antibody. The present review demonstrates some examples of proteomic and proteogenomic approaches and the methodological arsenal that proteomics can offer for studying antibodies, in particular, for identification of primary structure, evaluation of posttranslational modifications and application of bioinformatics tools for their decoding.