DLK1/DIO3 locus upregulation by a ß-catenin-dependent enhancer drives cell proliferation and liver tumorigenesis.

The CTNNB1 gene, encoding ß-catenin, is frequently mutated in hepatocellular carcinoma (HCC, ∼30%) and in hepatoblastoma (HB, >80%), in which DLK1/DIO3 locus induction is correlated with CTNNB1 mutations. Here, we aim to decipher how sustained ß-catenin activation regulates DLK1/DIO3 locus expression and the role this locus plays in HB and HCC development in mouse models deleted for Apc (ApcΔhep) or Ctnnb1-exon 3 (ß-cateninΔExon3) and in human CTNNB1-mutated hepatic cancer cells. We identified an enhancer site bound by TCF-4/ß-catenin complexes in an open conformation upon sustained ß-catenin activation (DLK1-Wnt responsive element [WRE]) and increasing DLK1/DIO3 locus transcription in ß-catenin-mutated human HB and mouse models. DLK1-WRE editing by CRISPR-Cas9 approach impaired DLK1/DIO3 locus expression and slowed tumor growth in subcutaneous CTNNB1-mutated tumor cell grafts, ApcΔhep HB and ß-cateninΔExon3 HCC. Tumor growth inhibition resulted either from increased FADD expression and subsequent caspase-3 cleavage in the first case or from decreased expression of cell cycle actors regulated by FoxM1 in the others. Therefore, the DLK1/DIO3 locus is an essential determinant of FoxM1-dependent cell proliferation during ß-catenin-driven liver tumorigenesis. Targeting the DLK1-WRE enhancer to silence the DLK1/DIO3 locus might thus represent an interesting therapeutic strategy to restrict tumor growth in primary liver cancers with CTNNB1 mutations.

DeSUMOylation of chromatin-bound proteins limits the rapid transcriptional reprogramming induced by daunorubicin in acute myeloid leukemias.

Genotoxicants have been used for decades as front-line therapies against cancer on the basis of their DNA-damaging actions. However, some of their non-DNA-damaging effects are also instrumental for killing dividing cells. We report here that the anthracycline Daunorubicin (DNR), one of the main drugs used to treat Acute Myeloid Leukemia (AML), induces rapid (3 h) and broad transcriptional changes in AML cells. The regulated genes are particularly enriched in genes controlling cell proliferation and death, as well as inflammation and immunity. These transcriptional changes are preceded by DNR-dependent deSUMOylation of chromatin proteins, in particular at active promoters and enhancers. Surprisingly, inhibition of SUMOylation with ML-792 (SUMO E1 inhibitor), dampens DNR-induced transcriptional reprogramming. Quantitative proteomics shows that the proteins deSUMOylated in response to DNR are mostly transcription factors, transcriptional co-regulators and chromatin organizers. Among them, the CCCTC-binding factor CTCF is highly enriched at SUMO-binding sites found in cis-regulatory regions. This is notably the case at the promoter of the DNR-induced NFKB2 gene. DNR leads to a reconfiguration of chromatin loops engaging CTCF- and SUMO-bound NFKB2 promoter with a distal cis-regulatory region and inhibition of SUMOylation with ML-792 prevents these changes.

Contribution of 3D genome topological domains to genetic risk of cancers: a genome-wide computational study.

