A systematic evaluation of normalization methods in quantitative label-free proteomics

CRABP1-complexes in exosome secretion

BACKGROUND

Cellular retinoic acid binding protein 1 (CRABP1) mediates rapid, non-canonical activity of retinoic acid (RA) by forming signalosomes via protein-protein interactions. Two signalosomes have been identified previously: CRABP1-MAPK and CRABP1-CaMKII. Crabp1 knockout (CKO) mice exhibited altered exosome profiles, but the mechanism of CRABP1 action was unclear. This study aimed to screen for and identify novel CRABP1 signalosomes that could modulate exosome secretion by using a combinatorial approach involving biochemical, bioinformatic and molecular studies.

METHODS

Immunoprecipitation coupled with mass spectrometry (IP-MS) identified candidate CRABP1-interacting proteins which were subsequently analyzed using GO Term Enrichment, Functional Annotation Clustering; and Pathway Analysis. Gene expression analysis of CKO samples revealed altered expression of genes related to exosome biogenesis and secretion. The effect of CRABP1 on exosome secretion was then experimentally validated using CKO mice and a Crabp1 knockdown P19 cell line.

RESULTS

IP-MS identified CRABP1-interacting targets. Bioinformatic analyses revealed significant association with actin cytoskeletal dynamics, kinases, and exosome secretion. The effect of CRABP1 on exosome secretion was experimentally validated by comparing circulating exosome numbers of CKO and wild type (WT) mice, and secreted exosomes from WT and siCRABP1-P19 cells. Pathway analysis identified kinase signaling and Arp2/3 complex as the major pathways where CRABP1-signalosomes modulate exosome secretion, which was validated in the P19 system.

CONCLUSION

The combinatorial approach allowed efficient screening for and identification of novel CRABP1-signalosomes. The results uncovered a novel function of CRABP1 in modulating exosome secretion, and suggested that CRABP1 could play roles in modulating intercellular communication and signal propagation.

MAPK14 /p38α Shapes the Molecular Landscape of Endometrial Cancer and promotes Tumorigenic Characteristics

The molecular underpinnings of H igh G rade E ndometrial C arcinoma (HGEC) metastatic growth and survival are poorly understood. Here we show that ascites-derived and primary tumor HGEC cell lines in 3D spheroid culture faithfully recapitulate key features of malignant peritoneal effusion and exhibit fundamentally distinct transcriptomic, proteomic and metabolomic landscapes when compared with conventional 2D monolayers. Using genetic screening platform we identify MAPK14 (which encodes the protein kinase p38α) as a specific requirement for HGEC in spheroid culture. MAPK14 /p38α has broad roles in programing the phosphoproteome, transcriptome and metabolome of HGEC spheroids, yet has negligible impact on monolayer cultures. MAPK14 promotes tumorigenicity in vivo and is specifically required to sustain a sub-population of spheroid cells that is enriched in cancer stemness markers. Therefore, spheroid growth of HGEC activates unique biological programs, including p38α signaling, that cannot be captured using 2D culture models and are highly relevant to malignant disease pathology.

Serum proteomics of mother-infant dyads carrying HLA-conferred type 1 diabetes risk

In-utero and dietary factors make important contributions toward health and development in early childhood. In this respect, serum proteomics of maturing infants can provide insights into studies of childhood diseases, which together with perinatal proteomes could reveal further biological perspectives. Accordingly, to determine differences between feeding groups and changes in infancy, serum proteomics analyses of mother-infant dyads with HLA-conferred susceptibility to type 1 diabetes (n = 22), weaned to either an extensively hydrolyzed or regular cow's milk formula, were made. The LC-MS/MS analyses included samples from the beginning of third trimester, the time of delivery, 3 months postpartum, cord blood, and samples from the infants at 3, 6, 9, and 12 months. Correlations between ranked protein intensities were detected within the dyads, together with perinatal and age-related changes. Comparison with intestinal permeability data revealed a number of significant correlations, which could merit further consideration in this context.

Untargeted metabolomics reveal signatures of a healthy lifestyle

This cross-sectional study investigated differences in the plasma metabolome in two groups of adults that were of similar age but varied markedly in body composition and dietary and physical activity patterns. Study participants included 52 adults in the lifestyle group (LIFE) (28 males, 24 females) and 52 in the control group (CON) (27 males, 25 females). The results using an extensive untargeted ultra high-performance liquid chromatography-high resolution mass spectrometry (UHPLC-HRMS) metabolomics analysis with 10,535 metabolite peaks identified 486 important metabolites (variable influence on projections scores of VIP ≥ 1) and 16 significantly enriched metabolic pathways that differentiated LIFE and CON groups. A novel metabolite signature of positive lifestyle habits emerged from this analysis highlighted by lower plasma levels of numerous bile acids, an amino acid profile characterized by higher histidine and lower glutamic acid, glutamine, β-alanine, phenylalanine, tyrosine, and proline, an elevated vitamin D status, higher levels of beneficial fatty acids and gut microbiome catabolism metabolites from plant substrates, and reduced levels of N-glycan degradation metabolites and environmental contaminants. This study established that the plasma metabolome is strongly associated with body composition and lifestyle habits. The robust lifestyle metabolite signature identified in this study is consistent with an improved life expectancy and a reduced risk for chronic disease.

High-throughput mass spectrometry maps the sepsis plasma proteome and differences in patient response

Sepsis, the dysregulated host response to infection causing life-threatening organ dysfunction, is a global health challenge requiring better understanding of pathophysiology and new therapeutic approaches. Here, we applied high-throughput tandem mass spectrometry to delineate the plasma proteome for sepsis and comparator groups (noninfected critical illness, postoperative inflammation, and healthy volunteers) involving 2612 samples (from 1611 patients) and 4553 liquid chromatography-mass spectrometry analyses acquired through a single batch of continuous measurements, with a throughput of 100 samples per day. We show how this scale of data can delineate proteins, pathways, and coexpression modules in sepsis and be integrated with paired leukocyte transcriptomic data (837 samples from n = 649 patients). We mapped the plasma proteomic landscape of the host response in sepsis, including changes over time, and identified features relating to etiology, clinical phenotypes (including organ failures), and severity. This work reveals subphenotypes informative for sepsis response state, disease processes, and outcome; identifies potential biomarkers; and advances opportunities for a precision medicine approach to sepsis.

Transkingdom Network Analysis (TkNA): a systems framework for inferring causal factors underlying host-microbiota and other multi-omic interactions

We present Transkingdom Network Analysis (TkNA), a unique causal-inference analytical framework that offers a holistic view of biological systems by integrating data from multiple cohorts and diverse omics types. TkNA helps to decipher key players and mechanisms governing host-microbiota (or any multi-omic data) interactions in specific conditions or diseases. TkNA reconstructs a network that represents a statistical model capturing the complex relationships between different omics in the biological system. It identifies robust and reproducible patterns of fold change direction and correlation sign across several cohorts to select differential features and their per-group correlations. The framework then uses causality-sensitive metrics, statistical thresholds and topological criteria to determine the final edges forming the transkingdom network. With the subsequent network's topological features, TkNA identifies nodes controlling a given subnetwork or governing communication between kingdoms and/or subnetworks. The computational time for the millions of correlations necessary for network reconstruction in TkNA typically takes only a few minutes, varying with the study design. Unlike most other multi-omics approaches that find only associations, TkNA focuses on establishing causality while accounting for the complex structure of multi-omic data. It achieves this without requiring huge sample sizes. Moreover, the TkNA protocol is user friendly, requiring minimal installation and basic familiarity with Unix. Researchers can access the TkNA software at https://github.com/CAnBioNet/TkNA/ .

Optimizing differential expression analysis for proteomics data via high-performing rules and ensemble inference

