LATS1 controls CTCF chromatin occupancy and hormonal response of 3D-grown breast cancer cells.

The cancer epigenome has been studied in cells cultured in two-dimensional (2D) monolayers, but recent studies highlight the impact of the extracellular matrix and the three-dimensional (3D) environment on multiple cellular functions. Here, we report the physical, biochemical, and genomic differences between T47D breast cancer cells cultured in 2D and as 3D spheroids. Cells within 3D spheroids exhibit a rounder nucleus with less accessible, more compacted chromatin, as well as altered expression of ~2000 genes, the majority of which become repressed. Hi-C analysis reveals that cells in 3D are enriched for regions belonging to the B compartment, have decreased chromatin-bound CTCF and increased fusion of topologically associating domains (TADs). Upregulation of the Hippo pathway in 3D spheroids results in the activation of the LATS1 kinase, which promotes phosphorylation and displacement of CTCF from DNA, thereby likely causing the observed TAD fusions. 3D cells show higher chromatin binding of progesterone receptor (PR), leading to an increase in the number of hormone-regulated genes. This effect is in part mediated by LATS1 activation, which favors cytoplasmic retention of YAP and CTCF removal.

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.

Chromatin topology defines estradiol-primed progesterone receptor and PAX2 binding in endometrial cancer cells.

Estrogen (E2) and Progesterone (Pg), via their specific receptors (ERalpha and PR), are major determinants in the development and progression of endometrial carcinomas, However, their precise mechanism of action and the role of other transcription factors involved are not entirely clear. Using Ishikawa endometrial cancer cells, we report that E2 treatment exposes a set of progestin-dependent PR binding sites which include both E2 and progestin target genes. ChIP-seq results from hormone-treated cells revealed a non-random distribution of PAX2 binding in the vicinity of these estrogen-promoted PR sites. Altered expression of hormone regulated genes in PAX2 knockdown cells suggests a role for PAX2 in fine-tuning ERalpha and PR interplay in transcriptional regulation. Analysis of long-range interactions by Hi-C coupled with ATAC-seq data showed that these regions, that we call 'progestin control regions' (PgCRs), exhibited an open chromatin state even before hormone exposure and were non-randomly associated with regulated genes. Nearly 20% of genes potentially influenced by PgCRs were found to be altered during progression of endometrial cancer. Our findings suggest that endometrial response to progestins in differentiated endometrial tumor cells results in part from binding of PR together with PAX2 to accessible chromatin regions. What maintains these regions open remains to be studied.

A set of accessible enhancers enables the initial response of breast cancer cells to physiological progestin concentrations.

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.

STAG2 loss-of-function affects short-range genomic contacts and modulates the basal-luminal transcriptional program of bladder cancer cells.

Cohesin exists in two variants containing STAG1 or STAG2. STAG2 is one of the most mutated genes in cancer and a major bladder tumor suppressor. Little is known about how its inactivation contributes to tumorigenesis. Here, we analyze the genomic distribution of STAG1 and STAG2 and perform STAG2 loss-of-function experiments using RT112 bladder cancer cells; we then analyze the genomic effects by integrating gene expression and chromatin interaction data. Functional compartmentalization exists between the cohesin complexes: cohesin-STAG2 displays a distinctive genomic distribution and mediates short and mid-ranged interactions that engage genes at higher frequency than those established by cohesin-STAG1. STAG2 knockdown results in down-regulation of the luminal urothelial signature and up-regulation of the basal transcriptional program, mirroring differences between STAG2-high and STAG2-low human bladder tumors. This is accompanied by rewiring of DNA contacts within topological domains, while compartments and domain boundaries remain refractive. Contacts lost upon depletion of STAG2 are assortative, preferentially occur within silent chromatin domains, and are associated with de-repression of lineage-specifying genes. Our findings indicate that STAG2 participates in the DNA looping that keeps the basal transcriptional program silent and thus sustains the luminal program. This mechanism may contribute to the tumor suppressor function of STAG2 in the urothelium.

MyoD induces ARTD1 and nucleoplasmic poly-ADP-ribosylation during fibroblast to myoblast transdifferentiation.

