Simulating the behaviour of glioblastoma multiforme based on patient MRI during treatments.

Glioblastoma multiforme is a brain cancer that still shows poor prognosis for patients despite the active research for new treatments. In this work, the goal is to model and simulate the evolution of tumour associated angiogenesis and the therapeutic response to glioblastoma multiforme. Multiple phenomena are modelled in order to fit different biological pathways, such as the cellular cycle, apoptosis, hypoxia or angiogenesis. This leads to a nonlinear system with 4 equations and 4 unknowns: the density of tumour cells, the [Formula: see text] concentration, the density of endothelial cells and the vascular endothelial growth factor concentration. This system is solved numerically on a mesh fitting the geometry of the brain and the tumour of a patient based on a 2D slice of MRI. We show that our numerical scheme is positive, and we give the energy estimates on the discrete solution to ensure its existence. The numerical scheme uses nonlinear control volume finite elements in space and is implicit in time. Numerical simulations have been done using the different standard treatments: surgery, chemotherapy and radiotherapy, in order to conform to the behaviour of a tumour in response to treatments according to empirical clinical knowledge. We find that our theoretical model exhibits realistic behaviours.

Reduction of RUNX1 transcription factor activity by a CBFA2T3-mimicking peptide: application to B cell precursor acute lymphoblastic leukemia.

BACKGROUND: B Cell Precursor Acute Lymphoblastic Leukemia (BCP-ALL) is the most common pediatric cancer. Identifying key players involved in proliferation of BCP-ALL cells is crucial to propose new therapeutic targets. Runt Related Transcription Factor 1 (RUNX1) and Core-Binding Factor Runt Domain Alpha Subunit 2 Translocated To 3 (CBFA2T3, ETO2, MTG16) are master regulators of hematopoiesis and are implicated in leukemia. METHODS: We worked with BCP-ALL mononuclear bone marrow patients' cells and BCP-ALL cell lines, and performed Chromatin Immunoprecipitations followed by Sequencing (ChIP-Seq), co-immunoprecipitations (co-IP), proximity ligation assays (PLA), luciferase reporter assays and mouse xenograft models. RESULTS: We demonstrated that CBFA2T3 transcript levels correlate with RUNX1 expression in the pediatric t(12;21) ETV6-RUNX1 BCP-ALL. By ChIP-Seq in BCP-ALL patients' cells and cell lines, we found that RUNX1 is recruited on its promoter and on an enhancer of CBFA2T3 located - 2 kb upstream CBFA2T3 promoter and that, subsequently, the transcription factor RUNX1 drives both RUNX1 and CBFA2T3 expression. We demonstrated that, mechanistically, RUNX1 and CBFA2T3 can be part of the same complex allowing CBFA2T3 to strongly potentiate the activity of the transcription factor RUNX1. Finally, we characterized a CBFA2T3-mimicking peptide that inhibits the interaction between RUNX1 and CBFA2T3, abrogating the activity of this transcription complex and reducing BCP-ALL lymphoblast proliferation. CONCLUSIONS: Altogether, our findings reveal a novel and important activation loop between the transcription regulator CBFA2T3 and the transcription factor RUNX1 that promotes BCP-ALL proliferation, supporting the development of an innovative therapeutic approach based on the NHR2 subdomain of CBFA2T3 protein.

DNA hydroxymethylation is associated with disease severity and persists at enhancers of oncogenic regions in multiple myeloma.