Identification of differentially expressed proteins in a proteomics workflow typically encompasses five key steps: raw data quantification, expression matrix construction, matrix normalization, missing value imputation (MVI), and differential expression analysis. The plethora of options in each step makes it challenging to identify optimal workflows that maximize the identification of differentially expressed proteins. To identify optimal workflows and their common properties, we conduct an extensive study involving 34,576 combinatoric experiments on 24 gold standard spike-in datasets. Applying frequent pattern mining techniques to top-ranked workflows, we uncover high-performing rules that demonstrate optimality has conserved properties. Via machine learning, we confirm optimal workflows are indeed predictable, with average cross-validation F1 scores and Matthew's correlation coefficients surpassing 0.84. We introduce an ensemble inference to integrate results from individual top-performing workflows for expanding differential proteome coverage and resolve inconsistencies. Ensemble inference provides gains in pAUC (up to 4.61%) and G-mean (up to 11.14%) and facilitates effective aggregation of information across varied quantification approaches such as topN, directLFQ, MaxLFQ intensities, and spectral counts. However, further development and evaluation are needed to establish acceptable frameworks for conducting ensemble inference on multiple proteomics workflows.

A framework for quality control in quantitative proteomics

A thorough evaluation of the quality, reproducibility, and variability of bottom-up proteomics data is necessary at every stage of a workflow from planning to analysis. We share real-world case studies applying adaptable quality control (QC) measures to assess sample preparation, system function, and quantitative analysis. System suitability samples are repeatedly measured longitudinally with targeted methods, and we share examples where they are used on three instrument platforms to identify severe system failures and track function over months to years. Internal QCs incorporated at protein and peptide-level allow our team to assess sample preparation issues and to differentiate system failures from sample-specific issues. External QC samples prepared alongside our experimental samples are used to verify the consistency and quantitative potential of our results during batch correction and normalization before assessing biological phenotypes. We combine these controls with rapid analysis using Skyline, longitudinal QC metrics using AutoQC, and server-based data deposition using PanoramaWeb. We propose that this integrated approach to QC be used as a starting point for groups to facilitate rapid quality control assessment to ensure that valuable instrument time is used to collect the best quality data possible.

FlexStat: combinatory differentially expressed protein extraction

Motivation

Mass spectrometry-based system proteomics allows identification of dysregulated protein hubs and associated disease-related features. Obtaining differentially expressed proteins (DEPs) is the most important step of downstream bioinformatics analysis. However, the extraction of statistically significant DEPs from datasets with multiple experimental conditions or disease types through currently available tools remains a laborious task. More often such an analysis requires considerable bioinformatics expertise, making it inaccessible to researchers with limited computational analytics experience.

Results

To uncover the differences among the many conditions within the data in a user-friendly manner, here we introduce FlexStat, a web-based interface that extracts DEPs through combinatory analysis. This tool accepts a protein expression matrix as input and systematically generates DEP results for every conceivable combination of various experimental conditions or disease types. FlexStat includes a suite of robust statistical tools for data preprocessing, in addition to DEP extraction, and publication-ready visualization, which are built on established R scientific libraries in an automated manner. This analytics suite was validated in diverse public proteomic datasets to showcase its high performance of rapid and simultaneous pairwise comparisons of comprehensive datasets.

Availability and implementation

FlexStat is implemented in R and is freely available at https://jglab.shinyapps.io/flexstatv1-pipeline-only/. The source code is accessible at https://github.com/kts-desilva/FlexStat/tree/main.

Matrix Metalloproteinase-7 in Urinary Extracellular Vesicles Identifies Rapid Disease Progression in Autosomal Dominant Polycystic Kidney Disease

SIGNIFICANCE STATEMENT

There is an unmet need for biomarkers of disease progression in autosomal dominant polycystic kidney disease (ADPKD). This study investigated urinary extracellular vesicles (uEVs) as a source of such biomarkers. Proteomic analysis of uEVs identified matrix metalloproteinase 7 (MMP-7) as a biomarker predictive of rapid disease progression. In validation studies, MMP-7 was predictive in uEVs but not in whole urine, possibly because uEVs are primarily secreted by tubular epithelial cells. Indeed, single-nucleus RNA sequencing showed that MMP-7 was especially increased in proximal tubule and thick ascending limb cells, which were further characterized by a profibrotic phenotype. Together, these data suggest that MMP-7 is a biologically plausible and promising uEV biomarker for rapid disease progression in ADPKD.

BACKGROUND

In ADPKD, there is an unmet need for early markers of rapid disease progression to facilitate counseling and selection for kidney-protective therapy. Our aim was to identify markers for rapid disease progression in uEVs.

METHODS

Six paired case-control groups ( n =10-59/group) of cases with rapid disease progression and controls with stable disease were formed from two independent ADPKD cohorts, with matching by age, sex, total kidney volume, and genetic variant. Candidate uEV biomarkers were identified by mass spectrometry and further analyzed using immunoblotting and an ELISA. Single-nucleus RNA sequencing of healthy and ADPKD tissue was used to identify the cellular origin of the uEV biomarker.

RESULTS

In the discovery proteomics experiments, the protein abundance of MMP-7 was significantly higher in uEVs of patients with rapid disease progression compared with stable disease. In the validation groups, a significant >2-fold increase in uEV-MMP-7 in patients with rapid disease progression was confirmed using immunoblotting. By contrast, no significant difference in MMP-7 was found in whole urine using ELISA. Compared with healthy kidney tissue, ADPKD tissue had significantly higher MMP-7 expression in proximal tubule and thick ascending limb cells with a profibrotic phenotype.

CONCLUSIONS

Among patients with ADPKD, rapid disease progressors have higher uEV-associated MMP-7. Our findings also suggest that MMP-7 is a biologically plausible biomarker for more rapid disease progression.

LFQRatio: A Normalization Method to Decipher Quantitative Proteome Changes in Microbial Coculture Systems

The value of synthetic microbial communities in biotechnology is gaining traction due to their ability to undertake more complex metabolic tasks than monocultures. However, a thorough understanding of strain interactions, productivity, and stability is often required to optimize growth and scale up cultivation. Quantitative proteomics can provide valuable insights into how microbial strains adapt to changing conditions in biomanufacturing. However, current workflows and methodologies are not suitable for simple artificial coculture systems where strain ratios are dynamic. Here, we established a workflow for coculture proteomics using an exemplar system containing two members, Azotobacter vinelandii and Synechococcus elongatus. Factors affecting the quantitative accuracy of coculture proteomics were investigated, including peptide physicochemical characteristics such as molecular weight, isoelectric point, hydrophobicity, and dynamic range as well as factors relating to protein identification such as varying proteome size and shared peptides between species. Different quantification methods based on spectral counts and intensity were evaluated at the protein and cell level. We propose a new normalization method, named "LFQRatio", to reflect the relative contributions of two distinct cell types emerging from cell ratio changes during cocultivation. LFQRatio can be applied to real coculture proteomics experiments, providing accurate insights into quantitative proteome changes in each strain.

Metaproteomics of Gut Fungi in Gnotobiotic Mice

Intestinal fungi are a fundamental component of the gut microbiome and play important roles in mammalian host biology. At the same time, the contribution of gut fungi to host health and disease remains understudied due to their low abundance. In that respect, gnotobiotic animals with defined microbial populations of reduced complexity represent a well-suited model system that highlights the effects of low abundant gut fungi on host physiology and other members of the microbial community. In this chapter, a label-free quantitative metaproteomic approach for the characterization of simplified microbial communities in gnotobiotic mice is presented. The model allows for exploring various research questions on the role of gut fungi in disease pathogenesis, microbial ecosystem maturation, or host-microbiome crosstalk.

Measurement and utilization of the proteomic reactivity by mass spectrometry

Chemical proteomics, which involves studying the covalent modifications of proteins by small molecules, has significantly contributed to our understanding of protein function and has become an essential tool in drug discovery. Mass spectrometry (MS) is the primary method for identifying and quantifying protein-small molecule adducts. In this review, we discuss various methods for measuring proteomic reactivity using MS and covalent proteomics probes that engage through reactivity-driven and proximity-driven mechanisms. We highlight the applications of these methods and probes in live-cell measurements, drug target identification and validation, and characterizing protein-small molecule interactions. We conclude the review with current developments and future opportunities in the field, providing our perspectives on analytical considerations for MS-based analysis of the proteomic reactivity landscape.

Label-Free Quantitation of Endogenous Peptides

