Trackplot: a fast and lightweight R script for epigenomic enrichment plots.

Motivation: BigWig files serve as essential inputs in epigenomic data visualization. However, current R packages for visualizing these files are limited, slow, and burdened by numerous dependencies. Results: We introduce trackplot, a minimal R script designed for the rapid generation of integrative genomics viewer (IGV) style track plots, profile plots, and heatmaps from bigWig files. This script offers speed, owing to its reliance on bwtool, resulting in performance gains of several magnitudes compared to equivalent packages. The script is lightweight, requiring only the data.table and bwtool packages as primary dependencies. Notably, the plots are generated in base R graphics, eliminating the need for additional packages. trackplot queries the University of California Santa Cruz (UCSC) genome browser for gene models thereby enhancing the reproducibility of analyses. The script extends its support to general transfer format (GTF) further enhancing its versatility. This tool addresses the gaps in existing bigWig visualization approaches by offering speed, simplicity, and minimal dependencies, thereby presenting a valuable asset to researchers in the fields of epigenomics. Availability and implementation: trackplot is implemented in R is made available under MIT license at https://github.com/PoisonAlien/trackplot.

The GAPDH redox switch safeguards reductive capacity and enables survival of stressed tumour cells.

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is known to contain an active-site cysteine residue undergoing oxidation in response to hydrogen peroxide, leading to rapid inactivation of the enzyme. Here we show that human and mouse cells expressing a GAPDH mutant lacking this redox switch retain catalytic activity but are unable to stimulate the oxidative pentose phosphate pathway and enhance their reductive capacity. Specifically, we find that anchorage-independent growth of cells and spheroids is limited by an elevation of endogenous peroxide levels and is largely dependent on a functional GAPDH redox switch. Likewise, tumour growth in vivo is limited by peroxide stress and suppressed when the GAPDH redox switch is disabled in tumour cells. The induction of additional intratumoural oxidative stress by chemo- or radiotherapy synergized with the deactivation of the GAPDH redox switch. Mice lacking the GAPDH redox switch exhibit altered fatty acid metabolism in kidney and heart, apparently in compensation for the lack of the redox switch. Together, our findings demonstrate the physiological and pathophysiological relevance of oxidative GAPDH inactivation in mammals.

TAL1 activation in T-cell acute lymphoblastic leukemia: a novel oncogenic 3' neo-enhancer.

T-cell acute lymphocytic leukemia protein 1 (TAL1) is one of the most frequently deregulated oncogenes in T-cell acute lymphoblastic leukemia (T-ALL). Its deregulation can occur through diverse cis-alterations, including SIL-TAL1 microdeletions, translocations with T-cell Receptor loci, and more recently described upstream intergenic non-coding mutations. These mutations consist of recurrent focal microinsertions that create an oncogenic neo-enhancer accompanied by activating epigenetic marks. This observation laid the groundwork for an innovative paradigm concerning the activation of proto-oncogenes via genomic alterations of non-coding intergenic regions. However, for the majority of T-ALL expressing TAL1 (TAL1+), the deregulation mechanism remains 'unresolved'. We took advantage of H3K27ac and H3K4me3 chromatin immunoprecipitation sequencing data of eight cases of T-ALL, including five TAL1+ cases. We identified a putative novel oncogenic neo-enhancer downstream of TAL1 in an unresolved monoallelic TAL1+ case. A rare but recurrent somatic heterozygous microinsertion within this region creates a de novo binding site for MYB transcription factor. Here we demonstrate that this mutation leads to increased enhancer activity, gain of active epigenetic marks, and TAL1 activation via recruitment of MYB. These results highlight the diversity of non-coding mutations that can drive oncogene activation.

Transcriptome analysis identifies TODL as a novel lncRNA associated with proliferation, differentiation, and tumorigenesis in liposarcoma through FOXM1.

