mapDIA: Preprocessing and statistical analysis of quantitative proteomics data from data independent acquisition mass spectrometry

High-Throughput Workflow for Detergent-free Cell-Based Proteomic Characterization

We have developed an automated cell-based workflow for the quantification of proteins by liquid chromatography-mass spectrometry (LC-MS) that facilitates large-scale perturbation studies carried out in a 96-well plate format and enables the preparation of one full plate in approximately 4 h, showcasing a high-throughput (HTP) concept. Cells were grown in a 96-well plate and lysed via ultrasonication. Proteins were subsequently solubilized, extracted, and processed into tryptic peptides for 2 h before being acquired by data-independent acquisition mass spectrometry (DIA-MS). This workflow leverages adaptive focused acoustics (AFA) technology for ultrasonication to aid cell lysis and protein solubilization on an automated liquid handling platform. As proof of principle, AC16 human cardiomyocyte-like cells were cultured in a 96-well plate under optimized conditions that were compatible with the downstream HTP pipeline. Over 30,000 peptides were identified, corresponding to the detection of 5100 unique proteins. 50% of measured proteins had an average coefficient of variation (CV) under 25% from approximately 30,000 cells. Our optimized detergent-free buffer consisting of ammonium bicarbonate yielded comparable findings. For the same number of cells, 5000 proteins were identified from 29,000 peptides, 40% of which demonstrated a CV under 25%.

Efficacy of IDOR-1117-2520, a novel, orally available CCR6 antagonist in preclinical models of skin dermatitis

BACKGROUND AND PURPOSE

The chemokine receptor CCR6 guides pathogenic T17 cells, implicated in autoimmune diseases including psoriasis, to sites of inflammation via the chemokine CCL20. Therefor, pharmacological inhibition of CCR6+ immune cell migration provides a novel therapeutic approach. Translatability of such an intervention has not yet been assessed in detail. We evaluated the translatability of the Aldara® mouse model induced skin inflammation to psoriasis, with particular focus on immune cell trafficking and assessed the efficacy of IDOR-1117-2520, a highly selective, potent and orally available CCR6 small inhibitor.

EXPERIMENTAL APPROACH

Effects of IDOR-1117-2520 were investigated in the Aldara® and IL23 mouse models of skin inflammation using flow cytometry, RNA sequencing and transcriptome-based cell type deconvolution approaches to characterise immune cell migration patterns. These results were compared to human psoriasis transcriptomics data.

KEY RESULTS

IDOR-1117-2520 dose dependently reduced infiltration of CCR6+ immune cells into inflamed skin, and was equally efficacious as IL-17 and IL-23 inhibition in models of skin inflammation. Pathway analysis showed molecular similarities in the immune response between human psoriasis and the Aldara® mouse model. IL-17/IL-23 pathway genes were expressed in both human psoriasis and the mouse model. CCR6 inhibition modulated multiple pathways associated with inflammation beyond the proximal IL-17/IL-23 pathway. A chemokine-chemokine receptor interaction map implicated CCL20-CCR6 as the dominant axis in recruiting pathogenic T17 cells in both the model and in human psoriasis.

CONCLUSION AND IMPLICATIONS

IDOR-1117-2520 could provide a promising novel targeted approach to treating psoriasis and, potentially, other autoimmune diseases involving the CCR6/CCL20 axis and the IL-17/IL-23 pathway. IDOR-1117-2520 is currently being evaluated in a clinical phase 1 trial (ISRCTN28892128).

Enhancing Proteomics Quality Control: Insights from the Visualization Tool QCeltis

Large-scale mass-spectrometry-based proteomics experiments are complex and prone to analytical variability, requiring rigorous quality checks across each step in the workflow: sample preparation, chromatography, mass spectrometry, and the bioinformatics stages. This includes quality control (QC) measures that address biological and technical variation. Most QC approaches involve detecting sample outliers and monitoring parameters related to sample preparation and mass spectrometer performance. Evaluating these parameters regularly is essential for reliable downstream analysis and proteomics research. Here, we introduce "QCeltis", a Python package designed to facilitate automated QC analysis across the proteomics workflow, aiding in the identification of technical biases and consistency verification. QCeltis is a versatile tool for detecting QC issues in large-scale data-independent acquisition proteomics experiments by not only identifying sample preparation and acquisition issues but also aiding in differentiating between QC issues vs batch effects. QCeltis is available for command-line use in Windows and Linux environments. We present three case studies showcasing QCeltis's capabilities across different data sets, including depleted plasma, whole blood vs plasma, and dried blood spot samples, emphasizing its potential impact on large-scale proteomics projects. This package can be used to enhance data reliability and enable nuanced downstream analysis and interpretation for proteomics studies.

Cross-Run Hybrid Features Improve the Identification of Data-Independent Acquisition Proteomics

The analysis of data-independent acquisition (DIA) mass spectrometry data is crucial for comprehensive proteomics studies. However, traditional single-run methods often fall short in terms of identification depth and consistency. We present HFDiscrim, a specialized multirun DIA analysis tool aimed at enhancing the depth and consistency of reliable peptide identifications of DIA analysis tools. HFDiscrim was extensively benchmarked on multiple data sets, including the MCB data set, the ccRCC data set, and a three-species benchmark mixture. Compared to PyProphet, HFDiscrim identified 22.04% more precursors, 19.1% more peptides, and 13.2% more proteins while maintaining a controllable false discovery rate. Furthermore, HFDiscrim demonstrated higher identification rates and improved reproducibility across multiple runs. HFDiscrim is publicly available at https://github.com/yachliu/HFDiscrim.

An Equine Protein Atlas Highlights Synovial Fluid Proteome Dynamics during Experimentally LPS-Induced Arthritis