Liquid chromatography-mass spectrometry (LC-MS)-based peptidomics methods allow for the detection and identification of many peptides in a complex biological mixture in an untargeted manner. Quantitative peptidomics approaches allow for comparisons of peptide abundance between different samples, allowing one to draw conclusions about peptide differences as a function of experimental treatment or physiology. While stable isotope labeling is a powerful approach for quantitative proteomics and peptidomics, advances in mass spectrometry instrumentation and analysis tools have allowed label-free methods to gain popularity in recent years. In a general label-free quantitative peptidomics experiment, peak intensity information for each peptide is compared across multiple LC-MS runs. Here, we outline a general approach for label-free quantitative peptidomics experiments, including steps for sample preparation, LC-MS data acquisition, data processing, and statistical analysis. Special attention is paid to address run-to-run variability, which can lead to several major problems in label-free experiments. Overall, our method provides researchers with a framework for the development of their own quantitative peptidomics workflows applicable to quantitation of peptides from a wide variety of different biological sources.

Evaluating Machine Learning Methods of Analyzing Multiclass Metabolomics

Multiclass metabolomic studies have become popular for revealing the differences in multiple stages of complex diseases, various lifestyles, or the effects of specific treatments. In multiclass metabolomics, there are multiple data manipulation steps for analyzing raw data, which consist of data filtering, the imputation of missing values, data normalization, marker identification, sample separation, classification, and so on. In each step, several to dozens of machine learning methods can be chosen for the given data set, with potentially hundreds or thousands of method combinations in the whole data processing chain. Therefore, a clear understanding of these machine learning methods is helpful for selecting an appropriate method combination for obtaining stable and reliable analytical results of specific data. However, there has rarely been an overall introduction or evaluation of these methods based on multiclass metabolomic data. Herein, detailed descriptions of these machine learning methods in multiple data manipulation steps are reviewed. Moreover, an assessment of these methods was performed using a benchmark data set for multiclass metabolomics. First, 12 imputation methods for imputing missing values were evaluated based on the PSS (Procrustes statistical shape analysis) and NRMSE (normalized root-mean-square error) values. Second, 17 normalization methods for processing multiclass metabolomic data were evaluated by applying the PMAD (pooled median absolute deviation) value. Third, different methods of identifying markers of multiclass metabolomics were evaluated based on the CWrel (relative weighted consistency) value. Fourth, nine classification methods for constructing multiclass models were assessed using the AUC (area under the curve) value. Performance evaluations of machine learning methods are highly recommended to select the most appropriate method combination before performing the final analysis of the given data. Overall, detailed descriptions and evaluation of various machine learning methods are expected to improve analyses of multiclass metabolomic data.

Comparative biodegradation analysis of three compostable polyesters by a marine microbial community

IMPORTANCE

Biodegradable plastics can be used in applications where the end product cannot be efficiently recycled due to high levels of contaminations, e.g., food or soil. Some of these plastics have a dedicated end of life, such as composting, but their degradation in the marine environment is poorly understood. In this study we showed that marine microbial communities can degrade a range of biodegradable polymers with different physical and chemical properties and use these as a sole carbon source for growth. We have also provided insights into the degradation mechanisms using a combined metagenomic and metaproteomic approach. In addition, we have identified three new enzymes that are capable of degrading both aliphatic polymers and aliphatic-aromatic copolymers, which can be used for biotechnological applications.

Workflow for Evaluating Normalization Tools for Omics Data Using Supervised and Unsupervised Machine Learning

To achieve high quality omics results, systematic variability in mass spectrometry (MS) data must be adequately addressed. Effective data normalization is essential for minimizing this variability. The abundance of approaches and the data-dependent nature of normalization have led some researchers to develop open-source academic software for choosing the best approach. While these tools are certainly beneficial to the community, none of them meet all of the needs of all users, particularly users who want to test new strategies that are not available in these products. Herein, we present a simple and straightforward workflow that facilitates the identification of optimal normalization strategies using straightforward evaluation metrics, employing both supervised and unsupervised machine learning. The workflow offers a "DIY" aspect, where the performance of any normalization strategy can be evaluated for any type of MS data. As a demonstration of its utility, we apply this workflow on two distinct datasets, an ESI-MS dataset of extracted lipids from latent fingerprints and a cancer spheroid dataset of metabolites ionized by MALDI-MSI, for which we identified the best-performing normalization strategies.

Proteomic advance of ischemic stroke: preclinical, clinical, and intervention

Ischemic stroke (IS) is the most common type of stroke and is characterized by high rates of mortality and long-term injury. The prediction and early diagnosis of IS are therefore crucial for optimal clinical intervention. Proteomics has provided important techniques for exploring protein markers associated with IS, but there has been no systematic evaluation and review of research that has used these techniques. Here, we review the differential proteins that have been found in cell- and animal- based studies and clinical trials of IS in the past 10 years; determine the key pathological proteins that have been identified in clinical trials; summarize the target proteins affected by interventions aimed at treating IS, with a focus on traditional Chinese medicine treatments. Overall, we clarify findings and problems that have been identified in recent proteomics research on IS and provide suggestions for improvements in this area. We also suggest areas that could be explored for determining the pathogenesis and developing interventions for IS.

Clostridioides difficile ferrosome organelles combat nutritional immunity

Iron is indispensable for almost all forms of life but toxic at elevated levels1-4. To survive within their hosts, bacterial pathogens have evolved iron uptake, storage and detoxification strategies to maintain iron homeostasis1,5,6. Recent studies showed that three Gram-negative environmental anaerobes produce iron-containing ferrosome granules7,8. However, it remains unclear whether ferrosomes are generated exclusively by Gram-negative bacteria. The Gram-positive bacterium Clostridioides difficile is the leading cause of nosocomial and antibiotic-associated infections in the USA9. Here we report that C. difficile undergoes an intracellular iron biomineralization process and stores iron in membrane-bound ferrosome organelles containing non-crystalline iron phosphate biominerals. We found that a membrane protein (FezA) and a P1B6-ATPase transporter (FezB), repressed by both iron and the ferric uptake regulator Fur, are required for ferrosome formation and play an important role in iron homeostasis during transition from iron deficiency to excess. Additionally, ferrosomes are often localized adjacent to cellular membranes as shown by cryo-electron tomography. Furthermore, using two mouse models of C. difficile infection, we demonstrated that the ferrosome system is activated in the inflamed gut to combat calprotectin-mediated iron sequestration and is important for bacterial colonization and survival during C. difficile infection.

Using ExpressAnalyst for Comprehensive Gene Expression Analysis in Model and Non-Model Organisms

ExpressAnalyst is a web-based platform that enables intuitive, end-to-end transcriptomics and proteomics data analysis. Users can start from FASTQ files, gene/protein abundance tables, or gene/protein lists. ExpressAnalyst will perform read quantification, gene expression table processing and normalization, differential expression analysis, or meta-analysis with complex study designs. The results are presented via various interactive visualizations such as volcano plots, heatmaps, networks, and ridgeline charts, with built-in functional enrichment analysis to allow flexible data exploration and understanding. ExpressAnalyst currently contains built-in support for 29 common organisms. For non-model organisms without good reference genomes, it can perform comprehensive transcriptome profiling directly from RNA-seq reads. These common tasks are covered in 11 Basic Protocols. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: RNA-seq count table uploading, processing, and normalization Basic Protocol 2: Differential expression analysis with linear models Basic Protocol 3: Functional analysis with volcano plot, enrichment network, and ridgeline visualization Basic Protocol 4: Hierarchical clustering analysis of transcriptomics data using interactive heatmaps Basic Protocol 5: Cross-species gene expression analysis based on ortholog mapping results Basic Protocol 6: Proteomics and microarray data processing and normalization Basic Protocol 7: Preparing multiple gene expression tables for meta-analysis Basic Protocol 8: Statistical and functional meta-analysis of gene expression data Basic Protocol 9: Functional analysis of transcriptomics signatures Basic Protocol 10: Dose-response and time-series data analysis Basic Protocol 11: RNA-seq reads processing and quantification with and without reference transcriptomes.

