Crosstalk between chromatin and Shavenbaby defines transcriptional output along the Drosophila intestinal stem cell lineage.

The transcription factor Shavenbaby (Svb), the only member of the OvoL family in Drosophila, controls the fate of various epithelial embryonic cells and adult stem cells. Post-translational modification of Svb produces two protein isoforms, Svb-ACT and Svb-REP, which promote adult intestinal stem cell renewal or differentiation, respectively. To define Svb mode of action, we used engineered cell lines and develop an unbiased method to identify Svb target genes across different contexts. Within a given cell type, Svb-ACT and Svb-REP antagonistically regulate the expression of a set of target genes, binding specific enhancers whose accessibility is constrained by chromatin landscape. Reciprocally, Svb-REP can influence local chromatin marks of active enhancers to help repressing target genes. Along the intestinal lineage, the set of Svb target genes progressively changes, together with chromatin accessibility. We propose that Svb-ACT-to-REP transition promotes enterocyte differentiation of intestinal stem cells through direct gene regulation and chromatin remodeling.

POLE3 is a repressor of unintegrated HIV-1 DNA required for efficient virus integration and escape from innate immune sensing.

Unintegrated retroviral DNA is transcriptionally silenced by host chromatin silencing factors. Here, we used the proteomics of isolated chromatin segments method to reveal viral and host factors associated with unintegrated HIV-1DNA involved in its silencing. By gene silencing using siRNAs, 46 factors were identified as potential repressors of unintegrated HIV-1DNA. Knockdown and knockout experiments revealed POLE3 as a transcriptional repressor of unintegrated HIV-1DNA. POLE3 maintains unintegrated HIV-1DNA in a repressive chromatin state, preventing RNAPII recruitment to the viral promoter. POLE3 and the recently identified host factors mediating unintegrated HIV-1 DNA silencing, CAF1 and SMC5/SMC6/SLF2, show specificity toward different forms of unintegrated HIV-1DNA. Loss of POLE3 impaired HIV-1 replication, suggesting that repression of unintegrated HIV-1DNA is important for optimal viral replication. POLE3 depletion reduces the integration efficiency of HIV-1. POLE3, by maintaining a repressive chromatin structure of unintegrated HIV-1DNA, ensures HIV-1 escape from innate immune sensing in primary CD4+ T cells.

Condensin positioning at telomeres by shelterin proteins drives sister-telomere disjunction in anaphase.

The localization of condensin along chromosomes is crucial for their accurate segregation in anaphase. Condensin is enriched at telomeres but how and for what purpose had remained elusive. Here, we show that fission yeast condensin accumulates at telomere repeats through the balancing acts of Taz1, a core component of the shelterin complex that ensures telomeric functions, and Mit1, a nucleosome remodeler associated with shelterin. We further show that condensin takes part in sister-telomere separation in anaphase, and that this event can be uncoupled from the prior separation of chromosome arms, implying a telomere-specific separation mechanism. Consistent with a cis-acting process, increasing or decreasing condensin occupancy specifically at telomeres modifies accordingly the efficiency of their separation in anaphase. Genetic evidence suggests that condensin promotes sister-telomere separation by counteracting cohesin. Thus, our results reveal a shelterin-based mechanism that enriches condensin at telomeres to drive in cis their separation during mitosis.

PAPγ associates with PAXT nuclear exosome to control the abundance of PROMPT ncRNAs.

Pervasive transcription of the human genome generates an abundance of RNAs that must be processed and degraded. The nuclear RNA exosome is the main RNA degradation machinery in the nucleus. However, nuclear exosome must be recruited to its substrates by targeting complexes, such as NEXT or PAXT. By proteomic analysis, we identify additional subunits of PAXT, including many orthologs of MTREC found in S. pombe. In particular, we show that polyA polymerase gamma (PAPγ) associates with PAXT. Genome-wide mapping of the binding sites of ZFC3H1, RBM27 and PAPγ shows that PAXT is recruited to the TSS of hundreds of genes. Loss of ZFC3H1 abolishes recruitment of PAXT subunits including PAPγ to TSSs and concomitantly increases the abundance of PROMPTs at the same sites. Moreover, PAPγ, as well as MTR4 and ZFC3H1, is implicated in the polyadenylation of PROMPTs. Our results thus provide key insights into the direct targeting of PROMPT ncRNAs by PAXT at their genomic sites.