BACKGROUND: Genome-wide association studies have identified statistical associations between various diseases, including cancers, and a large number of single-nucleotide polymorphisms (SNPs). However, they provide no direct explanation of the mechanisms underlying the association. Based on the recent discovery that changes in three-dimensional genome organization may have functional consequences on gene regulation favoring diseases, we investigated systematically the genome-wide distribution of disease-associated SNPs with respect to a specific feature of 3D genome organization: topologically associating domains (TADs) and their borders. RESULTS: For each of 449 diseases, we tested whether the associated SNPs are present in TAD borders more often than observed by chance, where chance (i.e., the null model in statistical terms) corresponds to the same number of pointwise loci drawn at random either in the entire genome, or in the entire set of disease-associated SNPs listed in the GWAS catalog. Our analysis shows that a fraction of diseases displays such a preferential localization of their risk loci. Moreover, cancers are relatively more frequent among these diseases, and this predominance is generally enhanced when considering only intergenic SNPs. The structure of SNP-based diseasome networks confirms that localization of risk loci in TAD borders differs between cancers and non-cancer diseases. Furthermore, different TAD border enrichments are observed in embryonic stem cells and differentiated cells, consistent with changes in topological domains along embryogenesis and delineating their contribution to disease risk. CONCLUSIONS: Our results suggest that, for certain diseases, part of the genetic risk lies in a local genetic variation affecting the genome partitioning in topologically insulated domains. Investigating this possible contribution to genetic risk is particularly relevant in cancers. This study thus opens a way of interpreting genome-wide association studies, by distinguishing two types of disease-associated SNPs: one with an effect on an individual gene, the other acting in interplay with 3D genome organization.

Cocaine-related DNA methylation in caudate neurons alters 3D chromatin structure of the IRXA gene cluster.

Epigenetic mechanisms, like those involving DNA methylation, are thought to mediate the relationship between chronic cocaine dependence and molecular changes in addiction-related neurocircuitry, but have been understudied in human brain. We initially used reduced representation bisulfite sequencing (RRBS) to generate a methylome-wide profile of cocaine dependence in human post-mortem caudate tissue. We focused on the Iroquois Homeobox A (IRXA) gene cluster, where hypomethylation in exon 3 of IRX2 in neuronal nuclei was associated with cocaine dependence. We replicated this finding in an independent cohort and found similar results in the dorsal striatum from cocaine self-administering mice. Using epigenome editing and 3C assays, we demonstrated a causal relationship between methylation within the IRX2 gene body, CTCF protein binding, three-dimensional (3D) chromatin interaction, and gene expression. Together, these findings suggest that cocaine-related hypomethylation of IRX2 contributes to the development and maintenance of cocaine dependence through alterations in 3D chromatin structure in the caudate nucleus.

High-salt-recovered sequences are associated with the active chromosomal compartment and with large ribonucleoprotein complexes including nuclear bodies.

The mammalian cell nucleus contains numerous discrete suborganelles named nuclear bodies. While recruitment of specific genomic regions into these large ribonucleoprotein (RNP) complexes critically contributes to higher-order functional chromatin organization, such regions remain ill-defined. We have developed the high-salt-recovered sequences-sequencing (HRS-seq) method, a straightforward genome-wide approach whereby we isolated and sequenced genomic regions associated with large high-salt insoluble RNP complexes. By using mouse embryonic stem cells (ESCs), we showed that these regions essentially correspond to the most highly expressed genes, and to cis-regulatory sequences like super-enhancers, that belong to the active A chromosomal compartment. They include both cell-type-specific genes, such as pluripotency genes in ESCs, and housekeeping genes associated with nuclear bodies, such as histone and snRNA genes that are central components of Histone Locus Bodies and Cajal bodies. We conclude that HRSs are associated with the active chromosomal compartment and with large RNP complexes including nuclear bodies. Association of such chromosomal regions with nuclear bodies is in agreement with the recently proposed phase separation model for transcription control and might thus play a central role in organizing the active chromosomal compartment in mammals.

Hypoxia differently modulates the release of mitochondrial and nuclear DNA.

