Publications

Spatial proteomics offers unprecedented insights into the localization, quantity, and interactions of proteins within cells and tissues, thereby enabling researchers and clinicians to map protein networks with high precision. This powerful approach is currently revolutionizing our understanding of cellular organization and function (Method of the Year 2024: spatial proteomics, 2024), and providing deeper insights into how protein distribution and dynamics contribute to physiological and pathological processes. Here, we discuss recent advances in targeted and untargeted (exploratory) spatial proteomics, shed light on emerging multiscale approaches for tissue profiling, and highlight their translational potential.

Multiplexed imaging is a powerful approach in spatial biology, although it is complex, expensive and labor-intensive. Here, we present the IBEX Knowledge-Base, a central resource for reagents, protocols and more, to enhance knowledge sharing, optimization and innovation of spatial proteomics techniques.

Human papilloma virus-negative head and neck squamous cell carcinoma (HNSCC) frequently harbors 11q13 amplifications. Among the oncogenes at this locus, CCND1 and ANO1 are linked to poor prognosis; however, their individual roles in treatment resistance remain unclear. The impact of Cyclin D1 and Ano1 overexpression on survival was analyzed using the TCGA HNSCC dataset and a Charité cohort treated with cisplatin (CDDP)-based radiochemotherapy. High Ano1 expression was primarily associated with poor overall survival in both datasets. The effects of CCND1 and ANO1 knockdown (KD) on radio- and drug sensitivity, along with changes in global protein expression, cell viability, growth, and DNA repair, were studied in an 11q13-amplified HNSCC cell line model of primary cisplatin resistance. Unique pathway alterations- VEGF in CCND1 KD and the Rho GTPase cycle in ANO1 KD- were observed, along with shared changes like DNA damage and cell cycle dysregulation. Silencing CCND1 or ANO1 increased CDDP sensitivity, while only ANO1 silencing increased radiosensitivity. Copanlisib and afatinib were identified as promising candidates for combination therapy of 11q13-amplified HNSCC tumors. We demonstrated a predominant role for Ano1 in treatment resistance in Cyclin D1highAno1high HNSCC tumors and identified novel potential treatment combinations for this high-risk patient group.

Tumors are complex ecosystems composed of malignant and non-malignant cells embedded in a dynamic extracellular matrix (ECM). In the tumor microenvironment, molecular phenotypes are controlled by cell-cell and ECM interactions in 3D cellular neighborhoods (CNs). While their inhibition can impede tumor progression, routine molecular tumor profiling fails to capture cellular interactions. Single-cell spatial transcriptomics (ST) maps receptor-ligand interactions but usually remains limited to 2D tissue sections and lacks ECM readouts. Here, we integrate 3D ST with ECM imaging in serial sections from one clinical lung carcinoma to systematically quantify molecular states, cell-cell interactions, and ECM remodeling in CN. Our integrative analysis pinpointed known immune escape and tumor invasion mechanisms, revealing several druggable drivers of tumor progression in the patient under study. This proof-of-principle study highlights the potential of in-depth CN profiling in routine clinical samples to inform microenvironment-directed therapies. A record of this paper's transparent peer review process is included in the supplemental information.

Muscle-invasive bladder cancer (MIBC) presents significant treatment challenges. Antibody-drug conjugates targeting human epidermal growth factor receptor 2 (HER2), trophoblast cell surface antigen 2 (TROP-2), and nectin cell adhesion molecule 4 (NECTIN4) offer promising therapeutic options. This study examined the spatial expression of HER2, TROP-2, and NECTIN4 in MIBC and metastases, their association with molecular subtypes, and clinical outcomes. Formalin-fixed, paraffin-embedded tissue samples from 251 patients with MIBC were analyzed using immunohistochemistry and tissue microarray analysis. Expression patterns between the tumor front and tumor center (TC) were compared, and statistical analyses assessed associations with molecular subtypes and clinical parameters. Additionally, 67 matched lymph node metastases and a secondary cohort comprising 16 distant metastases, including 7 matched primary tumors, were examined to explore the expression patterns in advanced tumor stages. In primary tumors, HER2 was predominantly negative (83%) but showed higher positivity in the TC. TROP-2 exhibited high overall expression (58% score 3+), whereas NECTIN4 displayed significant heterogeneity with stronger expression in the TC. Spatial overexpression of TROP-2 and NECTIN4 at the tumor front relative to the TC was associated with a better disease-free survival. Accurate assessment required 4 biopsies for HER2 and NECTIN4 and 3 for TROP-2. HER2 expression was associated with urothelial-like and genomically unstable molecular subtypes, whereas TROP-2 was widely expressed except in the mesenchymal-like subtype. NECTIN4 showed the absence of staining in basal, Mes-like, and Sc/NEC-like subtypes. Paired lymph node metastases showed higher expression scores for all 3 markers, whereas distant metastases showed reduced NECTIN4 expression. Additionally, lymph node metastases revealed a considerable heterogeneity for HER2 compared with their matched primary tumors. The spatial heterogeneity of HER2, TROP-2, and NECTIN4 expression necessitates multiple biopsies, particularly from the TC, for accurate evaluation. These findings underscore the need for personalized treatment strategies in MIBC, considering the increased risk of relapse associated with HER2 and NECTIN4 overexpression in the TC. Implementing a multibiopsy approach is critical to enhancing diagnostic accuracy.

Background: The study of circulating blood proteins in population cohorts offers new avenues to explore lifestyle-related and genetic influences describing and shaping human health.
Methods: Utilizing high-throughput mass spectrometry, we quantified 148 highly abundant proteins, functioning in the innate and adaptive immune system, coagulation and nutrient transport in 3632 blood plasma, and 500 serum samples from the CHRIS and BASE-II cross-sectional population studies, respectively. Through multiple regression analyses, we aimed to identify the main factors influencing the circulating proteome at population level.
Results: Many demographic covariates and common medications affect the concentration of high-abundant plasma proteins, but the most significant changes are linked to the use of hormonal contraceptives (HCU). HCU particularly alters amongst others the levels of Angiotensinogen and Transcortin. We robustly replicated these findings in the BASE-II cohort. Furthermore, our results indicate that combined hormonal contraceptives with ethinylestradiol have a stronger effect compared to bioidentical estrogens. Our analysis detects no lasting impact of hormonal contraceptives on the plasma proteome.
Conclusions: HCU is the dominant factor reshaping the high-abundant circulating blood proteome in two population studies. Given the high prevalence of HCU among young women, it is essential to account for this treatment in human proteome studies to avoid misinterpreting its impact as sex- or age-related effects. Although we did not investigate the influence of HCU-induced proteomic changes on human health, our data suggest that future studies on this topic are warranted.

Background: Biological age reflects inter-individual differences in biological function and capacity beyond chronological age. DNA methylation age (DNAmA) and its deviation from chronological age, DNAmA acceleration (DNAmAA), which was calculated as residuals of leukocyte cell count adjusted linear regression of DNAmA on chronological age, were used to estimate biological age in this study. Low levels of serum selenium, selenoprotein P (SELENOP), and the selenocysteine-containing glutathione peroxidase 3 (GPx3) are associated with adverse health outcomes and selenium supplementation is discussed as an anti-aging intervention.
Methods: In this study, we cross-sectionally analyzed 1568 older participants from the observational Berlin Aging Study II (mean age ± SD: 68.8 ± 3.7 years, 51% women). Serum selenium was measured by total reflection X-ray fluorescence (TXRF) spectroscopy and SELENOP was determined by sandwich ELISA. GPx3 was assessed as part of a proteomics dataset using liquid chromatography-mass spectrometry (LC-MS). The relationship between selenium biomarkers and epigenetic clock measures was analyzed using linear regression analyses. P values and 95% confidence intervals (not adjusted for multiple testing) are stated for each analysis.
Results: Participants with deficient serum selenium levels (< 90 μg/L) had a higher rate of biological aging (DunedinPACE, β = - 0.02, SE = 0.01, 95% CI - 0.033 to - 0.004, p = 0.010, n = 865). This association remained statistically significant after adjustment for age, sex, BMI, smoking, and first four genetic principal components (β = - 0.02, SE = 0.01, 95% CI - 0.034 to - 0.004, p = 0.012, n = 757). Compared to the highest quartile, participants in the lowest quartile of SELENOP levels showed an accelerated biological aging rate (DunedinPACE, β = - 0.03, SE = 0.01, 95% CI - 0.051 to - 0.008, p = 0.007, n = 740, fully adjusted model). Similarly, after adjustment for confounders, accelerated biological age was found in participants within the lowest GPx3 quartile compared to participants in the fourth quartile (DunedinPACE, β = - 0.04, SE = 0.01, 95% CI - 0.06 to - 0.02, p = 0.001, n = 674 and GrimAge, β = - 0.98, SE = 0.32, 95% CI - 1.6 to - 0.4, p = 0.002, n = 608). Only the association with GPx3 remained statistically significant after multiple testing correction.
Conclusion: Our study suggests that low levels of selenium biomarkers are associated with accelerated biological aging measured through epigenetic clocks. This effect was not substantially changed after adjustment for known confounders.

Advancing MS-based proteomics toward clinical applications evolves around developing standardized start-to-finish and fit-for-purpose workflows for clinical specimens. Steps along the method design involve the determination and optimization of several bioanalytical parameters such as selectivity, sensitivity, accuracy, and precision. In a joint effort, eight proteomics laboratories belonging to the MSCoreSys initiative including the CLINSPECT-M, MSTARS, DIASyM, and SMART-CARE consortia performed a longitudinal round-robin study to assess the analysis performance of plasma and serum as clinically relevant samples. A variety of LC-MS/MS setups including mass spectrometer models from ThermoFisher and Bruker as well as LC systems from ThermoFisher, Evosep, and Waters Corporation were used in this study. As key performance indicators, sensitivity, precision, and reproducibility were monitored over time. Protein identifications range between 300 and 400 IDs across different state-of-the-art MS instruments, with timsTOF Pro, Orbitrap Exploris 480, and Q Exactive HF-X being among the top performers. Overall, 71 proteins are reproducibly detectable in all setups in both serum and plasma samples, and 22 of these proteins are FDA-approved biomarkers, which are reproducibly quantified (CV < 20% with label-free quantification). In total, the round-robin study highlights a promising baseline for bringing MS-based measurements of serum and plasma samples closer to clinical utility.

Aims: The gastrointestinal (GI) tract is composed of distinct sub-regions, which exhibit segment-specific differences in microbial colonization and (patho)physiological characteristics. Gut microbes can be collectively considered as an active endocrine organ. Microbes produce metabolites, which can be taken up by the host and can actively communicate with the immune cells in the gut lamina propria with consequences for cardiovascular health. Variation in bacterial load and composition along the GI tract may influence the mucosal microenvironment and thus be reflected its interstitial fluid (IF). Characterization of the segment-specific microenvironment is challenging and largely unexplored because of lack of available tools.
Methods and results: Here, we developed methods, namely tissue centrifugation and elution, to collect IF from the mucosa of different intestinal segments. These methods were first validated in rats and mice, and the tissue elution method was subsequently translated for use in humans. These new methods allowed us to quantify microbiota-derived metabolites, mucosa-derived cytokines, and proteins at their site-of-action. Quantification of short-chain fatty acids showed enrichment in the colonic IF. Metabolite and cytokine analyses revealed differential abundances within segments, often significantly increased compared to plasma, and proteomics revealed that proteins annotated to the extracellular phase were site-specifically identifiable in IF. Lipopolysaccharide injections in rats showed significantly higher ileal IL-1β levels in IF compared to the systemic circulation, suggesting the potential of local as well as systemic effect.
Conclusion: Collection of IF from defined segments and the direct measurement of mediators at the site-of-action in rodents and humans bypasses the limitations of indirect analysis of faecal samples or serum, providing direct insight into this understudied compartment.

