Insulator function and topological domain border strength scale with architectural protein occupancy

Background

Chromosome conformation capture studies suggest that eukaryotic genomes are organized into structures called topologically associating domains. The borders of these domains are highly enriched for architectural proteins with characterized roles in insulator function. However, a majority of architectural protein binding sites localize within topological domains, suggesting sites associated with domain borders represent a functionally different subclass of these regulatory elements. How topologically associating domains are established and what differentiates border-associated from non-border architectural protein binding sites remain unanswered questions.

Results

By mapping the genome-wide target sites for several Drosophila architectural proteins, including previously uncharacterized profiles for TFIIIC and SMC-containing condensin complexes, we uncover an extensive pattern of colocalization in which architectural proteins establish dense clusters at the borders of topological domains. Reporter-based enhancer-blocking insulator activity as well as endogenous domain border strength scale with the occupancy level of architectural protein binding sites, suggesting co-binding by architectural proteins underlies the functional potential of these loci. Analyses in mouse and human stem cells suggest that clustering of architectural proteins is a general feature of genome organization, and conserved architectural protein binding sites may underlie the tissue-invariant nature of topologically associating domains observed in mammals.

Conclusions

We identify a spectrum of architectural protein occupancy that scales with the topological structure of chromosomes and the regulatory potential of these elements. Whereas high occupancy architectural protein binding sites associate with robust partitioning of topologically associating domains and robust insulator function, low occupancy sites appear reserved for gene-specific regulation within topological domains.     

A Comprehensive Map of Insulator Elements for the Drosophila Genome

Insulators are DNA sequences that control the interactions among genomic regulatory elements and act as chromatin boundaries. A thorough understanding of their location and function is necessary to address the complexities of metazoan gene regulation. We studied by ChIP–chip the genome-wide binding sites of 6 insulator-associated proteins—dCTCF, CP190, BEAF-32, Su(Hw), Mod(mdg4), and GAF—to obtain the first comprehensive map of insulator elements in Drosophila embryos. We identify over 14,000 putative insulators, including all classically defined insulators. We find two major classes of insulators defined by dCTCF/CP190/BEAF-32 and Su(Hw), respectively. Distributional analyses of insulators revealed that particular sub-classes of insulator elements are excluded between cis-regulatory elements and their target promoters; divide differentially expressed, alternative, and divergent promoters; act as chromatin boundaries; are associated with chromosomal breakpoints among species; and are embedded within active chromatin domains. Together, these results provide a map demarcating the boundaries of gene regulatory units and a framework for understanding insulator function during the development and evolution of Drosophila.     

Functional validation of mouse tyrosinase non-coding regulatory DNA elements by CRISPR–Cas9-mediated mutagenesis

Newly developed genome-editing tools, such as the clustered regularly interspaced short palindromic repeat (CRISPR)–Cas9 system, allow simple and rapid genetic modification in most model organisms and human cell lines. Here, we report the production and analysis of mice carrying the inactivation via deletion of a genomic insulator, a key non-coding regulatory DNA element found 5′ upstream of the mouse tyrosinase (Tyr) gene. Targeting sequences flanking this boundary in mouse fertilized eggs resulted in the efficient deletion or inversion of large intervening DNA fragments delineated by the RNA guides. The resulting genome-edited mice showed a dramatic decrease in Tyr gene expression as inferred from the evident decrease of coat pigmentation, thus supporting the functionality of this boundary sequence in vivo, at the endogenous locus. Several potential off-targets bearing sequence similarity with each of the two RNA guides used were analyzed and found to be largely intact. This study reports how non-coding DNA elements, even if located in repeat-rich genomic sequences, can be efficiently and functionally evaluated in vivo and, furthermore, it illustrates how the regulatory elements described by the ENCODE and EPIGENOME projects, in the mouse and human genomes, can be systematically validated.     

Boundaries of transcriptionally silent chromatin in Saccharomyces cerevisiae

In the budding yeast Saccharomyces cerevisiae, heterochromatic gene silencing has been found within HMR and HML silent mating type loci, the telomeres, and the rRNA-encoding DNA. There may be boundary elements that regulate the spread of silencing to protect genes adjacent to silenced domains from this epigenetic repressive effect. Many assays show that specific DNA regulatory elements separate a euchromatic locus from a neighboring heterochromatic domain and thereby function as a boundary. Alternatively, DNA-independent mechanisms such as competition between acetylated and deacetylated histones are also reported to contribute to gene insulation. However, the mechanism by which boundaries are formed is not clear. Here, the characteristics and functions of boundaries at silenced domains in S. cerevisiae are discussed.     

Human and mouse LSP1 genes code for highly conserved phosphoproteins

With use of the mouse LSP1 cDNA we isolated a human homologue of the mouse LSP1 gene from a human CTL cDNA library. The predicted protein sequence of human LSP1 is compared with the predicted mouse LSP1 protein sequence and regions of homology are identified in order to predict structural features of the LSP1 protein that might be important for its function. Both the human and mouse LSP1 proteins consist of two domains, an N-terminal acidic domain and a C-terminal basic domain. The C-terminal domains of the mouse and human LSP1 proteins are highly conserved and include several conserved, putative serine/threonine phosphorylation sites. Immunoprecipitation of LSP1 protein from 32P-orthophosphate-loaded cells show that both the mouse and human LSP1 proteins are phosphoproteins. The sequences of the putative Ca2(+)-binding sites present in the N-terminal domain of the mouse LSP1 protein are not conserved in the human LSP1 protein; however, a different Ca2(+)-binding site may exist in the human protein, indicating a functional conservation rather than a strict sequence conservation of the two proteins. The expression of the human LSP1 gene follows the same pattern as the expression of the mouse LSP1 gene. Southern analysis of human genomic DNA shows multiple LSP1-related fragments of varying intensity in contrast to the simple pattern found after similar analysis of mouse genomic DNA. By using different parts of the human LSP1 cDNA as a probe, we show that most of these multiple bands contain sequences homologous to the conserved C-terminal region of the LSP1 cDNA. This suggests that there are several LSP1-related genes present in the human genome.     

