ABSTRACT
The ultraviolet (UV) radiation triggers a pigmentation response in human skin, wherein, melanocytes rapidly activate divergent maturation and proliferation programs. Using single-cell sequencing, we demonstrate that these 2 programs are segregated in distinct subpopulations in melanocytes of human and zebrafish skin. The coexistence of these 2 cell states in cultured melanocytes suggests possible cell autonomy. Luria-Delbrück fluctuation test reveals that the initial establishment of these states is stochastic. Tracking of pigmenting cells ascertains that the stochastically acquired state is faithfully propagated in the progeny. A systemic approach combining single-cell multi-omics (RNA+ATAC) coupled to enhancer mapping with H3K27 acetylation successfully identified state-specific transcriptional networks. This comprehensive analysis led to the construction of a gene regulatory network (GRN) that under the influence of noise, establishes a bistable system of pigmentation and proliferation at the population level. This GRN recapitulates melanocyte behaviour in response to external cues that reinforce either of the states. Our work highlights that inherent stochasticity within melanocytes establishes dedicated states, and the mature state is sustained by selective enhancers mark through histone acetylation. While the initial cue triggers a proliferation response, the continued signal activates and maintains the pigmenting subpopulation via epigenetic imprinting. Thereby our study provides the basis of coexistence of distinct populations which ensures effective pigmentation response while preserving the self-renewal capacity.
Subject(s)
Cell Proliferation , Gene Regulatory Networks , Melanocytes , Skin Pigmentation , Zebrafish , Melanocytes/metabolism , Zebrafish/genetics , Animals , Humans , Skin Pigmentation/genetics , Skin Pigmentation/physiology , Stochastic Processes , Cell Differentiation/genetics , Histones/metabolism , Acetylation , Ultraviolet Rays , Single-Cell Analysis , Pigmentation/genetics , Enhancer Elements, Genetic/genetics , Epigenesis, Genetic , Skin/metabolism , Skin/cytologyABSTRACT
Intra-tumoral phenotypic heterogeneity promotes tumor relapse and therapeutic resistance and remains an unsolved clinical challenge. Decoding the interconnections among different biological axes of plasticity is crucial to understand the molecular origins of phenotypic heterogeneity. Here, we use multi-modal transcriptomic data-bulk, single-cell, and spatial transcriptomics-from breast cancer cell lines and primary tumor samples, to identify associations between epithelial-mesenchymal transition (EMT) and luminal-basal plasticity-two key processes that enable heterogeneity. We show that luminal breast cancer strongly associates with an epithelial cell state, but basal breast cancer is associated with hybrid epithelial/mesenchymal phenotype(s) and higher phenotypic heterogeneity. Mathematical modeling of core underlying gene regulatory networks representative of the crosstalk between the luminal-basal and epithelial-mesenchymal axes elucidate mechanistic underpinnings of the observed associations from transcriptomic data. Our systems-based approach integrating multi-modal data analysis with mechanism-based modeling offers a predictive framework to characterize intra-tumor heterogeneity and identify interventions to restrict it.
ABSTRACT
Epigenetic alterations that lead to differential expression of microRNAs (miRNAs/miR) are known to regulate tumour cell states, epithelial-mesenchymal transition (EMT) and the progression to metastasis in breast cancer. This study explores the key contribution of miRNA-18a in mediating a hybrid E/M cell state that is pivotal to the malignant transformation and tumour progression in the aggressive ER-negative subtype of breast cancer. The expression status and associated effects of miR-18a were evaluated in patient-derived breast tumour samples in combination with gene expression data from public datasets, and further validated in in vitro and in vivo breast cancer model systems. The clinical relevance of the study findings was corroborated against human breast tumour specimens (n = 446 patients). The down-regulated expression of miR-18a observed in ER-negative tumours was found to drive the enrichment of hybrid epithelial/mesenchymal (E/M) cells with luminal attributes, enhanced traits of migration, stemness, drug-resistance and immunosuppression. Further analysis of the miR-18a targets highlighted possible hypoxia-inducible factor 1-alpha (HIF-1α)-mediated signalling in these tumours. This is a foremost report that validates the dual role of miR-18a in breast cancer that is subtype-specific based on hormone receptor expression. The study also features a novel association of low miR-18a levels and subsequent enrichment of hybrid E/M cells, increased migration and stemness in a subgroup of ER-negative tumours that may be attributed to HIF-1α mediated signalling. The results highlight the possibility of stratifying the ER-negative disease into clinically relevant groups by analysing miRNA signatures.
Subject(s)
Breast Neoplasms , Epithelial-Mesenchymal Transition , Gene Expression Regulation, Neoplastic , MicroRNAs , MicroRNAs/genetics , MicroRNAs/metabolism , Humans , Epithelial-Mesenchymal Transition/genetics , Breast Neoplasms/genetics , Breast Neoplasms/pathology , Breast Neoplasms/metabolism , Female , Disease Progression , Receptors, Estrogen/metabolism , Receptors, Estrogen/genetics , Cell Line, Tumor , Hypoxia-Inducible Factor 1, alpha Subunit/metabolism , Hypoxia-Inducible Factor 1, alpha Subunit/genetics , Phenotype , Animals , Mice , Cell Movement/geneticsABSTRACT
Pierre Robin sequence is a rare congenital disorder with variable associations. A two-year-old female child was brought in with complaints of an abnormal gait and right lower limb shortening. A comprehensive clinical evaluation unveiled facial characteristics resembling those of Pierre Robin Sequence, including micrognathia, a cleft palate, and glossoptosis. The child also exhibited retrognathia, a low posterior hairline, bilateral developmental dysplasia of the hips, right foot syndactyly, and a left ectopic kidney as confirmed by ultrasonography. Radiological examination showed bilateral developmental dysplastic hips. The connection between Pierre Robin sequence and the presence of bilateral dysplastic hips and an ectopic kidney is an unusual and rare combination.
ABSTRACT
Neuroblastoma is the most frequent extracranial pediatric tumor and leads to 15% of all cancer-related deaths in children. Tumor relapse and therapy resistance in neuroblastoma are driven by phenotypic plasticity and heterogeneity between noradrenergic (NOR) and mesenchymal (MES) cell states. Despite the importance of this phenotypic plasticity, the dynamics and molecular patterns associated with these bidirectional cell-state transitions remain relatively poorly understood. Here, we analyze multiple RNA-seq datasets at both bulk and single-cell resolution, to understand the association between NOR- and MES-specific factors. We observed that NOR-specific and MES-specific expression patterns are largely mutually exclusive, exhibiting a "teams-like" behavior among the genes involved, reminiscent of our earlier observations in lung cancer and melanoma. This antagonism between NOR and MES phenotypes was also associated with metabolic reprogramming and with immunotherapy targets PD-L1 and GD2 as well as with experimental perturbations driving the NOR-MES and/or MES-NOR transition. Further, these "teams-like" patterns were seen only among the NOR- and MES-specific genes, but not in housekeeping genes, possibly highlighting a hallmark of network topology enabling cancer cell plasticity.
Subject(s)
Neoplasm Recurrence, Local , Neuroblastoma , Child , Humans , Neuroblastoma/genetics , Neuroblastoma/pathology , Gene Expression Regulation, Neoplastic , PhenotypeABSTRACT
Ectodermal appendages, such as the mammary gland (MG), are thought to have evolved from hair-associated apocrine glands to serve the function of milk secretion. Through the directed differentiation of mouse embryonic stem cells (mESCs), here, we report the generation of multilineage ESC-derived mammary organoids (MEMOs). We adapted the skin organoid model, inducing the dermal mesenchyme to transform into mammary-specific mesenchyme via the sequential activation of Bone Morphogenetic Protein 4 (BMP4) and Parathyroid Hormone-related Protein (PTHrP) and inhibition of hedgehog (HH) signaling. Using single-cell RNA sequencing, we identified gene expression profiles that demonstrate the presence of mammary-specific epithelial cells, fibroblasts, and adipocytes. MEMOs undergo ductal morphogenesis in Matrigel and can reconstitute the MG in vivo. Further, we demonstrate that the loss of function in placode regulators LEF1 and TBX3 in mESCs results in impaired skin and MEMO generation. In summary, our MEMO model is a robust tool for studying the development of ectodermal appendages, and it provides a foundation for regenerative medicine and disease modeling.
Subject(s)
Hedgehog Proteins , Mouse Embryonic Stem Cells , Mice , Animals , Hedgehog Proteins/metabolism , Mammary Glands, Animal , Epithelial Cells , Cell Differentiation , OrganoidsABSTRACT
Anoikis resistance or evasion of cell death triggered by cell detachment into suspension is a hallmark of cancer that is concurrent with cell survival and metastasis. The effects of frequent matrix detachment encounters on the development of anoikis resistance in cancer remains poorly defined. Here we show using a panel of ovarian cancer models, that repeated exposure to suspension stress in vitro followed by attached recovery growth leads to the development of anoikis resistance paralleling in vivo development of anoikis resistance in ovarian cancer ascites. This resistance is concurrent with enhanced invasion, chemoresistance and the ability of anoikis adapted cells to metastasize to distant sites. Adapted anoikis resistant cells show a heightened dependency on oxidative phosphorylation and can also evade immune surveillance. We find that such acquired anoikis resistance is not genetic, as acquired resistance persists for a finite duration in the absence of suspension stress. Transcriptional reprogramming is however essential to this process, as acquisition of adaptive anoikis resistance in vitro and in vivo is exquisitely sensitive to inhibition of CDK8/19 Mediator kinase, a pleiotropic regulator of transcriptional reprogramming. Our data demonstrate that growth after recovery from repeated exposure to suspension stress is a direct contributor to metastasis and that inhibition of CDK8/19 Mediator kinase during such adaptation provides a therapeutic opportunity to prevent both local and distant metastasis in cancer.
ABSTRACT
Intra-tumor heterogeneity contributes to treatment failure and poor survival in urothelial bladder carcinoma (UBC). Analyzing transcriptome from a UBC cohort, we report that intra-tumor transcriptomic heterogeneity indicates co-existence of tumor cells in epithelial and mesenchymal-like transcriptional states and bi-directional transition between them occurs within and between tumor subclones. We model spontaneous and reversible transition between these partially heritable states in cell lines and characterize their population dynamics. SMAD3, KLF4 and PPARG emerge as key regulatory markers of the transcriptional dynamics. Nutrient limitation, as in the core of large tumors, and radiation treatment perturb the dynamics, initially selecting for a transiently resistant phenotype and then reconstituting heterogeneity and growth potential, driving adaptive evolution. Dominance of transcriptional states with low PPARG expression indicates an aggressive phenotype in UBC patients. We propose that phenotypic plasticity and dynamic, non-genetic intra-tumor heterogeneity modulate both the trajectory of disease progression and adaptive treatment response in UBC.
Subject(s)
Carcinoma, Transitional Cell , Urinary Bladder Neoplasms , Humans , Urinary Bladder , PPAR gamma , Urinary Bladder Neoplasms/therapy , Carcinoma, Transitional Cell/pathology , Disease ProgressionABSTRACT
Intra-tumoral phenotypic heterogeneity promotes tumor relapse and therapeutic resistance and remains an unsolved clinical challenge. It manifests along multiple phenotypic axes and decoding the interconnections among these different axes is crucial to understand its molecular origins and to develop novel therapeutic strategies to control it. Here, we use multi-modal transcriptomic data analysis - bulk, single-cell and spatial transcriptomics - from breast cancer cell lines and primary tumor samples, to identify associations between epithelial-mesenchymal transition (EMT) and luminal-basal plasticity - two key processes that enable heterogeneity. We show that luminal breast cancer strongly associates with an epithelial cell state, but basal breast cancer is associated with hybrid epithelial/mesenchymal phenotype(s) and higher phenotypic heterogeneity. These patterns were inherent in methylation profiles, suggesting an epigenetic crosstalk between EMT and lineage plasticity in breast cancer. Mathematical modelling of core underlying gene regulatory networks representative of the crosstalk between the luminal-basal and epithelial-mesenchymal axes recapitulate and thus elucidate mechanistic underpinnings of the observed associations from transcriptomic data. Our systems-based approach integrating multi-modal data analysis with mechanism-based modeling offers a predictive framework to characterize intra-tumor heterogeneity and to identify possible interventions to restrict it.
ABSTRACT
BACKGROUND: Phenotypic heterogeneity of melanoma cells contributes to drug tolerance, increased metastasis, and immune evasion in patients with progressive disease. Diverse mechanisms have been individually reported to shape extensive intra-tumor and inter-tumor phenotypic heterogeneity, such as IFNγ signaling and proliferative to invasive transition, but how their crosstalk impacts tumor progression remains largely elusive. METHODS: Here, we integrate dynamical systems modeling with transcriptomic data analysis at bulk and single-cell levels to investigate underlying mechanisms behind phenotypic heterogeneity in melanoma and its impact on adaptation to targeted therapy and immune checkpoint inhibitors. We construct a minimal core regulatory network involving transcription factors implicated in this process and identify the multiple 'attractors' in the phenotypic landscape enabled by this network. Our model predictions about synergistic control of PD-L1 by IFNγ signaling and proliferative to invasive transition were validated experimentally in three melanoma cell lines-MALME3, SK-MEL-5 and A375. RESULTS: We demonstrate that the emergent dynamics of our regulatory network comprising MITF, SOX10, SOX9, JUN and ZEB1 can recapitulate experimental observations about the co-existence of diverse phenotypes (proliferative, neural crest-like, invasive) and reversible cell-state transitions among them, including in response to targeted therapy and immune checkpoint inhibitors. These phenotypes have varied levels of PD-L1, driving heterogeneity in immunosuppression. This heterogeneity in PD-L1 can be aggravated by combinatorial dynamics of these regulators with IFNγ signaling. Our model predictions about changes in proliferative to invasive transition and PD-L1 levels as melanoma cells evade targeted therapy and immune checkpoint inhibitors were validated in multiple RNA-seq data sets from in vitro and in vivo experiments. CONCLUSION: Our calibrated dynamical model offers a platform to test combinatorial therapies and provide rational avenues for the treatment of metastatic melanoma. This improved understanding of crosstalk among PD-L1 expression, proliferative to invasive transition and IFNγ signaling can be leveraged to improve the clinical management of therapy-resistant and metastatic melanoma.
Subject(s)
Melanoma , Neoplasms, Second Primary , Humans , B7-H1 Antigen/genetics , Immune Checkpoint Inhibitors , Melanoma/drug therapy , Melanoma/genetics , Cell LineABSTRACT
Termination codon readthrough (TCR) is a process in which ribosomes continue to translate an mRNA beyond a stop codon generating a C-terminally extended protein isoform. Here, we demonstrate TCR in mammalian NNAT mRNA, which encodes NNAT, a proteolipid important for neuronal differentiation. This is a programmed event driven by cis-acting RNA sequences present immediately upstream and downstream of the canonical stop codon and is negatively regulated by NONO, an RNA-binding protein known to promote neuronal differentiation. Unlike the canonical isoform NNAT, we determined that the TCR product (NNATx) does not show detectable interaction with the sarco/endoplasmic reticulum Ca2+-ATPase isoform 2 Ca2+ pump, cannot increase cytoplasmic Ca2+ levels, and therefore does not enhance neuronal differentiation in Neuro-2a cells. Additionally, an antisense oligonucleotide that targets a region downstream of the canonical stop codon reduced TCR of NNAT and enhanced the differentiation of Neuro-2a cells to cholinergic neurons. Furthermore, NNATx-deficient Neuro-2a cells, generated using CRISPR-Cas9, showed increased cytoplasmic Ca2+ levels and enhanced neuronal differentiation. Overall, these results demonstrate regulation of neuronal differentiation by TCR of NNAT. Importantly, this process can be modulated using a synthetic antisense oligonucleotide.
Subject(s)
Calcium , Neurons , Protein Biosynthesis , Animals , Calcium/metabolism , Cell Differentiation , Codon, Terminator , Mammals/metabolism , Oligonucleotides, Antisense/metabolism , Receptors, Antigen, T-Cell/metabolism , RNA, Messenger/genetics , RNA, Messenger/metabolism , Neurons/cytologyABSTRACT
Phenotypic heterogeneity of melanoma cells contributes to drug tolerance, increased metastasis, and immune evasion in patients with progressive disease. Diverse mechanisms have been individually reported to shape extensive intra- and inter-tumoral phenotypic heterogeneity, such as IFNγ signaling and proliferative to invasive transition, but how their crosstalk impacts tumor progression remains largely elusive. Here, we integrate dynamical systems modeling with transcriptomic data analysis at bulk and single-cell levels to investigate underlying mechanisms behind phenotypic heterogeneity in melanoma and its impact on adaptation to targeted therapy and immune checkpoint inhibitors. We construct a minimal core regulatory network involving transcription factors implicated in this process and identify the multiple "attractors" in the phenotypic landscape enabled by this network. Our model predictions about synergistic control of PD-L1 by IFNγ signaling and proliferative to invasive transition were validated experimentally in three melanoma cell lines - MALME3, SK-MEL-5 and A375. We demonstrate that the emergent dynamics of our regulatory network comprising MITF, SOX10, SOX9, JUN and ZEB1 can recapitulate experimental observations about the co-existence of diverse phenotypes (proliferative, neural crest-like, invasive) and reversible cell-state transitions among them, including in response to targeted therapy and immune checkpoint inhibitors. These phenotypes have varied levels of PD-L1, driving heterogeneity in immune-suppression. This heterogeneity in PD-L1 can be aggravated by combinatorial dynamics of these regulators with IFNγ signaling. Our model predictions about changes in proliferative to invasive transition and PD-L1 levels as melanoma cells evade targeted therapy and immune checkpoint inhibitors were validated in multiple data sets from in vitro and in vivo experiments. Our calibrated dynamical model offers a platform to test combinatorial therapies and provide rational avenues for the treatment of metastatic melanoma. This improved understanding of crosstalk among PD-L1 expression, proliferative to invasive transition and IFNγ signaling can be leveraged to improve the clinical management of therapy-resistant and metastatic melanoma.
ABSTRACT
(1) Background: Activation of the receptor tyrosine kinase Axl by Gas6 fosters oncogenic effects in hepatocellular carcinoma (HCC), associating with increased mortality of patients. The impact of Gas6/Axl signaling on the induction of individual target genes in HCC and its consequences is an open issue. (2) Methods: RNA-seq analysis of Gas6-stimulated Axl-proficient or Axl-deficient HCC cells was used to identify Gas6/Axl targets. Gain- and loss-of-function studies as well as proteomics were employed to characterize the role of PRAME (preferentially expressed antigen in melanoma). Expression of Axl/PRAME was assessed in publicly available HCC patient datasets and in 133 HCC cases. (3) Results: Exploitation of well-characterized HCC models expressing Axl or devoid of Axl allowed the identification of target genes including PRAME. Intervention with Axl signaling or MAPK/ERK1/2 resulted in reduced PRAME expression. PRAME levels were associated with a mesenchymal-like phenotype augmenting 2D cell migration and 3D cell invasion. Interactions with pro-oncogenic proteins such as CCAR1 suggested further tumor-promoting functions of PRAME in HCC. Moreover, PRAME showed elevated expression in Axl-stratified HCC patients, which correlates with vascular invasion and lowered patient survival. (4) Conclusions: PRAME is a bona fide target of Gas6/Axl/ERK signaling linked to EMT and cancer cell invasion in HCC.
ABSTRACT
Induction of immunoproteasome (IP) expression in tumour cells can enhance antigen presentation and immunogenicity. Recently, the overexpression of IP genes has been associated with better prognosis and response to immune checkpoint blockade (ICB) therapies in melanoma. However, the extent of this association in other solid tumours and how that is influenced by tumour cell-intrinsic and cell-extrinsic factors remain unclear. Here, we address this by exploring the gene expression patterns from available bulk and single-cell transcriptomic data of primary tumours. We find that tumours with high-IP expression exhibit cytotoxic immune cell infiltration and upregulation of IFN-γ and TNF-α pathways in tumour cells. However, the association of IP expression with overall survival (TCGA cohort) and response to ICB therapy (non-TCGA cohorts) is tumour-type specific (better in non-small-cell lung, breast, bladder and thymus; and worse in glioma and renal) and is greatly influenced by pro- or antitumourigenic immune cell infiltration patterns. This emphasises the need for considering immune cell infiltration patterns, along with IP expression, as a prognostic biomarker to predict overall survival or response to ICB therapies in solid tumours, besides melanoma.
Subject(s)
Carcinoma, Non-Small-Cell Lung , Lung Neoplasms , Melanoma , Humans , Prognosis , Melanoma/pathology , Gene Expression ProfilingABSTRACT
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.
ABSTRACT
Background: The initial insights from the studies on COVID-19 had been disappointing, indicating the necessity of an aggravated search for alternative strategies. In this regard, the adjunct potential of yoga has been proposed for enhancing the effectiveness of the standard of care with respect to COVID-19 management. We tested whether a telemodel of yoga intervention could aid in better clinical management for hospitalized patients with mild-to-moderate COVID-19 when complemented with the standard of care. Methods: This was a randomized controlled trial conducted at the Narayana Hrudyalaya, Bengaluru, India, on hospitalized patients with mild-to-moderate COVID-19 infection enrolled between 31 May and 22 July 2021. The patients (n = 225) were randomized in a 1:1 ratio [adjunct tele-yoga (n = 113) or standard of care]. The adjunct yoga group received intervention in tele-mode within 4-h post-randomization until 14 days along with the standard of care. The primary outcome was the clinical status on day 14 post-randomization, assessed with a seven-category ordinal scale. The secondary outcome set included scores on the COVID Outcomes Scale on day 7, follow-up for clinical status and all-cause mortality on day 28, post-randomization, duration of days at the hospital, 5th-day changes post-randomization for viral load expressed as cyclic threshold (Ct), and inflammatory markers and perceived stress scores on day 14. Results: As compared with the standard of care alone, the proportional odds of having a higher score on the 7-point ordinal scale on day 14 were ~1.8 for the adjunct tele-yoga group (OR = 1.83, 95% CI, 1.11-3.03). On day 5, there were significant reductions in CRP (P = 0.001) and LDH levels (P = 0.029) in the adjunct yoga group compared to the standard of care alone. CRP reduction was also observed as a potential mediator for the yoga-induced improvement of clinical outcomes. The Kaplan-Meier estimate of all-cause mortality on day 28 was the adjusted hazard ratio (HR) of 0.26 (95% CI, 0.05-1.30). Conclusion: The observed 1.8-fold improvement in the clinical status on day 14 of patients of COVID-19 with adjunct use of tele-yoga contests its use as a complementary treatment in hospital settings.
Subject(s)
COVID-19 , Yoga , Humans , COVID-19/therapy , SARS-CoV-2 , IndiaABSTRACT
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.
ABSTRACT
Epithelial-mesenchymal transition (EMT) and its reverse mesenchymal-epithelial transition (MET) are critical during embryonic development, wound healing and cancer metastasis. While phenotypic changes during short-term EMT induction are reversible, long-term EMT induction has been often associated with irreversibility. Here, we show that phenotypic changes seen in MCF10A cells upon long-term EMT induction by TGFß need not be irreversible, but have relatively longer time scales of reversibility than those seen in short-term induction. Next, using a phenomenological mathematical model to account for the chromatin-mediated epigenetic silencing of the miR-200 family by ZEB family, we highlight how the epigenetic memory gained during long-term EMT induction can slow the recovery to the epithelial state post-TGFß withdrawal. Our results suggest that epigenetic modifiers can govern the extent and time scale of EMT reversibility and advise caution against labelling phenotypic changes seen in long-term EMT induction as 'irreversible'.
Subject(s)
Epigenetic Memory , Epithelial-Mesenchymal Transition , Epigenesis, Genetic , Transforming Growth Factor betaABSTRACT
Epithelial-to-mesenchymal transition (EMT), a well-established phenomenon studied across pan-cancer types, has long been known to be a major player in driving tumor invasion and metastasis. Recent studies have highlighted the importance of partial EMT phenotypes in metastasis. Initially thought as a transitional state between epithelial and mesenchymal phenotypic states, partial EMT state is now widely recognized as a key driver of intra-tumoral heterogeneity and phenotypic plasticity, further accelerating tumor metastasis and therapeutic resistance. However, how tumor microenvironment regulates partial EMT phenotypes remains unclear. We have developed unique size-controlled three-dimensional microtumor models that recapitulate tumor-intrinsic hypoxia and the emergence of collectively migrating cells. In this study, we further interrogate these microtumor models to understand how tumor-intrinsic hypoxia regulates partial EMT and collective migration in hypoxic large microtumors fabricated from T47D breast cancer cells. We compared global gene expression profiles of hypoxic, migratory microtumors to that of non-hypoxic, non-migratory microtumors at early and late time-points. Using our microtumor models, we identified unique gene signatures for tumor-intrinsic hypoxia (early versus late), partial EMT and migration (pre-migratory versus migratory phenotype). Through differential gene expression analysis between the microtumor models with an overlap of hypoxia, partial EMT and migration signatures, we identified prolyl 4-hydroxylase subunit 2 (P4HA2), a hypoxia responsive gene, as a central regulator common to hypoxia, partial EMT and collective migration. Further, the inhibition of P4HA2 significantly blocked collective migration in hypoxic microtumors. Thus, using the integrated computational-experimental analysis, we identify the key role of P4HA2 in tumor-intrinsic hypoxia-driven partial EMT and collective migration.
ABSTRACT
The immense regenerative potential of the liver is attributed to the ability of its two key cell types - hepatocytes and cholangiocytes - to trans-differentiate to one another either directly or through intermediate progenitor states. However, the dynamic features of decision-making between these cell-fates during liver development and regeneration remains elusive. Here, we identify a core gene regulatory network comprising c/EBPα, TGFBR2, and SOX9 which is multistable in nature, enabling three distinct cell states - hepatocytes, cholangiocytes, and liver progenitor cells (hepatoblasts/oval cells) - and stochastic switching among them. Predicted expression signature for these three states are validated through multiple bulk and single-cell transcriptomic datasets collected across developmental stages and injury-induced liver repair. This network can also explain the experimentally observed spatial organization of phenotypes in liver parenchyma and predict strategies for efficient cellular reprogramming. Our analysis elucidates how the emergent dynamics of underlying regulatory networks drive diverse cell-fate decisions in liver development and regeneration.