In human proteomics, substantial efforts are ongoing to leverage large collections of mass spectrometry (MS) fragment ion spectra into extensive spectral libraries (SL) as a resource for data independent acquisition (DIA) analysis. Currently, such initiatives in equine research are still missing. Here we present a large-scale equine SL, comprising 6394 canonical proteins and 89,329 unique peptides, based on data dependent acquisition analysis of 75 tissue and body fluid samples from horses. The SL enabled large-scale DIA-MS based quantification of the same samples to generate a quantitative equine protein distribution atlas to infer dominant proteins in different organs and body fluids. Data mining revealed 163 proteins uniquely identified in a specific type of tissue or body fluid, serving as a starting point to determine tissue-specific or tissue-type-specific proteins. We showcase the SL by highlighting proteome dynamics in equine synovial fluid samples during experimental lipopolysaccharide-induced arthritis. A fuzzy c-means cluster analysis pinpointed SERPINB1, ATRN, NGAL, LTF, MMP1, and LBP as putative biomarkers for joint inflammation. This SL provides an extendable resource for future equine studies employing DIA-MS.

Analysis and Visualization of Quantitative Proteomics Data Using FragPipe-Analyst

The FragPipe computational proteomics platform is gaining widespread popularity among the proteomics research community because of its fast processing speed and user-friendly graphical interface. Although FragPipe produces well-formatted output tables that are ready for analysis, there is still a need for an easy-to-use and user-friendly downstream statistical analysis and visualization tool. FragPipe-Analyst addresses this need by providing an R shiny web server to assist FragPipe users in conducting downstream analyses of the resulting quantitative proteomics data. It supports major quantification workflows, including label-free quantification, tandem mass tags, and data-independent acquisition. FragPipe-Analyst offers a range of useful functionalities, such as various missing value imputation options, data quality control, unsupervised clustering, differential expression (DE) analysis using Limma, and gene ontology and pathway enrichment analysis using Enrichr. To support advanced analysis and customized visualizations, we also developed FragPipeAnalystR, an R package encompassing all FragPipe-Analyst functionalities that is extended to support site-specific analysis of post-translational modifications (PTMs). FragPipe-Analyst and FragPipeAnalystR are both open-source and freely available.

A proteomic approach supports the clinical relevance of TAT-Cx43266-283 in glioblastoma

Glioblastoma (GBM) is the most frequent and aggressive primary brain cancer. The Src inhibitor, TAT-Cx43266-283, exerts antitumor effects in in vitro and in vivo models of GBM. Because addressing the mechanism of action is essential to translate these results to a clinical setting, in this study we carried out an unbiased proteomic approach. Data-independent acquisition mass spectrometry proteomics allowed the identification of 190 proteins whose abundance was modified by TAT-Cx43266-283. Our results were consistent with the inhibition of Src as the mechanism of action of TAT-Cx43266-283 and unveiled antitumor effectors, such as p120 catenin. Changes in the abundance of several proteins suggested that TAT-Cx43266-283 may also impact the brain microenvironment. Importantly, the proteins whose abundance was reduced by TAT-Cx43266-283 correlated with an improved GBM patient survival in clinical datasets and none of the proteins whose abundance was increased by TAT-Cx43266-283 correlated with shorter survival, supporting its use in clinical trials.

Patient subtyping analysis of baseline multi-omic data reveals distinct pre-immune states associated with antibody response to seasonal influenza vaccination

Understanding the molecular mechanisms underpinning diverse vaccination responses is critical for developing efficient vaccines. Molecular subtyping can offer insights into heterogeneous nature of responses and aid in vaccine design. We analyzed multi-omic data from 62 haemagglutinin seasonal influenza vaccine recipients (2019-2020), including transcriptomics, proteomics, glycomics, and metabolomics data collected pre-vaccination. We performed a subtyping analysis on the integrated data revealing five subtypes with distinct molecular signatures. These subtypes differed in the expression of pre-existing adaptive or innate immunity signatures, which were linked to significant variation in baseline immunoglobulin A (IgA) and hemagglutination inhibition (HAI) titer levels. It is worth noting that these differences persisted through day 28 post-vaccination, indicating the effect of initial immune state on vaccination response. These findings highlight the significance of interpersonal variation in baseline immune status as a crucial factor in determining the effectiveness of seasonal vaccines. Ultimately, incorporating molecular profiling could enable personalized vaccine optimization.

Myocardial ultrastructure of human heart failure with preserved ejection fraction

Over half of patients with heart failure have a preserved ejection fraction (>50%, called HFpEF), a syndrome with substantial morbidity/mortality and few effective therapies1. Its dominant comorbidity is now obesity, which worsens disease and prognosis1-3. Myocardial data from patients with morbid obesity and HFpEF show depressed myocyte calcium-stimulated tension4 and disrupted gene expression of mitochondrial and lipid metabolic pathways5,6, abnormalities shared by human HF with a reduced EF but less so in HFpEF without severe obesity. The impact of severe obesity on human HFpEF myocardial ultrastructure remains unexplored. Here we assessed the myocardial ultrastructure in septal biopsies from patients with HFpEF using transmission electron microscopy. We observed sarcomere disruption and sarcolysis, mitochondrial swelling with cristae separation and dissolution and lipid droplet accumulation that was more prominent in the most obese patients with HFpEF and not dependent on comorbid diabetes. Myocardial proteomics revealed associated reduction in fatty acid uptake, processing and oxidation and mitochondrial respiration proteins, particularly in very obese patients with HFpEF.

Topical and oral peroxisome proliferator-activated receptor-α agonist ameliorates diabetic corneal neuropathy

Diabetic corneal neuropathy (DCN) is a common diabetic ocular complication with limited treatment options. In this study, we investigated the effects of topical and oral fenofibrate, a peroxisome proliferator-activated receptor-α agonist, on the amelioration of DCN using diabetic mice (n = 120). Ocular surface assessments, corneal nerve and cell imaging analysis, tear proteomics and its associated biological pathways, immuno-histochemistry and western blot on PPARα expression, were studied before and 12 weeks after treatment. At 12 weeks, PPARα expression markedly restored after topical and oral fenofibrate. Topical fenofibrate significantly improved corneal nerve fibre density (CNFD) and tortuosity coefficient. Likewise, oral fenofibrate significantly improved CNFD. Both topical and oral forms significantly improved corneal sensitivity. Additionally, topical and oral fenofibrate significantly alleviated diabetic keratopathy, with fenofibrate eye drops demonstrating earlier therapeutic effects. Both topical and oral fenofibrate significantly increased corneal β-III tubulin expression. Topical fenofibrate reduced neuroinflammation by significantly increasing the levels of nerve growth factor and substance P. It also significantly increased β-III-tubulin and reduced CDC42 mRNA expression in trigeminal ganglions. Proteomic analysis showed that neurotrophin signalling and anti-inflammation reactions were significantly up-regulated after fenofibrate treatment, whether applied topically or orally. This study concluded that both topical and oral fenofibrate ameliorate DCN, while topical fenofibrate significantly reduces neuroinflammation.

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.