Targeted serum proteomics of longitudinal samples from newly diagnosed youth with type 1 diabetes distinguishes markers of disease and C-peptide trajectory

AIMS/HYPOTHESIS

There is a growing need for markers that could help indicate the decline in beta cell function and recognise the need and efficacy of intervention in type 1 diabetes. Measurements of suitably selected serum markers could potentially provide a non-invasive and easily applicable solution to this challenge. Accordingly, we evaluated a broad panel of proteins previously associated with type 1 diabetes in serum from newly diagnosed individuals during the first year from diagnosis. To uncover associations with beta cell function, comparisons were made between these targeted proteomics measurements and changes in fasting C-peptide levels. To further distinguish proteins linked with the disease status, comparisons were made with measurements of the protein targets in age- and sex-matched autoantibody-negative unaffected family members (UFMs).

METHODS

Selected reaction monitoring (SRM) mass spectrometry analyses of serum, targeting 85 type 1 diabetes-associated proteins, were made. Sera from individuals diagnosed under 18 years (n=86) were drawn within 6 weeks of diagnosis and at 3, 6 and 12 months afterwards (288 samples in total). The SRM data were compared with fasting C-peptide/glucose data, which was interpreted as a measure of beta cell function. The protein data were further compared with cross-sectional SRM measurements from UFMs (n=194).

RESULTS

Eleven proteins had statistically significant associations with fasting C-peptide/glucose. Of these, apolipoprotein L1 and glutathione peroxidase 3 (GPX3) displayed the strongest positive and inverse associations, respectively. Changes in GPX3 levels during the first year after diagnosis indicated future fasting C-peptide/glucose levels. In addition, differences in the levels of 13 proteins were observed between the individuals with type 1 diabetes and the matched UFMs. These included GPX3, transthyretin, prothrombin, apolipoprotein C1 and members of the IGF family.

CONCLUSIONS/INTERPRETATION

The association of several targeted proteins with fasting C-peptide/glucose levels in the first year after diagnosis suggests their connection with the underlying changes accompanying alterations in beta cell function in type 1 diabetes. Moreover, the direction of change in GPX3 during the first year was indicative of subsequent fasting C-peptide/glucose levels, and supports further investigation of this and other serum protein measurements in future studies of beta cell function in type 1 diabetes.

Human cytomegalovirus seropositivity is associated with reduced patient survival during sepsis

BACKGROUND

Sepsis is one of the leading causes of death. Treatment attempts targeting the immune response regularly fail in clinical trials. As HCMV latency can modulate the immune response and changes the immune cell composition, we hypothesized that HCMV serostatus affects mortality in sepsis patients.

METHODS

We determined the HCMV serostatus (i.e., latency) of 410 prospectively enrolled patients of the multicenter SepsisDataNet.NRW study. Patients were recruited according to the SEPSIS-3 criteria and clinical data were recorded in an observational approach. We quantified 13 cytokines at Days 1, 4, and 8 after enrollment. Proteomics data were analyzed from the plasma samples of 171 patients.

RESULTS

The 30-day mortality was higher in HCMV-seropositive patients than in seronegative sepsis patients (38% vs. 25%, respectively; p = 0.008; HR, 1.656; 95% CI 1.135-2.417). This effect was observed independent of age (p = 0.010; HR, 1.673; 95% CI 1.131-2.477). The predictive value on the outcome of the increased concentrations of IL-6 was present only in the seropositive cohort (30-day mortality, 63% vs. 24%; HR 3.250; 95% CI 2.075-5.090; p < 0.001) with no significant differences in serum concentrations of IL-6 between the two groups. Procalcitonin and IL-10 exhibited the same behavior and were predictive of the outcome only in HCMV-seropositive patients.

CONCLUSION

We suggest that the predictive value of inflammation-associated biomarkers should be re-evaluated with regard to the HCMV serostatus. Targeting HCMV latency might open a new approach to selecting suitable patients for individualized treatment in sepsis.

Proteomic signaling of dual specificity phosphatase 4 (DUSP4) in Alzheimer's disease

DUSP4 is a member of the DUSP (Dual-Specificity Phosphatase) subfamily that is selective to the mitogen-activated protein kinases (MAPK) and has been implicated in a range of biological processes and functions in Alzheimer's disease (AD). In this study, we utilized stereotactic delivery of adeno-associated virus (AAV)-DUSP4 to overexpress DUSP4 in the dorsal hippocampus of 5xFAD and wildtype (WT) mice, then used mass spectrometry (MS)-based proteomics along with label-free quantification to profile the proteome and phosphoproteome in the hippocampus. We identified patterns of protein expression and phosphorylation that are modulated in 5xFAD mice and examined the sex-specific impact of DUSP4 overexpression on the 5xFAD proteome/phosphoproteome. In 5xFAD mice, a substantial number of proteins were up- or down-regulated in both male and female mice in comparison to age and sex-matched WT mice, many of which are involved in AD-related biological processes, such as the activated immune response or suppression of synaptic activities. Upon DUSP4 overexpression, significantly regulated proteins were found in pathways that were suppressed, such as the immune response, in male 5xFAD mice. In contrast, such a shift was absent in female mice. For the phosphoproteome, we detected an array of phosphorylation sites that are regulated in 5xFAD compared to WT, and are modulated by DUSP4 overexpression in each sex. Interestingly, the changes in 5xFAD- and DUSP4-associated phosphorylation occurred in opposite directions. Strikingly, both the 5xFAD- and DUSP4-associated phosphorylation changes were found for the most part in neurons, and play key roles in neuronal processes and synaptic function. Site-centric pathway analysis revealed that both the 5xFAD- and DUSP4-associated phosphorylation sites were enriched for a number of kinase sets in female, but only a limited number of sets of kinases in male mice. Taken together, our results suggest that male and female 5xFAD mice respond to DUSP4 overexpression via shared and sex-specific molecular mechanisms, which might underly similar reductions in amyloid pathology in both sexes, while learning deficits were reduced in only females with DUSP4 overexpression. Finally, we validated our findings with the sex-specific AD-associated proteomes in human cohorts and further developed DUSP4-centric proteomic network models and signaling maps for each sex.

Proteomic signaling of dual specificity phosphatase 4 (DUSP4) in Alzheimer's disease

DUSP4 is a member of the DUSP (Dual-Specificity Phosphatase) subfamily that is selective to the mitogen-activated protein kinases (MAPK) and has been implicated in a range of biological processes and functions in Alzheimer's disease (AD). In this study, we utilized stereotactic delivery of adeno-associated virus (AAV)-DUSP4 to overexpress DUSP4 in the dorsal hippocampus of 5×FAD and wildtype (WT) mice, then used mass spectrometry (MS)-based proteomics along with label-free quantification to profile the proteome and phosphoproteome in the hippocampus. We identified patterns of protein expression and phosphorylation that are modulated in 5×FAD mice and examined the sex-specific impact of DUSP4 overexpression on the 5×FAD proteome/phosphoproteome. In 5×FAD mice, a substantial number of proteins were up- or down-regulated in both male and female mice in comparison to age and sex-matched WT mice, many of which are involved in AD-related biological processes, such as the activated immune response or suppression of synaptic activities. Upon DUSP4 overexpression, significantly regulated proteins were found in pathways that were suppressed, such as the immune response, in male 5×FAD mice. In contrast, such a shift was absent in female mice. For the phosphoproteome, we detected an array of phosphorylation sites that are regulated in 5×FAD compared to WT, and are modulated by DUSP4 overexpression in each sex. Interestingly, the changes in 5×FAD- and DUSP4-associated phosphorylation occurred in opposite directions. Strikingly, both the 5×FAD- and DUSP4-associated phosphorylation changes were found for the most part in neurons, and play key roles in neuronal processes and synaptic function. Site-centric pathway analysis revealed that both the 5×FAD- and DUSP4-associated phosphorylation sites were enriched for a number of kinase sets in female, but only a limited number of sets of kinases in male mice. Taken together, our results suggest that male and female 5×FAD mice respond to DUSP4 overexpression via shared and sex-specific molecular mechanisms, which might underly similar reductions in amyloid pathology in both sexes, while learning deficits were reduced in only females with DUSP4 overexpression. Finally, we validated our findings with the sex-specific AD-associated proteomes in human cohorts and further developed DUSP4-centric proteomic network models and signaling maps for each sex.

