High-throughput identification of human SNPs affecting regulatory element activity.

Most of the millions of SNPs in the human genome are non-coding, and many overlap with putative regulatory elements. Genome-wide association studies (GWAS) have linked many of these SNPs to human traits or to gene expression levels, but rarely with sufficient resolution to identify the causal SNPs. Functional screens based on reporter assays have previously been of insufficient throughput to test the vast space of SNPs for possible effects on regulatory element activity. Here we leveraged the throughput and resolution of the survey of regulatory elements (SuRE) reporter technology to survey the effect of 5.9 million SNPs, including 57% of the known common SNPs, on enhancer and promoter activity. We identified more than 30,000 SNPs that alter the activity of putative regulatory elements, partially in a cell-type-specific manner. Integration of this dataset with GWAS results may help to pinpoint SNPs that underlie human traits.

Promoter-Intrinsic and Local Chromatin Features Determine Gene Repression in LADs.

It is largely unclear whether genes that are naturally embedded in lamina-associated domains (LADs) are inactive due to their chromatin environment or whether LADs are merely secondary to the lack of transcription. We show that hundreds of human promoters become active when moved from their native LAD position to a neutral context in the same cells, indicating that LADs form a repressive environment. Another set of promoters inside LADs is able to "escape" repression, although their transcription elongation is attenuated. By inserting reporters into thousands of genomic locations, we demonstrate that escaper promoters are intrinsically less sensitive to LAD repression. This is not simply explained by promoter strength but by the interplay between promoter sequence and local chromatin features that vary strongly across LADs. Enhancers also differ in their sensitivity to LAD chromatin. This work provides a general framework for the systematic understanding of gene regulation by repressive chromatin.

The large fraction of heterochromatin in Drosophila neurons is bound by both B-type lamin and HP1a.

BACKGROUND: In most mammalian cell lines, chromatin located at the nuclear periphery is represented by condensed heterochromatin, as evidenced by microscopy observations and DamID mapping of lamina-associated domains (LADs) enriched in dimethylated Lys9 of histone H3 (H3K9me2). However, in Kc167 cell culture, the only Drosophilla cell type where LADs have previously been mapped, they are neither H3K9me2-enriched nor overlapped with the domains of heterochromatin protein 1a (HP1a). RESULTS: Here, using cell type-specific DamID we mapped genome-wide LADs, HP1a and Polycomb (Pc) domains from the central brain, Repo-positive glia, Elav-positive neurons and the fat body of Drosophila third instar larvae. Strikingly, contrary to Kc167 cells of embryonic origin, in neurons and, to a lesser extent, in glia and the fat body, HP1a domains appear to overlap strongly with LADs in both the chromosome arms and pericentromeric regions. Accordingly, centromeres reside closer to the nuclear lamina in neurons than in Kc167 cells. As expected, active gene promoters are mostly not present in LADs, HP1a and Pc domains. These domains are occupied by silent or weakly expressed genes with genes residing in the HP1a-bound LADs expressed at the lowest level. CONCLUSIONS: In various differentiated Drosophila cell types, we discovered the existence of peripheral heterochromatin, similar to that observed in mammals. Our findings support the model that peripheral heterochromatin matures enhancing the repression of unwanted genes as cells terminally differentiate.

Massive reshaping of genome-nuclear lamina interactions during oncogene-induced senescence.

Cellular senescence is a mechanism that virtually irreversibly suppresses the proliferative capacity of cells in response to various stress signals. This includes the expression of activated oncogenes, which causes Oncogene-Induced Senescence (OIS). A body of evidence points to the involvement in OIS of chromatin reorganization, including the formation of senescence-associated heterochromatic foci (SAHF). The nuclear lamina (NL) is an important contributor to genome organization and has been implicated in cellular senescence and organismal aging. It interacts with multiple regions of the genome called lamina-associated domains (LADs). Some LADs are cell-type specific, whereas others are conserved between cell types and are referred to as constitutive LADs (cLADs). Here, we used DamID to investigate the changes in genome-NL interactions in a model of OIS triggered by the expression of the common BRAFV600E oncogene. We found that OIS cells lose most of their cLADS, suggesting the loss of a specific mechanism that targets cLADs to the NL. In addition, multiple genes relocated to the NL. Unexpectedly, they were not repressed, implying the abrogation of the repressive activity of the NL during OIS. Finally, OIS cells displayed an increased association of telomeres with the NL. Our study reveals that senescent cells acquire a new type of LAD organization and suggests the existence of as yet unknown mechanisms that tether cLADs to the NL and repress gene expression at the NL.

