Enrique Vidal

Analysis of the methylome of tumor cell-free deoxyribonucleic acid (DNA; cfDNA) has emerged as a powerful non-invasive technique for cancer subtyping and prognosis. However, its application is frequently hampered by the quality and total cfDNA yield. Here, we demonstrate the feasibility of very low-input cfDNA for whole-methylome and copy-number profiling studies using enzymatic conversion of unmethylated cysteines [enzymatic methyl-seq (EM-seq)] to better preserve DNA integrity. We created a model for predicting genomic subtyping and prognosis with high accuracy. We validated our tool by comparing whole-genome CpG sequencing with in situ cohorts generated with bisulfite conversion and array hybridization, demonstrating that, despite the different techniques and sample origins, information on cfDNA methylation is comparable with in situ cohorts. Our findings support use of liquid biopsy followed by EM-seq to assess methylome of cancer patients, enabling validation in external cohorts. This advance is particularly relevant for rare cancers like neuroblastomas where liquid-biopsy volume is restricted by ethical regulations in pediatric patients.

Genome organization plays a pivotal role in transcription, but how transcription factors (TFs) rewire the structure of the genome to initiate and maintain the programs that lead to oncogenic transformation remains poorly understood. Acute promyelocytic leukemia (APL) is a fatal subtype of leukemia driven by a chromosomal translocation between the promyelocytic leukemia (PML) and retinoic acid receptor α (RARα) genes. We used primary hematopoietic stem and progenitor cells (HSPCs) and leukemic blasts that express the fusion protein PML-RARα as a paradigm to temporally dissect the dynamic changes in the epigenome, transcriptome, and genome architecture induced during oncogenic transformation. We found that PML-RARα initiates a continuum of topologic alterations, including switches from A to B compartments, transcriptional repression, loss of active histone marks, and gain of repressive histone marks. Our multiomics-integrated analysis identifies Klf4 as an early down-regulated gene in PML-RARα-driven leukemogenesis. Furthermore, we characterized the dynamic alterations in the Klf4 cis-regulatory network during APL progression and demonstrated that ectopic Klf4 overexpression can suppress self-renewal and reverse the differentiation block induced by PML-RARα. Our study provides a comprehensive in vivo temporal dissection of the epigenomic and topological reprogramming induced by an oncogenic TF and illustrates how topological architecture can be used to identify new drivers of malignant transformation.

Here, we report that in T47D breast cancer cells 50 pM progestin is sufficient to activate cell cycle entry and the progesterone gene expression program. At this concentration, equivalent to the progesterone blood levels found around the menopause, progesterone receptor (PR) binds only to 2800 genomic sites, which are accessible to ATAC cleavage prior to hormone exposure. These highly accessible sites (HAs) are surrounded by well-organized nucleosomes and exhibit breast enhancer features, including estrogen receptor alpha (ERα), higher FOXA1 and BRD4 (bromodomain containing 4) occupancy. Although HAs are enriched in RAD21 and CTCF, PR binding is the driving force for the most robust interactions with hormone-regulated genes. HAs show higher frequency of 3D contacts among themselves than with other PR binding sites, indicating colocalization in similar compartments. Gene regulation via HAs is independent of classical coregulators and ATP-activated remodelers, relying mainly on MAP kinase activation that enables PR nuclear engagement. HAs are also preferentially occupied by PR and ERα in breast cancer xenografts derived from MCF-7 cells as well as from patients, indicating their potential usefulness as targets for therapeutic intervention.

The 11 zinc finger (ZF) protein CTCF regulates topologically associating domain formation and transcription through selective binding to thousands of genomic sites. Here, we replaced endogenous CTCF in mouse embryonic stem cells with green-fluorescent-protein-tagged wild-type or mutant proteins lacking individual ZFs to identify additional determinants of CTCF positioning and function. While ZF1 and ZF8-ZF11 are not essential for cell survival, ZF8 deletion strikingly increases the DNA binding off-rate of mutant CTCF, resulting in reduced CTCF chromatin residence time. Loss of ZF8 results in widespread weakening of topologically associating domains, aberrant gene expression and increased genome-wide DNA methylation. Thus, important chromatin-templated processes rely on accurate CTCF chromatin residence time, which we propose depends on local sequence and chromatin context as well as global CTCF protein concentration.

Summary: A hallmark of chromosome organization is the partition into transcriptionally active A and repressed B compartments and into topologically associating domains (TADs). Both structures were regarded absent from the inactive X chromosome, but to be re-established with transcriptional reactivation and chromatin opening during X-reactivation. Here, we combine a tailor-made mouse iPSC-reprogramming system and high-resolution Hi-C to produce the first time-course combining gene reactivation, chromatin opening and chromosome topology during X-reactivation. Contrary to previous observations, we uncover A/B-like compartments on the inactive X harboring multiple subcompartments. While partial X-reactivation initiates within a compartment rich in X-inactivation escapees, it then occurs rapidly along the chromosome, coinciding with acquisition of naive pluripotency, leading to downregulation of Xist . Importantly, we find that TAD formation precedes transcription, suggesting them to be causally independent. Instead, TADs form first within Xist-poor compartments, establishing Xist as common denominator, opposing both gene reactivation and TAD formation through separate mechanisms. Graphical Summary:

