An Ovol2-Zeb1 Mutual Inhibitory Circuit Governs Bidirectional and Multi-step Transition between Epithelial and Mesenchymal States

Energy Landscape Reveals the Underlying Mechanism of Cancer-Adipose Conversion in Gene Network Models

Cancer is a systemic heterogeneous disease involving complex molecular networks. Tumor formation involves an epithelial-mesenchymal transition (EMT), which promotes both metastasis and plasticity of cancer cells. Recent experiments have proposed that cancer cells can be transformed into adipocytes via a combination of drugs. However, the underlying mechanisms for how these drugs work, from a molecular network perspective, remain elusive. To reveal the mechanism of cancer-adipose conversion (CAC), this study adopts a systems biology approach by combing mathematical modeling and molecular experiments, based on underlying molecular regulatory networks. Four types of attractors are identified, corresponding to epithelial (E), mesenchymal (M), adipose (A) and partial/intermediate EMT (P) cell states on the CAC landscape. Landscape and transition path results illustrate that intermediate states play critical roles in the cancer to adipose transition. Through a landscape control approach, two new therapeutic strategies for drug combinations are identified, that promote CAC. These predictions are verified by molecular experiments in different cell lines. The combined computational and experimental approach provides a powerful tool to explore molecular mechanisms for cell fate transitions in cancer networks. The results reveal underlying mechanisms of intermediate cell states that govern the CAC, and identified new potential drug combinations to induce cancer adipogenesis.

Expression Patterns of Grainyhead-Like 2 and Ovo-Like 2 in Mouse Mammary Gland Development During Pregnancy, Lactation, and Weaning

Grainyhead-like 2 (Grhl2) is a transcription factor that regulates cell adhesion genes in mammary ductal development and serves as a repressor of the epithelial-mesenchymal transition. Conversely, Ovo-like2 (Ovol2) is a target gene of Grhl2 but functions as a substitute in Grhl2-deficient mice, facilitating successful epithelial barrier formation and lumen expansion in kidney-collecting ductal epithelial cells. Our objective was to examine the expression patterns of Grhl2, Ovol2, and their associated genes during the intricate phases of mouse mammary gland development. The mRNA expression of Grhl2 and Ovol2 increased after pregnancy. We observed Grhl2 protein presence in the epithelial cell's region, coinciding with acini formation, and its signal significantly correlated with E-cadherin (Cdh1) expression. However, Ovol2 was present in the epithelial region without a correlation with Cdh1. Similarly, Zeb1, a mesenchymal transcription factor, showed Cdh1-independent expression. Subsequently, we explored the interaction between Rab25, a small G protein, and Grhl2/Ovol2. The expressions of Grhl2 and Ovol2 exhibited a strong correlation with Rab25 and claudin-4, a tight junction protein. These findings suggest that Grhl2 and Ovol2 may collaborate to regulate genes associated with cell adhesion and are crucial for maintaining epithelial integrity during the different phases of mammary gland development.

RNA regulatory mechanisms controlling TGF-β signaling and EMT in cancer

Epithelial-mesenchymal transition (EMT) is a major contributor to metastatic progression and is prominently regulated by TGF-β signalling. Both EMT and TGF-β pathway components are tightly controlled by non-coding RNAs - including microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) - that collectively have major impacts on gene expression and resulting cellular states. While miRNAs are the best characterised regulators of EMT and TGF-β signaling and the miR-200-ZEB1/2 feedback loop plays a central role, important functions for lncRNAs and circRNAs are also now emerging. This review will summarise our current understanding of the roles of non-coding RNAs in EMT and TGF-β signaling with a focus on their functions in cancer progression.

ATF3 characterizes aggressive drug-tolerant persister cells in HGSOC

High-grade serous ovarian cancer (HGSOC) represents the most common and lethal subtype of ovarian cancer. Despite initial response to platinum-based standard therapy, patients commonly suffer from relapse that likely originates from drug-tolerant persister (DTP) cells. We generated isogenic clones of treatment-naïve and cisplatin-tolerant persister HGSOC cells. In addition, single-cell RNA sequencing of barcoded cells was performed in a xenograft model with HGSOC cell lines after platinum-based therapy. Published single-cell RNA-sequencing data from neo-adjuvant and non-treated HGSOC patients and patient data from TCGA were analyzed. DTP-derived cells exhibited morphological alterations and upregulation of epithelial-mesenchymal transition (EMT) markers. An aggressive subpopulation of DTP-derived cells showed high expression of the stress marker ATF3. Knockdown of ATF3 enhanced the sensitivity of aggressive DTP-derived cells to cisplatin-induced cell death, implying a role for ATF3 stress response in promoting a drug tolerant persister cell state. Furthermore, single cell lineage tracing to detect transcriptional changes in a HGSOC cell line-derived xenograft relapse model showed that cells derived from relapsed solid tumors express increased levels of EMT and multiple endoplasmic reticulum (ER) stress markers, including ATF3. Single cell RNA sequencing of epithelial cells from four HGSOC patients also identified a small cell population resembling DTP cells in all samples. Moreover, analysis of TCGA data from 259 HGSOC patients revealed a significant progression-free survival advantage for patients with low expression of the ATF3-associated partial EMT genes. These findings suggest that increased ATF3 expression together with partial EMT promote the development of aggressive DTP, and thereby relapse in HGSOC patients.

Data- and theory-driven approaches for understanding paths of epithelial-mesenchymal transition

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.

Cancer cell plasticity: from cellular, molecular, and genetic mechanisms to tumor heterogeneity and drug resistance

Cancer is a complex disease displaying a variety of cell states and phenotypes. This diversity, known as cancer cell plasticity, confers cancer cells the ability to change in response to their environment, leading to increased tumor diversity and drug resistance. This review explores the intricate landscape of cancer cell plasticity, offering a deep dive into the cellular, molecular, and genetic mechanisms that underlie this phenomenon. Cancer cell plasticity is intertwined with processes such as epithelial-mesenchymal transition and the acquisition of stem cell-like features. These processes are pivotal in the development and progression of tumors, contributing to the multifaceted nature of cancer and the challenges associated with its treatment. Despite significant advancements in targeted therapies, cancer cell adaptability and subsequent therapy-induced resistance remain persistent obstacles in achieving consistent, successful cancer treatment outcomes. Our review delves into the array of mechanisms cancer cells exploit to maintain plasticity, including epigenetic modifications, alterations in signaling pathways, and environmental interactions. We discuss strategies to counteract cancer cell plasticity, such as targeting specific cellular pathways and employing combination therapies. These strategies promise to enhance the efficacy of cancer treatments and mitigate therapy resistance. In conclusion, this review offers a holistic, detailed exploration of cancer cell plasticity, aiming to bolster the understanding and approach toward tackling the challenges posed by tumor heterogeneity and drug resistance. As articulated in this review, the delineation of cellular, molecular, and genetic mechanisms underlying tumor heterogeneity and drug resistance seeks to contribute substantially to the progress in cancer therapeutics and the advancement of precision medicine, ultimately enhancing the prospects for effective cancer treatment and patient outcomes.

