1
|
Copperman J, Mclean IC, Gross SM, Singh J, Chang YH, Zuckerman DM, Heiser LM. Single-cell morphodynamical trajectories enable prediction of gene expression accompanying cell state change. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.18.576248. [PMID: 38293173 PMCID: PMC10827140 DOI: 10.1101/2024.01.18.576248] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/01/2024]
Abstract
Extracellular signals induce changes to molecular programs that modulate multiple cellular phenotypes, including proliferation, motility, and differentiation status. The connection between dynamically adapting phenotypic states and the molecular programs that define them is not well understood. Here we develop data-driven models of single-cell phenotypic responses to extracellular stimuli by linking gene transcription levels to "morphodynamics" - changes in cell morphology and motility observable in time-lapse image data. We adopt a dynamics-first view of cell state by grouping single-cell trajectories into states with shared morphodynamic responses. The single-cell trajectories enable development of a first-of-its-kind computational approach to map live-cell dynamics to snapshot gene transcript levels, which we term MMIST, Molecular and Morphodynamics-Integrated Single-cell Trajectories. The key conceptual advance of MMIST is that cell behavior can be quantified based on dynamically defined states and that extracellular signals alter the overall distribution of cell states by altering rates of switching between states. We find a cell state landscape that is bound by epithelial and mesenchymal endpoints, with distinct sequences of epithelial to mesenchymal transition (EMT) and mesenchymal to epithelial transition (MET) intermediates. The analysis yields predictions for gene expression changes consistent with curated EMT gene sets and provides a prediction of thousands of RNA transcripts through extracellular signal-induced EMT and MET with near-continuous time resolution. The MMIST framework leverages true single-cell dynamical behavior to generate molecular-level omics inferences and is broadly applicable to other biological domains, time-lapse imaging approaches and molecular snapshot data.
Collapse
Affiliation(s)
- Jeremy Copperman
- Cancer Early Detection Advanced Research Center, Oregon Health and Science University, Portland OR 97239, U.S.A
| | - Ian C. Mclean
- Department of Biomedical Engineering, Oregon Health and Science University, Portland OR 97239, U.S.A
| | | | - Jalim Singh
- Knight Cancer Institute, Oregon Health and Science University, Portland OR 97239, U.S.A
| | - Young Hwan Chang
- Department of Biomedical Engineering, Oregon Health and Science University, Portland OR 97239, U.S.A
- Knight Cancer Institute, Oregon Health and Science University, Portland OR 97239, U.S.A
| | - Daniel M. Zuckerman
- Department of Biomedical Engineering, Oregon Health and Science University, Portland OR 97239, U.S.A
- Knight Cancer Institute, Oregon Health and Science University, Portland OR 97239, U.S.A
| | - Laura M. Heiser
- Department of Biomedical Engineering, Oregon Health and Science University, Portland OR 97239, U.S.A
- Knight Cancer Institute, Oregon Health and Science University, Portland OR 97239, U.S.A
- Cancer Early Detection Advanced Research Center, Oregon Health and Science University, Portland OR 97239, U.S.A
| |
Collapse
|
2
|
Calistri NL, Liby TA, Hu Z, Zhang H, Dane M, Gross SM, Heiser LM. TNBC response to paclitaxel phenocopies interferon response which reveals cell cycle-associated resistance mechanisms. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.06.04.596911. [PMID: 38895265 PMCID: PMC11185620 DOI: 10.1101/2024.06.04.596911] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
Paclitaxel is a standard of care neoadjuvant therapy for patients with triple negative breast cancer (TNBC); however, it shows limited benefit for locally advanced or metastatic disease. Here we used a coordinated experimental-computational approach to explore the influence of paclitaxel on the cellular and molecular responses of TNBC cells. We found that escalating doses of paclitaxel resulted in multinucleation, promotion of senescence, and initiation of DNA damage induced apoptosis. Single-cell RNA sequencing (scRNA-seq) of TNBC cells after paclitaxel treatment revealed upregulation of innate immune programs canonically associated with interferon response and downregulation of cell cycle progression programs. Systematic exploration of transcriptional responses to paclitaxel and cancer-associated microenvironmental factors revealed common gene programs induced by paclitaxel, IFNB, and IFNG. Transcription factor (TF) enrichment analysis identified 13 TFs that were both enriched based on activity of downstream targets and also significantly upregulated after paclitaxel treatment. Functional assessment with siRNA knockdown confirmed that the TFs FOSL1, NFE2L2 and ELF3 mediate cellular proliferation and also regulate nuclear structure. We further explored the influence of these TFs on paclitaxel-induced cell cycle behavior via live cell imaging, which revealed altered progression rates through G1, S/G2 and M phases. We found that ELF3 knockdown synergized with paclitaxel treatment to lock cells in a G1 state and prevent cell cycle progression. Analysis of publicly available breast cancer patient data showed that high ELF3 expression was associated with poor prognosis and enrichment programs associated with cell cycle progression. Together these analyses disentangle the diverse aspects of paclitaxel response and identify ELF3 upregulation as a putative biomarker of paclitaxel resistance in TNBC.
Collapse
Affiliation(s)
- Nicholas L Calistri
- Biomedical Engineering Department, Oregon Health & Science University, Portland Oregon
| | - Tiera A. Liby
- Biomedical Engineering Department, Oregon Health & Science University, Portland Oregon
| | - Zhi Hu
- Biomedical Engineering Department, Oregon Health & Science University, Portland Oregon
| | - Hongmei Zhang
- Biomedical Engineering Department, Oregon Health & Science University, Portland Oregon
| | - Mark Dane
- Biomedical Engineering Department, Oregon Health & Science University, Portland Oregon
| | - Sean M. Gross
- Biomedical Engineering Department, Oregon Health & Science University, Portland Oregon
| | - Laura M. Heiser
- Biomedical Engineering Department, Oregon Health & Science University, Portland Oregon
- Knight Cancer Institute, Oregon Health & Science University, Portland Oregon
| |
Collapse
|
3
|
Bennett JJR, Stern AD, Zhang X, Birtwistle MR, Pandey G. Low-frequency ERK and Akt activity dynamics are predictive of stochastic cell division events. NPJ Syst Biol Appl 2024; 10:65. [PMID: 38834572 PMCID: PMC11150372 DOI: 10.1038/s41540-024-00389-7] [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/29/2024] [Accepted: 05/20/2024] [Indexed: 06/06/2024] Open
Abstract
Understanding the dynamics of intracellular signaling pathways, such as ERK1/2 (ERK) and Akt1/2 (Akt), in the context of cell fate decisions is important for advancing our knowledge of cellular processes and diseases, particularly cancer. While previous studies have established associations between ERK and Akt activities and proliferative cell fate, the heterogeneity of single-cell responses adds complexity to this understanding. This study employed a data-driven approach to address this challenge, developing machine learning models trained on a dataset of growth factor-induced ERK and Akt activity time courses in single cells, to predict cell division events. The most predictive models were developed by applying discrete wavelet transforms (DWTs) to extract low-frequency features from the time courses, followed by using Ensemble Integration, a data integration and predictive modeling framework. The results demonstrated that these models effectively predicted cell division events in MCF10A cells (F-measure=0.524, AUC=0.726). ERK dynamics were found to be more predictive than Akt, but the combination of both measurements further enhanced predictive performance. The ERK model`s performance also generalized to predicting division events in RPE cells, indicating the potential applicability of these models and our data-driven methodology for predicting cell division across different biological contexts. Interpretation of these models suggested that ERK dynamics throughout the cell cycle, rather than immediately after growth factor stimulation, were associated with the likelihood of cell division. Overall, this work contributes insights into the predictive power of intra-cellular signaling dynamics for cell fate decisions, and highlights the potential of machine learning approaches in unraveling complex cellular behaviors.
Collapse
Affiliation(s)
- Jamie J R Bennett
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Alan D Stern
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Xiang Zhang
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC, USA
| | - Marc R Birtwistle
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC, USA.
| | - Gaurav Pandey
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA.
| |
Collapse
|
4
|
Richbourg NR, Irakoze N, Kim H, Peyton SR. Outlook and opportunities for engineered environments of breast cancer dormancy. SCIENCE ADVANCES 2024; 10:eadl0165. [PMID: 38457510 PMCID: PMC10923521 DOI: 10.1126/sciadv.adl0165] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/23/2023] [Accepted: 02/01/2024] [Indexed: 03/10/2024]
Abstract
Dormant, disseminated breast cancer cells resist treatment and may relapse into malignant metastases after decades of quiescence. Identifying how and why these dormant breast cancer cells are triggered into outgrowth is a key unsolved step in treating latent, metastatic breast cancer. However, our understanding of breast cancer dormancy in vivo is limited by technical challenges and ethical concerns with triggering the activation of dormant breast cancer. In vitro models avoid many of these challenges by simulating breast cancer dormancy and activation in well-controlled, bench-top conditions, creating opportunities for fundamental insights into breast cancer biology that complement what can be achieved through animal and clinical studies. In this review, we address clinical and preclinical approaches to treating breast cancer dormancy, how precisely controlled artificial environments reveal key interactions that regulate breast cancer dormancy, and how future generations of biomaterials could answer further questions about breast cancer dormancy.
Collapse
Affiliation(s)
- Nathan R. Richbourg
- Department of Chemical Engineering, University of Massachusetts Amherst, MA 01003, USA
| | - Ninette Irakoze
- Department of Chemical Engineering, University of Massachusetts Amherst, MA 01003, USA
| | - Hyuna Kim
- Molecular and Cellular Biology Graduate Program, University of Massachusetts Amherst, MA 01003, USA
| | - Shelly R. Peyton
- Department of Chemical Engineering, University of Massachusetts Amherst, MA 01003, USA
- Molecular and Cellular Biology Graduate Program, University of Massachusetts Amherst, MA 01003, USA
- Department of Biomedical Engineering, University of Massachusetts Amherst Amherst, MA 01003, USA
| |
Collapse
|
5
|
Golconda P, Andrade-Medina M, Oberstein A. Subconfluent ARPE-19 Cells Display Mesenchymal Cell-State Characteristics and Behave like Fibroblasts, Rather Than Epithelial Cells, in Experimental HCMV Infection Studies. Viruses 2023; 16:49. [PMID: 38257749 PMCID: PMC10821009 DOI: 10.3390/v16010049] [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: 11/02/2023] [Revised: 12/20/2023] [Accepted: 12/25/2023] [Indexed: 01/24/2024] Open
Abstract
Human cytomegalovirus (HCMV) has a broad cellular tropism and epithelial cells are important physiological targets during infection. The retinal pigment epithelial cell line ARPE-19 has been used to model HCMV infection in epithelial cells for decades and remains a commonly used cell type for studying viral entry, replication, and the cellular response to infection. We previously found that ARPE-19 cells, despite being derived from an epithelial cell explant, express extremely low levels of canonical epithelial proteins, such as E-cadherin and EpCAM. Here, we perform comparative studies of ARPE-19 and additional epithelial cell lines with strong epithelial characteristics. We find that ARPE-19 cells cultured under subconfluent conditions resemble mesenchymal fibroblasts, rather than epithelial cells; this is consistent with previous studies showing that ARPE-19 cultures require extended periods of high confluency culture to maintain epithelial characteristics. By reanalyzing public gene expression data and using machine learning, we find evidence that ARPE-19 cultures maintained across many labs exhibit mesenchymal characteristics and that the majority of studies employing ARPE-19 use them in a mesenchymal state. Lastly, by performing experimental HCMV infections across mesenchymal and epithelial cell lines, we find that ARPE-19 cells behave like mesenchymal fibroblasts, producing logarithmic yields of cell-free infectious progeny, while cell lines with strong epithelial character exhibit an atypical infectious cycle and naturally restrict the production of cell-free progeny. Our work highlights important characteristics of the ARPE-19 cell line and suggests that subconfluent ARPE-19 cells may not be optimal for modeling epithelial infection with HCMV or other human viruses. It also suggests that HCMV biosynthesis and/or spread may occur quite differently in epithelial cells compared to mesenchymal cells. These differences could contribute to viral persistence or pathogenesis in epithelial tissues.
Collapse
Affiliation(s)
| | | | - Adam Oberstein
- Department of Microbiology and Immunology, College of Medicine, University of Illinois at Chicago, 835 South Wolcott Ave., Chicago, IL 60612, USA; (P.G.); (M.A.-M.)
| |
Collapse
|
6
|
Erdem C, Gross SM, Heiser LM, Birtwistle MR. MOBILE pipeline enables identification of context-specific networks and regulatory mechanisms. Nat Commun 2023; 14:3991. [PMID: 37414767 PMCID: PMC10326020 DOI: 10.1038/s41467-023-39729-2] [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: 07/27/2022] [Accepted: 06/27/2023] [Indexed: 07/08/2023] Open
Abstract
Robust identification of context-specific network features that control cellular phenotypes remains a challenge. We here introduce MOBILE (Multi-Omics Binary Integration via Lasso Ensembles) to nominate molecular features associated with cellular phenotypes and pathways. First, we use MOBILE to nominate mechanisms of interferon-γ (IFNγ) regulated PD-L1 expression. Our analyses suggest that IFNγ-controlled PD-L1 expression involves BST2, CLIC2, FAM83D, ACSL5, and HIST2H2AA3 genes, which were supported by prior literature. We also compare networks activated by related family members transforming growth factor-beta 1 (TGFβ1) and bone morphogenetic protein 2 (BMP2) and find that differences in ligand-induced changes in cell size and clustering properties are related to differences in laminin/collagen pathway activity. Finally, we demonstrate the broad applicability and adaptability of MOBILE by analyzing publicly available molecular datasets to investigate breast cancer subtype specific networks. Given the ever-growing availability of multi-omics datasets, we envision that MOBILE will be broadly useful for identification of context-specific molecular features and pathways.
Collapse
Affiliation(s)
- Cemal Erdem
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC, USA
| | - Sean M Gross
- Department of Biomedical Engineering, Oregon Health & Science University, Portland, OR, USA
| | - Laura M Heiser
- Department of Biomedical Engineering, Oregon Health & Science University, Portland, OR, USA.
| | - Marc R Birtwistle
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC, USA.
- Department of Bioengineering, Clemson University, Clemson, SC, USA.
| |
Collapse
|
7
|
Copperman J, Gross SM, Chang YH, Heiser LM, Zuckerman DM. Morphodynamical cell state description via live-cell imaging trajectory embedding. Commun Biol 2023; 6:484. [PMID: 37142678 PMCID: PMC10160022 DOI: 10.1038/s42003-023-04837-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2022] [Accepted: 04/10/2023] [Indexed: 05/06/2023] Open
Abstract
Time-lapse imaging is a powerful approach to gain insight into the dynamic responses of cells, but the quantitative analysis of morphological changes over time remains challenging. Here, we exploit the concept of "trajectory embedding" to analyze cellular behavior using morphological feature trajectory histories-that is, multiple time points simultaneously, rather than the more common practice of examining morphological feature time courses in single timepoint (snapshot) morphological features. We apply this approach to analyze live-cell images of MCF10A mammary epithelial cells after treatment with a panel of microenvironmental perturbagens that strongly modulate cell motility, morphology, and cell cycle behavior. Our morphodynamical trajectory embedding analysis constructs a shared cell state landscape revealing ligand-specific regulation of cell state transitions and enables quantitative and descriptive models of single-cell trajectories. Additionally, we show that incorporation of trajectories into single-cell morphological analysis enables (i) systematic characterization of cell state trajectories, (ii) better separation of phenotypes, and (iii) more descriptive models of ligand-induced differences as compared to snapshot-based analysis. This morphodynamical trajectory embedding is broadly applicable to the quantitative analysis of cell responses via live-cell imaging across many biological and biomedical applications.
Collapse
Affiliation(s)
- Jeremy Copperman
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, OR, 97239, USA.
| | - Sean M Gross
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, OR, 97239, USA
| | - Young Hwan Chang
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, OR, 97239, USA
- Knight Cancer Institute, Oregon Health and Science University, Portland, OR, 97239, USA
| | - Laura M Heiser
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, OR, 97239, USA.
- Knight Cancer Institute, Oregon Health and Science University, Portland, OR, 97239, USA.
| | - Daniel M Zuckerman
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, OR, 97239, USA.
- Knight Cancer Institute, Oregon Health and Science University, Portland, OR, 97239, USA.
| |
Collapse
|
8
|
Erdem C, Birtwistle MR. MEMMAL: A tool for expanding large-scale mechanistic models with machine learned associations and big datasets. FRONTIERS IN SYSTEMS BIOLOGY 2023; 3:1099413. [PMID: 38269333 PMCID: PMC10807051 DOI: 10.3389/fsysb.2023.1099413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/26/2024]
Abstract
Computational models that can explain and predict complex sub-cellular, cellular, and tissue-level drug response mechanisms could speed drug discovery and prioritize patient-specific treatments (i.e., precision medicine). Some models are mechanistic with detailed equations describing known (or supposed) physicochemical processes, while some are statistical or machine learning-based approaches, that explain datasets but have no mechanistic or causal guarantees. These two types of modeling are rarely combined, missing the opportunity to explore possibly causal but data-driven new knowledge while explaining what is already known. Here, we explore combining machine learned associations with mechanistic models to develop computational models that could more fully represent cellular behavior. In this proposed MEMMAL (MEchanistic Modeling with MAchine Learning) framework, machine learning/statistical models built using omics datasets provide predictions for new interactions between genes and proteins where there is physicochemical uncertainty. These interactions are used as a basis for new reactions in mechanistic models. As a test case, we focused on incorporating novel IFNγ/PD-L1 related associations into a large-scale mechanistic model for cell proliferation and death to better recapitulate the recently released NIH LINCS Consortium MCF10A dataset and enable description of the cellular response to checkpoint inhibitor immunotherapies. This work is a template for combining big-data-inferred interactions with mechanistic models, which could be more broadly applicable for building multi-scale precision medicine and whole cell models.
Collapse
Affiliation(s)
- Cemal Erdem
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC, United States
| | - Marc R. Birtwistle
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC, United States
- Department of Bioengineering, Clemson University, Clemson, SC, United States
| |
Collapse
|