TRIC: an automated alignment strategy for reproducible protein quantification in targeted proteomics

Proteomics appending a complementary dimension to precision oncotherapy

Recent advances in high-throughput proteomic profiling technologies have facilitated the precise quantification of numerous proteins across multiple specimens concurrently. Researchers have the opportunity to comprehensively analyze the molecular signatures in plentiful medical specimens or disease pattern cell lines. Along with advances in data analysis and integration, proteomics data could be efficiently consolidated and employed to recognize precise elementary molecular mechanisms and decode individual biomarkers, guiding the precision treatment of tumors. Herein, we review a broad array of proteomics technologies and the progress and methods for the integration of proteomics data and further discuss how to better merge proteomics in precision medicine and clinical settings.

Targeting low levels of MIF expression as a potential therapeutic strategy for ALS

Mutations in SOD1 cause amyotrophic lateral sclerosis (ALS), a neurodegenerative disease characterized by motor neuron (MN) loss. We previously discovered that macrophage migration inhibitory factor (MIF), whose levels are extremely low in spinal MNs, inhibits mutant SOD1 misfolding and toxicity. In this study, we show that a single peripheral injection of adeno-associated virus (AAV) delivering MIF into adult SOD1G37R mice significantly improves their motor function, delays disease progression, and extends survival. Moreover, MIF treatment reduces neuroinflammation and misfolded SOD1 accumulation, rescues MNs, and corrects dysregulated pathways as observed by proteomics and transcriptomics. Furthermore, we reveal low MIF levels in human induced pluripotent stem cell-derived MNs from familial ALS patients with different genetic mutations, as well as in post mortem tissues of sporadic ALS patients. Our findings indicate that peripheral MIF administration may provide a potential therapeutic mechanism for modulating misfolded SOD1 in vivo and disease outcome in ALS patients.

Network-based elucidation of colon cancer drug resistance mechanisms by phosphoproteomic time-series analysis

Aberrant signaling pathway activity is a hallmark of tumorigenesis and progression, which has guided targeted inhibitor design for over 30 years. Yet, adaptive resistance mechanisms, induced by rapid, context-specific signaling network rewiring, continue to challenge therapeutic efficacy. Leveraging progress in proteomic technologies and network-based methodologies, we introduce Virtual Enrichment-based Signaling Protein-activity Analysis (VESPA)-an algorithm designed to elucidate mechanisms of cell response and adaptation to drug perturbations-and use it to analyze 7-point phosphoproteomic time series from colorectal cancer cells treated with clinically-relevant inhibitors and control media. Interrogating tumor-specific enzyme/substrate interactions accurately infers kinase and phosphatase activity, based on their substrate phosphorylation state, effectively accounting for signal crosstalk and sparse phosphoproteome coverage. The analysis elucidates time-dependent signaling pathway response to each drug perturbation and, more importantly, cell adaptive response and rewiring, experimentally confirmed by CRISPR knock-out assays, suggesting broad applicability to cancer and other diseases.

Quantification of Adaptive Immune Responses Against Protein-Binding Interfaces in the Streptococcal M1 Protein

Bacterial or viral antigens can contain subdominant protein regions that elicit weak antibody responses upon vaccination or infection although there is accumulating evidence that antibody responses against subdominant regions can enhance the protective immune response. One proposed mechanism for subdominant protein regions is the binding of host proteins that prevent antibody production against epitopes hidden within the protein binding interfaces. Here, we used affinity purification combined with quantitative mass spectrometry (AP-MS) to examine the level of competition between antigen-specific antibodies and host-pathogen protein interaction networks using the M1 protein from Streptococcus pyogenes as a model system. As most humans have circulating antibodies against the M1 protein, we first used AP-MS to show that the M1 protein interspecies protein network formed with human plasma proteins is largely conserved in naïve mice. Immunizing mice with the M1 protein generated a time-dependent increase of anti-M1 antibodies. AP-MS analysis comparing the composition of the M1-plasma protein network from naïve and immunized mice showed significant enrichment of 292 IgG peptides associated with 56 IgG chains in the immune mice. Despite the significant increase of bound IgGs, the levels of interacting plasma proteins were not significantly reduced in the immune mice. The results indicate that the antigen-specific polyclonal IgG against the M1 protein primarily targets epitopes outside the other plasma protein binding interfaces. In conclusion, this study demonstrates that AP-MS is a promising strategy to determine the relationship between antigen-specific antibodies and host-pathogen interaction networks that could be used to define subdominant protein regions of relevance for vaccine development.

Longitudinal profiling of the microbiome at four body sites reveals core stability and individualized dynamics during health and disease

To understand the dynamic interplay between the human microbiome and host during health and disease, we analyzed the microbial composition, temporal dynamics, and associations with host multi-omics, immune, and clinical markers of microbiomes from four body sites in 86 participants over 6 years. We found that microbiome stability and individuality are body-site specific and heavily influenced by the host. The stool and oral microbiome are more stable than the skin and nasal microbiomes, possibly due to their interaction with the host and environment. We identify individual-specific and commonly shared bacterial taxa, with individualized taxa showing greater stability. Interestingly, microbiome dynamics correlate across body sites, suggesting systemic dynamics influenced by host-microbial-environment interactions. Notably, insulin-resistant individuals show altered microbial stability and associations among microbiome, molecular markers, and clinical features, suggesting their disrupted interaction in metabolic disease. Our study offers comprehensive views of multi-site microbial dynamics and their relationship with host health and disease.

SeFilter-DIA: Squeeze-and-Excitation Network for Filtering High-Confidence Peptides of Data-Independent Acquisition Proteomics

Mass spectrometry is crucial in proteomics analysis, particularly using Data Independent Acquisition (DIA) for reliable and reproducible mass spectrometry data acquisition, enabling broad mass-to-charge ratio coverage and high throughput. DIA-NN, a prominent deep learning software in DIA proteome analysis, generates peptide results but may include low-confidence peptides. Conventionally, biologists have to manually screen peptide fragment ion chromatogram peaks (XIC) for identifying high-confidence peptides, a time-consuming and subjective process prone to variability. In this study, we introduce SeFilter-DIA, a deep learning algorithm, aiming at automating the identification of high-confidence peptides. Leveraging compressed excitation neural network and residual network models, SeFilter-DIA extracts XIC features and effectively discerns between high and low-confidence peptides. Evaluation of the benchmark datasets demonstrates SeFilter-DIA achieving 99.6% AUC on the test set and 97% for other performance indicators. Furthermore, SeFilter-DIA is applicable for screening peptides with phosphorylation modifications. These results demonstrate the potential of SeFilter-DIA to replace manual screening, providing an efficient and objective approach for high-confidence peptide identification while mitigating associated limitations.

Substantial downregulation of mitochondrial and peroxisomal proteins during acute kidney injury revealed by data-independent acquisition proteomics