While protein ADP-ribosylation was reported to regulate differentiation and dedifferentiation, it has so far not been studied during transdifferentiation. Here, we found that MyoD-induced transdifferentiation of fibroblasts to myoblasts promotes the expression of the ADP-ribosyltransferase ARTD1. Comprehensive analysis of the genome architecture by Hi-C and RNA-seq analysis during transdifferentiation indicated that ARTD1 locally contributed to A/B compartmentalization and coregulated a subset of MyoD target genes that were however not sufficient to alter transdifferentiation. Surprisingly, the expression of ARTD1 was accompanied by the continuous synthesis of nuclear ADP ribosylation that was neither dependent on the cell cycle nor induced by DNA damage. Conversely to the H2O2-induced ADP-ribosylation, the MyoD-dependent ADP-ribosylation was not associated to chromatin but rather localized to the nucleoplasm. Together, these data describe a MyoD-induced nucleoplasmic ADP-ribosylation that is observed particularly during transdifferentiation and thus potentially expands the plethora of cellular processes associated with ADP-ribosylation.

The impact of chromosomal fusions on 3D genome folding and recombination in the germ line.

The spatial folding of chromosomes inside the nucleus has regulatory effects on gene expression, yet the impact of genome reshuffling on this organization remains unclear. Here, we take advantage of chromosome conformation capture in combination with single-nucleotide polymorphism (SNP) genotyping and analysis of crossover events to study how the higher-order chromatin organization and recombination landscapes are affected by chromosomal fusions in the mammalian germ line. We demonstrate that chromosomal fusions alter the nuclear architecture during meiosis, including an increased rate of heterologous interactions in primary spermatocytes, and alterations in both chromosome synapsis and axis length. These disturbances in topology were associated with changes in genomic landscapes of recombination, resulting in detectable genomic footprints. Overall, we show that chromosomal fusions impact the dynamic genome topology of germ cells in two ways: (i) altering chromosomal nuclear occupancy and synapsis, and (ii) reshaping landscapes of recombination.

CTCF is dispensable for immune cell transdifferentiation but facilitates an acute inflammatory response.

Three-dimensional organization of the genome is important for transcriptional regulation1-7. In mammals, CTCF and the cohesin complex create submegabase structures with elevated internal chromatin contact frequencies, called topologically associating domains (TADs)8-12. Although TADs can contribute to transcriptional regulation, ablation of TAD organization by disrupting CTCF or the cohesin complex causes modest gene expression changes13-16. In contrast, CTCF is required for cell cycle regulation17, embryonic development and formation of various adult cell types18. 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.

C/EBPα mediates the growth inhibitory effect of progestins on breast cancer cells.

Steroid hormones are key gene regulators in breast cancer cells. While estrogens stimulate cell proliferation, progestins activate a single cell cycle followed by proliferation arrest. Here, we use biochemical and genome-wide approaches to show that progestins achieve this effect via a functional crosstalk with C/EBPα. Using ChIP-seq, we identify around 1,000 sites where C/EBPα binding precedes and helps binding of progesterone receptor (PR) in response to hormone. These regions exhibit epigenetic marks of active enhancers, and C/EBPα maintains an open chromatin conformation that facilitates loading of ligand-activated PR. Prior to hormone exposure, C/EBPα favors promoter-enhancer contacts that assure hormonal regulation of key genes involved in cell proliferation by facilitating binding of RAD21, YY1, and the Mediator complex. Knockdown of C/EBPα disrupts enhancer-promoter contacts and decreases the presence of these architectural proteins, highlighting its key role in 3D chromatin looping. Thus, C/EBPα fulfills a previously unknown function as a potential growth modulator in hormone-dependent breast cancer.

Three-Dimensional Genomic Structure and Cohesin Occupancy Correlate with Transcriptional Activity during Spermatogenesis.

