Shared enhancer gene regulatory networks between wound and oncogenic programs

  1. Swann Floc'hlay
  2. Ramya Balaji
  3. Dimitrije Stanković
  4. Valerie M Christiaens
  5. Carmen Bravo González-Blas
  6. Seppe De Winter
  7. Gert J Hulselmans
  8. Maxime De Waegeneer
  9. Xiaojiang Quan
  10. Duygu Koldere
  11. Mardelle Atkins
  12. Georg Halder
  13. Mirka Uhlirova
  14. Anne Classen  Is a corresponding author
  15. Stein Aerts  Is a corresponding author
  1. VIB-KU Leuven Center for Brain and Disease Research, Belgium
  2. University of Freiburg, Germany
  3. University of Colognee, Germany
  4. KU Leuven, Belgium
  5. Sam Houston State University, United States
  6. VIB-KU Leuven Center for Cancer Biology, Belgium
  7. University of Cologne, Germany

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)

The following data sets were generated
The following previously published data sets were used

Article and author information

Author details

  1. Swann Floc'hlay

    VIB-KU Leuven Center for Brain and Disease Research, Leuven, Belgium
    Competing interests
    The authors declare that no competing interests exist.
  2. Ramya Balaji

    Faculty of Biology, University of Freiburg, Freiburg, Germany
    Competing interests
    The authors declare that no competing interests exist.
  3. Dimitrije Stanković

    Institute for Genetics, University of Colognee, Cologne, Germany
    Competing interests
    The authors declare that no competing interests exist.
  4. Valerie M Christiaens

    Department of Human Genetics, KU Leuven, Leuven, Belgium
    Competing interests
    The authors declare that no competing interests exist.
  5. Carmen Bravo González-Blas

    VIB-KU Leuven Center for Brain and Disease Research, Leuven, Belgium
    Competing interests
    The authors declare that no competing interests exist.
  6. Seppe De Winter

    VIB-KU Leuven Center for Brain and Disease Research, Leuven, Belgium
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7907-1247
  7. Gert J Hulselmans

    Department of Human Genetics, KU Leuven, Leuven, Belgium
    Competing interests
    The authors declare that no competing interests exist.
  8. Maxime De Waegeneer

    VIB-KU Leuven Center for Brain and Disease Research, Leuven, Belgium
    Competing interests
    The authors declare that no competing interests exist.
  9. Xiaojiang Quan

    VIB-KU Leuven Center for Brain and Disease Research, Leuven, Belgium
    Competing interests
    The authors declare that no competing interests exist.
  10. Duygu Koldere

    VIB-KU Leuven Center for Brain and Disease Research, Leuven, Belgium
    Competing interests
    The authors declare that no competing interests exist.
  11. Mardelle Atkins

    Department of Biological Sciences, Sam Houston State University, Texas, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0245-2452
  12. Georg Halder

    VIB-KU Leuven Center for Cancer Biology, Leuven, Belgium
    Competing interests
    The authors declare that no competing interests exist.
  13. Mirka Uhlirova

    Institute for Genetics, University of Cologne, Cologne, Germany
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5735-8287
  14. Anne Classen

    Faculty of Biology, University of Freiburg, Freiburg, Germany
    For correspondence
    anne.classen@zbsa.uni-freiburg.de
    Competing interests
    The authors declare that no competing interests exist.
  15. Stein Aerts

    VIB-KU Leuven Center for Brain and Disease Research, Leuven, Belgium
    For correspondence
    stein.aerts@kuleuven.be
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8006-0315

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

  1. Utpal Banerjee, University of California, Los Angeles, United States

Version history

  1. Received: June 17, 2022
  2. Preprint posted: June 19, 2022 (view preprint)
  3. Accepted: May 2, 2023
  4. Accepted Manuscript published: May 3, 2023 (version 1)
  5. 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

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Swann Floc'hlay
  2. Ramya Balaji
  3. Dimitrije Stanković
  4. Valerie M Christiaens
  5. Carmen Bravo González-Blas
  6. Seppe De Winter
  7. Gert J Hulselmans
  8. Maxime De Waegeneer
  9. Xiaojiang Quan
  10. Duygu Koldere
  11. Mardelle Atkins
  12. Georg Halder
  13. Mirka Uhlirova
  14. Anne Classen
  15. Stein Aerts
(2023)
Shared enhancer gene regulatory networks between wound and oncogenic programs
eLife 12:e81173.
https://doi.org/10.7554/eLife.81173

Share this article

https://doi.org/10.7554/eLife.81173

Further reading

    1. Computational and Systems Biology
    2. Evolutionary Biology
    Kenya Hitomi, Yoichiro Ishii, Bei-Wen Ying
    Research Article

    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.

    1. Computational and Systems Biology
    2. Physics of Living Systems
    Taegon Chung, Iksoo Chang, Sangyeol Kim
    Research Article

    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.