BACKGROUND: Multiple myeloma (MM) is a heterogeneous plasma cell malignancy that remains challenging to cure. Global hypomethylation correlates with an aggressive phenotype of the disease, while hypermethylation is observed at particular regions of myeloma such as B cell-specific enhancers. The recently discovered active epigenetic mark 5-hydroxymethylCytosine (5hmC) may also play a role in tumor biology; however, little is known about its level and distribution in myeloma. In this study, we investigated the global level and the genomic localization of 5hmC in myeloma cells from 40 newly diagnosed patients, including paired relapses, and of control individuals. RESULTS: Compared to normal plasma cells, we found global 5hmC levels to be lower in myeloma (P < 0.001). Higher levels of 5hmC were found in lower grades of the International Staging System prognostic index (P < 0.05) and tend to associate with a longer overall survival (P < 0.1). From the hydroxymethylome data, we observed that the remaining 5hmC is organized in large domains overlapping with active chromatin marks and chromatin opening. We discovered that 5hmC strongly persists at key oncogenic genes such as CCND1, CCND2 and MMSET and characterized domains that are specifically hydroxymethylated in myeloma subgroups. Novel 5hmC-enriched domains were found at putative enhancers of CCND2 and MYC in newly diagnosed patients. CONCLUSIONS: 5hmC level is associated with clinical aspects of MM. Mapping 5hmC at a genome-wide level provides insights into the disease biology directly from genomic DNA, which makes it a potent mark to study epigenetics on large patient cohorts.

N6-Adenosine Methylation of miRNA-200b-3p Influences Its Functionality and Is a Theranostic Tool.

MicroRNAs (miRNAs or miRs) play crucial roles in biological and pathological processes. Some miRNAs also appear as promising biomarkers and therapeutic tools. However, the epitranscriptomic regulation of miRNAs is not yet fully elucidated in all of their fields of application. We report that adenosine methylation of miR-200b-3p inhibits its repressive function toward its mRNA targets such as XIAP by blocking the formation of the miRNA/3' UTRmRNA duplex. Our data indicate that the adenosine methylation of miR-200b-3p is associated with the survival of glioblastoma patients. Collectively, our data support the idea that the adenosine methylation of miR-200b-3p can be used as a prodrug having a selective cytotoxicity against cancer cells (while being harmless to peripheral blood mononuclear cells [PBMCs], astrocytes, neurons, and hepatocytes).

Proteogenomic analysis of Inhibitor of Differentiation 4 (ID4) in basal-like breast cancer.

BACKGROUND: Basal-like breast cancer (BLBC) is a poorly characterised, heterogeneous disease. Patients are diagnosed with aggressive, high-grade tumours and often relapse with chemotherapy resistance. Detailed understanding of the molecular underpinnings of this disease is essential to the development of personalised therapeutic strategies. Inhibitor of differentiation 4 (ID4) is a helix-loop-helix transcriptional regulator required for mammary gland development. ID4 is overexpressed in a subset of BLBC patients, associating with a stem-like poor prognosis phenotype, and is necessary for the growth of cell line models of BLBC through unknown mechanisms. METHODS: Here, we have defined unique molecular insights into the function of ID4 in BLBC and the related disease high-grade serous ovarian cancer (HGSOC), by combining RIME proteomic analysis, ChIP-seq mapping of genomic binding sites and RNA-seq. RESULTS: These studies reveal novel interactions with DNA damage response proteins, in particular, mediator of DNA damage checkpoint protein 1 (MDC1). Through MDC1, ID4 interacts with other DNA repair proteins (γH2AX and BRCA1) at fragile chromatin sites. ID4 does not affect transcription at these sites, instead binding to chromatin following DNA damage. Analysis of clinical samples demonstrates that ID4 is amplified and overexpressed at a higher frequency in BRCA1-mutant BLBC compared with sporadic BLBC, providing genetic evidence for an interaction between ID4 and DNA damage repair deficiency. CONCLUSIONS: These data link the interactions of ID4 with MDC1 to DNA damage repair in the aetiology of BLBC and HGSOC.

The hydroxymethylome of multiple myeloma identifies FAM72D as a 1q21 marker linked to proliferation.

Cell identity relies on the cross-talk between genetics and epigenetics and their impact on gene expression. Oxidation of 5-methylcytosine (5mC) into 5-hydroxymethylcytosine (5hmC) is the first step of an active DNA demethylation process occurring mainly at enhancers and gene bodies and, as such, participates in processes governing cell identity in normal and pathological conditions. Although genetic alterations are well documented in multiple myeloma (MM), epigenetic alterations associated with this disease have not yet been thoroughly analyzed. To gain insight into the biology of MM, genome-wide 5hmC profiles were obtained and showed that regions enriched in this modified base overlap with MM enhancers and super enhancers and are close to highly expressed genes. Through the definition of a MM-specific 5hmC signature, we identified FAM72D as a poor prognostic gene located on 1q21, a region amplified in high risk myeloma. We further uncovered that FAM72D functions as part of the FOXM1 transcription factor network controlling cell proliferation and survival and we evidenced an increased sensitivity of cells expressing high levels of FOXM1 and FAM72 to epigenetic drugs targeting histone deacetylases and DNA methyltransferases.