Phospho-Analyst: An Interactive, Easy-to-Use Web Platform To Analyze Quantitative Phosphoproteomics Data

Phosphoproteomics is nowadays the method of choice to comprehensively identify and quantify thousands of phosphorylated peptides and their associated proteins with the goal of interrogating changes in signal transduction pathways and other cellular processes. One of the most popular software suites to analyze phosphoproteomic data sets is MaxQuant, which converts mass spectrometric raw data into quantitative information on phosphopeptides and proteins. However, despite the increased utilization of phosphoproteomics in biomedical research, simple and user-friendly tools supporting downstream statistical analysis and interpretation of these highly complex outputs are still lacking. We have therefore developed Phospho-Analyst, which─similar to its sibling LFQ-Analyst─is an easy-to-use, interactive web application specifically designed to reproducibly perform differential expression analyses with "one click" and to visualize phosphoproteomic results in a meaningful and practical manner. Furthermore, if quantitative total proteomic information is available for the same samples, Phospho-Analyst automatically normalizes all phosphoproteomic results to underlying protein abundance levels, thereby ensuring that only genuine changes in phosphorylation events are considered. As such, Phospho-Analyst can not only be used by experienced proteomic veterans but also by researchers without any prior knowledge in (phospho)proteomics, statistics, or bioinformatics. Phospho-Analyst, including a detailed manual, is freely available at https://analyst-suites.org/apps/phospho-analyst/.

Proteomic investigation of neural stem cell to oligodendrocyte precursor cell differentiation reveals phosphorylation-dependent Dclk1 processing

Oligodendrocytes are generated via a two-step mechanism from pluripotent neural stem cells (NSCs): after differentiation of NSCs to oligodendrocyte precursor/NG2 cells (OPCs), they further develop into mature oligodendrocytes. The first step of this differentiation process is only incompletely understood. In this study, we utilized the neurosphere assay to investigate NSC to OPC differentiation in a time course-dependent manner by mass spectrometry-based (phospho-) proteomics. We identify doublecortin-like kinase 1 (Dclk1) as one of the most prominently regulated proteins in both datasets, and show that it undergoes a gradual transition between its short/long isoform during NSC to OPC differentiation. This is regulated by phosphorylation of its SP-rich region, resulting in inhibition of proteolytic Dclk1 long cleavage, and therefore Dclk1 short generation. Through interactome analyses of different Dclk1 isoforms by proximity biotinylation, we characterize their individual putative interaction partners and substrates. All data are available via ProteomeXchange with identifier PXD040652.

Multi-omics profiling of CSF from spinal muscular atrophy type 3 patients after nusinersen treatment: a 2-year follow-up multicenter retrospective study

Spinal muscular atrophy (SMA) is a neurodegenerative disorder caused by mutations in the SMN1 gene resulting in reduced levels of the SMN protein. Nusinersen, the first antisense oligonucleotide (ASO) approved for SMA treatment, binds to the SMN2 gene, paralogue to SMN1, and mediates the translation of a functional SMN protein. Here, we used longitudinal high-resolution mass spectrometry (MS) to assess both global proteome and metabolome in cerebrospinal fluid (CSF) from ten SMA type 3 patients, with the aim of identifying novel readouts of pharmacodynamic/response to treatment and predictive markers of treatment response. Patients had a median age of 33.5 [29.5; 38.25] years, and 80% of them were ambulant at time of the enrolment, with a median HFMSE score of 37.5 [25.75; 50.75]. Untargeted CSF proteome and metabolome were measured using high-resolution MS (nLC-HRMS) on CSF samples obtained before treatment (T0) and after 2 years of follow-up (T22). A total of 26 proteins were found to be differentially expressed between T0 and T22 upon VSN normalization and LIMMA differential analysis, accounting for paired replica. Notably, key markers of the insulin-growth factor signaling pathway were upregulated after treatment together with selective modulation of key transcription regulators. Using CombiROC multimarker signature analysis, we suggest that detecting a reduction of SEMA6A and an increase of COL1A2 and GRIA4 might reflect therapeutic efficacy of nusinersen. Longitudinal metabolome profiling, analyzed with paired t-Test, showed a significant shift for some aminoacid utilization induced by treatment, whereas other metabolites were largely unchanged. Together, these data suggest perturbation upon nusinersen treatment still sustained after 22 months of follow-up and confirm the utility of CSF multi-omic profiling as pharmacodynamic biomarker for SMA type 3. Nonetheless, validation studies are needed to confirm this evidence in a larger sample size and to further dissect combined markers of response to treatment.

Unraveling the Global Proteome and Phosphoproteome of Prostate Cancer Patient-Derived Xenografts

Resistance to androgen deprivation therapies leads to metastatic castration-resistant prostate cancer (mCRPC) of adenocarcinoma (AdCa) origin that can transform to emergent aggressive variant prostate cancer (AVPC) which has neuroendocrine (NE)-like features. To this end, we used LuCaP patient-derived xenograft (PDX) tumors, clinically relevant models that reflects and retains key features of the tumor from advanced prostate cancer patients. Here we performed proteome and phosphoproteome characterization of 48 LuCaP PDX tumors and identified over 94,000 peptides and 9,700 phosphopeptides corresponding to 7,738 proteins. When we compared 15 NE versus 33 AdCa PDX samples, we identified 309 unique proteins and 476 unique phosphopeptides that were significantly altered and corresponded to proteins that are known to distinguish these two phenotypes. Assessment of protein and RNA concordance from these tumors revealed increased dissonance in transcriptionally regulated proteins in NE and metabolite interconversion enzymes in AdCa.

Microbiome-derived cobalamin and succinyl-CoA as biomarkers for improved screening of anal cancer

Human papillomavirus can cause preinvasive, high-grade squamous intraepithelial lesions (HSILs) as precursors to cancer in the anogenital area, and the microbiome is suggested to be a contributing factor. Men who have sex with men (MSM) living with human immunodeficiency virus (HIV) have a high risk of anal cancer, but current screening strategies for HSIL detection lack specificity. Here, we investigated the anal microbiome to improve HSIL screening. We enrolled participants living with HIV, divided into a discovery (n = 167) and validation cohort (n = 46), and who were predominantly (93.9%) cisgender MSM undergoing HSIL screening with high-resolution anoscopy and anal biopsies. We identified no microbiome composition signatures associated with HSILs, but elevated levels of microbiome-encoded proteins producing succinyl coenzyme A and cobalamin were significantly associated with HSILs in both cohorts. Measurement of these candidate biomarkers alone in anal cytobrushes outperformed anal cytology as a diagnostic indicator for HSILs, increasing the sensitivity from 91.2% to 96.6%, the specificity from 34.1% to 81.8%, and reclassifying 82% of false-positive results as true negatives. We propose that these two microbiome-derived biomarkers may improve the current strategy of anal cancer screening.

Tidyproteomics: an open-source R package and data object for quantitative proteomics post analysis and visualization

BACKGROUND

The analysis of mass spectrometry-based quantitative proteomics data can be challenging given the variety of established analysis platforms, the differences in reporting formats, and a general lack of approachable standardized post-processing analyses such as sample group statistics, quantitative variation and even data filtering. We developed tidyproteomics to facilitate basic analysis, improve data interoperability and potentially ease the integration of new processing algorithms, mainly through the use of a simplified data-object.

RESULTS

The R package tidyproteomics was developed as both a framework for standardizing quantitative proteomics data and a platform for analysis workflows, containing discrete functions that can be connected end-to-end, thus making it easier to define complex analyses by breaking them into small stepwise units. Additionally, as with any analysis workflow, choices made during analysis can have large impacts on the results and as such, tidyproteomics allows researchers to string each function together in any order, select from a variety of options and in some cases develop and incorporate custom algorithms.

CONCLUSIONS