Acute kidney injury (AKI) manifests as a major health concern, particularly for the elderly. Understanding AKI-related proteome changes is critical for prevention and development of novel therapeutics to recover kidney function and to mitigate the susceptibility for recurrent AKI or development of chronic kidney disease. In this study, mouse kidneys were subjected to ischemia-reperfusion injury, and the contralateral kidneys remained uninjured to enable comparison and assess injury-induced changes in the kidney proteome. A ZenoTOF 7600 mass spectrometer was optimized for data-independent acquisition (DIA) to achieve comprehensive protein identification and quantification. Short microflow gradients and the generation of a deep kidney-specific spectral library allowed for high-throughput, comprehensive protein quantification. Upon AKI, the kidney proteome was completely remodeled, and over half of the 3945 quantified protein groups changed significantly. Downregulated proteins in the injured kidney were involved in energy production, including numerous peroxisomal matrix proteins that function in fatty acid oxidation, such as ACOX1, CAT, EHHADH, ACOT4, ACOT8, and Scp2. Injured kidneys exhibited severely damaged tissues and injury markers. The comprehensive and sensitive kidney-specific DIA-MS assays feature high-throughput analytical capabilities to achieve deep coverage of the kidney proteome, and will serve as useful tools for developing novel therapeutics to remediate kidney function.

Longitudinal profiling of the microbiome at four body sites reveals core stability and individualized dynamics during health and disease

To understand dynamic interplay between the human microbiome and host during health and disease, we analyzed the microbial composition, temporal dynamics, and associations with host multi-omics, immune and clinical markers of microbiomes from four body sites in 86 participants over six years. We found that microbiome stability and individuality are body-site-specific and heavily influenced by the host. The stool and oral microbiome were more stable than the skin and nasal microbiomes, possibly due to their interaction with the host and environment. Also, we identified individual-specific and commonly shared bacterial taxa, with individualized taxa showing greater stability. Interestingly, microbiome dynamics correlated across body sites, suggesting systemic coordination influenced by host-microbial-environment interactions. Notably, insulin-resistant individuals showed altered microbial stability and associations between microbiome, molecular markers, and clinical features, suggesting their disrupted interaction in metabolic disease. Our study offers comprehensive views of multi-site microbial dynamics and their relationship with host health and disease.

Study Highlights

The stability of the human microbiome varies among individuals and body sites.Highly individualized microbial genera are more stable over time.At each of the four body sites, systematic interactions between the environment, the host and bacteria can be detected.Individuals with insulin resistance have lower microbiome stability, a more diversified skin microbiome, and significantly altered host-microbiome interactions.

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

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

Data-Driven Tool for Cross-Run Ion Selection and Peak-Picking in Quantitative Proteomics with Data-Independent Acquisition LC-MS/MS

Proteomics provides molecular bases of biology and disease, and liquid chromatography-tandem mass spectrometry (LC-MS/MS) is a platform widely used for bottom-up proteomics. Data-independent acquisition (DIA) improves the run-to-run reproducibility of LC-MS/MS in proteomics research. However, the existing DIA data processing tools sometimes produce large deviations from true values for the peptides and proteins in quantification. Peak-picking error and incorrect ion selection are the two main causes of the deviations. We present a cross-run ion selection and peak-picking (CRISP) tool that utilizes the important advantage of run-to-run consistency of DIA and simultaneously examines the DIA data from the whole set of runs to filter out the interfering signals, instead of only looking at a single run at a time. Eight datasets acquired by mass spectrometers from different vendors with different types of mass analyzers were used to benchmark our CRISP-DIA against other currently available DIA tools. In the benchmark datasets, for analytes with large content variation among samples, CRISP-DIA generally resulted in 20 to 50% relative decrease in error rates compared to other DIA tools, at both the peptide precursor level and the protein level. CRISP-DIA detected differentially expressed proteins more efficiently, with 3.3 to 90.3% increases in the numbers of true positives and 12.3 to 35.3% decreases in the false positive rates, in some cases. In the real biological datasets, CRISP-DIA showed better consistencies of the quantification results. The advantages of assimilating DIA data in multiple runs for quantitative proteomics were demonstrated, which can significantly improve the quantification accuracy.

Advances in data-independent acquisition mass spectrometry towards comprehensive digital proteome landscape

The data-independent acquisition mass spectrometry (DIA-MS) has rapidly evolved as a powerful alternative for highly reproducible proteome profiling with a unique strength of generating permanent digital maps for retrospective analysis of biological systems. Recent advancements in data analysis software tools for the complex DIA-MS/MS spectra coupled to fast MS scanning speed and high mass accuracy have greatly expanded the sensitivity and coverage of DIA-based proteomics profiling. Here, we review the evolution of the DIA-MS techniques, from earlier proof-of-principle of parallel fragmentation of all-ions or ions in selected m/z range, the sequential window acquisition of all theoretical mass spectra (SWATH-MS) to latest innovations, recent development in computation algorithms for data informatics, and auxiliary tools and advanced instrumentation to enhance the performance of DIA-MS. We further summarize recent applications of DIA-MS and experimentally-derived as well as in silico spectra library resources for large-scale profiling to facilitate biomarker discovery and drug development in human diseases with emphasis on the proteomic profiling coverage. Toward next-generation DIA-MS for clinical proteomics, we outline the challenges in processing multi-dimensional DIA data set and large-scale clinical proteomics, and continuing need in higher profiling coverage and sensitivity.

Integration of data-independent acquisition (DIA) with co-fractionation mass spectrometry (CF-MS) to enhance interactome mapping capabilities

Proteomics technologies are continually advancing, providing opportunities to develop stronger and more robust protein interaction networks (PINs). In part, this is due to the ever-growing number of high-throughput proteomics methods that are available. This review discusses how data-independent acquisition (DIA) and co-fractionation mass spectrometry (CF-MS) can be integrated to enhance interactome mapping abilities. Furthermore, integrating these two techniques can improve data quality and network generation through extended protein coverage, less missing data, and reduced noise. CF-DIA-MS shows promise in expanding our knowledge of interactomes, notably for non-model organisms (NMOs). CF-MS is a valuable technique on its own, but upon the integration of DIA, the potential to develop robust PINs increases, offering a unique approach for researchers to gain an in-depth understanding into the dynamics of numerous biological processes.

Achieving quantitative reproducibility in label-free multisite DIA experiments through multirun alignment

DIA is a mainstream method for quantitative proteomics, but consistent quantification across multiple LC-MS/MS instruments remains a bottleneck in parallelizing data acquisition. One reason for this inconsistency and missing quantification is the retention time shift which current software does not adequately address for runs from multiple sites. We present multirun chromatogram alignment strategies to map peaks across columns, including the traditional reference-based Star method, and two novel approaches: MST and Progressive alignment. These reference-free strategies produce a quantitatively accurate data-matrix, even from heterogeneous multi-column studies. Progressive alignment also generates merged chromatograms from all runs which has not been previously achieved for LC-MS/MS data. First, we demonstrate the effectiveness of multirun alignment strategies on a gold-standard annotated dataset, resulting in a threefold reduction in quantitation error-rate compared to non-aligned DIA results. Subsequently, on a multi-species dataset that DIAlignR effectively controls the quantitative error rate, improves precision in protein measurements, and exhibits conservative peak alignment. We next show that the MST alignment reduces cross-site CV by 50% for highly abundant proteins when applied to a dataset from 11 different LC-MS/MS setups. Finally, the reanalysis of 949 plasma runs with multirun alignment revealed a more than 50% increase in insulin resistance (IR) and respiratory viral infection (RVI) proteins, identifying 11 and 13 proteins respectively, compared to prior analysis without it. The three strategies are implemented in our DIAlignR workflow (>2.3) and can be combined with linear, non-linear, or hybrid pairwise alignment.