BACKGROUND: We investigated the influence of hypoxia on the concentration of mitochondrial and nuclear cell-free DNA (McfDNA and NcfDNA, respectively). METHOD: By an ultra-sensitive quantitative PCR-based assay, McfDNA and NcfDNA were measured in the supernatants of different colorectal cell lines, and in the plasma of C57/Bl6 mice engrafted with TC1 tumour cells, in normoxic or hypoxic conditions. RESULTS: Our data when setting cell culture conditions highlighted the higher stability of McfDNA as compared to NcfDNA and revealed that cancer cells released amounts of nuclear DNA equivalent to the mass of a chromosome over a 6-h duration of incubation. In cell model, hypoxia induced a great increase in NcfDNA and McfDNA concentrations within the first 24 h. After this period, cfDNA total concentrations remained stable in hypoxia consecutive to a decrease of nuclear DNA release, and noteworthy, to a complete inhibition of daily mitochondrial DNA release. In TC1-engrafted mice submitted to intermittent hypoxia, plasma NcfDNA levels are much higher than in mice bred in normoxia, unlike plasma McfDNA concentration that is not impacted by hypoxia. CONCLUSION: This study suggests that hypoxia negatively modulates nuclear and, particularly, mitochondrial DNA releases in long-term hypoxia, and revealed that the underlying mechanisms are differently regulated.

EHMT2 directs DNA methylation for efficient gene silencing in mouse embryos.

The extent to which histone modifying enzymes contribute to DNA methylation in mammals remains unclear. Previous studies suggested a link between the lysine methyltransferase EHMT2 (also known as G9A and KMT1C) and DNA methylation in the mouse. Here, we used a model of knockout mice to explore the role of EHMT2 in DNA methylation during mouse embryogenesis. The Ehmt2 gene is expressed in epiblast cells but is dispensable for global DNA methylation in embryogenesis. In contrast, EHMT2 regulates DNA methylation at specific sequences that include CpG-rich promoters of germline-specific genes. These loci are bound by EHMT2 in embryonic cells, are marked by H3K9 dimethylation, and have strongly reduced DNA methylation in Ehmt2(-/-) embryos. EHMT2 also plays a role in the maintenance of germline-derived DNA methylation at one imprinted locus, the Slc38a4 gene. Finally, we show that DNA methylation is instrumental for EHMT2-mediated gene silencing in embryogenesis. Our findings identify EHMT2 as a critical factor that facilitates repressive DNA methylation at specific genomic loci during mammalian development.

Myod and H19-Igf2 locus interactions are required for diaphragm formation in the mouse.

The myogenic regulatory factor Myod and insulin-like growth factor 2 (Igf2) have been shown to interact in vitro during myogenic differentiation. In order to understand how they interact in vivo, we produced double-mutant mice lacking both the Myod and Igf2 genes. Surprisingly, these mice display neonatal lethality due to severe diaphragm atrophy. Alteration of diaphragm muscle development occurs as early as 15.5 days post-coitum in the double-mutant embryos and leads to a defect in the terminal differentiation of muscle progenitor cells. A negative-feedback loop was detected between Myod and Igf2 in embryonic muscles. Igf2 belongs to the imprinted H19-Igf2 locus. Molecular analyses show binding of Myod on a mesodermal enhancer (CS9) of the H19 gene. Chromatin conformation capture experiments reveal direct interaction of CS9 with the H19 promoter, leading to increased H19 expression in the presence of Myod. In turn, the non-coding H19 RNA represses Igf2 expression in trans. In addition, Igf2 also negatively regulates Myod expression, possibly by reducing the expression of the Srf transcription factor, a known Myod activator. In conclusion, Igf2 and Myod are tightly co-regulated in skeletal muscles and act in parallel pathways in the diaphragm, where they affect the progression of myogenic differentiation. Igf2 is therefore an essential player in the formation of a functional diaphragm in the absence of Myod.

Distinct polymer physics principles govern chromatin dynamics in mouse and Drosophila topological domains.

