Shared enhancer gene regulatory networks between wound and oncogenic programs
Abstract
Wound response programs are often activated during neoplastic growth in tumors. In both wound repair and tumor growth, cells respond to acute stress and balance the activation of multiple programs including apoptosis, proliferation, and cell migration. Central to those responses are the activation of the JNK/MAPK and JAK/STAT signaling pathways. Yet, to what extent these signaling cascades interact at the cis-regulatory level, and how they orchestrate different regulatory and phenotypic responses is still unclear. Here, we aim to characterize the regulatory states that emerge and cooperate in the wound response, using the Drosophila melanogaster wing disc as a model system, and compare these with cancer cell states induced by rasV12scrib-/- in the eye disc. We used single-cell multiome profiling to derive enhancer Gene Regulatory Networks (eGRNs) by integrating chromatin accessibility and gene expression signals. We identify a 'proliferative' eGRN, active in the majority of wounded cells and controlled by AP-1 and STAT. In a smaller, but distinct population of wound cells, a 'senescent' eGRN is activated and driven by C/EBP-like transcription factors (Irbp18, Xrp1, Slow border, and Vrille) and Scalloped. These two eGRN signatures are found to be active in tumor cells, at both gene expression and chromatin accessibility levels. Our single-cell multiome and eGRNs resource offers an in-depth characterisation of the senescence markers, together with a new perspective on the shared gene regulatory programs acting during wound response and oncogenesis.
Data availability
Single-cell sequencing data and aligned matrices have been deposited in GEO (accession code GSE205401)
-
Shared enhancer gene regulatory networks between wound and oncogenic programsNCBI Gene Expression Omnibus, GSE205401.
-
Single-cell transcriptomic analysis of 96hour AEL wild type wing imaginal discsNCBI Gene Expression Omnibus,GSE133204.
-
Gene expression atlas of a developing tissue by single cell expression correlation analysisNCBI Gene Expression Omnibus,GSE127832.
-
Single-cell transcriptomic analysis of the scrib mutant wing imaginal discsNCBI Gene Expression Omnibus,GSE130566.
Article and author information
Author details
Funding
European Research Council (724226_cisCONTROL)
- Valerie M Christiaens
- Gert J Hulselmans
- Stein Aerts
Fonds Wetenschappelijk Onderzoek (G0C0417N)
- Xiaojiang Quan
- Duygu Koldere
Fonds Wetenschappelijk Onderzoek (G094121N)
- Swann Floc'hlay
Deutsche Forschungsgemeinschaft (EXC-2189)
- Anne Classen
Deutsche Forschungsgemeinschaft (CL490/3-1)
- Anne Classen
Deutsche Forschungsgemeinschaft (EXC 2030)
- Mirka Uhlirova
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Reviewing Editor
- Utpal Banerjee, University of California, Los Angeles, United States
Version history
- Received: June 17, 2022
- Preprint posted: June 19, 2022 (view preprint)
- Accepted: May 2, 2023
- Accepted Manuscript published: May 3, 2023 (version 1)
- Version of Record published: June 2, 2023 (version 2)
Copyright
© 2023, Floc'hlay et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,680
- views
-
- 274
- downloads
-
- 4
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Evolutionary Biology
As the genome encodes the information crucial for cell growth, a sizeable genomic deficiency often causes a significant decrease in growth fitness. Whether and how the decreased growth fitness caused by genome reduction could be compensated by evolution was investigated here. Experimental evolution with an Escherichia coli strain carrying a reduced genome was conducted in multiple lineages for approximately 1000 generations. The growth rate, which largely declined due to genome reduction, was considerably recovered, associated with the improved carrying capacity. Genome mutations accumulated during evolution were significantly varied across the evolutionary lineages and were randomly localized on the reduced genome. Transcriptome reorganization showed a common evolutionary direction and conserved the chromosomal periodicity, regardless of highly diversified gene categories, regulons, and pathways enriched in the differentially expressed genes. Genome mutations and transcriptome reorganization caused by evolution, which were found to be dissimilar to those caused by genome reduction, must have followed divergent mechanisms in individual evolutionary lineages. Gene network reconstruction successfully identified three gene modules functionally differentiated, which were responsible for the evolutionary changes of the reduced genome in growth fitness, genome mutation, and gene expression, respectively. The diversity in evolutionary approaches improved the growth fitness associated with the homeostatic transcriptome architecture as if the evolutionary compensation for genome reduction was like all roads leading to Rome.
-
- Computational and Systems Biology
- Physics of Living Systems
Locomotion is a fundamental behavior of Caenorhabditis elegans (C. elegans). Previous works on kinetic simulations of animals helped researchers understand the physical mechanisms of locomotion and the muscle-controlling principles of neuronal circuits as an actuator part. It has yet to be understood how C. elegans utilizes the frictional forces caused by the tension of its muscles to perform sequenced locomotive behaviors. Here, we present a two-dimensional rigid body chain model for the locomotion of C. elegans by developing Newtonian equations of motion for each body segment of C. elegans. Having accounted for friction-coefficients of the surrounding environment, elastic constants of C. elegans, and its kymogram from experiments, our kinetic model (ElegansBot) reproduced various locomotion of C. elegans such as, but not limited to, forward-backward-(omega turn)-forward locomotion constituting escaping behavior and delta-turn navigation. Additionally, ElegansBot precisely quantified the forces acting on each body segment of C. elegans to allow investigation of the force distribution. This model will facilitate our understanding of the detailed mechanism of various locomotive behaviors at any given friction-coefficients of the surrounding environment. Furthermore, as the model ensures the performance of realistic behavior, it can be used to research actuator-controller interaction between muscles and neuronal circuits.