A recurrent chromosomal inversion suffices for driving escape from oncogene-induced senescence via subTAD reorganization.

Oncogene-induced senescence (OIS) is an inherent and important tumor suppressor mechanism. However, if not removed timely via immune surveillance, senescent cells also have detrimental effects. Although this has mostly been attributed to the senescence-associated secretory phenotype (SASP) of these cells, we recently proposed that "escape" from the senescent state is another unfavorable outcome. The mechanism underlying this phenomenon remains elusive. Here, we exploit genomic and functional data from a prototypical human epithelial cell model carrying an inducible CDC6 oncogene to identify an early-acquired recurrent chromosomal inversion that harbors a locus encoding the circadian transcription factor BHLHE40. This inversion alone suffices for BHLHE40 activation upon CDC6 induction and driving cell cycle re-entry of senescent cells, and malignant transformation. Ectopic overexpression of BHLHE40 prevented induction of CDC6-triggered senescence. We provide strong evidence in support of replication stress-induced genomic instability being a causative factor underlying "escape" from oncogene-induced senescence.

HMGB2 Loss upon Senescence Entry Disrupts Genomic Organization and Induces CTCF Clustering across Cell Types.

Processes like cellular senescence are characterized by complex events giving rise to heterogeneous cell populations. However, the early molecular events driving this cascade remain elusive. We hypothesized that senescence entry is triggered by an early disruption of the cells' three-dimensional (3D) genome organization. To test this, we combined Hi-C, single-cell and population transcriptomics, imaging, and in silico modeling of three distinct cells types entering senescence. Genes involved in DNA conformation maintenance are suppressed upon senescence entry across all cell types. We show that nuclear depletion of the abundant HMGB2 protein occurs early on the path to senescence and coincides with the dramatic spatial clustering of CTCF. Knocking down HMGB2 suffices for senescence-induced CTCF clustering and for loop reshuffling, while ectopically expressing HMGB2 rescues these effects. Our data suggest that HMGB2-mediated genomic reorganization constitutes a primer for the ensuing senescent program.

HMGB1 coordinates SASP-related chromatin folding and RNA homeostasis on the path to senescence.

Spatial organization and gene expression of mammalian chromosomes are maintained and regulated in conjunction with cell cycle progression. This is perturbed once cells enter senescence and the highly abundant HMGB1 protein is depleted from nuclei to act as an extracellular proinflammatory stimulus. Despite its physiological importance, we know little about the positioning of HMGB1 on chromatin and its nuclear roles. To address this, we mapped HMGB1 binding genome-wide in two primary cell lines. We integrated ChIP-seq and Hi-C with graph theory to uncover clustering of HMGB1-marked topological domains that harbor genes involved in paracrine senescence. Using simplified Cross-Linking and Immuno-Precipitation and functional tests, we show that HMGB1 is also a bona fide RNA-binding protein (RBP) binding hundreds of mRNAs. It presents an interactome rich in RBPs implicated in senescence regulation. The mRNAs of many of these RBPs are directly bound by HMGB1 and regulate availability of SASP-relevant transcripts. Our findings reveal a broader than hitherto assumed role for HMGB1 in coordinating chromatin folding and RNA homeostasis as part of a regulatory loop controlling cell-autonomous and paracrine senescence.

Binding of nuclear factor κB to noncanonical consensus sites reveals its multimodal role during the early inflammatory response.