Most high-grade serous ovarian cancers (OC) do not respond to current immunotherapies. To identify potential new actionable tumor antigens in OC, we performed immunopeptidomics on a human OC cell line expressing the HLA-A02:01 haplotype, which is commonly expressed across many racial and ethnic groups. From this dataset, we identified TTLL8, POTEE, and PKMYT1 peptides as candidate tumor antigens with low expression in normal tissues and upregulated expression in OC. Using tissue microarrays, we assessed the protein expression of TTLL8 and POTEE and their association with patient outcomes in a large cohort of OC patients. TTLL8 was found to be expressed in 56.7% of OC and was associated with a worse overall prognosis. POTEE was expressed in 97.2% of OC patients and had no significant association with survival. In patient TILs, increases in cytokine production and tetramer-positive populations identified antigen-specific CD8 T cell responses, which were dependent on antigen presentation by HLA class I. Antigen-specific T cells triggered cancer cell killing of antigen-pulsed OC cells. These findings suggest that TTLL8, POTEE, and PKMYT1 are potential targets for the development of antigen-targeted immunotherapy in OC.

Advances in single-cell proteomics allow quantification of half of the expressed proteome in an individual cell.

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

Spatial tissue proteomics integrating whole-slide imaging, laser microdissection, and ultrasensitive mass spectrometry is a powerful approach to link cellular phenotypes to functional proteome states in (patho)physiology. To be applicable to large patient cohorts and low sample input amounts, including single-cell applications, loss-minimized and streamlined end-to-end workflows are key. We here introduce an automated sample preparation protocol for laser microdissected samples utilizing the cellenONE robotic system, which has the capacity to process 192 samples in 3 h. Following laser microdissection collection directly into the proteoCHIP LF 48 or EVO 96 chip, our optimized protocol facilitates lysis, formalin de-crosslinking, and tryptic digest of low-input archival tissue samples. The seamless integration with the Evosep ONE LC system by centrifugation allows 'on-the-fly' sample clean-up, particularly pertinent for laser microdissection workflows. We validate our method in human tonsil archival tissue, where we profile proteomes of spatially-defined B-cell, T-cell, and epithelial microregions of 4000 μm2 to a depth of ∼2000 proteins and with high cell type specificity. We finally provide detailed equipment templates and experimental guidelines for broad accessibility.

Background: The aim of this study was to analyse transcriptomic differences between primary and recurrent high-grade serous ovarian carcinoma (HGSOC) to identify prognostic biomarkers.
Methods: We analysed 19 paired primary and recurrent HGSOC samples using targeted RNA sequencing. We selected the best candidates using in silico survival and pathway analysis and validated the biomarkers using immunohistochemistry on a cohort of 44 paired samples, an additional cohort of 504 primary HGSOCs and explored their function.
Results: We identified 233 differential expressed genes. Twenty-three showed a significant prognostic value for PFS and OS in silico. Seven markers (AHRR, COL5A2, FABP4, HMGCS2, ITGA5, SFRP2 and WNT9B) were chosen for validation at the protein level. AHRR expression was higher in primary tumours (p < 0.0001) and correlated with better patient survival (p < 0.05). Stromal SFRP2 expression was higher in recurrent samples (p = 0.009) and protein expression in primary tumours was associated with worse patient survival (p = 0.022). In multivariate analysis, tumour AHRR and SFRP2 remained independent prognostic markers. In vitro studies supported the anti-tumorigenic role of AHRR and the oncogenic function of SFRP2.
Conclusions: Our results underline the relevance of AHRR and SFRP2 proteins in aryl-hydrocarbon receptor and Wnt-signalling, respectively, and might lead to establishing them as biomarkers in HGSOC.

Objectives: Thyroid transcription factor 1 (TTF-1) is a well-established independent prognostic factor in lung adenocarcinoma (LUAD), irrespective of stage. This study aims to determine if TTF-1's prognostic impact is solely based on histomorphological differentiation (tumor grading) or if it independently relates to a biologically more aggressive phenotype. We analyzed a large bi-centric LUAD cohort to accurately assess TTF-1's prognostic value in relation to tumor grade.
Patients and methods: We studied 447 patients with resected LUAD from major German lung cancer centers (Berlin and Cologne), correlating TTF-1 status and grading with clinical, pathologic, and molecular data, alongside patient outcomes. TTF-1's impact was evaluated through univariate and multivariate Cox regression. Causal graph analysis was used to identify and account for potential confounders, improving the statistical estimation of TTF-1's predictive power for clinical outcomes.
Results: Univariate analysis revealed TTF-1 positivity associated with significantly longer disease-free survival (DFS) (median log HR -0.83; p = 0.018). Higher tumor grade showed a non-significant association with shorter DFS (median log HR 0.30; p = 0,62 for G1 to G2 and 0.68; p = 0,34 for G2 to G3). In multivariate analysis, TTF-1 positivity resulted in a significantly longer DFS (median log HR -0.65; p = 0.05) independent of all other parameters, including grading. Adjusting for potential confounders as indicated by the causal graph confirmed the superiority of TTF-1 over tumor grading in prognostics power.
Conclusions: TTF-1 status predicts relapse and survival in LUAD independently of tumor grading. The prognostic power of tumor grading is limited to TTF-1-positive patients, and the effect size of TTF-1 surpasses that of tumor grading. We recommend including TTF1 status as a prognostic factor in the diagnostic guidelines of LUAD.

Motivation: High dimensional single-cell mass cytometry data are confounded by unwanted covariance due to variations in cell size and staining efficiency, making analysis, and interpretation challenging.
Results: We present RUCova, a novel method designed to address confounding factors in mass cytometry data. RUCova removes unwanted covariance from measured markers applying multivariate linear regression based on surrogates of sources of unwanted covariance (SUCs) and principal component analysis (PCA). We exemplify the use of RUCova and show that it effectively removes unwanted covariance while preserving genuine biological signals. Our results demonstrate the efficacy of RUCova in elucidating complex data patterns, facilitating the identification of activated signalling pathways, and improving the classification of important cell populations such as apoptotic cells. By providing a robust framework for data normalization and interpretation, RUCova enhances the accuracy and reliability of mass cytometry analyses, contributing to advances in our understanding of cellular biology and disease mechanisms.

Cells that lack p53 signaling frequently occur in ulcerative colitis (UC) and are considered early drivers in UC-associated colorectal cancer (CRC). Epithelial injury during colitis is associated with transient stem cell reprogramming from the adult, homeostatic to a "fetal-like" regenerative state. Here, we use murine and organoid-based models to study the role of Trp53 during epithelial reprogramming. We find that p53 signaling is silent and dispensable during homeostasis but strongly up-regulated in the epithelium upon DSS-induced colitis. While in WT cells this causes termination of the regenerative state, crypts that lack Trp53 remain locked in the highly proliferative, regenerative state long-term. The regenerative state in WT cells requires high Wnt signaling to maintain elevated levels of glycolysis. Instead, Trp53 deficiency enables Wnt-independent glycolysis due to overexpression of rate-limiting enzyme PKM2. Our study reveals the context-dependent relevance of p53 signaling specifically in the injury-induced regenerative state, explaining the high abundance of clones lacking p53 signaling in UC and UC-associated CRC.

Acute myeloid leukemia (AML) is a heterogeneous disease characterized by genomic aberrations in oncogenes, cytogenetic abnormalities, and an aberrant epigenetic landscape. Nearly 50% of AML cases will relapse with current treatment. A major source of therapy resistance is the interaction of mesenchymal stroma with leukemic cells resulting in therapeutic protection. We aimed to determine pro-survival/anti-apoptotic protein networks involved in the stroma protection of leukemic cells. Proteomic profiling of cultured primary AML (n = 14) with Hs5 stroma cell line uncovered an up-regulation of energy-favorable metabolic proteins. Next, we modulated stroma-induced drug resistance with an epigenetic drug library, resulting in reduced apoptosis with histone deacetylase inhibitor (HDACi) treatment versus other epigenetic modifying compounds. Quantitative phosphoproteomic probing of this effect further revealed a metabolic-enriched phosphoproteome including significant up-regulation of acetyl-coenzyme A synthetase (ACSS2, S30) in leukemia-stroma HDACi treated cocultures compared with untreated monocultures. Validating these findings, we show ACSS2 substrate, acetate, promotes leukemic proliferation, ACSS2 knockout in leukemia cells inhibits leukemic proliferation and ACSS2 knockout in the stroma impairs leukemic metabolic fitness. Finally, we identify ACSS1/ACSS2-high expression AML subtype correlating with poor overall survival. Collectively, this study uncovers the leukemia-stroma phosphoproteome emphasizing a role for ACSS2 in mediating AML growth and drug resistance.

Salt sensitivity concerns blood pressure alterations after a change in salt intake (sodium chloride). The heart is a pump, and vessels are tubes; sodium can affect both. A high salt intake increases cardiac output, promotes vascular dysfunction and capillary rarefaction, and chronically leads to increased systemic vascular resistance. More recent findings suggest that sodium also acts as an important second messenger regulating energy metabolism and cellular functions. Besides endothelial cells and fibroblasts, sodium also affects innate and adaptive immunometabolism, immune cell function, and influences certain microbes and microbiota-derived metabolites. We propose the idea that the definition of salt sensitivity should be expanded beyond high blood pressure to cellular and molecular salt sensitivity.

Aneuploidy is a hallmark of human cancer, yet the molecular mechanisms to cope with aneuploidy-induced cellular stresses remain largely unknown. Here, we induce chromosome mis-segregation in non-transformed RPE1-hTERT cells and derive multiple stable clones with various degrees of aneuploidy. We perform a systematic genomic, transcriptomic and proteomic profiling of 6 isogenic clones, using whole-exome DNA, mRNA and miRNA sequencing, as well as proteomics. Concomitantly, we functionally interrogate their cellular vulnerabilities, using genome-wide CRISPR/Cas9 and large-scale drug screens. Aneuploid clones activate the DNA damage response and are more resistant to further DNA damage induction. Aneuploid cells also exhibit elevated RAF/MEK/ERK pathway activity and are more sensitive to clinically-relevant drugs targeting this pathway, and in particular to CRAF inhibition. Importantly, CRAF and MEK inhibition sensitize aneuploid cells to DNA damage-inducing chemotherapies and to PARP inhibitors. We validate these results in human cancer cell lines. Moreover, resistance of cancer patients to olaparib is associated with high levels of RAF/MEK/ERK signaling, specifically in highly-aneuploid tumors. Overall, our study provides a comprehensive resource for genetically-matched karyotypically-stable cells of various aneuploidy states, and reveals a therapeutically-relevant cellular dependency of aneuploid cells.