A uniquely efficacious type of CFTR corrector with complementary mode of action

Three distinct pharmacological corrector types (I, II, III) with different binding sites and additive behavior only partially rescue the F508del-cystic fibrosis transmembrane conductance regulator (CFTR) folding and trafficking defect observed in cystic fibrosis. We describe uniquely effective, macrocyclic CFTR correctors that were additive to the known corrector types, exerting a complementary "type IV" corrector mechanism. Macrocycles achieved wild-type-like folding efficiency of F508del-CFTR at the endoplasmic reticulum and normalized CFTR currents in reconstituted patient-derived bronchial epithelium. Using photo-activatable macrocycles, docking studies and site-directed mutagenesis a highly probable binding site and pose for type IV correctors was identified in a cavity between lasso helix-1 (Lh1) and transmembrane helix-1 of membrane spanning domain (MSD)-1, distinct from the known corrector binding sites. Since only F508del-CFTR fragments spanning from Lh1 until MSD2 responded to type IV correctors, these likely promote cotranslational assembly of Lh1, MSD1, and MSD2. Previously corrector-resistant CFTR folding mutants were also robustly rescued, suggesting substantial therapeutic potential for type IV correctors.

Complex changes in serum protein levels in COVID-19 convalescents

The COVID-19 pandemic, triggered by severe acute respiratory syndrome coronavirus 2, has affected millions of people worldwide. Much research has been dedicated to our understanding of COVID-19 disease heterogeneity and severity, but less is known about recovery associated changes. To address this gap in knowledge, we quantified the proteome from serum samples from 29 COVID-19 convalescents and 29 age-, race-, and sex-matched healthy controls. Samples were acquired within the first months of the pandemic. Many proteins from pathways known to change during acute COVID-19 illness, such as from the complement cascade, coagulation system, inflammation and adaptive immune system, had returned to levels seen in healthy controls. In comparison, we identified 22 and 15 proteins with significantly elevated and lowered levels, respectively, amongst COVID-19 convalescents compared to healthy controls. Some of the changes were similar to those observed for the acute phase of the disease, i.e. elevated levels of proteins from hemolysis, the adaptive immune systems, and inflammation. In contrast, some alterations opposed those in the acute phase, e.g. elevated levels of CETP and APOA1 which function in lipid/cholesterol metabolism, and decreased levels of proteins from the complement cascade (e.g. C1R, C1S, and VWF), the coagulation system (e.g. THBS1 and VWF), and the regulation of the actin cytoskeleton (e.g. PFN1 and CFL1) amongst COVID-19 convalescents. We speculate that some of these shifts might originate from a transient decrease in platelet counts upon recovery from the disease. Finally, we observed race-specific changes, e.g. with respect to immunoglobulins and proteins related to cholesterol metabolism.

The Molecular Twin artificial-intelligence platform integrates multi-omic data to predict outcomes for pancreatic adenocarcinoma patients

Contemporary analyses focused on a limited number of clinical and molecular biomarkers have been unable to accurately predict clinical outcomes in pancreatic ductal adenocarcinoma. Here we describe a precision medicine platform known as the Molecular Twin consisting of advanced machine-learning models and use it to analyze a dataset of 6,363 clinical and multi-omic molecular features from patients with resected pancreatic ductal adenocarcinoma to accurately predict disease survival (DS). We show that a full multi-omic model predicts DS with the highest accuracy and that plasma protein is the top single-omic predictor of DS. A parsimonious model learning only 589 multi-omic features demonstrated similar predictive performance as the full multi-omic model. Our platform enables discovery of parsimonious biomarker panels and performance assessment of outcome prediction models learning from resource-intensive panels. This approach has considerable potential to impact clinical care and democratize precision cancer medicine worldwide.

Patient Subtyping Analysis of Baseline Multi-omic Data Reveals Distinct Pre-immune States Predictive of Vaccination Responses

Understanding the molecular mechanisms that underpin diverse vaccination responses is a critical step toward developing efficient vaccines. Molecular subtyping approaches can offer valuable insights into the heterogeneous nature of responses and aid in the design of more effective vaccines. In order to explore the molecular signatures associated with the vaccine response, we analyzed baseline transcriptomics data from paired samples of whole blood, proteomics and glycomics data from serum, and metabolomics data from urine, obtained from influenza vaccine recipients (2019-2020 season) prior to vaccination. After integrating the data using a network-based model, we performed a subtyping analysis. The integration of multiple data modalities from 62 samples resulted in five baseline molecular subtypes with distinct molecular signatures. These baseline subtypes differed in the expression of pre-existing adaptive or innate immunity signatures, which were linked to significant variation across subtypes in baseline immunoglobulin A (IgA) and hemagglutination inhibition (HAI) titer levels. It is worth noting that these significant differences persisted through day 28 post-vaccination, indicating the effect of initial immune state on vaccination response. These findings highlight the significance of interpersonal variation in baseline immune status as a crucial factor in determining vaccine response and efficacy. Ultimately, incorporating molecular profiling could enable personalized vaccine optimization.

Proteomic Profiling of Advanced Melanoma Patients to Predict Therapeutic Response to Anti-PD-1 Therapy

PURPOSE

Despite high clinical need, there are no biomarkers that accurately predict the response of patients with metastatic melanoma to anti-PD-1 therapy.

EXPERIMENTAL DESIGN

In this multicenter study, we applied protein depletion and enrichment methods prior to various proteomic techniques to analyze a serum discovery cohort (n = 56) and three independent serum validation cohorts (n = 80, n = 12, n = 17). Further validation analyses by literature and survival analysis followed.

RESULTS