Mammalian cells have developed intricate mechanisms to interpret, integrate, and respond to extracellular stimuli. For example, tumor necrosis factor (TNF) rapidly activates proinflammatory genes, but our understanding of how this occurs against the ongoing transcriptional program of the cell is far from complete. Here, we monitor the early phase of this cascade at high spatiotemporal resolution in TNF-stimulated human endothelial cells. NF-κB, the transcription factor complex driving the response, interferes with the regulatory machinery by binding active enhancers already in interaction with gene promoters. Notably, >50% of these enhancers do not encode canonical NF-κB binding motifs. Using a combination of genomics tools, we find that binding site selection plays a key role in NF-κΒ-mediated transcriptional activation and repression. We demonstrate the latter by describing the synergy between NF-κΒ and the corepressor JDP2. Finally, detailed analysis of a 2.8-Mbp locus using sub-kbp-resolution targeted chromatin conformation capture and genome editing uncovers how NF-κΒ that has just entered the nucleus exploits pre-existing chromatin looping to exert its multimodal role. This work highlights the involvement of topology in cis-regulatory element function during acute transcriptional responses, where primary DNA sequence and its higher-order structure constitute a regulatory context leading to either gene activation or repression.

Analysis of computational footprinting methods for DNase sequencing experiments.

DNase-seq allows nucleotide-level identification of transcription factor binding sites on the basis of a computational search of footprint-like DNase I cleavage patterns on the DNA. Frequently in high-throughput methods, experimental artifacts such as DNase I cleavage bias affect the computational analysis of DNase-seq experiments. Here we performed a comprehensive and systematic study on the performance of computational footprinting methods. We evaluated ten footprinting methods in a panel of DNase-seq experiments for their ability to recover cell-specific transcription factor binding sites. We show that three methods--HINT, DNase2TF and PIQ--consistently outperformed the other evaluated methods and that correcting the DNase-seq signal for experimental artifacts significantly improved the accuracy of computational footprints. We also propose a score that can be used to detect footprints arising from transcription factors with potentially short residence times.

Epigenetic program and transcription factor circuitry of dendritic cell development.

Dendritic cells (DC) are professional antigen presenting cells that develop from hematopoietic stem cells through successive steps of lineage commitment and differentiation. Multipotent progenitors (MPP) are committed to DC restricted common DC progenitors (CDP), which differentiate into specific DC subsets, classical DC (cDC) and plasmacytoid DC (pDC). To determine epigenetic states and regulatory circuitries during DC differentiation, we measured consecutive changes of genome-wide gene expression, histone modification and transcription factor occupancy during the sequel MPP-CDP-cDC/pDC. Specific histone marks in CDP reveal a DC-primed epigenetic signature, which is maintained and reinforced during DC differentiation. Epigenetic marks and transcription factor PU.1 occupancy increasingly coincide upon DC differentiation. By integrating PU.1 occupancy and gene expression we devised a transcription factor regulatory circuitry for DC commitment and subset specification. The circuitry provides the transcription factor hierarchy that drives the sequel MPP-CDP-cDC/pDC, including Irf4, Irf8, Tcf4, Spib and Stat factors. The circuitry also includes feedback loops inferred for individual or multiple factors, which stabilize distinct stages of DC development and DC subsets. In summary, here we describe the basic regulatory circuitry of transcription factors that drives DC development.

The interaction of MYC with the trithorax protein ASH2L promotes gene transcription by regulating H3K27 modification.

The appropriate expression of the roughly 30,000 human genes requires multiple layers of control. The oncoprotein MYC, a transcriptional regulator, contributes to many of the identified control mechanisms, including the regulation of chromatin, RNA polymerases, and RNA processing. Moreover, MYC recruits core histone-modifying enzymes to DNA. We identified an additional transcriptional cofactor complex that interacts with MYC and that is important for gene transcription. We found that the trithorax protein ASH2L and MYC interact directly in vitro and co-localize in cells and on chromatin. ASH2L is a core subunit of KMT2 methyltransferase complexes that target histone H3 lysine 4 (H3K4), a mark associated with open chromatin. Indeed, MYC associates with H3K4 methyltransferase activity, dependent on the presence of ASH2L. MYC does not regulate this methyltransferase activity but stimulates demethylation and subsequently acetylation of H3K27. KMT2 complexes have been reported to associate with histone H3K27-specific demethylases, while CBP/p300, which interact with MYC, acetylate H3K27. Finally WDR5, another core subunit of KMT2 complexes, also binds directly to MYC and in genome-wide analyses MYC and WDR5 are associated with transcribed promoters. Thus, our findings suggest that MYC and ASH2L-KMT2 complexes cooperate in gene transcription by controlling H3K27 modifications and thereby regulate bivalent chromatin.