Liposarcoma, the most common soft tissue sarcoma, is a group of fat cell mesenchymal tumors with different histological subtypes. The dysregulation of long non-coding RNAs (lncRNAs) has been observed in human cancers including a few studies in sarcoma. However, the global transcriptome analysis and potential role of lncRNAs remain unexplored in liposarcoma. The present investigation uncovers the transcriptomic profile of liposarcoma by RNA sequencing to gain insight into the global transcriptional changes in liposarcoma. Our RNA sequencing analysis has identified that many oncogenic lncRNAs are differentially expressed in different subtypes of liposarcoma including MALAT1, PVT1, SNHG15, LINC00152, and MIR210HG. Importantly, we identified a highly overexpressed, unannotated, and novel lncRNA in dedifferentiated liposarcomas. We have named it TODL, transcript overexpressed in dedifferentiated liposarcoma. TODL lncRNA displayed significantly higher expression in dedifferentiated liposarcoma cell lines and patient samples. Interestingly, functional studies revealed that TODL lncRNA has an oncogenic function in liposarcoma cells by regulating proliferation, cell cycle, apoptosis, differentiation, and tumorigenesis in the murine model. Silencing of TODL lncRNA highlighted the enrichment of several key oncogenic signaling pathways including cell cycle, transcriptional misregulation, FOXM1 network, p53 signaling, PLK1 signaling, FoxO, and signaling Aurora signaling pathways. RNA pull-down assay revealed the binding of TODL lncRNA with FOXM1, an oncogenic transcription factor, and the key regulator of the cell cycle. Silencing of TODL lncRNA also induces adipogenesis in dedifferentiated liposarcomas. Altogether, our finding indicates that TODL could be utilized as a novel, specific diagnostic biomarker, and a pharmacological target for therapeutic development in controlling aggressive and metastatic dedifferentiated liposarcomas.

Genetic and Methylation Analysis of CTNNB1 in Benign and Malignant Melanocytic Lesions.

Melanocytic neoplasms have been genetically characterized in detail during the last decade. Recurrent CTNNB1 exon 3 mutations have been recognized in the distinct group of melanocytic tumors showing deep penetrating nevus-like morphology. In addition, they have been identified in 1-2% of advanced melanoma. Performing a detailed genetic analysis of difficult-to-classify nevi and melanomas with CTNNB1 mutations, we found that benign tumors (nevi) show characteristic morphological, genetic and epigenetic traits, which distinguish them from other nevi and melanoma. Malignant CTNNB1-mutant tumors (melanomas) demonstrated a different genetic profile, instead grouping clearly with other non-CTNNB1 melanomas in methylation assays. To further evaluate the role of CTNNB1 mutations in melanoma, we assessed a large cohort of clinically sequenced melanomas, identifying 38 tumors with CTNNB1 exon 3 mutations, including recurrent S45 (n = 13, 34%), G34 (n = 5, 13%), and S27 (n = 5, 13%) mutations. Locations and histological subtype of CTNNB1-mutated melanoma varied; none were reported as showing deep penetrating nevus-like morphology. The most frequent concurrent activating mutations were BRAF V600 (n = 21, 55%) and NRAS Q61 (n = 13, 34%). In our cohort, four of seven (58%) and one of nine (11%) patients treated with targeted therapy (BRAF and MEK Inhibitors) or immune-checkpoint therapy, respectively, showed disease control (partial response or stable disease). In summary, CTNNB1 mutations are associated with a unique melanocytic tumor type in benign tumors (nevi), which can be applied in a diagnostic setting. In advanced disease, no clear characteristics distinguishing CTNNB1-mutant from other melanomas were observed; however, studies of larger, optimally prospective, cohorts are warranted.

Genetic and methylation profiles distinguish benign, malignant and spitzoid melanocytic tumors.

