Title: Identification of EMT signaling cross-talk and gene regulatory networks by single-cell RNA sequencing
The epithelial-to-mesenchymal transition (EMT) plays a critical role during normal development and in cancer progression. EMT is induced by various signaling pathways, including TGF-β, BMP, Wnt–β-catenin, NOTCH, Shh, and receptor tyrosine kinases. In this study, we performed single-cell RNA sequencing on MCF10A cells undergoing EMT by TGF-β1 stimulation. Our comprehensive analysis revealed that cells progress through EMT at different paces. Using pseudotime clustering reconstruction of gene-expression profiles during EMT, we found sequential and parallel activation of EMT signaling pathways. We also observed various transitional cellular states during EMT. We identified regulatory signaling nodes that drive EMT with the expression of important microRNAs and transcription factors. Using a random circuit perturbation methodology, we demonstrate that the NOTCH signaling pathway acts as a key driver of TGF-β–induced EMT. Furthermore, we demonstrate that the gene signatures of pseudotime clusters corresponding to the intermediate hybrid EMT state are associated with poor patient outcome. Overall, this study provides insight into context-specific drivers of cancer progression and highlights the complexities of the EMT process.  more » « less
Award ID(s):
2019745
PAR ID:
10233654
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
118
Issue:
19
ISSN:
0027-8424
Page Range / eLocation ID:
e2102050118
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Epithelial-to-mesenchymal transition (EMT) plays an important role in many biological processes during development and cancer. The advent of single-cell transcriptome sequencing techniques allows the dissection of dynamical details underlying EMT with unprecedented resolution. Despite several single-cell data analysis on EMT, how cell communicates and regulates dynamics along the EMT trajectory remains elusive. Using single-cell transcriptomic datasets, here we infer the cell–cell communications and the multilayer gene–gene regulation networks to analyze and visualize the complex cellular crosstalk and the underlying gene regulatory dynamics along EMT. Combining with trajectory analysis, our approach reveals the existence of multiple intermediate cell states (ICSs) with hybrid epithelial and mesenchymal features. Analyses on the time-series datasets from cancer cell lines with different inducing factors show that the induced EMTs are context-specific: the EMT induced by transforming growth factor B1 (TGFB1) is synchronous, whereas the EMTs induced by epidermal growth factor and tumor necrosis factor are asynchronous, and the responses of TGF-β pathway in terms of gene expression regulations are heterogeneous under different treatments or among various cell states. Meanwhile, network topology analysis suggests that the ICSs during EMT serve as the signaling in cellular communication under different conditions. Interestingly, our analysis of a mouse skin squamous cell carcinoma dataset also suggests regardless of the significant discrepancy in concrete genes between in vitro and in vivo EMT systems, the ICSs play dominant role in the TGF-β signaling crosstalk. Overall, our approach reveals the multiscale mechanisms coupling cell–cell communications and gene–gene regulations responsible for complex cell-state transitions. 
    more » « less
  2. Hybrid epithelial/mesenchymal cells (E/M) are key players in aggressive cancer metastasis. It remains a challenge to understand how these cell states, which are mostly non-existent in healthy tissue, become stable phenotypes participating in collective cancer migration. The transcription factor Nrf2, which is associated with tumor progression and resistance to therapy, appears to be central to this process. Here, using a combination of immunocytochemistry, single cell biosensors, and computational modeling, we show that Nrf2 functions as a phenotypic stability factor for hybrid E/M cells by inhibiting a complete epithelial-mesenchymal transition (EMT) during collective cancer migration. We also demonstrate that Nrf2 and EMT signaling are spatially coordinated near the leading edge. In particular, computational analysis of an Nrf2-EMT-Notch network and experimental modulation of Nrf2 by pharmacological treatment or CRISPR/Cas9 gene editing reveal that Nrf2 stabilizes a hybrid E/M phenotype which is maximally observed in the interior region immediately behind the leading edge. We further demonstrate that the Nrf2-EMT-Notch network enhances Dll4 and Jagged1 expression at the leading edge, which correlates with the formation of leader cells and protruding tips. Altogether, our results provide direct evidence that Nrf2 acts as a phenotypic stability factor in restricting complete EMT and plays an important role in coordinating collective cancer migration. 
    more » « less
  3. ABSTRACT The methyltransferase enhancer of zeste homolog 2 (EZH2) regulates gene expression, and aberrant EZH2 expression and signaling can drive fibrosis and cancer. However, it is not clear how chemical and mechanical signals are integrated to regulate EZH2 and gene expression. We show that culture of cells on stiff matrices in concert with transforming growth factor (TGF)-β1 promotes nuclear localization of EZH2 and an increase in the levels of the corresponding histone modification, H3K27me3, thereby regulating gene expression. EZH2 activity and expression are required for TGFβ1- and stiffness-induced increases in H3K27me3 levels as well as for morphological and gene expression changes associated with epithelial–mesenchymal transition (EMT). Inhibition of Rho associated kinase (ROCK) proteins or myosin II signaling attenuates TGFβ1-induced nuclear localization of EZH2 and decreases H3K27me3 levels in cells cultured on stiff substrata, suggesting that cellular contractility, in concert with a major cancer signaling regulator TGFβ1, modulates EZH2 subcellular localization. These findings provide a contractility-dependent mechanism by which matrix stiffness and TGFβ1 together mediate EZH2 signaling to promote EMT. 
    more » « less
  4. Intermediate cell states (ICSs) during the epithelial–mesenchymal transition (EMT) are emerging as a driving force of cancer invasion and metastasis. ICSs typically exhibit hybrid epithelial/mesenchymal characteristics as well as cancer stem cell (CSC) traits including proliferation and drug resistance. Here, we analyze several single-cell RNA-seq (scRNA-seq) datasets to investigate the relation between several axes of cancer progression including EMT, CSC traits, and cell–cell signaling. To accomplish this task, we integrate computational methods for clustering and trajectory inference with analysis of EMT gene signatures, CSC markers, and cell–cell signaling pathways, and highlight conserved and specific processes across the datasets. Our analysis reveals that “standard” measures of pluripotency often used in developmental contexts do not necessarily correlate with EMT progression and expression of CSC-related markers. Conversely, an EMT circuit energy that quantifies the co-expression of epithelial and mesenchymal genes consistently increases along EMT trajectories across different cancer types and anatomical locations. Moreover, despite the high context specificity of signal transduction across different cell types, cells undergoing EMT always increased their potential to send and receive signals from other cells. 
    more » « less
  5. Metastatic melanoma presents a formidable challenge in oncology due to its high invasiveness and resistance to current treatments. Central to its ability to metastasize is the Notch signaling pathway, which, when activated through direct cell–cell interactions, propels cells into a metastatic state through mechanisms akin to the epithelial-mesenchymal transition (EMT). While the upregulation of miR-222 has been identified as a critical step in this metastatic progression, the mechanism through which this upregulation persists in the absence of active Notch signaling remains unclear. Here we introduce a dynamical system model that integrates miR-222 gene regulation with histone feedback mechanisms. Through computational analysis spanning both sustained and pulsatile ligand inputs, we delineate the non-linear decision boundaries that govern melanoma cell fate transitions, taking into account the dynamics of Notch signaling and the role of epigenetic modifications. Dimensional analysis reduces the 11-parameter system to three critical control groups governing chromatin modification rates and feedback strengths, providing a theoretical framework for parameter selection in the absence of complete kinetic measurements. Global sensitivity analysis identifies PRC2-mediated methylation and KDM5A-mediated demethylation as the dominant control parameters, while stochastic simulations show population heterogeneity consistent with the variable EMT responses observed in cancer cell populations. Our analysis examines the interplay between Notch signaling pathways and epigenetic regulation in dictating melanoma cell fate. 
    more » « less