Chromatin in 3D distinguishes dMes-4/NSD and Hypb/dSet2 in protecting genes from H3K27me3 silencing.

Cell type-specific barcoding of genomes requires the establishment of hundreds of heterochromatin domains where heterochromatin-associated repressive complexes hinder chromatin accessibility thereby silencing genes. At heterochromatin-euchromatin borders, regulation of accessibility not only depends on the delimitation of heterochromatin but may also involve interplays with nearby genes and their transcriptional activity, or alternatively on histone modifiers, chromatin barrier insulators, and more global demarcation of chromosomes into 3D compartmentalized domains and topological-associating domain (TADs). Here, we show that depletion of H3K36 di- or tri-methyl histone methyltransferases dMes-4/NSD or Hypb/dSet2 induces reproducible increasing levels of H3K27me3 at heterochromatin borders including in nearby promoters, thereby repressing hundreds of genes. Furthermore, dMes-4/NSD influences genes demarcated by insulators and TAD borders, within chromatin hubs, unlike transcription-coupled action of Hypb/dSet2 that protects genes independently of TADs. Insulator mutants recapitulate the increase of H3K27me3 upon dMes-4/NSD depletion unlike Hypb/dSet2. Hi-C data demonstrate how dMes-4/NSD blocks propagation of long-range interactions onto active regions. Our data highlight distinct mechanisms protecting genes from H3K27me3 silencing, highlighting a direct influence of H3K36me on repressive TADs.

The histone variant macroH2A1.1 regulates RNA polymerase II-paused genes within defined chromatin interaction landscapes.

The histone variant macroH2A1.1 plays a role in cancer development and metastasis. To determine the underlying molecular mechanisms, we mapped the genome-wide localization of endogenous macroH2A1.1 in the human breast cancer cell line MDA-MB-231. We demonstrate that macroH2A1.1 specifically binds to active promoters and enhancers in addition to facultative heterochromatin. Selective knock down of macroH2A1.1 deregulates the expression of hundreds of highly active genes. Depending on the chromatin landscape, macroH2A1.1 acts through two distinct molecular mechanisms. The first mitigates excessive transcription by binding over domains including the promoter and the gene body. The second stimulates expression of RNA polymerase II (Pol II)-paused genes, including genes regulating mammary tumor cell migration. In contrast to the first mechanism, macroH2A1.1 specifically associates with the transcription start site of Pol II-paused genes. These processes occur in a predefined local 3D genome landscape, but do not require rewiring of enhancer-promoter contacts. We thus propose that macroH2A1.1 serves as a transcriptional modulator with a potential role in assisting the conversion of promoter-locked Pol II into a productive, elongating Pol II.

Chromatin-associated MRN complex protects highly transcribing genes from genomic instability.

MRN-MDC1 plays a central role in the DNA damage response (DDR) and repair. Using proteomics of isolated chromatin fragments, we identified DDR factors, such as MDC1, among those highly associating with a genomic locus upon transcriptional activation. Purification of MDC1 in the absence of exogenous DNA damage revealed interactions with factors involved in gene expression and RNA processing, in addition to DDR factors. ChIP-seq showed that MRN subunits, MRE11 and NBS1, colocalized throughout the genome, notably at TSSs and bodies of actively transcribing genes, which was dependent on the RNAPII transcriptional complex rather than transcription per se. Depletion of MRN increased RNAPII abundance at MRE11/NBS1-bound genes. Prolonged MRE11 or NBS1 depletion induced single-nucleotide polymorphisms across actively transcribing MRN target genes. These data suggest that association of MRN with the transcriptional machinery constitutively scans active genes for transcription-induced DNA damage to preserve the integrity of the coding genome.

Exploring histone loading on HIV DNA reveals a dynamic nucleosome positioning between unintegrated and integrated viral genome.

