Termination codon readthrough of NNAT mRNA regulates calcium-mediated neuronal differentiation.

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.

Oncogenic gain of function due to p53 amyloids occurs through aberrant alteration of cell cycle and proliferation.

Transcription factor p53 (also known as TP53) has been shown to aggregate into cytoplasmic and nuclear inclusions, compromising its native tumor suppressive functions. Recently, p53 has been shown to form amyloids, which play a role in conferring cancerous properties to cells, leading to tumorigenesis. However, the exact pathways involved in p53 amyloid-mediated cellular transformations are unknown. Here, using an in cellulo model of full-length p53 amyloid formation, we demonstrate the mechanism of loss of p53 tumor-suppressive function with concomitant oncogenic gain of functions. Global gene expression profiling of cells suggests that p53 amyloid formation dysregulates genes associated with the cell cycle, proliferation, apoptosis and senescence along with major signaling pathways. This is further supported by a proteome analysis, showing a significant alteration in levels of p53 target proteins and enhanced metabolism, which enables the survival of cells. Our data indicate that specifically targeting the key molecules in pathways affected by p53 amyloid formation, such as cyclin-dependent kinase-1, leads to loss of the oncogenic phenotype and induces apoptosis of cells. Overall, our work establishes the mechanism of the transformation of cells due to p53 amyloids leading to cancer pathogenesis. This article has an associated First Person interview with the first author of the paper.

Identification and functional characterization of mRNAs that exhibit stop codon readthrough in Arabidopsis thaliana.

Stop codon readthrough (SCR) is the process of continuation of translation beyond the stop codon, generating protein isoforms with C-terminal extensions. SCR has been observed in viruses, fungi, and multicellular organisms, including mammals. However, SCR is largely unexplored in plants. In this study, we have analyzed ribosome profiling datasets to identify mRNAs that exhibit SCR in Arabidopsis thaliana. Analyses of the ribosome density, ribosome coverage, and three-nucleotide periodicity of the ribosome profiling reads in the mRNA region downstream of the stop codon provided strong evidence for SCR in mRNAs of 144 genes. We show that SCR generated putative evolutionarily conserved nuclear localization signals, transmembrane helices, and intrinsically disordered regions in the C-terminal extensions of several of these proteins. Furthermore, gene ontology functional enrichment analysis revealed that these 144 genes belong to three major functional groups-translation, photosynthesis, and abiotic stress tolerance. Using a luminescence-based readthrough assay, we experimentally demonstrated SCR in representative mRNAs belonging to each of these functional classes. Finally, using microscopy, we show that the SCR product of one gene that contains a nuclear localization signal at the C-terminal extension, CURT1B, localizes to the nucleus as predicted. Based on these observations, we propose that SCR plays an important role in plant physiology by regulating protein localization and function.

Let-7a-regulated translational readthrough of mammalian AGO1 generates a microRNA pathway inhibitor.

Translational readthrough generates proteins with extended C-termini, which often possess distinct properties. Here, we have used various reporter assays to demonstrate translational readthrough of AGO1 mRNA. Analysis of ribosome profiling data and mass spectrometry data provided additional evidence for translational readthrough of AGO1. The endogenous readthrough product, Ago1x, could be detected by a specific antibody both in vitro and in vivo. This readthrough process is directed by a cis sequence downstream of the canonical AGO1 stop codon, which is sufficient to drive readthrough even in a heterologous context. This cis sequence has a let-7a miRNA-binding site, and readthrough is promoted by let-7a miRNA. Interestingly, Ago1x can load miRNAs on target mRNAs without causing post-transcriptional gene silencing, due to its inability to interact with GW182. Because of these properties, Ago1x can serve as a competitive inhibitor of miRNA pathway. In support of this, we observed increased global translation in cells overexpressing Ago1x. Overall, our results reveal a negative feedback loop in the miRNA pathway mediated by the translational readthrough product of AGO1.

Lhx2 in germ cells suppresses endothelial cell migration in the developing ovary.