Accurate classification of melanocytic tumors is important for prognostic evaluation, treatment and follow-up protocols of patients. The majority of melanocytic proliferations can be classified solely based on clinical and pathological criteria, however in select cases a definitive diagnostic assessment remains challenging and additional diagnostic biomarkers would be advantageous. We analyzed melanomas, nevi, Spitz nevi and atypical spitzoid tumors using parallel sequencing (exons of 611 genes and 507 gene translocation analysis) and methylation arrays (850k Illumina EPIC). By combining detailed genetic and epigenetic analysis with reference-based and reference-free DNA methylome deconvolution we compared Spitz nevi to nevi and melanoma and assessed the potential for these methods in classifying challenging spitzoid tumors. Results were correlated with clinical and histologic features. Spitz nevi were found to cluster independently of nevi and melanoma and demonstrated a different mutation profile. Multiple copy number alterations and TERT promoter mutations were identified only in melanomas. Genome-wide methylation in Spitz nevi was comparable to benign nevi while the Leukocytes UnMethylation for Purity (LUMP) algorithm in Spitz nevi was comparable to melanoma. Histologically difficult to classify Spitz tumor cases were assessed which, based on methylation arrays, clustered between Spitz nevi and melanoma and in terms of genetic profile or copy number variations demonstrated worrisome features suggesting a malignant neoplasm. Comprehensive sequencing and methylation analysis verify Spitz nevi as an independent melanocytic entity distinct from both nevi and melanoma. Combined genetic and methylation assays can offer additional insights in diagnosing difficult to classify Spitzoid tumors.

Histone H3K27 demethylase KDM6A is an epigenetic gatekeeper of mTORC1 signalling in cancer.

OBJECTIVE: Large-scale genome sequencing efforts of human tumours identified epigenetic modifiers as one of the most frequently mutated gene class in human cancer. However, how these mutations drive tumour development and tumour progression are largely unknown. Here, we investigated the function of the histone demethylase KDM6A in gastrointestinal cancers, such as liver cancer and pancreatic cancer. DESIGN: Genetic alterations as well as expression analyses of KDM6A were performed in patients with liver cancer. Genetic mouse models of liver and pancreatic cancer coupled with Kdm6a-deficiency were investigated, transcriptomic and epigenetic profiling was performed, and in vivo and in vitro drug treatments were conducted. RESULTS: KDM6A expression was lost in 30% of patients with liver cancer. Kdm6a deletion significantly accelerated tumour development in murine liver and pancreatic cancer models. Kdm6a-deficient tumours showed hyperactivation of mTORC1 signalling, whereas endogenous Kdm6a re-expression by inducible RNA-interference in established Kdm6a-deficient tumours diminished mTORC1 activity resulting in attenuated tumour progression. Genome-wide transcriptional and epigenetic profiling revealed direct binding of Kdm6a to crucial negative regulators of mTORC1, such as Deptor, and subsequent transcriptional activation by epigenetic remodelling. Moreover, in vitro and in vivo genetic epistasis experiments illustrated a crucial function of Deptor and mTORC1 in Kdm6a-dependent tumour suppression. Importantly, KDM6A expression in human tumours correlates with mTORC1 activity and KDM6A-deficient tumours exhibit increased sensitivity to mTORC1 inhibition. CONCLUSION: KDM6A is an important tumour suppressor in gastrointestinal cancers and acts as an epigenetic toggle for mTORC1 signalling. Patients with KDM6A-deficient tumours could benefit of targeted therapy focusing on mTORC1 inhibition.

Epigenetic Modulation of Radiation-Induced Diacylglycerol Kinase Alpha Expression Prevents Pro-Fibrotic Fibroblast Response.

Radiotherapy, a common component in cancer treatment, can induce adverse effects including fibrosis in co-irradiated tissues. We previously showed that differential DNA methylation at an enhancer of diacylglycerol kinase alpha (DGKA) in normal dermal fibroblasts is associated with radiation-induced fibrosis. After irradiation, the transcription factor EGR1 is induced and binds to the hypomethylated enhancer, leading to increased DGKA and pro-fibrotic marker expression. We now modulated this DGKA induction by targeted epigenomic and genomic editing of the DGKA enhancer and administering epigenetic drugs. Targeted DNA demethylation of the DGKA enhancer in HEK293T cells resulted in enrichment of enhancer-related histone activation marks and radiation-induced DGKA expression. Mutations of the EGR1-binding motifs decreased radiation-induced DGKA expression in BJ fibroblasts and caused dysregulation of multiple fibrosis-related pathways. EZH2 inhibitors (GSK126, EPZ6438) did not change radiation-induced DGKA increase. Bromodomain inhibitors (CBP30, JQ1) suppressed radiation-induced DGKA and pro-fibrotic marker expression. Similar drug effects were observed in donor-derived fibroblasts with low DNA methylation. Overall, epigenomic manipulation of DGKA expression may offer novel options for a personalized treatment to prevent or attenuate radiotherapy-induced fibrosis.

