Single Nucleotide Variants (SNVs), including somatic point mutations and Single Nucleotide Polymorphisms (SNPs), in noncoding cis-regulatory elements (CREs) can affect gene regulation and lead to ...disease development. Several approaches have been developed to identify highly mutated regions, but these do not take into account the specific genomic context, and thus likelihood of mutation, of CREs.
Here, we present SMuRF (Significantly Mutated Region Finder), a user-friendly command-line tool to identify these significantly mutated regions from user-defined genomic intervals and SNVs. We demonstrate this using publicly available datasets in which SMuRF identifies 72 significantly mutated CREs in liver cancer, including known mutated gene promoters as well as previously unreported regions.
SMuRF is a helpful tool to allow the simple identification of significantly mutated regulatory elements. It is open-source and freely available on GitHub ( https://github.com/LupienLab/SMURF ).
Glioblastomas exhibit a hierarchical cellular organization, suggesting that they are driven by neoplastic stem cells that retain partial yet abnormal differentiation potential. Here, we show that a ...large subset of patient-derived glioblastoma stem cells (GSCs) express high levels of Achaete-scute homolog 1 (ASCL1), a proneural transcription factor involved in normal neurogenesis. ASCL1hi GSCs exhibit a latent capacity for terminal neuronal differentiation in response to inhibition of Notch signaling, whereas ASCL1lo GSCs do not. Increasing ASCL1 levels in ASCL1lo GSCs restores neuronal lineage potential, promotes terminal differentiation, and attenuates tumorigenicity. ASCL1 mediates these effects by functioning as a pioneer factor at closed chromatin, opening new sites to activate a neurogenic gene expression program. Directing GSCs toward terminal differentiation may provide therapeutic applications for a subset of GBM patients and strongly supports efforts to restore differentiation potential in GBM and other cancers.
Display omitted
•GSCs can be classified into two subgroups based on ASCL1 expression•ASCL1 is required for GSCs to undergo neuronal lineage differentiation•ASCL1hi GSCs are responsive to Notch pathway inhibitors•ASCL1 binds closed chromatin to activate target genes that drive neuronal fate
Glioblastoma is characterized by a block in cellular differentiation. Park et al. identify a subset of glioblastoma stem cells (GSCs) that express high levels of the proneural transcription factor ASCL1 and differentiate in response to Notch inhibition, effectively abrogating their stemness properties and tumorigenic potential.
Thousands of noncoding somatic single-nucleotide variants (SNVs) of unknown function are reported in tumors. Partitioning the genome according to cistromes reveals the enrichment of somatic SNVs in ...prostate tumors as opposed to adjacent normal tissue cistromes of master transcription regulators, including AR, FOXA1, and HOXB13. This parallels enrichment of prostate cancer genetic predispositions over these transcription regulators' tumor cistromes, exemplified at the 8q24 locus harboring both risk variants and somatic SNVs in cis-regulatory elements upregulating MYC expression. However, Massively Parallel Reporter Assays reveal that few SNVs can alter the transactivation potential of individual cis-regulatory elements. Instead, similar to inherited risk variants, SNVs accumulate in cistromes of master transcription regulators required for prostate cancer development.
Display omitted
•Cistrome partitioning reveals convergence of somatic and risk-variant SNVs•Somatic SNVs enrich within the tumor-defined as opposed to normal-defined cistromes•Somatic SNVs can modulate the output from CREs, such as seen at the MYC locus•Few SNVs alter CREs but collectively target essential transcription regulators
Mazrooei et al. show that somatic SNVs are enriched in cistromes of master transcription regulators in prostate cancers but not in adjacent normal tissues. Few SNVs alter CREs but collectively target essential transcription regulators.
Isocitrate dehydrogenase (IDH) genes 1 and 2 are frequently mutated in acute myeloid leukaemia (AML), low-grade glioma, cholangiocarcinoma (CC) and chondrosarcoma (CS). For AML, low-grade glioma and ...CC, mutant IDH status is associated with a DNA hypermethylation phenotype, implicating altered epigenome dynamics in the aetiology of these cancers. Here we show that the IDH variants in CS are also associated with a hypermethylation phenotype and display increased production of the oncometabolite 2-hydroxyglutarate, supporting the role of mutant IDH-produced 2-hydroxyglutarate as an inhibitor of TET-mediated DNA demethylation. Meta-analysis of the acute myeloid leukaemia, low-grade glioma, cholangiocarcinoma and CS methylation data identifies cancer-specific effectors within the retinoic acid receptor activation pathway among the hypermethylated targets. By analysing sequence motifs surrounding hypermethylated sites across the four cancer types, and using chromatin immunoprecipitation and western blotting, we identify the transcription factor EBF1 (early B-cell factor 1) as an interaction partner for TET2, suggesting a sequence-specific mechanism for regulating DNA methylation.
We investigated the role of 3D genome architecture in instructing functional properties of glioblastoma stem cells (GSCs) by generating sub-5-kb resolution 3D genome maps by in situ Hi-C. Contact ...maps at sub-5-kb resolution allow identification of individual DNA loops, domain organization, and large-scale genome compartmentalization. We observed differences in looping architectures among GSCs from different patients, suggesting that 3D genome architecture is a further layer of inter-patient heterogeneity for glioblastoma. Integration of DNA contact maps with chromatin and transcriptional profiles identified specific mechanisms of gene regulation, including the convergence of multiple super enhancers to individual stemness genes within individual cells. We show that the number of loops contacting a gene correlates with elevated transcription. These results indicate that stemness genes are hubs of interaction between multiple regulatory regions, likely to ensure their sustained expression. Regions of open chromatin common among the GSCs tested were poised for expression of immune-related genes, including
We demonstrate that this gene is co-expressed with stemness genes in GSCs and that CD276 can be targeted with an antibody-drug conjugate to eliminate self-renewing cells. Our results demonstrate that integrated structural genomics data sets can be employed to rationally identify therapeutic vulnerabilities in self-renewing cells.
We describe molecular convergence between BMI1 and CHD7 in the initiation of medulloblastoma. Identified in a functional genomic screen in mouse models, a BMI1High;CHD7Low expression signature within ...medulloblastoma characterizes patients with poor overall survival. We show that BMI1-mediated repression of the ERK1/2 pathway leads to increased proliferation and tumor burden in primary human MB cells and in a xenograft model, respectively. We provide evidence that repression of the ERK inhibitor DUSP4 by BMI1 is dependent on a more accessible chromatin configuration in G4 MB cells with low CHD7 expression. These findings extend current knowledge of the role of BMI1 and CHD7 in medulloblastoma pathogenesis, and they raise the possibility that pharmacological targeting of BMI1 or ERK may be particularly indicated in a subgroup of MB with low expression levels of CHD7.
Display omitted
•CHD7 inactivation induces MB in a forward genetic screen in Bmi1-overexpressing mice•A BMI1High;CHD7Low signature is found in MB patients with poor overall survival•CHD7Low favors chromatin accessibility at PcG target genes•ERK overactivity increases proliferation and tumor burden in BMI1High;CHD7Low MB
Badodi et al. find convergence of the chromatin modifiers BMI1 and CHD7 in medulloblastoma pathogenesis, and they show that this pathway regulates tumor proliferation and growth via ERK signaling.
Abstract
Motivation
The 3D genome architecture influences the regulation of genes by facilitating chromatin interactions between distal cis-regulatory elements and gene promoters. We implement Cross ...Cell-type Correlation based on DNA accessibility (C3D), a customizable computational tool that predicts chromatin interactions using an unsupervised algorithm that utilizes correlations in chromatin measurements, such as DNaseI hypersensitivity signals.
Results
C3D accurately predicts 32.7%, 18.3% and 24.1% of interactions, validated by ChIA-PET assays, between promoters and distal regions that overlie DNaseI hypersensitive sites in K562, MCF-7 and GM12878 cells, respectively.
Availability and implementation
Source code is open-source and freely available on GitHub (https://github.com/LupienLabOrganization/C3D) under the GNU GPLv3 license. C3D is implemented in Bash and R; it runs on any platform with Bash (≥4.0), R (≥3.1.1) and BEDTools (≥2.19.0). It requires the following R packages: GenomicRanges, Sushi, data.table, preprocessCore and dynamicTreeCut.
Supplementary information
Supplementary data are available at Bioinformatics online.
•A method is presented for 5hmC detection and analysis using Infinium 450K BeadChips.•The oxBS-450K method can discriminate between 5mC and 5hmC in human gDNA•5hmC levels were quantified genome-wide ...in 3 distinct biological samples.•The reported 5hmC signal was validated using mass spectrometry and pyrosequencing.•The effects of differing amounts of input DNA on final 5hmC call rate are discussed.
DNA methylation analysis has become an integral part of biomedical research. For high-throughput applications such as epigenome-wide association studies, the Infinium HumanMethylation450 (450K) BeadChip is currently the platform of choice. However, BeadChip processing relies on traditional bisulfite (BS) based protocols which cannot discriminate between 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC). Here, we report the adaptation of the recently developed oxidative bisulfite (oxBS) chemistry to specifically detect both 5mC and 5hmC in a single workflow using 450K BeadChips, termed oxBS-450K. Supported by validation using mass spectrometry and pyrosequencing, we demonstrate reproducible (R2>0.99) detection of 5hmC in human brain tissue using the optimised oxBS-450K protocol described here.
High-grade gliomas defined by histone 3 K27M driver mutations exhibit global loss of H3K27 trimethylation and reciprocal gain of H3K27 acetylation, respectively shaping repressive and active ...chromatin landscapes. We generated tumor-derived isogenic models bearing this mutation and show that it leads to pervasive H3K27ac deposition across the genome. In turn, active enhancers and promoters are not created de novo and instead reflect the epigenomic landscape of the cell of origin. H3K27ac is enriched at repeat elements, resulting in their increased expression, which in turn can be further amplified by DNA demethylation and histone deacetylase inhibitors providing an exquisite therapeutic vulnerability. These agents may therefore modulate anti-tumor immune responses as a therapeutic modality for this untreatable disease.
Display omitted
•H3K27M mutant gliomas exhibit an enhancer landscape reflecting lineage of origin•Elevated H3K27 acetylation is pervasively distributed across the epigenome•Repeat element expression is de-repressed by H3K27M mutation•H3K27M cells are more vulnerable to DNA demethylation triggering viral mimicry
Krug et al. reveal increased global H3K27ac deposition across the genome without creation of de novo active enhancers or promoters in high-grade glioma (HGG) with H3K27M mutations. H3K27ac enrichment at repeat elements in H3K27M HGG increases their expression, conferring sensitivity to epigenetic therapies.