Mammalian gametogenesis involves dramatic and tightly regulated chromatin remodeling, whose regulatory pathways remain largely unexplored. Here, we generate a comprehensive high-resolution structural and functional atlas of mouse spermatogenesis by combining in situ chromosome conformation capture sequencing (Hi-C), RNA sequencing (RNA-seq), and chromatin immunoprecipitation sequencing (ChIP-seq) of CCCTC-binding factor (CTCF) and meiotic cohesins, coupled with confocal and super-resolution microscopy. Spermatogonia presents well-defined compartment patterns and topological domains. However, chromosome occupancy and compartmentalization are highly re-arranged during prophase I, with cohesins bound to active promoters in DNA loops out of the chromosomal axes. Compartment patterns re-emerge in round spermatids, where cohesin occupancy correlates with transcriptional activity of key developmental genes. The compact sperm genome contains compartments with actively transcribed genes but no fine-scale topological domains, concomitant with the presence of protamines. Overall, we demonstrate how genome-wide cohesin occupancy and transcriptional activity is associated with three-dimensional (3D) remodeling during spermatogenesis, ultimately reprogramming the genome for the next generation.

ATP, Mg2+, Nuclear Phase Separation, and Genome Accessibility.

Misregulation of the processes controlling eukaryotic gene expression can result in disease. Gene expression is influenced by the surrounding chromatin; hence the nuclear environment is also of vital importance. Recently, understanding of chromatin hierarchical folding has increased together with the discovery of membrane-less organelles which are distinct, dynamic liquid droplets that merge and expand within the nucleus. These 'sieve'-like regions may compartmentalize and separate functionally distinct regions of chromatin. This article aims to discuss recent studies on nuclear phase within the context of poly(ADP-ribose), ATP, and Mg2+ levels, and we propose a combinatorial complex role for these molecules in phase separation and genome regulation. We also discuss the implications of this process for gene regulation and discuss possible strategies to test this.

Rapid reversible changes in compartments and local chromatin organization revealed by hyperosmotic shock.

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.

Arginine Citrullination at the C-Terminal Domain Controls RNA Polymerase II Transcription.

The post-translational modification of key residues at the C-terminal domain of RNA polymerase II (RNAP2-CTD) coordinates transcription, splicing, and RNA processing by modulating its capacity to act as a landing platform for a variety of protein complexes. Here, we identify a new modification at the CTD, the deimination of arginine and its conversion to citrulline by peptidyl arginine deiminase 2 (PADI2), an enzyme that has been associated with several diseases, including cancer. We show that, among PADI family members, only PADI2 citrullinates R1810 (Cit1810) at repeat 31 of the CTD. Depletion of PADI2 or loss of R1810 results in accumulation of RNAP2 at transcription start sites, reduced gene expression, and inhibition of cell proliferation. Cit1810 is needed for interaction with the P-TEFb (positive transcription elongation factor b) kinase complex and for its recruitment to chromatin. In this way, CTD-Cit1810 favors RNAP2 pause release and efficient transcription in breast cancer cells.

Hormone-control regions mediate steroid receptor-dependent genome organization.

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.

Unliganded Progesterone Receptor Governs Estrogen Receptor Gene Expression by Regulating DNA Methylation in Breast Cancer Cells.

Breast cancer prognosis and response to endocrine therapy strongly depends on the expression of the estrogen and progesterone receptors (ER and PR, respectively). Although much is known about ERα gene (ESR1) regulation after hormonal stimulation, how it is regulated in hormone-free condition is not fully understood. We used ER-/PR-positive breast cancer cells to investigate the role of PR in ESR1 regulation in the absence of hormones. We show that PR binds to the low-methylated ESR1 promoter and maintains both gene expression and DNA methylation of the ESR1 locus in hormone-deprived breast cancer cells. Depletion of PR reduces ESR1 expression, with a concomitant increase in gene promoter methylation. The high amount of methylation in the ESR1 promoter of PR-depleted cells persists after the stable re-expression of PR and inhibits PR binding to this genomic region. As a consequence, the rescue of PR expression in PR-depleted cells is insufficient to restore ESR1 expression. Consistently, DNA methylation impedes PR binding to consensus progesterone responsive elements. These findings contribute to understanding the complex crosstalk between PR and ER and suggest that the analysis of ESR1 promoter methylation in breast cancer cells can help to design more appropriate targeted therapies for breast cancer patients.

