Integrating single-cell multi-omics and prior biological knowledge for a functional characterization of the immune system.

The immune system comprises diverse specialized cell types that cooperate to defend the host against a wide range of pathogenic threats. Recent advancements in single-cell and spatial multi-omics technologies provide rich information about the molecular state of immune cells. Here, we review how the integration of single-cell and spatial multi-omics data with prior knowledge-gathered from decades of detailed biochemical studies-allows us to obtain functional insights, focusing on gene regulatory processes and cell-cell interactions. We present diverse applications in immunology and critically assess underlying assumptions and limitations. Finally, we offer a perspective on the ongoing technological and algorithmic developments that promise to get us closer to a systemic mechanistic understanding of the immune system.

Dynamic 3D proteomes reveal protein functional alterations at high resolution in situ.

Biological processes are regulated by intermolecular interactions and chemical modifications that do not affect protein levels, thus escaping detection in classical proteomic screens. We demonstrate here that a global protein structural readout based on limited proteolysis-mass spectrometry (LiP-MS) detects many such functional alterations, simultaneously and in situ, in bacteria undergoing nutrient adaptation and in yeast responding to acute stress. The structural readout, visualized as structural barcodes, captured enzyme activity changes, phosphorylation, protein aggregation, and complex formation, with the resolution of individual regulated functional sites such as binding and active sites. Comparison with prior knowledge, including other 'omics data, showed that LiP-MS detects many known functional alterations within well-studied pathways. It suggested distinct metabolite-protein interactions and enabled identification of a fructose-1,6-bisphosphate-based regulatory mechanism of glucose uptake in E. coli. The structural readout dramatically increases classical proteomics coverage, generates mechanistic hypotheses, and paves the way for in situ structural systems biology.

Dynamic chromatin accessibility licenses STAT5- and STAT6-dependent innate-like function of TH9 cells to promote allergic inflammation.

Allergic diseases are a major global health issue. Interleukin (IL)-9-producing helper T (TH9) cells promote allergic inflammation, yet TH9 cell effector functions are incompletely understood because their lineage instability makes them challenging to study. Here we found that resting TH9 cells produced IL-9 independently of T cell receptor (TCR) restimulation, due to STAT5- and STAT6-dependent bystander activation. This mechanism was seen in circulating cells from allergic patients and was restricted to recently activated cells. STAT5-dependent Il9/IL9 regulatory elements underwent remodeling over time, inactivating the locus. A broader 'allergic TH9' transcriptomic and epigenomic program was also unstable. In vivo, TH9 cells induced airway inflammation via TCR-independent, STAT-dependent mechanisms. In allergic patients, TH9 cell expansion was associated with responsiveness to JAK inhibitors. These findings suggest that TH9 cell instability is a negative checkpoint on bystander activation that breaks down in allergy and that JAK inhibitors should be considered for allergic patients with TH9 cell expansion.

B cell expansion hinders the stroma-epithelium regenerative cross talk during mucosal healing.

Therapeutic promotion of intestinal regeneration holds great promise, but defining the cellular mechanisms that influence tissue regeneration remains an unmet challenge. To gain insight into the process of mucosal healing, we longitudinally examined the immune cell composition during intestinal damage and regeneration. B cells were the dominant cell type in the healing colon, and single-cell RNA sequencing (scRNA-seq) revealed expansion of an IFN-induced B cell subset during experimental mucosal healing that predominantly located in damaged areas and associated with colitis severity. B cell depletion accelerated recovery upon injury, decreased epithelial ulceration, and enhanced gene expression programs associated with tissue remodeling. scRNA-seq from the epithelial and stromal compartments combined with spatial transcriptomics and multiplex immunostaining showed that B cells decreased interactions between stromal and epithelial cells during mucosal healing. Activated B cells disrupted the epithelial-stromal cross talk required for organoid survival. Thus, B cell expansion during injury impairs epithelial-stromal cell interactions required for mucosal healing, with implications for the treatment of IBD.

A Landscape of Pharmacogenomic Interactions in Cancer.