Mouse Hus1, a homolog of the Schizosaccharomyces pombe hus1+ cell cycle checkpoint gene.

Cell cycle checkpoints are regulatory mechanisms that arrest the cell cycle or initiate programmed cell death when critical events such as DNA replication fail to be completed or when DNA or spindle damage occurs. In fission yeast, cell cycle checkpoint responses to DNA replication blocks and DNA damage require the hus1+ gene. Mammalian homologs of hus1+ were recently identified, and here we report a detailed analysis of mouse Hus1. An approximately 4.2-kb full-length cDNA encoding the 32-kDa mouse Hus1 protein was isolated. The genomic structure and exon-intron boundary sequences of the gene were determined, and mouse Hus1 was found to consist of nine exons. Mouse Hus1 was mapped to the proximal end of chromosome 11 and is therefore a candidate gene for the mouse mutation germ cell deficient, which maps to the same genomic region. Finally, mouse Hus1 was found to be expressed in a variety of adult tissues and at several stages of embryonic development.     

Dynamic and facilitated binding of topoisomerase accelerates topological relaxation

How type 2 Topoisomerase (TopoII) proteins relax and simplify the topology of DNA molecules is one of the most intriguing open questions in genome and DNA biophysics. Most of the existing models neglect the dynamics of TopoII which is expected of proteins searching their targets via facilitated diffusion. Here, we show that dynamic binding of TopoII speeds up the topological relaxation of knotted substrates by enhancing the search of the knotted arc. Intriguingly, this in turn implies that the timescale of topological relaxation is virtually independent of the substrate length. We then discover that considering binding biases due to facilitated diffusion on looped substrates steers the sampling of the topological space closer to the boundaries between different topoisomers yielding an optimally fast topological relaxation. We discuss our findings in the context of topological simplification in vitro and in vivo.     

DNA loop domains in a 1.4-Mb region around the human hprt gene mapped by cleavage mediated by nuclear matrix-associated topoisomerase II

We have mapped the positions in a ～1.4-Mb region of genomic DNA around the human hprt gene which are accessible in vivo to cleavage by topoisomerase II associated with the nuclear matrix. These positions, which are interpreted as the boundaries of DNA loop domains, were mapped in K562 cells by examining the truncation of rare-cutter restriction fragments separated by pulsed field gel electrophoresis after topoisomerase II-mediated cleavage, using seven linked markers mapped in this region as probes for indirect end-labeling. Eleven cleavage positions were detected and were interpreted as defining ten loop domains of lengths between 70 and 210?kb (average ～135?kb); the hprt gene resides in a 150-kb loop domain. Loop domain boundaries coincided with three of the fifteen deletion breakpoints mapped in a 600-kb sector of this region in human lymphocytes, within the limits of resolution of pulsed field gel electrophoresis; this correlation was not statistically significant.     

Functional analysis of the homeodomain protein SIX5

SIX5 (previously known as myotonic dystrophy associated homeodomain protein-DMAHP) is a member of the SIX [sine oculis homeobox (Drosophila) homologue] gene family which encodes proteins containing a SIX domain adjacent to a homeodomain. To investigate the DNA binding specificities of these two domains in SIX5, they were expressed as GST fusion proteins, both separately and together. Affinity purified recombinant proteins and cell lysates from bacteria expressing the recombinant proteins were used in gel retardation assays with double stranded oligonucleotides representing putative DNA binding sites. The putative sites included two in the promoter region of DMPK (dystrophia myotonica protein kinase) and the previously characterised murine Six4 DNA binding site in the Na⁺/K⁺ ATPase α1 subunit gene (ATP1A1) regulatory element (ARE). None of the recombinant proteins showed any affinity for the two putative sites in DMPK. However, the two recombinant proteins containing the homeodomain both formed at least one specific complex with the ARE. The recombinant protein containing both domains formed a second specific complex with the ARE, assumed to be a dimer complex. Finally, a whole genome PCR-based screen was used to identify genomic DNA sequences to which SIX5 binds, as an initial stage in the identification of genes regulated by SIX5.     

Effects of core histone tail domains on the equilibrium constants for dynamic DNA site accessibility in nucleosomes

The N and C-terminal tail domains of the core histones play important roles in gene regulation, but the mechanisms through which they act are not known. These tail domains are highly positively charged and are the sites of numerous post-translational modifications, including many sites for lysine acetylation. Nucleosomes in which these tail domains have been removed by trypsin remain otherwise intact, and are used by many laboratories as a model system for highly acetylated nucleosomes. Here, we test the hypothesis that one role of the tail domains is to directly regulate the accessibility of nucleosomal DNA to other DNA-binding proteins. Three assays are used: equilibrium binding by a site-specific, DNA-binding protein, and dynamic accessibility to restriction enzymes or to a non-specific exonuclease. The effects of removal of the tail domains as monitored by each of these assays can be understood within the framework of the site exposure model for the dynamic equilibrium accessibility of target sites located within the nucleosomal DNA. Removal of the tail domains leads to a 1.5 to 14-fold increase in position-dependent equilibrium constants for site exposure. The smallness of the effect weighs against models for gene activation in which histone acetylation is a mandatory initial event, required to facilitate subsequent access of regulatory proteins to nucleosomal DNA target sites. Alternative roles for histone acetylation in gene regulation are discussed.