Epigenetic analysis of patients with T-ALL identifies poor outcomes and a hypomethylating agent-responsive subgroup.

Adult "T cell" acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy that is associated with poor outcomes, requiring additional therapeutic options. The DNA methylation landscapes of adult T-ALL remain undercharacterized. Here, we systematically analyzed the DNA methylation profiles of normal thymic-sorted T cell subpopulations and 143 primary adult T-ALLs as part of the French GRAALL 2003-2005 trial. Our results indicated that T-ALL is epigenetically heterogeneous consisting of five subtypes (C1-C5), which were either associated with co-occurring DNA methyltransferase 3 alpha (DNMT3A)/isocitrate dehydrogenase [NADP(+)] 2 (IDH2) mutations (C1), TAL bHLH transcription factor 1, erythroid differentiation factor (TAL1) deregulation (C2), T cell leukemia homeobox 3 (TLX3) (C3), TLX1/in cis-homeobox A9 (HOXA9) (C4), or in trans-HOXA9 overexpression (C5). Integrative analysis of DNA methylation and gene expression identified potential cluster-specific oncogenes and tumor suppressor genes. In addition to an aggressive hypomethylated subgroup (C1), our data identified an unexpected subset of hypermethylated T-ALL (C5) associated with poor outcome and primary therapeutic response. Using mouse xenografts, we demonstrated that hypermethylated T-ALL samples exhibited therapeutic responses to the DNA hypomethylating agent 5-azacytidine, which significantly (survival probability; P = 0.001 for C3, 0.01 for C4, and 0.0253 for C5) delayed tumor progression. These findings suggest that epigenetic-based therapies may provide an alternative treatment option in hypermethylated T-ALL.

Topography of transcriptionally active chromatin in glioblastoma.

Molecular profiling of the most aggressive brain tumor glioblastoma (GBM) on the basis of gene expression, DNA methylation, and genomic variations advances both cancer research and clinical diagnosis. The enhancer architectures and regulatory circuitries governing tumor-intrinsic transcriptional diversity and subtype identity are still elusive. Here, by mapping H3K27ac deposition, we analyze the active regulatory landscapes across 95 GBM biopsies, 12 normal brain tissues, and 38 cell line counterparts. Analyses of differentially regulated enhancers and super-enhancers uncovered previously unrecognized layers of intertumor heterogeneity. Integrative analysis of variant enhancer loci and transcriptome identified topographies of transcriptional enhancers and core regulatory circuitries in four molecular subtypes of primary tumors: AC1-mesenchymal, AC1-classical, AC2-proneural, and AC3-proneural. Moreover, this study reveals core oncogenic dependency on super-enhancer-driven transcriptional factors, long noncoding RNAs, and druggable targets in GBM. Through profiling of transcriptional enhancers, we provide clinically relevant insights into molecular classification, pathogenesis, and therapeutic intervention of GBM.

Methrix: an R/Bioconductor package for systematic aggregation and analysis of bisulfite sequencing data.

