Analysis of long and short enhancers in melanoma cell states
Abstract
Understanding how enhancers drive cell type specificity and efficiently identifying them is essential for the development of innovative therapeutic strategies. In melanoma, the melanocytic (MEL) and the mesenchymal-like (MES) states present themselves with different responses to therapy, making the identification of specific enhancers highly relevant. Using massively parallel reporter assays (MPRA) in a panel of patient-derived melanoma lines (MM lines), we set to identify and decipher melanoma enhancers by first focusing on regions with state specific H3K27 acetylation close to differentially expressed genes. An in-depth evaluation of those regions was then pursued by investigating the activity of overlapping ATAC-seq peaks along with a full tiling of the acetylated regions with 190 bp sequences. Activity was observed in more than 60% of the selected regions and we were able to precisely locate the active enhancers within ATAC-seq peaks. Comparison of sequence content with activity, using the deep learning model DeepMEL2, revealed that AP-1 alone is responsible for the MES enhancer activity. In contrast, SOX10 and MITF both influence MEL enhancer function with SOX10 being required to achieve high levels of activity. Overall, our MPRAs shed light on the relationship between long and short sequences in terms of their sequence content, enhancer activity, and specificity across melanoma cell states.
Data availability
Sequencing data have been deposited in GEO under accession codes GSE180879.Enhancer activity tables for each library is provided as source data.Scripts used for enhancer - barcode assignment, read processing and activity measurement and analysis are provided in the Scripts directory.
-
Analysis of long and short enhancers in melanoma cell statesNCBI Gene Expression Omnibus, GSE180879.
-
TFAP2A ChIP-seq in human primary melanocytesNCBI Gene Expression Omnibus, GSE67555.
-
Decoding the regulatory landscape of melanoma reveals TEADS as regulators of the invasive cell state.NCBI Gene Expression Omnibus, GSE60666.
-
cisTopic: cis-regulatory topic modelling on single-cell ATAC-seq dataNCBI Gene Expression Omnibus, GSE114557.
-
Cross-species analysis of melanoma enhancer logic using deep learningNCBI Gene Expression Omnibus, GSE142238.
-
Single-cell analysis of gene expression variation and phenotype switching in melanomaNCBI Gene Expression Omnibus, GSE134432.
Article and author information
Author details
Funding
H2020 European Research Council (724226_cis-CONTROL)
- Stein Aerts
KU Leuven (C14/18/092)
- Stein Aerts
Fonds Wetenschappelijk Onderzoek (1S03317N)
- Liesbeth Minnoye
Fonds Wetenschappelijk Onderzoek (12J6916N)
- Jonas Demeulemeester
Kom op tegen Kanker
- Jasper Wouters
H2020 European Research Council (724226_cis-CONTROL)
- David Mauduit
H2020 European Research Council (724226_cis-CONTROL)
- Valerie Christiaens
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2021, Mauduit 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
-
- 2,453
- views
-
- 322
- downloads
-
- 20
- 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
-
- Chromosomes and Gene Expression
RNA interference (RNAi) is a conserved pathway that utilizes Argonaute proteins and their associated small RNAs to exert gene regulatory function on complementary transcripts. While the majority of germline-expressed RNAi proteins reside in perinuclear germ granules, it is unknown whether and how RNAi pathways are spatially organized in other cell types. Here, we find that the small RNA biogenesis machinery is spatially and temporally organized during Caenorhabditis elegans embryogenesis. Specifically, the RNAi factor, SIMR-1, forms visible concentrates during mid-embryogenesis that contain an RNA-dependent RNA polymerase, a poly-UG polymerase, and the unloaded nuclear Argonaute protein, NRDE-3. Curiously, coincident with the appearance of the SIMR granules, the small RNAs bound to NRDE-3 switch from predominantly CSR-class 22G-RNAs to ERGO-dependent 22G-RNAs. NRDE-3 binds ERGO-dependent 22G-RNAs in the somatic cells of larvae and adults to silence ERGO-target genes; here we further demonstrate that NRDE-3-bound, CSR-class 22G-RNAs repress transcription in oocytes. Thus, our study defines two separable roles for NRDE-3, targeting germline-expressed genes during oogenesis to promote global transcriptional repression, and switching during embryogenesis to repress recently duplicated genes and retrotransposons in somatic cells, highlighting the plasticity of Argonaute proteins and the need for more precise temporal characterization of Argonaute-small RNA interactions.
-
- Chromosomes and Gene Expression
- Genetics and Genomics
A new method for mapping torsion provides insights into the ways that the genome responds to the torsion generated by RNA polymerase II.