BACKGROUND: In higher eukaryotes, the genome is partitioned into large "Topologically Associating Domains" (TADs) in which the chromatin displays favoured long-range contacts. While a crumpled/fractal globule organization has received experimental supports at higher-order levels, the organization principles that govern chromatin dynamics within these TADs remain unclear. Using simple polymer models, we previously showed that, in mouse liver cells, gene-rich domains tend to adopt a statistical helix shape when no significant locus-specific interaction takes place. RESULTS: Here, we use data from diverse 3C-derived methods to explore chromatin dynamics within mouse and Drosophila TADs. In mouse Embryonic Stem Cells (mESC), that possess large TADs (median size of 840 kb), we show that the statistical helix model, but not globule models, is relevant not only in gene-rich TADs, but also in gene-poor and gene-desert TADs. Interestingly, this statistical helix organization is considerably relaxed in mESC compared to liver cells, indicating that the impact of the constraints responsible for this organization is weaker in pluripotent cells. Finally, depletion of histone H1 in mESC alters local chromatin flexibility but not the statistical helix organization. In Drosophila, which possesses TADs of smaller sizes (median size of 70 kb), we show that, while chromatin compaction and flexibility are finely tuned according to the epigenetic landscape, chromatin dynamics within TADs is generally compatible with an unconstrained polymer configuration. CONCLUSIONS: Models issued from polymer physics can accurately describe the organization principles governing chromatin dynamics in both mouse and Drosophila TADs. However, constraints applied on this dynamics within mammalian TADs have a peculiar impact resulting in a statistical helix organization.

Global profiling of DNA methylation erasure in mouse primordial germ cells.

Epigenetic reprogramming, characterized by loss of cytosine methylation and histone modifications, occurs during mammalian development in primordial germ cells (PGCs), yet the targets and kinetics of this process are poorly characterized. Here we provide a map of cytosine methylation on a large portion of the genome in developing male and female PGCs isolated from mouse embryos. We show that DNA methylation erasure is global and affects genes of various biological functions. We also reveal complex kinetics of demethylation that are initiated at most genes in early PGC precursors around embryonic day 8.0-9.0. In addition, besides intracisternal A-particles (IAPs), we identify rare LTR-ERV1 retroelements and single-copy sequences that resist global methylation erasure in PGCs as well as in preimplantation embryos. Our data provide important insights into the targets and dynamics of DNA methylation reprogramming in mammalian germ cells.

Chromatin loop organization of the junb locus in mouse dendritic cells.

The junb gene behaves as an immediate early gene in bacterial lipopolysaccharide (LPS)-stimulated dendritic cells (DCs), where its transient transcriptional activation is necessary for the induction of inflammatory cytokines. junb is a short gene and its transcriptional activation by LPS depends on the binding of NF-κB to an enhancer located just downstream of its 3' UTR. Here, we have addressed the mechanisms underlying the transcriptional hyper-reactivity of junb. Using transfection and pharmacological assays to complement chromatin immunoprecipitation analyses addressing the localization of histones, polymerase II, negative elongation factor (NELF)-, DRB sensitivity-inducing factor (DSIF)- and Positive Transcription Factor b complexes, we demonstrate that junb is a RNA Pol II-paused gene where Pol II is loaded in the transcription start site domain but poorly active. Moreover, High salt-Recovered Sequence, chromosome conformation capture (3C)- and gene transfer experiments show that (i) junb is organized in a nuclear chromatin loop bringing into close spatial proximity the upstream promoter region and the downstream enhancer and (ii) this configuration permits immediate Pol II release on the junb body on binding of LPS-activated NF-κB to the enhancer. Thus, our work unveils a novel topological framework underlying fast junb transcriptional response in DCs. Moreover, it also points to a novel layer of complexity in the modes of action of NF-κB.

Chromatin and DNA Dynamics in Mouse Models of Liver Cancers.

In recent years, important efforts have been made to understand how the expression of a specific gene repertoire correlates with chromatin accessibility, histone mark deposition, as well as with chromatin looping establishing connectivity with regulatory regions. The emergence of new techniques for genome-wide analyses and their progressive optimization to work on low amounts of material allows the scientific community to obtain an integrated view of transcriptional landscapes in physiology and disease. Here, we describe our own experience aiming at correlating the TCF-4/ß-catenin cistrome during liver tumorigenesis with chromatin remodeling, histone mark modifications, and long-distance DNA looping.