Interplay between transcription regulators RUNX1 and FUBP1 activates an enhancer of the oncogene c-KIT and amplifies cell proliferation.

Runt-related transcription factor 1 (RUNX1) is a well-known master regulator of hematopoietic lineages but its mechanisms of action are still not fully understood. Here, we found that RUNX1 localizes on active chromatin together with Far Upstream Binding Protein 1 (FUBP1) in human B-cell precursor lymphoblasts, and that both factors interact in the same transcriptional regulatory complex. RUNX1 and FUBP1 chromatin localization identified c-KIT as a common target gene. We characterized two regulatory regions, at +700 bp and +30 kb within the first intron of c-KIT, bound by both RUNX1 and FUBP1, and that present active histone marks. Based on these regions, we proposed a novel FUBP1 FUSE-like DNA-binding sequence on the +30 kb enhancer. We demonstrated that FUBP1 and RUNX1 cooperate for the regulation of the expression of the oncogene c-KIT. Notably, upregulation of c-KIT expression by FUBP1 and RUNX1 promotes cell proliferation and renders cells more resistant to the c-KIT inhibitor imatinib mesylate, a common therapeutic drug. These results reveal a new mechanism of action of RUNX1 that implicates FUBP1, as a facilitator, to trigger transcriptional regulation of c-KIT and to regulate cell proliferation. Deregulation of this regulatory mechanism may explain some oncogenic function of RUNX1 and FUBP1.

Cytosine modifications modulate the chromatin architecture of transcriptional enhancers.

Epigenetic mechanisms are believed to play key roles in the establishment of cell-specific transcription programs. Accordingly, the modified bases 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) have been observed in DNA of genomic regulatory regions such as enhancers, and oxidation of 5mC into 5hmC by Ten-eleven translocation (TET) proteins correlates with enhancer activation. However, the functional relationship between cytosine modifications and the chromatin architecture of enhancers remains elusive. To gain insights into their function, 5mC and 5hmC levels were perturbed by inhibiting DNA methyltransferases and TETs during differentiation of mouse embryonal carcinoma cells into neural progenitors, and chromatin characteristics of enhancers bound by the pioneer transcription factors FOXA1, MEIS1, and PBX1 were interrogated. In a large fraction of the tested enhancers, inhibition of DNA methylation was associated with a significant increase in monomethylation of H3K4, a characteristic mark of enhancer priming. In addition, at some specific enhancers, 5mC oxidation by TETs facilitated chromatin opening, a process that may stabilize MEIS1 binding to these genomic regions.

ERRα induces H3K9 demethylation by LSD1 to promote cell invasion.

Lysine Specific Demethylase 1 (LSD1) removes mono- and dimethyl groups from lysine 4 of histone H3 (H3K4) or H3K9, resulting in repressive or activating (respectively) transcriptional histone marks. The mechanisms that control the balance between these two antagonist activities are not understood. We here show that LSD1 and the orphan nuclear receptor estrogen-related receptor α (ERRα) display commonly activated genes. Transcriptional activation by LSD1 and ERRα involves H3K9 demethylation at the transcriptional start site (TSS). Strikingly, ERRα is sufficient to induce LSD1 to demethylate H3K9 in vitro. The relevance of this mechanism is highlighted by functional data. LSD1 and ERRα coregulate several target genes involved in cell migration, including the MMP1 matrix metallo-protease, also activated through H3K9 demethylation at the TSS. Depletion of LSD1 or ERRα reduces the cellular capacity to invade the extracellular matrix, a phenomenon that is rescued by MMP1 reexpression. Altogether our results identify a regulatory network involving a direct switch in the biochemical activities of a histone demethylase, leading to increased cell invasion.