We identified several significantly regulated proteins as well as biological processes such as neutrophil degranulation, cell-substrate adhesion, and extracellular matrix organization. Analysis of the three independent serum validation cohorts confirmed the significant differences between responders (R) and nonresponders (NR) observed in the initial discovery cohort. In addition, literature-based validation highlighted 30 markers overlapping with previously published signatures. Survival analysis using the TCGA database showed that overexpression of 17 of the markers we identified correlated with lower overall survival in patients with melanoma.

CONCLUSIONS

Ultimately, this multilayered serum analysis led to a potential marker signature with 10 key markers significantly altered in at least two independent serum cohorts: CRP, LYVE1, SAA2, C1RL, CFHR3, LBP, LDHB, S100A8, S100A9, and SAA1, which will serve as the basis for further investigation. In addition to patient serum, we analyzed primary melanoma tumor cells from NR and found a potential marker signature with four key markers: LAMC1, PXDN, SERPINE1, and VCAN.

Functional and spatial proteomics profiling reveals intra- and intercellular signaling crosstalk in colorectal cancer

Precision oncology approaches for patients with colorectal cancer (CRC) continue to lag behind other solid cancers. Functional precision oncology-a strategy that is based on perturbing primary tumor cells from cancer patients-could provide a road forward to personalize treatment. We extend this paradigm to measuring proteome activity landscapes by acquiring quantitative phosphoproteomic data from patient-derived organoids (PDOs). We show that kinase inhibitors induce inhibitor- and patient-specific off-target effects and pathway crosstalk. Reconstruction of the kinase networks revealed that the signaling rewiring is modestly affected by mutations. We show non-genetic heterogeneity of the PDOs and upregulation of stemness and differentiation genes by kinase inhibitors. Using imaging mass-cytometry-based profiling of the primary tumors, we characterize the tumor microenvironment (TME) and determine spatial heterocellular crosstalk and tumor-immune cell interactions. Collectively, we provide a framework for inferring tumor cell intrinsic signaling and external signaling from the TME to inform precision (immuno-) oncology in CRC.

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.

Analysis of DIA proteomics data using MSFragger-DIA and FragPipe computational platform

Liquid chromatography (LC) coupled with data-independent acquisition (DIA) mass spectrometry (MS) has been increasingly used in quantitative proteomics studies. Here, we present a fast and sensitive approach for direct peptide identification from DIA data, MSFragger-DIA, which leverages the unmatched speed of the fragment ion indexing-based search engine MSFragger. Different from most existing methods, MSFragger-DIA conducts a database search of the DIA tandem mass (MS/MS) spectra prior to spectral feature detection and peak tracing across the LC dimension. To streamline the analysis of DIA data and enable easy reproducibility, we integrate MSFragger-DIA into the FragPipe computational platform for seamless support of peptide identification and spectral library building from DIA, data-dependent acquisition (DDA), or both data types combined. We compare MSFragger-DIA with other DIA tools, such as DIA-Umpire based workflow in FragPipe, Spectronaut, DIA-NN library-free, and MaxDIA. We demonstrate the fast, sensitive, and accurate performance of MSFragger-DIA across a variety of sample types and data acquisition schemes, including single-cell proteomics, phosphoproteomics, and large-scale tumor proteome profiling studies.

Accurate Label-Free Quantification by directLFQ to Compare Unlimited Numbers of Proteomes

Recent advances in mass spectrometry-based proteomics enable the acquisition of increasingly large datasets within relatively short times, which exposes bottlenecks in the bioinformatics pipeline. Although peptide identification is already scalable, most label-free quantification (LFQ) algorithms scale quadratic or cubic with the sample numbers, which may even preclude the analysis of large-scale data. Here we introduce directLFQ, a ratio-based approach for sample normalization and the calculation of protein intensities. It estimates quantities via aligning samples and ion traces by shifting them on top of each other in logarithmic space. Importantly, directLFQ scales linearly with the number of samples, allowing analyses of large studies to finish in minutes instead of days or months. We quantify 10,000 proteomes in 10 min and 100,000 proteomes in less than 2 h, a 1000-fold faster than some implementations of the popular LFQ algorithm MaxLFQ. In-depth characterization of directLFQ reveals excellent normalization properties and benchmark results, comparing favorably to MaxLFQ for both data-dependent acquisition and data-independent acquisition. In addition, directLFQ provides normalized peptide intensity estimates for peptide-level comparisons. It is an important part of an overall quantitative proteomic pipeline that also needs to include high sensitive statistical analysis leading to proteoform resolution. Available as an open-source Python package and a graphical user interface with a one-click installer, it can be used in the AlphaPept ecosystem as well as downstream of most common computational proteomics pipelines.

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).1.

Re-mining serum proteomics data reveals extensive post-translational modifications upon Zika and dengue infection

Zika virus (ZIKV) and dengue virus (DENV) are two closely related flaviviruses with similar symptoms. However, due to the implications of ZIKV infections for pregnancy outcomes, understanding differences in their molecular impact on the host is of high interest. Viral infections change the host proteome, including post-translational modifications. As modifications are diverse and of low abundance, they typically require additional sample processing which is not feasible for large cohort studies. Therefore, we tested the potential of next-generation proteomics data in its ability to prioritize specific modifications for later analysis. We re-mined published mass spectra from 122 serum samples from ZIKV and DENV patients for the presence of phosphorylated, methylated, oxidized, glycosylated/glycated, sulfated, and carboxylated peptides. We identified 246 modified peptides with significantly differential abundance in ZIKV and DENV patients. Amongst these, methionine-oxidized peptides from apolipoproteins and glycosylated peptides from immunoglobulin proteins were more abundant in ZIKV patient serum and generate hypotheses on the potential roles of the modification in the infection. The results demonstrate how data-independent acquisition techniques can help prioritize future analyses of peptide modifications.

A panel of blood biomarkers unique to sudden cardiac arrest

BACKGROUND

The identification of circulating biomarkers specific for sudden cardiac arrest (SCA) could enhance risk prediction. Of particular interest are biomarkers specific to SCA, independent of coronary artery disease (CAD).

OBJECTIVE

The purpose of this study was to identify biomarkers of SCA obtained close to the SCA event.

METHODS

