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.

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.

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.

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.

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.

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.

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.

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.

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.

Quantitative Analysis of Intra-chromosomal Contacts: The 3C-qPCR Method.

The chromosome conformation capture (3C) technique is fundamental to many population-based methods investigating chromatin dynamics and organization in eukaryotes. Here, we provide a modified quantitative 3C (3C-qPCR) protocol for improved quantitative analyses of intra-chromosomal contacts. We also describe an algorithm for data normalization which allows more accurate comparisons between contact profiles.

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.

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.

Contribution of Topological Domains and Loop Formation to 3D Chromatin Organization.

Recent investigations on 3D chromatin folding revealed that the eukaryote genomes are both highly compartmentalized and extremely dynamic. This review presents the most recent advances in topological domains' organization of the eukaryote genomes and discusses the relationship to chromatin loop formation. CTCF protein appears as a central factor of these two organization levels having either a strong insulating role at TAD borders, or a weaker architectural role in chromatin loop formation. TAD borders directly impact on chromatin dynamics by restricting contacts within specific genomic portions thus confining chromatin loop formation within TADs. We discuss how sub-TAD chromatin dynamics, constrained into a recently described statistical helix conformation, can produce functional interactions by contact stabilization.

Chromatin fiber allostery and the epigenetic code.

The notion of allostery introduced for proteins about fifty years ago has been extended since then to DNA allostery, where a locally triggered DNA structural transition remotely controls other DNA-binding events. We further extend this notion and propose that chromatin fiber allosteric transitions, induced by histone-tail covalent modifications, may play a key role in transcriptional regulation. We present an integrated scenario articulating allosteric mechanisms at different scales: allosteric transitions of the condensed chromatin fiber induced by histone-tail acetylation modify the mechanical constraints experienced by the embedded DNA, thus possibly controlling DNA-binding of allosteric transcription factors or further allosteric mechanisms at the linker DNA level. At a higher scale, different epigenetic constraints delineate different statistically dominant subsets of accessible chromatin fiber conformations, which each favors the assembly of dedicated regulatory complexes, as detailed on the emblematic example of the mouse Igf2-H19 gene locus and its parental imprinting. This physical view offers a mechanistic and spatially structured explanation of the observed correlation between transcriptional activity and histone modifications. The evolutionary origin of allosteric control supports to speak of an 'epigenetic code', by which events involved in transcriptional regulation are encoded in histone modifications in a context-dependent way.

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.

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.

H19 antisense RNA can up-regulate Igf2 transcription by activation of a novel promoter in mouse myoblasts.

It was recently shown that a long non-coding RNA (lncRNA), that we named the 91H RNA (i.e. antisense H19 transcript), is overexpressed in human breast tumours and contributes in trans to the expression of the Insulin-like Growth Factor 2 (IGF2) gene on the paternal chromosome. Our preliminary experiments suggested that an H19 antisense transcript having a similar function may also be conserved in the mouse. In the present work, we further characterise the mouse 91H RNA and, using a genetic complementation approach in H19 KO myoblast cells, we show that ectopic expression of the mouse 91H RNA can up-regulate Igf2 expression in trans despite almost complete unmethylation of the Imprinting-Control Region (ICR). We then demonstrate that this activation occurs at the transcriptional level by activation of a previously unknown Igf2 promoter which displays, in mouse tissues, a preferential mesodermic expression (Pm promoter). Finally, our experiments indicate that a large excess of the H19 transcript can counteract 91H-mediated Igf2 activation. Our work contributes, in conjunction with other recent findings, to open new horizons to our understanding of Igf2 gene regulation and functions of the 91H/H19 RNAs in normal and pathological conditions.

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.