Project description:Kilian2024 - Immune cell dynamics in Cue-Induced Extended Human Colitis Model
Single-cell technologies such as scRNA-seq and flow cytometry provide critical insights into immune cell behavior in inflammatory bowel disease (IBD). However, integrating these datasets into computational models for dynamic analysis remains challenging. Here, Kilian et al., (2024) developed a deterministic ODE-based model that incorporates these technologies to study immune cell population changes in murine colitis. The model parameters were optimized to fit experimental data, ensuring an accurate representation of immune cell behavior over time. It was then validated by comparing simulations with experimental data using Pearson’s correlation and further tested on independent datasets to confirm its robustness. Additionally, the model was applied to clinical bulk RNA-seq data from human IBD patients, providing valuable insights into immune system dynamics and potential therapeutic strategies.
Figure 4c, obtained from the simulation of human colitis model is highlighted here.
This model is described in the article:
Kilian, C., Ulrich, H., Zouboulis, V.A. et al. Longitudinal single-cell data informs deterministic modelling of inflammatory bowel disease. npj Syst Biol Appl 10, 69 (2024). https://doi.org/10.1038/s41540-024-00395-9
Abstract:
Single-cell-based methods such as flow cytometry or single-cell mRNA sequencing (scRNA-seq) allow deep molecular and cellular profiling of immunological processes. Despite their high throughput, however, these measurements represent only a snapshot in time. Here, we explore how longitudinal single-cell-based datasets can be used for deterministic ordinary differential equation (ODE)-based modelling to mechanistically describe immune dynamics. We derived longitudinal changes in cell numbers of colonic cell types during inflammatory bowel disease (IBD) from flow cytometry and scRNA-seq data of murine colitis using ODE-based models. Our mathematical model generalised well across different protocols and experimental techniques, and we hypothesised that the estimated model parameters reflect biological processes. We validated this prediction of cellular turnover rates with KI-67 staining and with gene expression information from the scRNA-seq data not used for model fitting. Finally, we tested the translational relevance of the mathematical model by deconvolution of longitudinal bulk mRNA-sequencing data from a cohort of human IBD patients treated with olamkicept. We found that neutrophil depletion may contribute to IBD patients entering remission. The predictive power of IBD deterministic modelling highlights its potential to advance our understanding of immune dynamics in health and disease.
This model was curated during the Hackathon hosted by BioMed X GmbH in 2024.
Project description:Pancreatic cancer is a complex disease with a desmoplastic stroma, extreme hypoxia, and inherent resistance to therapy. Understanding the signaling and adaptive response of such an aggressive cancer is key to making advances in therapeutic efficacy and understanding disease progression. Redox factor-1 (Ref-1), a redox signaling protein, regulates the DNA binding activity of several transcription factors, including HIF-1. The conversion of HIF-1 from an oxidized to reduced state leads to enhancement of its DNA binding. In our previously published work, knockdown of Ref-1 under normoxia resulted in altered gene expression patterns on pathways including EIF2, protein kinase A, and mTOR. In this study, single cell RNA sequencing (scRNA-seq) and proteomics were used to explore the effects of Ref-1 on metabolic pathways under hypoxia.Results: We also integrated the scRNA data analysis with the proteomic analysis and found that the differentially expressed genes and pathways identified from the scRNA-seq data are highly consistent to the significant proteins observed in the proteomics data, especially for the upregulated cell cycle and transcription pathways and downregulated metabolic, apoptosis and signaling pathways under hypoxia. Conclusion: The scRNA-seq and proteomics data consistently demonstrated down-regulated central metabolism pathways in APE1/Ref-1 knockdown vs scrambled control under both normoxia and hypoxia conditions. Experimental Methods: scRNA-seq comparing pancreatic cancer cells expressing less than 20% of the Ref-1 protein was analyzed using left truncated mixture Gaussian model. Matched samples were also collected for bulk proteomic analysis of the four conditions. scRNA-seq data was validated using proteomics and qRT-PCR. Ref-1’s role in mitochondrial function was confirmed using mitochondrial function assays and qRT-PCR. Results: We also integrated the scRNA data analysis with the proteomic analysis and found that the differentially expressed genes and pathways identified from the scRNA-seq data are highly consistent to the significant proteins observed in the proteomics data, especially for the upregulated cell cycle and transcription pathways and downregulated metabolic, apoptosis and signaling pathways under hypoxia. Conclusion: The scRNA-seq and proteomics data consistently demonstrated down-regulated central metabolism pathways in APE1/Ref-1 knockdown vs scrambled control under both normoxia and hypoxia conditions.
Project description:Idiopathic pulmonary fibrosis (IPF) is a progressive lung disease characterized by repetitive alveolar injuries with excessive deposition of extracellular matrix (ECM) proteins. A crucial need in understanding IPF pathogenesis is identifying cell types associated with histopathological regions, particularly local fibrosis centers known as fibroblast foci. To address this, we integrated published spatial transcriptomics and single-cell RNA sequencing (scRNA-seq) transcriptomics and adopted the Query method and the Overlap method to determine cell type enrichments in histopathological regions. Distinct fibroblast cell types are highly associated with fibroblast foci, and transitional alveolar type 2 and aberrant KRT5-/KRT17+ epithelial cells are associated with morphologically normal alveoli in human IPF lungs. Furthermore, we employed laser capture microdissection directed mass spectrometry to profile proteins. By comparing with another published similar dataset, common differentially expressed proteins and enriched pathways related to ECM structure organization and collagen processing were identified in fibroblast foci. Importantly, cell type enrichment results from innovative spatial proteomics and scRNA-seq data integration accord with those from spatial transcriptomics and scRNA-seq data integration, supporting the capability and versatility of the entire approach. In summary, we integrated spatial multi-omics with scRNA-seq data to identify disease-associated cell types and potential targets for novel therapies in IPF intervention. The approach can be further applied to other disease areas characterized by spatial heterogeneity.
Project description:Energy metabolism and extracellular matrix function are closely connected to orchestrate and maintain tissue organization, but the crosstalk is poorly understood. Here, we used scRNA-seq analysis to uncover the importance of respiration for extracellular matrix homeostasis in mature cartilage. Genetic inhibition of respiration in cartilage results in the expansion of a central area of 1-month-old mouse femur head cartilage showing disorganized chondrocytes and increased deposition of extracellular matrix material. scRNA-seq analysis identified a cluster-specific decrease in mitochondrial DNA-encoded respiratory chain genes and a unique regulation of extracellular matrix-related genes in nonarticular chondrocyte clusters. These changes were associated with alterations in extracellular matrix composition, a shift in the collagen/non-collagen protein content and an increase of collagen crosslinking and ECM stiffness. The results demonstrate, based on findings of the scRNA-seq analysis, that respiration is a key factor contributing to ECM integrity and mechanostability in cartilage and presumably also in many other tissues.
Project description:This study delves into the proteomic intricacies of drug-resistant cells (DRCs) within prostate cancer, which are known for their pivotal roles in therapeutic resistance, relapse, and metastasis. Utilizing single-cell proteomics (SCP) with an optimized high-throughput Data Independent Acquisition (DIA) approach with the throughput of 60 sample per day, we characterized the proteomic landscape of DRCs in comparison to parental PC3 cells. This optimized DIA method allowed for robust and reproducible protein quantification at the single-cell level, enabling the identification and quantification of over 1,300 proteins per cell on average. Distinct proteomic sub-clusters within the DRC population were identified, closely linked to variations in cell size. The study uncovered novel protein signatures, including the regulation of proteins critical for cell adhesion and metabolic processes, as well as the upregulation of surface proteins and transcription factors pivotal for cancer progression. Furthermore, by integrating SCP and single-cell RNA-seq (scRNA-seq) data, we identified six upregulated and ten downregulated genes consistently altered in drug-treated cells across both SCP and scRNA-seq platforms. These findings underscore the heterogeneity of DRCs and their unique molecular signatures, providing valuable insights into their biological behavior and potential therapeutic targets.
Project description:MAP4K family kinases are key kinases for T-cell-mediated immune responses; however, in vivo roles of MAP4K2 in immune regulation remain unclear. Using T-cell-specific Map4k2 conditional knockout (T-Map4k2 cKO) mice, single-cell RNA sequencing (scRNA-seq), and mass spectrometry analysis, we found that MAP4K2 interacted with DDX39B, induced FOXP3 gene expression, and promoted Treg differentiation. Mechanistically, MAP4K2 directly phosphorylated the DEAD box protein DDX39B, leading to DDX39B nuclear translocation and subsequent Foxp3 RNA splicing. MAP4K2-induced FOXP3 mRNA levels were abolished in DDX39B knockout T cells. Furthermore, T-Map4k2 cKO mice displayed the reduction of Treg population and sustained inflammation during remission phase of EAE autoimmune disease model. Remarkably, the anti-PD-1 immunotherapeutic effect on pancreatic cancer was drastically improved in T-Map4k2 cKO mice, Treg-specific Map4k2-deficient mice, adaptively transferred mice, or MAP4K2-inhibitor-treated mice. Consistently, scRNA-seq analysis of human pancreatic patients showed increased MAP4K2 levels in infiltrating Treg cells. Collectively, MAP4K2 promotes Treg differentiation by inducing DDX39B nuclear translocation, leading to the attenuation of tumor immunity.
Project description:MAP4K family kinases are key kinases for T-cell-mediated immune responses; however, in vivo roles of MAP4K2 in immune regulation remain unclear. Using T- cell-specific Map4k2 conditional knockout (T-Map4k2 cKO) mice, single-cell RNA sequencing (scRNA-seq), and mass spectrometry analysis, we found that MAP4K2 interacted with DDX39B, induced FOXP3 gene expression, and promoted Treg differentiation. Mechanistically, MAP4K2 directly phosphorylated the DEAD box protein DDX39B, leading to DDX39B nuclear translocation and subsequent Foxp3 RNA splicing. MAP4K2-induced FOXP3 mRNA levels were abolished in DDX39B knockout T cells. Furthermore, T-Map4k2 cKO mice displayed the reduction of Treg population and sustained inflammation during remission phase of EAE autoimmune disease model. Remarkably, the anti-PD-1 immunotherapeutic effect on pancreatic cancer was drastically improved in T-Map4k2 cKO mice, Treg-specific Map4k2- deficient mice, or MAP4K2-inhibitor-treated mice. Consistently, scRNA-seq analysis of human pancreatic or lung cancer patients showed increased MAP4K2 levels in infiltrating Treg cells. Collectively, MAP4K2 promotes Treg differentiation by inducing DDX39B nuclear translocation, leading to the attenuation of tumor immunity.
Project description:Aging is a universal biological phenomenon linked to many diseases, such as cancer or neurodegeneration. However, the molecular mechanisms underlying aging, or how lifestyle interventions such as cognitive stimulation can ameliorate this process, are yet to be clarified. Here, we performed a multi-omic profiling, including RNA-seq, ATAC-seq, ChIP-seq, EM-seq, SWATH-MS and single cell Multiome scRNA and scATAC-seq, in the dorsal hippocampus of young and old mouse subjects which were subject to cognitive stimulation using the paradigm of environmental enrichment. In this study we were able to describe the epigenomic landscape of aging and cognitive stimulation.
Project description:Aneuploidy is a frequent feature of human tumors. Germline mutations leading to aneuploidy are very rare in humans, and their tumor-promoting properties are mostly unknown at the molecular level. We report here novel germline biallelic mutations in MAD1L1, the gene encoding the Spindle Assembly Checkpoint (SAC) protein MAD1, in a 36-year-old female with a dozen of neoplasias, including five malignant tumors. Functional studies in peripheral blood cells demonstrated lack of full-length protein and deficient SAC response, resulting in ~30-40% of aneuploid cells as detected by cytogenetic and single-cell (sc) DNA analysis. scRNA-seq analysis of proband blood cells identified mitochondrial stress accompanied by systemic inflammation with enhanced interferon and NFkB signaling. The inference of chromosomal aberrations from scRNA-seq analysis detected inflammatory signals both in aneuploid and euploid cells, suggesting a non-cell autonomous response to aneuploidy. In addition to random aneuploidies, MAD1L1 mutations resulted in specific clonal expansions of T-cells with chromosome 18 gains and enhanced cytotoxic profile, as well as intermediate B-cells with chromosome 12 gains and transcriptomic signatures characteristic of chronic lymphocytic leukemia cells. These data point to MAD1L1 mutations as the cause of a new variant of mosaic variegated aneuploidy syndrome (MVA) with systemic inflammation and unprecedented tumor susceptibility.
Project description:Genetic and non-genetic heterogeneity within cancer cell populations represents a major challenge to anti-cancer therapies. We currently lack robust methods to determine how pre-existing and adaptive features affect cellular responses to therapies. Here, by conducting clonal fitness mapping and transcriptional characterization using expressed barcodes and single-cell RNA-sequencing, we have developed TraCe-seq, a method that captures at clonal resolution the origin, fate, and differential early adaptive transcriptional programs of cells in a complex population in response to distinct treatments. We used TraCe-seq to benchmark how next-generation dual EGFR inhibitors-degraders compare to standard EGFR kinase inhibitors in EGFR-mutant lung cancer cells. To QC the TraCe-seq strategy, single-cell RNA-seq libraries were generated from a variety of human cancer cell lines transduced with the TraCe-seq library to validate the TraCe-seq strategy. Specifically, 5 different cell lines (PC9, MCF-10A, MDA-MB-231, NCI-H358, and NCI-H1373) were each transduced with a unique TraCe-seq barcode. The transduced cells were selected with puromycin only, dissociated to single cell suspensions, and then mixed together. The complex mixture of the 5 cell lines was profiled by 10X scRNA-seq. Furthermore, transduced NCI-H1373 cells were sorted by FACS to enrich for the top 50% of eGFP positive cells, and sorted cells were cultured briefly and used to construct scRNA-seq libraries and profiled by 10x scRNA-seq. To carry out the full TraCe-seq experiment, ~600 PC9 cells carrying unique TraCe-seq barcodes were expanded over 12 doublings to establish the barcoded population. A subset of the barcoded PC9 population was used to generate scRNA-seq libraries and profiled by 10x scRNA-seq prior to treatment to establish a baseline transcription profile for each barcoded clone. The rest of the cells were then treated for four days with 1 µM erlotinib, 1 µM GNE-069, or 1 µM GNE-104 respectively. scRNA-seq libraries were then generated form the treated cells and profiled by 10x scRNA-seq.