B cell receptor (BCR) signaling is required for the survival and maturation of B cells and is deregulated in B cell lymphomas. While proximal BCR signaling is well studied, little is known about the crosstalk of downstream effector pathways, and a comprehensive quantitative network analysis of BCR signaling is missing. Here, we semi-quantitatively modelled BCR signaling in Burkitt lymphoma (BL) cells using systematically perturbed phosphorylation data of BL-2 and BL-41 cells. The models unveiled feedback and crosstalk structures in the BCR signaling network, including a negative crosstalk from p38 to MEK/ERK. The relevance of the crosstalk was verified for BCR and CD40 signaling in different BL cells and confirmed by global phosphoproteomics on ERK itself and known ERK target sites. Compared to the starting network, the trained network for BL-2 cells was better transferable to BL-41 cells. Moreover, the BL-2 network was also suited to model BCR signaling in Diffuse large B cell lymphoma cells lines with aberrant BCR signaling (HBL-1, OCI-LY3), indicating that BCR aberration does not cause a major downstream rewiring.

Multiple myeloma (MM) is a plasma cell malignancy of the bone marrow. Despite therapeutic advances, MM remains incurable, and better risk stratification as well as new therapies are therefore highly needed. The proteome of MM has not been systematically assessed before and holds the potential to uncover insight into disease biology and improved prognostication in addition to genetic and transcriptomic studies. Here we provide a comprehensive multiomics analysis including deep tandem mass tag-based quantitative global (phospho)proteomics, RNA sequencing, and nanopore DNA sequencing of 138 primary patient-derived plasma cell malignancies encompassing treatment-naive MM, plasma cell leukemia and the premalignancy monoclonal gammopathy of undetermined significance, as well as healthy controls. We found that the (phospho)proteome of malignant plasma cells are highly deregulated as compared with healthy plasma cells and is both defined by chromosomal alterations as well as posttranscriptional regulation. A prognostic protein signature was identified that is associated with aggressive disease independent of established risk factors in MM. Integration with functional genetics and single-cell RNA sequencing revealed general and genetic subtype-specific deregulated proteins and pathways in plasma cell malignancies that include potential targets for (immuno)therapies. Our study demonstrates the potential of proteogenomics in cancer and provides an easily accessible resource for investigating protein regulation and new therapeutic approaches in MM.

Accessing the natural genetic diversity of species unveils hidden genetic traits, clarifies gene functions and allows the generalizability of laboratory findings to be assessed. One notable discovery made in natural isolates of Saccharomyces cerevisiae is that aneuploidy—an imbalance in chromosome copy numbers—is frequent (around 20%), which seems to contradict the substantial fitness costs and transient nature of aneuploidy when it is engineered in the laboratory. Here we generate a proteomic resource and merge it with genomic and transcriptomic data for 796 euploid and aneuploid natural isolates. We find that natural and lab-generated aneuploids differ specifically at the proteome. In lab-generated aneuploids, some proteins—especially subunits of protein complexes—show reduced expression, but the overall protein levels correspond to the aneuploid gene dosage. By contrast, in natural isolates, more than 70% of proteins encoded on aneuploid chromosomes are dosage compensated, and average protein levels are shifted towards the euploid state chromosome-wide. At the molecular level, we detect an induction of structural components of the proteasome, increased levels of ubiquitination, and reveal an interdependency of protein turnover rates and attenuation. Our study thus highlights the role of protein turnover in mediating aneuploidy tolerance, and shows the utility of exploiting the natural diversity of species to attain generalizable molecular insights into complex biological processes.

Increased throughput in proteomic experiments can improve accessibility of proteomic platforms, reduce costs, and facilitate new approaches in systems biology and biomedical research. Here we propose combination of analytical flow rate chromatography with ion mobility separation of peptide ions, data-independent acquisition, and data analysis with the DIA-NN software suite, to achieve high-quality proteomic experiments from limited sample amounts, at a throughput of up to 400 samples per day. For instance, when benchmarking our workflow using a 500-μL/min flow rate and 3-min chromatographic gradients, we report the quantification of 5211 proteins from 2 μg of a mammalian cell-line standard at high quantitative accuracy and precision. We further used this platform to analyze blood plasma samples from a cohort of COVID-19 inpatients, using a 3-min chromatographic gradient and alternating column regeneration on a dual pump system. The method delivered a comprehensive view of the COVID-19 plasma proteome, allowing classification of the patients according to disease severity and revealing plasma biomarker candidates.

Protein glycosylation, a complex and heterogeneous post-translational modification that is frequently dysregulated in disease, has been difficult to analyse at scale. Here we report a data-independent acquisition technique for the large-scale mass-spectrometric quantification of glycopeptides in plasma samples. The technique, which we named ‘OxoScan-MS’, identifies oxonium ions as glycopeptide fragments and exploits a sliding-quadrupole dimension to generate comprehensive and untargeted oxonium ion maps of precursor masses assigned to fragment ions from non-enriched plasma samples. By applying OxoScan-MS to quantify 1,002 glycopeptide features in the plasma glycoproteomes from patients with COVID-19 and healthy controls, we found that severe COVID-19 induces differential glycosylation in IgA, haptoglobin, transferrin and other disease-relevant plasma glycoproteins. OxoScan-MS may allow for the quantitative mapping of glycoproteomes at the scale of hundreds to thousands of samples.

Despite abundant genomic and phenotypic data across individuals and environments, the functional impact of most mutations on phenotype remains unclear. Here, we bridge this gap by linking genome to proteome in 800 meiotic progeny from an intercross between two closely related Saccharomyces cerevisiae isolates adapted to distinct niches. Modest genetic distance between the parents generated remarkable proteomic diversity that was amplified in the progeny and captured by 6,476 genotype-protein associations, over 1,600 of which we resolved to single variants. Proteomic adaptation emerged through the combined action of numerous cis- and trans-regulatory mutations, a regulatory architecture that was conserved across the species. Notably, trans-regulatory variants often arose in proteins not traditionally associated with gene regulation, such as enzymes. Moreover, the proteomic consequences of mutations predicted fitness under various stresses. Our study demonstrates that the collective action of natural genetic variants drives dramatic proteome diversification, with molecular consequences that forecast phenotypic outcomes.

Data-independent acquisition (DIA) is increasingly preferred over data-dependent acquisition due to its higher throughput and fewer missing values. Whereas data-dependent acquisition often uses stable isotope labeling to improve quantification, DIA mostly relies on label-free approaches. Efforts to integrate DIA with isotope labeling include chemical methods like mass differential tags for relative and absolute quantification and dimethyl labeling, which, while effective, complicate sample preparation. Stable isotope labeling by amino acids in cell culture (SILAC) achieves high labeling efficiency through the metabolic incorporation of heavy labels into proteins in vivo. However, the need for metabolic incorporation limits the direct use in clinical scenarios and certain high-throughput experiments. Spike-in SILAC (SiS) methods use an externally generated heavy sample as an internal reference, enabling SILAC-based quantification even for samples that cannot be directly labeled. Here, we combine DIA-SiS, leveraging the robust quantification of SILAC without the complexities associated with chemical labeling. We developed DIA-SiS and rigorously assessed its performance with mixed-species benchmark samples on bulk and single cell-like amount level. We demonstrate that DIA-SiS substantially improves proteome coverage and quantification compared to label-free approaches and reduces incorrectly quantified proteins. Additionally, DIA-SiS proves effective in analyzing proteins in low-input formalin-fixed paraffin-embedded tissue sections. DIA-SiS combines the precision of stable isotope-based quantification with the simplicity of label-free sample preparation, facilitating simple, accurate, and comprehensive proteome profiling.

Glyceraldehyde (GA) is a three-carbon monosaccharide that can be present in cells as a by-product of fructose metabolism. Bruno Mendel and Otto Warburg showed that the application of GA to cancer cells inhibits glycolysis and their growth. However, the molecular mechanism by which this occurred was not clarified. We describe a novel multi-modal mechanism by which the L-isomer of GA (L-GA) inhibits neuroblastoma cell growth. L-GA induces significant changes in the metabolic profile, promotes oxidative stress and hinders nucleotide biosynthesis. GC-MS and 13C-labeling was employed to measure the flow of carbon through glycolytic intermediates under L-GA treatment. It was found that L-GA is a potent inhibitor of glycolysis due to its proposed targeting of NAD(H)-dependent reactions. This results in growth inhibition, apoptosis and a redox crisis in neuroblastoma cells. It was confirmed that the redox mechanisms were modulated via L-GA by proteomic analysis. Analysis of nucleotide pools in L-GA-treated cells depicted a previously unreported observation, in which nucleotide biosynthesis is significantly inhibited. The inhibitory action of L-GA was partially relieved with the co-application of the antioxidant N-acetyl-cysteine. We present novel evidence for a simple sugar that inhibits cancer cell proliferation via dysregulating its fragile homeostatic environment.

Background
Current treatment for head and neck squamous cell carcinoma (HNSCC) involves surgery, radiotherapy, and chemotherapy. Despite aggressive multimodal approaches, tumour recurrence occurs in 40–60 % of cases, leading to poor survival outcomes. HNSCC lacks common genetic drivers for tailored therapies, and reliable biomarkers for treatment selection are scarce. We investigated the procedural requirements for incorporating drug- and radiosensitivity screens in patient-derived organoids (PDOs) within a clinical trial framework.
Patients and methods
Fresh tumour samples (N = 198) from 186 HNSCC patients were included. Success rates of organoid establishment were correlated with clinical and procedural parameters. Timelines for establishment of PDO cultures were determined, and their long-term growth potential assessed by serial passaging. Additionally, we conducted whole exome sequencing on matched tumour-organoid pairs. Three PDO models were employed to establish radiosensitivity assays.
Results
In total, PDO models displaying histomorphological features and genomic alterations of parental tumours were successfully established for 35 % of patient tumours. Success rates rose to 77 % for samples with a tumour cell content of 30 % or higher. Advanced patient age, prior radiotherapy, and delays in tissue processing were identified as negative predictors for engraftment. The estimated time interval needed for screens was compatible with PDO-guided selection of curative-intent radiotherapy regimens.
Conclusions
Our findings suggest that with high-quality samples and efficient tissue processing, PDO screens can be successfully performed in 77 % of HNSCC patients. Given the procedural challenges involved, future clinical trials aiming to the utility of PDOs for guiding treatment decisions should consider implementing centralised PDO screening.

Spatial proteomics combining microscopy-based cell phenotyping with ultrasensitive mass-spectrometry-based proteomics is an emerging and powerful concept to study cell function and heterogeneity in (patho)physiology. However, optimized workflows that preserve morphological information for phenotype discovery and maximize proteome coverage of few or even single cells from laser microdissected tissue are currently lacking. Here, we report a robust and scalable workflow for the proteomic analysis of ultra-low-input archival material. Benchmarking in murine liver resulted in up to 2,000 quantified proteins from single hepatocyte contours and nearly 5,000 proteins from 50-cell regions. Applied to human tonsil, we profiled 146 microregions including T and B lymphocyte niches and quantified cell-type-specific markers, cytokines, and transcription factors. These data also highlighted proteome dynamics within activated germinal centers, illuminating sites undergoing B cell proliferation and somatic hypermutation. This approach has broad implications in biomedicine, including early disease profiling and drug target and biomarker discovery. A record of this paper’s transparent peer review process is included in the supplemental information.