The aim of the present study was to understand the biology of unintegrated HIV-1 DNA and reveal the mechanisms involved in its transcriptional silencing. We found that histones are loaded on HIV-1 DNA after its nuclear import and before its integration in the host genome. Nucleosome positioning analysis along the unintegrated and integrated viral genomes revealed major differences in nucleosome density and position. Indeed, in addition to the well-known nucleosomes Nuc0, Nuc1, and Nuc2 loaded on integrated HIV-1 DNA, we also found NucDHS, a nucleosome that covers the DNase hypersensitive site, in unintegrated viral DNA. In addition, unintegrated viral DNA-associated Nuc0 and Nuc2 were positioned slightly more to the 5' end relative to their position in integrated DNA. The presence of NucDHS in the proximal region of the long terminal repeat (LTR) promoter was associated with the absence of RNAPII and of the active histone marks H3K4me3 and H3ac at the LTR. Conversely, analysis of integrated HIV-1 DNA showed a loss of NucDHS, loading of RNAPII, and enrichment in active histone marks within the LTR. We propose that unintegrated HIV-1 DNA adopts a repressive chromatin structure that competes with the transcription machinery, leading to its silencing.

Dense neural networks for predicting chromatin conformation.

BACKGROUND: DNA inside eukaryotic cells wraps around histones to form the 11nm chromatin fiber that can further fold into higher-order DNA loops, which may depend on the binding of architectural factors. Predicting how the DNA will fold given a distribution of bound factors, here viewed as a type of sequence, is currently an unsolved problem and several heterogeneous polymer models have shown that many features of the measured structure can be reproduced from simulations. However a model that determines the optimal connection between sequence and structure and that can rapidly assess the effects of varying either one is still lacking. RESULTS: Here we train a dense neural network to solve for the local folding of chromatin, connecting structure, represented as a contact map, to a sequence of bound chromatin factors. The network includes a convolutional filter that compresses the large number of bound chromatin factors into a single 1D sequence representation that is optimized for predicting structure. We also train a network to solve the inverse problem, namely given only structural information in the form of a contact map, predict the likely sequence of chromatin states that generated it. CONCLUSIONS: By carrying out sensitivity analysis on both networks, we are able to highlight the importance of chromatin contexts and neighborhoods for regulating long-range contacts, along with critical alterations that affect contact formation. Our analysis shows that the networks have learned physical insights that are informative and intuitive about this complex polymer problem.

Predicting double-strand DNA breaks using epigenome marks or DNA at kilobase resolution.

Double-strand breaks (DSBs) result from the attack of both DNA strands by multiple sources, including radiation and chemicals. DSBs can cause the abnormal chromosomal rearrangements associated with cancer. Recent techniques allow the genome-wide mapping of DSBs at high resolution, enabling the comprehensive study of their origins. However, these techniques are costly and challenging. Hence, we devise a computational approach to predict DSBs using the epigenomic and chromatin context, for which public data are readily available from the ENCODE project. We achieve excellent prediction accuracy at high resolution. We identify chromatin accessibility, activity, and long-range contacts as the best predictors.

TAD-free analysis of architectural proteins and insulators.

The three-dimensional (3D) organization of the genome is intimately related to numerous key biological functions including gene expression and DNA replication regulations. The mechanisms by which molecular drivers functionally organize the 3D genome, such as topologically associating domains (TADs), remain to be explored. Current approaches consist in assessing the enrichments or influences of proteins at TAD borders. Here, we propose a TAD-free model to directly estimate the blocking effects of architectural proteins, insulators and DNA motifs on long-range contacts, making the model intuitive and biologically meaningful. In addition, the model allows analyzing the whole Hi-C information content (2D information) instead of only focusing on TAD borders (1D information). The model outperforms multiple logistic regression at TAD borders in terms of parameter estimation accuracy and is validated by enhancer-blocking assays. In Drosophila, the results support the insulating role of simple sequence repeats and suggest that the blocking effects depend on the number of repeats. Motif analysis uncovered the roles of the transcriptional factors pannier and tramtrack in blocking long-range contacts. In human, the results suggest that the blocking effects of the well-known architectural proteins CTCF, cohesin and ZNF143 depend on the distance between loci, where each protein may participate at different scales of the 3D chromatin organization.

Dynamic chromatin technologies: from individual molecules to epigenomic regulation in cells.

The establishment and maintenance of chromatin states involves multiscale dynamic processes integrating transcription factor and multiprotein effector dynamics, cycles of chemical chromatin modifications, and chromatin structural organization. Recent developments in genomic technologies are emerging that are enabling a view beyond ensemble- and time-averaged properties and are revealing the importance of dynamic chromatin states for cell fate decisions, differentiation and reprogramming at the single-cell level. Concurrently, biochemical and single-molecule methodologies are providing key insights into the underlying molecular mechanisms. Combining results from defined in vitro and single-molecule studies with single-cell genomic approaches thus holds great promise for understanding chromatin-based transcriptional memory and cell fate. In this Review, we discuss recent developments in biochemical, single-molecule biophysical and single-cell genomic technologies and review how the findings from these approaches can be integrated to paint a comprehensive picture of dynamic chromatin states.