Ovol2 promoter mutations in mice and human illuminate species-specific phenotypic divergence

Pathogenic variants in the highly conserved OVOL2 promoter region cause posterior polymorphous corneal dystrophy (PPCD) 1 by inducing an ectopic expression of the endothelial OVOL2 mRNA. Here we produced an allelic series of Ovol2 promoter mutations in the mouse model including the heterozygous c.-307T>C variant (RefSeq NM_021220.4) causing PPCD1 in humans. Despite the high evolutionary conservation of the Ovol2 promoter, only some alterations of its sequence had phenotypic consequences in mice. Four independent sequence variants in the distal part of the Ovol2 promoter had no significant effect on endothelial Ovol2 mRNA level or caused any ocular phenotype. In contrast, the mutation c.-307T>C resulted in increased Ovol2 expression in the corneal endothelium. However, only a small fraction of adult mice c.-307T>C heterozygotes developed ocular phenotypes such as irido-corneal adhesions, and corneal opacity. Interestingly, phenotypic penetrance was increased at embryonic stages. Notably, c.-307T>C mutation is located next to the Ovol1/Ovol2 transcription factor binding site. Mice carrying an allele with a deletion encompassing the Ovol2 binding site c.-307_-320del showed significant Ovol2 gene upregulation in the cornea endothelium and exhibited phenotypes similar to the c.-307T>C mutation. In conclusion, although the mutations c.-307T>C and -307_-320del lead to a comparably strong increase in endothelial Ovol2 expression as seen in PPCD1 patients, endothelial dystrophy was not observed in the mouse model, implicating species-specific differences in endothelial cell biology. Nonetheless, the emergence of dominant ocular phenotypes associated with Ovol2 promoter variants in mice implies a potential role of this gene in eye development and disease.

Identifying critical regulatory interactions in cell fate decision and transition by systematic perturbation analysis

One of the most significant challenges in biology is to elucidate the roles of various regulatory interactions in cell fate decision and transition. However, it remains to be fully clarified how they cooperate and determine fate transition. Here, a general framework based on statistical analysis and bifurcation theory is proposed to identify crucial regulatory interactions and how they play decisive roles in fate transition. More exactly, specific feedback loops determine occurrence of bifurcations by which cell fate transition can be realized. While regulatory interactions in the feedback loops determine the direction of transition. In addition, two-parameter bifurcation analysis further provides detailed understanding of how the fate transition based on statistical analysis occurs. Statistical analysis can also be used to reveal synergistic combinatorial perturbations by which fate transition can be more efficiently realized. The integrative analysis approach can be used to identify critical regulatory interactions in cell fate transition and reveal how specific cell fate transition occurs. To verify feasibility of the approach, the epithelial to mesenchymal transition (EMT) network is chosen as an illustrative example. In agreement with experimental observations, the approach reveals some critical regulatory interactions and underlying mechanisms in cell fate determination and transitions between three states. The approach can also be applied to analyze other regulatory networks related to cell fate decision and transition.

Crosstalk between chromatin and Shavenbaby defines transcriptional output along the Drosophila intestinal stem cell lineage

The transcription factor Shavenbaby (Svb), the only member of the OvoL family in Drosophila, controls the fate of various epithelial embryonic cells and adult stem cells. Post-translational modification of Svb produces two protein isoforms, Svb-ACT and Svb-REP, which promote adult intestinal stem cell renewal or differentiation, respectively. To define Svb mode of action, we used engineered cell lines and develop an unbiased method to identify Svb target genes across different contexts. Within a given cell type, Svb-ACT and Svb-REP antagonistically regulate the expression of a set of target genes, binding specific enhancers whose accessibility is constrained by chromatin landscape. Reciprocally, Svb-REP can influence local chromatin marks of active enhancers to help repressing target genes. Along the intestinal lineage, the set of Svb target genes progressively changes, together with chromatin accessibility. We propose that Svb-ACT-to-REP transition promotes enterocyte differentiation of intestinal stem cells through direct gene regulation and chromatin remodeling.

The role of flavonoids in the regulation of epithelial-mesenchymal transition in cancer: A review on targeting signaling pathways and metastasis

One of the hallmarks of cancer is metastasis, a process that entails the spread of cancer cells to distant regions in the body, culminating in tumor formation in secondary organs. Importantly, the proinflammatory environment surrounding cancer cells further contributes to cancer cell transformation and extracellular matrix destruction. During metastasis, front-rear polarity and emergence of migratory and invasive features are manifestations of epithelial-mesenchymal transition (EMT). A variety of transcription factors (TFs) are implicated in the execution of EMT, the most prominent belonging to the Snail Family Transcriptional Repressor (SNAI) and Zinc Finger E-Box Binding Homeobox (ZEB) families of TFs. These TFs are regulated by interaction with specific microRNAs (miRNAs), as miR34 and miR200. Among the several secondary metabolites produced in plants, flavonoids constitute a major group of bioactive molecules, with several described effects including antioxidant, antiinflammatory, antidiabetic, antiobesogenic, and anticancer effects. This review scrutinizes the modulatory role of flavonoids on the activity of SNAI/ZEB TFs and on their regulatory miRNAs, miR-34, and miR-200. The modulatory role of flavonoids can attenuate mesenchymal features and stimulate epithelial features, thereby inhibiting and reversing EMT. Moreover, this modulation is concomitant with the attenuation of signaling pathways involved in diverse processes as cell proliferation, cell growth, cell cycle progression, apoptosis inhibition, morphogenesis, cell fate, cell migration, cell polarity, and wound healing. The antimetastatic potential of these versatile compounds is emerging and represents an opportunity for the synthesis of more specific and potent agents.

Dynamical hallmarks of cancer: Phenotypic switching in melanoma and epithelial-mesenchymal plasticity

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach.

