1
|
Zhu L, Yang S, Zhang K, Wang H, Fang X, Wang J. Uncovering underlying physical principles and driving forces of cell differentiation and reprogramming from single-cell transcriptomics. Proc Natl Acad Sci U S A 2024; 121:e2401540121. [PMID: 39150785 DOI: 10.1073/pnas.2401540121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Accepted: 06/28/2024] [Indexed: 08/18/2024] Open
Abstract
Recent advances in single-cell sequencing technology have revolutionized our ability to acquire whole transcriptome data. However, uncovering the underlying transcriptional drivers and nonequilibrium driving forces of cell function directly from these data remains challenging. We address this by learning cell state vector fields from discrete single-cell RNA velocity to quantify the single-cell global nonequilibrium driving forces as landscape and flux. From single-cell data, we quantified the Waddington landscape, showing that optimal paths for differentiation and reprogramming deviate from the naively expected landscape gradient paths and may not pass through landscape saddles at finite fluctuations, challenging conventional transition state estimation of kinetic rate for cell fate decisions due to the presence of the flux. A key insight from our study is that stem/progenitor cells necessitate greater energy dissipation for rapid cell cycles and self-renewal, maintaining pluripotency. We predict optimal developmental pathways and elucidate the nucleation mechanism of cell fate decisions, with transition states as nucleation sites and pioneer genes as nucleation seeds. The concept of loop flux quantifies the contributions of each cycle flux to cell state transitions, facilitating the understanding of cell dynamics and thermodynamic cost, and providing insights into optimizing biological functions. We also infer cell-cell interactions and cell-type-specific gene regulatory networks, encompassing feedback mechanisms and interaction intensities, predicting genetic perturbation effects on cell fate decisions from single-cell omics data. Essentially, our methodology validates the landscape and flux theory, along with its associated quantifications, offering a framework for exploring the physical principles underlying cellular differentiation and reprogramming and broader biological processes through high-throughput single-cell sequencing experiments.
Collapse
Affiliation(s)
- Ligang Zhu
- College of Physics, Jilin University, Changchun 130021, China
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, China
| | - Songlin Yang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, China
| | - Kun Zhang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, China
| | - Hong Wang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, China
| | - Xiaona Fang
- College of Chemistry, Northeast Normal University, Changchun 130024, China
| | - Jin Wang
- Center for Theoretical Interdisciplinary Sciences, Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou 325001, China
- Department of Chemistry, Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794
| |
Collapse
|
2
|
Bonham-Carter B, Schiebinger G. Cellular proliferation biases clonal lineage tracing and trajectory inference. Bioinformatics 2024; 40:btae483. [PMID: 39102821 PMCID: PMC11316616 DOI: 10.1093/bioinformatics/btae483] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2023] [Revised: 03/28/2024] [Accepted: 07/30/2024] [Indexed: 08/07/2024] Open
Abstract
MOTIVATION Lineage tracing and trajectory inference from single-cell RNA-sequencing data hold tremendous potential for uncovering the genetic programs driving development and disease. Single cell datasets are thought to provide an unbiased view on the diverse cellular architecture of tissues. Sampling bias, however, can skew single cell datasets away from the cellular composition they are meant to represent. RESULTS We demonstrate a novel form of sampling bias, caused by a statistical phenomenon related to repeated sampling from a growing, heterogeneous population. Relative growth rates of cells influence the probability that they will be sampled in clones observed across multiple time points. We support our probabilistic derivations with a simulation study and an analysis of a real time-course of T-cell development. We find that this bias can impact fate probability predictions, and we explore how to develop trajectory inference methods which are robust to this bias. AVAILABILITY AND IMPLEMENTATION Source code for the simulated datasets and to create the figures in this manuscript is freely available in python at https://github.com/rbonhamcarter/simulate-clones. A python implementation of the extension of the LineageOT method is freely available at https://github.com/rbonhamcarter/LineageOT/tree/multi-time-clones.
Collapse
Affiliation(s)
- Becca Bonham-Carter
- Department of Mathematics, University of British Columbia, Vancouver, BC V6T 1Z4, Canada
| | - Geoffrey Schiebinger
- Department of Mathematics, University of British Columbia, Vancouver, BC V6T 1Z4, Canada
| |
Collapse
|
3
|
Li K, Wang R, Liu GW, Peng ZY, Wang JC, Xiao GD, Tang SC, Du N, Zhang J, Zhang J, Ren H, Sun X, Yang YP, Liu DP. Refining the optimal CAF cluster marker for predicting TME-dependent survival expectancy and treatment benefits in NSCLC patients. Sci Rep 2024; 14:16766. [PMID: 39034310 PMCID: PMC11271481 DOI: 10.1038/s41598-024-55375-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Accepted: 02/22/2024] [Indexed: 07/23/2024] Open
Abstract
The tumor microenvironment (TME) plays a pivotal role in the onset, progression, and treatment response of cancer. Among the various components of the TME, cancer-associated fibroblasts (CAFs) are key regulators of both immune and non-immune cellular functions. Leveraging single-cell RNA sequencing (scRNA) data, we have uncovered previously hidden and promising roles within this specific CAF subgroup, paving the way for its clinical application. However, several critical questions persist, primarily stemming from the heterogeneous nature of CAFs and the use of different fibroblast markers in various sample analyses, causing confusion and hindrance in their clinical implementation. In this groundbreaking study, we have systematically screened multiple databases to identify the most robust marker for distinguishing CAFs in lung cancer, with a particular focus on their potential use in early diagnosis, staging, and treatment response evaluation. Our investigation revealed that COL1A1, COL1A2, FAP, and PDGFRA are effective markers for characterizing CAF subgroups in most lung adenocarcinoma datasets. Through comprehensive analysis of treatment responses, we determined that COL1A1 stands out as the most effective indicator among all CAF markers. COL1A1 not only deciphers the TME signatures related to CAFs but also demonstrates a highly sensitive and specific correlation with treatment responses and multiple survival outcomes. For the first time, we have unveiled the distinct roles played by clusters of CAF markers in differentiating various TME groups. Our findings confirm the sensitive and unique contributions of CAFs to the responses of multiple lung cancer therapies. These insights significantly enhance our understanding of TME functions and drive the translational application of extensive scRNA sequence results. COL1A1 emerges as the most sensitive and specific marker for defining CAF subgroups in scRNA analysis. The CAF ratios represented by COL1A1 can potentially serve as a reliable predictor of treatment responses in clinical practice, thus providing valuable insights into the influential roles of TME components. This research marks a crucial step forward in revolutionizing our approach to cancer diagnosis and treatment.
Collapse
Affiliation(s)
- Kai Li
- Department of Otorhinolaryngology‑Head and Neck Surgery, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, 710061, Shaanxi, China
| | - Rui Wang
- Department of Thoracic Surgery and Oncology, Cancer Centre, The First Affiliated Hospital of Xi'an Jiaotong University, 277 Yanta West Road, Xi'an, 710061, Shaanxi, China
| | - Guo-Wei Liu
- Department of Thoracic Surgery, Qinghai Provincial People's Hospital, Gonghe Road No. 2, Chengdong District, Xining, 810007, Qinghai, China
| | - Zi-Yang Peng
- School of Future Technology, National Local Joint Engineering Research Center for Precision Surgery and Regenerative Medicine, Xi'an Jiaotong University, Xi'an, 710061, Shaanxi, China
| | - Ji-Chang Wang
- Department of Vascular Surgery, First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, 710061, Shaanxi, China
| | - Guo-Dong Xiao
- Oncology Department, The First Affiliated Hospital of Zhengzhou University, Zheng Zhou, 450052, Henan, China
| | - Shou-Ching Tang
- Section of Hematology Oncology, Department of Internal Medicine, LSUHSC Cancer Center, School of Medicine, 1700 Tulane Avenue, New Orleans, LA, 70112, USA
| | - Ning Du
- Department of Thoracic Surgery and Oncology, Cancer Centre, The First Affiliated Hospital of Xi'an Jiaotong University, 277 Yanta West Road, Xi'an, 710061, Shaanxi, China
| | - Jia Zhang
- Department of Thoracic Surgery and Oncology, Cancer Centre, The First Affiliated Hospital of Xi'an Jiaotong University, 277 Yanta West Road, Xi'an, 710061, Shaanxi, China
| | - Jing Zhang
- Department of Thoracic Surgery and Oncology, Cancer Centre, The First Affiliated Hospital of Xi'an Jiaotong University, 277 Yanta West Road, Xi'an, 710061, Shaanxi, China
| | - Hong Ren
- Department of Thoracic Surgery and Oncology, Cancer Centre, The First Affiliated Hospital of Xi'an Jiaotong University, 277 Yanta West Road, Xi'an, 710061, Shaanxi, China
| | - Xin Sun
- Department of Thoracic Surgery and Oncology, Cancer Centre, The First Affiliated Hospital of Xi'an Jiaotong University, 277 Yanta West Road, Xi'an, 710061, Shaanxi, China
- Department of Pathology, Anatomy and Cell Biology, Sidney Kimmel Cancer Center, Thomas Jefferson University, Philadelphia, PA, 19107, USA
| | - Yi-Ping Yang
- Department of Radiotherapy, Shaanxi Provincial Tumor Hospital, 309 Yanta W Rd, Yanta District, Xi'an, 710063, Shaanxi, China.
| | - Da-Peng Liu
- Department of Thoracic Surgery and Oncology, Cancer Centre, The First Affiliated Hospital of Xi'an Jiaotong University, 277 Yanta West Road, Xi'an, 710061, Shaanxi, China.
| |
Collapse
|
4
|
Coulombe B, Chapleau A, Macintosh J, Durcan TM, Poitras C, Moursli YA, Faubert D, Pinard M, Bernard G. Towards a Treatment for Leukodystrophy Using Cell-Based Interception and Precision Medicine. Biomolecules 2024; 14:857. [PMID: 39062571 PMCID: PMC11274857 DOI: 10.3390/biom14070857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2024] [Revised: 07/12/2024] [Accepted: 07/14/2024] [Indexed: 07/28/2024] Open
Abstract
Cell-based interception and precision medicine is a novel approach aimed at improving healthcare through the early detection and treatment of diseased cells. Here, we describe our recent progress towards developing cell-based interception and precision medicine to detect, understand, and advance the development of novel therapeutic approaches through a single-cell omics and drug screening platform, as part of a multi-laboratory collaborative effort, for a group of neurodegenerative disorders named leukodystrophies. Our strategy aims at the identification of diseased cells as early as possible to intercept progression of the disease prior to severe clinical impairment and irreversible tissue damage.
Collapse
Affiliation(s)
- Benoit Coulombe
- Translational Proteomics Laboratory, Institut de Recherches Cliniques de Montréal, Montréal, QC H2W 1R7, Canada; (C.P.); (Y.A.M.); (M.P.)
- Department of Biochemistry and Molecular Medicine, Université de Montréal, Montréal, QC H3T 1A8, Canada
| | - Alexandra Chapleau
- Department of Neurology and Neurosurgery, Pediatrics and Human Genetics, McGill University, Montréal, QC H9X 3V9, Canada; (A.C.); (J.M.); (G.B.)
- Child Health and Human Development Program, Research Institute of the McGill University Health Centre, Montréal, QC H4A 3J1, Canada
- The Neuro’s Early Drug Discovery Unit (EDDU), McGill University, Montréal, QC H9X 3V9, Canada;
| | - Julia Macintosh
- Department of Neurology and Neurosurgery, Pediatrics and Human Genetics, McGill University, Montréal, QC H9X 3V9, Canada; (A.C.); (J.M.); (G.B.)
- Child Health and Human Development Program, Research Institute of the McGill University Health Centre, Montréal, QC H4A 3J1, Canada
| | - Thomas M. Durcan
- The Neuro’s Early Drug Discovery Unit (EDDU), McGill University, Montréal, QC H9X 3V9, Canada;
| | - Christian Poitras
- Translational Proteomics Laboratory, Institut de Recherches Cliniques de Montréal, Montréal, QC H2W 1R7, Canada; (C.P.); (Y.A.M.); (M.P.)
| | - Yena A. Moursli
- Translational Proteomics Laboratory, Institut de Recherches Cliniques de Montréal, Montréal, QC H2W 1R7, Canada; (C.P.); (Y.A.M.); (M.P.)
| | - Denis Faubert
- Mass Spectrometry and Proteomics Platform, Institut de Recherches Cliniques de Montréal, Montréal, QC H2W 1R7, Canada;
| | - Maxime Pinard
- Translational Proteomics Laboratory, Institut de Recherches Cliniques de Montréal, Montréal, QC H2W 1R7, Canada; (C.P.); (Y.A.M.); (M.P.)
| | - Geneviève Bernard
- Department of Neurology and Neurosurgery, Pediatrics and Human Genetics, McGill University, Montréal, QC H9X 3V9, Canada; (A.C.); (J.M.); (G.B.)
- Child Health and Human Development Program, Research Institute of the McGill University Health Centre, Montréal, QC H4A 3J1, Canada
- Department Specialized Medicine, Division of Medical Genetics, McGill University Health Centre, Montréal, QC H4A 3J1, Canada
| |
Collapse
|
5
|
Rukhlenko OS, Imoto H, Tambde A, McGillycuddy A, Junk P, Tuliakova A, Kolch W, Kholodenko BN. Cell State Transition Models Stratify Breast Cancer Cell Phenotypes and Reveal New Therapeutic Targets. Cancers (Basel) 2024; 16:2354. [PMID: 39001416 PMCID: PMC11240448 DOI: 10.3390/cancers16132354] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2024] [Revised: 06/17/2024] [Accepted: 06/23/2024] [Indexed: 07/16/2024] Open
Abstract
Understanding signaling patterns of transformation and controlling cell phenotypes is a challenge of current biology. Here we applied a cell State Transition Assessment and Regulation (cSTAR) approach to a perturbation dataset of single cell phosphoproteomic patterns of multiple breast cancer (BC) and normal breast tissue-derived cell lines. Following a separation of luminal, basal, and normal cell states, we identified signaling nodes within core control networks, delineated causal connections, and determined the primary drivers underlying oncogenic transformation and transitions across distinct BC subtypes. Whereas cell lines within the same BC subtype have different mutational and expression profiles, the architecture of the core network was similar for all luminal BC cells, and mTOR was a main oncogenic driver. In contrast, core networks of basal BC were heterogeneous and segregated into roughly four major subclasses with distinct oncogenic and BC subtype drivers. Likewise, normal breast tissue cells were separated into two different subclasses. Based on the data and quantified network topologies, we derived mechanistic cSTAR models that serve as digital cell twins and allow the deliberate control of cell movements within a Waddington landscape across different cell states. These cSTAR models suggested strategies of normalizing phosphorylation networks of BC cell lines using small molecule inhibitors.
Collapse
Affiliation(s)
- Oleksii S Rukhlenko
- Systems Biology Ireland, School of Medicine, University College Dublin, D04 V1W8 Dublin, Ireland
| | - Hiroaki Imoto
- Systems Biology Ireland, School of Medicine, University College Dublin, D04 V1W8 Dublin, Ireland
| | - Ayush Tambde
- Systems Biology Ireland, School of Medicine, University College Dublin, D04 V1W8 Dublin, Ireland
- Stratford College, D06 T9V3 Dublin, Ireland
| | - Amy McGillycuddy
- Systems Biology Ireland, School of Medicine, University College Dublin, D04 V1W8 Dublin, Ireland
- School of Biological, Health and Sports Sciences, Technological University, D07 H6K8 Dublin, Ireland
| | - Philipp Junk
- Systems Biology Ireland, School of Medicine, University College Dublin, D04 V1W8 Dublin, Ireland
| | - Anna Tuliakova
- Systems Biology Ireland, School of Medicine, University College Dublin, D04 V1W8 Dublin, Ireland
| | - Walter Kolch
- Systems Biology Ireland, School of Medicine, University College Dublin, D04 V1W8 Dublin, Ireland
- Conway Institute of Biomolecular and Biomedical Research, University College Dublin, D04 V1W8 Dublin, Ireland
| | - Boris N Kholodenko
- Systems Biology Ireland, School of Medicine, University College Dublin, D04 V1W8 Dublin, Ireland
- Conway Institute of Biomolecular and Biomedical Research, University College Dublin, D04 V1W8 Dublin, Ireland
- Department of Pharmacology, Yale University School of Medicine, New Haven, CT 06520, USA
| |
Collapse
|
6
|
Peng D, Cahan P. OneSC: A computational platform for recapitulating cell state transitions. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.31.596831. [PMID: 38895453 PMCID: PMC11185539 DOI: 10.1101/2024.05.31.596831] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
Computational modelling of cell state transitions has been a great interest of many in the field of developmental biology, cancer biology and cell fate engineering because it enables performing perturbation experiments in silico more rapidly and cheaply than could be achieved in a wet lab. Recent advancements in single-cell RNA sequencing (scRNA-seq) allow the capture of high-resolution snapshots of cell states as they transition along temporal trajectories. Using these high-throughput datasets, we can train computational models to generate in silico 'synthetic' cells that faithfully mimic the temporal trajectories. Here we present OneSC, a platform that can simulate synthetic cells across developmental trajectories using systems of stochastic differential equations govern by a core transcription factors (TFs) regulatory network. Different from the current network inference methods, OneSC prioritizes on generating Boolean network that produces faithful cell state transitions and steady cell states that mimic real biological systems. Applying OneSC to real data, we inferred a core TF network using a mouse myeloid progenitor scRNA-seq dataset and showed that the dynamical simulations of that network generate synthetic single-cell expression profiles that faithfully recapitulate the four myeloid differentiation trajectories going into differentiated cell states (erythrocytes, megakaryocytes, granulocytes and monocytes). Finally, through the in-silico perturbations of the mouse myeloid progenitor core network, we showed that OneSC can accurately predict cell fate decision biases of TF perturbations that closely match with previous experimental observations.
Collapse
Affiliation(s)
- Da Peng
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, 21205, USA
| | - Patrick Cahan
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, 21205, USA
- Institute for Cell Engineering, Johns Hopkins University, Baltimore, Maryland, 21205, USA
- Department of Molecular Biology and Genetics, Johns Hopkins University, Baltimore, Maryland, 21205, USA
| |
Collapse
|
7
|
Rivero-Garcia I, Torres M, Sánchez-Cabo F. Deep generative models in single-cell omics. Comput Biol Med 2024; 176:108561. [PMID: 38749321 DOI: 10.1016/j.compbiomed.2024.108561] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2024] [Revised: 04/30/2024] [Accepted: 05/05/2024] [Indexed: 05/31/2024]
Abstract
Deep Generative Models (DGMs) are becoming instrumental for inferring probability distributions inherent to complex processes, such as most questions in biomedical research. For many years, there was a lack of mathematical methods that would allow this inference in the scarce data scenario of biomedical research. The advent of single-cell omics has finally made square the so-called "skinny matrix", allowing to apply mathematical methods already extensively used in other areas. Moreover, it is now possible to integrate data at different molecular levels in thousands or even millions of samples, thanks to the number of single-cell atlases being collaboratively generated. Additionally, DGMs have proven useful in other frequent tasks in single-cell analysis pipelines, from dimensionality reduction, cell type annotation to RNA velocity inference. In spite of its promise, DGMs need to be used with caution in biomedical research, paying special attention to its use to answer the right questions and the definition of appropriate error metrics and validation check points that confirm not only its correct use but also its relevance. All in all, DGMs provide an exciting tool that opens a bright future for the integrative analysis of single-cell -omics to understand health and disease.
Collapse
Affiliation(s)
- Inés Rivero-Garcia
- Universidad Politécnica de Madrid, Madrid, 28040, Spain; Centro Nacional de Investigaciones Cardiovasculares (CNIC), Madrid, 28029, Spain
| | - Miguel Torres
- Centro Nacional de Investigaciones Cardiovasculares (CNIC), Madrid, 28029, Spain
| | - Fátima Sánchez-Cabo
- Centro Nacional de Investigaciones Cardiovasculares (CNIC), Madrid, 28029, Spain.
| |
Collapse
|
8
|
Hossain I, Fanfani V, Fischer J, Quackenbush J, Burkholz R. Biologically informed NeuralODEs for genome-wide regulatory dynamics. Genome Biol 2024; 25:127. [PMID: 38773638 PMCID: PMC11106922 DOI: 10.1186/s13059-024-03264-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2023] [Accepted: 04/30/2024] [Indexed: 05/24/2024] Open
Abstract
BACKGROUND Gene regulatory network (GRN) models that are formulated as ordinary differential equations (ODEs) can accurately explain temporal gene expression patterns and promise to yield new insights into important cellular processes, disease progression, and intervention design. Learning such gene regulatory ODEs is challenging, since we want to predict the evolution of gene expression in a way that accurately encodes the underlying GRN governing the dynamics and the nonlinear functional relationships between genes. Most widely used ODE estimation methods either impose too many parametric restrictions or are not guided by meaningful biological insights, both of which impede either scalability, explainability, or both. RESULTS We developed PHOENIX, a modeling framework based on neural ordinary differential equations (NeuralODEs) and Hill-Langmuir kinetics, that overcomes limitations of other methods by flexibly incorporating prior domain knowledge and biological constraints to promote sparse, biologically interpretable representations of GRN ODEs. We tested the accuracy of PHOENIX in a series of in silico experiments, benchmarking it against several currently used tools. We demonstrated PHOENIX's flexibility by modeling regulation of oscillating expression profiles obtained from synchronized yeast cells. We also assessed the scalability of PHOENIX by modeling genome-scale GRNs for breast cancer samples ordered in pseudotime and for B cells treated with Rituximab. CONCLUSIONS PHOENIX uses a combination of user-defined prior knowledge and functional forms from systems biology to encode biological "first principles" as soft constraints on the GRN allowing us to predict subsequent gene expression patterns in a biologically explainable manner.
Collapse
Affiliation(s)
| | - Viola Fanfani
- Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Jonas Fischer
- Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | | | - Rebekka Burkholz
- CISPA Helmholtz Center for Information Security, Saarbrücken, Germany
| |
Collapse
|
9
|
Hong T, Xing J. Data- and theory-driven approaches for understanding paths of epithelial-mesenchymal transition. Genesis 2024; 62:e23591. [PMID: 38553870 PMCID: PMC11017362 DOI: 10.1002/dvg.23591] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2024] [Revised: 02/16/2024] [Accepted: 03/16/2024] [Indexed: 04/02/2024]
Abstract
Reversible transitions between epithelial and mesenchymal cell states are a crucial form of epithelial plasticity for development and disease progression. Recent experimental data and mechanistic models showed multiple intermediate epithelial-mesenchymal transition (EMT) states as well as trajectories of EMT underpinned by complex gene regulatory networks. In this review, we summarize recent progress in quantifying EMT and characterizing EMT paths with computational methods and quantitative experiments including omics-level measurements. We provide perspectives on how these studies can help relating fundamental cell biology to physiological and pathological outcomes of EMT.
Collapse
Affiliation(s)
- Tian Hong
- Department of Biochemistry & Cellular and Molecular Biology, The University of Tennessee, Knoxville, Knoxville TN, USA
| | - Jianhua Xing
- Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, PA, USA
- UPMC-Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA, USA
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA, USA
| |
Collapse
|
10
|
Bailly R, Malfante M, Allier C, Paviolo C, Ghenim L, Padmanabhan K, Bardin S, Mars J. Detecting abnormal cell behaviors from dry mass time series. Sci Rep 2024; 14:7053. [PMID: 38528035 DOI: 10.1038/s41598-024-57684-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Accepted: 03/20/2024] [Indexed: 03/27/2024] Open
Abstract
The prediction of pathological changes on single cell behaviour is a challenging task for deep learning models. Indeed, in self-supervised learning methods, no prior labels are used for the training and all of the information for event predictions are extracted from the data themselves. We present here a novel self-supervised learning model for the detection of anomalies in a given cell population, StArDusTS. Cells are monitored over time, and analysed to extract time-series of dry mass values. We assessed its performances on different cell lines, showing a precision of 96% in the automatic detection of anomalies. Additionally, anomaly detection was also associated with cell measurement errors inherent to the acquisition or analysis pipelines, leading to an improvement of the upstream methods for feature extraction. Our results pave the way to novel architectures for the continuous monitoring of cell cultures in applied research or bioproduction applications, and for the prediction of pathological cellular changes.
Collapse
Affiliation(s)
- Romain Bailly
- Univ. Grenoble Alpes, CEA, List, F-38000, Grenoble, France
- Univ. Grenoble Alpes, CNRS, Grenoble-INP, GIPSA-Lab, 38000, Grenoble, France
| | | | - Cédric Allier
- Univ. Grenoble Alpes, CEA, Leti, F-38000, Grenoble, France
- Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, VA, USA
| | - Chiara Paviolo
- Univ. Grenoble Alpes, CEA, Leti, F-38000, Grenoble, France
| | - Lamya Ghenim
- Univ. Grenoble Alpes, INSERM, CEA-IRIG, BGE, Biomics, F-38000, Grenoble, France
| | - Kiran Padmanabhan
- Institut de Génomique Fonctionnelle de Lyon, Univ. Lyon, CNRS/ENS, UMR 5242, Lyon, France
| | - Sabine Bardin
- Institut Curie, PSL Research University, CNRS, UMR 144, Molecular Mechanisms of Intracellular Transport, F-75005, Paris, France
| | - Jérôme Mars
- Univ. Grenoble Alpes, CNRS, Grenoble-INP, GIPSA-Lab, 38000, Grenoble, France
| |
Collapse
|
11
|
Baghdassarian HM, Lewis NE. Resource allocation in mammalian systems. Biotechnol Adv 2024; 71:108305. [PMID: 38215956 PMCID: PMC11182366 DOI: 10.1016/j.biotechadv.2023.108305] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Revised: 12/17/2023] [Accepted: 12/18/2023] [Indexed: 01/14/2024]
Abstract
Cells execute biological functions to support phenotypes such as growth, migration, and secretion. Complementarily, each function of a cell has resource costs that constrain phenotype. Resource allocation by a cell allows it to manage these costs and optimize their phenotypes. In fact, the management of resource constraints (e.g., nutrient availability, bioenergetic capacity, and macromolecular machinery production) shape activity and ultimately impact phenotype. In mammalian systems, quantification of resource allocation provides important insights into higher-order multicellular functions; it shapes intercellular interactions and relays environmental cues for tissues to coordinate individual cells to overcome resource constraints and achieve population-level behavior. Furthermore, these constraints, objectives, and phenotypes are context-dependent, with cells adapting their behavior according to their microenvironment, resulting in distinct steady-states. This review will highlight the biological insights gained from probing resource allocation in mammalian cells and tissues.
Collapse
Affiliation(s)
- Hratch M Baghdassarian
- Bioinformatics and Systems Biology Graduate Program, University of California, San Diego, La Jolla, CA 92093, USA; Department of Pediatrics, University of California, San Diego, La Jolla, CA 92093, USA
| | - Nathan E Lewis
- Department of Pediatrics, University of California, San Diego, La Jolla, CA 92093, USA; Department of Bioengineering, University of California, San Diego, La Jolla, CA 92093, USA.
| |
Collapse
|
12
|
Kucinski I, Campos J, Barile M, Severi F, Bohin N, Moreira PN, Allen L, Lawson H, Haltalli MLR, Kinston SJ, O'Carroll D, Kranc KR, Göttgens B. A time- and single-cell-resolved model of murine bone marrow hematopoiesis. Cell Stem Cell 2024; 31:244-259.e10. [PMID: 38183977 PMCID: PMC7615671 DOI: 10.1016/j.stem.2023.12.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2023] [Revised: 09/25/2023] [Accepted: 12/04/2023] [Indexed: 01/08/2024]
Abstract
The paradigmatic hematopoietic tree model is increasingly recognized to be limited, as it is based on heterogeneous populations largely defined by non-homeostatic assays testing cell fate potentials. Here, we combine persistent labeling with time-series single-cell RNA sequencing to build a real-time, quantitative model of in vivo tissue dynamics for murine bone marrow hematopoiesis. We couple cascading single-cell expression patterns with dynamic changes in differentiation and growth speeds. The resulting explicit linkage between molecular states and cellular behavior reveals widely varying self-renewal and differentiation properties across distinct lineages. Transplanted stem cells show strong acceleration of differentiation at specific stages of erythroid and neutrophil production, illustrating how the model can quantify the impact of perturbations. Our reconstruction of dynamic behavior from snapshot measurements is akin to how a kinetoscope allows sequential images to merge into a movie. We posit that this approach is generally applicable to understanding tissue-scale dynamics at high resolution.
Collapse
Affiliation(s)
- Iwo Kucinski
- Wellcome-MRC Cambridge Stem Cell Institute, Department of Haematology, Jeffrey Cheah Biomedical Centre, University of Cambridge, Cambridge, UK
| | - Joana Campos
- Centre for Haemato-Oncology, Barts Cancer Institute, Queen Mary University of London, London EC1M 6BQ, UK; Institute of Cancer Research, London SM2 5NG, UK
| | - Melania Barile
- Wellcome-MRC Cambridge Stem Cell Institute, Department of Haematology, Jeffrey Cheah Biomedical Centre, University of Cambridge, Cambridge, UK; Centre for Translational Stem Cell Biology, Hong Kong SAR, China
| | - Francesco Severi
- Centre for Regenerative Medicine, University of Edinburgh, Edinburgh EH16 4UU, UK; Wellcome Centre for Cell Biology, University of Edinburgh, Edinburgh EH9 3BF, UK
| | - Natacha Bohin
- Centre for Haemato-Oncology, Barts Cancer Institute, Queen Mary University of London, London EC1M 6BQ, UK
| | - Pedro N Moreira
- Centre for Regenerative Medicine, University of Edinburgh, Edinburgh EH16 4UU, UK; Wellcome Centre for Cell Biology, University of Edinburgh, Edinburgh EH9 3BF, UK
| | - Lewis Allen
- Centre for Haemato-Oncology, Barts Cancer Institute, Queen Mary University of London, London EC1M 6BQ, UK; Institute of Cancer Research, London SM2 5NG, UK
| | - Hannah Lawson
- Centre for Haemato-Oncology, Barts Cancer Institute, Queen Mary University of London, London EC1M 6BQ, UK; Institute of Cancer Research, London SM2 5NG, UK
| | - Myriam L R Haltalli
- Wellcome-MRC Cambridge Stem Cell Institute, Department of Haematology, Jeffrey Cheah Biomedical Centre, University of Cambridge, Cambridge, UK
| | - Sarah J Kinston
- Wellcome-MRC Cambridge Stem Cell Institute, Department of Haematology, Jeffrey Cheah Biomedical Centre, University of Cambridge, Cambridge, UK
| | - Dónal O'Carroll
- Centre for Regenerative Medicine, University of Edinburgh, Edinburgh EH16 4UU, UK; Wellcome Centre for Cell Biology, University of Edinburgh, Edinburgh EH9 3BF, UK.
| | - Kamil R Kranc
- Centre for Haemato-Oncology, Barts Cancer Institute, Queen Mary University of London, London EC1M 6BQ, UK; Institute of Cancer Research, London SM2 5NG, UK.
| | - Berthold Göttgens
- Wellcome-MRC Cambridge Stem Cell Institute, Department of Haematology, Jeffrey Cheah Biomedical Centre, University of Cambridge, Cambridge, UK.
| |
Collapse
|
13
|
Jumde G, Spanjaard B, Junker JP. Inference of differentiation trajectories by transfer learning across biological processes. Cell Syst 2024; 15:75-82.e5. [PMID: 38128536 DOI: 10.1016/j.cels.2023.12.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2023] [Revised: 07/28/2023] [Accepted: 12/06/2023] [Indexed: 12/23/2023]
Abstract
Stem cells differentiate into distinct fates by transitioning through a series of transcriptional states. Current computational approaches allow reconstruction of differentiation trajectories from single-cell transcriptomics data, but it remains unknown to what degree differentiation can be predicted across biological processes. Here, we use transfer learning to infer differentiation processes and quantify predictability in early embryonic development and adult hematopoiesis. Overall, we find that non-linear methods outperform linear approaches, and we achieved the best predictions with a custom variational autoencoder that explicitly models changes in transcriptional variance. We observed a high accuracy of predictions in embryonic development, but we found somewhat lower agreement with the real data in adult hematopoiesis. We demonstrate that this discrepancy can be explained by a higher degree of concordant transcriptional processes along embryonic differentiation compared with adult homeostasis. In summary, we establish a framework for quantifying and exploiting predictability of cellular differentiation trajectories.
Collapse
Affiliation(s)
- Gaurav Jumde
- Max Delbrück Center for Molecular Medicine, Berlin Institute for Medical Systems Biology, 10115 Berlin, Germany; Humboldt Universität zu Berlin, Faculty of Life Sciences, Department of Biology, 10115 Berlin, Germany
| | - Bastiaan Spanjaard
- Max Delbrück Center for Molecular Medicine, Berlin Institute for Medical Systems Biology, 10115 Berlin, Germany; Charité Universitätsmedizin Berlin, 10117 Berlin, Germany.
| | - Jan Philipp Junker
- Max Delbrück Center for Molecular Medicine, Berlin Institute for Medical Systems Biology, 10115 Berlin, Germany; Charité Universitätsmedizin Berlin, 10117 Berlin, Germany.
| |
Collapse
|
14
|
Sha Y, Qiu Y, Zhou P, Nie Q. Reconstructing growth and dynamic trajectories from single-cell transcriptomics data. NAT MACH INTELL 2023; 6:25-39. [PMID: 38274364 PMCID: PMC10805654 DOI: 10.1038/s42256-023-00763-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2023] [Accepted: 10/25/2023] [Indexed: 01/27/2024]
Abstract
Time-series single-cell RNA sequencing (scRNA-seq) datasets provide unprecedented opportunities to learn dynamic processes of cellular systems. Due to the destructive nature of sequencing, it remains challenging to link the scRNA-seq snapshots sampled at different time points. Here we present TIGON, a dynamic, unbalanced optimal transport algorithm that reconstructs dynamic trajectories and population growth simultaneously as well as the underlying gene regulatory network from multiple snapshots. To tackle the high-dimensional optimal transport problem, we introduce a deep learning method using a dimensionless formulation based on the Wasserstein-Fisher-Rao (WFR) distance. TIGON is evaluated on simulated data and compared with existing methods for its robustness and accuracy in predicting cell state transition and cell population growth. Using three scRNA-seq datasets, we show the importance of growth in the temporal inference, TIGON's capability in reconstructing gene expression at unmeasured time points and its applications to temporal gene regulatory networks and cell-cell communication inference.
Collapse
Affiliation(s)
- Yutong Sha
- Department of Mathematics, University of California, Irvine, Irvine, CA USA
| | - Yuchi Qiu
- Department of Mathematics, Michigan State University, East Lansing, MI USA
| | - Peijie Zhou
- Department of Mathematics, University of California, Irvine, Irvine, CA USA
| | - Qing Nie
- Department of Mathematics, University of California, Irvine, Irvine, CA USA
- Department of Developmental and Cell Biology, University of California, Irvine, Irvine, CA USA
- The NSF-Simons Center for Multiscale Cell Fate Research, University of California, Irvine, Irvine, CA USA
| |
Collapse
|
15
|
Deng H, Rukhlendo OS, Joshi D, Hu X, Junk P, Tuliakova A, Kholodenko BN, Schwartz MA. cSTAR analysis identifies endothelial cell cycle as a key regulator of flow-dependent artery remodeling. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.10.24.563764. [PMID: 37961694 PMCID: PMC10634797 DOI: 10.1101/2023.10.24.563764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/15/2023]
Abstract
Fluid shear stress (FSS) from blood flow is sensed by vascular endothelial cells (ECs) to determine vessel stability, remodeling and susceptibility to atherosclerosis and other inflammatory diseases but the regulatory networks that govern these behaviors are only partially understood. We used cSTAR, a powerful new computational method, to define EC transcriptomic states under low shear stress (LSS) that triggers vessel inward remodeling, physiological shear stress (PSS) that stabilizes vessels, high shear stress (HSS) that triggers outward remodeling, and oscillatory shear stress (OSS) that confers disease susceptibility, all in comparison to cells under static conditions (STAT). We combined these results with the LINCS database where EC transcriptomic responses to drug treatments to define a preliminary regulatory network in which the cyclin-dependent kinases CDK1/2 play a central role in promoting vessel stability. Experimental analysis showed that PSS induced a strong late G1 cell cycle arrest in which CDK2 was activated. EC deletion of CDK2 in mice resulted in inward artery remodeling and both pulmonary and systemic hypertension. These results validate use of cSTAR to determine EC state and in vivo vessel behavior, reveal unexpected features of EC phenotype under different FSS conditions, and identify CDK2 as a key element within the EC regulatory network that governs artery remodeling.
Collapse
|
16
|
Kana O, Nault R, Filipovic D, Marri D, Zacharewski T, Bhattacharya S. Generative modeling of single-cell gene expression for dose-dependent chemical perturbations. PATTERNS (NEW YORK, N.Y.) 2023; 4:100817. [PMID: 37602218 PMCID: PMC10436058 DOI: 10.1016/j.patter.2023.100817] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Revised: 12/07/2022] [Accepted: 07/14/2023] [Indexed: 08/22/2023]
Abstract
Single-cell sequencing reveals the heterogeneity of cellular response to chemical perturbations. However, testing all relevant combinations of cell types, chemicals, and doses is a daunting task. A deep generative learning formalism called variational autoencoders (VAEs) has been effective in predicting single-cell gene expression perturbations for single doses. Here, we introduce single-cell variational inference of dose-response (scVIDR), a VAE-based model that predicts both single-dose and multiple-dose cellular responses better than existing models. We show that scVIDR can predict dose-dependent gene expression across mouse hepatocytes, human blood cells, and cancer cell lines. We biologically interpret the latent space of scVIDR using a regression model and use scVIDR to order individual cells based on their sensitivity to chemical perturbation by assigning each cell a "pseudo-dose" value. We envision that scVIDR can help reduce the need for repeated animal testing across tissues, chemicals, and doses.
Collapse
Affiliation(s)
- Omar Kana
- Department of Pharmacology and Toxicology, Michigan State University, East Lansing, MI 48824, USA
- Institute for Integrative Toxicology, Michigan State University, East Lansing, MI 48824, USA
- Institute for Quantitative Health Science & Engineering, Michigan State University, East Lansing, MI 48824, USA
| | - Rance Nault
- Institute for Integrative Toxicology, Michigan State University, East Lansing, MI 48824, USA
- Department of Biochemistry and Molecular Biology Michigan State University, Michigan State University, East Lansing, MI 48824, USA
| | - David Filipovic
- Institute for Quantitative Health Science & Engineering, Michigan State University, East Lansing, MI 48824, USA
- Department of Biomedical Engineering, Michigan State University, East Lansing, MI 48824, USA
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, MI 48824, USA
| | - Daniel Marri
- Institute for Quantitative Health Science & Engineering, Michigan State University, East Lansing, MI 48824, USA
- Department of Biomedical Engineering, Michigan State University, East Lansing, MI 48824, USA
| | - Tim Zacharewski
- Institute for Integrative Toxicology, Michigan State University, East Lansing, MI 48824, USA
- Department of Biochemistry and Molecular Biology Michigan State University, Michigan State University, East Lansing, MI 48824, USA
| | - Sudin Bhattacharya
- Department of Pharmacology and Toxicology, Michigan State University, East Lansing, MI 48824, USA
- Institute for Integrative Toxicology, Michigan State University, East Lansing, MI 48824, USA
- Institute for Quantitative Health Science & Engineering, Michigan State University, East Lansing, MI 48824, USA
- Department of Biomedical Engineering, Michigan State University, East Lansing, MI 48824, USA
| |
Collapse
|
17
|
Kernfeld E, Yang Y, Weinstock JS, Battle A, Cahan P. A systematic comparison of computational methods for expression forecasting. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.28.551039. [PMID: 37577640 PMCID: PMC10418073 DOI: 10.1101/2023.07.28.551039] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/15/2023]
Abstract
Due to the abundance of single cell RNA-seq data, a number of methods for predicting expression after perturbation have recently been published. Expression prediction methods are enticing because they promise to answer pressing questions in fields ranging from developmental genetics to cell fate engineering and because they are faster, cheaper, and higher-throughput than their experimental counterparts. However, the absolute and relative accuracy of these methods is poorly characterized, limiting their informed use, their improvement, and the interpretation of their predictions. To address these issues, we created a benchmarking platform that combines a panel of large-scale perturbation datasets with an expression forecasting software engine that encompasses or interfaces to current methods. We used our platform to systematically assess methods, parameters, and sources of auxiliary data. We found that uninformed baseline predictions, which were not always included in prior evaluations, yielded the same or better mean absolute error than benchmarked methods in all test cases. These results cast doubt on the ability of current expression forecasting methods to provide mechanistic insights or to rank hypotheses for experimental follow-up. However, given the rapid pace of innovation in the field, new approaches may yield more accurate expression predictions. Our platform will serve as a neutral benchmark to improve methods and to identify contexts in which expression prediction can succeed.
Collapse
|
18
|
Hossain I, Fanfani V, Quackenbush J, Burkholz R. Biologically informed NeuralODEs for genome-wide regulatory dynamics. RESEARCH SQUARE 2023:rs.3.rs-2675584. [PMID: 36993392 PMCID: PMC10055646 DOI: 10.21203/rs.3.rs-2675584/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Models that are formulated as ordinary differential equations (ODEs) can accurately explain temporal gene expression patterns and promise to yield new insights into important cellular processes, disease progression, and intervention design. Learning such ODEs is challenging, since we want to predict the evolution of gene expression in a way that accurately encodes the causal gene-regulatory network (GRN) governing the dynamics and the nonlinear functional relationships between genes. Most widely used ODE estimation methods either impose too many parametric restrictions or are not guided by meaningful biological insights, both of which impedes scalability and/or explainability. To overcome these limitations, we developed PHOENIX, a modeling framework based on neural ordinary differential equations (NeuralODEs) and Hill-Langmuir kinetics, that can flexibly incorporate prior domain knowledge and biological constraints to promote sparse, biologically interpretable representations of ODEs. We test accuracy of PHOENIX in a series of in silico experiments benchmarking it against several currently used tools for ODE estimation. We also demonstrate PHOENIX's flexibility by studying oscillating expression data from synchronized yeast cells and assess its scalability by modelling genome-scale breast cancer expression for samples ordered in pseudotime. Finally, we show how the combination of user-defined prior knowledge and functional forms from systems biology allows PHOENIX to encode key properties of the underlying GRN, and subsequently predict expression patterns in a biologically explainable way.
Collapse
Affiliation(s)
- Intekhab Hossain
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Viola Fanfani
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Rebekka Burkholz
- Helmholtz Center for Information Security (CISPA), Saarbrücken, Germany
| |
Collapse
|
19
|
Allesøe RL, Lundgaard AT, Hernández Medina R, Aguayo-Orozco A, Johansen J, Nissen JN, Brorsson C, Mazzoni G, Niu L, Biel JH, Brasas V, Webel H, Benros ME, Pedersen AG, Chmura PJ, Jacobsen UP, Mari A, Koivula R, Mahajan A, Vinuela A, Tajes JF, Sharma S, Haid M, Hong MG, Musholt PB, De Masi F, Vogt J, Pedersen HK, Gudmundsdottir V, Jones A, Kennedy G, Bell J, Thomas EL, Frost G, Thomsen H, Hansen E, Hansen TH, Vestergaard H, Muilwijk M, Blom MT, 't Hart LM, Pattou F, Raverdy V, Brage S, Kokkola T, Heggie A, McEvoy D, Mourby M, Kaye J, Hattersley A, McDonald T, Ridderstråle M, Walker M, Forgie I, Giordano GN, Pavo I, Ruetten H, Pedersen O, Hansen T, Dermitzakis E, Franks PW, Schwenk JM, Adamski J, McCarthy MI, Pearson E, Banasik K, Rasmussen S, Brunak S, Thomas CE, Haussler R, Beulens J, Rutters F, Nijpels G, van Oort S, Groeneveld L, Elders P, Giorgino T, Rodriquez M, Nice R, Perry M, Bianzano S, Graefe-Mody U, Hennige A, Grempler R, Baum P, Stærfeldt HH, Shah N, Teare H, Ehrhardt B, Tillner J, Dings C, Lehr T, Scherer N, Sihinevich I, Cabrelli L, Loftus H, Bizzotto R, Tura A, Dekkers K, van Leeuwen N, Groop L, Slieker R, Ramisch A, Jennison C, McVittie I, Frau F, Steckel-Hamann B, Adragni K, Thomas M, Pasdar NA, Fitipaldi H, Kurbasic A, Mutie P, Pomares-Millan H, Bonnefond A, Canouil M, Caiazzo R, Verkindt H, Holl R, Kuulasmaa T, Deshmukh H, Cederberg H, Laakso M, Vangipurapu J, Dale M, Thorand B, Nicolay C, Fritsche A, Hill A, Hudson M, Thorne C, Allin K, Arumugam M, Jonsson A, Engelbrechtsen L, Forman A, Dutta A, Sondertoft N, Fan Y, Gough S, Robertson N, McRobert N, Wesolowska-Andersen A, Brown A, Davtian D, Dawed A, Donnelly L, Palmer C, White M, Ferrer J, Whitcher B, Artati A, Prehn C, Adam J, Grallert H, Gupta R, Sackett PW, Nilsson B, Tsirigos K, Eriksen R, Jablonka B, Uhlen M, Gassenhuber J, Baltauss T, de Preville N, Klintenberg M, Abdalla M. Discovery of drug-omics associations in type 2 diabetes with generative deep-learning models. Nat Biotechnol 2023; 41:399-408. [PMID: 36593394 PMCID: PMC10017515 DOI: 10.1038/s41587-022-01520-x] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2022] [Accepted: 09/20/2022] [Indexed: 01/03/2023]
Abstract
The application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale and heterogeneous nature of multi-modal data makes integration and inference a non-trivial task. We developed a deep-learning-based framework, multi-omics variational autoencoders (MOVE), to integrate such data and applied it to a cohort of 789 people with newly diagnosed type 2 diabetes with deep multi-omics phenotyping from the DIRECT consortium. Using in silico perturbations, we identified drug-omics associations across the multi-modal datasets for the 20 most prevalent drugs given to people with type 2 diabetes with substantially higher sensitivity than univariate statistical tests. From these, we among others, identified novel associations between metformin and the gut microbiota as well as opposite molecular responses for the two statins, simvastatin and atorvastatin. We used the associations to quantify drug-drug similarities, assess the degree of polypharmacy and conclude that drug effects are distributed across the multi-omics modalities.
Collapse
Affiliation(s)
- Rosa Lundbye Allesøe
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark.,Copenhagen Research Centre for Mental Health, Mental Health Centre Copenhagen, Copenhagen University Hospital, Copenhagen, Denmark
| | - Agnete Troen Lundgaard
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Ricardo Hernández Medina
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Alejandro Aguayo-Orozco
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Joachim Johansen
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Jakob Nybo Nissen
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Caroline Brorsson
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Gianluca Mazzoni
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Lili Niu
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Jorge Hernansanz Biel
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Valentas Brasas
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Henry Webel
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Michael Eriksen Benros
- Copenhagen Research Centre for Mental Health, Mental Health Centre Copenhagen, Copenhagen University Hospital, Copenhagen, Denmark.,Department of Immunology and Microbiology, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Anders Gorm Pedersen
- Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Piotr Jaroslaw Chmura
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Ulrik Plesner Jacobsen
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Andrea Mari
- C.N.R. Institute of Neuroscience, Padova, Italy
| | - Robert Koivula
- Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK
| | - Anubha Mahajan
- Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK
| | - Ana Vinuela
- Department of Genetic Medicine and Development, University of Geneva Medical School, Geneva, Switzerland.,Biosciences Institute, Faculty of Medical Sciences, Newcastle University, Newcastle, UK
| | | | - Sapna Sharma
- Research Unit of Molecular Epidemiology, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Bavaria, Germany.,Institute of Epidemiology, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Bavaria, Germany.,Chair of Food Chemistry and Molecular and Sensory Science, Technical University of Munich, Freising, Germany
| | - Mark Haid
- Metabolomics and Proteomics Core, Helmholtz Zentrum Muenchen, German Research Center for Environmental Health, Neuherberg, Germany
| | - Mun-Gwan Hong
- Affinity Proteomics, Science for Life Laboratory, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, Solna, Sweden
| | - Petra B Musholt
- Research and Development Global Development, Translational Medicine and Clinical Pharmacology, Sanofi-Aventis Deutschland, Frankfurt, Germany
| | - Federico De Masi
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Josef Vogt
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Helle Krogh Pedersen
- Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark.,Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Valborg Gudmundsdottir
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Angus Jones
- University of Exeter Medical School, Exeter, UK
| | - Gwen Kennedy
- The Immunoassay Biomarker Core Laboratory, School of Medicine, University of Dundee, Dundee, UK
| | - Jimmy Bell
- Research Centre for Optimal Health, Department of Life Sciences, University of Westminster, London, UK
| | - E Louise Thomas
- Research Centre for Optimal Health, Department of Life Sciences, University of Westminster, London, UK
| | - Gary Frost
- Section for Nutrition Research, Faculty of Medicine, Imperial College London, London, UK
| | - Henrik Thomsen
- Department of Radiology, Copenhagen University Hospital Herlev-Gentofte, Herlev, Denmark
| | - Elizaveta Hansen
- Department of Radiology, Copenhagen University Hospital Herlev-Gentofte, Herlev, Denmark
| | - Tue Haldor Hansen
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Henrik Vestergaard
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Mirthe Muilwijk
- Department of Epidemiology and Data Science, Amsterdam Public Health Research Institute, Amsterdam UMC, Vrije Universiteit Amsterdam, Amsterdam, the Netherlands
| | - Marieke T Blom
- Department of General Practice, Amsterdam Public Health Research Institute, Amsterdam UMC, Vrije Universiteit Amsterdam, Amsterdam, the Netherlands
| | - Leen M 't Hart
- Department of Epidemiology and Data Science, Amsterdam Public Health Research Institute, Amsterdam UMC, Vrije Universiteit Amsterdam, Amsterdam, the Netherlands.,Department of Biomedical Data Science, Section Molecular Epidemiology, Leiden University Medical Center, Leiden, the Netherlands.,Department of Cell and Chemical Biology, Leiden University Medical Center, Leiden, the Netherlands
| | - Francois Pattou
- Inserm, Univ Lille, CHU Lille, Lille Pasteur Institute, EGID, Lille, France
| | - Violeta Raverdy
- Inserm, Univ Lille, CHU Lille, Lille Pasteur Institute, EGID, Lille, France
| | - Soren Brage
- MRC Epidemiology Unit, University of Cambridge School of Clinical Medicine, Cambridge, UK
| | - Tarja Kokkola
- Department of Medicine, University of Eastern Finland, Kuopio, Finland
| | - Alison Heggie
- Institute of Cellular Medicine, Newcastle University, Newcastle, UK
| | - Donna McEvoy
- Diabetes Research Network, Royal Victoria Infirmary, Newcastle, UK
| | - Miranda Mourby
- Centre for Health, Law and Emerging Technologies (HeLEX), Faculty of Law, University of Oxford, Oxford, UK
| | - Jane Kaye
- Centre for Health, Law and Emerging Technologies (HeLEX), Faculty of Law, University of Oxford, Oxford, UK
| | | | | | - Martin Ridderstråle
- Lund University Diabetes Centre, Department of Clinical Sciences, Lund University, Malmö, Sweden
| | - Mark Walker
- Translational and Clinical Research Institute, Faculty of Medical Sciences, Newcastle University, Newcastle, UK
| | - Ian Forgie
- Division of Population Health & Genomics, School of Medicine, University of Dundee, Dundee, UK
| | - Giuseppe N Giordano
- Genetic and Molecular Epidemiology Unit, Lund University Diabetes Centre, Department of Clinical Sciences, CRC, Lund University, SUS, Malmö, Sweden
| | - Imre Pavo
- Eli Lilly Regional Operations, Vienna, Austria
| | - Hartmut Ruetten
- Research and Development Global Development, Translational Medicine and Clinical Pharmacology, Sanofi-Aventis Deutschland, Frankfurt, Germany
| | - Oluf Pedersen
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Torben Hansen
- Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
| | - Emmanouil Dermitzakis
- Department of Genetic Medicine and Development, University of Geneva Medical School, Geneva, Switzerland
| | - Paul W Franks
- Lund University Diabetes Centre, Department of Clinical Sciences, Lund University, Malmö, Sweden.,Harvard T.H. Chan School of Public Health, Boston, MA, USA.,OCDEM, Radcliffe Department of Medicine, University of Oxford, Oxford, UK
| | - Jochen M Schwenk
- Affinity Proteomics, Science for Life Laboratory, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, Solna, Sweden
| | - Jerzy Adamski
- Institute of Experimental Genetics, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany.,Department of Biochemistry, Yong Loo Lin School of Medicine, National University of Singapore, Singapore, Singapore.,Institute of Biochemistry, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia
| | - Mark I McCarthy
- Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK.,Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Oxford, UK.,Genentech, South San Francisco, CA, USA
| | - Ewan Pearson
- Division of Population Health & Genomics, School of Medicine, University of Dundee, Dundee, UK
| | - Karina Banasik
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Simon Rasmussen
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark.
| | - Søren Brunak
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark. .,Department of Health Technology, Technical University of Denmark, Kongens Lyngby, Denmark.
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
20
|
Ouyang JF, Chothani S, Rackham OJ. Deep learning models will shape the future of stem cell research. Stem Cell Reports 2023; 18:6-12. [PMID: 36630908 PMCID: PMC9860061 DOI: 10.1016/j.stemcr.2022.11.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2022] [Revised: 11/10/2022] [Accepted: 11/11/2022] [Indexed: 01/12/2023] Open
Abstract
Our ability to understand and control stem cell biology is being augmented by developments on two fronts, our ability to collect more data describing cell state and our capability to comprehend these data using deep learning models. Here we consider the impact deep learning will have in the future of stem cell research. We explore the importance of generating data suitable for these methods, the requirement for close collaboration between experimental and computational researchers, and the challenges we face to do this fairly and effectively. Achieving this will ensure that the resulting deep learning models are biologically meaningful and computationally tractable.
Collapse
Affiliation(s)
- John F. Ouyang
- Duke-NUS Medical School, Program in Cardiovascular and Metabolic Disorders (CVMD) and Centre for Computational Biology (CCB), Singapore, Singapore
| | - Sonia Chothani
- Duke-NUS Medical School, Program in Cardiovascular and Metabolic Disorders (CVMD) and Centre for Computational Biology (CCB), Singapore, Singapore
| | - Owen J.L. Rackham
- Duke-NUS Medical School, Program in Cardiovascular and Metabolic Disorders (CVMD) and Centre for Computational Biology (CCB), Singapore, Singapore,School of Biological Sciences, University of Southampton, Southampton, UK,The Alan Turing Institute, The British Library, London, UK,Corresponding author
| |
Collapse
|
21
|
Xing J. Reconstructing data-driven governing equations for cell phenotypic transitions: integration of data science and systems biology. Phys Biol 2022; 19:10.1088/1478-3975/ac8c16. [PMID: 35998617 PMCID: PMC9585661 DOI: 10.1088/1478-3975/ac8c16] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 08/23/2022] [Indexed: 11/11/2022]
Abstract
Cells with the same genome can exist in different phenotypes and can change between distinct phenotypes when subject to specific stimuli and microenvironments. Some examples include cell differentiation during development, reprogramming for induced pluripotent stem cells and transdifferentiation, cancer metastasis and fibrosis progression. The regulation and dynamics of cell phenotypic conversion is a fundamental problem in biology, and has a long history of being studied within the formalism of dynamical systems. A main challenge for mechanism-driven modeling studies is acquiring sufficient amount of quantitative information for constraining model parameters. Advances in quantitative experimental approaches, especially high throughput single-cell techniques, have accelerated the emergence of a new direction for reconstructing the governing dynamical equations of a cellular system from quantitative single-cell data, beyond the dominant statistical approaches. Here I review a selected number of recent studies using live- and fixed-cell data and provide my perspective on future development.
Collapse
Affiliation(s)
- Jianhua Xing
- Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, PA 15232, USA
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15232, USA
- UPMC-Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA, USA
| |
Collapse
|
22
|
Xing J. Reconstructing data-driven governing equations for cell phenotypic transitions: integration of data science and systems biology. Phys Biol 2022. [PMID: 35998617 DOI: 10.48550/arxiv.2203.14964] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
Abstract
Cells with the same genome can exist in different phenotypes and can change between distinct phenotypes when subject to specific stimuli and microenvironments. Some examples include cell differentiation during development, reprogramming for induced pluripotent stem cells and transdifferentiation, cancer metastasis and fibrosis progression. The regulation and dynamics of cell phenotypic conversion is a fundamental problem in biology, and has a long history of being studied within the formalism of dynamical systems. A main challenge for mechanism-driven modeling studies is acquiring sufficient amount of quantitative information for constraining model parameters. Advances in quantitative experimental approaches, especially high throughput single-cell techniques, have accelerated the emergence of a new direction for reconstructing the governing dynamical equations of a cellular system from quantitative single-cell data, beyond the dominant statistical approaches. Here I review a selected number of recent studies using live- and fixed-cell data and provide my perspective on future development.
Collapse
Affiliation(s)
- Jianhua Xing
- Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, PA 15232, United States of America.,Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15232, United States of America.,UPMC-Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA, United States of America
| |
Collapse
|
23
|
Cho H, Kuo YH, Rockne RC. Comparison of cell state models derived from single-cell RNA sequencing data: graph versus multi-dimensional space. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2022; 19:8505-8536. [PMID: 35801475 PMCID: PMC9308174 DOI: 10.3934/mbe.2022395] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Single-cell sequencing technologies have revolutionized molecular and cellular biology and stimulated the development of computational tools to analyze the data generated from these technology platforms. However, despite the recent explosion of computational analysis tools, relatively few mathematical models have been developed to utilize these data. Here we compare and contrast two cell state geometries for building mathematical models of cell state-transitions with single-cell RNA-sequencing data with hematopoeisis as a model system; (i) by using partial differential equations on a graph representing intermediate cell states between known cell types, and (ii) by using the equations on a multi-dimensional continuous cell state-space. As an application of our approach, we demonstrate how the calibrated models may be used to mathematically perturb normal hematopoeisis to simulate, predict, and study the emergence of novel cell states during the pathogenesis of acute myeloid leukemia. We particularly focus on comparing the strength and weakness of the graph model and multi-dimensional model.
Collapse
Affiliation(s)
- Heyrim Cho
- Department of Mathematics, University of California Riverside, Riverside, CA, USA
- Interdisciplinary Center for Quantitative Modeling in Biology, University of California Riverside, Riverside, CA, USA
| | - Ya-Huei Kuo
- Department of Hematologic Malignancies Translational Science, City of Hope, Duarte, CA, USA
| | - Russell C. Rockne
- Department of Computational and Quantitative Medicine, Division of Mathematical Oncology, City of Hope, Duarte, CA, USA
- Interdisciplinary Center for Quantitative Modeling in Biology, University of California Riverside, Riverside, CA, USA
| |
Collapse
|
24
|
Fang D, Tan XH, Song WP, Gu YY, Pan JC, Yang XQ, Song WD, Yuan YM, Peng J, Zhang ZC, Xin ZC, Li XS, Guan RL. Single-Cell RNA Sequencing of Human Corpus Cavernosum Reveals Cellular Heterogeneity Landscapes in Erectile Dysfunction. Front Endocrinol (Lausanne) 2022; 13:874915. [PMID: 35518933 PMCID: PMC9066803 DOI: 10.3389/fendo.2022.874915] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/13/2022] [Accepted: 03/16/2022] [Indexed: 11/26/2022] Open
Abstract
PURPOSE To assess the diverse cell populations of human corpus cavernosum in patients with severe erectile dysfunction (ED) at the single-cell level. METHODS Penile tissues collected from three patients were subjected to single-cell RNA sequencing using the BD Rhapsody™ platform. Common bioinformatics tools were used to analyze cellular heterogeneity and gene expression profiles from generated raw data, including the packages Seurat, Monocle, and CellPhoneDB. RESULTS Disease-related heterogeneity of cell types was determined in the cavernous tissue such as endothelial cells (ECs), smooth muscle cells, fibroblasts, and immune cells. Reclustering analysis of ECs identified an arteriole ECs subcluster and another one with gene signatures of fibroblasts. The proportion of fibroblasts was higher than the other cell populations and had the most significant cellular heterogeneity, in which a distinct subcluster co-expressed endothelial markers. The transition trajectory of differentiation from smooth muscle cells into fibroblasts was depicted using the pseudotime analysis, suggesting that the expansion of corpus cavernosum is possibly compromised as a result of fibrosis. Cell-cell communications among ECs, smooth muscle cells, fibroblasts, and macrophages were robust, which indicated that inflammation may also have a crucial role in the development of ED. CONCLUSIONS Our study has demonstrated a comprehensive single-cell atlas of cellular components in human corpus cavernosum of ED, providing in-depth insights into the pathogenesis. Future research is warranted to explore disease-specific alterations for individualized treatment of ED.
Collapse
Affiliation(s)
- Dong Fang
- Department of Urology, Peking University First Hospital, Beijing, China
- Institute of Urology, Peking University, Beijing, China
- Beijing Key Laboratory of Urogenital Diseases (male) Molecular Diagnosis and Treatment Center, Beijing, China
| | - Xiao-Hui Tan
- Department of Urology, Peking University First Hospital, Beijing, China
- Institute of Urology, Peking University, Beijing, China
- Beijing Key Laboratory of Urogenital Diseases (male) Molecular Diagnosis and Treatment Center, Beijing, China
| | - Wen-Peng Song
- Institute of Urology, Peking University, Beijing, China
- Beijing Key Laboratory of Urogenital Diseases (male) Molecular Diagnosis and Treatment Center, Beijing, China
- Department of Dental Implant Center, Beijing Stomatological Hospital, School of Stomatology, Capital Medical University, Beijing, China
| | - Yang-Yang Gu
- Institute of Urology, Peking University, Beijing, China
- Beijing Key Laboratory of Urogenital Diseases (male) Molecular Diagnosis and Treatment Center, Beijing, China
- Department of Radiation Medicine, Institute of Systems Biomedicine, School of Basic Medical Sciences, Peking University Health Science Center, Beijing, China
| | - Jian-Cheng Pan
- Male Reproductive and Sexual Medicine, Department of Urology, The Second Hospital of Tianjin Medical University, Tianjin, China
- Institute of Urology, Tianjin Medical University, Tianjin, China
| | - Xiao-Qing Yang
- Male Reproductive and Sexual Medicine, Department of Urology, The Second Hospital of Tianjin Medical University, Tianjin, China
- Institute of Urology, Tianjin Medical University, Tianjin, China
| | - Wei-Dong Song
- Department of Urology, Peking University First Hospital, Beijing, China
- Institute of Urology, Peking University, Beijing, China
- Beijing Key Laboratory of Urogenital Diseases (male) Molecular Diagnosis and Treatment Center, Beijing, China
| | - Yi-Ming Yuan
- Department of Urology, Peking University First Hospital, Beijing, China
- Institute of Urology, Peking University, Beijing, China
- Beijing Key Laboratory of Urogenital Diseases (male) Molecular Diagnosis and Treatment Center, Beijing, China
| | - Jing Peng
- Department of Urology, Peking University First Hospital, Beijing, China
- Institute of Urology, Peking University, Beijing, China
- Beijing Key Laboratory of Urogenital Diseases (male) Molecular Diagnosis and Treatment Center, Beijing, China
| | - Zhi-Chao Zhang
- Department of Urology, Peking University First Hospital, Beijing, China
- Institute of Urology, Peking University, Beijing, China
- Beijing Key Laboratory of Urogenital Diseases (male) Molecular Diagnosis and Treatment Center, Beijing, China
| | - Zhong-Cheng Xin
- Male Reproductive and Sexual Medicine, Department of Urology, The Second Hospital of Tianjin Medical University, Tianjin, China
- Institute of Urology, Tianjin Medical University, Tianjin, China
| | - Xue-Song Li
- Department of Urology, Peking University First Hospital, Beijing, China
- Institute of Urology, Peking University, Beijing, China
- Beijing Key Laboratory of Urogenital Diseases (male) Molecular Diagnosis and Treatment Center, Beijing, China
- *Correspondence: Rui-Li Guan, ; Xue-Song Li,
| | - Rui-Li Guan
- Department of Urology, Peking University First Hospital, Beijing, China
- Institute of Urology, Peking University, Beijing, China
- Beijing Key Laboratory of Urogenital Diseases (male) Molecular Diagnosis and Treatment Center, Beijing, China
- *Correspondence: Rui-Li Guan, ; Xue-Song Li,
| |
Collapse
|
25
|
Jiang Q, Zhang S, Wan L. Dynamic inference of cell developmental complex energy landscape from time series single-cell transcriptomic data. PLoS Comput Biol 2022; 18:e1009821. [PMID: 35073331 PMCID: PMC8812873 DOI: 10.1371/journal.pcbi.1009821] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2021] [Revised: 02/03/2022] [Accepted: 01/10/2022] [Indexed: 12/27/2022] Open
Abstract
Time series single-cell RNA sequencing (scRNA-seq) data are emerging. However, dynamic inference of an evolving cell population from time series scRNA-seq data is challenging owing to the stochasticity and nonlinearity of the underlying biological processes. This calls for the development of mathematical models and methods capable of reconstructing cellular dynamic transition processes and uncovering the nonlinear cell-cell interactions. In this study, we present GraphFP, a nonlinear Fokker-Planck equation on graph based model and dynamic inference framework, with the aim of reconstructing the cell state-transition complex potential energy landscape from time series single-cell transcriptomic data. The free energy of our model explicitly takes into account of the cell-cell interactions in a nonlinear quadratic term. We then recast the model inference problem in the form of a dynamic optimal transport framework and solve it efficiently with the adjoint method of optimal control. We evaluated GraphFP on the time series scRNA-seq data set of embryonic murine cerebral cortex development. We illustrated that it 1) reconstructs cell state potential energy, which is a measure of cellular differentiation potency, 2) faithfully charts the probability flows between paired cell states over the dynamic processes of cell differentiation, and 3) accurately quantifies the stochastic dynamics of cell type frequencies on probability simplex in continuous time. We also illustrated that GraphFP is robust in terms of cluster labelling with different resolutions, as well as parameter choices. Meanwhile, GraphFP provides a model-based approach to delineate the cell-cell interactions that drive cell differentiation. GraphFP software is available at https://github.com/QiJiang-QJ/GraphFP. Dynamic inference of cell development processes from time series scRNA-seq data is a major challenge. Here, we present GraphFP, a coherent computational framework that simultaneously reconstructs the cell state-transition complex potential energy landscape and infers cell-cell interactions from time series single-cell transcriptomic data. Based on the mathematical framework of nonlinear Fokker-Planck equation on graph, GraphFP models the stochastic dynamics of the cell state/type frequencies on probability simplex in continuous time, where the free energy with a nonlinear quadratic interaction term is employed to characterize cell-cell interactions. We formulate the model inference problem in the form of a dynamic optimal transport framework and solve it efficiently with the celebrated adjoint method. GraphFP allows for 1) reconstructing cell state potential energy, which is a measure of cellular differentiation potency, 2) charting the probability flows between paired cell states over dynamic processes, 3) quantifying the stochastic dynamics of cell type frequencies on probability simplex in continuous time, and 4) delineating cell-cell interactions that drive cell differentiation. We show how GraphFP can be used to faithfully reveal and accurately quantify the cell development processes using the embryonic murine cerebral cortex development time series scRNA-seq dataset.
Collapse
Affiliation(s)
- Qi Jiang
- NCMIS, LSC, LSEC, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing, China
- School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing, China
| | - Shuo Zhang
- NCMIS, LSC, LSEC, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing, China
- School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing, China
| | - Lin Wan
- NCMIS, LSC, LSEC, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing, China
- School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing, China
- * E-mail:
| |
Collapse
|
26
|
Xu Y, Liu X, Cao X, Huang C, Liu E, Qian S, Liu X, Wu Y, Dong F, Qiu CW, Qiu J, Hua K, Su W, Wu J, Xu H, Han Y, Fu C, Yin Z, Liu M, Roepman R, Dietmann S, Virta M, Kengara F, Zhang Z, Zhang L, Zhao T, Dai J, Yang J, Lan L, Luo M, Liu Z, An T, Zhang B, He X, Cong S, Liu X, Zhang W, Lewis JP, Tiedje JM, Wang Q, An Z, Wang F, Zhang L, Huang T, Lu C, Cai Z, Wang F, Zhang J. Artificial intelligence: A powerful paradigm for scientific research. Innovation (N Y) 2021; 2:100179. [PMID: 34877560 PMCID: PMC8633405 DOI: 10.1016/j.xinn.2021.100179] [Citation(s) in RCA: 78] [Impact Index Per Article: 26.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2021] [Accepted: 10/26/2021] [Indexed: 12/18/2022] Open
Abstract
Artificial intelligence (AI) coupled with promising machine learning (ML) techniques well known from computer science is broadly affecting many aspects of various fields including science and technology, industry, and even our day-to-day life. The ML techniques have been developed to analyze high-throughput data with a view to obtaining useful insights, categorizing, predicting, and making evidence-based decisions in novel ways, which will promote the growth of novel applications and fuel the sustainable booming of AI. This paper undertakes a comprehensive survey on the development and application of AI in different aspects of fundamental sciences, including information science, mathematics, medical science, materials science, geoscience, life science, physics, and chemistry. The challenges that each discipline of science meets, and the potentials of AI techniques to handle these challenges, are discussed in detail. Moreover, we shed light on new research trends entailing the integration of AI into each scientific discipline. The aim of this paper is to provide a broad research guideline on fundamental sciences with potential infusion of AI, to help motivate researchers to deeply understand the state-of-the-art applications of AI-based fundamental sciences, and thereby to help promote the continuous development of these fundamental sciences.
Collapse
Affiliation(s)
- Yongjun Xu
- Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Xin Liu
- Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Xin Cao
- Zhongshan Hospital Institute of Clinical Science, Fudan University, Shanghai 200032, China
| | - Changping Huang
- Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Enke Liu
- Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Songshan Lake Materials Laboratory, Dongguan, Guangdong 523808, China
| | - Sen Qian
- Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China
| | - Xingchen Liu
- Institute of Coal Chemistry, Chinese Academy of Sciences, Taiyuan 030001, China
| | - Yanjun Wu
- Institute of Software, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Fengliang Dong
- National Center for Nanoscience and Technology, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Cheng-Wei Qiu
- Department of Electrical and Computer Engineering, National University of Singapore, Singapore 117583, Singapore
| | - Junjun Qiu
- Department of Gynaecology, Obstetrics and Gynaecology Hospital, Fudan University, Shanghai 200011, China
- Shanghai Key Laboratory of Female Reproductive Endocrine-Related Diseases, Shanghai 200011, China
| | - Keqin Hua
- Department of Gynaecology, Obstetrics and Gynaecology Hospital, Fudan University, Shanghai 200011, China
- Shanghai Key Laboratory of Female Reproductive Endocrine-Related Diseases, Shanghai 200011, China
| | - Wentao Su
- School of Food Science and Technology, Dalian Polytechnic University, Dalian 116034, China
| | - Jian Wu
- Second Affiliated Hospital School of Medicine, and School of Public Health, Zhejiang University, Hangzhou 310058, China
| | - Huiyu Xu
- Department of Obstetrics and Gynecology, Peking University Third Hospital, Beijing 100191, China
| | - Yong Han
- Zhejiang Provincial People’s Hospital, Hangzhou 310014, China
| | - Chenguang Fu
- School of Materials Science and Engineering, Zhejiang University, Hangzhou 310027, China
| | - Zhigang Yin
- Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, China
| | - Miao Liu
- Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Songshan Lake Materials Laboratory, Dongguan, Guangdong 523808, China
| | - Ronald Roepman
- Medical Center, Radboud University, 6500 Nijmegen, the Netherlands
| | - Sabine Dietmann
- Institute for Informatics, Washington University School of Medicine, St. Louis, MO 63110, USA
| | - Marko Virta
- Department of Microbiology, University of Helsinki, 00014 Helsinki, Finland
| | - Fredrick Kengara
- School of Pure and Applied Sciences, Bomet University College, Bomet 20400, Kenya
| | - Ze Zhang
- Agriculture College of Shihezi University, Xinjiang 832000, China
| | - Lifu Zhang
- Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
- Agriculture College of Shihezi University, Xinjiang 832000, China
| | - Taolan Zhao
- Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
| | - Ji Dai
- The Brain Cognition and Brain Disease Institute, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
- University of Chinese Academy of Sciences, Beijing 100049, China
- Shenzhen-Hong Kong Institute of Brain Science-Shenzhen Fundamental Research Institutions, Shenzhen 518055, China
| | | | - Liang Lan
- Department of Communication Studies, Hong Kong Baptist University, Hong Kong, China
| | - Ming Luo
- South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China
- Center of Economic Botany, Core Botanical Gardens, Chinese Academy of Sciences, Guangzhou 510650, China
| | - Zhaofeng Liu
- Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Tao An
- Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China
| | - Bin Zhang
- Institute of Coal Chemistry, Chinese Academy of Sciences, Taiyuan 030001, China
| | - Xiao He
- Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China
| | - Shan Cong
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, China
| | - Xiaohong Liu
- Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences, Chongqing 400714, China
| | - Wei Zhang
- Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences, Chongqing 400714, China
| | - James P. Lewis
- Institute of Coal Chemistry, Chinese Academy of Sciences, Taiyuan 030001, China
| | - James M. Tiedje
- Center for Microbial Ecology, Department of Plant, Soil and Microbial Sciences, Michigan State University, East Lansing, MI 48824, USA
| | - Qi Wang
- Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
- Zhejiang Lab, Hangzhou 311121, China
| | - Zhulin An
- Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Fei Wang
- Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Libo Zhang
- Institute of Software, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Tao Huang
- Shanghai Institute of Nutrition and Health, Chinese Academy of Sciences, Shanghai 200031, China
| | - Chuan Lu
- Department of Computer Science, Aberystwyth University, Aberystwyth, Ceredigion SY23 3FL, UK
| | - Zhipeng Cai
- Department of Computer Science, Georgia State University, Atlanta, GA 30303, USA
| | - Fang Wang
- Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Jiabao Zhang
- Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|