In vivo temporal resolution of acute promyelocytic leukemia progression reveals a role of Klf4 in suppressing early leukemic transformation.

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.

Coordinated changes in gene expression, H1 variant distribution and genome 3D conformation in response to H1 depletion.

Up to seven members of the histone H1 family may contribute to chromatin compaction and its regulation in human somatic cells. In breast cancer cells, knock-down of multiple H1 variants deregulates many genes, promotes the appearance of genome-wide accessibility sites and triggers an interferon response via activation of heterochromatic repeats. However, how these changes in the expression profile relate to the re-distribution of H1 variants as well as to genome conformational changes have not been yet studied. Here, we combined ChIP-seq of five endogenous H1 variants with Chromosome Conformation Capture analysis in wild-type and H1.2/H1.4 knock-down T47D cells. The results indicate that H1 variants coexist in the genome in two large groups depending on the local GC content and that their distribution is robust with respect to H1 depletion. Despite the small changes in H1 variants distribution, knock-down of H1 translated into more isolated but de-compacted chromatin structures at the scale of topologically associating domains (TADs). Such changes in TAD structure correlated with a coordinated gene expression response of their resident genes. This is the first report describing simultaneous profiling of five endogenous H1 variants and giving functional evidence of genome topology alterations upon H1 depletion in human cancer cells.

TADs enriched in histone H1.2 strongly overlap with the B compartment, inaccessible chromatin, and AT-rich Giemsa bands.

Giemsa staining of metaphase chromosomes results in a characteristic banding useful for identification of chromosomes and its alterations. We have investigated in silico whether Giemsa bands (G bands) correlate with epigenetic and topological features of the interphase genome. Staining of G-positive bands decreases with GC content; nonetheless, G-negative bands are GC heterogeneous. High GC bands are enriched in active histone marks, RNA polymerase II, and SINEs and associate with gene richness, gene expression, and early replication. Low GC bands are enriched in repressive marks, lamina-associated domains, and LINEs. Histone H1 variants distribute heterogeneously among G bands: H1X is enriched at high GC bands and H1.2 is abundant at low GC, compacted bands. According to epigenetic features and H1 content, G bands can be organized in clusters useful to compartmentalize the genome. Indeed, we have obtained Hi-C chromosome interaction maps and compared topologically associating domains (TADs) and A/B compartments to G banding. TADs with high H1.2/H1X ratio strongly overlap with B compartment, late replicating, and inaccessible chromatin and low GC bands. We propose that GC content is a strong driver of chromatin compaction and 3D genome organization, that Giemsa staining recapitulates this organization denoted by high-throughput techniques, and that H1 variants distribute at distinct chromatin domains. DATABASES: Hi-C data on T47D breast cancer cells have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE147627.