Twenty cases (survivors of SCA) and 40 age- and sex-matched controls were compared, with a replication analysis of 29 cases matched to 57 controls. A secondary analysis compared 20 SCA cases to 20 controls with CAD. Blood samples were obtained from SCA survivors at a median of 11 months after the SCA event. Proteins were analyzed on a mass spectrometer using data-independent acquisition; a subset of cytokines were analyzed using immunoassays; and 1153 lipids (13 classes) were analyzed. A false discovery rate P value of <.05 identified associated proteins.

RESULTS

Patients had a mean age of 58 years (range 25-87 years), and 70% were male. A total of 26 protein biomarkers associated with SCA when cases were compared with controls, of which 20 differentiated SCA from CAD. The replication analysis identified 8 of 26 biomarkers, of which 6 were not overlapping with CAD. The top identified biological processes involved the extracellular matrix, coagulation cascades, and platelet activation. Lipids in the lysophosphatidylcholine class were implicated in SCA through the CAD pathway.

CONCLUSION

We identified a panel of novel blood biomarkers specifically associated with SCA, including several that may be involved outside the CAD pathway. These biomarkers could have mechanistic significance and the potential to enhance clinical prediction of SCA.

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.

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.

NeuroLINCS Proteomics: Defining human-derived iPSC proteomes and protein signatures of pluripotency

The National Institute of Health (NIH) Library of integrated network-based cellular signatures (LINCS) program is premised on the generation of a publicly available data resource of cell-based biochemical responses or "signatures" to genetic or environmental perturbations. NeuroLINCS uses human inducible pluripotent stem cells (hiPSCs), derived from patients and healthy controls, and differentiated into motor neuron cell cultures. This multi-laboratory effort strives to establish i) robust multi-omic workflows for hiPSC and differentiated neuronal cultures, ii) public annotated data sets and iii) relevant and targetable biological pathways of spinal muscular atrophy (SMA) and amyotrophic lateral sclerosis (ALS). Here, we focus on the proteomics and the quality of the developed workflow of hiPSC lines from 6 individuals, though epigenomics and transcriptomics data are also publicly available. Known and commonly used markers representing 73 proteins were reproducibly quantified with consistent expression levels across all hiPSC lines. Data quality assessments, data levels and metadata of all 6 genetically diverse human iPSCs analysed by DIA-MS are parsable and available as a high-quality resource to the public.

Digoxin Induces Human Astrocyte Reaction In Vitro

Astrocyte reaction is a complex cellular process involving astrocytes in response to various types of CNS injury and a marker of neurotoxicity. It has been abundantly studied in rodents but relatively poorly in human cells due to limited access to the brain. Astrocytes play important roles in cerebral energy metabolism and are also key players in neuroinflammation. Astroglial metabolic and inflammatory changes have been reported with age, leading to the hypothesis that mitochondrial metabolism and inflammatory responses are interconnected. However, the relationship between energy metabolism and astrocyte reactivity in the context of neurotoxicity is not known. We hypothesized that changes in energy metabolism of astrocytes will be coupled to their activation by xenobiotics. Astrocyte reaction and associated energy metabolic changes were assessed by immunostaining, gene expression, proteomics, metabolomics, and extracellular flux analyses after 24 h of exposure of human ReN-derived astrocytes to digoxin (1-10 µM) or TNFα (30 ng/ml) used as a positive control. Strong astrocytic reaction was observed, accompanied by increased glycolysis at low concentrations of digoxin (0.1 and 0.5 µM) and after TNFα exposure, suggesting that increased glycolysis may be a common feature of reactive astrocytes, independent of the triggering molecule. In conclusion, whether astrocyte activation is triggered by cytokines or a xenobiotic, it is strongly tied to energy metabolism in human ReN-derived astrocytes. Increased glycolysis might be considered as an endpoint to detect astrocyte activation by potentially neurotoxic compounds in vitro. Finally, ReN-derived astrocytes may help to decipher mechanisms of neurotoxicity in ascertaining the ability of chemicals to directly target astrocytes.

Oxidative Stress and Extracellular Matrix Remodeling Are Signature Pathways of Extracellular Vesicles Released upon Morphine Exposure on Human Brain Microvascular Endothelial Cells

Morphine, a commonly used antinociceptive drug in hospitals, is known to cross the blood-brain barrier (BBB) by first passing through brain endothelial cells. Despite its pain-relieving effect, morphine also has detrimental effects, such as the potential induction of redox imbalance in the brain. However, there is still insufficient evidence of these effects on the brain, particularly on the brain endothelial cells and the extracellular vesicles that they naturally release. Indeed, extracellular vesicles (EVs) are nanosized bioparticles produced by almost all cell types and are currently thought to reflect the physiological state of their parent cells. These vesicles have emerged as a promising source of biomarkers by indicating the functional or dysfunctional state of their parent cells and, thus, allowing a better understanding of the biological processes involved in an adverse state. However, there is very little information on the morphine effect on human brain microvascular endothelial cells (HBMECs), and even less on their released EVs. Therefore, the current study aimed at unraveling the detrimental mechanisms of morphine exposure (at 1, 10, 25, 50 and 100 µM) for 24 h on human brain microvascular endothelial cells as well as on their associated EVs. Isolation of EVs was carried out using an affinity-based method. Several orthogonal techniques (NTA, western blotting and proteomics analysis) were used to validate the EVs enrichment, quality and concentration. Data-independent mass spectrometry (DIA-MS)-based proteomics was applied in order to analyze the proteome modulations induced by morphine on HBMECs and EVs. We were able to quantify almost 5500 proteins in HBMECs and 1500 proteins in EVs, of which 256 and 148, respectively, were found to be differentially expressed in at least one condition. Pathway enrichment analysis revealed that the "cell adhesion and extracellular matrix remodeling" process and the "HIF1 pathway", a pathway related to oxidative stress responses, were significantly modulated upon morphine exposure in HBMECs and EVs. Altogether, the combination of proteomics and bioinformatics findings highlighted shared pathways between HBMECs exposed to morphine and their released EVs. These results put forward molecular signatures of morphine-induced toxicity in HBMECs that were also carried by EVs. Therefore, EVs could potentially be regarded as a useful tool to investigate brain endothelial cells dysfunction, and to a different extent, the BBB dysfunction in patient circulation using these "signature pathways".