Systematic studies of cancer genomes have provided unprecedented insights into the molecular nature of cancer. Using this information to guide the development and application of therapies in the clinic is challenging. Here, we report how cancer-driven alterations identified in 11,289 tumors from 29 tissues (integrating somatic mutations, copy number alterations, DNA methylation, and gene expression) can be mapped onto 1,001 molecularly annotated human cancer cell lines and correlated with sensitivity to 265 drugs. We find that cell lines faithfully recapitulate oncogenic alterations identified in tumors, find that many of these associate with drug sensitivity/resistance, and highlight the importance of tissue lineage in mediating drug response. Logic-based modeling uncovers combinations of alterations that sensitize to drugs, while machine learning demonstrates the relative importance of different data types in predicting drug response. Our analysis and datasets are rich resources to link genotypes with cellular phenotypes and to identify therapeutic options for selected cancer sub-populations.

Hyperactivated PI3Kδ promotes self and commensal reactivity at the expense of optimal humoral immunity.

Gain-of-function mutations in the gene encoding the phosphatidylinositol-3-OH kinase catalytic subunit p110δ (PI3Kδ) result in a human primary immunodeficiency characterized by lymphoproliferation, respiratory infections and inefficient responses to vaccines. However, what promotes these immunological disturbances at the cellular and molecular level remains unknown. We generated a mouse model that recapitulated major features of this disease and used this model and patient samples to probe how hyperactive PI3Kδ fosters aberrant humoral immunity. We found that mutant PI3Kδ led to co-stimulatory receptor ICOS-independent increases in the abundance of follicular helper T cells (TFH cells) and germinal-center (GC) B cells, disorganized GCs and poor class-switched antigen-specific responses to immunization, associated with altered regulation of the transcription factor FOXO1 and pro-apoptotic and anti-apoptotic members of the BCL-2 family. Notably, aberrant responses were accompanied by increased reactivity to gut bacteria and a broad increase in autoantibodies that were dependent on stimulation by commensal microbes. Our findings suggest that proper regulation of PI3Kδ is critical for ensuring optimal host-protective humoral immunity despite tonic stimulation from the commensal microbiome.

Deep spatial profiling of human COVID-19 brains reveals neuroinflammation with distinct microanatomical microglia-T-cell interactions.

COVID-19 can cause severe neurological symptoms, but the underlying pathophysiological mechanisms are unclear. Here, we interrogated the brain stems and olfactory bulbs in postmortem patients who had COVID-19 using imaging mass cytometry to understand the local immune response at a spatially resolved, high-dimensional, single-cell level and compared their immune map to non-COVID respiratory failure, multiple sclerosis, and control patients. We observed substantial immune activation in the central nervous system with pronounced neuropathology (astrocytosis, axonal damage, and blood-brain-barrier leakage) and detected viral antigen in ACE2-receptor-positive cells enriched in the vascular compartment. Microglial nodules and the perivascular compartment represented COVID-19-specific, microanatomic-immune niches with context-specific cellular interactions enriched for activated CD8+ T cells. Altered brain T-cell-microglial interactions were linked to clinical measures of systemic inflammation and disturbed hemostasis. This study identifies profound neuroinflammation with activation of innate and adaptive immune cells as correlates of COVID-19 neuropathology, with implications for potential therapeutic strategies.

Prospective derivation of a living organoid biobank of colorectal cancer patients.

In Rspondin-based 3D cultures, Lgr5 stem cells from multiple organs form ever-expanding epithelial organoids that retain their tissue identity. We report the establishment of tumor organoid cultures from 20 consecutive colorectal carcinoma (CRC) patients. For most, organoids were also generated from adjacent normal tissue. Organoids closely recapitulate several properties of the original tumor. The spectrum of genetic changes within the "living biobank" agrees well with previous large-scale mutational analyses of CRC. Gene expression analysis indicates that the major CRC molecular subtypes are represented. Tumor organoids are amenable to high-throughput drug screens allowing detection of gene-drug associations. As an example, a single organoid culture was exquisitely sensitive to Wnt secretion (porcupine) inhibitors and carried a mutation in the negative Wnt feedback regulator RNF43, rather than in APC. Organoid technology may fill the gap between cancer genetics and patient trials, complement cell-line- and xenograft-based drug studies, and allow personalized therapy design. PAPERCLIP.