LIM-homeobox genes play multiple roles in developmental processes, but their roles in gonad development are not completely understood. Herein, we report that Lhx2, Ils2, Lmx1a, and Lmx1b are expressed in a sexually dimorphic manner in mouse, rat, and human gonads during sex determination. Amongst these, Lhx2 has female biased expression in the developing gonads of species with environmental and genetic modes of sex determination. Single-cell RNAseq analysis revealed that Lhx2 is exclusively expressed in the germ cells of the developing mouse ovaries. To elucidate the roles of Lhx2 in the germ cells, we analyzed the phenotypes of Lhx2 knockout XX gonads. While the gonads developed appropriately in Lhx2 knockout mice and the somatic cells were correctly specified in the developing ovaries, transcriptome analysis revealed enrichment of genes in the angiogenesis pathway. There was an elevated expression of several pro-angiogenic factors in the Lhx2 knockout ovaries. The elevated expression of pro-angiogenic factors was associated with an increase in numbers of endothelial cells in the Lhx2-/- ovaries at E13.5. Gonad recombination assays revealed that the increased numbers of endothelial cells in the XX gonads in absence of Lhx2 was due to ectopic migration of endothelial cells in a cell non-autonomous manner. We also found that, there was increased expression of several endothelial cell-enriched male-biased genes in Lhx2 knockout ovaries. Also, in absence of Lhx2, the migrated endothelial cells formed an angiogenic network similar to that of the wild type testis, although the coelomic blood vessel did not form. Together, our results suggest that Lhx2 in the germ cells is required to suppress vascularization in the developing ovary. These results suggest a need to explore the roles of germ cells in the control of vascularization in developing gonads. Preprint version of the article is available on BioRxiv at https://doi.org/10.1101/2022.03.07.483280.

A Computational Systems Biology Approach Identifies SLUG as a Mediator of Partial Epithelial-Mesenchymal Transition (EMT).

Epithelial-mesenchymal plasticity comprises reversible transitions among epithelial, hybrid epithelial/mesenchymal (E/M) and mesenchymal phenotypes, and underlies various aspects of aggressive tumor progression such as metastasis, therapy resistance, and immune evasion. The process of cells attaining one or more hybrid E/M phenotypes is termed as partial epithelial mesenchymal transition (EMT). Cells in hybrid E/M phenotype(s) can be more aggressive than those in either fully epithelial or mesenchymal state. Thus, identifying regulators of hybrid E/M phenotypes is essential to decipher the rheostats of phenotypic plasticity and consequent accelerators of metastasis. Here, using a computational systems biology approach, we demonstrate that SLUG (SNAIL2) - an EMT-inducing transcription factor - can inhibit cells from undergoing a complete EMT and thus stabilize them in hybrid E/M phenotype(s). It expands the parametric range enabling the existence of a hybrid E/M phenotype, thereby behaving as a phenotypic stability factor. Our simulations suggest that this specific property of SLUG emerges from the topology of the regulatory network it forms with other key regulators of epithelial-mesenchymal plasticity. Clinical data suggest that SLUG associates with worse patient prognosis across multiple carcinomas. Together, our results indicate that SLUG can stabilize hybrid E/M phenotype(s).

Prioritizing and characterizing functionally relevant genes across human tissues.

Knowledge of genes that are critical to a tissue's function remains difficult to ascertain and presents a major bottleneck toward a mechanistic understanding of genotype-phenotype links. Here, we present the first machine learning model-FUGUE-combining transcriptional and network features, to predict tissue-relevant genes across 30 human tissues. FUGUE achieves an average cross-validation auROC of 0.86 and auPRC of 0.50 (expected 0.09). In independent datasets, FUGUE accurately distinguishes tissue or cell type-specific genes, significantly outperforming the conventional metric based on tissue-specific expression alone. Comparison of tissue-relevant transcription factors across tissue recapitulate their developmental relationships. Interestingly, the tissue-relevant genes cluster on the genome within topologically associated domains and furthermore, are highly enriched for differentially expressed genes in the corresponding cancer type. We provide the prioritized gene lists in 30 human tissues and an open-source software to prioritize genes in a novel context given multi-sample transcriptomic data.