Protein cysteinyl thiols are susceptible to reduction-oxidation reactions that can influence protein function. Accurate quantification of cysteine oxidation is therefore crucial for decoding protein redox regulation. Here, we present CysQuant, a novel approach for simultaneous quantification of cysteine oxidation degrees and protein abundancies. CysQuant involves light/heavy iodoacetamide isotopologues for differential labeling of reduced and reversibly oxidized cysteines analyzed by data-dependent acquisition (DDA) or data-independent acquisition mass spectrometry (DIA-MS). Using plexDIA with in silico predicted spectral libraries, we quantified an average of 18% cysteine oxidation in Arabidopsis thaliana by DIA-MS, including a subset of highly oxidized cysteines forming disulfide bridges in AlphaFold2 predicted structures. Applying CysQuant to Arabidopsis seedlings exposed to excessive light, we successfully quantified the well-established increased reduction of Calvin-Benson cycle enzymes and discovered yet uncharacterized redox-sensitive disulfides in chloroplastic enzymes. Overall, CysQuant is a highly versatile tool for assessing the cysteine modification status that can be widely applied across various mass spectrometry platforms and organisms.

Scientists employing omics in life science studies face challenges such as the modeling of multiassay studies, recording of all relevant parameters, and managing many samples with their metadata. They must manage many large files that are the results of the assays or subsequent computation. Users with diverse backgrounds, ranging from computational scientists to wet-lab scientists, have dissimilar needs when it comes to data access, with programmatic interfaces being favored by the former and graphical ones by the latter.
We introduce SODAR, the system for omics data access and retrieval. SODAR is a software package that addresses these challenges by providing a web-based graphical user interface for managing multiassay studies and describing them using the ISA (Investigation, Study, Assay) data model and the ISA-Tab file format. Data storage is handled using the iRODS data management system, which handles large quantities of files and substantial amounts of data. SODAR also offers programmable APIs and command-line access for metadata and file storage.
SODAR supports complex omics integration studies and can be easily installed. The software is written in Python 3 and freely available at https://github.com/bihealth/sodar-server under the MIT license.

Mass spectrometry-based proteomics has been rapidly gaining traction as a powerful analytical method both in basic research and translation. While the problem of error control in peptide and protein identification has been addressed extensively, the quality of the resulting quantities remains challenging to evaluate. Here we introduce QuantUMS (Quantification using an Uncertainty Minimising Solution), a machine learning-based method which minimises errors and eliminates bias in peptide and protein quantification by integrating multiple sources of quantitative information. In combination with data-independent acquisition proteomics, QuantUMS boosts accuracy and precision of quantities, as well as reports an uncertainty metric, enabling effective filtering of data for downstream analysis. The algorithm has linear complexity with respect to the number of mass spectrometry acquisitions in the experiment and is thus scalable to infinitely large proteomic experiments. For an easy implementation in a proteomics laboratory, we integrate QuantUMS in our automated DIA-NN software suite.

Functional genomic strategies have become fundamental for annotating gene function and regulatory networks. Here, we combined functional genomics with proteomics by quantifying protein abundances in a genome-scale knockout library in Saccharomyces cerevisiae, using data-independent acquisition mass spectrometry. We find that global protein expression is driven by a complex interplay of (1) general biological properties, including translation rate, protein turnover, the formation of protein complexes, growth rate, and genome architecture, followed by (2) functional properties, such as the connectivity of a protein in genetic, metabolic, and physical interaction networks. Moreover, we show that functional proteomics complements current gene annotation strategies through the assessment of proteome profile similarity, protein covariation, and reverse proteome profiling. Thus, our study reveals principles that govern protein expression and provides a genome-spanning resource for functional annotation.

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is known to contain an active-site cysteine residue undergoing oxidation in response to hydrogen peroxide, leading to rapid inactivation of the enzyme. Here we show that human and mouse cells expressing a GAPDH mutant lacking this redox switch retain catalytic activity but are unable to stimulate the oxidative pentose phosphate pathway and enhance their reductive capacity. Specifically, we find that anchorage-independent growth of cells and spheroids is limited by an elevation of endogenous peroxide levels and is largely dependent on a functional GAPDH redox switch. Likewise, tumour growth in vivo is limited by peroxide stress and suppressed when the GAPDH redox switch is disabled in tumour cells. The induction of additional intratumoural oxidative stress by chemo- or radiotherapy synergized with the deactivation of the GAPDH redox switch. Mice lacking the GAPDH redox switch exhibit altered fatty acid metabolism in kidney and heart, apparently in compensation for the lack of the redox switch. Together, our findings demonstrate the physiological and pathophysiological relevance of oxidative GAPDH inactivation in mammals.

Colorectal cancer progression is intrinsically linked to stepwise deregulation of the intestinal differentiation trajectory. In this process, sequential mutations of APC, KRAS, TP53, and SMAD4 enable oncogenic signaling and establish the hallmarks of cancer. Here, we use mass cytometry of isogenic human colon organoids and patient-derived cancer organoids to capture oncogenic signaling, cell phenotypes, and differentiation states in a high-dimensional single-cell map. We define a differentiation axis in all tumor progression states from normal to cancer. Our data show that colorectal cancer driver mutations shape the distribution of cells along the differentiation axis. In this regard, subsequent mutations can have stem cell promoting or restricting effects. Individual nodes of the cancer cell signaling network remain coupled to the differentiation state, regardless of the presence of driver mutations. We use single-cell RNA sequencing to link the (phospho-)protein signaling network to transcriptomic states with biological and clinical relevance. Our work highlights how oncogenes gradually shape signaling and transcriptomes during tumor progression.

The molecular heterogeneity of cancer cells contributes to the often partial response to targeted therapies and relapse of disease due to the escape of resistant cell populations. While single-cell sequencing has started to improve our understanding of this heterogeneity, it offers a mostly descriptive view on cellular types and states. To obtain more functional insights, we propose scGeneRAI, an explainable deep learning approach that uses layer-wise relevance propagation (LRP) to infer gene regulatory networks from static single-cell RNA sequencing data for individual cells. We benchmark our method with synthetic data and apply it to single-cell RNA sequencing data of a cohort of human lung cancers. From the predicted single-cell networks our approach reveals characteristic network patterns for tumor cells and normal epithelial cells and identifies subnetworks that are observed only in (subgroups of) tumor cells of certain patients. While current state-of-the-art methods are limited by their ability to only predict average networks for cell populations, our approach facilitates the reconstruction of networks down to the level of single cells which can be utilized to characterize the heterogeneity of gene regulation within and across tumors.

Genetically identical cells are known to differ in many physiological parameters such as growth rate and drug tolerance. Metabolic specialization is believed to be a cause of such phenotypic heterogeneity, but detection of metabolically divergent subpopulations remains technically challenging. We developed a proteomics-based technology, termed differential isotope labelling by amino acids (DILAC), that can detect producer and consumer subpopulations of a particular amino acid within an isogenic cell population by monitoring peptides with multiple occurrences of the amino acid. We reveal that young, morphologically undifferentiated yeast colonies contain subpopulations of lysine producers and consumers that emerge due to nutrient gradients. Deconvoluting their proteomes using DILAC, we find evidence for in situ cross-feeding where rapidly growing cells ferment and provide the more slowly growing, respiring cells with ethanol. Finally, by combining DILAC with fluorescence-activated cell sorting, we show that the metabolic subpopulations diverge phenotypically, as exemplified by a different tolerance to the antifungal drug amphotericin B. Overall, DILAC captures previously unnoticed metabolic heterogeneity and provides experimental evidence for the role of metabolic specialization and cross-feeding interactions as a source of phenotypic heterogeneity in isogenic cell populations.

Despite numerous diagnostic and therapeutic advances, pancreatic ductal adenocarcinoma (PDAC) has a high mortality rate, and is the fourth leading cause of cancer death in developing countries. Besides its increasing prevalence, pancreatic malignancies are characterized by poor prognosis. Omics technologies have potential relevance for PDAC assessment but are time-intensive and relatively cost-intensive and limited by tissue heterogeneity. Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) can obtain spatially distinct peptide-signatures and enables tumor classification within a feasible time with relatively low cost. While MALDI-MSI data sets are inherently large, machine learning methods have the potential to greatly decrease processing time. We present a pilot study investigating the potential of MALDI-MSI in combination with neural networks, for classification of pancreatic ductal adenocarcinoma. Neural-network models were trained to distinguish between pancreatic ductal adenocarcinoma and other pancreatic cancer types. The proposed methods are able to correctly classify the PDAC types with an accuracy of up to 86% and a sensitivity of 82%. This study demonstrates that machine learning tools are able to identify different pancreatic carcinoma from complex MALDI data, enabling fast prediction of large data sets. Our results encourage a more frequent use of MALDI-MSI and machine learning in histopathological studies in the future.

Context
Humans respond profoundly to changes in diet, while nutrition and environment have a great impact on population health. It is therefore important to deeply characterize the human nutritional responses.
Objective
Endocrine parameters and the metabolome of human plasma are rapidly responding to acute nutritional interventions such as caloric restriction or a glucose challenge. It is less well understood whether the plasma proteome would be equally dynamic, and whether it could be a source of corresponding biomarkers.
Methods
We used high-throughput mass spectrometry to determine changes in the plasma proteome of i) 10 healthy, young, male individuals in response to 2 days of acute caloric restriction followed by refeeding; ii) 200 individuals of the Ely epidemiological study before and after a glucose tolerance test at 4 time points (0, 30, 60, 120 minutes); and iii) 200 random individuals from the Generation Scotland study. We compared the proteomic changes detected with metabolome data and endocrine parameters.
Results
Both caloric restriction and the glucose challenge substantially impacted the plasma proteome. Proteins responded across individuals or in an individual-specific manner. We identified nutrient-responsive plasma proteins that correlate with changes in the metabolome, as well as with endocrine parameters. In particular, our study highlights the role of apolipoprotein C1 (APOC1), a small, understudied apolipoprotein that was affected by caloric restriction and dominated the response to glucose consumption and differed in abundance between individuals with and without type 2 diabetes.
Conclusion
Our study identifies APOC1 as a dominant nutritional responder in humans and highlights the interdependency of acute nutritional response proteins and the endocrine system.

FOXP3+ regulatory T cells (Tregs) are central for peripheral tolerance, and their deregulation is associated with autoimmunity. Dysfunctional autoimmune Tregs display pro-inflammatory features and altered mitochondrial metabolism, but contributing factors remain elusive. High salt (HS) has been identified to alter immune function and to promote autoimmunity. By investigating longitudinal transcriptional changes of human Tregs, we identified that HS induces metabolic reprogramming, recapitulating features of autoimmune Tregs. Mechanistically, extracellular HS raises intracellular Na+, perturbing mitochondrial respiration by interfering with the electron transport chain (ETC). Metabolic disturbance by a temporary HS encounter or complex III blockade rapidly induces a pro-inflammatory signature and FOXP3 downregulation, leading to long-term dysfunction in vitro and in vivo. The HS-induced effect could be reversed by inhibition of mitochondrial Na+/Ca2+ exchanger (NCLX). Our results indicate that salt could contribute to metabolic reprogramming and that short-term HS encounter perturb metabolic fitness and long-term function of human Tregs with important implications for autoimmunity.