Gene regulatory network inference in the era of single-cell multi-omics.

The interplay between chromatin, transcription factors and genes generates complex regulatory circuits that can be represented as gene regulatory networks (GRNs). The study of GRNs is useful to understand how cellular identity is established, maintained and disrupted in disease. GRNs can be inferred from experimental data - historically, bulk omics data - and/or from the literature. The advent of single-cell multi-omics technologies has led to the development of novel computational methods that leverage genomic, transcriptomic and chromatin accessibility information to infer GRNs at an unprecedented resolution. Here, we review the key principles of inferring GRNs that encompass transcription factor-gene interactions from transcriptomics and chromatin accessibility data. We focus on the comparison and classification of methods that use single-cell multimodal data. We highlight challenges in GRN inference, in particular with respect to benchmarking, and potential further developments using additional data modalities.

Spatial multi-omic map of human myocardial infarction.

Myocardial infarction is a leading cause of death worldwide1. Although advances have been made in acute treatment, an incomplete understanding of remodelling processes has limited the effectiveness of therapies to reduce late-stage mortality2. Here we generate an integrative high-resolution map of human cardiac remodelling after myocardial infarction using single-cell gene expression, chromatin accessibility and spatial transcriptomic profiling of multiple physiological zones at distinct time points in myocardium from patients with myocardial infarction and controls. Multi-modal data integration enabled us to evaluate cardiac cell-type compositions at increased resolution, yielding insights into changes of the cardiac transcriptome and epigenome through the identification of distinct tissue structures of injury, repair and remodelling. We identified and validated disease-specific cardiac cell states of major cell types and analysed them in their spatial context, evaluating their dependency on other cell types. Our data elucidate the molecular principles of human myocardial tissue organization, recapitulating a gradual cardiomyocyte and myeloid continuum following ischaemic injury. In sum, our study provides an integrative molecular map of human myocardial infarction, represents an essential reference for the field and paves the way for advanced mechanistic and therapeutic studies of cardiac disease.

Decoding myofibroblast origins in human kidney fibrosis.

Kidney fibrosis is the hallmark of chronic kidney disease progression; however, at present no antifibrotic therapies exist1-3. The origin, functional heterogeneity and regulation of scar-forming cells that occur during human kidney fibrosis remain poorly understood1,2,4. Here, using single-cell RNA sequencing, we profiled the transcriptomes of cells from the proximal and non-proximal tubules of healthy and fibrotic human kidneys to map the entire human kidney. This analysis enabled us to map all matrix-producing cells at high resolution, and to identify distinct subpopulations of pericytes and fibroblasts as the main cellular sources of scar-forming myofibroblasts during human kidney fibrosis. We used genetic fate-tracing, time-course single-cell RNA sequencing and ATAC-seq (assay for transposase-accessible chromatin using sequencing) experiments in mice, and spatial transcriptomics in human kidney fibrosis, to shed light on the cellular origins and differentiation of human kidney myofibroblasts and their precursors at high resolution. Finally, we used this strategy to detect potential therapeutic targets, and identified NKD2 as a myofibroblast-specific target in human kidney fibrosis.

Analysis of the Human Kinome and Phosphatome by Mass Cytometry Reveals Overexpression-Induced Effects on Cancer-Related Signaling.

Kinase and phosphatase overexpression drives tumorigenesis and drug resistance. We previously developed a mass-cytometry-based single-cell proteomics approach that enables quantitative assessment of overexpression effects on cell signaling. Here, we applied this approach in a human kinome- and phosphatome-wide study to assess how 649 individually overexpressed proteins modulated cancer-related signaling in HEK293T cells in an abundance-dependent manner. Based on these data, we expanded the functional classification of human kinases and phosphatases and showed that the overexpression effects include non-catalytic roles. We detected 208 previously unreported signaling relationships. The signaling dynamics analysis indicated that the overexpression of ERK-specific phosphatases sustains proliferative signaling. This suggests a phosphatase-driven mechanism of cancer progression. Moreover, our analysis revealed a drug-resistant mechanism through which overexpression of tyrosine kinases, including SRC, FES, YES1, and BLK, induced MEK-independent ERK activation in melanoma A375 cells. These proteins could predict drug sensitivity to BRAF-MEK concurrent inhibition in cells carrying BRAF mutations.