MicroRNA and Protein Cargos of Human Limbal Epithelial Cell-Derived Exosomes and Their Regulatory Roles in Limbal Stromal Cells of Diabetic and Non-Diabetic Corneas

Epithelial and stromal/mesenchymal limbal stem cells contribute to corneal homeostasis and cell renewal. Extracellular vesicles (EVs), including exosomes (Exos), can be paracrine mediators of intercellular communication. Previously, we described cargos and regulatory roles of limbal stromal cell (LSC)-derived Exos in non-diabetic (N) and diabetic (DM) limbal epithelial cells (LECs). Presently, we quantify the miRNA and proteome profiles of human LEC-derived Exos and their regulatory roles in N- and DM-LSC. We revealed some miRNA and protein differences in DM vs. N-LEC-derived Exos' cargos, including proteins involved in Exo biogenesis and packaging that may affect Exo production and ultimately cellular crosstalk and corneal function. Treatment by N-Exos, but not by DM-Exos, enhanced wound healing in cultured N-LSCs and increased proliferation rates in N and DM LSCs vs. corresponding untreated (control) cells. N-Exos-treated LSCs reduced the keratocyte markers ALDH3A1 and lumican and increased the MSC markers CD73, CD90, and CD105 vs. control LSCs. These being opposite to the changes quantified in wounded LSCs. Overall, N-LEC Exos have a more pronounced effect on LSC wound healing, proliferation, and stem cell marker expression than DM-LEC Exos. This suggests that regulatory miRNA and protein cargo differences in DM- vs. N-LEC-derived Exos could contribute to the disease state.

Proteome-wide copy-number estimation from transcriptomics

