Comprehensive molecular interaction map of TGFβ induced epithelial to mesenchymal transition in breast cancer

Breast cancer is one of the prevailing cancers globally, with a high mortality rate. Metastatic breast cancer (MBC) is an advanced stage of cancer, characterised by a highly nonlinear, heterogeneous process involving numerous singling pathways and regulatory interactions. Epithelial-mesenchymal transition (EMT) emerges as a key mechanism exploited by cancer cells. Transforming Growth Factor-β (TGFβ)-dependent signalling is attributed to promote EMT in advanced stages of breast cancer. A comprehensive regulatory map of TGFβ induced EMT was developed through an extensive literature survey. The network assembled comprises of 312 distinct species (proteins, genes, RNAs, complexes), and 426 reactions (state transitions, nuclear translocations, complex associations, and dissociations). The map was developed by following Systems Biology Graphical Notation (SBGN) using Cell Designer and made publicly available using MINERVA ( http://35.174.227.105:8080/minerva/?id=Metastatic_Breast_Cancer_1 ). While the complete molecular mechanism of MBC is still not known, the map captures the elaborate signalling interplay of TGFβ induced EMT-promoting MBC. Subsequently, the disease map assembled was translated into a Boolean model utilising CaSQ and analysed using Cell Collective. Simulations of these have captured the known experimental outcomes of TGFβ induced EMT in MBC. Hub regulators of the assembled map were identified, and their transcriptome-based analysis confirmed their role in cancer metastasis. Elaborate analysis of this map may help in gaining additional insights into the development and progression of metastatic breast cancer.

Mesenchymal stem cells derived from CHIR99021 and TGF‑β induction remained on the colicomentum and improved cardiac function of a rat model of acute myocardium infarction

Human induced pluripotent stem cells (hiPSCs) have been regarded as a potential stem cell source for cell therapy. However, the production of cells with mesenchymal potential from hiPSCs through spontaneous differentiation is time consuming and laborious. In the present study, the combined use of the GSK-3 inhibitor CHIR99021 and TGF-β was used to obtain mesenchymal stem cell (MSC)-like cells from hiPSCs. During the induction process, the transcription of epithelial-mesenchymal transition (EMT)-related genes N-cadherin and Vimentin in the transformed cells was upregulated, whereas the transcription of E-cadherin and pluripotency-related transcription factors SOX2, OCT4 and NANOG did not change significantly. This indicated that whilst cells were pluripotent, EMT was initiated by the upregulation of transcription of EMT promoting genes. Both SMAD-dependent and independent signalling pathways were significantly activated by the combined induction treatment compared with the single factor induction. The hiPSC-derived MSC-like cells (hiPSC-MSCs) expressed MSC-related markers and acquired osteogenic, chondrogenic and adipogenic differentiation potentials. After being injected into the peritoneal cavity of rats, the hiPSC-MSCs secreted angiogenic and immune-regulatory factors and remained on the colicomentum for 3 weeks. Within an 11-week period, four intraperitoneal hiPSC-MSC injections (1x107 cells/injection) into acute myocardial infarction (AMI) model rats significantly increased the left ventricular ejection fraction, left ventricular fractional shortening and angiogenesis and significantly reduced scar size and the extent of apoptosis in the infarcted area compared with that of the control PBS injection. Symptoms of hiPSC-MSC-induced immune reaction or tumour formation were not observed over the course of the experiment in the hiSPC-MSC treated rats. In conclusion, the CHIR99021 and TGF-β combined induction was a rapid and effective method to obtain MSC-like cells from hiPSCs and multiple high dose intraperitoneal injections of hiPSC-derived MSCs were safe and effective at restoring cardiac function in an AMI rat model.

DNA-PKcs/AKT1 inhibits epithelial-mesenchymal transition during radiation-induced pulmonary fibrosis by inducing ubiquitination and degradation of Twist1

INTRODUCTION

Radiation-induced pulmonary fibrosis (RIPF) is a chronic, progressive, irreversible lung interstitial disease that develops after radiotherapy. Although several previous studies have focused on the mechanism of epithelial-mesenchymal transition (EMT) in lung epithelial cells, the essential factors involved in this process remain poorly understood. The DNA-dependent protein kinase catalytic subunit (DNA-PKcs) exhibits strong repair capacity when cells undergo radiation-induced damage; whether DNA-PKcs regulates EMT during RIPF remains unclear.

OBJECTIVES

To investigate the role and molecular mechanism of DNA-PKcs in RIPF and provide an important theoretical basis for utilising DNA-PKcs-targeted drugs for preventing RIPF.

METHODS

DNA-PKcs knockout (DPK-/-) mice were generated via the Cas9/sgRNA technique and subjected to whole chest ionizing radiation (IR) at a 20 Gy dose. Before whole chest IR, the mice were intragastrically administered the DNA-PKcs-targeted drug VND3207. Lung tissues were collected at 1 and 5 months after IR.

RESULTS

The expression of DNA-PKcs is low in pulmonary fibrosis (PF) patients. DNA-PKcs deficiency significantly exacerbated RIPF by promoting EMT in lung epithelial cells. Mechanistically, DNA-PKcs deletion by shRNA or inhibitor NU7441 maintained the protein stability of Twist1. Furthermore, AKT1 mediated the interaction between DNA-PKcs and Twist1. High Twist1 expression and EMT-associated changes caused by DNA-PKcs deletion were blocked by insulin-like growth factor-1 (IGF-1), an AKT1 agonist. The radioprotective drug VND3207 prevented IR-induced EMT and alleviated RIPF in mice by stimulating the kinase activity of DNA-PKcs.

CONCLUSION

Our study clarified the critical role and mechanism of DNA-PKcs in RIPF and showed that it could be a potential target for preventing RIPF.

Investigating the Role of SNAI1 and ZEB1 Expression in Prostate Cancer Progression and Immune Modulation of the Tumor Microenvironment

Prostate cancer (PCa) is an immunologically cold tumor and the molecular processes that underlie this behavior are poorly understood. In this study, we investigated a primary cohort of intermediate-risk PCa (n = 51) using two NanoString profiling panels designed to study cancer progression and immune response. We identified differentially expressed genes (DEGs) and pathways associated with biochemical recurrence (BCR) and clinical risk. Confirmatory analysis was performed using the TCGA-PRAD cohort. Noteworthy DEGs included collagens such as COL1A1, COL1A2, and COL3A1. Changes in the distribution of collagens may influence the immune activity in the tumor microenvironment (TME). In addition, immune-related DEGs such as THY1, IRF5, and HLA-DRA were also identified. Enrichment analysis highlighted pathways such as those associated with angiogenesis, TGF-beta, UV response, and EMT. Among the 39 significant DEGs, 11 (28%) were identified as EMT target genes for ZEB1 using the Harmonizome database. Elevated ZEB1 expression correlated with reduced BCR risk. Immune landscape analysis revealed that ZEB1 was associated with increased immunosuppressive cell types in the TME, such as naïve B cells and M2 macrophages. Increased expression of both ZEB1 and SNAI1 was associated with elevated immune checkpoint expression. In the future, modulation of EMT could be beneficial for overcoming immunotherapy resistance in a cold tumor, such as PCa.

Cancer spreading patterns based on epithelial-mesenchymal plasticity

Introduction: Metastasis is a major cause of cancer-related deaths, underscoring the necessity to discern the rules and patterns of cancer cell spreading. Epithelial-mesenchymal plasticity contributes to cancer aggressiveness and metastasis. Despite establishing key determinants of cancer aggressiveness and metastatic ability, a comprehensive understanding of the underlying mechanism is unknown. We aimed to propose a classification system for cancer cells based on epithelial-mesenchymal plasticity, focusing on hysteresis of the epithelial-mesenchymal transition and the hybrid epithelial/mesenchymal phenotype. Methods: We extensively reviewed the concept of epithelial-mesenchymal plasticity, specifically considering the hysteresis of the epithelial-mesenchymal transition and the hybrid epithelial/mesenchymal phenotype. Results: In this review and hypothesis article, based on epithelial-mesenchymal plasticity, especially the hysteresis of epithelial-mesenchymal transition and the hybrid epithelial/mesenchymal phenotype, we proposed a classification of cancer cells, indicating that cancer cells with epithelial-mesenchymal plasticity potential could be classified into four types: irreversible hysteresis, weak hysteresis, strong hysteresis, and hybrid epithelial/mesenchymal phenotype. These four types of cancer cells had varied biology, spreading features, and prognoses. Discussion: Our results highlight that the proposed classification system offers insights into the diverse behaviors of cancer cells, providing implications for cancer aggressiveness and metastasis.

Escargot a Snail superfamily member and its multiple roles in Drosophila melanogaster development

The Snail superfamily of transcription factors plays a crucial role in metazoan development; one of the most important vertebrate members of this family is Snai1 which is orthologous to the Drosophila melanogaster esg gene. This review offers a comprehensive examination of the roles of the esg gene in Drosophila development, covering its expression pattern and downstream targets, and draws parallels between the vertebrate Snai1 family proteins on controlling the epithelial-to-mesenchymal transition and esg. This gene regulates stemness, ploidy, and pluripontency. esg is expressed in various tissues during development, including the gut, imaginal discs, and neuroblasts. The functions of the esg include the suppression of differentiation in intestinal stem cells and the preservation of diploidy in imaginal cells. In the nervous system development, esg expression also inhibits neuroblast differentiation, thus regulating the number of neurons and the moment in development of neuronal differentiation. Loss of esg function results in diverse developmental defects, including defects in intestinal stem cell maintenance and differentiation, and alters imaginal disc and nervous system development. Expression levels of esg also play a role in regulating longevity and metabolism in adult stages. This review provides an overview of the current understanding of esg's developmental role, emphasizing cellular and tissue effects that arise from its loss of function. The insights gained may contribute to a better understanding of evolutionary conserved developmental mechanisms and certain metabolic diseases.

CD276 promotes epithelial-mesenchymal transition in esophageal squamous cell carcinoma through the TGF-β/SMAD signaling

OBJECTIVE

Aberrant expression of CD276 has been reported in malignant tumors. However, the exact role and mechanisms of CD276 influence the progression of esophageal squamous cell carcinoma (ESCC) still need to be understood.

METHODS

Bioinformatics analysis of data from The Cancer Genome Atlas and Gene Expression Omnibus databases, along with immunohistochemistry staining, was used to explore the expression patterns of CD276 in ESCC. Cell counting kit-8 and Transwell assays were employed to evaluate the effects of CD276 expression on tumor cell proliferation and motility. Western blotting and Transwell assays were used to explore the potential pathways through which CD276 mediates the progression of ESCC. Moreover, the in vivo role of CD276 in tumor progression was investigated by establishing a lung metastasis mouse model.

RESULTS

A significant upregulation of CD276 was observed in ESCC tissues compared to adjacent tissues. The inhibition of CD276 had no evident impact on ESCC cell proliferation but notably hindered their migratory and invasive properties and the expression of epithelial-mesenchymal transition (EMT) markers. Inversely, overexpressing CD276 led to an upregulation of EMT markers, underscoring the capacity of CD276 to amplify the motility of ESCC cells. Furthermore, CD276 was found to enhance the migratory and invasive abilities of ESCC cells by activating the TGF-β/SMAD signaling but not the PI3K/AKT pathway. In vivo studies demonstrated that CD276 facilitates pulmonary metastasis.

CONCLUSION

CD276 is significant upregulation in ESCC tissues and facilitates the EMT process in ESCC cells via the TGF-β/SMAD signaling, thus promoting the progression of ESCC.

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.

Tipping points in epithelial-mesenchymal lineages from single-cell transcriptomics data

Understanding cell fate decision-making during complex biological processes is an open challenge that is now aided by high-resolution single-cell sequencing technologies. Specifically, it remains challenging to identify and characterize transition states corresponding to "tipping points" whereby cells commit to new cell states. Here, we present a computational method that takes advantage of single-cell transcriptomics data to infer the stability and gene regulatory networks (GRNs) along cell lineages. Our method uses the unspliced and spliced counts from single-cell RNA sequencing data and cell ordering along lineage trajectories to train an RNA splicing multivariate model, from which cell-state stability along the lineage is inferred based on spectral analysis of the model's Jacobian matrix. Moreover, the model infers the RNA cross-species interactions resulting in GRNs and their variation along the cell lineage. When applied to epithelial-mesenchymal transition in ovarian and lung cancer-derived cell lines, our model predicts a saddle-node transition between the epithelial and mesenchymal states passing through an unstable, intermediate cell state. Furthermore, we show that the underlying GRN controlling epithelial-mesenchymal transition rearranges during the transition, resulting in denser and less modular networks in the intermediate state. Overall, our method represents a flexible tool to study cell lineages with a combination of theory-driven modeling and single-cell transcriptomics data.

High-throughput and proteome-wide discovery of endogenous biomolecular condensates

Phase separation inside mammalian cells regulates the formation of the biomolecular condensates that are related to gene expression, signalling, development and disease. However, a large population of endogenous condensates and their candidate phase-separating proteins have yet to be discovered in a quantitative and high-throughput manner. Here we demonstrate that endogenously expressed biomolecular condensates can be identified across a cell's proteome by sorting proteins across varying oligomeric states. We employ volumetric compression to modulate the concentrations of intracellular proteins and the degree of crowdedness, which are physical regulators of cellular biomolecular condensates. The changes in degree of the partition of proteins into condensates or phase separation led to varying oligomeric states of the proteins, which can be detected by coupling density gradient ultracentrifugation and quantitative mass spectrometry. In total, we identified 1,518 endogenous condensate proteins, of which 538 have not been reported before. Furthermore, we demonstrate that our strategy can identify condensate proteins that respond to specific biological processes.

Increased heterogeneity in expression of genes associated with cancer progression and drug resistance

Fluctuations in the number of regulatory molecules and differences in timings of molecular events can generate variation in gene expression among genetically identical cells in the same environmental condition. This variation, termed as expression noise, can create differences in metabolic state and cellular functions, leading to phenotypic heterogeneity. Expression noise and phenotypic heterogeneity have been recognized as important contributors to intra-tumor heterogeneity, and have been associated with cancer growth, progression, and therapy resistance. However, how expression noise changes with cancer progression in actual cancer patients has remained poorly explored. Such an analysis, through identification of genes with increasing expression noise, can provide valuable insights into generation of intra-tumor heterogeneity, and could have important implications for understanding immune-suppression, drug tolerance and therapy resistance. In this work, we performed a genome-wide identification of changes in gene expression noise with cancer progression using single-cell RNA-seq data of lung adenocarcinoma patients at different stages of cancer. We identified 37 genes in epithelial cells that showed an increasing noise trend with cancer progression, many of which were also associated with cancer growth, EMT and therapy resistance. We found that expression of several of these genes was positively associated with expression of mitochondrial genes, suggesting an important role of mitochondria in generation of heterogeneity. In addition, we uncovered substantial differences in sample-specific noise profiles which could have implications for personalized prognosis and treatment.

Cooperative pro-tumorigenic adaptation to oncogenic RAS through epithelial-to-mesenchymal plasticity

In breast cancers, aberrant activation of the RAS/MAPK pathway is strongly associated with mesenchymal features and stemness traits, suggesting an interplay between this mitogenic signaling pathway and epithelial-to-mesenchymal plasticity (EMP). By using inducible models of human mammary epithelial cells, we demonstrate herein that the oncogenic activation of RAS promotes ZEB1-dependent EMP, which is necessary for malignant transformation. Notably, EMP is triggered by the secretion of pro-inflammatory cytokines from neighboring RAS-activated senescent cells, with a prominent role for IL-6 and IL-1α. Our data contrast with the common view of cellular senescence as a tumor-suppressive mechanism and EMP as a process promoting late stages of tumor progression in response to signals from the tumor microenvironment. We highlighted here a pro-tumorigenic cooperation of RAS-activated mammary epithelial cells, which leverages on oncogene-induced senescence and EMP to trigger cellular reprogramming and malignant transformation.

Cationic solid lipid nanoparticles (SLN) complexed with plasmid DNA enhance prostate cancer cells (PC-3) migration

Nanotechnology applications in biomedicine have increased in recent decades, primarily as therapeutic agents, drugs, and gene delivery systems. Among the nanoparticles used in medicine, we highlight cationic solid lipid nanoparticles (SLN). Given their nontoxic properties, much research has focused on the beneficial effects of SLN for drug or gene delivery system. However, little attention has been paid to the adverse impacts of SLN on the cellular environment, particularly their influence on intracellular signaling pathways. In this work, we investigate the effects triggered by cationic SLN on human prostate non-tumor cells (PNT1A) and tumor cells (PC-3). Our results demonstrate that cationic SLN enhances the migration of PC-3 prostate cancer cells but not PNT1A non-tumor prostate cells, an unexpected and unprecedented development. Furthermore, we observed that the enhanced cell migration velocity is a concentration-dependent and nanoparticle-dependent effect, and not related to any individual nanoparticle component. Moreover, cationic SLN increased vimentin expression (p < 0.05) but SLN did not affect Smad2 nuclear translocation. Meanwhile, EMT-related (epithelial-to-mesenchymal transition) proteins, such as ZEB1, underwent nuclear translocation when treated with cationic SLN, thereby affecting PC-3 cell motility through ZEB1 and vimentin modulation. From a therapeutic perspective, cationic SLN could potentially worsen a patient's condition if these results were reproduced in vivo. Understanding the in vitro molecular mechanisms triggered by nanomaterials and their implications for cell function is crucial for defining their safe and effective use.

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.

Single-cell morphodynamical trajectories enable prediction of gene expression accompanying cell state change

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.

Endothelial Protein C Receptor and 3K3A-Activated Protein C Protect Mice from Allergic Contact Dermatitis in a Contact Hypersensitivity Model

Endothelial protein C receptor (EPCR) is a receptor for the natural anti-coagulant activated protein C (aPC). It mediates the anti-inflammatory and barrier-protective functions of aPC through the cleavage of protease-activated receptor (PAR)1/2. Allergic contact dermatitis is a common skin disease characterized by inflammation and defective skin barrier. This study investigated the effect of EPCR and 3K3A-aPC on allergic contact dermatitis using a contact hypersensitivity (CHS) model. CHS was induced using 1-Fluoro-2,4-dinitrobenzene in EPCR-deficient (KO) and matched wild-type mice and mice treated with 3K3A-aPC, a mutant form of aPC with diminished anti-coagulant activity. Changes in clinical and histological features, cytokines, and immune cells were examined. EPCRKO mice displayed more severe CHS, with increased immune cell infiltration in the skin and higher levels of inflammatory cytokines and IgE than wild-type mice. EPCR, aPC, and PAR1/2 were expressed by the skin epidermis, with EPCR presenting almost exclusively in the basal layer. EPCRKO increased the epidermal expression of aPC and PAR1, whereas in CHS, their expression was reduced compared to wild-type mice. 3K3A-aPC reduced CHS severity in wild-type and EPCRKO mice by suppressing immune cell infiltration/activation and inflammatory cytokines. In summary, EPCRKO exacerbated CHS, whereas 3K3A-aPC could reduce the severity of CHS in both EPCRKO and wild-type mice.

Dual role of CASP8AP2/FLASH in regulating epithelial-to-mesenchymal transition plasticity (EMP)

BACKGROUND

Epithelial-to-mesenchymal transition (EMT) is a developmental program that consists of the loss of epithelial features concomitant with the acquisition of mesenchymal features. Activation of EMT in cancer facilitates the acquisition of aggressive traits and cancer invasion. EMT plasticity (EMP), the dynamic transition between multiple hybrid states in which cancer cells display both epithelial and mesenchymal markers, confers survival advantages for cancer cells in constantly changing environments during metastasis.

METHODS

RNAseq analysis was performed to assess genome-wide transcriptional changes in cancer cells depleted for histone regulators FLASH, NPAT, and SLBP. Quantitative PCR and Western blot were used for the detection of mRNA and protein levels. Computational analysis was performed on distinct sets of genes to determine the epithelial and mesenchymal score in cancer cells and to correlate FLASH expression with EMT markers in the CCLE collection.

RESULTS

We demonstrate that loss of FLASH in cancer cells gives rise to a hybrid E/M phenotype with high epithelial scores even in the presence of TGFβ, as determined by computational methods using expression of predetermined sets of epithelial and mesenchymal genes. Multiple genes involved in cell-cell junction formation are similarly specifically upregulated in FLASH-depleted cells, suggesting that FLASH acts as a repressor of the epithelial phenotype. Further, FLASH expression in cancer lines is inversely correlated with the epithelial score. Nonetheless, subsets of mesenchymal markers were distinctly up-regulated in FLASH, NPAT, or SLBP-depleted cells.

CONCLUSIONS

The ZEB1low/SNAILhigh/E-cadherinhigh phenotype described in FLASH-depleted cancer cells is driving a hybrid E/M phenotype in which epithelial and mesenchymal markers coexist.

Stabilizing vimentin phosphorylation inhibits stem-like cell properties and metastasis of hybrid epithelial/mesenchymal carcinomas

Epithelial-mesenchymal transition (EMT) empowers epithelial cells with mesenchymal and stem-like attributes, facilitating metastasis, a leading cause of cancer-related mortality. Hybrid epithelial-mesenchymal (E/M) cells, retaining both epithelial and mesenchymal traits, exhibit heightened metastatic potential and stemness. The mesenchymal intermediate filament, vimentin, is upregulated during EMT, enhancing the resilience and invasiveness of carcinoma cells. The phosphorylation of vimentin is critical to its structure and function. Here, we identify that stabilizing vimentin phosphorylation at serine 56 induces multinucleation, specifically in hybrid E/M cells with stemness properties but not epithelial or mesenchymal cells. Cancer stem-like cells are especially susceptible to vimentin-induced multinucleation relative to differentiated cells, leading to a reduction in self-renewal and stemness. As a result, vimentin-induced multinucleation leads to sustained inhibition of stemness properties, tumor initiation, and metastasis. These observations indicate that a single, targetable phosphorylation event in vimentin is critical for stemness and metastasis in carcinomas with hybrid E/M properties.

MicroRNA-200b-mediated reversion of a spectrum of epithelial-to-mesenchymal transition states in recessive dystrophic epidermolysis bullosa squamous cell carcinomas

BACKGROUND

Cutaneous squamous cell carcinoma (SCC) is the leading cause of death in patients with recessive dystrophic epidermolysis bullosa (RDEB). However, the survival time from first diagnosis differs between patients; some tumours spread particularly fast, while others may remain localized for years. As treatment options are limited, there is an urgent need for further insights into the pathomechanisms of RDEB tumours, to foster therapy development and support clinical decision-making.

OBJECTIVES

To investigate differences in RDEB tumours of diverging aggressiveness at the molecular and phenotypic level, with a particular focus on epithelial-to-mesenchymal (EMT) transition states and thus microRNA-200b (miR-200b) as a regulator.

METHODS

Primary RDEB-SCC keratinocyte lines were characterized with respect to their EMT state. For this purpose, cell morphology was classified and the expression of EMT markers analysed using immunofluorescence, flow cytometry, semi-quantitative reverse transcriptase polymerase chain reaction and Western blotting. The motility of RDEB-SCC cells was determined and conditioned medium of RDEB-SCC cells was used to treat endothelial cells in an angiogenesis assay. In addition, we mined previously generated microRNA (miRNA) profiling data to identify a candidate with potential therapeutic relevance and performed transient miRNA transfection studies to investigate the candidate's ability to reverse EMT characteristics.

RESULTS

We observed high variability in EMT state in the RDEB-SCC cell lines, which correlated with in situ analysis of two available patient biopsies and respective clinical disease course. Furthermore, we identified miR-200b-3p to be downregulated in RDEB-SCCs, and the extent of deregulation significantly correlated with the EMT features of the various tumour lines. miR-200b-3p was reintroduced into RDEB-SCC cell lines with pronounced EMT features, which resulted in a significant increase in epithelial characteristics, including cell morphology, EMT marker expression, migration and angiogenic potential.

CONCLUSIONS

RDEB-SCCs exist in different EMT states and the level of miR-200b is indicative of how far an RDEB-SCC has gone down the EMT path. Moreover, the reintroduction of miR-200b significantly reduced mesenchymal features.

Prostate transmembrane androgen inducible protein 1 (PMEPA1): regulation and clinical implications

Prostate transmembrane androgen inducible protein 1 (PMEPA1) can promote or inhibit prostate cancer cell growth based on the cancer cell response to the androgen receptor (AR). Further, it can be upregulated by transforming growth factor (TGF), which downregulates transforming growth factor-β (TGF-β) signaling by interfering with R-Smad phosphorylation to facilitate TGF-β receptor degradation. Studies have indicated the increased expression of PMEPA1 in some solid tumors and its functioning as a regulator of multiple signaling pathways. This review highlights the multiple potential signaling pathways associated with PMEPA1 and the role of the PMEPA1 gene in regulating prognosis, including transcriptional regulation and epithelial mesenchymal transition (EMT). Moreover, the relevant implications in and outside tumors, for example, as a biomarker and its potential functions in lysosomes have also been discussed.

TGF-β, EMT, and resistance to anti-cancer treatment

Transforming growth factor-β (TGF-β) signaling regulates cell-specific programs involved in embryonic development, wound-healing, and immune homeostasis. Yet, during tumor progression, these TGF-β-mediated programs are altered, leading to epithelial cell plasticity and a reprogramming of epithelial cells into mesenchymal lineages through epithelial-to-mesenchymal transition (EMT), a critical developmental program in morphogenesis and organogenesis. These changes, in turn, lead to enhanced carcinoma cell invasion, metastasis, immune cell differentiation, immune evasion, and chemotherapy resistance. Here, we discuss EMT as one of the critical programs associated with carcinoma cell plasticity and the influence exerted by TGF-β on carcinoma status and function. We further explore the composition of carcinoma and other cell populations within the tumor microenvironment, and consider the relevant outcomes related to the programs associated with cancer treatment resistance.

Identifying key factors in cell fate decisions by machine learning interpretable strategies

Cell fate decisions and transitions are common in almost all developmental processes. Therefore, it is important to identify the decision-making mechanisms and important individual molecules behind the fate decision processes. In this paper, we propose an interpretable strategy based on systematic perturbation, unsupervised hierarchical cluster analysis (HCA), machine learning (ML), and Shapley additive explanation (SHAP) analysis for inferring the contribution and importance of individual molecules in cell fate decision and transition processes. In order to verify feasibility of the approach, we apply it to the core epithelial to mesenchymal transition (EMT)-metastasis network. The key factors identified in EMT-metastasis are consistent with relevant experimental observations. The approach presented here can be applied to other biological networks to identify important factors related to cell fate decisions and transitions.

Cusp bifurcation in a metastatic regulatory network

Understanding the potential for cancers to metastasize is still relatively unknown. While many predictive methods may use deep learning or stochastic processes, we highlight a long standing mathematical concept that may be useful for modeling metastatic breast cancer systems. Ordinary differential equations (ODEs) can model cell state transitions by considering the pertinent environmental variables as well as the paths systems take over time. Bifurcation theory is a branch of dynamical systems which studies changes in the behavior of an ODE system while one or more parameters are varied. Many studies have applied concepts in one-parameter bifurcation theory to model biological network dynamics, and cell division. However, studies of two-parameter bifurcations are much more rare. Two-parameter bifurcations have not been studied in metastatic systems. Here we show how a specific two-parameter bifurcation phenomenon called a cusp bifurcation separates two qualitatively different metastatic cell state transitions modalities and propose a new perspective on defining such transitions based on mathematical theory. We hope the observations and verification methods detailed here may help in the understanding of metastatic potential from a basic biological perspective and in clinical settings.

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.

A logical model of Ewing sarcoma cell epithelial-to-mesenchymal transition supports the existence of hybrid cellular phenotypes

Ewing sarcoma (ES) is a highly aggressive pediatric tumor driven by the RNA-binding protein EWS (EWS)/friend leukemia integration 1 transcription factor (FLI1) chimeric transcription factor, which is involved in epithelial-mesenchymal transition (EMT). EMT stabilizes a hybrid cell state, boosting metastatic potential and drug resistance. Nevertheless, the mechanisms underlying the maintenance of this hybrid phenotype in ES remain elusive. Our study proposes a logical EMT model for ES, highlighting zinc finger E-box-binding homeobox 2 (ZEB2), miR-145, and miR-200 circuits that maintain hybrid states. The model aligns with experimental findings and reveals a previously unknown circuit supporting the mesenchymal phenotype. These insights emphasize the role of ZEB2 in the maintenance of the hybrid state in ES.

The epithelial-mesenchymal plasticity landscape: principles of design and mechanisms of regulation

Epithelial-mesenchymal plasticity (EMP) enables cells to interconvert between several states across the epithelial-mesenchymal landscape, thereby acquiring hybrid epithelial/mesenchymal phenotypic features. This plasticity is crucial for embryonic development and wound healing, but also underlies the acquisition of several malignant traits during cancer progression. Recent research using systems biology and single-cell profiling methods has provided novel insights into the main forces that shape EMP, which include the microenvironment, lineage specification and cell identity, and the genome. Additionally, key roles have emerged for hysteresis (cell memory) and cellular noise, which can drive stochastic transitions between cell states. Here, we review these forces and the distinct but interwoven layers of regulatory control that stabilize EMP states or facilitate epithelial-mesenchymal transitions (EMTs) and discuss the therapeutic potential of manipulating the EMP landscape.

Theoretical and computational tools to model multistable gene regulatory networks

The last decade has witnessed a surge of theoretical and computational models to describe the dynamics of complex gene regulatory networks, and how these interactions can give rise to multistable and heterogeneous cell populations. As the use of theoretical modeling to describe genetic and biochemical circuits becomes more widespread, theoreticians with mathematical and physical backgrounds routinely apply concepts from statistical physics, non-linear dynamics, and network theory to biological systems. This review aims at providing a clear overview of the most important methodologies applied in the field while highlighting current and future challenges. It also includes hands-on tutorials to solve and simulate some of the archetypical biological system models used in the field. Furthermore, we provide concrete examples from the existing literature for theoreticians that wish to explore this fast-developing field. Whenever possible, we highlight the similarities and differences between biochemical and regulatory networks and 'classical' systems typically studied in non-equilibrium statistical and quantum mechanics.

Tumor-mediated immunosuppression and cytokine spreading affects the relation between EMT and PD-L1 status

Epithelial-mesenchymal transition (EMT) and immune resistance mediated by Programmed Death-Ligand 1 (PD-L1) upregulation are established drivers of tumor progression. Their bi-directional crosstalk has been proposed to facilitate tumor immunoevasion, yet the impact of immunosuppression and spatial heterogeneity on the interplay between these processes remains to be characterized. Here we study the role of these factors using mathematical and spatial models. We first designed models incorporating immunosuppressive effects on T cells mediated via PD-L1 and the EMT-inducing cytokine Transforming Growth Factor beta (TGFβ). Our models predict that PD-L1-mediated immunosuppression merely reduces the difference in PD-L1 levels between EMT states, while TGFβ-mediated suppression also causes PD-L1 expression to correlate negatively with TGFβ within each EMT phenotype. We subsequently embedded the models in multi-scale spatial simulations to explicitly describe heterogeneity in cytokine levels and intratumoral heterogeneity. Our multi-scale models show that Interferon gamma (IFNγ)-induced partial EMT of a tumor cell subpopulation can provide some, albeit limited protection to bystander tumor cells. Moreover, our simulations show that the true relationship between EMT status and PD-L1 expression may be hidden at the population level, highlighting the importance of studying EMT and PD-L1 status at the single-cell level. Our findings deepen the understanding of the interactions between EMT and the immune response, which is crucial for developing novel diagnostics and therapeutics for cancer patients.

Noncoding RNA circuitry in melanoma onset, plasticity, and therapeutic response

Melanoma, the cancer of the melanocyte, is the deadliest form of skin cancer with an aggressive nature, propensity to metastasize and tendency to resist therapeutic intervention. Studies have identified that the re-emergence of developmental pathways in melanoma contributes to melanoma onset, plasticity, and therapeutic response. Notably, it is well known that noncoding RNAs play a critical role in the development and stress response of tissues. In this review, we focus on the noncoding RNAs, including microRNAs, long non-coding RNAs, circular RNAs, and other small RNAs, for their functions in developmental mechanisms and plasticity, which drive onset, progression, therapeutic response and resistance in melanoma. Going forward, elucidation of noncoding RNA-mediated mechanisms may provide insights that accelerate development of novel melanoma therapies.

Population dynamics of EMT elucidates the timing and distribution of phenotypic intra-tumoral heterogeneity

The Epithelial-to-Mesenchymal Transition (EMT) is a hallmark of cancer metastasis and morbidity. EMT is a non-binary process, and cells can be stably arrested en route to EMT in an intermediate hybrid state associated with enhanced tumor aggressiveness and worse patient outcomes. Understanding EMT progression in detail will provide fundamental insights into the mechanisms underlying metastasis. Despite increasingly available single-cell RNA sequencing (scRNA-seq) data that enable in-depth analyses of EMT at the single-cell resolution, current inferential approaches are limited to bulk microarray data. There is thus a great need for computational frameworks to systematically infer and predict the timing and distribution of EMT-related states at single-cell resolution. Here, we develop a computational framework for reliable inference and prediction of EMT-related trajectories from scRNA-seq data. Our model can be utilized across a variety of applications to predict the timing and distribution of EMT from single-cell sequencing data.

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.

Fiber density and matrix stiffness modulate distinct cell migration modes in a 3D stroma mimetic composite hydrogel

The peritumoral stroma is a complex 3D tissue that provides cells with myriad biophysical and biochemical cues. Histologic observations suggest that during metastatic spread of carcinomas, these cues influence transformed epithelial cells, prompting a diversity of migration modes spanning single cell and multicellular phenotypes. Purported consequences of these variations in tumor escape strategies include differential metastatic capability and therapy resistance. Therefore, understanding how cues from the peritumoral stromal microenvironment regulate migration mode has both prognostic and therapeutic value. Here, we utilize a synthetic stromal mimetic in which matrix fiber density and bulk hydrogel mechanics can be orthogonally tuned to investigate the contribution of these two key matrix attributes on MCF10A migration mode phenotypes, epithelial-mesenchymal transition (EMT), and invasive potential. We develop an automated computational image analysis framework to extract migratory phenotypes from fluorescent images and determine 3D migration metrics relevant to metastatic spread. Using this analysis, we find that matrix fiber density and bulk hydrogel mechanics distinctly contribute to a variety of MCF10A migration modes including amoeboid, single mesenchymal, clusters, and strands. We identify combinations of physical and soluble cues that induce a variety of migration modes originating from the same MCF10A spheroid and use these settings to examine a functional consequence of migration mode -resistance to apoptosis. We find that cells migrating as strands are more resistant to staurosporine-induced apoptosis than either disconnected clusters or individual invading cells. Improved models of the peritumoral stromal microenvironment and understanding of the relationships between matrix attributes and cell migration mode can aid ongoing efforts to identify effective cancer therapeutics that address cell plasticity-based therapy resistances. STATEMENT OF SIGNIFICANCE: Stromal extracellular matrix structure dictates both cell homeostasis and activation towards migratory phenotypes. However decoupling the effects of myriad biophysical cues has been difficult to achieve. Here, we encapsulate electrospun fiber segments within an amorphous hydrogel to create a fiber-reinforced hydrogel composite in which fiber density and hydrogel stiffness can be orthogonally tuned. Quantification of 3D cell migration reveal these two parameters uniquely contribute to a diversity of migration phenotypes spanning amoeboid, single mesenchymal, multicellular cluster, and collective strand. By tuning biophysical and biochemical cues to elicit heterogeneous migration phenotypes, we find that collective strands best resist apoptosis. This work establishes a composite approach to modulate fibrous topography and bulk hydrogel mechanics and identified biomaterial parameters to direct distinct 3D cell migration phenotypes.

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.

Klotho prevents transforming growth factor-β2-induced senescent-like morphological changes in the retinal pigment epithelium

Degenerative changes of the retinal pigment epithelium (RPE) triggered by transforming growth factor-β2 (TGF-β2) and oxidative stress play a critical role in the progression of age-related macular degeneration (AMD). The expression of α-klotho, an antiaging protein, declines with age, increasing the risk factors for age-related diseases. Here, we investigated the protective effects of soluble α-klotho on TGF-β2-induced RPE degeneration. The morphological changes induced by TGF-β2, including epithelial-mesenchymal transition (EMT), were attenuated in the mouse RPE by the intravitreal injection (IVT) of α-klotho. In ARPE19 cells, EMT and morphological alterations by TGF-β2 were attenuated by co-incubation with α-klotho. TGF-β2 decreased miR-200a accompanied by zinc finger e-box binding homeobox1 (ZEB1) upregulation and EMT, all of which were prevented by α-klotho co-treatment. Inhibitor of miR-200a mimicked TGF-β2-induced morphological changes, which were recovered by ZEP1 silencing, but not by α-klotho, implying the upstream regulation of α-klotho on miR-200a-ZEP1-EMT axis. α-Klotho inhibited receptor binding of TGF-β2, Smad2/3 phosphorylation, extracellular signal-regulated protein kinase 1/2 (ERK1/2)-a mechanistic target of rapamycin (mTOR) activation and oxidative stress via NADPH oxidase 4 (NOX4) upregulation. Furthermore, α-klotho recovered the TGF-β2-induced mitochondrial activation and superoxide generation. Interestingly, TGF-β2 upregulated α-klotho expression in the RPE cells, and genetic suppression of endogenous α-klotho aggravated TGF-β2-induced oxidative stress and EMT. Lastly, α-klotho abrogated senescence-associated signaling molecules and phenotypes induced by long-term incubation with TGF-β2. Hence, our findings indicate that the antiaging α-klotho plays a protective role against EMT and degeneration of the RPE, demonstrating the therapeutic potential for age-related retinal diseases, including the dry type of AMD.

Noncanonical Wnt signaling promotes colon tumor growth, chemoresistance and tumor fibroblast activation

Colon tumors of the mesenchymal subtype have the lowest overall survival. Snail1 is essential for the acquisition of this phenotype, characterized by increased tumor stemness and invasion, and high resistance to chemotherapy. Here, we find that Snail1 expression in colon tumor cells is dependent on an autocrine noncanonical Wnt pathway. Accordingly, depletion of Ror2, the co-receptor for noncanonical Wnts such as Wnt5a, potently decreases Snail1 expression. Wnt5a, Ror2, and Snail1 participate in a self-stimulatory feedback loop since Wnt5a increases its own synthesis in a Ror2- and Snail1-dependent fashion. This Wnt5a/Ror2/Snail1 axis controls tumor invasion, chemoresistance, and formation of tumor spheres. It also stimulates TGFβ synthesis; consequently, tumor cells expressing Snail1 are more efficient in activating cancer-associated fibroblasts than the corresponding controls. Ror2 downmodulation or inhibition of the Wnt5a pathway decreases Snail1 expression in primary colon tumor cells and their ability to form tumors and liver metastases. Finally, the expression of SNAI1, ROR2, and WNT5A correlates in human colon and other tumors. These results identify inhibition of the noncanonical Wnt pathway as a putative colon tumor therapy.

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.

RSV infection does not induce EMT

Respiratory syncytial virus (RSV) infection does not cause severe disease in most of us despite suffering from multiple RSV infections in our lives. However, infants, young children, older adults, and immunocompromised patients are unfortunately vulnerable to RSV-associated severe diseases. A recent study suggested that RSV infection causes cell expansion, resulting in bronchial wall thickening in vitro. Whether the virus-induced changes in the lung airway resemble epithelial-mesenchymal transition (EMT) is still unknown. Here, we report that RSV does not induce EMT in three different in vitro lung models: the epithelial A549 cell line, primary normal human bronchial epithelial cells, and pseudostratified airway epithelium. We found that RSV increases the cell surface area and perimeter in the infected airway epithelium, which is distinct from the effects of a potent EMT inducer, TGF-β1-driven cell elongation-indicative of cell motility. A genome-wide transcriptome analysis revealed that both RSV and TGF-β1 have distinct modulation patterns of the transcriptome, which suggests that RSV-induced changes are distinct from EMT.

The ELF3 transcription factor is associated with an epithelial phenotype and represses epithelial-mesenchymal transition

BACKGROUND

Epithelial-mesenchymal plasticity (EMP) involves bidirectional transitions between epithelial, mesenchymal and multiple intermediary hybrid epithelial/mesenchymal phenotypes. While the process of epithelial-mesenchymal transition (EMT) and its associated transcription factors are well-characterised, the transcription factors that promote mesenchymal-epithelial transition (MET) and stabilise hybrid E/M phenotypes are less well understood.

RESULTS

Here, we analyse multiple publicly-available transcriptomic datasets at bulk and single-cell level and pinpoint ELF3 as a factor that is strongly associated with an epithelial phenotype and is inhibited during EMT. Using mechanism-based mathematical modelling, we also show that ELF3 inhibits the progression of EMT. This behaviour was also observed in the presence of an EMT inducing factor WT1. Our model predicts that the MET induction capacity of ELF3 is stronger than that of KLF4, but weaker than that of GRHL2. Finally, we show that ELF3 levels correlates with worse patient survival in a subset of solid tumour types.

CONCLUSION

ELF3 is shown to be inhibited during EMT progression and is also found to inhibit the progression of complete EMT suggesting that ELF3 may be able to counteract EMT induction, including in the presence of EMT-inducing factors, such as WT1. The analysis of patient survival data indicates that the prognostic capacity of ELF3 is specific to cell-of-origin or lineage.

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.

A quantitative evaluation of topological motifs and their coupling in gene circuit state distributions

One of the major challenges in biology is to understand how gene interactions collaborate to determine overall functions of biological systems. Here, we present a new computational framework that enables systematic, high-throughput, and quantitative evaluation of how small transcriptional regulatory circuit motifs, and their coupling, contribute to functions of a dynamical biological system. We illustrate how this approach can be applied to identify four-node gene circuits, circuit motifs, and motif coupling responsible for various gene expression state distributions, including those derived from single-cell RNA sequencing data. We also identify seven major classes of four-node circuits from clustering analysis of state distributions. The method is applied to establish phenomenological models of gene circuits driving human neuron differentiation, revealing important biologically relevant regulatory interactions. Our study will shed light on a better understanding of gene regulatory mechanisms in creating and maintaining cellular states.

TGFβ1-Induced EMT in the MCF10A Mammary Epithelial Cell Line Model Is Executed Independently of SNAIL1 and ZEB1 but Relies on JUNB-Coordinated Transcriptional Regulation

Epithelial-mesenchymal transition (EMT) fosters cancer cell invasion and metastasis, the main cause of cancer-related mortality. Growing evidence that SNAIL and ZEB transcription factors, typically portrayed as master regulators of EMT, may be dispensable for this process, led us to re-investigate its mechanistic underpinnings. For this, we used an unbiased computational approach that integrated time-resolved analyses of chromatin structure and differential gene expression, to predict transcriptional regulators of TGFβ1-inducible EMT in the MCF10A mammary epithelial cell line model. Bioinformatic analyses indicated comparatively minor contributions of SNAIL proteins and ZEB1 to TGFβ1-induced EMT, whereas the AP-1 subunit JUNB was anticipated to have a much larger impact. CRISPR/Cas9-mediated loss-of-function studies confirmed that TGFβ1-induced EMT proceeded independently of SNAIL proteins and ZEB1. In contrast, JUNB was necessary and sufficient for EMT in MCF10A cells, but not in A549 lung cancer cells, indicating cell-type-specificity of JUNB EMT-regulatory capacity. Nonetheless, the JUNB-dependence of EMT-associated transcriptional reprogramming in MCF10A cells allowed to define a gene expression signature which was regulated by TGFβ1 in diverse cellular backgrounds, showed positively correlated expression with TGFβ signaling in multiple cancer transcriptomes, and was predictive of patient survival in several cancer types. Altogether, our findings provide novel mechanistic insights into the context-dependent control of TGFβ1-driven EMT and thereby may lead to improved diagnostic and therapeutic options.

A data-driven optimization method for coarse-graining gene regulatory networks

One major challenge in systems biology is to understand how various genes in a gene regulatory network (GRN) collectively perform their functions and control network dynamics. This task becomes extremely hard to tackle in the case of large networks with hundreds of genes and edges, many of which have redundant regulatory roles and functions. The existing methods for model reduction usually require the detailed mathematical description of dynamical systems and their corresponding kinetic parameters, which are often not available. Here, we present a data-driven method for coarse-graining large GRNs, named SacoGraci, using ensemble-based mathematical modeling, dimensionality reduction, and gene circuit optimization by Markov Chain Monte Carlo methods. SacoGraci requires network topology as the only input and is robust against errors in GRNs. We benchmark and demonstrate its usage with synthetic, literature-based, and bioinformatics-derived GRNs. We hope SacoGraci will enhance our ability to model the gene regulation of complex biological systems.

The versatile roles of ADAM8 in cancer cell migration, mechanics, and extracellular matrix remodeling

The posttranslational proteolytic cleavage is a unique and irreversible process that governs the function and half-life of numerous proteins. Thereby the role of the family of A disintegrin and metalloproteases (ADAMs) plays a leading part. A member of this family, ADAM8, has gained attention in regulating disorders, such as neurogenerative diseases, immune function and cancer, by attenuating the function of proteins nearby the extracellular membrane leaflet. This process of "ectodomain shedding" can alter the turnover rate of a number of transmembrane proteins that function in cell adhesion and receptor signal transduction. In the past, the major focus of research about ADAMs have been on neurogenerative diseases, such as Alzheimer, however, there seems to be evidence for a connection between ADAM8 and cancer. The role of ADAMs in the field of cancer research has gained recent attention, but it has been not yet been extensively addressed. Thus, this review article highlights the various roles of ADAM8 with particular emphasis on pathological conditions, such as cancer and malignant cancer progression. Here, the shedding function, direct and indirect matrix degradation, effects on cancer cell mobility and transmigration, and the interplay of ADAM8 with matrix-embedded neighboring cells are presented and discussed. Moreover, the most probable mechanical impact of ADAM8 on cancer cells and their matrix environment is addressed and debated. In summary, this review presents recent advances in substrates/ligands and functions of ADAM8 in its new role in cancer and its potential link to cell mechanical properties and discusses matrix mechanics modifying properties. A deeper comprehension of the regulatory mechanisms governing the expression, subcellular localization, and activity of ADAM8 is expected to reveal appropriate drug targets that will permit a more tailored and fine-tuned modification of its proteolytic activity in cancer development and metastasis.

Dissected subgroups predict the risk of recurrence of stage II colorectal cancer and select rational treatment

Background

Stage II colorectal cancer(CRC) patients after surgery alone have a five-year survival rate of ~60-80%; the incremental benefit of adjuvant chemotherapy is <5%. Predicting risk of recurrence and selecting effective personalized adjuvant drugs for stage II CRC using formalin-fixed, paraffin-embedded(FFPE) samples is a major challenge.

Methods

1319 stage II CRC patients who enrolled in 2011-2019 at Sun Yat-sen University Cancer Center were screened. RNAseq data of FFPE tumor samples of 222 stage II microsatellite stable(MSS) CRC patients(recurrence (n=47), norecurrence (n=175), median follow-up=41 months) were used to develop a method TFunctionalProg for dissecting heterogeneous subgroups of recurrence and predicting risk of recurrence.

Results

TFunctionalProg showed significant predictive values in 222 stage II MSS CRCs. The TFunctionalProg low-risk group had significantly better recurrence free survival (validation set: HR=4.78, p-value=1e-4, low-risk group three-year recurrence free survival=92.6%, high-risk group three-year recurrence free survival=59.7%). TFunctionalProg dissected two subgroups of transition states of stage II MSS CRCs at a high risk of recurrence; each state displays distinct levels of hybrid epithelial-mesenchymal traits, CD8+ T cell suppression mechanisms and FOLFOX resistance. Based on mechanisms in two subgroups, TFunctionalProg proposed personalized rational adjuvant drug combinations of immunotherapy, chemotherapy and repurposed CNS drugs. TFunctionalProg provides different utilities from ctDNA-based prognostic biomarkers.

Conclusion

TFunctionalProg was validated using FFPE samples to predict the risk of recurrence and propose rational adjuvant drug combinations for stage II CRC.

Cancer plasticity: Investigating the causes for this agility

Cancer is not a hard-wired phenomenon but an evolutionary disease. From the onset of carcinogenesis, cancer cells continuously adapt and evolve to satiate their ever-growing proliferation demands. This results in the formation of multiple subtypes of cancer cells with different phenotypes, cellular compositions, and consequently displaying varying degrees of tumorigenic identity and function. This phenomenon is referred to as cancer plasticity, during which the cancer cells exist in a plethora of cellular states having distinct phenotypes. With the advent of modern technologies equipped with enhanced resolution and depth, for example, single-cell RNA-sequencing and advanced computational tools, unbiased cancer profiling at a single-cell resolution are leading the way in understanding cancer cell rewiring both spatially and temporally. In this review, the processes and mechanisms that give rise to cancer plasticity include both intrinsic genetic factors such as epigenetic changes, differential expression due to changes in DNA, RNA, or protein content within the cancer cell, as well as extrinsic environmental factors such as tissue perfusion, extracellular milieu are detailed and their influence on key cancer plasticity hallmarks such as epithelial-mesenchymal transition (EMT) and cancer cell stemness (CSCs) are discussed. Due to therapy evasion and drug resistance, tumor heterogeneity caused by cancer plasticity has major therapeutic ramifications. Hence, it is crucial to comprehend all the cellular and molecular mechanisms that control cellular plasticity. How this process evades therapy, and the therapeutic avenue of targeting cancer plasticity must be diligently investigated.

Regulation of epithelial-mesenchymal transition by tumor microenvironmental signals and its implication in cancer therapeutics

Epithelial-mesenchymal transition (EMT) has been implicated in various aspects of tumor development, including tumor invasion and metastasis, cancer stemness, and therapy resistance. Diverse stroma cell types along with biochemical and biophysical factors in the tumor microenvironment impinge on the EMT program to impact tumor progression. Here we provide an in-depth review of various tumor microenvironmental signals that regulate EMT in cancer. We discuss the molecular mechanisms underlying the role of EMT in therapy resistance and highlight new therapeutic approaches targeting the tumor microenvironment to impact EMT and tumor progression.

"In medio stat virtus": Insights into hybrid E/M phenotype attitudes

Epithelial-mesenchymal plasticity (EMP) refers to the ability of cells to dynamically interconvert between epithelial (E) and mesenchymal (M) phenotypes, thus generating an array of hybrid E/M intermediates with mixed E and M features. Recent findings have demonstrated how these hybrid E/M rather than fully M cells play key roles in most of physiological and pathological processes involving EMT. To this regard, the onset of hybrid E/M state coincides with the highest stemness gene expression and is involved in differentiation of either normal and cancer stem cells. Moreover, hybrid E/M cells are responsible for wound healing and create a favorable immunosuppressive environment for tissue regeneration. Nevertheless, hybrid state is responsible of metastatic process and of the increasing of survival, apoptosis and therapy resistance in cancer cells. The present review aims to describe the main features and the emerging concepts regulating EMP and the formation of E/M hybrid intermediates by describing differences and similarities between cancer and normal hybrid stem cells. In particular, the comprehension of hybrid E/M cells biology will surely advance our understanding of their features and how they could be exploited to improve tissue regeneration and repair.

Bidirectional crosstalk between epithelial-mesenchymal plasticity and IFNγ-induced PD-L1 expression promotes tumour progression

Epithelial-mesenchymal transition (EMT) and immunoevasion through upregulation of programmed death-ligand 1 (PD-L1) are important drivers of cancer progression. While EMT has been proposed to facilitate PD-L1-mediated immunosuppression, molecular mechanisms of their interaction remain obscure. Here, we provide insight into these mechanisms by proposing a mathematical model that describes the crosstalk between EMT and interferon gamma (IFNγ)-induced PD-L1 expression. Our model shows that via interaction with microRNA-200 (miR-200), the multi-stability of the EMT regulatory circuit is mirrored in PD-L1 levels, which are further amplified by IFNγ stimulation. This IFNγ-mediated effect is most prominent for cells in a fully mesenchymal state and less strong for those in an epithelial or partially mesenchymal state. In addition, bidirectional crosstalk between miR-200 and PD-L1 implies that IFNγ stimulation allows cells to undergo EMT for lower amounts of inducing signal, and the presence of IFNγ accelerates EMT and decelerates mesenchymal-epithelial transition (MET). Overall, our model agrees with published findings and provides insight into possible mechanisms behind EMT-mediated immune evasion, and primary, adaptive, or acquired resistance to immunotherapy. Our model can be used as a starting point to explore additional crosstalk mechanisms, as an improved understanding of these mechanisms is indispensable for developing better diagnostic and therapeutic options for cancer patients.

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.

Interleukin-34 promotes the proliferation and epithelial-mesenchymal transition of gastric cancer cells

BACKGROUND

Interleukin (IL)-34 is a pro-inflammatory cytokine involved in tumor development. The role of IL-34 in the proliferation and epithelial-mesenchymal transition (EMT) of gastric cancer (GC) remains to be investigated.

AIM

To investigate whether and how IL-34 affects the proliferation of GC cells and EMT.

METHODS

Using immunohistochemical staining, the expression of IL-34 protein was detected in 60 paired GC and normal paracancerous tissues and the relationship between IL-34 and clinicopathological factors was analyzed. The expression of IL-34 mRNA and protein in normal gastric epithelial cell lines and GC was detected using quantitative real-time polymerase chain reaction (qRT-PCR) and western blotting, respectively. Stable IL-34 knockdown and overexpression in AGS cell lines were established by lentiviral infection and validated by qRT-PCR and western blotting. The cholecystokinin-8 assay, clone formation assay, cell scratch assay, and transwell system were used to detect GC cell proliferation, clone formation, migration, and invasion capacity, respectively. The effects of IL-34 on the growth of GC transplant tumors were assessed using a subcutaneous transplant tumor assay in nude mice. The effects of IL-34 on the expression level of EMT-associated proteins in AGS cells were examined by western blotting.

RESULTS

Expression of IL-34 protein and mRNA was higher in GC cell lines than in GES-1 cells. Compared to matched normal paraneoplastic tissues, the expression of IL-34 protein was higher in 60 GC tissues, which was correlated with tumor size, T-stage, N-stage, tumor, node and metastasis stage, and degree of differentiation. Knockdown of IL-34 expression inhibited the proliferation, clone formation, migration, and invasion of AGS cells, while overexpression of IL-34 promoted cell proliferation, clone formation, migration, and invasion. Furthermore, the reduction of IL-34 promoted the expression of E-cadherin in AGS cells but inhibited the expression of vimentin and N-cadherin. Overexpression of IL-34 inhibited E-cadherin expression but promoted expression of vimentin and N-cadherin in AGS cells. Overexpression of IL-34 promoted the growth of subcutaneous transplanted tumors in nude mice.

CONCLUSION

IL-34 expression is increased in GC tissues and cell lines compared to normal gastric tissues or cell lines. In GC cells, IL-34 promoted proliferation, clone formation, migration, and invasion by regulating EMT-related protein expression cells. Interference with IL-34 may represent a novel strategy for diagnosis and targeted therapy of GC.