MetalinksDB: a flexible and contextualizable resource of metabolite-protein interactions.

From the catalytic breakdown of nutrients to signaling, interactions between metabolites and proteins play an essential role in cellular function. An important case is cell-cell communication, where metabolites, secreted into the microenvironment, initiate signaling cascades by binding to intra- or extracellular receptors of neighboring cells. Protein-protein cell-cell communication interactions are routinely predicted from transcriptomic data. However, inferring metabolite-mediated intercellular signaling remains challenging, partially due to the limited size of intercellular prior knowledge resources focused on metabolites. Here, we leverage knowledge-graph infrastructure to integrate generalistic metabolite-protein with curated metabolite-receptor resources to create MetalinksDB. MetalinksDB is an order of magnitude larger than existing metabolite-receptor resources and can be tailored to specific biological contexts, such as diseases, pathways, or tissue/cellular locations. We demonstrate MetalinksDB's utility in identifying deregulated processes in renal cancer using multi-omics bulk data. Furthermore, we infer metabolite-driven intercellular signaling in acute kidney injury using spatial transcriptomics data. MetalinksDB is a comprehensive and customizable database of intercellular metabolite-protein interactions, accessible via a web interface (https://metalinks.omnipathdb.org/) and programmatically as a knowledge graph (https://github.com/biocypher/metalinks). We anticipate that by enabling diverse analyses tailored to specific biological contexts, MetalinksDB will facilitate the discovery of disease-relevant metabolite-mediated intercellular signaling processes.

NADH Shuttling Couples Cytosolic Reductive Carboxylation of Glutamine with Glycolysis in Cells with Mitochondrial Dysfunction.

The bioenergetics and molecular determinants of the metabolic response to mitochondrial dysfunction are incompletely understood, in part due to a lack of appropriate isogenic cellular models of primary mitochondrial defects. Here, we capitalize on a recently developed cell model with defined levels of m.8993T>G mutation heteroplasmy, mTUNE, to investigate the metabolic underpinnings of mitochondrial dysfunction. We found that impaired utilization of reduced nicotinamide adenine dinucleotide (NADH) by the mitochondrial respiratory chain leads to cytosolic reductive carboxylation of glutamine as a new mechanism for cytosol-confined NADH recycling supported by malate dehydrogenase 1 (MDH1). We also observed that increased glycolysis in cells with mitochondrial dysfunction is associated with increased cell migration in an MDH1-dependent fashion. Our results describe a novel link between glycolysis and mitochondrial dysfunction mediated by reductive carboxylation of glutamine.

Assessing the impact of transcriptomics data analysis pipelines on downstream functional enrichment results.

Transcriptomics is widely used to assess the state of biological systems. There are many tools for the different steps, such as normalization, differential expression, and enrichment. While numerous studies have examined the impact of method choices on differential expression results, little attention has been paid to their effects on further downstream functional analysis, which typically provides the basis for interpretation and follow-up experiments. To address this, we introduce FLOP, a comprehensive nextflow-based workflow combining methods to perform end-to-end analyses of transcriptomics data. We illustrate FLOP on datasets ranging from end-stage heart failure patients to cancer cell lines. We discovered effects not noticeable at the gene-level, and observed that not filtering the data had the highest impact on the correlation between pipelines in the gene set space. Moreover, we performed three benchmarks to evaluate the 12 pipelines included in FLOP, and confirmed that filtering is essential in scenarios of expected moderate-to-low biological signal. Overall, our results underscore the impact of carefully evaluating the consequences of the choice of preprocessing methods on downstream enrichment analyses. We envision FLOP as a valuable tool to measure the robustness of functional analyses, ultimately leading to more reliable and conclusive biological findings.

Drugst.One - a plug-and-play solution for online systems medicine and network-based drug repurposing.