Detection of active transcription factor binding sites with the combination of DNase hypersensitivity and histone modifications.

MOTIVATION: The identification of active transcriptional regulatory elements is crucial to understand regulatory networks driving cellular processes such as cell development and the onset of diseases. It has recently been shown that chromatin structure information, such as DNase I hypersensitivity (DHS) or histone modifications, significantly improves cell-specific predictions of transcription factor binding sites. However, no method has so far successfully combined both DHS and histone modification data to perform active binding site prediction. RESULTS: We propose here a method based on hidden Markov models to integrate DHS and histone modifications occupancy for the detection of open chromatin regions and active binding sites. We have created a framework that includes treatment of genomic signals, model training and genome-wide application. In a comparative analysis, our method obtained a good trade-off between sensitivity versus specificity and superior area under the curve statistics than competing methods. Moreover, our technique does not require further training or sequence information to generate binding location predictions. Therefore, the method can be easily applied on new cell types and allow flexible downstream analysis such as de novo motif finding. AVAILABILITY AND IMPLEMENTATION: Our framework is available as part of the Regulatory Genomics Toolbox. The software information and all benchmarking data are available at http://costalab.org/wp/dh-hmm. CONTACT: ivan.costa@rwth-aachen.de or eduardo.gusmao@rwth-aachen.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A ONECUT1 regulatory, non-coding region in pancreatic development and diabetes.

In a patient with permanent neonatal syndromic diabetes clinically similar to cases with ONECUT1 biallelic mutations, we identified a disease-causing deletion located upstream of ONECUT1. Through genetic, genomic, and functional studies, we identified a crucial regulatory region acting as an enhancer of ONECUT1 specifically during pancreatic development. This enhancer region contains a low-frequency variant showing a strong association with type 2 diabetes and other glycemic traits, thus extending the contribution of this region to common forms of diabetes. Clinical relevance is provided by experimentally tailored therapy options for patients carrying ONECUT1 coding or regulatory mutations.

Replicative senescence is associated with nuclear reorganization and with DNA methylation at specific transcription factor binding sites.

BACKGROUND: Primary cells enter replicative senescence after a limited number of cell divisions. This process needs to be considered in cell culture experiments, and it is particularly important for regenerative medicine. Replicative senescence is associated with reproducible changes in DNA methylation (DNAm) at specific sites in the genome. The mechanism that drives senescence-associated DNAm changes remains unknown - it may involve stochastic DNAm drift due to imperfect maintenance of epigenetic marks or it is directly regulated at specific sites in the genome. RESULTS: In this study, we analyzed the reorganization of nuclear architecture and DNAm changes during long-term culture of human fibroblasts and mesenchymal stromal cells (MSCs). We demonstrate that telomeres shorten and shift towards the nuclear center at later passages. In addition, DNAm profiles, either analyzed by MethylCap-seq or by 450k IlluminaBeadChip technology, revealed consistent senescence-associated hypermethylation in regions associated with H3K27me3, H3K4me3, and H3K4me1 histone marks, whereas hypomethylation was associated with chromatin containing H3K9me3 and lamina-associated domains (LADs). DNA hypermethylation was significantly enriched in the vicinity of genes that are either up- or downregulated at later passages. Furthermore, specific transcription factor binding motifs (e.g. EGR1, TFAP2A, and ETS1) were significantly enriched in differentially methylated regions and in the promoters of differentially expressed genes. CONCLUSIONS: Senescence-associated DNA hypermethylation occurs at specific sites in the genome and reflects functional changes in the course of replicative senescence. These results indicate that tightly regulated epigenetic modifications during long-term culture contribute to changes in nuclear organization and gene expression.