Small chromosomal regions position themselves autonomously according to their chromatin class.

The spatial arrangement of chromatin is linked to the regulation of nuclear processes. One striking aspect of nuclear organization is the spatial segregation of heterochromatic and euchromatic domains. The mechanisms of this chromatin segregation are still poorly understood. In this work, we investigated the link between the primary genomic sequence and chromatin domains. We analyzed the spatial intranuclear arrangement of a human artificial chromosome (HAC) in a xenospecific mouse background in comparison to an orthologous region of native mouse chromosome. The two orthologous regions include segments that can be assigned to three major chromatin classes according to their gene abundance and repeat repertoire: (1) gene-rich and SINE-rich euchromatin; (2) gene-poor and LINE/LTR-rich heterochromatin; and (3) gene-depleted and satellite DNA-containing constitutive heterochromatin. We show, using fluorescence in situ hybridization (FISH) and 4C-seq technologies, that chromatin segments ranging from 0.6 to 3 Mb cluster with segments of the same chromatin class. As a consequence, the chromatin segments acquire corresponding positions in the nucleus irrespective of their chromosomal context, thereby strongly suggesting that this is their autonomous property. Interactions with the nuclear lamina, although largely retained in the HAC, reveal less autonomy. Taken together, our results suggest that building of a functional nucleus is largely a self-organizing process based on mutual recognition of chromosome segments belonging to the major chromatin classes.

Genome-wide mapping of autonomous promoter activity in human cells.

Previous methods to systematically characterize sequence-intrinsic activity of promoters have been limited by relatively low throughput and the length of the sequences that could be tested. Here we present 'survey of regulatory elements' (SuRE), a method that assays more than 108 DNA fragments, each 0.2-2 kb in size, for their ability to drive transcription autonomously. In SuRE, a plasmid library of random genomic fragments upstream of a 20-bp barcode is constructed, and decoded by paired-end sequencing. This library is used to transfect cells, and barcodes in transcribed RNA are quantified by high-throughput sequencing. When applied to the human genome, we achieve 55-fold genome coverage, allowing us to map autonomous promoter activity genome-wide in K562 cells. By computational modeling we delineate subregions within promoters that are relevant for their activity. We show that antisense promoter transcription is generally dependent on the sense core promoter sequences, and that most enhancers and several families of repetitive elements act as autonomous transcription initiation sites.

High-throughput assessment of context-dependent effects of chromatin proteins.

BACKGROUND: Chromatin proteins control gene activity in a concerted manner. We developed a high-throughput assay to study the effects of the local chromatin environment on the regulatory activity of a protein of interest. The assay combines a previously reported multiplexing strategy based on barcoded randomly integrated reporters with Gal4-mediated tethering. We applied the assay to Drosophila heterochromatin protein 1a (HP1a), which is mostly known as a repressive protein but has also been linked to transcriptional activation. RESULTS: Recruitment to over 1000 genomic locations revealed that HP1a is a potent repressor able to silence even highly expressing reporter genes. However, the local chromatin context can modulate HP1a function. In pericentromeric regions, HP1a-induced repression was enhanced by twofold. In regions marked by a H3K36me3-rich chromatin signature, HP1a-dependent silencing was significantly decreased. We found no evidence for an activating function of HP1a in our experimental system. Furthermore, we did not observe stable transmission of repression over mitotic divisions after loss of targeted HP1a. CONCLUSIONS: The multiplexed tethered reporter assay should be applicable to a large number of chromatin proteins and will be a useful tool to dissect combinatorial regulatory interactions in chromatin.