Metabolism is deeply intertwined with aging. Effects of metabolic interventions on aging have been explained with intracellular metabolism, growth control, and signaling. Studying chronological aging in yeast, we reveal a so far overlooked metabolic property that influences aging via the exchange of metabolites. We observed that metabolites exported by young cells are re-imported by chronologically aging cells, resulting in cross-generational metabolic interactions. Then, we used self-establishing metabolically cooperating communities (SeMeCo) as a tool to increase metabolite exchange and observed significant lifespan extensions. The longevity of the SeMeCo was attributable to metabolic reconfigurations in methionine consumer cells. These obtained a more glycolytic metabolism and increased the export of protective metabolites that in turn extended the lifespan of cells that supplied them with methionine. Our results establish metabolite exchange interactions as a determinant of cellular aging and show that metabolically cooperating cells can shape the metabolic environment to extend their lifespan.

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.

Current mass spectrometry methods enable high-throughput proteomics of large sample amounts, but proteomics of low sample amounts remains limited in depth and throughput. To increase the throughput of sensitive proteomics, we developed an experimental and computational framework, called plexDIA, for simultaneously multiplexing the analysis of peptides and samples. Multiplexed analysis with plexDIA increases throughput multiplicatively with the number of labels without reducing proteome coverage or quantitative accuracy. By using three-plex non-isobaric mass tags, plexDIA enables quantification of threefold more protein ratios among nanogram-level samples. Using 1-hour active gradients, plexDIA quantified ~8,000 proteins in each sample of labeled three-plex sets and increased data completeness, reducing missing data more than twofold across samples. Applied to single human cells, plexDIA quantified ~1,000 proteins per cell and achieved 98% data completeness within a plexDIA set while using ~5 minutes of active chromatography per cell. These results establish a general framework for increasing the throughput of sensitive and quantitative protein analysis.

The assimilation, incorporation, and metabolism of sulfur is a fundamental process across all domains of life, yet how cells deal with varying sulfur availability is not well understood. We studied an unresolved conundrum of sulfur fixation in yeast, in which organosulfur auxotrophy caused by deletion of the homocysteine synthase Met17p is overcome when cells are inoculated at high cell density. In combining the use of self-establishing metabolically cooperating (SeMeCo) communities with proteomic, genetic, and biochemical approaches, we discovered an uncharacterized gene product YLL058Wp, herein named Hydrogen Sulfide Utilizing-1 (HSU1). Hsu1p acts as a homocysteine synthase and allows the cells to substitute for Met17p by reassimilating hydrosulfide ions leaked from met17Δ cells into O-acetyl-homoserine and forming homocysteine. Our results show that cells can cooperate to achieve sulfur fixation, indicating that the collective properties of microbial communities facilitate their basic metabolic capacity to overcome sulfur limitation.

Pre-eclampsia (PE) is a syndrome that affects multiple organ systems and is the most severe hypertensive disorder in pregnancy. It frequently leads to preterm delivery, maternal and fetal morbidity and mortality and life-long complications1. We currently lack efficient screening tools2, 3 and early therapies4, 5 to address PE. To investigate the early stages of early onset PE, and identify candidate markers and pathways, we performed spatio-temporal multi-omics profiling of human PE placentae and healthy controls and validated targets in early gestation in a longitudinal clinical cohort. We used a single-nuclei RNA-seq approach combined with spatial proteo- and transcriptomics and mechanistic in vitro signalling analyses to bridge the gap from late pregnancy disease to early pregnancy pathomechanisms. We discovered a key disruption in villous trophoblast differentiation, which is driven by the increase of transcriptional coactivator p300, that ultimately ends with a senescence-associated secretory phenotype (SASP) of trophoblasts. We found a significant increase in the senescence marker activin A in preeclamptic maternal serum in early gestation, before the development of clinical symptoms, indicating a translation of the placental syndrome to the maternal side. Our work describes a new disease progression, starting with a disturbed transition in villous trophoblast differentiation. Our study identifies potential pathophysiology-relevant biomarkers for the early diagnosis of the disease as well as possible targets for interventions, which would be crucial steps toward protecting the mother and child from gestational mortality and morbidity and an increased risk of cardiovascular disease later in life.

Despite the availabilty of imaging-based and mass-spectrometry-based methods for spatial proteomics, a key challenge remains connecting images with single-cell-resolution protein abundance measurements. Here, we introduce Deep Visual Proteomics (DVP), which combines artificial-intelligence-driven image analysis of cellular phenotypes with automated single-cell or single-nucleus laser microdissection and ultra-high-sensitivity mass spectrometry. DVP links protein abundance to complex cellular or subcellular phenotypes while preserving spatial context. By individually excising nuclei from cell culture, we classified distinct cell states with proteomic profiles defined by known and uncharacterized proteins. In an archived primary melanoma tissue, DVP identified spatially resolved proteome changes as normal melanocytes transition to fully invasive melanoma, revealing pathways that change in a spatial manner as cancer progresses, such as mRNA splicing dysregulation in metastatic vertical growth that coincides with reduced interferon signaling and antigen presentation. The ability of DVP to retain precise spatial proteomic information in the tissue context has implications for the molecular profiling of clinical samples.

Microbial communities are composed of cells of varying metabolic capacity, and regularly include auxotrophs that lack essential metabolic pathways. Through analysis of auxotrophs for amino acid biosynthesis pathways in microbiome data derived from >12,000 natural microbial communities obtained as part of the Earth Microbiome Project (EMP), and study of auxotrophic–prototrophic interactions in self-establishing metabolically cooperating yeast communities (SeMeCos), we reveal a metabolically imprinted mechanism that links the presence of auxotrophs to an increase in metabolic interactions and gains in antimicrobial drug tolerance. As a consequence of the metabolic adaptations necessary to uptake specific metabolites, auxotrophs obtain altered metabolic flux distributions, export more metabolites and, in this way, enrich community environments in metabolites. Moreover, increased efflux activities reduce intracellular drug concentrations, allowing cells to grow in the presence of drug levels above minimal inhibitory concentrations. For example, we show that the antifungal action of azoles is greatly diminished in yeast cells that uptake metabolites from a metabolically enriched environment. Our results hence provide a mechanism that explains why cells are more robust to drug exposure when they interact metabolically.

The ubiquitin-proteasome system (UPS) was discovered about 40 years ago and is known to regulate a multitude of cellular processes including protein homeostasis. Ubiquitylated proteins are recognized by downstream effectors, resulting in alterations of protein abundance, activity, or localization. Not surprisingly, the ubiquitylation machinery is dysregulated in numerous diseases, including cancers and neurodegeneration. Mass spectrometry (MS)-based proteomics has emerged as a transformative technology for characterizing protein ubiquitylation in an unbiased fashion. Here, we provide an overview of the different MS-based approaches for studying protein ubiquitylation. We review various methods for enriching and quantifying ubiquitin modifications at the peptide or protein level, outline MS acquisition, and data processing approaches and discuss key challenges. Finally, we examine how MS-based ubiquitinomics can aid both basic biology and drug discovery research.

Background: COVID-19 has highly variable clinical courses. The search for prognostic host factors for COVID-19 outcome is a priority.
Methods: We performed logistic regression for ICU admission against a polygenic score (PGS) for Cystatin C (CyC) production in patients with COVID-19. We analyzed the predictive value of longitudinal plasma CyC levels in an independent cohort of patients hospitalized with COVID-19.
Findings: In four cohorts spanning European and African ancestry populations, we identified a significant association between CyC-production PGS and odds of critical illness (n cases=2,319), with the strongest association captured in the UKB cohort (OR 2.13, 95% CI 1.58-2.87, p=7.12e-7). Plasma proteomics from an independent cohort of hospitalized COVID-19 patients (n cases = 131) demonstrated that CyC production was associated with COVID-specific mortality (p=0.0007).
Interpretation: Our findings suggest that CyC may be useful for stratification of patients and it has functional role in the host response to COVID-19.

The outbreak of a novel coronavirus (SARS-CoV-2) in 2019 led to a worldwide pandemic, which remains an integral part of our lives to this day. Coronavirus disease (COVID-19) is a flu like condition, often accompanied by high fever and respiratory distress. In some cases, conjointly with other co-morbidities, COVID-19 can become severe, leading to lung arrest and even death. Although well-known from a clinical standpoint, the mechanistic understanding of lethal COVID-19 is still rudimentary. Studying the pathology and changes on a molecular level associated with the resulting COVID-19 disease is impeded by the highly infectious nature of the virus and the concomitant sampling challenges. We were able to procure COVID-19 post-mortem lung tissue specimens by our collaboration with the BSL-3 laboratory of the Biobanking and BioMolecular resources Research Infrastructure Austria which we subjected to state-of-the-art quantitative proteomic analysis to better understand the pulmonary manifestations of lethal COVID-19. Lung tissue samples from age-matched non-COVID-19 patients who died within the same period were used as controls. Samples were subjected to parallel accumulation-serial fragmentation combined with data-independent acquisition (diaPASEF) on a timsTOF Pro and obtained raw data was processed using DIA-NN software. Here we report that terminal COVID-19 patients display an increase in inflammation, acute immune response and blood clot formation (with concomitant triggering of fibrinolysis). Furthermore, we describe that COVID-19 diseased lungs undergo severe extracellular matrix restructuring, which was corroborated on the histopathological level. However, although undergoing an injury, diseased lungs seem to have impaired proliferative and tissue repair signalling, with several key kinase-mediated signalling pathways being less active. This might provide a mechanistic link to post-acute sequelae of COVID-19 (PASC; "Long COVID"). Overall, we emphasize the importance of histopathological patient stratification when interpreting molecular COVID-19 data.

Scanning SWATH coupled with normal-flow LC has been recently introduced for high-content, high-throughput proteomics analysis, which requires a relatively large amount of sample injection. Here we established the microflow LC coupled with Scanning SWATH for samples with relatively small quantities. First, we optimized several key parameters of the LC and MS settings, including C18 particle size for the analytical column, LC gradient and flow rate, as well as effective ion accumulation time and isolation window width for MS acquisition. We then compared the optimized Scanning SWATH method with the conventional variable window SWATH (referred to as SWATH) method. Results showed that the total ion chromatogram signals in Scanning SWATH were 10 times higher than that of SWATH, and Scanning SWATH identified 12.2–22.2% more peptides than SWATH. Finally, we employed 120 min Scanning SWATH to acquire the proteomes of 62 formalin-fixed, paraffin-embedded (FFPE) tissue samples from 31 patients with hepatocellular carcinoma (HCC). Altogether, 92 334 peptides and 8516 proteins were quantified. Besides the reported biomarkers, including ANXA2, MCM7, SUOX, and AKR1B10, we identified new potential HCC biomarkers such as CST5, TP53, CEBPB, and E2F4. Taken together, we present an optimal workflow integrating microflow LC and Scanning SWATH that effectively improves the protein identification and quantitation.

Phosphoproteomics routinely quantifies changes in the levels of thousands of phosphorylation sites, but functional analysis of such data remains a major challenge. While databases like PhosphoSitePlus contain information about many phosphorylation sites, the vast majority of known sites is not assigned to any protein kinase. Assigning changes in the phosphoproteome to the activity of individual kinases therefore remains a key challenge. A recent large-scale study systematically identified in vitro substrates for most human protein kinases. Here, we reprocessed and filtered these data to generate an in vitro Kinase-to-Phosphosite database (iKiP-DB). We show that iKiP-DB can accurately predict changes in kinase activity in published phosphoproteomic data sets for both well-studied and poorly characterized kinases. We apply iKiP-DB to a newly generated phosphoproteomic analysis of SARS-CoV-2 infected human lung epithelial cells and provide evidence for coronavirus-induced changes in host cell kinase activity. In summary, we show that iKiP-DB is widely applicable to facilitate the functional analysis of phosphoproteomic data sets.