Tidyproteomics aims to simplify data exploration from multiple platforms, provide control over individual functions and analysis order, and serve as a tool to assemble complex repeatable processing workflows in a logical flow. Datasets in tidyproteomics are easy to work with, have a structure that allows for biological annotations to be added, and come with a framework for developing additional analysis tools. The consistent data structure and accessible analysis and plotting tools also offers a way for researchers to save time on mundane data manipulation tasks.

The Edge Effect in High-Throughput Proteomics: A Cautionary Tale

In order for mass spectrometry to continue to grow as a platform for high-throughput clinical and translational research, careful consideration must be given to quality control by ensuring that the assay performs reproducibly and accurately and precisely. In particular, the throughput required for large cohort clinical validation in biomarker discovery and diagnostic screening has driven the growth of multiplexed targeted liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) assays paired with sample preparation and analysis in multiwell plates. However, large scale MS-based proteomics studies are often plagued by batch effects: sources of technical variation in the data, which can arise from a diverse array of sources such as sample preparation batches, different reagent lots, or indeed MS signal drift. These batch effects can confound the detection of true signal differences, resulting in incorrect conclusions being drawn about significant biological effects or lack thereof. Here, we present an intraplate batch effect termed the edge effect arising from temperature gradients in multiwell plates, commonly reported in preclinical cell culture studies but not yet reported in a clinical proteomics setting. We present methods herein to ameliorate the phenomenon including proper assessment of heating techniques for multiwell plates and incorporation of surrogate standards, which can normalize for intraplate variation.

Multi-omics analyses reveal ClpP activators disrupt essential mitochondrial pathways in triple-negative breast cancer

ClpP activators ONC201 and related small molecules (TR compounds, Madera Therapeutics), have demonstrated significant anti-cancer potential in vitro and in vivo studies, including clinical trials for refractory solid tumors. Though progress has been made in identifying specific phenotypic outcomes following ClpP activation, the exact mechanism by which ClpP activation leads to broad anti-cancer activity has yet to be fully elucidated. In this study, we utilized a multi-omics approach to identify the ClpP-dependent proteomic, transcriptomic, and metabolomic changes resulting from ONC201 or the TR compound TR-57 in triple-negative breast cancer cells. Applying mass spectrometry-based methods of proteomics and metabolomics, we identified ∼8,000 proteins and 588 metabolites, respectively. From proteomics data, 113 (ONC201) and 191 (TR-57) proteins significantly increased and 572 (ONC201) and 686 (TR-57) proteins significantly decreased in this study. Gene ontological (GO) analysis revealed strong similarities between proteins up- or downregulated by ONC201 or TR-57 treatment. Notably, this included the downregulation of many mitochondrial processes and proteins, including mitochondrial translation and mitochondrial matrix proteins. We performed a large-scale transcriptomic analysis of WT SUM159 cells, identifying ∼7,700 transcripts (746 and 1,100 significantly increasing, 795 and 1,013 significantly decreasing in ONC201 and TR-57 treated cells, respectively). Less than 21% of these genes were affected by these compounds in ClpP null cells. GO analysis of these data demonstrated additional similarity of response to ONC201 and TR-57, including a decrease in transcripts related to the mitochondrial inner membrane and matrix, cell cycle, and nucleus, and increases in other nuclear transcripts and transcripts related to metal-ion binding. Comparison of response between both compounds demonstrated a highly similar response in all -omics datasets. Analysis of metabolites also revealed significant similarities between ONC201 and TR-57 with increases in α-ketoglutarate and 2-hydroxyglutaric acid and decreased ureidosuccinic acid, L-ascorbic acid, L-serine, and cytidine observed following ClpP activation in TNBC cells. Further analysis identified multiple pathways that were specifically impacted by ClpP activation, including ATF4 activation, heme biosynthesis, and the citrulline/urea cycle. In summary the results of our studies demonstrate that ONC201 and TR-57 induce highly similar and broad effects against multiple mitochondrial processes required for cell proliferation.

Transkingdom Network Analysis (TkNA): a systems framework for inferring causal factors underlying host-microbiota and other multi-omic interactions

Technological advances have generated tremendous amounts of high-throughput omics data. Integrating data from multiple cohorts and diverse omics types from new and previously published studies can offer a holistic view of a biological system and aid in deciphering its critical players and key mechanisms. In this protocol, we describe how to use Transkingdom Network Analysis (TkNA), a unique causal-inference analytical framework that can perform meta-analysis of cohorts and detect master regulators among measured parameters that govern pathological or physiological responses of host-microbiota (or any multi-omic data) interactions in a particular condition or disease. TkNA first reconstructs the network that represents a statistical model capturing the complex relationships between the different omics of the biological system. Here, it selects differential features and their per-group correlations by identifying robust and reproducible patterns of fold change direction and sign of correlation across several cohorts. Next, a causality-sensitive metric, statistical thresholds, and a set of topological criteria are used to select the final edges that form the transkingdom network. The second part of the analysis involves interrogating the network. Using the network's local and global topology metrics, it detects nodes that are responsible for control of given subnetwork or control of communication between kingdoms and/or subnetworks. The underlying basis of the TkNA approach involves fundamental principles including laws of causality, graph theory and information theory. Hence, TkNA can be used for causal inference via network analysis of any host and/or microbiota multi-omics data. This quick and easy-to-run protocol requires very basic familiarity with the Unix command-line environment.

prolfqua: A Comprehensive R-Package for Proteomics Differential Expression Analysis

Mass spectrometry is widely used for quantitative proteomics studies, relative protein quantification, and differential expression analysis of proteins. There is a large variety of quantification software and analysis tools. Nevertheless, there is a need for a modular, easy-to-use application programming interface in R that transparently supports a variety of well principled statistical procedures to make applying them to proteomics data, comparing and understanding their differences easy. The prolfqua package integrates essential steps of the mass spectrometry-based differential expression analysis workflow: quality control, data normalization, protein aggregation, statistical modeling, hypothesis testing, and sample size estimation. The package makes integrating new data formats easy. It can be used to model simple experimental designs with a single explanatory variable and complex experiments with multiple factors and hypothesis testing. The implemented methods allow sensitive and specific differential expression analysis. Furthermore, the package implements benchmark functionality that can help to compare data acquisition, data preprocessing, or data modeling methods using a gold standard data set. The application programmer interface of prolfqua strives to be clear, predictable, discoverable, and consistent to make proteomics data analysis application development easy and exciting. Finally, the prolfqua R-package is available on GitHub https://github.com/fgcz/prolfqua, distributed under the MIT license. It runs on all platforms supported by the R free software environment for statistical computing and graphics.

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.

Intratympanic steroid treatment can reduce ROS and immune response in human perilymph investigated by in-depth proteome analysis

Intratympanic (IT) steroid treatment is one of the most widely used and effective treatments for inner ear disorders such as sudden sensorineural hearing loss (SNHL). However, a clear mechanism of IT steroids in inner ear recovery has not yet been revealed. Therefore, we investigated proteome changes in extracted human perilymph after steroid treatment. In this study, we applied a tandem mass spectrometry (MS/MS)-based proteomics approach to discover global proteome changes by comparing human perilymph after steroid treatment with non-treated perilymph group. Using liquid chromatography-MS/MS analysis, we selected 156 differentially expressed proteins (DEPs) that were statistically significant according to Student's t-test. Functional annotation analysis showed that upregulated proteins after steroid treatment are related to apoptosis signaling, as well as reactive oxygen species (ROS) and immune responses. The protein-protein interaction (PPI) clusters the proteins associated with these processes and attempts to observe signaling circuitry, which mediates cellular response after IT steroid treatments. Moreover, we also considered the interactome analysis of DEPs and observed that those with high interaction scores were categorized as having equivalent molecular functions (MFs). Collectively, we suggest that DEPs and interacting proteins in human perilymph after steroid treatment would inhibit the apoptotic and adaptive immune processes that may lead to anti-inflammatory effects.

Multivariate Analysis with the R Package mixOmics