OVOL2: an epithelial lineage determiner with emerging roles in energy homeostasis

Ovo like zinc finger 2 (OVOL2) is an evolutionarily conserved regulator of epithelial lineage determination and differentiation during embryogenesis. OVOL2 binds to DNA using zinc-finger domains to suppress epithelial-mesenchymal transition (EMT), which is critical for tumor metastasis. However, recent studies have suggested some noncanonical roles of OVOL2 that do not rely on the DNA binding of zinc-finger domains or regulation of EMT. OVOL2 and EMT regulators have emerging roles in adipogenesis, thermogenesis, and lipid metabolism. Here, we review different roles of OVOL2 from embryo development to adult tissue homeostasis, and discuss how OVOL2 and other EMT regulators orchestrate a regulatory network to control energy homeostasis. Last, we propose potential applications of targeting OVOL2 to reduce human obesity.

Unraveling the Molecular Puzzle: Exploring Gene Networks across Diverse EMT Status of Cell Lines

Understanding complex disease mechanisms requires a comprehensive understanding of the gene regulatory networks, as complex diseases are often characterized by the dysregulation and dysfunction of molecular networks, rather than abnormalities in single genes. Specifically, the exploration of cell line-specific gene networks can provide essential clues for precision medicine, as this methodology can uncover molecular interplays specific to particular cell line statuses, such as drug sensitivity, cancer progression, etc. In this article, we provide a comprehensive review of computational strategies for cell line-specific gene network analysis: (1) cell line-specific gene regulatory network estimation and analysis of gene networks under varying epithelial-mesenchymal transition (EMT) statuses of cell lines; and (2) an explainable artificial intelligence approach for interpreting the estimated massive multiple EMT-status-specific gene networks. The objective of this review is to help readers grasp the concept of computational network biology, which holds significant implications for precision medicine by offering crucial clues.

ZHX2 deficiency enriches hybrid MET cells through regulating E-cadherin expression

Growing evidence indicates that the epithelial to mesenchymal (E/M) hybrid state plays a key role in tumorigenesis. Importantly, a hybrid mesenchymal to epithelial transition (MET) state in which individual cells express both epithelial and mesenchymal markers was recently identified in vivo, further strengthening the bonds between the hybrid EMT state and cancer progression. However, the role and the molecular mechanisms by which the hybrid MET state is maintained in triple-negative breast cancer cells (TNBC) remain elusive. Here, we find that loss of ZHX2 expression results in the hybrid MET phenotype in mesenchymal TNBC cells. Mechanistically, through directly binding to the CDH1 promoter, depletion of ZHX2 specifically reactivates expression of CDH1 encoding E-cadherin, an epithelial marker that is crucial for maintaining epithelial phenotype. Functionally, loss of ZHX2 expression enriches the hybrid MET cells and inhibits the migration and dissemination of TNBC cells or organoids, which could be reversed by restoration of E-cadherin. Moreover, depletion of ZHX2 suppresses lung metastasis in preclinical models of TNBC. In patients with TNBC, ZHX2 expression was amplified and negatively correlated with the expression of E-cadherin. These findings suggest that loss of ZHX2 promotes the hybrid MET state to impair TNBC progression.

Role of epithelial-mesenchymal transition factor SNAI1 and its targets in ovarian cancer aggressiveness

Ovarian cancer remains the most lethal gynecologic malignancy in the USA. For over twenty years, epithelial-mesenchymal transition (EMT) has been characterized extensively in development and disease. The dysregulation of this process in cancer has been identified as a mechanism by which epithelial tumors become more aggressive, allowing them to survive and invade distant tissues. This occurs in part due to the increased expression of the EMT transcription factor, SNAI1 (Snail). In the case of epithelial ovarian cancer, Snail has been shown to contribute to cancer invasion, stemness, chemoresistance, and metabolic changes. Thus, in this review, we focus on summarizing current findings on the role of EMT (specifically, factors downstream of Snail) in determining ovarian cancer aggressiveness.

Resistance to mesenchymal reprogramming sustains clonal propagation in metastatic breast cancer

The acquisition of mesenchymal traits is considered a hallmark of breast cancer progression. However, the functional relevance of epithelial-to-mesenchymal transition (EMT) remains controversial and context dependent. Here, we isolate epithelial and mesenchymal populations from human breast cancer metastatic biopsies and assess their functional potential in vivo. Strikingly, progressively decreasing epithelial cell adhesion molecule (EPCAM) levels correlate with declining disease propagation. Mechanistically, we find that persistent EPCAM expression marks epithelial clones that resist EMT induction and propagate competitively. In contrast, loss of EPCAM defines clones arrested in a mesenchymal state, with concomitant suppression of tumorigenicity and metastatic potential. This dichotomy results from distinct clonal trajectories impacting global epigenetic programs that are determined by the interplay between human ZEB1 and its target GRHL2. Collectively, our results indicate that susceptibility to irreversible EMT restrains clonal propagation, whereas resistance to mesenchymal reprogramming sustains disease spread in multiple models of human metastatic breast cancer, including patient-derived cells in vivo.

EMT/MET plasticity in cancer and Go-or-Grow decisions in quiescence: the two sides of the same coin?

Epithelial mesenchymal transition (EMT) and mesenchymal epithelial transition (MET) are genetic determinants of cellular plasticity. These programs operate in physiological (embryonic development, wound healing) and pathological (organ fibrosis, cancer) conditions. In cancer, EMT and MET interfere with various signalling pathways at different levels. This results in gross alterations in the gene expression programs, which affect most, if not all hallmarks of cancer, such as response to proliferative and death-inducing signals, tumorigenicity, and cell stemness. EMT in cancer cells involves large scale reorganisation of the cytoskeleton, loss of epithelial integrity, and gain of mesenchymal traits, such as mesenchymal type of cell migration. In this regard, EMT/MET plasticity is highly relevant to the Go-or-Grow concept, which postulates the dichotomous relationship between cell motility and proliferation. The Go-or-Grow decisions are critically important in the processes in which EMT/MET plasticity takes the central stage, mobilisation of stem cells during wound healing, cancer relapse, and metastasis. Here we outline the maintenance of quiescence in stem cell and metastatic niches, focusing on the implication of EMT/MET regulatory networks in Go-or-Grow switches. In particular, we discuss the analogy between cells residing in hybrid quasi-mesenchymal states and G_Alert, an intermediate phase allowing quiescent stem cells to enter the cell cycle rapidly.