Cellular glutamine synthesis is thought to be an important resistance factor in protecting cells from nutrient deprivation and may also contribute to drug resistance. The application of ‟targeted stable isotope resolved metabolomics” allowed to directly measure the activity of glutamine synthetase in the cell. With the help of this method, the fate of glutamine derived nitrogen within the biochemical network of the cells was traced. The application of stable isotope labelled substrates and analyses of isotope enrichment in metabolic intermediates allows the determination of metabolic activity and flux in biological systems. In our study we used stable isotope labelled substrates of glutamine synthetase to demonstrate its role in the starvation response of cancer cells. We applied 13C labelled glutamate and 15N labelled ammonium and determined the enrichment of both isotopes in glutamine and nucleotide species. Our results show that the metabolic compensatory pathways to overcome glutamine depletion depend on the ability to synthesise glutamine via glutamine synthetase. We demonstrate that the application of dual-isotope tracing can be used to address specific reactions within the biochemical network directly. Our study highlights the potential of concurrent isotope tracing methods in medical research.

The rapid rise of monkeypox (MPX) cases outside previously endemic areas prompts for a better understanding of the disease. We studied the plasma proteome of a group of MPX patients with a similar infection history and clinical manifestation typical for the current outbreak. We report that MPX in this case series is associated with a strong plasma proteomic response among nutritional and acute phase response proteins. Moreover, we report a correlation between plasma proteins and disease severity. Contrasting the MPX host response with that of COVID-19, we find a range of similarities, but also important differences. For instance, CFHR1 is induced in COVID-19, but suppressed in MPX, reflecting the different roles of the complement system in the two infectious diseases. Of note, the spatial overlap in response proteins suggested that a COVID-19 biomarker panel assay could be repurposed for MPX. Applying a targeted protein panel assay provided encouraging results and distinguished MPX cases from healthy controls. Hence, our results provide a first proteomic characterization of the MPX human host response and encourage further research on protein-panel assays in emerging infectious diseases.

The possibility to record proteomes in high throughput and at high quality has opened new avenues for biomedical research, drug discovery, systems biology, and clinical translation. However, high-throughput proteomic experiments often require high sample amounts and can be less sensitive compared to conventional proteomic experiments. Here, we introduce and benchmark Zeno SWATH MS, a data-independent acquisition technique that employs a linear ion trap pulsing (Zeno trap pulsing) to increase the sensitivity in high-throughput proteomic experiments. We demonstrate that when combined with fast micro- or analytical flow-rate chromatography, Zeno SWATH MS increases protein identification with low sample amounts. For instance, using 20 min micro-flow-rate chromatography, Zeno SWATH MS identified more than 5000 proteins consistently, and with a coefficient of variation of 6%, from a 62.5 ng load of human cell line tryptic digest. Using 5 min analytical flow-rate chromatography (800 µl/min), Zeno SWATH MS identified 4907 proteins from a triplicate injection of 2 µg of a human cell lysate, or more than 3000 proteins from a 250 ng tryptic digest. Zeno SWATH MS hence facilitates sensitive high-throughput proteomic experiments with low sample amounts, mitigating the current bottlenecks of high-throughput proteomics.

Background
Global healthcare systems continue to be challenged by the COVID-19 pandemic, and there is a need for clinical assays that can help optimise resource allocation, support treatment decisions, and accelerate the development and evaluation of new therapies.
Methods
We developed a multiplexed proteomics assay for determining disease severity and prognosis in COVID-19. The assay quantifies up to 50 peptides, derived from 30 known and newly introduced COVID-19-related protein markers, in a single measurement using routine-lab compatible analytical flow rate liquid chromatography and multiple reaction monitoring (LC-MRM). We conducted two observational studies in patients with COVID-19 hospitalised at Charité – Universitätsmedizin Berlin, Germany before (from March 1 to 26, 2020, n=30) and after (from April 4 to November 19, 2020, n=164) dexamethasone became standard of care. The study is registered in the German and the WHO International Clinical Trials Registry (DRKS00021688).
Findings
The assay produces reproducible (median inter-batch CV of 10.9%) absolute quantification of 47 peptides with high sensitivity (median LLOQ of 143 ng/ml) and accuracy (median 96.8%). In both studies, the assay reproducibly captured hallmarks of COVID-19 infection and severity, as it distinguished healthy individuals, mild, moderate, and severe COVID-19. In the post-dexamethasone cohort, the assay predicted survival with an accuracy of 0.83 (108/130), and death with an accuracy of 0.76 (26/34) in the median 2.5 weeks before the outcome, thereby outperforming compound clinical risk assessments such as SOFA, APACHE II, and ABCS scores.
Interpretation
Disease severity and clinical outcomes of patients with COVID-19 can be stratified and predicted by the routine-applicable panel assay that combines known and novel COVID-19 biomarkers. The prognostic value of this assay should be prospectively assessed in larger patient cohorts for future support of clinical decisions, including evaluation of sample flow in routine setting. The possibility to objectively classify COVID-19 severity can be helpful for monitoring of novel therapies, especially in early clinical trials.

Understanding the pathological properties of dysregulated protein networks in individual patients’ tumors is the basis for precision therapy. Functional experiments are commonly used, but cover only parts of the oncogenic signaling networks, whereas methods that reconstruct networks from omics data usually only predict average network features across tumors. Here, we show that the explainable AI method layer-wise relevance propagation (LRP) can infer protein interaction networks for individual patients from proteomic profiling data. LRP reconstructs average and individual interaction networks with an AUC of 0.99 and 0.93, respectively, and outperforms state-of-the-art network prediction methods for individual tumors. Using data from The Cancer Proteome Atlas, we identify known and potentially novel oncogenic network features, among which some are cancer-type specific and show only minor variation among patients, while others are present across certain tumor types but differ among individual patients. Our approach may therefore support predictive diagnostics in precision oncology by inferring “patient-level” oncogenic mechanisms.

The diagnosis of sinonasal tumors is challenging due to a heterogeneous spectrum of various differential diagnoses as well as poorly defined, disputed entities such as sinonasal undifferentiated carcinomas (SNUCs). In this study, we apply a machine learning algorithm based on DNA methylation patterns to classify sinonasal tumors with clinical-grade reliability. We further show that sinonasal tumors with SNUC morphology are not as undifferentiated as their current terminology suggests but rather reassigned to four distinct molecular classes defined by epigenetic, mutational and proteomic profiles. This includes two classes with neuroendocrine differentiation, characterized by IDH2 or SMARCA4/ARID1A mutations with an overall favorable clinical course, one class composed of highly aggressive SMARCB1-deficient carcinomas and another class with tumors that represent potentially previously misclassified adenoid cystic carcinomas. Our findings can aid in improving the diagnostic classification of sinonasal tumors and could help to change the current perception of SNUCs.

Severe COVID-19 is linked to both dysfunctional immune response and unrestrained immunopathology, and it remains unclear whether T cells contribute to disease pathology. Here, we combined single-cell transcriptomics and single-cell proteomics with mechanistic studies to assess pathogenic T cell functions and inducing signals. We identified highly activated CD16+ T cells with increased cytotoxic functions in severe COVID-19. CD16 expression enabled immune-complex-mediated, T cell receptor-independent degranulation and cytotoxicity not found in other diseases. CD16+ T cells from COVID-19 patients promoted microvascular endothelial cell injury and release of neutrophil and monocyte chemoattractants. CD16+ T cell clones persisted beyond acute disease maintaining their cytotoxic phenotype. Increased generation of C3a in severe COVID-19 induced activated CD16+ cytotoxic T cells. Proportions of activated CD16+ T cells and plasma levels of complement proteins upstream of C3a were associated with fatal outcome of COVID-19, supporting a pathological role of exacerbated cytotoxicity and complement activation in COVID-19.

The dia-PASEF technology uses ion mobility separation to reduce signal interferences and increase sensitivity in proteomic experiments. Here we present a two-dimensional peak-picking algorithm and generation of optimized spectral libraries, as well as take advantage of neural network-based processing of dia-PASEF data. Our computational platform boosts proteomic depth by up to 83% compared to previous work, and is specifically beneficial for fast proteomic experiments and those with low sample amounts. It quantifies over 5300 proteins in single injections recorded at 200 samples per day throughput using Evosep One chromatography system on a timsTOF Pro mass spectrometer and almost 9000 proteins in single injections recorded with a 93-min nanoflow gradient on timsTOF Pro 2, from 200 ng of HeLa peptides. A user-friendly implementation is provided through the incorporation of the algorithms in the DIA-NN software and by the FragPipe workflow for spectral library generation.

The immunomodulatory drugs (IMiDs) lenalidomide and pomalidomide are highly effective treatments for multiple myeloma. However, virtually all patients eventually relapse due to acquired drug resistance with resistance-causing genetic alterations being found only in a small subset of cases. To identify non-genetic mechanisms of drug resistance, we here perform integrated global quantitative tandem mass tag (TMT)-based proteomic and phosphoproteomic analyses and RNA sequencing in five paired pre-treatment and relapse samples from multiple myeloma patients. These analyses reveal a CDK6-governed protein resistance signature that includes myeloma high-risk factors such as TRIP13 and RRM1. Overexpression of CDK6 in multiple myeloma cell lines reduces sensitivity to IMiDs while CDK6 inhibition by palbociclib or CDK6 degradation by proteolysis targeting chimeras (PROTACs) is highly synergistic with IMiDs in vitro and in vivo. This work identifies CDK6 upregulation as a druggable target in IMiD-resistant multiple myeloma and highlights the use of proteomic studies to uncover non-genetic resistance mechanisms in cancer.

Global healthcare systems are challenged by the COVID-19 pandemic. There is a need to optimize allocation of treatment and resources in intensive care, as clinically established risk assessments such as SOFA and APACHE II scores show only limited performance for predicting the survival of severely ill COVID-19 patients. Additional tools are also needed to monitor treatment, including experimental therapies in clinical trials. Comprehensively capturing human physiology, we speculated that proteomics in combination with new data-driven analysis strategies could produce a new generation of prognostic discriminators. We studied two independent cohorts of patients with severe COVID-19 who required intensive care and invasive mechanical ventilation. SOFA score, Charlson comorbidity index, and APACHE II score showed limited performance in predicting the COVID-19 outcome. Instead, the quantification of 321 plasma protein groups at 349 timepoints in 50 critically ill patients receiving invasive mechanical ventilation revealed 14 proteins that showed trajectories different between survivors and non-survivors. A predictor trained on proteomic measurements obtained at the first time point at maximum treatment level (i.e. WHO grade 7), which was weeks before the outcome, achieved accurate classification of survivors (AUROC 0.81). We tested the established predictor on an independent validation cohort (AUROC 1.0). The majority of proteins with high relevance in the prediction model belong to the coagulation system and complement cascade. Our study demonstrates that plasma proteomics can give rise to prognostic predictors substantially outperforming current prognostic markers in intensive care.