FOXA1 Directs H3K4 Monomethylation at Enhancers via Recruitment of the Methyltransferase MLL3.

FOXA1 is a pioneer factor that binds to enhancer regions that are enriched in H3K4 mono- and dimethylation (H3K4me1 and H3K4me2). We performed a FOXA1 rapid immunoprecipitation mass spectrometry of endogenous proteins (RIME) screen in ERα-positive MCF-7 breast cancer cells and found histone-lysine N-methyltransferase (MLL3) as the top FOXA1-interacting protein. MLL3 is typically thought to induce H3K4me3 at promoter regions, but recent findings suggest it may contribute to H3K4me1 deposition. We performed MLL3 chromatin immunoprecipitation sequencing (ChIP-seq) in breast cancer cells, and MLL3 was shown to occupy regions marked by FOXA1 occupancy and H3K4me1 and H3K4me2. MLL3 binding was dependent on FOXA1, indicating that FOXA1 recruits MLL3 to chromatin. MLL3 silencing decreased H3K4me1 at enhancer elements but had no appreciable impact on H3K4me3 at enhancer elements. We propose a mechanism whereby the pioneer factor FOXA1 recruits the chromatin modifier MLL3 to facilitate the deposition of H3K4me1 histone marks, subsequently demarcating active enhancer elements.

Progesterone receptor modulates ERα action in breast cancer.

Progesterone receptor (PR) expression is used as a biomarker of oestrogen receptor-α (ERα) function and breast cancer prognosis. Here we show that PR is not merely an ERα-induced gene target, but is also an ERα-associated protein that modulates its behaviour. In the presence of agonist ligands, PR associates with ERα to direct ERα chromatin binding events within breast cancer cells, resulting in a unique gene expression programme that is associated with good clinical outcome. Progesterone inhibited oestrogen-mediated growth of ERα(+) cell line xenografts and primary ERα(+) breast tumour explants, and had increased anti-proliferative effects when coupled with an ERα antagonist. Copy number loss of PGR, the gene coding for PR, is a common feature in ERα(+) breast cancers, explaining lower PR levels in a subset of cases. Our findings indicate that PR functions as a molecular rheostat to control ERα chromatin binding and transcriptional activity, which has important implications for prognosis and therapeutic interventions.

Dynamic estrogen receptor interactomes control estrogen-responsive trefoil Factor (TFF) locus cell-specific activities.

Estradiol signaling is ideally suited for analyzing the molecular and functional linkages between the different layers of information directing transcriptional regulations: the DNA sequence, chromatin modifications, and the spatial organization of the genome. Hence, the estrogen receptor (ER) can bind at a distance from its target genes and engages timely and spatially coordinated processes to regulate their expression. In the context of the coordinated regulation of colinear genes, identifying which ER binding sites (ERBSs) regulate a given gene still remains a challenge. Here, we investigated the coordination of such regulatory events at a 2-Mb genomic locus containing the estrogen-sensitive trefoil factor (TFF) cluster of genes in breast cancer cells. We demonstrate that this locus exhibits a hormone- and cohesin-dependent reduction in the plasticity of its three-dimensional organization that allows multiple ERBSs to be dynamically brought to the vicinity of estrogen-sensitive genes. Additionally, by using triplex-forming oligonucleotides, we could precisely document the functional links between ER engagement at given ERBSs and the regulation of particular genes. Hence, our data provide evidence of a formerly suggested cooperation of enhancers toward gene regulation and also show that redundancy between ERBSs can occur.

Development of an Illumina-based ChIP-exonuclease method provides insight into FoxA1-DNA binding properties.

ChIP-exonuclease (ChIP-exo) is a modified ChIP-seq approach for high resolution mapping of transcription factor DNA sites. We describe an Illumina-based ChIP-exo method which provides a global improvement of the data quality of estrogen receptor (ER) ChIP and insights into the motif structure for key ER-associated factors. ChIP-exo of the ER pioneer factor FoxA1 identifies protected DNA with a predictable 8 bp overhang from the Forkhead motif, which we term mesas. We show that mesas occur in multiple cellular contexts and exist as single or overlapping motifs. Our Illumina-based ChIP-exo provides high resolution mapping of transcription factor binding sites.