Proteomic Signatures of the Serological Response to Influenza Vaccination in a Large Human Cohort Study

The serological response to the influenza virus vaccine is highly heterogeneous for reasons that are not entirely clear. While the impact of demographic factors such as age, body mass index (BMI), sex, prior vaccination and titer levels are known to impact seroconversion, they only explain a fraction of the response. To identify signatures of the vaccine response, we analyzed 273 protein levels from 138 serum samples of influenza vaccine recipients (2019-2020 season). We found that levels of proteins functioning in cholesterol transport were positively associated with seroconversion, likely linking to the known impact of BMI. When adjusting seroconversion for the demographic factors, we identified additional, unexpected signatures: proteins regulating actin cytoskeleton dynamics were significantly elevated in participants with high adjusted seroconversion. Viral strain specific analysis showed that this trend was largely driven by the H3N2 strain. Further, we identified complex associations between adjusted seroconversion and other factors: levels of proteins of the complement system associated positively with adjusted seroconversion in younger participants, while they were associated negatively in the older population. We observed the opposite trends for proteins of high density lipoprotein remodeling, transcription, and hemostasis. In sum, careful integrative modeling can extract new signatures of seroconversion from highly variable data that suggest links between the humoral response as well as immune cell communication and migration.

Global, in situ analysis of the structural proteome in individuals with Parkinson's disease to identify a new class of biomarker

Parkinson's disease (PD) is a prevalent neurodegenerative disease for which robust biomarkers are needed. Because protein structure reflects function, we tested whether global, in situ analysis of protein structural changes provides insight into PD pathophysiology and could inform a new concept of structural disease biomarkers. Using limited proteolysis-mass spectrometry (LiP-MS), we identified 76 structurally altered proteins in cerebrospinal fluid (CSF) of individuals with PD relative to healthy donors. These proteins were enriched in processes misregulated in PD, and some proteins also showed structural changes in PD brain samples. CSF protein structural information outperformed abundance information in discriminating between healthy participants and those with PD and improved the discriminatory performance of CSF measures of the hallmark PD protein α-synuclein. We also present the first analysis of inter-individual variability of a structural proteome in healthy individuals, identifying biophysical features of variable protein regions. Although independent validation is needed, our data suggest that global analyses of the human structural proteome will guide the development of novel structural biomarkers of disease and enable hypothesis generation about underlying disease processes.

Data-independent acquisition proteomics methods for analyzing post-translational modifications

Protein post-translational modifications (PTMs) increase the functional diversity of the cellular proteome. Accurate and high throughput identification and quantification of protein PTMs is a key task in proteomics research. Recent advancements in data-independent acquisition (DIA) mass spectrometry (MS) technology have achieved deep coverage and accurate quantification of proteins and PTMs. This review provides an overview of DIA data processing methods that cover three aspects of PTMs analysis, i.e., detection of PTMs, site localization, and characterization of complex modification moieties, such as glycosylation. In addition, a survey of deep learning methods that boost DIA-based PTMs analysis is presented, including in silico spectral library generation, as well as feature scoring and error rate control. The limitations and future directions of DIA methods for PTMs analysis are also discussed. Novel data analysis methods will take advantage of advanced MS instrumentation techniques to empower DIA MS for in-depth and accurate PTMs measurements. This article is protected by copyright. All rights reserved.

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.

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.

Morphine-induced modulation of Nrf2-antioxidant response element signaling pathway in primary human brain microvascular endothelial cells

Morphine is one of the most potent opioid analgesic used for pain treatment. Morphine action in the central nervous system requires crossing the blood-brain barrier. Due to the controversial relationship between morphine and oxidative stress, the potential pro- or antioxidant effects of morphine in the blood-brain barrier is important to be understood, as oxidative stress could cause its disruption and predispose to neurodegenerative diseases. However, investigation is scarce in human brain endothelial cells. Therefore, the present study evaluated the impact of morphine exposure at three different concentrations (1, 10 and 100 µM) for 24 h and 48 h on primary human brain microvascular endothelial cells. A quantitative data-independent acquisition mass spectrometry strategy was used to analyze proteome modulations. Almost 3000 proteins were quantified of which 217 were reported to be significantly regulated in at least one condition versus untreated control. Pathway enrichment analysis unveiled dysregulation of the Nrf2 pathway involved in oxidative stress response. Seahorse assay underlined mitochondria dysfunctions, which were supported by significant expression modulations of relevant mitochondrial proteins. In conclusion, our study revealed the dysregulation of the Nrf2 pathway and mitochondria dysfunctions after morphine exposure, highlighting a potential redox imbalance in human brain endothelial cells.

Answer ALS, a large-scale resource for sporadic and familial ALS combining clinical and multi-omics data from induced pluripotent cell lines

Answer ALS is a biological and clinical resource of patient-derived, induced pluripotent stem (iPS) cell lines, multi-omic data derived from iPS neurons and longitudinal clinical and smartphone data from over 1,000 patients with ALS. This resource provides population-level biological and clinical data that may be employed to identify clinical-molecular-biochemical subtypes of amyotrophic lateral sclerosis (ALS). A unique smartphone-based system was employed to collect deep clinical data, including fine motor activity, speech, breathing and linguistics/cognition. The iPS spinal neurons were blood derived from each patient and these cells underwent multi-omic analytics including whole-genome sequencing, RNA transcriptomics, ATAC-sequencing and proteomics. The intent of these data is for the generation of integrated clinical and biological signatures using bioinformatics, statistics and computational biology to establish patterns that may lead to a better understanding of the underlying mechanisms of disease, including subgroup identification. A web portal for open-source sharing of all data was developed for widespread community-based data analytics.

S-adenosylmethionine inhibits la ribonucleoprotein domain family member 1 in murine liver and human liver cancer cells

BACKGROUND AND AIMS

Methionine adenosyltransferase 1A (MAT1A) is responsible for S-adenosylmethionine (SAMe) biosynthesis in the liver. Mice lacking Mat1a have hepatic SAMe depletion and develop NASH and HCC spontaneously. Several kinases are activated in Mat1a knockout (KO) mice livers. However, characterizing the phospho-proteome and determining whether they contribute to liver pathology remain open for study. Our study aimed to provide this knowledge.