MOTIVATION: Whole-genome bisulfite sequencing (WGBS) measures DNA methylation at base pair resolution resulting in large bedGraph like coverage files. Current options for processing such files are hindered by discrepancies in file format specification, speed, and memory requirements. RESULTS: We developed methrix, an R package, which provides a toolset for systematic analysis of large datasets. Core functionality of the package includes a comprehensive bedGraph or similar tab-separated text file reader-which summarizes methylation calls based on annotated reference indices, infers and collapses strands and handles uncovered reference CpG sites while facilitating a flexible input file format specification. Additional optimized functions for quality control filtering, subsetting and visualization allow user-friendly and effective processing of WGBS results. Easy integration with tools for differentially methylated region (DMR) calling and annotation further eases the analysis of genome-wide methylation data. Overall, methrix enriches established WGBS workflows by bringing together computational efficiency and versatile functionality. AVAILABILITY AND IMPLEMENTATION: Methrix is implemented as an R package, made available under MIT license at https://github.com/CompEpigen/methrix and can be installed from the Bioconductor repository. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Globally altered epigenetic landscape and delayed osteogenic differentiation in H3.3-G34W-mutant giant cell tumor of bone.

The neoplastic stromal cells of giant cell tumor of bone (GCTB) carry a mutation in H3F3A, leading to a mutant histone variant, H3.3-G34W, as a sole recurrent genetic alteration. We show that in patient-derived stromal cells H3.3-G34W is incorporated into the chromatin and associates with massive epigenetic alterations on the DNA methylation, chromatin accessibility and histone modification level, that can be partially recapitulated in an orthogonal cell line system by the introduction of H3.3-G34W. These epigenetic alterations affect mainly heterochromatic and bivalent regions and provide possible explanations for the genomic instability, as well as the osteolytic phenotype of GCTB. The mutation occurs in differentiating mesenchymal stem cells and associates with an impaired osteogenic differentiation. We propose that the observed epigenetic alterations reflect distinct differentiation stages of H3.3 WT and H3.3 MUT stromal cells and add to H3.3-G34W-associated changes.

TP63, SOX2, and KLF5 Establish a Core Regulatory Circuitry That Controls Epigenetic and Transcription Patterns in Esophageal Squamous Cell Carcinoma Cell Lines.

BACKGROUND & AIMS: We investigated the transcriptome of esophageal squamous cell carcinoma (ESCC) cells, activity of gene regulatory (enhancer and promoter regions), and the effects of blocking epigenetic regulatory proteins. METHODS: We performed chromatin immunoprecipitation sequencing with antibodies against H3K4me1, H3K4me3, and H3K27ac and an assay for transposase-accessible chromatin to map the enhancer regions and accessible chromatin in 8 ESCC cell lines. We used the CRC_Mapper algorithm to identify core regulatory circuitry transcription factors in ESCC cell lines, and determined genome occupancy profiles for 3 of these factors. In ESCC cell lines, expression of transcription factors was knocked down with small hairpin RNAs, promoter and enhancer regions were disrupted by CRISPR/Cas9 genome editing, or bromodomains and extraterminal (BET) family proteins and histone deacetylases (HDACs) were inhibited with ARV-771 and romidepsin, respectively. ESCC cell lines were then analyzed by whole-transcriptome sequencing, immunoprecipitation, immunoblots, immunohistochemistry, and viability assays. Interactions between distal enhancers and promoters were identified and verified with circular chromosome conformation capture sequencing. NOD-SCID mice were given injections of modified ESCC cells, some mice where given injections of HDAC or BET inhibitors, and growth of xenograft tumors was measured. RESULTS: We identified super-enhancer-regulated circuits and transcription factors TP63, SOX2, and KLF5 as core regulatory factors in ESCC cells. Super-enhancer regulation of ALDH3A1 mediated by core regulatory factors was required for ESCC viability. We observed direct interactions between the promoter region of TP63 and functional enhancers, mediated by the core regulatory circuitry transcription factors. Deletion of enhancer regions from ESCC cells decreased expression of the core regulatory circuitry transcription factors and reduced cell viability; these same results were observed with knockdown of each core regulatory circuitry transcription factor. Incubation of ESCC cells with BET and HDAC disrupted the core regulatory circuitry program and the epigenetic modifications observed in these cells; mice given injections of HDAC or BET inhibitors developed smaller xenograft tumors from the ESCC cell lines. Xenograft tumors grew more slowly in mice given the combination of ARV-771 and romidepsin than mice given either agent alone. CONCLUSIONS: In epigenetic and transcriptional analyses of ESCC cell lines, we found the transcription factors TP63, SOX2, and KLF5 to be part of a core regulatory network that determines chromatin accessibility, epigenetic modifications, and gene expression patterns in these cells. A combination of epigenetic inhibitors slowed growth of xenograft tumors derived from ESCC cells in mice.