In recent decades, the development of new drugs has become increasingly expensive and inefficient, and the molecular mechanisms of most pharmaceuticals remain poorly understood. In response, computational systems and network medicine tools have emerged to identify potential drug repurposing candidates. However, these tools often require complex installation and lack intuitive visual network mining capabilities. To tackle these challenges, we introduce Drugst.One, a platform that assists specialized computational medicine tools in becoming user-friendly, web-based utilities for drug repurposing. With just three lines of code, Drugst.One turns any systems biology software into an interactive web tool for modeling and analyzing complex protein-drug-disease networks. Demonstrating its broad adaptability, Drugst.One has been successfully integrated with 21 computational systems medicine tools. Available at https://drugst.one, Drugst.One has significant potential for streamlining the drug discovery process, allowing researchers to focus on essential aspects of pharmaceutical treatment research.

Metabolic Communication by SGLT2 Inhibition.

BACKGROUND: SGLT2 (sodium-glucose cotransporter 2) inhibitors (SGLT2i) can protect the kidneys and heart, but the underlying mechanism remains poorly understood. METHODS: To gain insights on primary effects of SGLT2i that are not confounded by pathophysiologic processes or are secondary to improvement by SGLT2i, we performed an in-depth proteomics, phosphoproteomics, and metabolomics analysis by integrating signatures from multiple metabolic organs and body fluids after 1 week of SGLT2i treatment of nondiabetic as well as diabetic mice with early and uncomplicated hyperglycemia. RESULTS: Kidneys of nondiabetic mice reacted most strongly to SGLT2i in terms of proteomic reconfiguration, including evidence for less early proximal tubule glucotoxicity and a broad downregulation of the apical uptake transport machinery (including sodium, glucose, urate, purine bases, and amino acids), supported by mouse and human SGLT2 interactome studies. SGLT2i affected heart and liver signaling, but more reactive organs included the white adipose tissue, showing more lipolysis, and, particularly, the gut microbiome, with a lower relative abundance of bacteria taxa capable of fermenting phenylalanine and tryptophan to cardiovascular uremic toxins, resulting in lower plasma levels of these compounds (including p-cresol sulfate). SGLT2i was detectable in murine stool samples and its addition to human stool microbiota fermentation recapitulated some murine microbiome findings, suggesting direct inhibition of fermentation of aromatic amino acids and tryptophan. In mice lacking SGLT2 and in patients with decompensated heart failure or diabetes, the SGLT2i likewise reduced circulating p-cresol sulfate, and p-cresol impaired contractility and rhythm in human induced pluripotent stem cell-derived engineered heart tissue. CONCLUSIONS: SGLT2i reduced microbiome formation of uremic toxins such as p-cresol sulfate and thereby their body exposure and need for renal detoxification, which, combined with direct kidney effects of SGLT2i, including less proximal tubule glucotoxicity and a broad downregulation of apical transporters (including sodium, amino acid, and urate uptake), provides a metabolic foundation for kidney and cardiovascular protection.

Complementing Cell Taxonomies with a Multicellular Analysis of Tissues.

The application of single-cell molecular profiling coupled with spatial technologies has enabled charting of cellular heterogeneity in reference tissues and in disease. This new wave of molecular data has highlighted the expected diversity of single-cell dynamics upon shared external queues and spatial organizations. However, little is known about the relationship between single-cell heterogeneity and the emergence and maintenance of robust multicellular processes in developed tissues and its role in (patho)physiology. Here, we present emerging computational modeling strategies that use increasingly available large-scale cross-condition single-cell and spatial datasets to study multicellular organization in tissues and complement cell taxonomies. This perspective should enable us to better understand how cells within tissues collectively process information and adapt synchronized responses in disease contexts and to bridge the gap between structural changes and functions in tissues.

Molecular causality in the advent of foundation models.

Correlation is not causation: this simple and uncontroversial statement has far-reaching implications. Defining and applying causality in biomedical research has posed significant challenges to the scientific community. In this perspective, we attempt to connect the partly disparate fields of systems biology, causal reasoning, and machine learning to inform future approaches in the field of systems biology and molecular medicine.