Endogenous purification reveals GREB1 as a key estrogen receptor regulatory factor.

Estrogen receptor-α (ER) is the driving transcription factor in most breast cancers, and its associated proteins can influence drug response, but direct methods for identifying interacting proteins have been limited. We purified endogenous ER using an approach termed RIME (rapid immunoprecipitation mass spectrometry of endogenous proteins) and discovered the interactome under agonist- and antagonist-liganded conditions in breast cancer cells, revealing transcriptional networks in breast cancer. The most estrogen-enriched ER interactor is GREB1, a potential clinical biomarker with no known function. GREB1 is shown to be a chromatin-bound ER coactivator and is essential for ER-mediated transcription, because it stabilizes interactions between ER and additional cofactors. We show a GREB1-ER interaction in three xenograft tumors, and using a directed protein-protein approach, we find GREB1-ER interactions in half of ER(+) primary breast cancers. This finding is supported by histological expression of GREB1, which shows that GREB1 is expressed in half of ER(+) cancers, and predicts good clinical outcome. These findings reveal an unexpected role for GREB1 as an estrogen-specific ER cofactor that is expressed in drug-sensitive contexts.

Dynamic hydroxymethylation of deoxyribonucleic acid marks differentiation-associated enhancers.

Enhancers are developmentally controlled transcriptional regulatory regions whose activities are modulated through histone modifications or histone variant deposition. In this study, we show by genome-wide mapping that the newly discovered deoxyribonucleic acid (DNA) modification 5-hydroxymethylcytosine (5hmC) is dynamically associated with transcription factor binding to distal regulatory sites during neural differentiation of mouse P19 cells and during adipocyte differentiation of mouse 3T3-L1 cells. Functional annotation reveals that regions gaining 5hmC are associated with genes expressed either in neural tissues when P19 cells undergo neural differentiation or in adipose tissue when 3T3-L1 cells undergo adipocyte differentiation. Furthermore, distal regions gaining 5hmC together with H3K4me2 and H3K27ac in P19 cells behave as differentiation-dependent transcriptional enhancers. Identified regions are enriched in motifs for transcription factors regulating specific cell fates such as Meis1 in P19 cells and PPARγ in 3T3-L1 cells. Accordingly, a fraction of hydroxymethylated Meis1 sites were associated with a dynamic engagement of the 5-methylcytosine hydroxylase Tet1. In addition, kinetic studies of cytosine hydroxymethylation of selected enhancers indicated that DNA hydroxymethylation is an early event of enhancer activation. Hence, acquisition of 5hmC in cell-specific distal regulatory regions may represent a major event of enhancer progression toward an active state and participate in selective activation of tissue-specific genes.

Epigenetic switch involved in activation of pioneer factor FOXA1-dependent enhancers.

Transcription factors (TFs) bind specifically to discrete regions of mammalian genomes called cis-regulatory elements. Among those are enhancers, which play key roles in regulation of gene expression during development and differentiation. Despite the recognized central regulatory role exerted by chromatin in control of TF functions, much remains to be learned regarding the chromatin structure of enhancers and how it is established. Here, we have analyzed on a genomic-scale enhancers that recruit FOXA1, a pioneer transcription factor that triggers transcriptional competency of these cis-regulatory sites. Importantly, we found that FOXA1 binds to genomic regions showing local DNA hypomethylation and that its cell-type-specific recruitment to chromatin is linked to differential DNA methylation levels of its binding sites. Using neural differentiation as a model, we showed that induction of FOXA1 expression and its subsequent recruitment to enhancers is associated with DNA demethylation. Concomitantly, histone H3 lysine 4 methylation is induced at these enhancers. These epigenetic changes may both stabilize FOXA1 binding and allow for subsequent recruitment of transcriptional regulatory effectors. Interestingly, when cloned into reporter constructs, FOXA1-dependent enhancers were able to recapitulate their cell type specificity. However, their activities were inhibited by DNA methylation. Hence, these enhancers are intrinsic cell-type-specific regulatory regions of which activities have to be potentiated by FOXA1 through induction of an epigenetic switch that includes notably DNA demethylation.