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

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.

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.

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 (MS)-based proteomics enable the acquisition of increasingly large datasets within relatively short times, which exposes bottlenecks in the bioinformatics pipeline. Whereas 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 minutes and 100,000 proteomes in less than two hours - a thousand-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 (DDA) and data-independent acquisition (DIA). Additionally, directFQ 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 GUI 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.

•

Multi-omic analysis of differentiated C9ORF72 iPSC-derived motor neurons

•

Network-based integrative computational analysis

•

Pathogenic versus compensatory pathways elucidated using C9ORF72 Drosophila model

•

Pathways confirmed with alternative differentiation protocol and postmortem data

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.

Molecular pathways behind acquired obesity: Adipose tissue and skeletal muscle multiomics in monozygotic twin pairs discordant for BMI

Tissue-specific mechanisms prompting obesity-related development complications in humans remain unclear. We apply multiomics analyses of subcutaneous adipose tissue and skeletal muscle to examine the effects of acquired obesity among 49 BMI-discordant monozygotic twin pairs. Overall, adipose tissue appears to be more affected by excess body weight than skeletal muscle. In heavier co-twins, we observe a transcriptional pattern of downregulated mitochondrial pathways in both tissues and upregulated inflammatory pathways in adipose tissue. In adipose tissue, heavier co-twins exhibit lower creatine levels; in skeletal muscle, glycolysis- and redox stress-related protein and metabolite levels remain higher. Furthermore, metabolomics analyses in both tissues reveal that several proinflammatory lipids are higher and six of the same lipid derivatives are lower in acquired obesity. Finally, in adipose tissue, but not in skeletal muscle, mitochondrial downregulation and upregulated inflammation are associated with a fatty liver, insulin resistance, and dyslipidemia, suggesting that adipose tissue dominates in acquired obesity.

Multiomics analyses of adipose tissue and skeletal muscle in BMI-discordant twins

Excess body weight downregulates mitochondrial pathways in both tissues

Excess body weight upregulates proinflammatory pathways in both tissues

Adipose tissue alterations are associated with metabolic health in acquired obesity

New Proteomic Signatures to Distinguish Between Zika and Dengue Infections

Distinguishing between Zika and dengue virus infections is critical for accurate treatment, but we still lack detailed understanding of their impact on their host. To identify new protein signatures of the two infections, we used next-generation proteomics to profile 122 serum samples from 62 Zika and dengue patients. We quantified >500 proteins and identified 13 proteins that were significantly differentially expressed (adjusted p-value < 0.05). These proteins typically function in infection and wound healing, with several also linked to pregnancy and brain function. We successfully validated expression differences with Carbonic Anhydrase 2 in both the original and an independent sample set. Three of the differentially expressed proteins, i.e., Fibrinogen Alpha, Platelet Factor 4 Variant 1, and Pro-Platelet Basic Protein, predicted Zika virus infection at a ∼70% true-positive and 6% false-positive rate. Further, we showed that intraindividual temporal changes in protein signatures can disambiguate diagnoses and serve as indicators for past infections. Taken together, we demonstrate that serum proteomics can provide new resources that serve to distinguish between different viral infections.

•

Analysis of human serum samples with extreme protein abundance ranges

•

Unique protein signatures for Zika and dengue virus infection

•

Temporal changes in protein signatures as indicators for past infections

•

Machine learning to account for confounding factors

Isoform-resolved correlation analysis between mRNA abundance regulation and protein level degradation

Profiling of biological relationships between different molecular layers dissects regulatory mechanisms that ultimately determine cellular function. To thoroughly assess the role of protein post‐translational turnover, we devised a strategy combining pulse stable isotope‐labeled amino acids in cells (pSILAC), data‐independent acquisition mass spectrometry (DIA‐MS), and a novel data analysis framework that resolves protein degradation rate on the level of mRNA alternative splicing isoforms and isoform groups. We demonstrated our approach by the genome‐wide correlation analysis between mRNA amounts and protein degradation across different strains of HeLa cells that harbor a high grade of gene dosage variation. The dataset revealed that specific biological processes, cellular organelles, spatial compartments of organelles, and individual protein isoforms of the same genes could have distinctive degradation rate. The protein degradation diversity thus dissects the corresponding buffering or concerting protein turnover control across cancer cell lines. The data further indicate that specific mRNA splicing events such as intron retention significantly impact the protein abundance levels. Our findings support the tight association between transcriptome variability and proteostasis and provide a methodological foundation for studying functional protein degradation.