Inducible DamID systems for genomic mapping of chromatin proteins in Drosophila.

Dam identification (DamID) is a powerful technique to generate genome-wide maps of chromatin protein binding. Due to its high sensitivity, it is particularly suited to study the genome interactions of chromatin proteins in small tissue samples in model organisms such as Drosophila Here, we report an intein-based approach to tune the expression level of Dam and Dam-fusion proteins in Drosophila by addition of a ligand to fly food. This helps to suppress possible toxic effects of Dam. In addition, we describe a strategy for genetically controlled expression of Dam in a specific cell type in complex tissues. We demonstrate the utility of the latter by generating a glia-specific map of Polycomb in small samples of brain tissue. These new DamID tools will be valuable for the mapping of binding patterns of chromatin proteins in Drosophila tissues and especially in cell lineages.

Genome-wide maps of nuclear lamina interactions in single human cells.

Mammalian interphase chromosomes interact with the nuclear lamina (NL) through hundreds of large lamina-associated domains (LADs). We report a method to map NL contacts genome-wide in single human cells. Analysis of nearly 400 maps reveals a core architecture consisting of gene-poor LADs that contact the NL with high cell-to-cell consistency, interspersed by LADs with more variable NL interactions. The variable contacts tend to be cell-type specific and are more sensitive to changes in genome ploidy than the consistent contacts. Single-cell maps indicate that NL contacts involve multivalent interactions over hundreds of kilobases. Moreover, we observe extensive intra-chromosomal coordination of NL contacts, even over tens of megabases. Such coordinated loci exhibit preferential interactions as detected by Hi-C. Finally, the consistency of NL contacts is inversely linked to gene activity in single cells and correlates positively with the heterochromatic histone modification H3K9me3. These results highlight fundamental principles of single-cell chromatin organization. VIDEO ABSTRACT.

Using TRIP for genome-wide position effect analysis in cultured cells.

The influence of local chromatin context on gene expression can be explored by integrating a transcription reporter at different locations in the genome as a sensor. Here we provide a detailed protocol for analyzing thousands of reporters integrated in parallel (TRIP) at a genome-wide level. TRIP is based on tagging each reporter with a unique barcode, which is used for independent reporter expression analysis and integration site mapping. Compared with previous methods for studying position effects, TRIP offers a 100-1,000-fold higher throughput in a faster and less-labor-intensive manner. The entire experimental protocol takes ∼42 d to complete, with high-throughput sequencing and data analysis requiring an additional ∼11 d. TRIP was developed by using transcription reporters in mouse embryonic stem (mES) cells, but because of its flexibility the method can be used to probe the influence of chromatin context on a variety of molecular processes in any transfectable cell line.

Chromatin position effects assayed by thousands of reporters integrated in parallel.

Reporter genes integrated into the genome are a powerful tool to reveal effects of regulatory elements and local chromatin context on gene expression. However, so far such reporter assays have been of low throughput. Here, we describe a multiplexing approach for the parallel monitoring of transcriptional activity of thousands of randomly integrated reporters. More than 27,000 distinct reporter integrations in mouse embryonic stem cells, obtained with two different promoters, show ∼1,000-fold variation in expression levels. Data analysis indicates that lamina-associated domains act as attenuators of transcription, likely by reducing access of transcription factors to binding sites. Furthermore, chromatin compaction is predictive of reporter activity. We also found evidence for crosstalk between neighboring genes and estimate that enhancers can influence gene expression on average over ∼20 kb. The multiplexed reporter assay is highly flexible in design and can be modified to query a wide range of aspects of gene regulation.

Single-cell dynamics of genome-nuclear lamina interactions.