Protein copy numbers constrain systems-level properties of regulatory networks, but absolute proteomic data remain scarce compared to transcriptomics obtained by RNA sequencing. We addressed this persistent gap by relating mRNA to protein statistically using best-available data from quantitative proteomics-transcriptomics for 4366 genes in 369 cell lines. The approach starts with a central estimate of protein copy number and hierarchically appends mRNA-protein and mRNA-mRNA dependencies to define an optimal gene-specific model that links mRNAs to protein. For dozens of independent cell lines and primary prostate samples, these protein inferences from mRNA outmatch stringent null models, a count-based protein-abundance repository, and empirical protein-to-mRNA ratios. The optimal mRNA-to-protein relationships capture biological processes along with hundreds of known protein-protein interaction complexes, suggesting mechanistic relationships are embedded. We use the method to estimate viral-receptor abundances of CD55-CXADR from human heart transcriptomes and build 1489 systems-biology models of coxsackievirus B3 infection susceptibility. When applied to 796 RNA sequencing profiles of breast cancer from The Cancer Genome Atlas, inferred copy-number estimates collectively reclassify 26% of Luminal A and 29% of Luminal B tumors. Protein-based reassignments strongly involve a pharmacologic target for luminal breast cancer (CDK4) and an α-catenin that is often undetectable at the mRNA level (CTTNA2). Thus, by adopting a gene-centered perspective of mRNA-protein covariation across different biological contexts, we achieve accuracies comparable to the technical reproducibility limits of contemporary proteomics. The collection of gene-specific models is assembled as a web tool for users seeking mRNA-guided predictions of absolute protein abundance (http://janeslab.shinyapps.io/Pinferna).

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.

Phenotypic Characterization of Acoustically Enriched Extracellular Vesicles from Pathogen-Activated Platelets

Extracellular vesicles (EVs) are derived from the membrane of platelets and released into the circulation upon activation or injury. Analogous to the parent cell, platelet-derived EVs play an important role in hemostasis and immune responses by transfer of bioactive cargo from the parent cells. Platelet activation and release of EVs increase in several pathological inflammatory diseases, such as sepsis. We have previously reported that the M1 protein released from the bacterial pathogen Streptococcus pyogenes directly mediates platelet activation. In this study, EVs were isolated from these pathogen-activated platelets using acoustic trapping, and their inflammation phenotype was characterized using quantitative mass spectrometry-based proteomics and cell-based models of inflammation. We determined that M1 protein mediated release of platelet-derived EVs that contained the M1 protein. The isolated EVs derived from pathogen-activated platelets contained a similar protein cargo to those from physiologically activated platelets (thrombin) and included platelet membrane proteins, granule proteins, cytoskeletal proteins, coagulation factors, and immune mediators. Immunomodulatory cargo, complement proteins, and IgG3 were significantly enriched in EVs isolated from M1 protein-stimulated platelets. Acoustically enriched EVs were functionally intact and exhibited pro-inflammatory effects on addition to blood, including platelet-neutrophil complex formation, neutrophil activation, and cytokine release. Collectively, our findings reveal novel aspects of pathogen-mediated platelet activation during invasive streptococcal infection.

Network-centric analysis of co-fractionated protein complex profiles using SECAT

The Size-Exclusion Chromatography Analysis Toolkit (SECAT) elucidates protein complex dynamics using co-fractionated bottom-up mass spectrometry (CF-MS) data. Here, we present a protocol for the network-centric analysis and interpretation of CF-MS profiles using SECAT. We describe the technical steps for preprocessing, scoring, semi-supervised machine learning, and quantification, including common pitfalls and their solutions. We further provide guidance for data export, visualization, and the interpretation of SECAT results to discover dysregulated proteins and interactions, supporting new hypotheses and biological insights. For complete details on the use and execution of this protocol, please refer to Rosenberger et al. (2020).¹.

M. tuberculosis relies on trace oxygen to maintain energy homeostasis and survive in hypoxic environments

The bioenergetic mechanisms by which Mycobacterium tuberculosis survives hypoxia are poorly understood. Current models assume that the bacterium shifts to an alternate electron acceptor or fermentation to maintain membrane potential and ATP synthesis. Counterintuitively, we find here that oxygen itself is the principal terminal electron acceptor during hypoxic dormancy. M. tuberculosis can metabolize oxygen efficiently at least two orders of magnitude below the concentration predicted to occur in hypoxic lung granulomas. Despite a difference in apparent affinity for oxygen, both the cytochrome bcc:aa₃ and cytochrome bd oxidase respiratory branches are required for hypoxic respiration. Simultaneous inhibition of both oxidases blocks oxygen consumption, reduces ATP levels, and kills M. tuberculosis under hypoxia. The capacity of mycobacteria to scavenge trace levels of oxygen, coupled with the absence of complex regulatory mechanisms to achieve hierarchal control of the terminal oxidases, may be a key determinant of long-term M. tuberculosis survival in hypoxic lung granulomas.

Strategies for Increasing the Depth and Throughput of Protein Analysis by plexDIA

Accurate protein quantification is key to identifying protein markers, regulatory relationships between proteins, and pathophysiological mechanisms. Realizing this potential requires sensitive and deep protein analysis of a large number of samples. Toward this goal, proteomics throughput can be increased by parallelizing the analysis of both precursors and samples using multiplexed data independent acquisition (DIA) implemented by the plexDIA framework: https://plexDIA.slavovlab.net. Here we demonstrate the improved precisions of retention time estimates within plexDIA and how this enables more accurate protein quantification. plexDIA has demonstrated multiplicative gains in throughput, and these gains may be substantially amplified by improving the multiplexing reagents, data acquisition, and interpretation. We discuss future directions for advancing plexDIA, which include engineering optimized mass-tags for high-plexDIA, introducing isotopologous carriers, and developing algorithms that utilize the regular structures of plexDIA data to improve sensitivity, proteome coverage, and quantitative accuracy. These advances in plexDIA will increase the throughput of functional proteomic assays, including quantifying protein conformations, turnover dynamics, modifications states and activities. The sensitivity of these assays will extend to single-cell analysis, thus enabling functional single-cell protein analysis.

Substantial Downregulation of Mitochondrial and Peroxisomal Proteins during Acute Kidney Injury revealed by Data-Independent Acquisition Proteomics

Acute kidney injury (AKI) manifests as a major health concern, particularly for the elderly. Understanding AKI-related proteome changes is critical for prevention and development of novel therapeutics to recover kidney function and to mitigate the susceptibility for recurrent AKI or development of chronic kidney disease. In this study, mouse kidneys were subjected to ischemia-reperfusion injury, and the contralateral kidneys remained uninjured to enable comparison and assess injury-induced changes in the kidney proteome. A fast-acquisition rate ZenoTOF 7600 mass spectrometer was introduced for data-independent acquisition (DIA) for comprehensive protein identification and quantification. Short microflow gradients and the generation of a deep kidney-specific spectral library allowed for high-throughput, comprehensive protein quantification. Upon AKI, the kidney proteome was completely remodeled, and over half of the 3,945 quantified protein groups changed significantly. Downregulated proteins in the injured kidney were involved in energy production, including numerous peroxisomal matrix proteins that function in fatty acid oxidation, such as ACOX1, CAT, EHHADH, ACOT4, ACOT8, and Scp2. Injured mice exhibited severely declined health. The comprehensive and sensitive kidney-specific DIA assays highlighted here feature high-throughput analytical capabilities to achieve deep coverage of the kidney proteome and will serve as useful tools for developing novel therapeutics to remediate kidney function.

Network-based elucidation of colon cancer drug resistance by phosphoproteomic time-series analysis

Aberrant signaling pathway activity is a hallmark of tumorigenesis and progression, which has guided targeted inhibitor design for over 30 years. Yet, adaptive resistance mechanisms, induced by rapid, context-specific signaling network rewiring, continue to challenge therapeutic efficacy. By leveraging progress in proteomic technologies and network-based methodologies, over the past decade, we developed VESPA-an algorithm designed to elucidate mechanisms of cell response and adaptation to drug perturbations-and used it to analyze 7-point phosphoproteomic time series from colorectal cancer cells treated with clinically-relevant inhibitors and control media. Interrogation of tumor-specific enzyme/substrate interactions accurately inferred kinase and phosphatase activity, based on their inferred substrate phosphorylation state, effectively accounting for signal cross-talk and sparse phosphoproteome coverage. The analysis elucidated time-dependent signaling pathway response to each drug perturbation and, more importantly, cell adaptive response and rewiring that was experimentally confirmed by CRISPRko assays, suggesting broad applicability to cancer and other diseases.

Profiling Serum Intact N-Glycopeptides Using Data-Independent Acquisition Mass Spectrometry

Large-scale profiling of intact glycopeptides is critical but challenging in glycoproteomics. Data-independent acquisition (DIA) mass spectrometry is an emerging technology with deep proteome coverage as well as accurate quantitative capability for large-scale proteomics studies and has also been applied to the field of glycoproteomics. In this protocol, we describe how to analyze data from a DIA experiment for profiling serum intact N-glycopeptides. We present a comprehensive data analysis workflow using GproDIA, including glycopeptide spectral library building, chromatographic feature extraction from the DIA data, and feature scoring with appropriate statistical control of error rates. We anticipate that this method could provide a powerful tool to explore the serum glycoproteome.

Meiotic nuclear pore complex remodeling provides key insights into nuclear basket organization

Nuclear pore complexes (NPCs) are large proteinaceous assemblies that mediate nuclear compartmentalization. NPCs undergo large-scale structural rearrangements during mitosis in metazoans and some fungi. However, our understanding of NPC remodeling beyond mitosis remains limited. Using time-lapse fluorescence microscopy, we discovered that NPCs undergo two mechanistically separable remodeling events during budding yeast meiosis in which parts or all of the nuclear basket transiently dissociate from the NPC core during meiosis I and II, respectively. Meiosis I detachment, observed for Nup60 and Nup2, is driven by Polo kinase-mediated phosphorylation of Nup60 at its interface with the Y-complex. Subsequent reattachment of Nup60-Nup2 to the NPC core is facilitated by a lipid-binding amphipathic helix in Nup60. Preventing Nup60-Nup2 reattachment causes misorganization of the entire nuclear basket in gametes. Strikingly, meiotic nuclear basket remodeling also occurs in the distantly related fission yeast, Schizosaccharomyces pombe. Our study reveals a conserved and developmentally programmed aspect of NPC plasticity, providing key mechanistic insights into the nuclear basket organization.

BIRCH: An Automated Workflow for Evaluation, Correction, and Visualization of Batch Effect in Bottom-Up Mass Spectrometry-Based Proteomics Data

Recent surges in large-scale mass spectrometry (MS)-based proteomics studies demand a concurrent rise in methods to facilitate reliable and reproducible data analysis. Quantification of proteins in MS analysis can be affected by variations in technical factors such as sample preparation and data acquisition conditions leading to batch effects, which adds to noise in the data set. This may in turn affect the effectiveness of any biological conclusions derived from the data. Here we present Batch-effect Identification, Representation, and Correction of Heterogeneous data (BIRCH), a workflow for analysis and correction of batch effect through an automated, versatile, and easy to use web-based tool with the goal of eliminating technical variation. BIRCH also supports diagnosis of the data to check for the presence of batch effects, feasibility of batch correction, and imputation to deal with missing values in the data set. To illustrate the relevance of the tool, we explore two case studies, including an iPSC-derived cell study and a Covid vaccine study to show different context-specific use cases. Ultimately this tool can be used as an extremely powerful approach for eliminating technical bias while retaining biological bias, toward understanding disease mechanisms and potential therapeutics.

Blood proteome of acute intracranial hemorrhage in infant victims of abusive head trauma

Abusive head trauma (AHT) is a leading cause of mortality and morbidity in infants. While the reported incidence is close to 40 cases per 100'000 births/year, misdiagnoses are commonly observed in cases with atypical, subacute, or chronic presentation. Currently, standard clinical evaluation of inflicted intracranial hemorrhagic injury (ICH) in infants urgently requires a screening test able to identify infants who need additional investigations. Blood biomarkers characteristic of AHT may assist in detecting these infants, improving prognosis through early medical care. To date, the application of innovative omics technologies in retrospective studies of AHT in infants is rare, due also to the blood serum and cerebrospinal fluid of AHT cases being scarce and not systematically accessible. Here, we explored the circulating blood proteomes of infants with severe AHT and their atraumatic controls. We discovered 165 circulating serum proteins that display differential changes in AHT cases compared with atraumatic controls. The peripheral blood proteomes of pediatric AHT commonly reflect: (i) potentially secreted proteome from injured brain, and (ii) proteome dysregulated in the system's circulation by successive biological events following acute ICH. This study opens up a novel opportunity for research efforts in clinical screening of AHT cases.

Hyperphosphorylation of hepatic proteome characterizes nonalcoholic fatty liver disease in S-adenosylmethionine deficiency

Methionine adenosyltransferase 1a (MAT1A) is responsible for hepatic S-adenosyl-L-methionine (SAMe) biosynthesis. Mat1a ^-/- mice have hepatic SAMe depletion, develop nonalcoholic steatohepatitis (NASH) which is reversed with SAMe administration. We examined temporal alterations in the proteome/phosphoproteome in pre-disease and NASH Mat1a ^-/- mice, effects of SAMe administration, and compared to human nonalcoholic fatty liver disease (NAFLD). Mitochondrial and peroxisomal lipid metabolism proteins were altered in pre-disease mice and persisted in NASH Mat1a ^-/- mice, which exhibited more progressive alterations in cytoplasmic ribosomes, ER, and nuclear proteins. A common mechanism found in both pre-disease and NASH livers was a hyperphosphorylation signature consistent with casein kinase 2α (CK2α) and AKT1 activation, which was normalized by SAMe administration. This was mimicked in human NAFLD with a metabolomic signature (M-subtype) resembling Mat1a ^-/- mice. In conclusion, we have identified a common proteome/phosphoproteome signature between Mat1a ^-/- mice and human NAFLD M-subtype that may have pathophysiological and therapeutic implications.

Introducing untargeted data-independent acquisition for metaproteomics of complex microbial samples

Mass spectrometry-based metaproteomics is a relatively new field of research that enables the characterization of the functionality of microbiota. Recently, we demonstrated the applicability of data-independent acquisition (DIA) mass spectrometry to the analysis of complex metaproteomic samples. This allowed us to circumvent many of the drawbacks of the previously used data-dependent acquisition (DDA) mass spectrometry, mainly the limited reproducibility when analyzing samples with complex microbial composition. However, the DDA-assisted DIA approach still required additional DDA data on the samples to assist the analysis. Here, we introduce, for the first time, an untargeted DIA metaproteomics tool that does not require any DDA data, but instead generates a pseudospectral library directly from the DIA data. This reduces the amount of required mass spectrometry data to a single DIA run per sample. The new DIA-only metaproteomics approach is implemented as a new open-source software package named glaDIAtor, including a modern web-based graphical user interface to facilitate wide use of the tool by the community.

Implementing the reuse of public DIA proteomics datasets: from the PRIDE database to Expression Atlas

The number of mass spectrometry (MS)-based proteomics datasets in the public domain keeps increasing, particularly those generated by Data Independent Acquisition (DIA) approaches such as SWATH-MS. Unlike Data Dependent Acquisition datasets, the re-use of DIA datasets has been rather limited to date, despite its high potential, due to the technical challenges involved. We introduce a (re-)analysis pipeline for public SWATH-MS datasets which includes a combination of metadata annotation protocols, automated workflows for MS data analysis, statistical analysis, and the integration of the results into the Expression Atlas resource. Automation is orchestrated with Nextflow, using containerised open analysis software tools, rendering the pipeline readily available and reproducible. To demonstrate its utility, we reanalysed 10 public DIA datasets from the PRIDE database, comprising 1,278 SWATH-MS runs. The robustness of the analysis was evaluated, and the results compared to those obtained in the original publications. The final expression values were integrated into Expression Atlas, making SWATH-MS experiments more widely available and combining them with expression data originating from other proteomics and transcriptomics datasets.

Global, distinctive, and personal changes in molecular and microbial profiles by specific fibers in humans

Dietary fibers act through the microbiome to improve cardiovascular health and prevent metabolic disorders and cancer. To understand the health benefits of dietary fiber supplementation, we investigated two popular purified fibers, arabinoxylan (AX) and long-chain inulin (LCI), and a mixture of five fibers. We present multiomic signatures of metabolomics, lipidomics, proteomics, metagenomics, a cytokine panel, and clinical measurements on healthy and insulin-resistant participants. Each fiber is associated with fiber-dependent biochemical and microbial responses. AX consumption associates with a significant reduction in LDL and an increase in bile acids, contributing to its observed cholesterol reduction. LCI is associated with an increase in Bifidobacterium. However, at the highest LCI dose, there is increased inflammation and elevation in the liver enzyme alanine aminotransferase. This study yields insights into the effects of fiber supplementation and the mechanisms behind fiber-induced cholesterol reduction, and it shows effects of individual, purified fibers on the microbiome.

Precision environmental health monitoring by longitudinal exposome and multi-omics profiling

Conventional environmental health studies have primarily focused on limited environmental stressors at the population level, which lacks the power to dissect the complexity and heterogeneity of individualized environmental exposures. Here, as a pilot case study, we integrated deep-profiled longitudinal personal exposome and internal multi-omics to systematically investigate how the exposome shapes a single individual's phenome. We annotated thousands of chemical and biological components in the personal exposome cloud and found they were significantly correlated with thousands of internal biomolecules, which was further cross-validated using corresponding clinical data. Our results showed that agrochemicals and fungi predominated in the highly diverse and dynamic personal exposome, and the biomolecules and pathways related to the individual's immune system, kidney, and liver were highly associated with the personal external exposome. Overall, this data-driven longitudinal monitoring study shows the potential dynamic interactions between the personal exposome and internal multi-omics, as well as the impact of the exposome on precision health by producing abundant testable hypotheses.

The impact of genomic variation on protein phosphorylation states and regulatory networks

Genomic variation impacts on cellular networks by affecting the abundance (e.g., protein levels) and the functional states (e.g., protein phosphorylation) of their components. Previous work has focused on the former, while in this context, the functional states of proteins have largely remained neglected. Here, we generated high‐quality transcriptome, proteome, and phosphoproteome data for a panel of 112 genomically well‐defined yeast strains. Genetic effects on transcripts were generally transmitted to the protein layer, but specific gene groups, such as ribosomal proteins, showed diverging effects on protein levels compared with RNA levels. Phosphorylation states proved crucial to unravel genetic effects on signaling networks. Correspondingly, genetic variants that cause phosphorylation changes were mostly different from those causing abundance changes in the respective proteins. Underscoring their relevance for cell physiology, phosphorylation traits were more strongly correlated with cell physiological traits such as chemical compound resistance or cell morphology, compared with transcript or protein abundance. This study demonstrates how molecular networks mediate the effects of genomic variants to cellular traits and highlights the particular importance of protein phosphorylation.

Benchmarking of analysis strategies for data-independent acquisition proteomics using a large-scale dataset comprising inter-patient heterogeneity

Numerous software tools exist for data-independent acquisition (DIA) analysis of clinical samples, necessitating their comprehensive benchmarking. We present a benchmark dataset comprising real-world inter-patient heterogeneity, which we use for in-depth benchmarking of DIA data analysis workflows for clinical settings. Combining spectral libraries, DIA software, sparsity reduction, normalization, and statistical tests results in 1428 distinct data analysis workflows, which we evaluate based on their ability to correctly identify differentially abundant proteins. From our dataset, we derive bootstrap datasets of varying sample sizes and use the whole range of bootstrap datasets to robustly evaluate each workflow. We find that all DIA software suites benefit from using a gas-phase fractionated spectral library, irrespective of the library refinement used. Gas-phase fractionation-based libraries perform best against two out of three reference protein lists. Among all investigated statistical tests non-parametric permutation-based statistical tests consistently perform best.

Data independent acquisition (DIA) has been gaining momentum in clinical proteomics. Here, the authors create a benchmark dataset comprising inter-patient heterogeneity to compare popular DIA data analysis workflows for identifying differentially abundant proteins.

Interaction of Bartonella henselae with Fibronectin Represents the Molecular Basis for Adhesion to Host Cells

Deciphering the mechanisms of bacterial host cell adhesion is a clue for preventing infections. We describe the underestimated role that the extracellular matrix protein fibronectin plays in the adhesion of human-pathogenic Bartonella henselae to host cells.

Proteomic analysis of descending thoracic aorta identifies unique and universal signatures of aneurysm and dissection

Objective

Methods

Results

Conclusions

Diseases of the descending thoracic aorta such as aneurysms and dissections carry a high degree of morbidity and mortality. At present, a complete understanding is still lacking of the genetics that drive these diseases and why some aortic segments dissect in the presence or absence of an aneurysm. We compared and contrasted the whole proteome expression of descending aortas from patients with normal, dissected, aneurysmal, and aneurysmal with dissected pathology aortic tissue. We uncovered potential tissue markers that might serve as future targets for therapy or predictors of disease progression.

Distinct serotypes of streptococcal M proteins mediate fibrinogen-dependent platelet activation and pro-inflammatory effects

Sepsis is a life-threatening complication of infection that is characterized by a dysregulated inflammatory state and disturbed hemostasis. Platelets are the main regulators of hemostasis, and they also respond to inflammation. The human pathogen Streptococcus pyogenes can cause local infection that may progress to sepsis. There are more than 200 serotypes of S. pyogenes defined according to sequence variations in the M protein. The M1 serotype is among 10 serotypes that are predominant in invasive infection. M1 protein can be released from the surface and has previously been shown to generate platelet, neutrophil, and monocyte activation. The platelet-dependent proinflammatory effects of other serotypes of M protein associated with invasive infection (M3, M5, M28, M49, and M89) are now investigated using a combination of multiparameter flow cytometry, enzyme-linked immunosorbent assay (ELISA), aggregometry, and quantitative mass spectrometry. We demonstrate that only M1, M3, and M5 protein serotypes can bind fibrinogen in plasma and mediate fibrinogen- and IgG-dependent platelet activation and aggregation, release of granule proteins, upregulation of CD62P to the platelet surface, and complex formation with neutrophils and monocytes. Neutrophil and monocyte activation, determined as upregulation of surface CD11b, is also mediated by M1, M3, and M5 protein serotypes, while M28, M49, and M89 proteins failed to mediate activation of platelets or leukocytes. Collectively, our findings reveal novel aspects of the immunomodulatory role of fibrinogen acquisition and platelet activation during streptococcal infections.

Comparative Proteomic Analysis of HPV(+) Oropharyngeal Squamous Cell Carcinoma Recurrence

Deintensification therapy for human papillomavirus-related oropharyngeal squamous cell carcinoma (HPV(+) OPSCC) is under active investigation. An adaptive treatment approach based on molecular stratification could identify high-risk patients predisposed to recurrence and better select for appropriate treatment regimens. Collectively, 40 HPV(+) OPSCC FFPE samples (20 disease-free, 20 recurrent) were surveyed using mass spectrometry-based proteomic analysis via data-independent acquisition to obtain fold change and false discovery differences. Ten-year overall survival was 100.0 and 27.7% for HPV(+) disease-free and recurrent cohorts, respectively. Of 1414 quantified proteins, 77 demonstrated significant differential expression. Top enriched functional pathways included those involved in programmed cell death (73 proteins, p = 7.43 × 10^-30), apoptosis (73 proteins, p = 5.56 × 10^-9), β-catenin independent WNT signaling (47 proteins, p = 1.45 × 10^-15), and Rho GTPase signaling (69 proteins, p = 1.09 × 10^-5). PFN1 (p = 1.0 × 10^-3), RAD23B (p = 2.9 × 10^-4), LDHB (p = 1.0 × 10^-3), and HINT1 (p = 3.8 × 10^-3) pathways were significantly downregulated in the recurrent cohort. On functional validation via immunohistochemistry (IHC) staining, 46.9% (PFN1), 71.9% (RAD23B), 59.4% (LDHB), and 84.4% (HINT1) of cases were corroborated with mass spectrometry findings. Development of a multilateral molecular signature incorporating these targets may characterize high-risk disease, predict treatment response, and augment current management paradigms in head and neck cancer.

Quantitative proteomic analysis of human peripheral nerves from subjects with type 2 diabetes

AIMS

Diabetic peripheral neuropathy (DPN) is a common and severe complication to type 2 diabetes. The pathogenesis of DPN is not fully known, but several pathways and gene polymorphisms contributing to DPN are described. DPN can be studied using nerve biopsies, but studies on the proteome of the nerve itself, and its surrounding tissue as a whole, are lacking. Studies on the posterior interosseous nerve (PIN) have proposed PIN a useful indicator of DPN.

METHODS

A quantitative mass spectrometry-based proteomics analysis was made of peripheral nerves from age- and gender-matched living human male tissue donors; nine type 2 diabetes subjects, with decreased sural nerve action potentials indicating DPN, and six controls without type 2 diabetes, with normal electrophysiology results.

RESULTS

A total of 2617 proteins were identified. Linear regression was used to discover which proteins were differentially expressed between type 2 diabetes and controls. Only soft signals were found. Therefore, clustering of the 500 most variable proteins was made to find clusters of similar proteins in type 2 diabetes subjects and healthy controls.

CONCLUSIONS

This feasibility study shows, for the first time, that the use of quantitative mass spectrometry enables quantification of proteins from nerve biopsies from subjects with and without type 2 diabetes, which may aid in finding biomarkers of importance to DPN development.

Streptococcus pyogenes Forms Serotype- and Local Environment-Dependent Interspecies Protein Complexes

Streptococcus pyogenes is known to cause both mucosal and systemic infections in humans. In this study, we used a combination of quantitative and structural mass spectrometry techniques to determine the composition and structure of the interaction network formed between human plasma proteins and the surfaces of different S. pyogenes serotypes. Quantitative network analysis revealed that S. pyogenes forms serotype-specific interaction networks that are highly dependent on the domain arrangement of the surface-attached M protein. Subsequent structural mass spectrometry analysis and computational modeling of one of the M proteins, M28, revealed that the network structure changes across different host microenvironments. We report that M28 binds secretory IgA via two separate binding sites with high affinity in saliva. During vascular leakage mimicked by increasing plasma concentrations in saliva, the binding of secretory IgA was replaced by the binding of monomeric IgA and C4b-binding protein (C4BP). This indicates that an upsurge of C4BP in the local microenvironment due to damage to the mucosal membrane drives the binding of C4BP and monomeric IgA to M28. These results suggest that S. pyogenes has evolved to form microenvironment-dependent host-pathogen protein complexes to combat human immune surveillance during both mucosal and systemic infections.

IMPORTANCEStreptococcus pyogenes (group A Streptococcus [GAS]), is a human-specific Gram-positive bacterium. Each year, the bacterium affects 700 million people globally, leading to 160,000 deaths. The clinical manifestations of S. pyogenes are diverse, ranging from mild and common infections like tonsillitis and impetigo to life-threatening systemic conditions such as sepsis and necrotizing fasciitis. S. pyogenes expresses multiple virulence factors on its surface to localize and initiate infections in humans. Among all these expressed virulence factors, the M protein is the most important antigen. In this study, we perform an in-depth characterization of the human protein interactions formed around one of the foremost human pathogens. This strategy allowed us to decipher the protein interaction networks around different S. pyogenes strains on a global scale and to compare and visualize how such interactions are mediated by M proteins.

GproDIA enables data-independent acquisition glycoproteomics with comprehensive statistical control

Large-scale profiling of intact glycopeptides is critical but challenging in glycoproteomics. Data independent acquisition (DIA) is an emerging technology with deep proteome coverage and accurate quantitative capability in proteomics studies, but is still in the early stage of development in the field of glycoproteomics. We propose GproDIA, a framework for the proteome-wide characterization of intact glycopeptides from DIA data with comprehensive statistical control by a 2-dimentional false discovery rate approach and a glycoform inference algorithm, enabling accurate identification of intact glycopeptides using wide isolation windows. We further utilize a semi-empirical spectrum prediction strategy to expand the coverage of spectral libraries of glycopeptides. We benchmark our method for N-glycopeptide profiling on DIA data of yeast and human serum samples, demonstrating that DIA with GproDIA outperforms the data-dependent acquisition-based methods for glycoproteomics in terms of capacity and data completeness of identification, as well as accuracy and precision of quantification. We expect that this work can provide a powerful tool for glycoproteomic studies.

The autoimmune signature of hyperinflammatory multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) manifests as a severe and uncontrolled inflammatory response with multiorgan involvement, occurring weeks after SARS-CoV-2 infection. Here, we utilized proteomics, RNA sequencing, autoantibody arrays, and B cell receptor (BCR) repertoire analysis to characterize MIS-C immunopathogenesis and identify factors contributing to severe manifestations and intensive care unit admission. Inflammation markers, humoral immune responses, neutrophil activation, and complement and coagulation pathways were highly enriched in MIS-C patient serum, with a more hyperinflammatory profile in severe than in mild MIS-C cases. We identified a strong autoimmune signature in MIS-C, with autoantibodies targeted to both ubiquitously expressed and tissue-specific antigens, suggesting autoantigen release and excessive antigenic drive may result from systemic tissue damage. We further identified a cluster of patients with enhanced neutrophil responses as well as high anti-Spike IgG and autoantibody titers. BCR sequencing of these patients identified a strong imprint of antigenic drive with substantial BCR sequence connectivity and usage of autoimmunity-associated immunoglobulin heavy chain variable region (IGHV) genes. This cluster was linked to a TRBV11-2 expanded T cell receptor (TCR) repertoire, consistent with previous studies indicating a superantigen-driven pathogenic process. Overall, we identify a combination of pathogenic pathways that culminate in MIS-C and may inform treatment.