Nα-terminal acetylation of proteins by NatA and NatB serves distinct physiological roles in Saccharomyces cerevisiae

N-terminal (Nt) acetylation is a highly prevalent co-translational protein modification in eukaryotes, catalyzed by at least five Nt acetyltransferases (Nats) with differing specificities. Nt acetylation has been implicated in protein quality control, but its broad biological significance remains elusive. We investigate the roles of the two major Nats of S. cerevisiae, NatA and NatB, by performing transcriptome, translatome, and proteome profiling of natAΔ and natBΔ mutants. Our results reveal a range of NatA- and NatB-specific phenotypes. NatA is implicated in systemic adaptation control, because natAΔ mutants display altered expression of transposons, sub-telomeric genes, pheromone response genes, and nuclear genes encoding mitochondrial ribosomal proteins. NatB predominantly affects protein folding, because natBΔ mutants, to a greater extent than natA mutants, accumulate protein aggregates, induce stress responses, and display reduced fitness in the absence of the ribosome-associated chaperone Ssb. These phenotypic differences indicate that controlling Nat activities may serve to elicit distinct cellular responses.

Complement pathway changes at age 12 are associated with psychotic experiences at age 18 in a longitudinal population-based study: evidence for a role of stress

The complement cascade is a major component of the immune defence against infection, and there is increasing evidence for a role of dysregulated complement in major psychiatric disorders. We undertook a directed proteomic analysis of the complement signalling pathway (n = 29 proteins) using data-independent acquisition. Participants were recruited from the UK avon longitudinal study of parents and children (ALSPAC) cohort who participated in psychiatric assessment interviews at ages 12 and 18. Protein expression levels at age 12 among individuals who reported psychotic experiences (PEs) at age 18 (n = 64) were compared with age-matched controls (n = 67). Six out of the 29 targeted complement proteins or protein subcomponents were significantly upregulated following correction for multiple comparisons (VTN↑, C1RL↑, C8B↑, C8A↑, CFH↑, and C5↑). We then undertook an unbiased plasma proteomic analysis of mice exposed to chronic social stress and observed dysregulation of 11 complement proteins, including three that were altered in the same direction in individuals with PE (C1R↑, CFH↑, and C5↑). Our findings indicate that dysregulation of the complement protein pathway in blood is associated with incidence of psychotic experiences and that these changes may reflect exposure to stress.

ProteomeExpert: a docker image based web-server for exploring, modeling, visualizing, and mining quantitative proteomic data sets

Summary

The rapid progresses of high-throughput sequencing technology-based omics and mass spectrometry-based proteomics, such as data-independent acquisition and its penetration to clinical studies have generated increasing number of proteomic datasets containing hundreds to thousands of samples. To analyze these quantitative proteomic datasets and other omics (e.g. transcriptomics and metabolomics) datasets more efficiently and conveniently, we present a web server-based software tool ProteomeExpert implemented in Docker, which offers various analysis tools for experimental design, data mining, interpretation and visualization of quantitative proteomic datasets. ProteomeExpert can be deployed on an operating system with Docker installed or with R language environment.

Availability and implementation

The Docker image of ProteomeExpert is freely available from https://hub.docker.com/r/lifeinfo/proteomeexpert. The source code of ProteomeExpert is also openly accessible at http://www.github.com/guomics-lab/ProteomeExpert/. In addition, a demo server is provided at https://proteomic.shinyapps.io/peserver/.

Supplementary information

Supplementary data are available at Bioinformatics online.

MaxDIA enables library-based and library-free data-independent acquisition proteomics

MaxDIA is a software platform for analyzing data-independent acquisition (DIA) proteomics data within the MaxQuant software environment. Using spectral libraries, MaxDIA achieves deep proteome coverage with substantially better coefficients of variation in protein quantification than other software. MaxDIA is equipped with accurate false discovery rate (FDR) estimates on both library-to-DIA match and protein levels, including when using whole-proteome predicted spectral libraries. This is the foundation of discovery DIA-hypothesis-free analysis of DIA samples without library and with reliable FDR control. MaxDIA performs three- or four-dimensional feature detection of fragment data, and scoring of matches is augmented by machine learning on the features of an identification. MaxDIA's bootstrap DIA workflow performs multiple rounds of matching with increasing quality of recalibration and stringency of matching to the library. Combining MaxDIA with two new technologies-BoxCar acquisition and trapped ion mobility spectrometry-both lead to deep and accurate proteome quantification.