A role for partial epithelial-to-mesenchymal transition in enabling stemness in homeostasis and cancer

Stem cells have self-renewal capacities and the ability to give rise to differentiated cells thereby sustaining tissues during homeostasis and injury. This structural hierarchy extends to tumours which harbor stem-like cells deemed cancer stem cells that propagate the tumour and drive metastasis and relapse. The process of epithelial-to-mesenchymal transition (EMT), which plays an important role in development and cancer cell migration, was shown to be correlated with stemness in both homeostasis and cancer indicating that stemness can be acquired and is not necessarily an intrinsic trait. Nowadays it is experimentally proven that the activation of an EMT program does not necessarily drive cells towards a fully mesenchymal phenotype but rather to hybrid E/M states. This review offers the latest advances in connecting the EMT status and stem-cell state of both non-transformed and cancer cells. Recent literature clearly shows that hybrid EMT states have a higher probability of acquiring stem cell traits. The position of a cell along the EMT-axis which coincides with a stem cell-like state is known as the stemness window. We show how the original EMT-state of a cell dictates the EMT/MET inducing programmes required to reach stemness. Lastly we present the mechanism of stemness regulation and the regulatory feedback loops which position cells at a certain EMT state along the EMT axis.

Dynamics of Epithelial-Mesenchymal Plasticity: What Have Single-Cell Investigations Elucidated So Far?

Epithelial-mesenchymal plasticity (EMP) is a key driver of cancer metastasis and therapeutic resistance, through which cancer cells can reversibly and dynamically alter their molecular and functional traits along the epithelial-mesenchymal spectrum. While cells in the epithelial phenotype are usually tightly adherent, less metastatic, and drug-sensitive, those in the hybrid epithelial/mesenchymal and/or mesenchymal state are more invasive, migratory, drug-resistant, and immune-evasive. Single-cell studies have emerged as a powerful tool in gaining new insights into the dynamics of EMP across various cancer types. Here, we review many recent studies that employ single-cell analysis techniques to better understand the dynamics of EMP in cancer both in vitro and in vivo. These single-cell studies have underlined the plurality of trajectories cells can traverse during EMP and the consequent heterogeneity of hybrid epithelial/mesenchymal phenotypes seen at both preclinical and clinical levels. They also demonstrate how diverse EMP trajectories may exhibit hysteretic behavior and how the rate of such cell-state transitions depends on the genetic/epigenetic background of recipient cells, as well as the dose and/or duration of EMP-inducing growth factors. Finally, we discuss the relationship between EMP and patient survival across many cancer types. We also present a next set of questions related to EMP that could benefit much from single-cell observations and pave the way to better tackle phenotypic switching and heterogeneity in clinic.

Involvement of Epithelial-Mesenchymal Transition Genes in Small Cell Lung Cancer Phenotypic Plasticity

Small cell lung cancer (SCLC) is an aggressive cancer recalcitrant to treatment, arising predominantly from epithelial pulmonary neuroendocrine (NE) cells. Intratumor heterogeneity plays critical roles in SCLC disease progression, metastasis, and treatment resistance. At least five transcriptional SCLC NE and non-NE cell subtypes were recently defined by gene expression signatures. Transition from NE to non-NE cell states and cooperation between subtypes within a tumor likely contribute to SCLC progression by mechanisms of adaptation to perturbations. Therefore, gene regulatory programs distinguishing SCLC subtypes or promoting transitions are of great interest. Here, we systematically analyze the relationship between SCLC NE/non-NE transition and epithelial to mesenchymal transition (EMT)-a well-studied cellular process contributing to cancer invasiveness and resistance-using multiple transcriptome datasets from SCLC mouse tumor models, human cancer cell lines, and tumor samples. The NE SCLC-A2 subtype maps to the epithelial state. In contrast, SCLC-A and SCLC-N (NE) map to a partial mesenchymal state (M1) that is distinct from the non-NE, partial mesenchymal state (M2). The correspondence between SCLC subtypes and the EMT program paves the way for further work to understand gene regulatory mechanisms of SCLC tumor plasticity with applicability to other cancer types.

Heterogeneity and plasticity of epithelial-mesenchymal transition (EMT) in cancer metastasis: Focusing on partial EMT and regulatory mechanisms

Epithelial-mesenchymal transition (EMT) or mesenchymal-epithelial transition (MET) plays critical roles in cancer metastasis. Recent studies, especially those based on single-cell sequencing, have revealed that EMT is not a binary process, but a heterogeneous and dynamic disposition with intermediary or partial EMT states. Multiple double-negative feedback loops involved by EMT-related transcription factors (EMT-TFs) have been identified. These feedback loops between EMT drivers and MET drivers finely regulate the EMT transition state of the cell. In this review, the general characteristics, biomarkers and molecular mechanisms of different EMT transition states were summarized. We additionally discussed the direct and indirect roles of EMT transition state in tumour metastasis. More importantly, this article provides direct evidence that the heterogeneity of EMT is closely related to the poor prognosis in gastric cancer. Notably, a seesaw model was proposed to explain how tumour cells regulate themselves to remain in specific EMT transition states, including epithelial state, hybrid/intermediate state and mesenchymal state. Additionally, this article also provides a review of the current status, limitations and future perspectives of EMT signalling in clinical applications.

Genomic and microenvironmental heterogeneity shaping epithelial-to-mesenchymal trajectories in cancer

The epithelial to mesenchymal transition (EMT) is a key cellular process underlying cancer progression, with multiple intermediate states whose molecular hallmarks remain poorly characterised. To fill this gap, we present a method to robustly evaluate EMT transformation in individual tumours based on transcriptomic signals. We apply this approach to explore EMT trajectories in 7180 tumours of epithelial origin and identify three macro-states with prognostic and therapeutic value, attributable to epithelial, hybrid E/M and mesenchymal phenotypes. We show that the hybrid state is relatively stable and linked with increased aneuploidy. We further employ spatial transcriptomics and single cell datasets to explore the spatial heterogeneity of EMT transformation and distinct interaction patterns with cytotoxic, NK cells and fibroblasts in the tumour microenvironment. Additionally, we provide a catalogue of genomic events underlying distinct evolutionary constraints on EMT transformation. This study sheds light on the aetiology of distinct stages along the EMT trajectory, and highlights broader genomic and environmental hallmarks shaping the mesenchymal transformation of primary tumours.