RNA-Binding Protein ZFP36L1 Suppresses Hypoxia and Cell-Cycle Signaling.

ZFP36L1 is a tandem zinc-finger RNA-binding protein that recognizes conserved adenylate-uridylate-rich elements (ARE) located in 3'untranslated regions (UTR) to mediate mRNA decay. We hypothesized that ZFP36L1 is a negative regulator of a posttranscriptional hub involved in mRNA half-life regulation of cancer-related transcripts. Analysis of in silico data revealed that ZFP36L1 was significantly mutated, epigenetically silenced, and downregulated in a variety of cancers. Forced expression of ZFP36L1 in cancer cells markedly reduced cell proliferation in vitro and in vivo, whereas silencing of ZFP36L1 enhanced tumor cell growth. To identify direct downstream targets of ZFP36L1, systematic screening using RNA pull-down of wild-type and mutant ZFP36L1 as well as whole transcriptome sequencing of bladder cancer cells {plus minus} tet-on ZFP36L1 was performed. A network of 1,410 genes was identified as potential direct targets of ZFP36L1. These targets included a number of key oncogenic transcripts such as HIF1A, CCND1, and E2F1. ZFP36L1 specifically bound to the 3'UTRs of these targets for mRNA degradation, thus suppressing their expression. Dual luciferase reporter assays and RNA electrophoretic mobility shift assays showed that wild-type, but not zinc-finger mutant ZFP36L1, bound to HIF1A 3'UTR and mediated HIF1A mRNA degradation, leading to reduced expression of HIF1A and its downstream targets. Collectively, our findings reveal an indispensable role of ZFP36L1 as a posttranscriptional safeguard against aberrant hypoxic signaling and abnormal cell-cycle progression. SIGNIFICANCE: RNA-binding protein ZFP36L1 functions as a tumor suppressor by regulating the mRNA stability of a number of mRNAs involved in hypoxia and cell-cycle signaling.

Identification of somatic alterations in lipoma using whole exome sequencing.

Lipomas are benign fatty tumors with a high prevalence rate, mostly found in adults but have a good prognosis. Until now, reason for lipoma occurrence not been identified. We performed whole exome sequencing to define the mutational spectrum in ten lipoma patients along with their matching control samples. We presented genomic insight into the development of lipomas, the most common benign tumor of soft tissue. Our analysis identified 412 somatic variants including missense mutations, splice site variants, frameshift indels, and stop gain/lost. Copy number variation analysis highlighted minor aberrations in patients. Kinase genes and transcriptions factors were among the validated mutated genes critical for cell proliferation and survival. Pathway analysis revealed enrichment of calcium, Wnt and phospholipase D signaling in patients. In conclusion, whole exome sequencing in lipomas identified mutations in genes with a possible role in development and progression of lipomas.

Identification of a novel enhancer of CEBPE essential for granulocytic differentiation.

CCAAT/enhancer binding protein ε (CEBPE) is an essential transcription factor for granulocytic differentiation. Mutations of CEBPE occur in individuals with neutrophil-specific granule deficiency (SGD), which is characterized by defects in neutrophil maturation. Cebpe-knockout mice also exhibit defects in terminal differentiation of granulocytes, a phenotype reminiscent of SGD. Analysis of DNase I hypersensitive sites sequencing data revealed an open chromatin region 6 kb downstream of the transcriptional start site of Cebpe in murine myeloid cells. We identified an interaction between this +6-kb region and the core promoter of Cebpe using circular chromosome conformation capture sequencing (4C-seq). To understand the role of this putative enhancer in transcriptional regulation of Cebpe, we targeted it using catalytically inactive Cas9 fused to Krüppel-associated box (KRAB) domain and observed a significant downregulation of transcript and protein levels of CEBPE in cells expressing guide RNA targeting the +6-kb region. To further investigate the role of this novel enhancer further in myelopoiesis, we generated mice with deletion of this region using CRISPR/Cas9 technology. Germline deletion of the +6-kb enhancer resulted in reduced levels of CEBPE and its target genes and caused a severe block in granulocytic differentiation. We also identified binding of CEBPA and CEBPE to the +6-kb enhancer, which suggests their role in regulating the expression of Cebpe In summary, we have identified a novel enhancer crucial for regulating expression of Cebpe and required for normal granulocytic differentiation.