Uncovering direct and indirect molecular determinants of chromatin loops using a computational integrative approach.

Chromosomal organization in 3D plays a central role in regulating cell-type specific transcriptional and DNA replication timing programs. Yet it remains unclear to what extent the resulting long-range contacts depend on specific molecular drivers. Here we propose a model that comprehensively assesses the influence on contacts of DNA-binding proteins, cis-regulatory elements and DNA consensus motifs. Using real data, we validate a large number of predictions for long-range contacts involving known architectural proteins and DNA motifs. Our model outperforms existing approaches including enrichment test, random forests and correlation, and it uncovers numerous novel long-range contacts in Drosophila and human. The model uncovers the orientation-dependent specificity for long-range contacts between CTCF motifs in Drosophila, highlighting its conserved property in 3D organization of metazoan genomes. Our model further unravels long-range contacts depending on co-factors recruited to DNA indirectly, as illustrated by the influence of cohesin in stabilizing long-range contacts between CTCF sites. It also reveals asymmetric contacts such as enhancer-promoter contacts that highlight opposite influences of the transcription factors EBF1, EGR1 or MEF2C depending on RNA Polymerase II pausing.

Computational Identification of Genomic Features That Influence 3D Chromatin Domain Formation.

Recent advances in long-range Hi-C contact mapping have revealed the importance of the 3D structure of chromosomes in gene expression. A current challenge is to identify the key molecular drivers of this 3D structure. Several genomic features, such as architectural proteins and functional elements, were shown to be enriched at topological domain borders using classical enrichment tests. Here we propose multiple logistic regression to identify those genomic features that positively or negatively influence domain border establishment or maintenance. The model is flexible, and can account for statistical interactions among multiple genomic features. Using both simulated and real data, we show that our model outperforms enrichment test and non-parametric models, such as random forests, for the identification of genomic features that influence domain borders. Using Drosophila Hi-C data at a very high resolution of 1 kb, our model suggests that, among architectural proteins, BEAF-32 and CP190 are the main positive drivers of 3D domain borders. In humans, our model identifies well-known architectural proteins CTCF and cohesin, as well as ZNF143 and Polycomb group proteins as positive drivers of domain borders. The model also reveals the existence of several negative drivers that counteract the presence of domain borders including P300, RXRA, BCL11A and ELK1.

Predicting the spatial organization of chromosomes using epigenetic data.

Chromosome folding can reinforce the demarcation between euchromatin and heterochromatin. Two new studies show how epigenetic data, including DNA methylation, can accurately predict chromosome folding in three dimensions. Such computational approaches reinforce the idea of a linkage between epigenetically marked chromatin domains and their segregation into distinct compartments at the megabase scale or topological domains at a higher resolution. Please see related articles: http://dx.doi.org/10.1186/s13059-015-0741-y and http://dx.doi.org/10.1186/s13059-015-0740-z.

Probing long-range interactions by extracting free energies from genome-wide chromosome conformation capture data.

BACKGROUND: A variety of DNA binding proteins are involved in regulating and shaping the packing of chromatin. They aid the formation of loops in the DNA that function to isolate different structural domains. A recent experimental technique, Hi-C, provides a method for determining the frequency of such looping between all distant parts of the genome. Given that the binding locations of many chromatin associated proteins have also been measured, it has been possible to make estimates for their influence on the long-range interactions as measured by Hi-C. However, a challenge in this analysis is the predominance of non-specific contacts that mask out the specific interactions of interest. RESULTS: We show that transforming the Hi-C contact frequencies into free energies gives a natural method for separating out the distance dependent non-specific interactions. In particular we apply Principal Component Analysis (PCA) to the transformed free energy matrix to identify the dominant modes of interaction. PCA identifies systematic effects as well as high frequency spatial noise in the Hi-C data which can be filtered out. Thus it can be used as a data driven approach for normalizing Hi-C data. We assess this PCA based normalization approach, along with several other normalization schemes, by fitting the transformed Hi-C data using a pairwise interaction model that takes as input the known locations of bound chromatin factors. The result of fitting is a set of predictions for the coupling energies between the various chromatin factors and their effect on the energetics of looping. We show that the quality of the fit can be used as a means to determine how much PCA filtering should be applied to the Hi-C data. CONCLUSIONS: We find that the different normalizations of the Hi-C data vary in the quality of fit to the pairwise interaction model. PCA filtering can improve the fit, and the predicted coupling energies lead to biologically meaningful insights for how various chromatin bound factors influence the stability of DNA loops in chromatin.