Stop codon read-through of mammalian MTCH2 leading to an unstable isoform regulates mitochondrial membrane potential.

Stop codon read-through (SCR) is a process of continuation of translation beyond a stop codon. This phenomenon, which occurs only in certain mRNAs under specific conditions, leads to a longer isoform with properties different from that of the canonical isoform. MTCH2, which encodes a mitochondrial protein that regulates mitochondrial metabolism, was selected as a potential read-through candidate based on evolutionary conservation observed in the proximal region of its 3' UTR. Here, we demonstrate translational read-through across two evolutionarily conserved, in-frame stop codons of MTCH2 using luminescence- and fluorescence-based assays, and by analyzing ribosome-profiling and mass spectrometry (MS) data. This phenomenon generates two isoforms, MTCH2x and MTCH2xx (single- and double-SCR products, respectively), in addition to the canonical isoform MTCH2, from the same mRNA. Our experiments revealed that a cis-acting 12-nucleotide sequence in the proximal 3' UTR of MTCH2 is the necessary signal for SCR. Functional characterization showed that MTCH2 and MTCH2x were localized to mitochondria with a long t1/2 (>36 h). However, MTCH2xx was found predominantly in the cytoplasm. This mislocalization and its unique C terminus led to increased degradation, as shown by greatly reduced t1/2 (<1 h). MTCH2 read-through-deficient cells, generated using CRISPR-Cas9, showed increased MTCH2 expression and, consistent with this, decreased mitochondrial membrane potential. Thus, double-SCR of MTCH2 regulates its own expression levels contributing toward the maintenance of normal mitochondrial membrane potential.

Metabolic origin of the fused aminoacyl-tRNA synthetase, glutamyl-prolyl-tRNA synthetase.

About 1 billion years ago, in a single-celled holozoan ancestor of all animals, a gene fusion of two tRNA synthetases formed the bifunctional enzyme, glutamyl-prolyl-tRNA synthetase (EPRS). We propose here that a confluence of metabolic, biochemical, and environmental factors contributed to the specific fusion of glutamyl- (ERS) and prolyl- (PRS) tRNA synthetases. To test this idea, we developed a mathematical model that centers on the precursor-product relationship of glutamic acid and proline, as well as metabolic constraints on free glutamic acid availability near the time of the fusion event. Our findings indicate that proline content increased in the proteome during the emergence of animals, thereby increasing demand for free proline. Together, these constraints contributed to a marked cellular depletion of glutamic acid and its products, with potentially catastrophic consequences. In response, an ancient organism invented an elegant solution in which genes encoding ERS and PRS fused to form EPRS, forcing coexpression of the two enzymes and preventing lethal dysregulation. The substantial evolutionary advantage of this coregulatory mechanism is evidenced by the persistence of EPRS in nearly all extant animals.

Pierre Robin sequence coexisting with ectopic kidney, bilateral developmental hip dysplasia in a child - A rare case report.

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.

Mutually exclusive teams-like patterns of gene regulation characterize phenotypic heterogeneity along the noradrenergic-mesenchymal axis in neuroblastoma.

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.

Elucidating the Role of MicroRNA-18a in Propelling a Hybrid Epithelial-Mesenchymal Phenotype and Driving Malignant Progression in ER-Negative Breast Cancer.

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.

Spatiotemporal modulation of growth factors directs the generation of multilineage mouse embryonic stem cell-derived mammary organoids.

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.

Increased prevalence of hybrid epithelial/mesenchymal state and enhanced phenotypic heterogeneity in basal breast cancer.

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.

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.

Dynamical modelling of proliferative-invasive plasticity and IFNγ signaling in melanoma reveals mechanisms of PD-L1 expression heterogeneity.

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.

Dynamical modeling of proliferative-invasive plasticity and IFNγ signaling in melanoma reveals mechanisms of PD-L1 expression heterogeneity.

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.

Prognostic association of immunoproteasome expression in solid tumours is governed by the immediate immune environment.

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.

Transcriptional state dynamics lead to heterogeneity and adaptive tumor evolution in urothelial bladder carcinoma.

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.

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.