Ex vivo mapping of enhancer networks that define the transcriptional program driving melanoma metastasis.

Mortality from vmelanoma is associated with metastatic disease, but the mechanisms leading to spreading of the cancer cells remain obscure. Spatial profiling revealed that melanoma is characterized by a high degree of heterogeneity, which is established by the ability of melanoma cells to switch between different phenotypical stages. This plasticity, likely a heritage from embryonic pathways, accounts for a relevant part of the metastatic potential of these lesions, and requires the rapid and efficient reorganization of the transcriptional landscape of melanoma cells. A large part of the non-coding genome cooperates to control gene expression, specifically through the activity of enhancers (ENHs). In this study, we aimed to identify ex vivo the network of active ENHs and to outline their cooperative interactions in supporting transcriptional adaptation during melanoma metastatic progression. We conducted a genome-wide analysis to map active ENHs distribution in a retrospective cohort of 39 melanoma patients, comparing the profiles obtained in primary (N = 19) and metastatic (N = 20) melanoma lesions. Unsupervised clustering showed that the profile for acetylated histone H3 at lysine 27 (H3K27ac) efficiently segregates lesions into three different clusters corresponding to progressive stages of the disease. We reconstructed the map of super-ENHs (SEs) and cooperative ENHs that associate with metastatic progression in melanoma, which showed that cooperation among regulatory elements is a mandatory requirement for transcriptional plasticity. We also showed that these elements carry out specialized and non-redundant functions, and indicated the existence of a hierarchical organization, with SEs on top as masterminds of the entire transcriptional program and classical ENHs as executors. By providing an innovative vision of how the chromatin landscape of melanoma works during metastatic spreading, our data also point out the need to integrate functional profiling in the analysis of cancer lesions to increase definition and improve interpretation of tumor heterogeneity.