Long-range chromatin interactions at the mouse Igf2/H19 locus reveal a novel paternally expressed long non-coding RNA.

Parental genomic imprinting at the Igf2/H19 locus is controlled by a methylation-sensitive CTCF insulator that prevents the access of downstream enhancers to the Igf2 gene on the maternal chromosome. However, on the paternal chromosome, it remains unclear whether long-range interactions with the enhancers are restricted to the Igf2 promoters or whether they encompass the entire gene body. Here, using the quantitative chromosome conformation capture assay, we show that, in the mouse liver, the endodermal enhancers have low contact frequencies with the Igf2 promoters but display, on the paternal chromosome, strong interactions with the intragenic differentially methylated regions 1 and 2. Interestingly, we found that enhancers also interact with a so-far poorly characterized intergenic region of the locus that produces a novel imprinted long non-coding transcript that we named the paternally expressed Igf2/H19 intergenic transcript (PIHit) RNA. PIHit is expressed exclusively from the paternal chromosome, contains a novel discrete differentially methylated region in a highly conserved sequence and, surprisingly, does not require an intact ICR/H19 gene region for its imprinting. Altogether, our data reveal a novel imprinted domain in the Igf2/H19 locus and lead us to propose a model for chromatin folding of this locus on the paternal chromosome.

H19 acts as a trans regulator of the imprinted gene network controlling growth in mice.

The imprinted H19 gene produces a non-coding RNA of unknown function. Mice lacking H19 show an overgrowth phenotype, due to a cis effect of the H19 locus on the adjacent Igf2 gene. To explore the function of the RNA itself, we produced transgenic mice overexpressing H19. We observed postnatal growth reduction in two independent transgenic lines and detected a decrease of Igf2 expression in embryos. An extensive analysis of several other genes from the newly described imprinted gene network (IGN) was performed in both loss- and gain-of-function animals. We found that H19 deletion leads to the upregulation of several genes of the IGN. This overexpression is restored to the wild-type level by transgenic expression of H19. We therefore propose that the H19 gene participates as a trans regulator in the fine-tuning of this IGN in the mouse embryo. This is the first in vivo evidence of a functional role for the H19 RNA. Our results also bring further experimental evidence for the existence of the IGN and open new perspectives in the comprehension of the role of genomic imprinting in embryonic growth and in human imprinting pathologies.

Quantitative Chromosome Conformation Capture (3C-qPCR).

Many population-based methods investigating chromatin dynamics and organization in eukaryotes are based on the chromosome conformation capture (3C) method. Here, we provide an updated version of the quantitative 3C (3C-qPCR) protocol for improved and simplified quantitative analyses of intra-chromosomal contacts.

The High-Salt Recovered Sequence-Sequencing (HRS-seq) Method: Exploring Genome Association with Nuclear Bodies.

Recent works indicate that, at specific loci, interactions of chromatin with membrane-less organelles self-assembled through mechanisms of phase separation, like nuclear bodies, are crucial to regulate genome functions, and in particular transcription. Here we describe the protocol of the high-salt recovered sequence sequencing method whose principle relies on high-throughput sequencing of genomic DNA trapped into large RNP complexes that are made insoluble by high-salt treatments.

Endogenous Retroviral Sequences Behave as Putative Enhancers Controlling Gene Expression through HP1-Regulated Long-Range Chromatin Interactions.