APPROACH AND RESULTS

We performed phospho-proteomics in Mat1a KO mice livers with and without SAMe treatment to identify SAMe-dependent changes that may contribute to liver pathology. Our studies used Mat1a KO mice at different ages treated with and without SAMe, cell lines, in vitro translation and kinase assays, and human liver specimens. We found that the most striking change was hyperphosphorylation and increased content of La-related protein 1 (LARP1), which, in the unphosphorylated form, negatively regulates translation of 5'-terminal oligopyrimidine (TOP)-containing mRNAs. Consistently, multiple TOP proteins are induced in KO livers. Translation of TOP mRNAs ribosomal protein S3 and ribosomal protein L18 was enhanced by LARP1 overexpression in liver cancer cells. We identified LARP1-T449 as a SAMe-sensitive phospho-site of cyclin-dependent kinase 2 (CDK2). Knocking down CDK2 lowered LARP1 phosphorylation and prevented LARP1-overexpression-mediated increase in translation. LARP1-T449 phosphorylation induced global translation, cell growth, migration, invasion, and expression of oncogenic TOP-ribosomal proteins in HCC cells. LARP1 expression is increased in human NASH and HCC.

CONCLUSIONS

Our results reveal a SAMe-sensitive mechanism of LARP1 phosphorylation that may be involved in the progression of NASH to HCC.

Comparative assessment and novel strategy on methods for imputing proteomics data

Missing values are a major issue in quantitative proteomics analysis. While many methods have been developed for imputing missing values in high-throughput proteomics data, a comparative assessment of imputation accuracy remains inconclusive, mainly because mechanisms contributing to true missing values are complex and existing evaluation methodologies are imperfect. Moreover, few studies have provided an outlook of future methodological development. We first re-evaluate the performance of eight representative methods targeting three typical missing mechanisms. These methods are compared on both simulated and masked missing values embedded within real proteomics datasets, and performance is evaluated using three quantitative measures. We then introduce fused regularization matrix factorization, a low-rank global matrix factorization framework, capable of integrating local similarity derived from additional data types. We also explore a biologically-inspired latent variable modeling strategy—convex analysis of mixtures—for missing value imputation and present preliminary experimental results. While some winners emerged from our comparative assessment, the evaluation is intrinsically imperfect because performance is evaluated indirectly on artificial missing or masked values not authentic missing values. Nevertheless, we show that our fused regularization matrix factorization provides a novel incorporation of external and local information, and the exploratory implementation of convex analysis of mixtures presents a biologically plausible new approach.

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.

An integrated multi-omic analysis of iPSC-derived motor neurons from C9ORF72 ALS patients

Neurodegenerative diseases are challenging for systems biology because of the lack of reliable animal models or patient samples at early disease stages. Induced pluripotent stem cells (iPSCs) could address these challenges. We investigated DNA, RNA, epigenetics, and proteins in iPSC-derived motor neurons from patients with ALS carrying hexanucleotide expansions in C9ORF72. Using integrative computational methods combining all omics datasets, we identified novel and known dysregulated pathways. We used a C9ORF72 Drosophila model to distinguish pathways contributing to disease phenotypes from compensatory ones and confirmed alterations in some pathways in postmortem spinal cord tissue of patients with ALS. A different differentiation protocol was used to derive a separate set of C9ORF72 and control motor neurons. Many individual -omics differed by protocol, but some core dysregulated pathways were consistent. This strategy of analyzing patient-specific neurons provides disease-related outcomes with small numbers of heterogeneous lines and reduces variation from single-omics to elucidate network-based signatures.

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.

Paraquat-induced cholesterol biosynthesis proteins dysregulation in human brain microvascular endothelial cells

Despite Paraquat (PQ) being banned in several countries, it is still one of the most commonly used herbicides in agriculture. This compound is known to induce damaging effects on human and animal brain cells by generating Reactive Oxygen Species (ROS). However, there is few evidence of PQ effect on Human Brain Microvascular Endothelial Cells (HBMECs), one of the major component of the Blood–Brain Barrier (BBB). The present study aimed at unraveling biological mechanisms associated to the exposure of 1, 10 and 100 µM of PQ for 24 h on HBMECs. High-throughput mass spectrometry-based proteomics using data-independent acquisition (DIA) was applied. Biological pathway enrichment and cellular assays such as mitochondrial respiration and cholesterol level were performed to verify proteomics results. A total of 3753 proteins were quantified out of which 419 were significantly modulated by paraquat exposure. Biological pathway enrichment revealed the ubiquinone metabolism, a pathway directly linked to mitochondrial complex I proteins, confirming the well-known mechanism of PQ inducing oxidative stress. Additionally, this study also described the cholesterol biosynthesis modulation on HBMECs not yet described. In conclusion, our data indicate the toxic effect of PQ on HBMECs by downregulating proteins involved in mitochondrial complex I and cholesterol pathways.

Application of Mass Spectrometry-Based Proteomics to Barley Research

Barley (Hordeum vulgare) is the fourth most cultivated crop in the world in terms of production volume, and it is also the most important raw material of the malting and brewing industries. Barley belongs to the grass (Poaceae) family and plays an important role in food security and food safety for both humans and livestock. With the global population set to reach 9.7 billion by 2050, but with less available and/or suitable land for agriculture, the use of biotechnology tools in breeding programs are of considerable importance in the quest to meet the growing food gap. Proteomics as a member of the "omics" technologies has become popular for the investigation of proteins in cereal crops and particularly barley and its related products such as malt and beer. This technology has been applied to study how proteins in barley respond to adverse environmental conditions including abiotic and/or biotic stresses, how they are impacted during food processing including malting and brewing, and the presence of proteins implicated in celiac disease. Moreover, proteomics can be used in the future to inform breeding programs that aim to enhance the nutritional value and broaden the application of this crop in new food and beverage products. Mass spectrometry analysis is a valuable tool that, along with genomics and transcriptomics, can inform plant breeding strategies that aim to produce superior barley varieties. In this review, recent studies employing both qualitative and quantitative mass spectrometry approaches are explored with a focus on their application in cultivation, manufacturing, processing, quality, and the safety of barley and its related products.