Parallelized multidimensional analytic framework applied to mammary epithelial cells uncovers regulatory principles in EMT

A proper understanding of disease etiology will require longitudinal systems-scale reconstruction of the multitiered architecture of eukaryotic signaling. Here we combine state-of-the-art data acquisition platforms and bioinformatics tools to devise PAMAF, a workflow that simultaneously examines twelve omics modalities, i.e., protein abundance from whole-cells, nucleus, exosomes, secretome and membrane; N-glycosylation, phosphorylation; metabolites; mRNA, miRNA; and, in parallel, single-cell transcriptomes. We apply PAMAF in an established in vitro model of TGFβ-induced epithelial to mesenchymal transition (EMT) to quantify >61,000 molecules from 12 omics and 10 timepoints over 12 days. Bioinformatics analysis of this EMT-ExMap resource allowed us to identify; -topological coupling between omics, -four distinct cell states during EMT, -omics-specific kinetic paths, -stage-specific multi-omics characteristics, -distinct regulatory classes of genes, -ligand-receptor mediated intercellular crosstalk by integrating scRNAseq and subcellular proteomics, and -combinatorial drug targets (e.g., Hedgehog signaling and CAMK-II) to inhibit EMT, which we validate using a 3D mammary duct-on-a-chip platform. Overall, this study provides a resource on TGFβ signaling and EMT.

Network topology metrics explaining enrichment of hybrid epithelial mesenchymal phenotypes in metastasis

Epithelial to Mesenchymal Transition (EMT) and its reverse—Mesenchymal to Epithelial Transition (MET) are hallmarks of metastasis. Cancer cells use this reversible cellular programming to switch among Epithelial (E), Mesenchymal (M), and hybrid Epithelial/Mesenchymal (hybrid E/M) state(s) and seed tumors at distant sites. Hybrid E/M cells are often more aggressive and metastatic than the “pure” E and M cells. Thus, identifying mechanisms to inhibit hybrid E/M cells can be promising in curtailing metastasis. While multiple gene regulatory networks (GRNs) based mathematical models for EMT/MET have been developed recently, identifying topological signatures enriching hybrid E/M phenotypes remains to be done. Here, we investigate the dynamics of 13 different GRNs and report an interesting association between “hybridness” and the number of negative/positive feedback loops across the networks. While networks having more negative feedback loops favor hybrid phenotype(s), networks having more positive feedback loops (PFLs) or many HiLoops–specific combinations of PFLs, support terminal (E and M) phenotypes. We also establish a connection between “hybridness” and network-frustration by showing that hybrid phenotypes likely result from non-reinforcing interactions among network nodes (genes) and therefore tend to be more frustrated (less stable). Our analysis, thus, identifies network topology-based signatures that can give rise to, as well as prevent, the emergence of hybrid E/M phenotype in GRNs underlying EMP. Our results can have implications in terms of targeting specific interactions in GRNs as a potent way to restrict switching to the hybrid E/M phenotype(s) to curtail metastasis.

Landscape of epithelial-mesenchymal plasticity as an emergent property of coordinated teams in regulatory networks

Elucidating the design principles of regulatory networks driving cellular decision-making has fundamental implications in mapping and eventually controlling cell-fate decisions. Despite being complex, these regulatory networks often only give rise to a few phenotypes. Previously, we identified two 'teams' of nodes in a small cell lung cancer regulatory network that constrained the phenotypic repertoire and aligned strongly with the dominant phenotypes obtained from network simulations (Chauhan et al., 2021). However, it remained elusive whether these 'teams' exist in other networks, and how do they shape the phenotypic landscape. Here, we demonstrate that five different networks of varying sizes governing epithelial-mesenchymal plasticity comprised of two 'teams' of players - one comprised of canonical drivers of epithelial phenotype and the other containing the mesenchymal inducers. These 'teams' are specific to the topology of these regulatory networks and orchestrate a bimodal phenotypic landscape with the epithelial and mesenchymal phenotypes being more frequent and dynamically robust to perturbations, relative to the intermediary/hybrid epithelial/mesenchymal ones. Our analysis reveals that network topology alone can contain information about corresponding phenotypic distributions, thus obviating the need to simulate them. We propose 'teams' of nodes as a network design principle that can drive cell-fate canalization in diverse decision-making processes.

Cooperative RNA degradation stabilizes intermediate epithelial-mesenchymal states and supports a phenotypic continuum

Multiple intermediate epithelial-mesenchymal transition (EMT) states reflecting hybrid epithelial and mesenchymal phenotypes were observed in physiological and pathological conditions. Previous theoretical models explaining multiple EMT states rely on regulatory loops involving transcriptional feedback, which produce three or four attractors. This is incompatible with the observed continuum-like EMT spectrum. Here, we used mass-action-based models to describe post-transcriptional regulations, finding that cooperative RNA degradation via multiple microRNA binding sites can generate four-attractor systems without transcriptional feedback. Furthermore, the newly identified intermediates-enabling circuits are common in the EMT regulatory network, and they can synergize with transcriptional feedback to support phenotypic continuum. Finally, our model predicted a role of miR-101 in multistate EMT, and we identified evidence from single-cell RNA-sequencing data that support the prediction. Our work reveals a previously unknown role of cooperative RNA degradation and microRNAs in EMT, providing a framework that can bridge the gap between mechanistic models and single-cell experiments.

Stochastic population dynamics of cancer stemness and adaptive response to therapies

Intratumoral heterogeneity can exist along multiple axes: Cancer stem cells (CSCs)/non-CSCs, drug-sensitive/drug-tolerant states, and a spectrum of epithelial-hybrid-mesenchymal phenotypes. Further, these diverse cell-states can switch reversibly among one another, thereby posing a major challenge to therapeutic efficacy. Therefore, understanding the origins of phenotypic plasticity and heterogeneity remains an active area of investigation. While genomic components (mutations, chromosomal instability) driving heterogeneity have been well-studied, recent reports highlight the role of non-genetic mechanisms in enabling both phenotypic plasticity and heterogeneity. Here, we discuss various processes underlying phenotypic plasticity such as stochastic gene expression, chromatin reprogramming, asymmetric cell division and the presence of multiple stable gene expression patterns ('attractors'). These processes can facilitate a dynamically evolving cell population such that a subpopulation of (drug-tolerant) cells can survive lethal drug exposure and recapitulate population heterogeneity on drug withdrawal, leading to relapse. These drug-tolerant cells can be both pre-existing and also induced by the drug itself through cell-state reprogramming. The dynamics of cell-state transitions both in absence and presence of the drug can be quantified through mathematical models. Such a dynamical systems approach to elucidating patterns of intratumoral heterogeneity by integrating longitudinal experimental data with mathematical models can help design effective combinatorial and/or sequential therapies for better clinical outcomes.