The high-dimensional nature of proteomics data presents challenges for statistical analysis and biological interpretation. Multivariate analysis, combined with insightful visualization can help to reveal the underlying patterns in complex biological data. This chapter introduces the R package mixOmics which focuses on data exploration and integration. We first introduce methods for single data sets: both Principal Component Analysis, which can identify the patterns of variance present in data, and sparse Partial Least Squares Discriminant Analysis, which aims to identify variables that can classify samples into known groups. We then present integrative methods with Projection to Latent Structures and further extensions for discriminant analysis. We illustrate each technique on a breast cancer multi-omics study and provide the R code and data as online supplementary material for readers interested in reproducing these analyses.

MS-DAP Platform for Downstream Data Analysis of Label-Free Proteomics Uncovers Optimal Workflows in Benchmark Data Sets and Increased Sensitivity in Analysis of Alzheimer's Biomarker Data

In the rapidly moving proteomics field, a diverse patchwork of data analysis pipelines and algorithms for data normalization and differential expression analysis is used by the community. We generated a mass spectrometry downstream analysis pipeline (MS-DAP) that integrates both popular and recently developed algorithms for normalization and statistical analyses. Additional algorithms can be easily added in the future as plugins. MS-DAP is open-source and facilitates transparent and reproducible proteome science by generating extensive data visualizations and quality reporting, provided as standardized PDF reports. Second, we performed a systematic evaluation of methods for normalization and statistical analysis on a large variety of data sets, including additional data generated in this study, which revealed key differences. Commonly used approaches for differential testing based on moderated t-statistics were consistently outperformed by more recent statistical models, all integrated in MS-DAP. Third, we introduced a novel normalization algorithm that rescues deficiencies observed in commonly used normalization methods. Finally, we used the MS-DAP platform to reanalyze a recently published large-scale proteomics data set of CSF from AD patients. This revealed increased sensitivity, resulting in additional significant target proteins which improved overlap with results reported in related studies and includes a large set of new potential AD biomarkers in addition to previously reported.

Translation factor eIF5a is essential for IFNγ production and cell cycle regulation in primary CD8+ T lymphocytes

Control of mRNA translation adjusts protein production rapidly and facilitates local cellular responses to environmental conditions. Traditionally initiation of translation is considered to be a major translational control point, however, control of peptide elongation is also important. Here we show that the function of the elongation factor, eIF5a, is regulated dynamically in naïve CD8+ T cells upon activation by post-translational modification, whereupon it facilitates translation of specific subsets of proteins. eIF5a is essential for long-term survival of effector CD8+ T cells and sequencing of nascent polypeptides indicates that the production of proteins which regulate proliferation and key effector functions, particularly the production of IFNγ and less acutely TNF production and cytotoxicity, is dependent on the presence of functional eIF5a. Control of translation in multiple immune cell lineages is required to co-ordinate immune responses and these data illustrate that translational elongation contributes to post-transcriptional regulons important for the control of inflammation.

Multi-omic interactions in the gut of children at the onset of islet autoimmunity

BACKGROUND

The gastrointestinal ecosystem is a highly complex environment with a profound influence on human health. Inflammation in the gut, linked to an altered gut microbiome, has been associated with the development of multiple human conditions including type 1 diabetes (T1D). Viruses infecting the gastrointestinal tract, especially enteroviruses, are also thought to play an important role in T1D pathogenesis possibly via overlapping mechanisms. However, it is not known whether the microbiome and virome act together or which risk factor may be of greater importance at the time when islet autoimmunity is initiated.

RESULTS

Here, we apply an integrative approach to combine comprehensive fecal virome, microbiome, and metaproteome data sampled before and at the onset of islet autoimmunity in 40 children at increased risk of T1D. We show strong age-related effects, with microbial and metaproteome diversity increasing with age while host antibody number and abundance declined with age. Mastadenovirus, which has been associated with a reduced risk of T1D, was associated with profound changes in the metaproteome indicating a functional shift in the microbiota. Multi-omic factor analysis modeling revealed a cluster of proteins associated with carbohydrate transport from the genus Faecalibacterium were associated with islet autoimmunity.

CONCLUSIONS

These findings demonstrate the interrelatedness of the gut microbiota, metaproteome and virome in young children. We show a functional remodeling of the gut microbiota accompanies both islet autoimmunity and viral infection with a switch in function in Faecalibacterium occurring at the onset of islet autoimmunity. Video Abstract.

Exploratory Metabolomics Underscores the Folate Enzyme ALDH1L1 as a Regulator of Glycine and Methylation Reactions

Folate (vitamin B9) is involved in one-carbon transfer reactions and plays a significant role in nucleic acid synthesis and control of cellular proliferation, among other key cellular processes. It is now recognized that the role of folates in different stages of carcinogenesis is complex, and more research is needed to understand how folate reactions become dysregulated in cancers and the metabolic consequences that occur as a result. ALDH1L1 (cytosolic 10-formyltetrahydrofolate dehydrogenase), an enzyme of folate metabolism expressed in many tissues, is ubiquitously downregulated in cancers and is not expressed in cancer cell lines. The RT4 cell line (derived from papillary bladder cancer) which expresses high levels of ALDH1L1 represents an exception, providing an opportunity to explore the metabolic consequences of the loss of this enzyme. We have downregulated this protein in RT4 cells (shRNA driven knockdown or CRISPR driven knockout) and compared metabolomes of ALDH1L1-expressing and -deficient cells to determine if metabolic changes linked to the loss of this enzyme might provide proliferative and/or survival advantages for cancer cells. In this study, cell extracts were analyzed using Ultra High Performance Liquid Chromatography High Resolution Mass Spectrometry (UHPLC-HR-MS). A total of 13,339 signals were identified or annotated using an in-house library and public databases. Supervised and unsupervised multivariate analysis revealed metabolic differences between RT4 cells and ALDH1L1-deficient clones. Glycine (8-fold decrease) and metabolites derived from S-adenosylmethionine utilizing pathways were significantly decreased in the ALDH1L1-deficient clones, compared with RT4 cells. Other changes linked to ALDH1L1 downregulation include decreased levels of amino acids, Krebs cycle intermediates, and ribose-5-phosphate, and increased nicotinic acid. While the ALDH1L1-catalyzed reaction is directly linked to glycine biosynthesis and methyl group flux, its overall effect on cellular metabolism extends beyond immediate metabolic pathways controlled by this enzyme.

Does filter-aided sample preparation provide sufficient method linearity for quantitative plant shotgun proteomics?

Due to its outstanding throughput and analytical resolution, gel-free LC-based shotgun proteomics represents the gold standard of proteome analysis. Thereby, the efficiency of sample preparation dramatically affects the correctness and reliability of protein quantification. Thus, the steps of protein isolation, solubilization, and proteolysis represent the principal bottleneck of shotgun proteomics. The desired performance of the sample preparation protocols can be achieved by the application of detergents. However, these compounds ultimately compromise reverse-phase chromatographic separation and disrupt electrospray ionization. Filter-aided sample preparation (FASP) represents an elegant approach to overcome these limitations. Although this method is comprehensively validated for cell proteomics, its applicability to plants and compatibility with plant-specific protein isolation protocols remain to be confirmed. Thereby, the most important gap is the absence of the data on the linearity of underlying protein quantification methods for plant matrices. To fill this gap, we address here the potential of FASP in combination with two protein isolation protocols for quantitative analysis of pea (Pisum sativum) seed and Arabidopsis thaliana leaf proteomes by the shotgun approach. For this aim, in comprehensive spiking experiments with bovine serum albumin (BSA), we evaluated the linear dynamic range (LDR) of protein quantification in the presence of plant matrices. Furthermore, we addressed the interference of two different plant matrices in quantitative experiments, accomplished with two alternative sample preparation workflows in comparison to conventional FASP-based digestion of cell lysates, considered here as a reference. The spiking experiments revealed high sensitivities (LODs of up to 4 fmol) for spiked BSA and LDRs of at least 0.6 × 10². Thereby, phenol extraction yielded slightly better recoveries, whereas the detergent-based method showed better linearity. Thus, our results indicate the very good applicability of FASP to quantitative plant proteomics with only limited impact of the protein isolation technique on the method's overall performance.

Microbiota independent effects of oligosaccharides on Caco-2 cells -A semi-targeted metabolomics approach using DI-FT-ICR-MS coupled with pathway enrichment analysis