ANPELA: analysis and performance assessment of the label-free quantification workflow for metaproteomic studies

Label-free quantification (LFQ) with a specific and sequentially integrated workflow of acquisition technique, quantification tool and processing method has emerged as the popular technique employed in metaproteomic research to provide a comprehensive landscape of the adaptive response of microbes to external stimuli and their interactions with other organisms or host cells. The performance of a specific LFQ workflow is highly dependent on the studied data. Hence, it is essential to discover the most appropriate one for a specific data set. However, it is challenging to perform such discovery due to the large number of possible workflows and the multifaceted nature of the evaluation criteria. Herein, a web server ANPELA (https://idrblab.org/anpela/) was developed and validated as the first tool enabling performance assessment of whole LFQ workflow (collective assessment by five well-established criteria with distinct underlying theories), and it enabled the identification of the optimal LFQ workflow(s) by a comprehensive performance ranking. ANPELA not only automatically detects the diverse formats of data generated by all quantification tools but also provides the most complete set of processing methods among the available web servers and stand-alone tools. Systematic validation using metaproteomic benchmarks revealed ANPELA's capabilities in 1 discovering well-performing workflow(s), (2) enabling assessment from multiple perspectives and (3) validating LFQ accuracy using spiked proteins. ANPELA has a unique ability to evaluate the performance of whole LFQ workflow and enables the discovery of the optimal LFQs by the comprehensive performance ranking of all 560 workflows. Therefore, it has great potential for applications in metaproteomic and other studies requiring LFQ techniques, as many features are shared among proteomic studies.

Robust optimization of SWATH-MS workflow for human blood serum proteome analysis using a quality by design approach

Background

It is not enough to optimize proteomics assays. It is critical those assays are robust to operating conditions. Without robust assays, proteomic biomarkers are unlikely to translate readily into the clinic. This study outlines a structured approach to the identification of a robust operating window for proteomics assays and applies that method to Sequential Window Acquisition of all Theoretical Spectra Mass Spectroscopy (SWATH-MS).

Methods

We used a sequential quality by design approach exploiting a fractional screening design to first identify critical SWATH-MS parameters, then using response surface methods to identify a robust operating window with good reproducibility, before validating those settings in a separate validation study.

Results

The screening experiment identified two critical SWATH-MS parameters. We modelled the number of proteins and reproducibility as a function of those parameters identifying an operating window permitting robust maximization of the number of proteins quantified in human serum. In a separate validation study, these settings were shown to give good proteome-wide coverage and high quantification reproducibility.

Conclusions

Using design of experiments permits identification of a robust operating window for SWATH-MS. The method gives a good understanding of proteomics assays and greater data-driven confidence in SWATH-MS performance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12014-021-09323-z.

Diagnostics and correction of batch effects in large-scale proteomic studies: a tutorial

Advancements in mass spectrometry-based proteomics have enabled experiments encompassing hundreds of samples. While these large sample sets deliver much-needed statistical power, handling them introduces technical variability known as batch effects. Here, we present a step-by-step protocol for the assessment, normalization, and batch correction of proteomic data. We review established methodologies from related fields and describe solutions specific to proteomic challenges, such as ion intensity drift and missing values in quantitative feature matrices. Finally, we compile a set of techniques that enable control of batch effect adjustment quality. We provide an R package, "proBatch", containing functions required for each step of the protocol. We demonstrate the utility of this methodology on five proteomic datasets each encompassing hundreds of samples and consisting of multiple experimental designs. In conclusion, we provide guidelines and tools to make the extraction of true biological signal from large proteomic studies more robust and transparent, ultimately facilitating reliable and reproducible research in clinical proteomics and systems biology.

Alendronate-induced Perturbation of the Bone Proteome and Microenvironmental Pathophysiology

Objectives: Bisphosphonates (BPs) are powerful inhibitors of osteoclastogenesis and are used to prevent osteoporotic bone loss and reduce the risk of osteoporotic fracture in patients suffering from postmenopausal osteoporosis. Patients with breast cancer or gynecological malignancies being treated with BPs or those receiving bone-targeted therapy for metastatic prostate cancer are at increased risk of bisphosphonate-related osteonecrosis of the jaw (BRONJ). Although BPs markedly ameliorate osteoporosis, their adverse effects largely limit the clinical application of these drugs. This study focused on providing a deeper understanding of one of the most popular BPs, the alendronate (ALN)-induced perturbation of the bone proteome and microenvironmental pathophysiology. Methods: To understand the molecular mechanisms underlying ALN-induced side-effects, an unbiased and global proteomics approach combined with big data bioinformatics was applied. This was followed by biochemical and functional analyses to determine the clinicopathological mechanisms affected by ALN. Results: The findings from this proteomics study suggest that the RIPK3/Wnt/GSK3/β-catenin signaling pathway is significantly perturbed upon ALN treatment, resulting in abnormal angiogenesis, inflammation, anabolism, remodeling, and mineralization in bone cells in an in vitro cell culture system. Conclusion: Our investigation into potential key signaling mechanisms in response to ALN provides a rational basis for suppressing BP-induced adverse effect and presents various therapeutic strategies.

DIAproteomics: A Multifunctional Data Analysis Pipeline for Data-Independent Acquisition Proteomics and Peptidomics

Data-independent acquisition (DIA) is becoming a leading analysis method in biomedical mass spectrometry. The main advantages include greater reproducibility and sensitivity and a greater dynamic range compared with data-dependent acquisition (DDA). However, the data analysis is complex and often requires expert knowledge when dealing with large-scale data sets. Here we present DIAproteomics, a multifunctional, automated, high-throughput pipeline implemented in the Nextflow workflow management system that allows one to easily process proteomics and peptidomics DIA data sets on diverse compute infrastructures. The central components are well-established tools such as the OpenSwathWorkflow for the DIA spectral library search and PyProphet for the false discovery rate assessment. In addition, it provides options to generate spectral libraries from existing DDA data and to carry out the retention time and chromatogram alignment. The output includes annotated tables and diagnostic visualizations from the statistical postprocessing and computation of fold-changes across pairwise conditions, predefined in an experimental design. DIAproteomics is well documented open-source software and is available under a permissive license to the scientific community at https://www.openms.de/diaproteomics/.

Proteogenomic Workflow Reveals Molecular Phenotypes Related to Breast Cancer Mammographic Appearance

Proteogenomic approaches have enabled the generat̲ion of novel information levels when compared to single omics studies although burdened by extensive experimental efforts. Here, we improved a data-independent acquisition mass spectrometry proteogenomic workflow to reveal distinct molecular features related to mammographic appearances in breast cancer. Our results reveal splicing processes detectable at the protein level and highlight quantitation and pathway complementarity between RNA and protein data. Furthermore, we confirm previously detected enrichments of molecular pathways associated with estrogen receptor-dependent activity and provide novel evidence of epithelial-to-mesenchymal activity in mammography-detected spiculated tumors. Several transcript-protein pairs displayed radically different abundances depending on the overall clinical properties of the tumor. These results demonstrate that there are differentially regulated protein networks in clinically relevant tumor subgroups, which in turn alter both cancer biology and the abundance of biomarker candidates and drug targets.

A Simple Optimization Workflow to Enable Precise and Accurate Imputation of Missing Values in Proteomic Data Sets

Missing values in proteomic data sets have real consequences on downstream data analysis and reproducibility. Although several imputation methods exist to handle missing values, no single imputation method is best suited for a diverse range of data sets, and no clear strategy exists for evaluating imputation methods for clinical DIA-MS data sets, especially at different levels of protein quantification. To navigate through the different imputation strategies available in the literature, we have established a strategy to assess imputation methods on clinical label-free DIA-MS data sets. We used three DIA-MS data sets with real missing values to evaluate eight imputation methods with multiple parameters at different levels of protein quantification: a dilution series data set, a small pilot data set, and a clinical proteomic data set comparing paired tumor and stroma tissue. We found that imputation methods based on local structures within the data, like local least-squares (LLS) and random forest (RF), worked well in our dilution series data set, whereas imputation methods based on global structures within the data, like BPCA, performed well in the other two data sets. We also found that imputation at the most basic protein quantification level-fragment level-improved accuracy and the number of proteins quantified. With this analytical framework, we quickly and cost-effectively evaluated different imputation methods using two smaller complementary data sets to narrow down to the larger proteomic data set's most accurate methods. This acquisition strategy allowed us to provide reproducible evidence of the accuracy of the imputation method, even in the absence of a ground truth. Overall, this study indicates that the most suitable imputation method relies on the overall structure of the data set and provides an example of an analytic framework that may assist in identifying the most appropriate imputation strategies for the differential analysis of proteins.

From coarse to fine: the absolute Escherichia coli proteome under diverse growth conditions

Accurate measurements of cellular protein concentrations are invaluable to quantitative studies of gene expression and physiology in living cells. Here, we developed a versatile mass spectrometric workflow based on data-independent acquisition proteomics (DIA/SWATH) together with a novel protein inference algorithm (xTop). We used this workflow to accurately quantify absolute protein abundances in Escherichia coli for > 2,000 proteins over > 60 growth conditions, including nutrient limitations, non-metabolic stresses, and non-planktonic states. The resulting high-quality dataset of protein mass fractions allowed us to characterize proteome responses from a coarse (groups of related proteins) to a fine (individual) protein level. Hereby, a plethora of novel biological findings could be elucidated, including the generic upregulation of low-abundant proteins under various metabolic limitations, the non-specificity of catabolic enzymes upregulated under carbon limitation, the lack of large-scale proteome reallocation under stress compared to nutrient limitations, as well as surprising strain-dependent effects important for biofilm formation. These results present valuable resources for the systems biology community and can be used for future multi-omics studies of gene regulation and metabolic control in E. coli.