The role of tumor-associated macrophages in oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is a common head and neck cancer with a high recurrence rate and a low 5-year survival rate. Tumor-associated macrophages (TAMs) are important immune cells in the tumor microenvironment, which play an important role in the progression of many tumors. This article reviews the origin, and the role of TAMs in the invasion, metastasis, angiogenesis and immunosuppression of OSCC. Therapeutic strategies targeting TAMs are also discussed in hopes of providing new ideas for the treatment of OSCC.

A data-assimilation approach to predict population dynamics during epithelial-mesenchymal transition

Epithelial-mesenchymal transition (EMT) is a biological process that plays a central role in embryonic development, tissue regeneration, and cancer metastasis. Transforming growth factor-β (TGFβ) is a potent inducer of this cellular transition, comprising transitions from an epithelial state to partial or hybrid EMT state(s), to a mesenchymal state. Recent experimental studies have shown that, within a population of epithelial cells, heterogeneous phenotypical profiles arise in response to different time- and TGFβ dose-dependent stimuli. This offers a challenge for computational models, as most model parameters are generally obtained to represent typical cell responses, not necessarily specific responses nor to capture population variability. In this study, we applied a data-assimilation approach that combines limited noisy observations with predictions from a computational model, paired with parameter estimation. Synthetic experiments mimic the biological heterogeneity in cell states that is observed in epithelial cell populations by generating a large population of model parameter sets. Analysis of the parameters for virtual epithelial cells with biologically significant characteristics (e.g., EMT prone or resistant) illustrates that these sub-populations have identifiable critical model parameters. We perform a series of in silico experiments in which a forecasting system reconstructs the EMT dynamics of each virtual cell within a heterogeneous population exposed to time-dependent exogenous TGFβ dose and either an EMT-suppressing or EMT-promoting perturbation. We find that estimating population-specific critical parameters significantly improved the prediction accuracy of cell responses. Thus, with appropriate protocol design, we demonstrate that a data-assimilation approach successfully reconstructs and predicts the dynamics of a heterogeneous virtual epithelial cell population in the presence of physiological model error and parameter uncertainty.

Comparative single-cell transcriptomes of dose and time dependent epithelial–mesenchymal spectrums

Epithelial–mesenchymal transition (EMT) is a cellular process involved in development and disease progression. Intermediate EMT states were observed in tumors and fibrotic tissues, but previous in vitro studies focused on time-dependent responses with single doses of signals; it was unclear whether single-cell transcriptomes support stable intermediates observed in diseases. Here, we performed single-cell RNA-sequencing with human mammary epithelial cells treated with multiple doses of TGF-β. We found that dose-dependent EMT harbors multiple intermediate states at nearly steady state. Comparisons of dose- and time-dependent EMT transcriptomes revealed that the dose-dependent data enable higher sensitivity to detect genes associated with EMT. We identified cell clusters unique to time-dependent EMT, reflecting cells en route to stable states. Combining dose- and time-dependent cell clusters gave rise to accurate prognosis for cancer patients. Our transcriptomic data and analyses uncover a stable EMT continuum at the single-cell resolution, and complementary information of two types of single-cell experiments.

Epithelial-Mesenchymal Plasticity in Tumor Immune Evasion

Epithelial-mesenchymal transition (EMT) is a fundamental process that occurs during embryogenesis and tissue repair. However, EMT can be hijacked by malignant cells, where it may promote immune evasion and metastasis. Classically considered a dichotomous transition, EMT in cancer has recently been considered a plastic process whereby malignant cells display and interconvert among hybrid epithelial/mesenchymal (E/M) states. Epithelial-mesenchymal plasticity (EMP) and associated hybrid E/M states are divergent from classical EMT, with unique immunomodulatory effects. Here, we review recent insights into the EMP-immune cross-talk, highlighting possible mechanisms of immune evasion conferred by hybrid E/M states and roles of immune cells in EMP.

Sirtuins and Hypoxia in EMT Control

Epithelial–mesenchymal transition (EMT), a physiological process during embryogenesis, can become pathological in the presence of different driving forces. Reduced oxygen tension or hypoxia is one of these forces, triggering a large number of molecular pathways with aberrant EMT induction, resulting in cancer and fibrosis onset. Both hypoxia-induced factors, HIF-1α and HIF-2α, act as master transcription factors implicated in EMT. On the other hand, hypoxia-dependent HIF-independent EMT has also been described. Recently, a new class of seven proteins with deacylase activity, called sirtuins, have been implicated in the control of both hypoxia responses, HIF-1α and HIF-2α activation, as well as EMT induction. Intriguingly, different sirtuins have different effects on hypoxia and EMT, acting as either activators or inhibitors, depending on the tissue and cell type. Interestingly, sirtuins and HIF can be activated or inhibited with natural or synthetic molecules. Moreover, recent studies have shown that these natural or synthetic molecules can be better conveyed using nanoparticles, representing a valid strategy for EMT modulation. The following review, by detailing the aspects listed above, summarizes the interplay between hypoxia, sirtuins, and EMT, as well as the possible strategies to modulate them by using a nanoparticle-based approach.

Revisiting the miR-200 Family: A Clan of Five Siblings with Essential Roles in Development and Disease

Over two decades of studies on small noncoding RNA molecules illustrate the significance of microRNAs (miRNAs/miRs) in controlling multiple physiological and pathological functions through post-transcriptional and spatiotemporal gene expression. Among the plethora of miRs that are essential during animal embryonic development, in this review, we elaborate the indispensable role of the miR-200 family (comprising miR-200a, -200b, 200c, -141, and -429) in governing the cellular functions associated with epithelial homeostasis, such as epithelial differentiation and neurogenesis. Additionally, in pathological contexts, miR-200 family members are primarily involved in tumor-suppressive roles, including the reversal of the cancer-associated epithelial–mesenchymal transition dedifferentiation process, and are dysregulated during organ fibrosis. Moreover, recent eminent studies have elucidated the crucial roles of miR-200s in the pathophysiology of multiple neurodegenerative diseases and tissue fibrosis. Lastly, we summarize the key studies that have recognized the potential use of miR-200 members as biomarkers for the diagnosis and prognosis of cancers, elaborating the application of these small biomolecules in aiding early cancer detection and intervention.