Severe COVID-19 is linked to both dysfunctional immune response and unrestrained immunopathology, and it remains unclear whether T cells contribute to disease pathology. Here, we combined single-cell transcriptomics and single-cell proteomics with mechanistic studies to assess pathogenic T cell functions and inducing signals. We identified highly activated CD16+ T cells with increased cytotoxic functions in severe COVID-19. CD16 expression enabled immune-complex-mediated, T cell receptor-independent degranulation and cytotoxicity not found in other diseases. CD16+ T cells from COVID-19 patients promoted microvascular endothelial cell injury and release of neutrophil and monocyte chemoattractants. CD16+ T cell clones persisted beyond acute disease maintaining their cytotoxic phenotype. Increased generation of C3a in severe COVID-19 induced activated CD16+ cytotoxic T cells. Proportions of activated CD16+ T cells and plasma levels of complement proteins upstream of C3a were associated with fatal outcome of COVID-19, supporting a pathological role of exacerbated cytotoxicity and complement activation in COVID-19.

The particularly high mortality of epithelial ovarian cancer (EOC) is in part linked to limited understanding of its molecular signatures. Although there are data available on in situ N-glycosylation in EOC tissue, previous studies focused primarily on neutral N-glycan species and, hence, still little is known regarding EOC tissue-specific sialylation. In this proof-of-concept study, we implemented MALDI mass spectrometry imaging (MALDI-MSI) in combination with sialic acid derivatization to simultaneously investigate neutral and sialylated N-glycans in formalin-fixed paraffin-embedded tissue microarray specimens of less common EOC histotypes and non-malignant borderline ovarian tumor (BOT). The applied protocol allowed detecting over 50 m/z species, many of which showed differential tissue distribution. Most importantly, it could be demonstrated that α2,6- and α2,3-sialylated N-glycans are enriched in tissue regions corresponding to tumor and adjacent tumor-stroma, respectively. Interestingly, analogous N-glycosylation patterns were observed in tissue cores of BOT, suggesting that regio-specific N-glycan distribution might occur already in non-malignant ovarian pathologies. All in all, our data provide proof that the combination of MALDI-MSI and sialic acid derivatization is suitable for delineating regio-specific N-glycan distribution in EOC and BOT tissues and might serve as a promising strategy for future glycosylation-based biomarker discovery studies.

Available automated methods for peak detection in untargeted metabolomics suffer from poor precision. We present NeatMS, which uses machine learning based on a convoluted neural network to reduce the number and fraction of false peaks. NeatMS comes with a pre-trained model representing expert knowledge in the differentiation of true chemical signal from noise. Furthermore, it provides all necessary functions to easily train new models or improve existing ones by transfer learning. Thus, the tool improves peak curation and contributes to the robust and scalable analysis of large-scale experiments. We show how to integrate it into different liquid chromatography–mass spectrometry (LC-MS) analysis workflows, quantify its performance, and compare it to various other approaches. NeatMS software is available as open source on github under permissive MIT license and is also provided as easy-to-install PyPi and Bioconda packages.

The coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has been affecting the world since January 2020 and has caused millions of deaths. To gain a better insight into molecular changes underlying the COVID-19 disease, we investigated here the N-glycosylation of three immunoglobulin G (IgG) fractions isolated from plasma of 35 severe COVID-19 patients, namely total IgG1, total IgG2, and anti-Spike IgG, by means of MALDI-TOF-MS. All analyses were performed at the glycopeptide level to assure subclass- and site-specific information. For each COVID-19 patient, the analyses included three blood withdrawals at different time-points of hospitalization, which allowed profiling longitudinal alterations in IgG glycosylation. The COVID-19 patients presented altered IgG N-glycosylation profiles in all investigated IgG fractions. The most pronounced COVID-19-related changes were observed in the glycosylation profiles of antigen-specific anti-Spike IgG1. Anti-Spike IgG1 fucosylation and galactosylation showed the strongest variation during the disease course, with the difference in anti-Spike IgG1 fucosylation being significantly correlated with patients’ age. Decreases in anti-Spike IgG1 galactosylation and sialylation in the course of the disease were found to be significantly correlated with the difference in anti-Spike IgG plasma concentration. The present findings suggest that patients’ age and anti-S IgG abundance might influence IgG N-glycosylation alterations occurring in COVID-19.

Treating refractory or relapsed neuroblastoma remains challenging. Monitoring body fluids for tumor-derived molecular information indicating minimal residual disease supports more frequent diagnostic surveillance and may have the power to detect resistant subclones before they give rise to relapses. If actionable targets are identified from liquid biopsies, targeted treatment options can be considered earlier.
Total cfDNA concentrations in blood plasma from patients with high-risk neuroblastoma were higher than in healthy controls and consistently correlated with neuron-specific enolase levels and lactate dehydrogenase activity but not with 123I-meta-iodobenzylguanidine scores at relapse diagnosis. Targeted cfDNA diagnostics proved superior for early relapse detection to all current diagnostics in 2 patients. Marker analysis in cfDNA indicated intratumor heterogeneity for cell clones harboring MYCN amplifications and druggable ALK alterations that were not detectable in matched tumor tissue samples in 17 patients from our cohort. Proof of concept is provided for molecular target detection in cerebrospinal fluid from patients with isolated central nervous system relapses.
Tumor-specific alterations can be identified and monitored during disease course in liquid biopsies from pediatric patients with high-risk neuroblastoma. This approach to cfDNA surveillance warrants further prospective validation and exploitation for diagnostic purposes and to guide therapeutic decisions.

COVID-19 is highly variable in its clinical presentation, ranging from asymptomatic infection to severe organ damage and death. We characterized the time-dependent progression of the disease in 139 COVID-19 inpatients by measuring 86 accredited diagnostic parameters, such as blood cell counts and enzyme activities, as well as untargeted plasma proteomes at 687 sampling points. We report an initial spike in a systemic inflammatory response, which is gradually alleviated and followed by a protein signature indicative of tissue repair, metabolic reconstitution, and immunomodulation. We identify prognostic marker signatures for devising risk-adapted treatment strategies and use machine learning to classify therapeutic needs. We show that the machine learning models based on the proteome are transferable to an independent cohort. Our study presents a map linking routinely used clinical diagnostic parameters to plasma proteomes and their dynamics in an infectious disease.

Neutrophils provide a critical line of defense in immune responses to various pathogens, but also inflict self-damage upon transition to a hyperactivated, procoagulant state. Recent work has highlighted proinflammatory neutrophil phenotypes contributing to lung injury and acute respiratory distress syndrome (ARDS) in patients suffering from COVID-19. Here, we utilize state-of-the art mass spectrometry-based proteomics, transcriptomic and correlative analyses as well as functional in vitro and in vivo studies to dissect how neutrophils contribute to the progression to severe COVID-19. We identify a reinforcing loop of both systemic and neutrophil intrinsic interleukin-8 (CXCL8/IL-8) dysregulation, which initiates and perpetuates neutrophil-driven immunopathology. This positive feedback loop of systemic and neutrophil autocrine IL-8 production leads to an activated, prothrombotic neutrophil phenotype characterized by degranulation and neutrophil extracellular trap (NET) formation. In severe COVID-19, neutrophils directly initiate the coagulation and complement cascade, highlighting a link to the immunothrombotic state observed in these patients. Targeting the IL-8-CXCR-1/-2 axis interferes with this vicious cycle and attenuates neutrophil activation, degranulation, NETosis, and IL-8 release. Finally, we show that blocking IL-8-like signaling reduces SARS-CoV-2 spike protein-induced, hACE2-dependent pulmonary microthrombosis in mice. In summary, our data provide comprehensive insights into the activation mechanisms of neutrophils in COVID-19 and uncover a self-sustaining neutrophil-IL-8-axis as promising therapeutic target in severe SARS-CoV-2 infection.

Protein-protein interactions (PPIs) mediated by intrinsically disordered regions (IDRs) are often based on short linear motifs (SLiM). SLiMs are implicated in signal transduction and gene regulation, yet remain technically laborious and notoriously challenging to study. Here, we present an optimized method for a PRotein Interaction Screen on a peptide MAtrix (PRISMA) in combination with quantitative mass spectrometry. The protocol was benchmarked with previously described SLiM based PPIs using peptides derived from EGFR, SOS1, GLUT1 and CEBPB and extended to map binding partners of kinase activation loops. The detailed protocol provides practical considerations for setting up a PRISMA screen and subsequently implementing PRISMA on a liquid handling robotic platform as a cost effective high-throughput method. Optimized PRISMA can be universally applied to systematically study SLiM based interactions and associated post translational modifications (PTMs) or mutations to advance our understanding of the largely uncharacterized interactomes of intrinsically disordered protein regions.

Formalin-fixed paraffin-embedded (FFPE) tissues are a valuable resource for retrospective clinical studies. Here, we evaluate the feasibility of (phospho-)proteomics on FFPE lung tissue regarding protein extraction, quantification, pre-analytics, and sample size. After comparing protein extraction protocols, we use the best-performing protocol for the acquisition of deep (phospho-)proteomes from lung squamous cell and adenocarcinoma with >8,000 quantified proteins and >14,000 phosphosites with a tandem mass tag (TMT) approach. With a microscaled approach, we quantify 7,000 phosphosites, enabling the analysis of FFPE biopsies with limited tissue amounts. We also investigate the influence of pre-analytical variables including fixation time and heat-assisted de-crosslinking on protein extraction efficiency and proteome coverage. Our improved workflows provide quantitative information on protein abundance and phosphosite regulation for the most relevant oncogenes, tumor suppressors, and signaling pathways in lung cancer. Finally, we present general guidelines to which methods are best suited for different applications, highlighting TMT methods for comprehensive (phospho-)proteome profiling for focused clinical studies and label-free methods for large cohorts.

Accurate quantification of the proteome remains challenging for large sample series and longitudinal experiments. We report a data-independent acquisition method, Scanning SWATH, that accelerates mass spectrometric (MS) duty cycles, yielding quantitative proteomes in combination with short gradients and high-flow (800 µl min-1) chromatography. Exploiting a continuous movement of the precursor isolation window to assign precursor masses to tandem mass spectrometry (MS/MS) fragment traces, Scanning SWATH increases precursor identifications by ~70% compared to conventional data-independent acquisition (DIA) methods on 0.5-5-min chromatographic gradients. We demonstrate the application of ultra-fast proteomics in drug mode-of-action screening and plasma proteomics. Scanning SWATH proteomes capture the mode of action of fungistatic azoles and statins. Moreover, we confirm 43 and identify 11 new plasma proteome biomarkers of COVID-19 severity, advancing patient classification and biomarker discovery. Thus, our results demonstrate a substantial acceleration and increased depth in fast proteomic experiments that facilitate proteomic drug screens and clinical studies.

Risk classification plays a crucial role in clinical management and therapy decisions in children with neuroblastoma. Risk assessment is currently based on patient criteria and molecular factors in single tumor biopsies at diagnosis. Growing evidence of extensive neuroblastoma intratumor heterogeneity drives the need for novel diagnostics to assess molecular profiles more comprehensively in spatial resolution to better predict risk for tumor progression and therapy resistance. We present a pilot study investigating the feasibility and potential of matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) to identify spatial peptide heterogeneity in neuroblastoma tissues of divergent current risk classification: high versus low/intermediate risk. Univariate (receiver operating characteristic analysis) and multivariate (segmentation, principal component analysis) statistical strategies identified spatially discriminative risk-associated MALDI-based peptide signatures. The AHNAK nucleoprotein and collapsin response mediator protein 1 (CRMP1) were identified as proteins associated with these peptide signatures, and their differential expression in the neuroblastomas of divergent risk was immunohistochemically validated. This proof-of-concept study demonstrates that MALDI-MSI combined with univariate and multivariate analysis strategies can identify spatially discriminative risk-associated peptide signatures in neuroblastoma tissues. These results suggest a promising new analytical strategy improving risk classification and providing new biological insights into neuroblastoma intratumor heterogeneity.