About half of the mammalian genome is constituted of repeated elements, among which endogenous retroviruses (ERVs) are known to influence gene expression and cancer development. The HP1 (Heterochromatin Protein 1) proteins are known to be essential for heterochromatin establishment and function and its loss in hepatocytes leads to the reactivation of specific ERVs and to liver tumorigenesis. Here, by studying two ERVs located upstream of genes upregulated upon loss of HP1, Mbd1 and Trim24, we show that these HP1-dependent ERVs behave as either alternative promoters or as putative enhancers forming a loop with promoters of endogenous genes depending on the genomic context and HP1 expression level. These ERVs are characterised by a specific HP1-independent enrichment in heterochromatin-associated marks H3K9me3 and H4K20me3 as well as in the enhancer-specific mark H3K4me1, a combination that might represent a bookmark of putative ERV-derived enhancers. These ERVs are further enriched in a HP1-dependent manner in H3K27me3, suggesting a critical role of this mark together with HP1 in the silencing of the ERVs, as well as for the repression of the associated genes. Altogether, these results lead to the identification of a new regulatory hub involving the HP1-dependent formation of a physical loop between specific ERVs and endogenous genes.

CTCF is a DNA methylation-sensitive positive regulator of the INK/ARF locus.

The INK4B-ARF-INK4A (INK/ARF) locus is composed of three tumor suppressor genes, which are kept silenced by DNA methylation in different cancer types. In addition, a non-coding RNA (ANRIL) is transcribed in the anti-sense orientation upstream of the ARF gene. The resulting divergent promoter region is bound by the chromatin insulator protein CTCF in association with histone H3 tri-methylated on lysine 4, irrespective of transcription of ANRIL and ARF. Methylation of the overlapping CpG island abolishes CTCF binding and the associated modification, which can be restored by 5-Aza-2'-deoxycytidine (5-Aza-dC) treatment. shRNA knock down of CTCF expression dramatically reduces the induction of ANRIL and ARF, but also that of INK4A and INK4B expression by 5-Aza-dC. We propose that CTCF is an essential factor for transcription of the INK/ARF locus and that abrogation of its binding by DNA methylation contributes to the permanent silencing of several genes of the locus in tumors.

Exploring Mammalian Genome within Phase-Separated Nuclear Bodies: Experimental Methods and Implications for Gene Expression.

The importance of genome organization at the supranucleosomal scale in the control of gene expression is increasingly recognized today. In mammals, Topologically Associating Domains (TADs) and the active/inactive chromosomal compartments are two of the main nuclear structures that contribute to this organization level. However, recent works reviewed here indicate that, at specific loci, chromatin interactions with nuclear bodies could also be crucial to regulate genome functions, in particular transcription. They moreover suggest that these nuclear bodies are membrane-less organelles dynamically self-assembled and disassembled through mechanisms of phase separation. We have recently developed a novel genome-wide experimental method, High-salt Recovered Sequences sequencing (HRS-seq), which allows the identification of chromatin regions associated with large ribonucleoprotein (RNP) complexes and nuclear bodies. We argue that the physical nature of such RNP complexes and nuclear bodies appears to be central in their ability to promote efficient interactions between distant genomic regions. The development of novel experimental approaches, including our HRS-seq method, is opening new avenues to understand how self-assembly of phase-separated nuclear bodies possibly contributes to mammalian genome organization and gene expression.

Genomic imprinting controls matrix attachment regions in the Igf2 gene.

Genomic imprinting at the Igf2/H19 locus originates from allele-specific DNA methylation, which modifies the affinity of some proteins for their target sequences. Here, we show that AT-rich DNA sequences located in the vicinity of previously characterized differentially methylated regions (DMRs) of the imprinted Igf2 gene are conserved between mouse and human. These sequences have all the characteristics of matrix attachment regions (MARs), which are known as versatile regulatory elements involved in chromatin structure and gene expression. Combining allele-specific nuclear matrix binding assays and real-time PCR quantification, we show that retention of two of these Igf2 MARs (MAR0 and MAR2) in the nuclear matrix fraction depends on the tissue and is specific to the paternal allele. Furthermore, on this allele, the Igf2 MAR2 is functionally linked to the neighboring DMR2 while, on the maternal allele, it is controlled by the imprinting-control region. Our work clearly demonstrates that genomic imprinting controls matrix attachment regions in the Igf2 gene.