Distinct cytosine modification profiles define epithelial-to-mesenchymal cell-state transitions

Background: Epithelial-to-mesenchymal transition (EMT) is an early step in the invasion-metastasis cascade, involving progression through intermediate cell states. Due to challenges with isolating intermediate cell states, genome-wide cytosine modifications that define transition are not completely understood. Methods: The authors measured multiple DNA cytosine modification marks and chromatin accessibility across clonal populations residing in specific EMT states. Results: Clones exhibiting more intermediate EMT phenotypes demonstrated increased 5-hydroxymethylcytosine and decreased 5-methylcytosine. Open chromatin regions containing increased 5-hydroxymethylcytosine CpG loci were enriched in EMT transcription factor motifs and were associated with Rho GTPases. Conclusion: The results indicate the importance of both distinct and shared epigenetic profiles associated with EMT processes that may be targeted to prevent EMT progression.

Capybara: A computational tool to measure cell identity and fate transitions

Measuring cell identity in development, disease, and reprogramming is challenging as cell types and states are in continual transition. Here, we present Capybara, a computational tool to classify discrete cell identity and intermediate "hybrid" cell states, supporting a metric to quantify cell fate transition dynamics. We validate hybrid cells using experimental lineage tracing data to demonstrate the multi-lineage potential of these intermediate cell states. We apply Capybara to diagnose shortcomings in several cell engineering protocols, identifying hybrid states in cardiac reprogramming and off-target identities in motor neuron programming, which we alleviate by adding exogenous signaling factors. Further, we establish a putative in vivo correlate for induced endoderm progenitors. Together, these results showcase the utility of Capybara to dissect cell identity and fate transitions, prioritizing interventions to enhance the efficiency and fidelity of stem cell engineering.

Nonmodular oscillator and switch based on RNA decay drive regeneration of multimodal gene expression

Periodic gene expression dynamics are key to cell and organism physiology. Studies of oscillatory expression have focused on networks with intuitive regulatory negative feedback loops, leaving unknown whether other common biochemical reactions can produce oscillations. Oscillation and noise have been proposed to support mammalian progenitor cells’ capacity to restore heterogenous, multimodal expression from extreme subpopulations, but underlying networks and specific roles of noise remained elusive. We use mass-action-based models to show that regulated RNA degradation involving as few as two RNA species—applicable to nearly half of human protein-coding genes—can generate sustained oscillations without explicit feedback. Diverging oscillation periods synergize with noise to robustly restore cell populations’ bimodal expression on timescales of days. The global bifurcation organizing this divergence relies on an oscillator and bistable switch which cannot be decomposed into two structural modules. Our work reveals surprisingly rich dynamics of post-transcriptional reactions and a potentially widespread mechanism underlying development, tissue regeneration, and cancer cell heterogeneity.

Heterogeneity and dynamic of EMT through the plasticity of ribosome and mRNA translation

Growing evidence exposes translation and its translational machinery as key players in establishing and maintaining physiological and pathological biological processes. Examining translation may not only provide new biological insight but also identify novel innovative therapeutic targets in several fields of biology, including that of epithelial-to-mesenchymal transition (EMT). EMT is currently considered as a dynamic and reversible transdifferentiation process sustaining the transition from an epithelial to mesenchymal phenotype, known to be mainly driven by transcriptional reprogramming. However, it seems that the characterization of EMT plasticity is challenging, relying exclusively on transcriptomic and epigenetic approaches. Indeed, heterogeneity in EMT programs was reported to depend on the biological context. Here, by reviewing the involvement of translational control, translational machinery and ribosome biogenesis characterizing the different types of EMT, from embryonic and adult physiological to pathological contexts, we discuss the added value of integrating translational control and its machinery to depict the heterogeneity and dynamics of EMT programs.

Transcriptional and post-transcriptional control of epithelial-mesenchymal plasticity: why so many regulators?

The dynamic transition between epithelial-like and mesenchymal-like cell states has been a focus for extensive investigation for decades, reflective of the importance of Epithelial-Mesenchymal Transition (EMT) through development, in the adult, and the contributing role EMT has to pathologies including metastasis and fibrosis. Not surprisingly, regulation of the complex genetic networks that underlie EMT have been attributed to multiple transcription factors and microRNAs. What is surprising, however, are the sheer number of different regulators (hundreds of transcription factors and microRNAs) for which critical roles have been described. This review seeks not to collate these studies, but to provide a perspective on the fundamental question of whether it is really feasible that so many regulators play important roles and if so, what does this tell us about EMT and more generally, the genetic machinery that controls complex biological processes.

Ferritinophagy-Mediated Ferroptosis and Activation of Keap1/Nrf2/HO-1 Pathway Were Conducive to EMT Inhibition of Gastric Cancer Cells in Action of 2,2'-Di-pyridineketone Hydrazone Dithiocarbamate Butyric Acid Ester