Promoter bivalency favors an open chromatin architecture in embryonic stem cells.

In embryonic stem cells (ESCs), developmental gene promoters are characterized by their bivalent chromatin state, with simultaneous modification by MLL2 and Polycomb complexes. Although essential for embryogenesis, bivalency is functionally not well understood. Here, we show that MLL2 plays a central role in ESC genome organization. We generate a catalog of bona fide bivalent genes in ESCs and demonstrate that loss of MLL2 leads to increased Polycomb occupancy. Consequently, promoters lose accessibility, long-range interactions are redistributed, and ESCs fail to differentiate. We pose that bivalency balances accessibility and long-range connectivity of promoters, allowing developmental gene expression to be properly modulated.

Lamin B1 mapping reveals the existence of dynamic and functional euchromatin lamin B1 domains.

Lamins (A/C and B) are major constituents of the nuclear lamina (NL). Structurally conserved lamina-associated domains (LADs) are formed by genomic regions that contact the NL. Lamins are also found in the nucleoplasm, with a yet unknown function. Here we map the genome-wide localization of lamin B1 in an euchromatin-enriched fraction of the mouse genome and follow its dynamics during the epithelial-to-mesenchymal transition (EMT). Lamin B1 associates with actively expressed and open euchromatin regions, forming dynamic euchromatin lamin B1-associated domains (eLADs) of about 0.3 Mb. Hi-C data link eLADs to the 3D organization of the mouse genome during EMT and correlate lamin B1 enrichment at topologically associating domain (TAD) borders with increased border strength. Having reduced levels of lamin B1 alters the EMT transcriptional signature and compromises the acquisition of mesenchymal traits. Thus, during EMT, the process of genome reorganization in mouse involves dynamic changes in eLADs.

Distinct roles of cohesin-SA1 and cohesin-SA2 in 3D chromosome organization.

Two variant cohesin complexes containing SMC1, SMC3, RAD21 and either SA1 (also known as STAG1) or SA2 (also known as STAG2) are present in all cell types. We report here their genomic distribution and specific contributions to genome organization in human cells. Although both variants are found at CCCTC-binding factor (CTCF) sites, a distinct population of the SA2-containing cohesin complexes (hereafter referred to as cohesin-SA2) localize to enhancers lacking CTCF, are linked to tissue-specific transcription and cannot be replaced by the SA1-containing cohesin complex (cohesin-SA1) when SA2 is absent, a condition that has been observed in several tumors. Downregulation of each of these variants has different consequences for gene expression and genome architecture. Our results suggest that cohesin-SA1 preferentially contributes to the stabilization of topologically associating domain boundaries together with CTCF, whereas cohesin-SA2 promotes cell-type-specific contacts between enhancers and promoters independently of CTCF. Loss of cohesin-SA2 rewires local chromatin contacts and alters gene expression. These findings provide insights into how cohesin mediates chromosome folding and establish a novel framework to address the consequences of mutations in cohesin genes in cancer.

OneD: increasing reproducibility of Hi-C samples with abnormal karyotypes.

The three-dimensional conformation of genomes is an essential component of their biological activity. The advent of the Hi-C technology enabled an unprecedented progress in our understanding of genome structures. However, Hi-C is subject to systematic biases that can compromise downstream analyses. Several strategies have been proposed to remove those biases, but the issue of abnormal karyotypes received little attention. Many experiments are performed in cancer cell lines, which typically harbor large-scale copy number variations that create visible defects on the raw Hi-C maps. The consequences of these widespread artifacts on the normalized maps are mostly unexplored. We observed that current normalization methods are not robust to the presence of large-scale copy number variations, potentially obscuring biological differences and enhancing batch effects. To address this issue, we developed an alternative approach designed to take into account chromosomal abnormalities. The method, called OneD, increases reproducibility among replicates of Hi-C samples with abnormal karyotype, outperforming previous methods significantly. On normal karyotypes, OneD fared equally well as state-of-the-art methods, making it a safe choice for Hi-C normalization. OneD is fast and scales well in terms of computing resources for resolutions up to 5 kb.