The nuclear lamina (NL) interacts with hundreds of large genomic regions termed lamina associated domains (LADs). The dynamics of these interactions and the relation to epigenetic modifications are poorly understood. We visualized the fate of LADs in single cells using a "molecular contact memory" approach. In each nucleus, only ~30% of LADs are positioned at the periphery; these LADs are in intermittent molecular contact with the NL but remain constrained to the periphery. Upon mitosis, LAD positioning is not detectably inherited but instead is stochastically reshuffled. Contact of individual LADs with the NL is linked to transcriptional repression and H3K9 dimethylation in single cells. Furthermore, we identify the H3K9 methyltransferase G9a as a regulator of NL contacts. Collectively, these results highlight principles of the dynamic spatial architecture of chromosomes in relation to gene regulation.

Constitutive nuclear lamina-genome interactions are highly conserved and associated with A/T-rich sequence.

In metazoans, the nuclear lamina is thought to play an important role in the spatial organization of interphase chromosomes, by providing anchoring sites for large genomic segments named lamina-associated domains (LADs). Some of these LADs are cell-type specific, while many others appear constitutively associated with the lamina. Constitutive LADs (cLADs) may contribute to a basal chromosome architecture. By comparison of mouse and human lamina interaction maps, we find that the sizes and genomic positions of cLADs are strongly conserved. Moreover, cLADs are depleted of synteny breakpoints, pointing to evolutionary selective pressure to keep cLADs intact. Paradoxically, the overall sequence conservation is low for cLADs. Instead, cLADs are universally characterized by long stretches of DNA of high A/T content. Cell-type specific LADs also tend to adhere to this "A/T rule" in embryonic stem cells, but not in differentiated cells. This suggests that the A/T rule represents a default positioning mechanism that is locally overruled during lineage commitment. Analysis of paralogs suggests that during evolution changes in A/T content have driven the relocation of genes to and from the nuclear lamina, in tight association with changes in expression level. Taken together, these results reveal that the spatial organization of mammalian genomes is highly conserved and tightly linked to local nucleotide composition.

Glucosylated hydroxymethyluracil, DNA base J, prevents transcriptional readthrough in Leishmania.

Some Ts in nuclear DNA of trypanosomes and Leishmania are hydroxylated and glucosylated to yield base J (ß-D-glucosyl-hydroxymethyluracil). In Leishmania, about 99% of J is located in telomeric repeats. We show here that most of the remaining J is located at chromosome-internal RNA polymerase II termination sites. This internal J and telomeric J can be reduced by a knockout of J-binding protein 2 (JBP2), an enzyme involved in the first step of J biosynthesis. J levels are further reduced by growing Leishmania JBP2 knockout cells in BrdU-containing medium, resulting in cell death. The loss of internal J in JBP2 knockout cells is accompanied by massive readthrough at RNA polymerase II termination sites. The readthrough varies between transcription units but may extend over 100 kb. We conclude that J is required for proper transcription termination and infer that the absence of internal J kills Leishmania by massive readthrough of transcriptional stops.

Interactions among Polycomb domains are guided by chromosome architecture.

Polycomb group (PcG) proteins bind and regulate hundreds of genes. Previous evidence has suggested that long-range chromatin interactions may contribute to the regulation of PcG target genes. Here, we adapted the Chromosome Conformation Capture on Chip (4C) assay to systematically map chromosomal interactions in Drosophila melanogaster larval brain tissue. Our results demonstrate that PcG target genes interact extensively with each other in nuclear space. These interactions are highly specific for PcG target genes, because non-target genes with either low or high expression show distinct interactions. Notably, interactions are mostly limited to genes on the same chromosome arm, and we demonstrate that a topological rather than a sequence-based mechanism is responsible for this constraint. Our results demonstrate that many interactions among PcG target genes exist and that these interactions are guided by overall chromosome architecture.

The insulator protein SU(HW) fine-tunes nuclear lamina interactions of the Drosophila genome.