In metastasis of cancer cells, the epithelial-mesenchymal transition (EMT) is prerequired. Ferroptosis is an iron-mediated cellular death process, but whether it involves EMT regulation remains elusive. In addition, how stress responders (Nrf2) respond to the redox alteration and cross-talking between them needs to be determined. Our data revealed that DpdtbA (2,2'-di-pyridineketone hydrazone dithiocarbamate butyric acid ester) resisted TGF-β1-induced EMT in gastric cancer lines (SGC-7901 and MGC-823) through ferritinophagy-mediated ROS production. Furthermore, the depletion of Gpx4 and xCT as well as enhanced lipid peroxidation indicated that DpdtbA acted as Erastin did in ferroptosis induction, which thus provided chance to explore the causal relationship between ferroptosis and EMT. Our data illustrated that ferritinophagy-mediated ferroptosis promoted the EMT inhibition. In addition, activated Nrf2 involved the regulation on both ferroptosis and EMT in response to the alteration in the cellular redox environment. In brief, ferritinophagy-mediated ferroptosis and activation of the Keap1/Nrf2/HO-1 pathway were conducive to the EMT inhibition.

Dual role of Ovol2 on the germ cell lineage segregation during gastrulation in mouse embryogenesis

In mammals, primordial germ cells (PGCs), the origin of the germ line, are specified from the epiblast at the posterior region where gastrulation simultaneously occurs, yet the functional relationship between PGC specification and gastrulation remains unclear. Here, we show that OVOL2, a transcription factor conserved across the animal kingdom, balances these major developmental processes by repressing the epithelial-to-mesenchymal transition (EMT) that drives gastrulation and the upregulation of genes associated with PGC specification. Ovol2a, a splice variant encoding a repressor domain, directly regulates EMT-related genes and, consequently, induces re-acquisition of potential pluripotency during PGC specification, whereas Ovol2b, another splice variant missing the repressor domain, directly upregulates genes associated with PGC specification. Taken together, these results elucidate the molecular mechanism underlying allocation of the germ line among epiblast cells differentiating into somatic cells through gastrulation.

This article has an associated ‘The people behind the papers’ interview.

Regulation of epithelial-to-mesenchymal transition in hypoxia by the HIF-1α network

Epithelial-to-mesenchymal transition (EMT) plays a significant role in cancer metastasis. A series of models have focused on EMT regulation by TGF-β network. However, how EMT is regulated under hypoxia is less understood. We developed a model of HIF-1α network to explore the potential link between EMT and the network topology. Our results revealed that three positive feedback loops, composed of HIF-1α and its three targets SNAIL, TWIST, and miR-210, should be sequentially activated to induce EMT under aggravating hypoxia. We suggested that the number of the positive feedback loops is critical for determining the number of stable states in EMT. Our work may advance the understanding of the significance of network topology in the regulation of EMT.

NRF2-dependent Epigenetic Regulation can Promote the Hybrid Epithelial/Mesenchymal Phenotype

The epithelial-mesenchymal transition (EMT) is a cellular process critical for wound healing, cancer metastasis and embryonic development. Recent efforts have identified the role of hybrid epithelial/mesenchymal states, having both epithelial and mesehncymal traits, in enabling cancer metastasis and resistance to various therapies. Also, previous work has suggested that NRF2 can act as phenotypic stability factor to help stablize such hybrid states. Here, we incorporate a phenomenological epigenetic feedback effect into our previous computational model for EMT signaling. We show that this type of feedback can stabilize the hybrid state as compared to the fully mesenchymal phenotype if NRF2 can influence SNAIL at an epigenetic level, as this link makes transitions out of hybrid state more difficult. However, epigenetic regulation on other NRF2-related links do not significantly change the EMT dynamics. Finally, we considered possible cell division effects in our epigenetic regulation model, and our results indicate that the degree of epigenetic inheritance does not appear to be a critical factor for the hybrid E/M state stabilizing behavior of NRF2.

Coordinate control of basal epithelial cell fate and stem cell maintenance by core EMT transcription factor Zeb1

Maintenance of undifferentiated, long-lived, and often quiescent stem cells in the basal compartment is important for homeostasis and regeneration of multiple epithelial tissues, but the molecular mechanisms that coordinately control basal cell fate and stem cell quiescence are elusive. Here, we report an epithelium-intrinsic requirement for Zeb1, a core transcriptional inducer of epithelial-to-mesenchymal transition, for mammary epithelial ductal side branching and for basal cell regenerative capacity. Our findings uncover an evolutionarily conserved role of Zeb1 in promoting basal cell fate over luminal differentiation. We show that Zeb1 loss results in increased basal cell proliferation at the expense of quiescence and self-renewal. Moreover, Zeb1 cooperates with YAP to activate Axin2 expression, and inhibition of Wnt signaling partially restores stem cell function to Zeb1-deficient basal cells. Thus, Zeb1 is a transcriptional regulator that maintains both basal cell fate and stem cell quiescence, and it functions in part through suppressing Wnt signaling.

Interconnected high-dimensional landscapes of epithelial-mesenchymal plasticity and stemness in cancer

Establishing macrometastases at distant organs is a highly challenging process for cancer cells, with extremely high attrition rates. A very small percentage of disseminated cells have the ability to dynamically adapt to their changing micro-environments through reversibly switching to another phenotype, aiding metastasis. Such plasticity can be exhibited along one or more axes-epithelial-mesenchymal plasticity (EMP) and cancer stem cells (CSCs) being the two most studied, and often tacitly assumed to be synonymous. Here, we review the emerging concepts related to EMP and CSCs across multiple cancers. Both processes are multi-dimensional in nature; for instance, EMP can be defined on morphological, molecular and functional changes, which may or may not be synchronized. Similarly, self-renewal, multi-lineage potential, and resistance to anoikis and/or therapy may not all occur simultaneously in CSCs. Thus, understanding the complexity in defining EMP and CSCs is essential if we are to understand their contribution to cancer metastasis. This will require a more comprehensive understanding of the non-linearity of these processes. These processes are dynamic, reversible, and semi-independent in nature; cells traverse the inter-connected high-dimensional EMP and CSC landscapes in diverse paths, each of which may exhibit a distinct EMP-CSC coupling. Our proposed model offers a potential unifying framework for elucidating the coupled decision-making along these dimensions and highlights a key set of open questions to be answered.

Immunosuppressive Traits of the Hybrid Epithelial/Mesenchymal Phenotype

Recent preclinical and clinical data suggests enhanced metastatic fitness of hybrid epithelial/mesenchymal (E/M) phenotypes, but mechanistic details regarding their survival strategies during metastasis remain unclear. Here, we investigate immune-evasive strategies of hybrid E/M states. We construct and simulate the dynamics of a minimalistic regulatory network encompassing the known associations among regulators of EMT (epithelial-mesenchymal transition) and PD-L1, an established immune-suppressor. Our simulations for the network consisting of SLUG, ZEB1, miR-200, CDH1 and PD-L1, integrated with single-cell and bulk RNA-seq data analysis, elucidate that hybrid E/M cells can have high levels of PD-L1, similar to those seen in cells with a full EMT phenotype, thus obviating the need for cancer cells to undergo a full EMT to be immune-evasive. Specifically, in breast cancer, we show the co-existence of hybrid E/M phenotypes, enhanced resistance to anti-estrogen therapy and increased PD-L1 levels. Our results underscore how the emergent dynamics of interconnected regulatory networks can coordinate different axes of cellular fitness during metastasis.

Unraveling the stochastic transition mechanism between oscillation states by landscape and minimum action path theory

Cell fate transitions have been studied from various perspectives, such as the transition between stable states, or the transition between stable states and oscillation states. However, there is a lack...