Despite the correlation of clinical outcome and molecular subtypes of high-grade serous ovarian cancer (HGSOC), contemporary gene expression signatures have not been implemented in clinical practice to stratify patients for targeted therapy. Hence, we aimed to examine the potential of unsupervised matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI-IMS) to stratify patients who might benefit from targeted therapeutic strategies. Molecular subtyping of paraffin-embedded tissue samples from 279 HGSOC patients was performed by NanoString analysis (ground truth labeling). Next, we applied MALDI-IMS paired with machine-learning algorithms to identify distinct mass profiles on the same paraffin-embedded tissue sections and distinguish HGSOC subtypes by proteomic signature. Finally, we devised a novel approach to annotate spectra of stromal origin. We elucidated a MALDI-derived proteomic signature (135 peptides) able to classify HGSOC subtypes. Random forest classifiers achieved an area under the curve (AUC) of 0.983. Furthermore, we demonstrated that the exclusion of stroma-associated spectra provides tangible improvements to classification quality (AUC = 0.988). Moreover, novel MALDI-based stroma annotation achieved near-perfect classifications (AUC = 0.999). Here, we present a concept integrating MALDI-IMS with machine-learning algorithms to classify patients according to distinct molecular subtypes of HGSOC. This has great potential to assign patients for personalized treatment.

Dietary high salt (HS) is a leading risk factor for mortality and morbidity. Serum sodium transiently increases postprandially but can also accumulate at sites of inflammation affecting differentiation and function of innate and adaptive immune cells. Here, we focus on how changes in extracellular sodium, mimicking alterations in the circulation and tissues, affect the early metabolic, transcriptional, and functional adaption of human and murine mononuclear phagocytes.
Extracellular sodium was taken up into the intracellular compartment, followed by the inhibition of mitochondrial respiration in murine and human macrophages. Mechanistically, HS reduces mitochondrial membrane potential, electron transport chain complex II activity, oxygen consumption, and ATP production independently of the polarization status of macrophages. Subsequently, cell activation is altered with improved bactericidal function in HS-treated M1-like macrophages and diminished CD4+ T cell migration in HS-treated M2-like macrophages. Pharmacological uncoupling of the electron transport chain under normal salt phenocopies HS-induced transcriptional changes and bactericidal function of human and murine mononuclear phagocytes. Clinically, also in vivo, rise in plasma sodium concentration within the physiological range reversibly reduces mitochondrial function in human monocytes. In both a 14-day and single meal HS challenge, healthy volunteers displayed a plasma sodium increase of and respectively, that correlated with decreased monocytic mitochondrial oxygen consumption.

COVID-19-induced “acute respiratory distress syndrome” (ARDS) is associated with prolonged respiratory failure and high mortality, but the mechanistic basis of lung injury remains incompletely understood. Here, we analyze pulmonary immune responses and lung pathology in two cohorts of patients with COVID-19 ARDS using functional single-cell genomics, immunohistology, and electron microscopy. We describe an accumulation of CD163-expressing monocyte-derived macrophages that acquired a profibrotic transcriptional phenotype during COVID-19 ARDS. Gene set enrichment and computational data integration revealed a significant similarity between COVID-19-associated macrophages and profibrotic macrophage populations identified in idiopathic pulmonary fibrosis. COVID-19 ARDS was associated with clinical, radiographic, histopathological, and ultrastructural hallmarks of pulmonary fibrosis. Exposure of human monocytes to SARS-CoV-2, but not influenza A virus or viral RNA analogs, was sufficient to induce a similar profibrotic phenotype in vitro. In conclusion, we demonstrate that SARS-CoV-2 triggers profibrotic macrophage responses and pronounced fibroproliferative ARDS.

Longitudinal analyses of the innate immune system, including the earliest time points, are essential to understand the immunopathogenesis and clinical course of coronavirus disease (COVID-19). Here, we performed a detailed characterization of natural killer (NK) cells in 205 patients (403 samples; days 2 to 41 after symptom onset) from four independent cohorts using single-cell transcriptomics and proteomics together with functional studies. We found elevated interferon (IFN)-α plasma levels in early severe COVD-19 alongside increased NK cell expression of IFN-stimulated genes (ISGs) and genes involved in IFN-α signaling, while upregulation of tumor necrosis factor (TNF)-induced genes was observed in moderate diseases. NK cells exert anti-SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) activity but are functionally impaired in severe COVID-19. Further, NK cell dysfunction may be relevant for the development of fibrotic lung disease in severe COVID-19, as NK cells exhibited impaired anti-fibrotic activity. Our study indicates preferential IFN-α and TNF responses in severe and moderate COVID-19, respectively, and associates a prolonged IFN-α-induced NK cell response with poorer disease outcome.

In COVID-19, immune responses are key in determining disease severity. However, cellular mechanisms at the onset of inflammatory lung injury in SARS-CoV-2 infection, particularly involving endothelial cells, remain ill-defined. Using Syrian hamsters as a model for moderate COVID-19, we conduct a detailed longitudinal analysis of systemic and pulmonary cellular responses, and corroborate it with datasets from COVID-19 patients. Monocyte-derived macrophages in lungs exert the earliest and strongest transcriptional response to infection, including induction of pro-inflammatory genes, while epithelial cells show weak alterations. Without evidence for productive infection, endothelial cells react, depending on cell subtypes, by strong and early expression of anti-viral, pro-inflammatory, and T cell recruiting genes. Recruitment of cytotoxic T cells as well as emergence of IgM antibodies precede viral clearance at day 5 post infection. Investigating SARS-CoV-2 infected Syrian hamsters thus identifies cell type-specific effector functions, providing detailed insights into pathomechanisms of COVID-19 and informing therapeutic strategies.

Despite the overall poor prognosis of pancreatic cancer there is heterogeneity in clinical courses of tumors not assessed by conventional risk stratification. This yields the need of additional markers for proper assessment of prognosis and multimodal clinical management. We provide a proof of concept study evaluating the feasibility of Matrix-assisted laser desorption/ionization (MALDI) mass spectrometry imaging (MSI) to identify specific peptide signatures linked to prognostic parameters of pancreatic cancer. On 18 patients with exocrine pancreatic cancer after tumor resection, MALDI imaging analysis was performed additional to histopathological assessment. Principal component analysis (PCA) was used to explore discrimination of peptide signatures of prognostic histopathological features and receiver operator characteristic (ROC) to identify which specific m/z values are the most discriminative between the prognostic subgroups of patients. Out of 557 aligned m/z values discriminate peptide signatures for the prognostic histopathological features lymphatic vessel invasion (pL, 16 m/z values, eight proteins), nodal metastasis (pN, two m/z values, one protein) and angioinvasion (pV, 4 m/z values, two proteins) were identified. These results yield proof of concept that MALDI-MSI of pancreatic cancer tissue is feasible to identify peptide signatures of prognostic relevance and can augment risk assessment.

Very high risk neuroblastoma is characterised by increased MAPK signalling, and targeting MAPK signalling is a promising therapeutic strategy. We used a deeply characterised panel of neuroblastoma cell lines and found that the sensitivity to MEK inhibitors varied drastically between these cell lines. By generating quantitative perturbation data and mathematical modelling, we determined potential resistance mechanisms. We found that negative feedbacks within MAPK signalling and via the IGF receptor mediate re-activation of MAPK signalling upon treatment in resistant cell lines. By using cell-line specific models, we predict that combinations of MEK inhibitors with RAF or IGFR inhibitors can overcome resistance, and tested these predictions experimentally. In addition, phospho-proteomic profiling confirmed the cell-specific feedback effects and synergy of MEK and IGFR targeted treatment. Our study shows that a quantitative understanding of signalling and feedback mechanisms facilitated by models can help to develop and optimise therapeutic strategies. Our findings should be considered for the planning of future clinical trials introducing MEKi in the treatment of neuroblastoma.

Recent advances in cancer research and diagnostics largely rely on new developments in microscopic or molecular profiling techniques, offering high levels of detail with respect to either spatial or molecular features, but usually not both. Here, we present an explainable machine-learning approach for the integrated profiling of morphological, molecular and clinical features from breast cancer histology. First, our approach allows for the robust detection of cancer cells and tumour-infiltrating lymphocytes in histological images, providing precise heatmap visualizations explaining the classifier decisions. Second, molecular features, including DNA methylation, gene expression, copy number variations, somatic mutations and proteins are predicted from histology. Molecular predictions reach balanced accuracies up to 78%, whereas accuracies of over 95% can be achieved for subgroups of patients. Finally, our explainable AI approach allows assessment of the link between morphological and molecular cancer properties. The resulting computational multiplex-histology analysis can help promote basic cancer research and precision medicine through an integrated diagnostic scoring of histological, clinical and molecular features.

The COVID-19 pandemic is an unprecedented global challenge, and point-of-care diagnostic classifiers are urgently required. Here, we present a platform for ultra-high-throughput serum and plasma proteomics that builds on ISO13485 standardization to facilitate simple implementation in regulated clinical laboratories. Our low-cost workflow handles up to 180 samples per day, enables high precision quantification, and reduces batch effects for large-scale and longitudinal studies. We use our platform on samples collected from a cohort of early hospitalized cases of the SARS-CoV-2 pandemic and identify 27 potential biomarkers that are differentially expressed depending on the WHO severity grade of COVID-19. They include complement factors, the coagulation system, inflammation modulators, and pro-inflammatory factors upstream and downstream of interleukin 6. All protocols and software for implementing our approach are freely available. In total, this work supports the development of routine proteomic assays to aid clinical decision making and generate hypotheses about potential COVID-19 therapeutic targets.

We present an easy-to-use integrated software suite, DIA-NN, that exploits deep neural networks and new quantification and signal correction strategies for the processing of data-independent acquisition (DIA) proteomics experiments. DIA-NN improves the identification and quantification performance in conventional DIA proteomic applications, and is particularly beneficial for high-throughput applications, as it is fast and enables deep and confident proteome coverage when used in combination with fast chromatographic methods.

Modern life science is generating large data sets at a unprecedented speed. A major source of data are so-called omics (e.g., genomics, metabolomics, or proteomics) experiments. Consequently, management and analysis of scientific data has become a major challenge. Further, the heterogeneity of projects makes “one size fits all” data analysis systems infeasible and calls for specialized data analysis platforms. The authors are actively developing applications for the FAIR (findable, accessible, interoperable, and reuseable, cf. Wilkinson et al. (2016) data management of omics data and their analysis. In order to prevent duplication of work, we have extracted the commonly useful components into SODAR Core, a Python framework on top of Django. It is used in our actively developed applications and also proved useful in internal web app prototypes. Examples for using SODAR Core in the development of scientific data management web apps is Digestiflow (Holtgrewe et al., 2019) and the Filesfolders module shipping with SODAR Core. An example for using SODAR Core in the development of scientific data analysis web apps is VarFish (Holtgrewe et al., 2020).