Three-dimensional organization of the genome is important for transcriptional regulation^{1,2,3,4,5,6,7}. In mammals, CTCF and the cohesin complex create submegabase structures with elevated internal chromatin contact frequencies, called topologically associating domains (TADs)^8,9,10,11,12. Although TADs can contribute to transcriptional regulation, ablation of TAD organization by disrupting CTCF or the cohesin complex causes modest gene expression changes^13,14,15,16. In contrast, CTCF is required for cell cycle regulation¹⁷, embryonic development and formation of various adult cell types¹⁸. To uncouple the role of CTCF in cell-state transitions and cell proliferation, we studied the effect of CTCF depletion during the conversion of human leukemic B cells into macrophages with minimal cell division. CTCF depletion disrupts TAD organization but not cell transdifferentiation. In contrast, CTCF depletion in induced macrophages impairs the full-blown upregulation of inflammatory genes after exposure to endotoxin. Our results demonstrate that CTCF-dependent genome topology is not strictly required for a functional cell-fate conversion but facilitates a rapid and efficient response to an external stimulus.

How pluripotent stem cells differentiate into the main germ layers is a key question of developmental biology. Here, we show that the chromatin-related factor Whsc1 (also known as Nsd2 and MMSET) has a dual role in pluripotency exit and germ layer specification of embryonic stem cells. On induction of differentiation, a proportion of Whsc1-depleted embryonic stem cells remain entrapped in a pluripotent state and fail to form mesendoderm, although they are still capable of generating neuroectoderm. These functions of Whsc1 are independent of its methyltransferase activity. Whsc1 binds to enhancers of the mesendodermal regulators Gata4, T (Brachyury), Gata6 and Foxa2, together with Brd4, and activates the expression of these genes. Depleting each of these regulators also delays pluripotency exit, suggesting that they mediate the effects observed with Whsc1. Our data indicate that Whsc1 links silencing of the pluripotency regulatory network with activation of mesendoderm lineages.

Chromosome conformation capture techniques, such as Hi-C, are fundamental in characterizing genome organization. These methods have revealed several genomic features, such as chromatin loops, whose disruption can have dramatic effects in gene regulation. Unfortunately, their detection is difficult; current methods require that the users choose the resolution of interaction maps based on dataset quality and sequencing depth. Here, we introduce Binless, a resolution-agnostic method that adapts to the quality and quantity of available data, to detect both interactions and differences. Binless relies on an alternate representation of Hi-C data, which leads to a more detailed classification of paired-end reads. Using a large-scale benchmark, we demonstrate that Binless is able to call interactions with higher reproducibility than other existing methods. Binless, which is freely available, can thus reliably be used to identify chromatin loops as well as for differential analysis of chromatin interaction maps.

In breast cancer cells, some topologically associating domains (TADs) behave as hormonal gene regulation units, within which gene transcription is coordinately regulated in response to steroid hormones. Here we further describe that responsive TADs contain 20- to 100-kb-long clusters of intermingled estrogen receptor (ESR1) and progesterone receptor (PGR) binding sites, hereafter called hormone-control regions (HCRs). In T47D cells, we identified more than 200 HCRs, which are frequently bound by unliganded ESR1 and PGR. These HCRs establish steady long-distance inter-TAD interactions between them and organize characteristic looping structures with promoters in their TADs even in the absence of hormones in ESR1⁺-PGR⁺ cells. This organization is dependent on the expression of the receptors and is further dynamically modulated in response to steroid hormones. HCRs function as platforms that integrate different signals, resulting in some cases in opposite transcriptional responses to estrogens or progestins. Altogether, these results suggest that steroid hormone receptors act not only as hormone-regulated sequence-specific transcription factors but also as local and global genome organizers.

Nuclear architecture is decisive for the assembly of transcriptional responses. However, how chromosome organization is dynamically modulated to permit rapid and transient transcriptional changes in response to environmental challenges remains unclear. Here we show that hyperosmotic stress disrupts different levels of chromosome organization, ranging from A/B compartment changes to reduction in the number and insulation of topologically associating domains (TADs). Concomitantly, transcription is greatly affected, TAD borders weaken, and RNA Polymerase II runs off from hundreds of transcription end sites. Stress alters the binding profiles of architectural proteins, which explains the disappearance of local chromatin organization. These processes are dynamic, and cells rapidly reconstitute their default chromatin conformation after stress removal, uncovering an intrinsic organization. Transcription is not required for local chromatin reorganization, while compartment recovery is partially transcription-dependent. Thus, nuclear organization in mammalian cells can be rapidly modulated by environmental changes in a reversible manner.