Integrator complex regulates NELF-mediated RNA polymerase II pause/release and processivity at coding genes.

RNA polymerase II (RNAPII) pausing/termination shortly after initiation is a hallmark of gene regulation. Here, we show that negative elongation factor (NELF) interacts with Integrator complex subunits (INTScom), RNAPII and Spt5. The interaction between NELF and INTScom subunits is RNA and DNA independent. Using both human immunodeficiency virus type 1 promoter and genome-wide analyses, we demonstrate that Integrator subunits specifically control NELF-mediated RNAPII pause/release at coding genes. The strength of RNAPII pausing is determined by the nature of the NELF-associated INTScom subunits. Interestingly, in addition to controlling RNAPII pause-release INTS11 catalytic subunit of the INTScom is required for RNAPII processivity. Finally, INTScom target genes are enriched in human immunodeficiency virus type 1 transactivation response element/NELF binding element and in a 3' box sequence required for small nuclear RNA biogenesis. Revealing these unexpected functions of INTScom in regulating RNAPII pause-release and completion of mRNA synthesis of NELF-target genes will contribute to our understanding of the gene expression cycle.

Chromatin insulator factors involved in long-range DNA interactions and their role in the folding of the Drosophila genome.

Chromatin insulators are genetic elements implicated in the organization of chromatin and the regulation of transcription. In Drosophila, different insulator types were characterized by their locus-specific composition of insulator proteins and co-factors. Insulators mediate specific long-range DNA contacts required for the three dimensional organization of the interphase nucleus and for transcription regulation, but the mechanisms underlying the formation of these contacts is currently unknown. Here, we investigate the molecular associations between different components of insulator complexes (BEAF32, CP190 and Chromator) by biochemical and biophysical means, and develop a novel single-molecule assay to determine what factors are necessary and essential for the formation of long-range DNA interactions. We show that BEAF32 is able to bind DNA specifically and with high affinity, but not to bridge long-range interactions (LRI). In contrast, we show that CP190 and Chromator are able to mediate LRI between specifically-bound BEAF32 nucleoprotein complexes in vitro. This ability of CP190 and Chromator to establish LRI requires specific contacts between BEAF32 and their C-terminal domains, and dimerization through their N-terminal domains. In particular, the BTB/POZ domains of CP190 form a strict homodimer, and its C-terminal domain interacts with several insulator binding proteins. We propose a general model for insulator function in which BEAF32/dCTCF/Su(HW) provide DNA specificity (first layer proteins) whereas CP190/Chromator are responsible for the physical interactions required for long-range contacts (second layer). This network of organized, multi-layer interactions could explain the different activities of insulators as chromatin barriers, enhancer blockers, and transcriptional regulators, and suggest a general mechanism for how insulators may shape the organization of higher-order chromatin during cell division.

Insulators recruit histone methyltransferase dMes4 to regulate chromatin of flanking genes.

Chromosomal domains in Drosophila are marked by the insulator-binding proteins (IBPs) dCTCF/Beaf32 and cofactors that participate in regulating long-range interactions. Chromosomal borders are further enriched in specific histone modifications, yet the role of histone modifiers and nucleosome dynamics in this context remains largely unknown. Here, we show that IBP depletion impairs nucleosome dynamics specifically at the promoters and coding sequence of genes flanked by IBP binding sites. Biochemical purification identifies the H3K36 histone methyltransferase NSD/dMes-4 as a novel IBP cofactor, which specifically co-regulates the chromatin accessibility of hundreds of genes flanked by dCTCF/Beaf32. NSD/dMes-4 presets chromatin before the recruitment of transcriptional activators including DREF that triggers Set2/Hypb-dependent H3K36 trimethylation, nucleosome positioning, and RNA splicing. Our results unveil a model for how IBPs regulate nucleosome dynamics and gene expression through NSD/dMes-4, which may regulate H3K27me3 spreading. Our data uncover how IBPs dynamically regulate chromatin organization depending on distinct cofactors.