Chromatin remodeling mediated by ARID1A is indispensable for normal hematopoiesis in mice.

Precise regulation of chromatin architecture is vital to physiological processes including hematopoiesis. ARID1A is a core component of the mammalian SWI/SNF complex, which is one of the ATP-dependent chromatin remodeling complexes. To uncover the role of ARID1A in hematopoietic development, we utilized hematopoietic cell-specific deletion of Arid1a in mice. We demonstrate that ARID1A is essential for maintaining the frequency and function of hematopoietic stem cells and its loss impairs the differentiation of both myeloid and lymphoid lineages. ARID1A deficiency led to a global reduction in open chromatin and ensuing transcriptional changes affected key genes involved in hematopoietic development. We also observed that silencing of ARID1A affected ATRA-induced differentiation of NB4 cells, suggesting its role in granulocytic differentiation of human leukemic cells. Overall, our study provides a comprehensive elucidation of the function of ARID1A in hematopoiesis and highlights the central role of ARID1A-containing SWI/SNF complex in maintaining chromatin dynamics in hematopoietic cells.

Bromodomain and extraterminal proteins foster the core transcriptional regulatory programs and confer vulnerability in liposarcoma.

Liposarcomas (LPSs) are a group of malignant mesenchymal tumors showing adipocytic differentiation. Here, to gain insight into the enhancer dysregulation and transcriptional addiction in this disease, we chart super-enhancer structures in both LPS tissues and cell lines. We identify a bromodomain and extraterminal (BET) protein-cooperated FUS-DDIT3 function in myxoid LPS and a BET protein-dependent core transcriptional regulatory circuitry consisting of FOSL2, MYC, and RUNX1 in de-differentiated LPS. Additionally, SNAI2 is identified as a crucial downstream target that enforces both proliferative and metastatic potentials to de-differentiated LPS cells. Genetic depletion of BET genes, core transcriptional factors, or SNAI2 mitigates consistently LPS malignancy. We also reveal a compelling susceptibility of LPS cells to BET protein degrader ARV-825. BET protein depletion confers additional advantages to circumvent acquired resistance to Trabectedin, a chemotherapy drug for LPS. Moreover, this study provides a framework for discovering and targeting of core oncogenic transcriptional programs in human cancers.

Super-enhancer-associated MEIS1 promotes transcriptional dysregulation in Ewing sarcoma in co-operation with EWS-FLI1.

As the second most common malignant bone tumor in children and adolescents, Ewing sarcoma is initiated and exacerbated by a chimeric oncoprotein, most commonly, EWS-FLI1. In this study, we apply epigenomic analysis to characterize the transcription dysregulation in this cancer, focusing on the investigation of super-enhancer and its associated transcriptional regulatory mechanisms. We demonstrate that super-enhancer-associated transcripts are significantly enriched in EWS-FLI1 target genes, contribute to the aberrant transcriptional network of the disease, and mediate the exceptional sensitivity of Ewing sarcoma to transcriptional inhibition. Through integrative analysis, we identify MEIS1 as a super-enhancer-driven oncogene, which co-operates with EWS-FLI1 in transcriptional regulation, and plays a key pro-survival role in Ewing sarcoma. Moreover, APCDD1, another super-enhancer-associated gene, acting as a downstream target of both MEIS1 and EWS-FLI1, is also characterized as a novel tumor-promoting factor in this malignancy. These data delineate super-enhancer-mediated transcriptional deregulation in Ewing sarcoma, and uncover numerous candidate oncogenes which can be exploited for further understanding of the molecular pathogenesis for this disease.