Specific interactions of the genome with the nuclear lamina (NL) are thought to assist chromosome folding inside the nucleus and to contribute to the regulation of gene expression. High-resolution mapping has recently identified hundreds of large, sharply defined lamina-associated domains (LADs) in the human genome, and suggested that the insulator protein CTCF may help to demarcate these domains. Here, we report the detailed structure of LADs in Drosophila cells, and investigate the putative roles of five insulator proteins in LAD organization. We found that the Drosophila genome is also organized in discrete LADs, which are about five times smaller than human LADs but contain on average a similar number of genes. Systematic comparison to new and published insulator binding maps shows that only SU(HW) binds preferentially at LAD borders and at specific positions inside LADs, while GAF, CTCF, BEAF-32 and DWG are mostly absent from these regions. By knockdown and overexpression studies we demonstrate that SU(HW) weakens genome - NL interactions through a local antagonistic effect, but we did not obtain evidence that it is essential for border formation. Our results provide insights into the evolution of LAD organization and identify SU(HW) as a fine-tuner of genome - NL interactions.

Molecular maps of the reorganization of genome-nuclear lamina interactions during differentiation.

The three-dimensional organization of chromosomes within the nucleus and its dynamics during differentiation are largely unknown. To visualize this process in molecular detail, we generated high-resolution maps of genome-nuclear lamina interactions during subsequent differentiation of mouse embryonic stem cells via lineage-committed neural precursor cells into terminally differentiated astrocytes. This reveals that a basal chromosome architecture present in embryonic stem cells is cumulatively altered at hundreds of sites during lineage commitment and subsequent terminal differentiation. This remodeling involves both individual transcription units and multigene regions and affects many genes that determine cellular identity. Often, genes that move away from the lamina are concomitantly activated; many others, however, remain inactive yet become unlocked for activation in a next differentiation step. These results suggest that lamina-genome interactions are widely involved in the control of gene expression programs during lineage commitment and terminal differentiation.

Histone H1 binding is inhibited by histone variant H3.3.

Linker histones are involved in the formation of higher-order chromatin structure and the regulation of specific genes, yet it remains unclear what their principal binding determinants are. We generated a genome-wide high-resolution binding map for linker histone H1 in Drosophila cells, using DamID. H1 binds at similar levels across much of the genome, both in classic euchromatin and heterochromatin. Strikingly, there are pronounced dips of low H1 occupancy around transcription start sites for active genes and at many distant cis-regulatory sites. H1 dips are not due to lack of nucleosomes; rather, all regions with low binding of H1 show enrichment of the histone variant H3.3. Knockdown of H3.3 causes H1 levels to increase at these sites, with a concomitant increase in nucleosome repeat length. These changes are independent of transcriptional changes. Our results show that the H3.3 protein counteracts association of H1, providing a mechanism to keep diverse genomic sites in an open chromatin conformation.

High-resolution mapping of heterochromatin redistribution in a Drosophila position-effect variegation model.

BACKGROUND: Position-effect variegation (PEV) is the stochastic transcriptional silencing of a gene positioned adjacent to heterochromatin. white-mottled X-chromosomal inversions in Drosophila are classic PEV models that show variegation of the eye color gene white due to its relocation next to pericentric heterochromatin. It has been suggested that in these models the spreading of heterochromatin across the rearrangement breakpoint causes the silencing of white. However, the extent of this spreading and the precise pattern of heterochromatin redistribution have remained unclear. To obtain insight into the mechanism of PEV, we constructed high-resolution binding maps of Heterochromatin Protein 1 (HP1) on white-mottled chromosomes. RESULTS: We find that HP1 invades euchromatin across the inversion breakpoints over approximately 175 kb and approximately 30 kb, causing de novo association of HP1 with 20 genes. However, HP1 binding levels in these regions show substantial local variation, and white is the most strongly bound gene. Remarkably, white is also the only gene that is detectably repressed by heterochromatin. Furthermore, we find that HP1 binding to the invaded region is particularly sensitive to the dosage of the histone methyltransferase Su(var)3-9, indicating that the de novo formed heterochromatin is less stable than naturally occurring constitutive heterochromatin. CONCLUSION: Our molecular maps demonstrate that heterochromatin can invade a normally euchromatic region, yet the strength of HP1 binding and effects on gene expression are highly dependent on local context. Our data suggest that the white gene has an unusual intrinsic affinity for heterochromatin, which may cause this gene to be more sensitive to PEV than most other genes.