obaDIA: one-step biological analysis pipeline for data-independent acquisition and other quantitative proteomics data

MOTIVATION

Data mining and data quality evaluation are indispensable constituents of quantitative proteomics, but few integrated tools available.

RESULTS

We introduced obaDIA, a one-step pipeline to generate visualizable and comprehensive results for quantitative proteomics data. obaDIA supports fragment-level, peptide-level and protein-level abundance matrices from DIA technique, as well as protein-level abundance matrices from other quantitative proteomic techniques. The result contains abundance matrix statistics, differential expression analysis, protein functional annotation and enrichment analysis. Additionally, enrichment strategies which use total proteins or expressed proteins as background are optional, and HTML based interactive visualization for differentially expressed proteins in the KEGG pathway is offered, which helps biological significance mining. In short, obaDIA is an automatic tool for bioinformatics analysis for quantitative proteomics.

AVAILABILITY AND IMPLEMENTATION

obaDIA is freely available from https://github.com/yjthu/obaDIA.git.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Comparison of the rabbit and human corneal endothelial proteomes regarding proliferative capacity

The shortage of human donor corneas has raised important concerns about engineering of corneal endothelial cells (CECs) for clinical use. However, due to the limited proliferative capacity of human CECs, driving them into proliferation and regeneration may be difficult. Unlike human CECs, rabbit CECs have a marked proliferative capacity. To clarify the potential reason for this difference, we analysed the proteomes of four human corneal endothelium samples and four rabbit corneal endothelium samples with quantitative label-free proteomics and downstream analysis. We discovered that vitamin and selenocompound metabolism and some signaling pathways such as NF-kappa B signaling pathway differed between the samples. Moreover, TGFβ, PITX2 and keratocan were distinctively expressed in rabbit samples, which might be associated with active proliferation in rabbit CECs. This study illustrates the proteomic differences between human and rabbit CECs and might promote CEC engineering strategies.

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.

Gene and protein expression in human megakaryocytes derived from induced pluripotent stem cells

BACKGROUND

There is interest in deriving megakaryocytes (MKs) from pluripotent stem cells (iPSC) for biological studies. We previously found that genomic structural integrity and genotype concordance is maintained in iPSC-derived MKs.

OBJECTIVE

To establish a comprehensive dataset of genes and proteins expressed in iPSC-derived MKs.

METHODS

iPSCs were reprogrammed from peripheral blood mononuclear cells (MNCs) and MKs were derived from the iPSCs in 194 healthy European American and African American subjects. mRNA was isolated and gene expression measured by RNA sequencing. Protein expression was measured in 62 of the subjects using mass spectrometry.

RESULTS AND CONCLUSIONS

MKs expressed genes and proteins known to be important in MK and platelet function and demonstrated good agreement with previous studies in human MKs derived from CD34+ progenitor cells. The percent of cells expressing the MK markers CD41 and CD42a was consistent in biological replicates, but variable across subjects, suggesting that unidentified subject-specific factors determine differentiation of MKs from iPSCs. Gene and protein sets important in platelet function were associated with increasing expression of CD41/42a, while those related to more basic cellular functions were associated with lower CD41/42a expression. There was differential gene expression by the sex and race (but not age) of the subject. Numerous genes and proteins were highly expressed in MKs but not known to play a role in MK or platelet function; these represent excellent candidates for future study of hematopoiesis, platelet formation, and/or platelet function.

MetaboKit: a comprehensive data extraction tool for untargeted metabolomics

We have developed MetaboKit, a comprehensive software package for compound identification and relative quantification in mass spectrometry-based untargeted metabolomics analysis. In data dependent acquisition (DDA) analysis, MetaboKit constructs a customized spectral library with compound identities from reference spectral libraries, adducts, dimers, in-source fragments (ISF), MS/MS fragmentation spectra, and more importantly the retention time information unique to the chromatography system used in the experiment. Using the customized library, the software performs targeted peak integration for precursor ions in DDA analysis and for precursor and product ions in data independent acquisition (DIA) analysis. With its stringent identification algorithm requiring matches by both MS and MS/MS data, MetaboKit provides identification results with significantly greater specificity than the competing software packages without loss in sensitivity. The proposed MS/MS-based screening of ISFs also reduces the chance of unverifiable identification of ISFs considerably. MetaboKit's quantification module produced peak area values highly correlated with known concentrations in a DIA analysis of the metabolite standards at both MS1 and MS2 levels. Moreover, the analysis of Cdk1Liv-/- mouse livers showed that MetaboKit can identify a wide range of lipid species and their ISFs, and quantitatively reconstitute the well-characterized fatty liver phenotype in these mice. In DIA data, the MS1-level and MS2-level peak area data produced similar fold change estimates in the differential abundance analysis, and the MS2-level peak area data allowed for quantitative comparisons in compounds whose precursor ion chromatogram was too noisy for peak integration.

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.

Plasma Metabolome and Lipidome Associations with Type 2 Diabetes and Diabetic Nephropathy

We conducted untargeted metabolomics analysis of plasma samples from a cross-sectional case–control study with 30 healthy controls, 30 patients with diabetes mellitus and normal renal function (DM-N), and 30 early diabetic nephropathy (DKD) patients using liquid chromatography-mass spectrometry (LC-MS). We employed two different modes of MS acquisition on a high-resolution MS instrument for identification and semi-quantification, and analyzed data using an advanced multivariate method for prioritizing differentially abundant metabolites. We obtained semi-quantification data for 1088 unique compounds (~55% lipids), excluding compounds that may be either exogenous compounds or treated as medication. Supervised classification analysis over a confounding-free partial correlation network shows that prostaglandins, phospholipids, nucleotides, sugars, and glycans are elevated in the DM-N and DKD patients, whereas glutamine, phenylacetylglutamine, 3-indoxyl sulfate, acetylphenylalanine, xanthine, dimethyluric acid, and asymmetric dimethylarginine are increased in DKD compared to DM-N. The data recapitulate the well-established plasma metabolome changes associated with DM-N and suggest uremic solutes and oxidative stress markers as the compounds indicating early renal function decline in DM patients.