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.
Reviewing Editor
- Jian Xu, University of Texas Southwestern Medical Center, United States
Publication history
- Received: June 28, 2021
- Preprint posted: July 27, 2021 (view preprint)
- Accepted: December 6, 2021
- Accepted Manuscript published: December 7, 2021 (version 1)
- Version of Record published: December 21, 2021 (version 2)
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
-
- 947
- Page views
-
- 158
- Downloads
-
- 2
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
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)
Further reading
-
- Chromosomes and Gene Expression
- Microbiology and Infectious Disease
Mosquitoes transmit numerous pathogens, but large gaps remain in our understanding of their physiology. To facilitate explorations of mosquito biology, we have created Aegypti-Atlas (http://aegyptiatlas.buchonlab.com/), an online resource hosting RNAseq profiles of Ae. aegypti body parts (head, thorax, abdomen, gut, Malpighian tubules, ovaries), gut regions (crop, proventriculus, anterior and posterior midgut, hindgut), and a gut time course of blood meal digestion. Using Aegypti-Atlas, we provide insights into regionalization of gut function, blood feeding response, and immune defenses. We find that the anterior and posterior midgut possess digestive specializations which are preserved in the blood-fed state. Blood feeding initiates the sequential induction and repression/depletion of multiple cohorts of peptidases. With respect to defense, immune signaling components, but not recognition or effector molecules, show enrichment in ovaries. Basal expression of antimicrobial peptides is dominated by holotricin and gambicin, which are expressed in carcass and digestive tissues, respectively, in a mutually exclusive manner. In the midgut, gambicin and other effectors are almost exclusively expressed in the anterior regions, while the posterior midgut exhibits hallmarks of immune tolerance. Finally, in a cross-species comparison between Ae. aegypti and Anopheles gambiae midguts, we observe that regional digestive and immune specializations are conserved, indicating that our dataset may be broadly relevant to multiple mosquito species. We demonstrate that the expression of orthologous genes is highly correlated, with the exception of a ‘species signature’ comprising a few highly/disparately expressed genes. With this work, we show the potential of Aegypti-Atlas to unlock a more complete understanding of mosquito biology.
-
- Chromosomes and Gene Expression
- Genetics and Genomics
For some inducible genes, the rate and molecular mechanism of transcriptional activation depends on the prior experiences of the cell. This phenomenon, called epigenetic transcriptional memory, accelerates reactivation and requires both changes in chromatin structure and recruitment of poised RNA Polymerase II (RNAPII) to the promoter. Memory of inositol starvation in budding yeast involves a positive feedback loop between transcription factor-dependent interaction with the nuclear pore complex and histone H3 lysine 4 dimethylation (H3K4me2). While H3K4me2 is essential for recruitment of RNAPII and faster reactivation, RNAPII is not required for H3K4me2. Unlike RNAPII-dependent H3K4me2 associated with transcription, RNAPII-independent H3K4me2 requires Nup100, SET3C, the Leo1 subunit of the Paf1 complex and, upon degradation of an essential transcription factor, is inherited through multiple cell cycles. The writer of this mark (COMPASS) physically interacts with the potential reader (SET3C), suggesting a molecular mechanism for the spreading and re-incorporation of H3K4me2 following DNA replication.