Milk oligosaccharides (MOS) and galactooligosaccharides (GOS) are associated with many benefits, including anti-microbial effects and immune-modulating properties. However, the cellular mechanisms of these are largely unknown. In this study, the effects of enriched GOS and MOS mixtures from caprine and bovine milk consisting mainly 6'-galactosyllactose, 3'-sialyllactose, and 6'-sialyllactose on Caco-2 cells were investigated, and the treatment-specific metabolomes were described. In the control, the cells were treated with a sugar mix consisting of one-third each of glucose, galactose and lactose. A local metabolomics workflow with pathway enrichment was established, which specifically addresses DI-FT-ICR-MS analyses and includes adaptations in terms of measurement technology and sample matrices. By including quality parameters, especially the isotope pattern, we increased the precision of annotation. The independence from online tools, the fast adaptability to changes in databases, and the specific adjustment to the measurement technology and biomaterial used, proved to be a great advantage. For the first time it was possible to find 71 active pathways in a Caco-2 cell experiment. These pathways were assigned to 12 main categories, with amino acid metabolism and carbohydrate metabolism being the most dominant categories in terms of the number of metabolites and metabolic pathways. Treatment of Caco-2 cells with high GOS and glucose contents resulted in significant effects on several metabolic pathways, whereas the MOS containing treatments resulted only for individual metabolites in significant changes. An effect based on bovine or caprine origin alone could not be observed. Thus, it was shown that MOS and GOS containing treatments can exert microbiome-independent effects on the metabolome of Caco-2 cells.

OmicsEV: a tool for comprehensive quality evaluation of omics data tables

SUMMARY

RNA-Seq and mass spectrometry-based studies generate omics data tables with measurements for tens of thousands of genes across all samples in a study. The success of a study relies on the quality of these data tables, which is determined by both experimental data generation and computational methods used to process raw experimental data into quantitative data tables. We present OmicsEV, an R package for the quality evaluation of omics data tables. For each data table, OmicsEV uses a series of methods to evaluate data depth, data normalization, batch effect, biological signal, platform reproducibility and multi-omics concordance, producing comprehensive visual and quantitative evaluation results that help assess the data quality of individual data tables and facilitate the identification of the optimal data processing method and parameters for the omics study under investigation.

AVAILABILITY AND IMPLEMENTATION

The source code and the user manual of OmicsEV are available at https://github.com/bzhanglab/OmicsEV, and the source code is released under the GPL-3 license.

Advances, obstacles, and opportunities for machine learning in proteomics

The fields of proteomics and machine learning are both large disciplines, each producing well over 5,000 publications per year. However, studies combining both fields are still relatively rare, with only about 2% of recent proteomics papers including machine learning. This review, which focuses on the intersection of the fields, is intended to inspire proteomics researchers to develop skills and knowledge in the application of machine learning. A brief tutorial introduction to machine learning is provided, and research advances that rely on both fields, particularly as they relate to proteomics tools development and biomarker discovery, are highlighted. Key knowledge gaps and opportunities for scientific advancement are also enumerated.

A system suitability testing platform for untargeted, high-resolution mass spectrometry

The broad coverage of untargeted metabolomics poses fundamental challenges for the harmonization of measurements along time, even if they originate from the very same instrument. Internal isotopic standards can hardly cover the chemical complexity of study samples. Therefore, they are insufficient for normalizing data a posteriori as done for targeted metabolomics. Instead, it is crucial to verify instrument’s performance a priori, that is, before samples are injected. Here, we propose a system suitability testing platform for time-of-flight mass spectrometers independent of liquid chromatography. It includes a chemically defined quality control mixture, a fast acquisition method, software for extracting ca. 3,000 numerical features from profile data, and a simple web service for monitoring. We ran a pilot for 21 months and present illustrative results for anomaly detection or learning causal relationships between the spectral features and machine settings. Beyond mere detection of anomalies, our results highlight several future applications such as 1) recommending instrument retuning strategies to achieve desired values of quality indicators, 2) driving preventive maintenance, and 3) using the obtained, detailed spectral features for posterior data harmonization.

Plasma Proteomics Enable Differentiation of Lung Adenocarcinoma from Chronic Obstructive Pulmonary Disease (COPD)

Chronic obstructive pulmonary disease (COPD) is a major risk factor for the development of lung adenocarcinoma (AC). AC often develops on underlying COPD; thus, the differentiation of both entities by biomarker is challenging. Although survival of AC patients strongly depends on early diagnosis, a biomarker panel for AC detection and differentiation from COPD is still missing. Plasma samples from 176 patients with AC with or without underlying COPD, COPD patients, and hospital controls were analyzed using mass-spectrometry-based proteomics. We performed univariate statistics and additionally evaluated machine learning algorithms regarding the differentiation of AC vs. COPD and AC with COPD vs. COPD. Univariate statistics revealed significantly regulated proteins that were significantly regulated between the patient groups. Furthermore, random forest classification yielded the best performance for differentiation of AC vs. COPD (area under the curve (AUC) 0.935) and AC with COPD vs. COPD (AUC 0.916). The most influential proteins were identified by permutation feature importance and compared to those identified by univariate testing. We demonstrate the great potential of machine learning for differentiation of highly similar disease entities and present a panel of biomarker candidates that should be considered for the development of a future biomarker panel.

Anti-TNF Therapies Suppress Adipose Tissue Inflammation in Crohn’s Disease

Anti-TNF biologics have been shown to markedly improve the quality of life for patients with Crohn’s disease (CD), yet one-third of patients fail to benefit from this treatment. Patients with CD develop a characteristic wrapping of visceral adipose tissue (VAT) in the inflamed intestinal area, termed creeping fat, and it is known that adipose tissue expansion influences the efficacy of anti-TNF drugs. We questioned whether anti-TNF therapies impact the creeping fat in CD, which might affect the outcome of the disease. Adipose tissue biopsies were obtained from a cohort of 14 patients with CD that received anti-TNF drugs and from 29 non-anti-TNF-treated patients (control group) matched by sex, age, and body mass index undergoing surgical interventions for symptomatic complications. We found that anti-TNF therapies restored adipose tissue morphology and suppressed immune cell infiltration in the creeping fat. Additionally, anti-TNF treatments appeared to markedly improve the pro-inflammatory phenotype of adipose-tissue macrophages and adipose-tissue-derived stem cells. Our study provides evidence that anti-TNF medications influence immune cells and progenitor cells in the creeping of patients with CD, suppressing inflammation. We propose that perilesional VAT should be considered when administering anti-TNF therapy in patients with CD.

Elevated ganglioside GM2 activator (GM2A) in human brain tissue reduces neurite integrity and spontaneous neuronal activity

BACKGROUND

Alzheimer's Disease (AD) affects millions globally, but therapy development is lagging. New experimental systems that monitor neuronal functions in conditions approximating the AD brain may be beneficial for identifying new therapeutic strategies.

METHODS

We expose cultured neurons to aqueous-soluble human brain extract from 43 individuals across a spectrum of AD pathology. Multi-electrode arrays (MEAs) and live-cell imaging were used to assess neuronal firing and neurite integrity (NI), respectively, following treatments of rat cortical neurons (MEA) and human iPSC-derived neurons (iN) with human brain extracts.

RESULTS

We observe associations between spontaneous activity and Aβ42:40 levels, between neurite integrity and oligomeric Aβ, and between neurite integrity and tau levels present in the brain extracts. However, these associations with Aβ and tau do not fully account for the effects observed. Proteomic profiling of the brain extracts revealed additional candidates correlated with neuronal structure and activity. Neurotoxicity in MEA and NI assays was associated with proteins implicated in lysosomal storage disorders, while neuroprotection was associated with proteins of the WAVE regulatory complex controlling actin cytoskeleton dynamics. Elevated ganglioside GM2 activator (GM2A) associates with reductions in both NI and MEA activity, and cell-derived GM2A alone is sufficient to induce a loss of neurite integrity and a reduction in neuronal firing.

CONCLUSIONS

The techniques and data herein introduce a system for modeling neuronal vulnerability in response to factors in the human brain and provide insights into proteins potentially contributing to AD pathogenesis.