KRAB–Zinc Finger Proteins and KAP1 Can Mediate Long-Range Transcriptional Repression through Heterochromatin Spreading

doi:10.1371/journal.pgen.1000869

PLoS Genet. 2010 March; 6(3): e1000869.

Published online 2010 March 5. doi: 10.1371/journal.pgen.1000869

PMCID: PMC2832679

Copyright Groner et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

KRAB–Zinc Finger Proteins and KAP1 Can Mediate Long-Range Transcriptional Repression through Heterochromatin Spreading

Anna C. Groner,^1,² Sylvain Meylan,^1,² Angela Ciuffi,³ Nadine Zangger,^1,² Giovanna Ambrosini,^1,⁴ Nicolas Dénervaud,¹ Philipp Bucher,^1,⁴ and Didier Trono^1,²^*

¹School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland

²Frontiers-in-Genetics National Center of Competence in Research, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland

³Institute of Microbiology, University Hospital Center and University of Lausanne, Lausanne, Switzerland

⁴Swiss Institute of Bioinformatics, Lausanne, Switzerland

Hiten D. Madhani, Editor

University of California San Francisco, United States of America

* E-mail: didier.trono/at/epfl.ch

Conceived and designed the experiments: ACG SM DT. Performed the experiments: ACG SM. Analyzed the data: ACG SM AC NZ GA ND PB. Contributed reagents/materials/analysis tools: AC NZ GA ND PB. Wrote the paper: ACG DT.

Received October 20, 2009; Accepted February 2, 2010.

This article has been cited by other articles in PMC.

Abstract

Krüppel-associated box domain-zinc finger proteins (KRAB–ZFPs) are tetrapod-specific transcriptional repressors encoded in the hundreds by the human genome. In order to explore their as yet ill-defined impact on gene expression, we developed an ectopic repressor assay, allowing the study of KRAB–mediated transcriptional regulation at hundreds of different transcriptional units. By targeting a drug-controllable KRAB–containing repressor to gene-trapping lentiviral vectors, we demonstrate that KRAB and its corepressor KAP1 can silence promoters located several tens of kilobases (kb) away from their DNA binding sites, with an efficiency which is generally higher for promoters located within 15 kb or less. Silenced promoters exhibit a loss of histone H3-acetylation, an increase in H3 lysine 9 trimethylation (H3K9me3), and a drop in RNA Pol II recruitment, consistent with a block of transcriptional initiation following the establishment of silencing marks. Furthermore, we reveal that KRAB–mediated repression is established by the long-range spreading of H3K9me3 and heterochromatin protein 1 β (HP1β) between the repressor binding site and the promoter. We confirm the biological relevance of this phenomenon by documenting KAP1–dependent transcriptional repression at an endogenous KRAB–ZFP gene cluster, where KAP1 binds to the 3′ end of genes and mediates propagation of H3K9me3 and HP1β towards their 5′ end. Together, our data support a model in which KRAB/KAP1 recruitment induces long-range repression through the spread of heterochromatin. This finding not only suggests auto-regulatory mechanisms in the control of KRAB–ZFP gene clusters, but also provides important cues for interpreting future genome-wide DNA binding data of KRAB–ZFPs and KAP1.

Author Summary

The regulation of gene activity by transcription factors is crucial to the function of all cells. Here, we studied the mechanisms of action of the largest family of gene regulators encoded by the human genome, the so-called KRAB–containing zinc finger proteins (KRAB–ZFPs), which in concert with their universal cofactor KAP1 act as transcriptional repressors. For this, we used two parallel approaches. First, by targeting an ectopic KRAB domain to hundreds of different genes, we found that KRAB/KAP1 can repress promoters located several tens of kilobases from the repressor DNA docking site. We further could show that KRAB induces such long-range effects by mediating the spread of repressive chromatin marks along the body of the gene, resulting in a block of transcriptional initiation at the promoter. In a second set of experiments, we analyzed an endogenous KRAB–ZFP gene cluster, where we could also document KAP1–dependent heterochromatin spreading and transcriptional repression. Together, these results support a model whereby KRAB–ZFPs and KAP1 can mediate long-range transcriptional repression through the spread of silencing chromatin marks. This study thus provides insight into KRAB/KAP1–induced gene regulation at KRAB–ZFP gene clusters, and will further help interpret genome-wide studies of KRAB–ZFPs and KAP1 DNA binding patterns.

Introduction

The proper control of gene expression is paramount to all cellular events, and is orchestrated through a sophisticated balance of activating and repressing influences. Krüppel-associated box domain zinc finger proteins (KRAB-ZFP) constitute the single largest group of transcriptional repressors encoded by the genomes of higher organisms. After appearing in early tetrapods, the KRAB-ZFP family has rapidly expanded and diverged through multiple rounds of gene and segmental duplications, to give rise to more than three hundred and fifty members annotated in both mice and humans [1]–[6]. In spite of their numerical abundance, wide range of tissue-specific expression and dynamic evolutionary history, the physiological functions of KRAB-ZFPs collectively remain ill-defined, and few of their targets have been identified [3],[7]. However, emerging evidence links KRAB/KAP1-mediated regulation to processes as essential and diverse as stem cell pluripotency, early embryonic development and differentiation, genomic imprinting, response to DNA damage and control of behavioral stress [8]–[14]. Furthermore, KAP1 controls endogenous retroviruses in embryonic stem cells, a process crucial for the maintenance of genomic stability [15].

KRAB-ZFPs all harbor a so-called KRAB domain situated upstream of an array of two to forty C2H2 zinc fingers, which provide sequence-specific DNA binding ability [2]. KRAB recruits KAP1 (KRAB-associated protein 1, also known as TRIM28, Tif1β or KRIP-1) [16]–[19], which acts as a scaffold for various heterochromatin-inducing factors, such as heterochromatin protein 1 (HP1), the histone methyltransferase SETDB1, the nucleosome-remodeling and histone deacetylation (NuRD) complex, the nuclear receptor corepressor complex 1 (N-CoR1) and, at least during early embryonic development, de novo DNA methyltransferases [20]–[27].

The phylogenetically conserved family of HP1 proteins is implicated in a variety of nuclear events, such as transcriptional repression and maintenance of chromosome structure [28]. HP1 harbors two major regions: the chromo domain, which binds to repressive di- and trimethylated histone 3 lysine 9 (H3K9me2 and H3K9me3, respectively) residues, and the chromo shadow domain, involved in HP1 homodimerization and recruitment of other partners [29]–[31]. Therefore, HP1 bridges histones with other chromatin-associated proteins, that promote heterochromatin spreading [32]. An advancing front of heterochromatization can thus be propagated by the creation of HP1 binding sites, through histone methyltransferase-mediated H3K9 methylation, followed by the HP1-mediated recruitment of more histone methyltransferase for another round of H3K9 methylation/HP1 binding [33]–[36]. Little is known about the kinetics, efficiency, self-perpetuating ability and action-range of this process. Tethering of a regulated KRAB repressor domain suggested that KRAB/KAP1 induced heterochromatin formation has a rather limited spreading potential in euchromatin, as silencing could be exerted no farther than 2–3 kilobases (kb) away from the repressor binding site [37]. Yet a more recent Chromatin IP (ChIP)-on-chip analysis performed in the human testicular carcinoma cell line Ntera2 revealed close to 7,000 KAP1 binding sites, a number of which were located at the 3′end of KRAB-ZFP genes [38]. While this suggests auto-regulatory negative feedback loops for these genes, such a process would imply that KRAB/KAP1 binding can affect promoters situated at very significant distances. In agreement with such a model, large heterochromatin domains associated with both HP1β and SUV39h1 were found on chromosome 19, where most of the KRAB-ZFP gene clusters reside [39].

To examine the genomic features of KRAB/KAP1-mediated transcriptional regulation, we developed an ectopic repressor assay. In our system, promoterless lentiviral vectors serve as gene traps to drive reporter expression from cellular promoters. Drug-controllable docking of an ectopic KRAB-based repressor then allows an assessment of the effects of KRAB/KAP1 recruitment at these loci. Using this system, we found that KRAB-induced silencing can repress promoters situated several tens of kilobases away from the repressor primary docking site through reduced RNA Pol II binding. Furthermore, we observed that this phenomenon is independent of promoter strength, facilitated by HP1 and associated with spreading of heterochromatin marks between repressor binding site and targeted promoters. Finally, we could document KAP1-mediated transcriptional repression at an endogenous KRAB-ZFP gene cluster by propagation of HP1β and H3K9me3 from the 3′ end of these genes to their transcriptional start site. Our results indicate that KRAB/KAP1 induce long-range repression through the spread of heterochromatin.

Results

KRAB–mediated silencing can act over several tens of kilobases

In order to study the features of KRAB/KAP1-induced silencing, we exposed lentivirally trapped cellular promoters to a drug-regulated KRAB-containing repressor. The tTRKRAB protein contains the KRAB domain of the human KOX1 ZFP fused to the E. coli tetracycline repressor (tTR), and binds to Tet operator sequences (TetO) in a doxycycline (dox)-controllable fashion [40],[41]. We engineered human immunodeficiency virus (HIV)-derived lentiviral (LV) gene trap vectors carrying tandem TetO repeats and a promoter-less puromycin resistance-GFP fusion reporter (puro^R-GFP) downstream of a potent adenoviral splice acceptor. This design predicted that i) reporter expression would occur from the promoters of active genes targeted by the integrants, and ii) dox withdrawal would result in tTRKRAB binding to the TetO sites present in the proviruses, thus exposing the trapped promoters to potential KRAB/KAP1-mediated repression (Figure 1A).

Figure 1

KRAB–mediated silencing can act over several tens of kilobases.

To study the effects of tTRKRAB recruitment to endogenous genes, we engineered HeLa cells stably expressing tTRKRAB and infected those with the TetO-containing LV trapping/silencing (“TrapSil”) vectors at a low multiplicity of infection to limit integrants to one per cell. Clones expressing the trap reporter puro^R-GFP were then isolated and subjected to FACS analysis in the presence (-dox) or absence (+dox) of tTRKRAB binding. Comparing the mean fluorescence intensity values of GFP (MFI) in both conditions revealed heterogeneous effects of tTRKRAB recruitment, with clones exhibiting either high or low levels of reporter silencing (Figure 1B and Table S1). We saw no correlation between the strength of the trapped promoters and their susceptibility to tTRKRAB-mediated repression (Figure 1B). We could also verify proper recruitment of KRAB/KAP1 in different clones by ChIP targeting KAP1 (Figure S1 and data not shown). In addition, we analyzed the effects of tTRKRAB recruitment on populations of LV TrapSil clones and found that the level of GFP was decreased on both the protein and mRNA level upon repressor binding. We also found that the average silencing efficiency of the LV TrapSil population was at around 60%, which is consistent with repressible clones being more prevalent in this population than irrepressible clones (Figure S2).

HIV-derived vectors preferentially integrate in downstream regions of transcribed genes [42],[43]. Therefore, we evaluated the impact of distance between the tTRKRAB docking site and the trapped promoter on KRAB-induced repression in cell populations with distinct LV TrapSil integrants. For this, we used GFP-based cell sorting to purify subpopulations exhibiting trapped promoters, which were silenced by either more than 90% (“repressible”) or less than 10% (“irrepressible”). Proviral integration sites were mapped by ligation-mediated PCR on genomic DNA, followed by large scale pyrosequencing [44],[45]. To ensure an unequivocal analysis, we chose to consider only trapping events that mapped to characterized UCSC known genes with single transcriptional start sites. Using these criteria we were able to obtain 484 repressible integrants in 221 genes and 699 irrepressible integrants in 371 genes (Table S2). When we compared the distance between the proviral integration sites and their trapped promoters in both populations, we found that irrepressible promoters were on average located further away from the repressor-docking site than their repressible counterparts (Figure 1C). We observed a similar pattern when analyzing trapped genes harboring multiple integration events from both the repressible and the irrepressible groups. By analyzing six of the most targeted genes, we found that the provirus integrated closest to the promoter systematically induced a repressible phenotype, while integrants located further downstream often fell into the irrepressible group (Figure S3). The comparison between the effect of tTRKRAB-mediated silencing and the range of its action also revealed that repression could be exerted at several tens of kilobases away from the tTRKRAB binding site, albeit with a sharp decrease at around 15 kb (Figure 1C). Finally, the overlap between the two histograms indicates that distance between repressor-binding site and the targeted promoter was not the sole determinant of susceptibility to KRAB-mediated repression, but that sequence- or locus-dependent factors must also play a role.

Silenced promoters exhibit increased H3K9me3 and decreased H3-acetylation and RNA Pol II recruitment

To investigate the molecular bases of repressible and irrepressible phenotypes we characterized four clones (I, IX, XI and XVI) in more detail using 5′RACE to map the trapped fusion transcripts (Table S1) and by performing ChIP studies. All clones had the integrant within 10–20 kb from the trapped promoter, yet three (clones I, IX and XI) were repressible and one (clone XVI) was irrepressible (Figure 1B). The level of histone modifications known to be affected by KRAB/KAP1-mediated repression, namely histone H3-acetylation (H3Ac) and H3K9me3 was evaluated at each of the four trapped promoters. In all three repressible clones, H3Ac levels dropped significantly upon tTRKRAB-binding, consistent with a close to complete deacetylation of the trapped allele (Figure 2A). In contrast, the irrepressible promoter of the casein kinase 1d (CSNK1D) gene, trapped in clone XVI, exhibited constant levels of this chromatin mark (Figure 2A). An inverse image was obtained for H3K9me3, with increases of different magnitudes at the three silenced promoters upon dox removal. No such change in H3K9me3 levels was observed at the irrepressible promoter (Figure 2B). To exclude the possibility that increased H3K9me3 was due to chromatin condensation, we also monitored levels of total histone H3 at these promoters. Comparison of the dox+ and dox- conditions indicated that the observed H3K9me3 increases could not be solely attributed to changes in total H3 (Figure 2C). In order to probe the impact of KAP1-induced repression on transcriptional initiation, we measured the levels of hypophosphorylated RNA Pol II present at trapped promoters (Figure 2D). We found that tTRKRAB binding resulted in a significant decrease of the recruitment of this enzyme in repressible but not irrepressible clones. Noteworthy, the remaining levels of Pol II binding detected in the absence of dox most likely originated from the sister allele of the trapped promoter.

Figure 2

Silenced promoters exhibit increased H3K9me3 and decreased H3-acetylation and RNA Pol II recruitment.

HP1β and H3K9me3 spread from the KRAB–binding site to the promoter

To investigate further the mechanisms of this long-range repression, we measured levels of both H3K9me3 and HP1β enrichment between the promoters and the TrapSil integrants of the loci under study. It was previously demonstrated that all HP1 isoforms can interact with KAP1 [20],[22],[29],[46],[47]. However, recent FRET results revealed distinct physical interactions between KAP1 and HP1β or HP1γ during differentiation of embryonal carcinoma F9 cells to primitive endoderm-like cells [48]. Interestingly, only the KAP1-HP1β interaction induced a relocalization of genes from euchromatin to heterochromatin upon differentiation, suggesting a specific role for this isoform during KAP1-mediated gene silencing. In addition, both HP1β and KAP1 co-localize at the 3′end of KRAB-ZFP genes [38],[39], making HP1β a likely candidate to be involved in propagating KRAB/KAP1-mediated long-range repression.

Indeed, global levels of H3K9me3 and HP1β increased along the body of trapped genes upon tTRKRAB recruitment, and the magnitude and efficiency of propagation of this phenomenon correlated with its phenotypic outcome (Figure 3, Figure S4, Figure S5). Clone I, where the prostaglandin E synthase 3 (PTGES3) gene was trapped, exhibited high reporter GFP expression, which was completely silenced upon tTRKRAB binding (Figure 1B). Correspondingly, KRAB repressor binding led to a massive increase of both HP1β and H3K9me3 along this gene (Figure 3B). In comparison, the irrepressible CSNK1D gene body, showed smaller changes in HP1β and H3K9me3 loading, which gradually decreased towards the promoter (Figure 3A). The ChIP efficiencies between the +dox and the –dox samples were comparable for both factors, as shown by analysis of several control loci (Figure S6). We additionally monitored the levels of H3K9me3 and HP1β at the trapped and repressible KRAB-ZFP gene ZNF77. More specifically, tTRKRAB binding to the 3′end of ZNF77-1 led to an increase of both factors at the promoters and gene bodies of ZNF77-1, as well as the neighboring ZNF57 gene, which is located approximately 40 kb away (Figure S5). Therefore, we conclude that KRAB/KAP1 recruitment triggers the long-range spreading of silent chromatin marks. Efficient transcriptional silencing, however, is only established when this process extends over the promoter itself, interfering with proper RNA Pol II recruitment and transcriptional initiation (Figure 2D).

Figure 3

HP1β and H3K9me3 spread from the KRAB–binding site to the promoter.

HP1 facilitates KRAB/KAP1–mediated long-range repression

We further explored the role of HP1 recruitment to KAP1 in KRAB/KAP1-mediated long-range repression by performing complementation studies in a KAP1 knockout mouse embryonic fibroblast (MEF) cell line stably expressing tTRKRAB. For this, we isolated two LV TrapSil clones, and used 5′RACE to map the fusion transcripts resulting from promoter trapping. Since both proviruses were intronic, we assigned integration to the middle of the targeted intron. This revealed that clone 1 had integrated in the gene Sestd1 (SEC14 and spectrin domains 1) at a distance of 15 kb from the promoter and that clone 2 had trapped the gene Prcp (prolylcarboxypeptidase) and was situated 25 kb away from its TSS (Table S1). Both clones were therefore suitable to serve in our study of KRAB/KAP1-mediated long-range repression. We found that GFP expression was independent of the presence of dox (Figure 4A), as predicted from the knockout of KAP1 (Figure 4B). However, silencing could be significantly restored in both clones by lentiviral vector-mediated complementation with wild type KAP1 (Figure 4A and 4B). By comparison, the transduction with a vector expressing a KAP1 mutant defective for HP1 binding (KAP1^{R487E, V488E}) [22] only rescued tTRKRAB-mediated silencing to levels that were about six times lower than observed with wild type KAP1 in both clones (Figure 4A and 4B). This strongly suggests that HP1 recruitment by KAP1 facilitates the long-range impact of the KRAB/KAP1 silencing complex.

Figure 4

HP1 facilitates KRAB/KAP1–mediated long-range repression.

KAP1 represses KRAB-ZFP genes via 3′ to 5′ long-range heterochromatin spreading

Recent genome-wide studies have identified KAP1 and HP1β binding sites at the 3′end of KRAB-ZFP genes, suggesting long-range auto-regulatory feedback loops [38],[39]. To investigate this possibility directly, we generated a KAP1 knockdown HeLa cell line by transduction with an shRNA-expressing lentiviral vector. KAP1 mRNA and protein levels were about ten-fold lower in these cells than in their wild type control (Figure S7). We then analyzed a group of five KRAB-ZFPs located on human chromosome 19, where most members of the family reside in closely packed clusters (Figure 5, Figure 6, and Figure 7). In KAP1 knockdown cells, four out of five analyzed cluster-specific transcripts, namely ZNF554, ZNF555, ZNF556 and ZNF57, showed slight upregulation, which reached significance in the latter three genes, when compared to a cell line expressing an shRNA directed against GFP (Figure 5A). This is consistent with a de-repression of these transcripts following KAP1 depletion. The ZNF77 gene is predicted to produce two partly overlapping transcripts, ZNF77-1 and ZNF77-2. However, only ZNF77-1 was reliably detected in our cells and its levels were identical in KAP1 knockdown and control cells (Figure 5A and data not shown). We further verified that KAP1 depletion did not lead to a global decrease in heterochromatin protein expression, since four out of five analyzed genes displayed similar RNA levels in cells expressing either shKAP1 or shGFP (Figure 5B).

Figure 5

KAP1 mediates transcriptional repression of KRAB–ZFPs.

Figure 6

KAP1 is enriched at the 3′end of KRAB–ZFPs.

Figure 7

KAP1 mediates the spread of HP1β and H3K9me3 at KRAB–ZFP clusters.

KAP1-specific ChIP identified KAP1 binding sites at the 3′end of ZNF554, ZNF555 and ZNF556, which are located more than 10 kb from the corresponding promoters, as well as within the region of ZNF77-1/2 overlap (Figure 6B–6E). We then examined the presence of HP1β and H3K9me3 throughout the five genes present in the cluster. In ZNF554, ZNF555, ZNF556 and ZNF77, we found peaks of HP1β and H3K9me3 enrichment at the KAP1 binding site, the intensity of which progressively decreased towards the promoter (Figure 7A–7D). We also observed elevated levels of the two silencing marks in the middle of ZNF57, for which no KAP1 association had been detected, possibly because it was not covered by our qPCR primer sets (Figure 7D). When KAP1 was targeted by shRNA, levels of HP1β and H3K9me3 decreased, consistent with KAP1 mediating the spread of silent chromatin at the affected gene bodies (Figure 7A–7D). Therefore, we conclude that KAP1 represses KRAB-ZFP gene clusters by mediating 3′ to 5′ spreading of HP1β and H3K9me3 from endogenous binding sites located at the distal end of these genes. This closely resembles results obtained with the ectopic TrapSil vectors, demonstrating the biological significance of KRAB/KAP1 in long-range repression through the propagation of heterochromatin-associated modifications.

Discussion

In the present work, we demonstrate that KRAB-ZFPs and their co-repressor KAP1 can induce long-range transcriptional repression through heterochromatin spreading. For this, we studied lentivirally-trapped promoters silenced by an ectopic KRAB-containing repressor, as well as the KAP1-mediated modulation of the chromatin landscape at an endogenous KRAB-ZFP gene cluster. We could document the propagation of HP1β and the silencing mark H3K9me3 from repressor binding sites to targeted promoters, and showed that gene silencing results from a block in transcriptional initiation. Furthermore, we could demonstrate that this silencing process is KAP1-dependent and facilitated by HP1.

Our results demonstrate that KRAB/KAP1 have the ability to repress promoters at far greater distances than suspected based on earlier publications. Indeed, it was previously found that recruitment of a regulatable KRAB-based repressor to an array of integrated reporter genes led to transcriptional repression of the promoter adjacent to the repressor binding site, but did not affect a transcriptional unit located 2–3 kb away [37]. Here, we frequently observed repression of promoters located several tens of kilobases away from the tTRKRAB binding site, although this process was most efficient for promoters located within 15 kb or less. These results are consistent with a recent KAP1 ChIP-chip analysis, which documented the association of this co-repressor with promoters and the 3′end of KRAB-ZFP genes, and corroborate data indicating that the KRAB-ZFP ZNF263 is often found within the body of genes [38],[49]. The novel approach used here, which allows targeting of chromatin-modifying proteins to hundreds of active genes at once, provides a functional correlation to these mapping studies.

Heterochromatin domains at KRAB–ZFP genes

The large family of KRAB-ZFP genes is the result of a recent evolutionary phenomenon, with iterative gene duplication events that have led to large sequence homologies between its different members. Therefore, it was proposed that the large heterochromatin domains containing HP1β and SUV39h1 at KRAB-ZFP genes act as a protection from deleterious recombination events, that such a homology would enhance [39]. We show here that KAP1, very likely in cooperation with specific KRAB-ZFPs, binds to the 3′end of KRAB-ZFP genes and mediates the spread of heterochromatin through their gene bodies, thus providing a molecular mechanism for the formation of these heterochromatic domains at KRAB-ZFP gene clusters. Furthermore, by examining a specific KRAB-ZFP gene cluster, we documented the KAP1-mediated control of its transcriptional units, by demonstrating increased expression of KRAB-ZFP genes upon KAP1 depletion. Importantly, the absence of silent chromatin marks at the promoters of these KRAB-ZFP genes, as well as the range of the transcriptional modulation observed, may be consistent with auto-regulatory feedback loops between different KRAB-ZFPs. We further propose that heterochromatic modifications concentrated in gene bodies of KRAB-ZFPs rather than in their promoter, ensure protection from recombination events, while preserving the transcriptional potential of KRAB-ZFP genes. Further analysis of specific KRAB-ZFPs and their binding profiles will provide more functional insight into this type of regulation.

Mechanism of KRAB/KAP1–induced heterochromatin spreading

Spreading is considered to be an intrinsic property of heterochromatin [50]. For mammals it has been proposed that initial heterochromatin assembly at the centromere is followed by a spreading process involving HP1 and H3K9me3 in a self-perpetuating loop. It remains unclear what brings about the initial H3K9me3 that nucleates heterochromatin formation. In addition, in mammalian cells spreading of heterochromatin relies on HP1 in cooperation with additional factors such as SUV39H, since neither H3K9me2 nor H3K9me3 alone is sufficient to recruit HP1 to chromatin [51]. By analogy, KRAB/KAP1-mediated silencing may depend on similar interactions between KAP1, HP1 and H3K9me3, since long-range repression is abrogated when KAP1-HP1 binding is lost, and effective silencing is only seen if levels of H3K9me3 increase at the affected promoters themselves. In agreement with this model, recent studies on the MEST promoter, a primary KAP1 target, showed that repression is rapidly lost upon disruption of the interaction between KAP1 and HP1, and is accompanied by a decrease in H3K9me3 [52].

Remarkably, we documented a decreasing gradient of HP1β and H3K9me3 from the KAP1 docking site to the promoter in both ectopically KRAB-repressed genes and at an endogenous KRAB-ZFP gene cluster. Although it is difficult to compare the binding curves of these different factors due to the non-linear binding kinetics of the antibodies used for their capture, this may suggest that KAP1 nucleates the spread of heterochromatin. Notably, KAP1 is enriched along the bodies of all KRAB-ZFP genes (Figure 6B–6E), consistent with a model of low levels of KAP1 spreading away from KAP1 primary binding sites. It has been proposed that the propagation of heterochromatin could occur by diffusion of an enzymatically-competent macromolecular complex from a source site to nearby lower affinity sites [53]. Our data suggest that KAP1, HP1 and specific histone methyltransferases may form such a complex to mediate the spread of silent chromatin.

Consequences of KRAB/KAP1–induced heterochromatin spreading

Our analyses of TrapSil-targeted genes and of an endogenous KRAB-ZFP gene cluster indicate not only that in both settings silent chromatin marks spread from sites of KAP1 recruitment, but also reveal that the transcriptional consequences of this event vary. Our finding that KAP1 depletion upregulated several transcripts within a KRAB-ZFP cluster, even in the absence of promoter silencing marks suggests that expression was modulated primarily by chromatin changes over the transcribed region of these genes. Therefore, our results imply broad functions for KRAB/KAP1-mediated long-range repression, encompassing not only transcriptional promoter silencing but also transcriptional dampening, as has previously been described for HP1 binding within the coding regions of genes [54],[55].

Our analysis of trapped cellular genes exposed to the tTRKRAB ectopic repressor revealed outcomes largely dictated by whether or not silencing marks reached their promoters. It is likely that the cumulative effect of increased H3K9me3 and decreased H3Ac interferes with binding of transcriptional activator proteins and the ultimate recruitment of RNA Pol II. This hypothesis is not only consistent with our finding of less hypophosphorylated RNA Pol II at trapped and silenced promoters, but was also recently put forward to explain how histone modifications regulate the expression of the L-type pyruvate kinase gene [56].

Analysis of trapped populations with different silencing phenotypes showed that the distance between sites of KRAB/KAP1 recruitment and trapped promoters had a strong influence on KRAB-mediated silencing, but was not the sole factor at play. Indeed, some promoters located at great distances from the tTRKRAB-binding site were efficiently shut off, while others situated close-by escaped this process. Several possibilities come to mind to explain this variation. First, distinct subgroups of promoters may exhibit intrinsically different susceptibility to KRAB/KAP1-mediated silencing. Second, counteracting influences may be exerted at some promoters, for instance downregulation of H3K9me3 levels by histone demethylases, as exemplified by the proposed action of JMJD1A and JMJD2C on pluripotency-associated genes in mouse ES cells [57],[58]. Third, an active nucleosome remodeling mechanism may protect a promoter from heterochromatization. This is consistent with our observation that the irrepressible trapped promoter of the casein kinase 1d locus exhibited lower levels of histone 3 upon tTRKRAB recruitment (Figure 2C). By analogy, the induction of the mammalian interferon-β promoter through external cues is also accompanied by nucleosome remodeling [59].

More generally, there may be genomic loci that exhibit chromatin environments permissive or restrictive to KRAB/KAP1-mediated long-range repression. Examples of both situations in this study include the long-range spreading of heterochromatin detected over a 50 kb long region harboring several KRAB-ZFP genes and the impaired spreading observed in the vicinity of the HNRNPA2B1-CBX3 ubiquitously-acting chromatin opening element (A2UCOE) [60]. In the first of these two cases the site of tTRKRAB recruitment at the 3′end of the ZNF77-1 gene coincided with an endogenous KAP1 binding site, so that the impact of tTRKRAB and KAP1 recruitment to this locus could be directly compared (Figure 6E and Figure S5). Their effect on promoters seemed to differ, as low basal levels of H3K9me3 at the ZNF77-1 and ZNF57 promoters dramatically increased upon tTRKRAB recruitment, which correlated with transcriptional repression. This result is consistent with our model that elevated H3K9me3 at active promoters leads to efficient transcriptional repression by KRAB/KAP1. If levels of this histone modification, however, are already low, then depletion of KAP1 will not affect the promoter as strongly. In contrast, we found that increases of both H3K9me3 and HP1β along the gene body were dependent on both KAP1 and tTRKRAB, consistent with the idea that both factors mediate the spread of heterochromatin from an initial binding site. We further speculate that specific regions, such as KRAB-ZFP gene clusters, are primed for KRAB/KAP1-mediated long-range repression, where this mechanism contributes to heterochromatin structure and long-range transcriptional repression.

In a second example we made use of trapped genes containing several TrapSil integrations with distinct silencing phenotypes (Figure S3). One of these genes is HNRNPA2B1, the promoter of which was shown to be part of a UCOE exhibiting a methylation-free CpG island (Figure S3E). This island encompasses the two divergently transcribing promoters of the housekeeping genes HNRNPA2B1 and CBX3 and allows for stable transgene expression in centromeric heterochromatin [60]. Although in our study the integrant closest to the promoter exhibited a repressible phenotype, the first irrepressible integrant was located at a distance of 3.5 kb, which is considerably lower than the observed average of 15 kb. Moreover, this gene contained a very high ratio of irrepressible/repressible integrants. It was recently shown that the DNA CpG-methylation free region spans 3 kb into the HNRNPA2B1 gene body, probably allowing for very efficient transcriptional initiation [61]. Therefore, it is conceivable that the high frequency of irrepressible clones found here is linked to the specific architecture of this promoter. We are currently trying to establish more detailed correlations between the effect of tTRKRAB-mediated silencing on genes and their individual chromatin environment.

In some cases, long-range acting elements such as enhancers might be the primary targets of KRAB/KAP1-induced chromatin modifications, with a secondary impact on promoters. More generally, specific cis-acting sequences or chromatin features located between KAP1 recruitment sites and transcriptional start sites may alter the spread of heterochromatin, either positively or negatively. Conversely, while the hereby observed long-range spreading of silencing heterochromatin suggests that some clusters of genes might be co-regulated through common KRAB-ZFP binding sites, it also implies the existence of limiting mechanisms to avoid that it systematically be the case.

Perspectives

In conclusion, our results reveal that KRAB-ZFPs, which constitute the largest family of transcription factors encoded by the human genome, mediate long-range repression through the spread of heterochromatin in a process influenced by the genomic context. Therefore, it is tempting to speculate that the long-range epigenetic silencing (LRES) of large chromosomal regions observed in cancer cells [62] may be due to aberrant spreading of heterochromatin stemming from the dysregulation of KRAB/KAP1-mediated epigenetic regulation. Indeed, the expression of KAP1 has been found to be altered in several types of malignancies, for instance lung cancer [63],[64]. Furthermore, since tTRKRAB recruitment results in long-range repression at permissive endogenous loci, its targeting to specific chromosomal regions, through the site-specific introduction of TetO sequences, could be a useful tool to elucidate functions of interesting regulatory landscapes. Finally, the controlled recruitment of KRAB-containing repressors should facilitate studies on the kinetics of heterochromatin spreading in different chromosomal environments, at distinct developmental stages, and in various physiopathogical conditions, thereby shedding light on the intricacies of transcriptional repression.

Materials and Methods

Cell lines and culture

HeLa, 293T and HCT116 cells (from ATCC, Manassas, VA) were cultured using standard methods. Mouse embryonic fibroblasts from homozygous KAP1 floxed mice were derived as described [27]. Doxycycline (Sigma Aldrich) was used at a concentration of 1 µg/mL. GFP expression was determined by flow cytometry using the Beckton Dickinson FACScan.

Vectors

LV- based gene trap vectors [65] were modified by introducing 7 TetO repeats. pLV-tTR-KRAB-Red was described [66]. LV-based plasmids expressing human wild type or mutant (R487E, V488E) KAP1 were a gift from Frank Rauscher (Wistar Institute). LV-based plasmids expressing shRNAs targeting KAP1 and GFP were purchased from Sigma and Addgene, respectively (for sequences see Table S3). LV- and MLV-based particles were produced as described [67].

Western blot

Cells were lysed using standard procedures and subjected to SDS-PAGE. Antibodies used were mouse monoclonal anti-PCNA (Oncogene Research Products, Boston, MA) and mouse monoclonal anti-KAP1 (MAB 3662, Chemicon), followed by secondary goat anti-mouse IRDye 800CW antibody (Li-cor). Imaging and quantification of the Western blot were performed on ODYSSEY infrared imager using the machine-specific software (Li-cor).

RNA procedures

Total RNA was isolated using an RNeasy plus kit (Qiagen).

For Northern blot 15 µg of RNA were run on a denaturing agarose gel and transferred onto a nylon membrane, which was pre-hybridized in ULTRAhyb buffer (Ambion) and after addition of a GFP-specific DIG-labelled probe (Roche, PCR DIG probe synthesis) hybridized overnight at 42°C. The membrane was then subjected to washing and detection procedures as recommended by the manufacturer (Roche, DIG Northern kit). For 5′RACE a previously described protocol was used [68]. Briefly, RNA was converted into cDNA with Super Script II Reverse Transcriptase (Invitrogen) and the trap-specific reporter primer RSP-6 or RSP-8 (for sequences see Table S3). Primary PCR reactions were performed using 2 µl of dA-tailed cDNA products and the reporter specific primers IRES(−) or RSP-6 and anchoring primers QA and QT. Nested PCR reactions were done using the reporter specific primers IRES1b or IRES(−) and anchoring primer QB. PCR products were gel purified and cloned into pCR4 by using a TOPO TA cloning kit (Invitrogen) and subsequently sequenced. Sequences were extracted from chromatograms using phred [69]. pCR4 and virus specific sequences were identified with BLAT [70] and removed prior to annotation. The location of each amplified cDNA fragment was determined with a blastn search for exact matches on the human genome (ensembl version, August 2007). Since most integrations were intronic and thus not present in the amplified cDNA, we arbitrarily assigned the proviral integration to the middle of the intron.

For quantitative Real-Time PCR (qPCR), RNA was converted into cDNA as in 5′RACE, however using random hexamer primers. For each sample a control reaction lacking RT was included. Diluted cDNA was used to assay the expression of each gene with 1x Power Sybr (Applied Biosystems), and a 400 nM concentration of each gene-specific primer. To verify specificity, each PCR was followed by a melting curve analysis. The increase in fluorescence was analysed with the SDS software, version 2.2.2 (Applied Biosystems). For all amplification plots the baseline data were set with the automatic cycle threshold function. A mean quantity was calculated from duplicate or triplicate PCRs for each sample, and a normalization factor was obtained by calculating the geometric mean of the values of the selected housekeeping genes, which was subsequently used to normalize the relative amounts of RNAs of interest [71].

Ligation-mediated PCR (LM–PCR)

Integration sites were mapped by LM-PCR using a previously described protocol [44],[45]. Integration site determination was processed by selecting integrations site amplicons satisfying all quality controls (4/5/4 primer sequence and co-amplified viral genome fragment) and using FetchGWI [72] for alignment on the March 2006 Human Genome assembly. The gene list used for integrant mapping was downloaded from the UCSC table browser by selecting “UCSC known genes”, further including SpID and genesymbol fields. TxStart and TxEnd were downloaded to determine the transcriptional start site and 3′ end for corresponding genes.

Chromatin immunoprecipitation (ChIP)

ChIP was performed according to the Upstate protocol with minor modifications (http://www.millipore.com/techpublications/tech1/mcproto407). About 10⁷ cells were crosslinked with 1% formaldehyde for 15 min at RT, quenched by adding glycine and rinsed with PBS. Cells were then resuspended in lysis buffer and sheared by sonication with a Branson digital sonicator (model 250) on ice four times for 20 s. 100 µl of sonicated chromatin was directly de-crosslinked and used as the total input (TI) reference in quantitative PCR analysis at a dilution of 1 [ratio]

350 or 1

700. 100 µl of sonicated chromatin was used for each ChIP reaction and was diluted in 900 µl dilution buffer and precleared with 80 µl salmon-sperm DNA- proteinA agarose beads (Upstate). Precipitating antibodies H3Ac (Upstate, 06-599), H3K9me3 (Abcam ab 8898), total H3 (Abcam ab 1791), HP1β (Abcam, ab 49938), RNA Pol II 8WG16 (Abcam, ab 817) and rabbit polyclonal KAP1 (gift from D. Schultz) were added overnight (IP). Chromatin-antibody complexes were then captured, washed and eluted with 100 mM NaHCO₃, 1% SDS. Cross-links between DNA and proteins were reversed by addition of NaCl (for H3Ac, H3K9me3, total H3, HP1β) and incubation at 65°C. DNA was precipitated after incubation with RNase A and Proteinase K and resuspended in 50 µl H₂O and subjected to quantitative PCR analysis at a dilution of either 1 [ratio]

4 or 1

8. Quantitative PCR was done as stated in the RT-PCR section and primer sequences are available in Table S3. Negative control reactions without antibody were run for each sample and in all cases gave negligible results. To quantify the relative enrichment of proteins or specific histone modifications at a given sequence a ratio between the relative quantities of IP and TI was calculated with the help of a standard curve and expressed as the % of input value. In the case of H3Ac, total H3 and RNA Pol II these values were further normalized to the relative enrichment value of GAPDH. In all the other cases ChIP efficiency was comparable between samples as measured at specific control loci, including p53BP2 for H3K9me3 [73] and ZNF556 for HP1β (Figure S6 and Figure S8). Fold changes were calculated as ratios of -dox/+dox enrichments, with the ratios of the respective positive controls set as 1, the controls consisted of p53BP2 for H3K9me3 and of ZNF556 for HP1β (Table S4).

Supporting Information

Figure S1

KAP1 is recruited to the provirus of repressible and irrepressible clones. The relative enrichment (% of input) of KAP1 was quantified by ChIP analyses on the provirus of the irrepressible casein kinase 1d (CSNK1D) gene, the repressible prostaglandin E synthase 3 (PTGES3) gene and at the repressible ephrin receptor B4 (EPHB4) gene. All values are expressed as means +SEM of duplicate experiments.

(0.06 MB TIF)

Click here for additional data file.^{(55K, tif)}

Figure S2

tTRKRAB mediates silencing of retrovirally trapped cellular promoters in a selected HeLa cell population. (A) Mechanism of promoter trapping and dox-controllable silencing following TrapSil-mediated transduction of tTRKRAB-expressing cells. (B,C) tTRKRAB-expressing HeLa cells infected with LV- or MLV-derived TrapSil vectors were selected in puromycin before dox withdrawal and analysis by (B) FACS (GFP positive cells are in upper left quadrant) and (C) Northern blot (using a GFP-specific probe). The lower panels show a longer exposure of the blot and the 28S rRNA loading control.

(2.17 MB TIF)

Click here for additional data file.^{(2.0M, tif)}

Figure S3

Trapped genes carrying multiple TrapSil integrants of different silencing phenotypes. (A–F) Graphic depiction of six trapped genes carrying multiple LV-based TrapSil integrants with different silencing phenotypes: (A) calnexin precursor (CANX), (B) nucleophosmin 1 (NPM1), (C) sodium-coupled neutral amino acid transporter A2 (SLC38A2), (D) Y-box binding protein 1 (YBX1), (E) heterogeneous nuclear ribonucleoproteins A2/B1 (HNRNPA2B1), and (F) HNRNP C1/C2 (HNRNPC). The proviruses exhibiting a repressible phenotype (>90% silencing) are depicted as black triangles above the baseline, whereas the irrepressible counterparts (<10% silencing) are depicted as light grey triangles below the baseline.

(0.27 MB TIF)

Click here for additional data file.^{(268K, tif)}

Figure S4

HP1β and H3K9me3 spread from the KRAB-binding site to the promoter. ChIP analyses quantifying the relative enrichment (% of input) of both HP1β and H3K9me3 were performed for the repressible clone IX in the presence and absence of tTRKRAB binding. The interrogated sequence at the ephrin receptor B4 (EPHB4) locus spanned from the proviral tTRKRAB binding sites (light grey circles) to the trapped promoter. qPCR amplicons are depicted as letters and are not drawn to scale. All values are expressed as means +SEM of triplicate experiments. Fold changes were calculated as ratios of -dox/+dox enrichments, with the ratios of the respective positive controls set as 1. The controls consisted of p53BP2 for H3K9me3 and of ZNF556 for HP1β (Table S4).

(0.12 MB TIF)

Click here for additional data file.^{(118K, tif)}

Figure S5

HP1β and H3K9me3 spread along the 50 kb-spanning ZNF77/57 locus upon tTRKRAB binding. ChIP analyses quantifying the relative enrichment (% of input) of both HP1β and H3K9me3 were performed for the repressible clone XI in the presence and absence of dox. The ZNF77-1 promoter drives the expression of the integrated TrapSil provirus (depicted as light grey circles). The interrogated sequence spanned the whole ZNF77/57 locus and the respective qPCR amplicons are depicted as letters and are not drawn to scale. All values are expressed as means +SEM of triplicate experiments. Fold changes were calculated as ratios of -dox/+dox enrichments, with the ratios of the respective positive controls set as 1. The controls consisted of p53BP2 for H3K9me3 and of ZNF556 for HP1β (Table S4).

(0.30 MB TIF)

Click here for additional data file.^{(296K, tif)}

Figure S6

ChIP analyses of different TrapSil clones at control loci. (A–D) We ensured that there were similar amounts of ChIP material in the +dox compared to the -dox samples for each TrapSil clone by analyzing HP1β or H3K9me3 relative enrichment levels at control loci. These control loci are: p53BP2, ZNF554, ZNF555, ZNF556, ZNF77-1, and ZNF77-2 and were probed in the TrapSil clones (A) XVI, (B) I, (C) XI, and (D) IX.

(0.37 MB TIF)

Click here for additional data file.^{(366K, tif)}

Figure S7

Quantification of the KAP1 knockdown efficiency in a stable HeLa cell line. HeLa cells were stably transduced with lentiviruses expressing shRNA targeting either KAP1 or GFP. The levels of knockdown were quantified by using (A) qPCR measurements normalized to EEF1α and (B) western blot analyses for KAP1 levels with PCNA as a loading control. The qPCR values are expressed as means +SEM of triplicate experiments.

(0.12 MB TIF)

Click here for additional data file.^{(120K, tif)}

Figure S8

ChIP analyses of knockdown HeLa cell lines at control loci. We ensured that there were similar amounts of ChIP material in the shKAP1, compared to the shGFP cell lines, by analyzing HP1β or H3K9me3 relative enrichment levels at the control gene p53BP2 and at the human satellite 2 repeats (Sat2).

(0.08 MB TIF)

Click here for additional data file.^{(79K, tif)}

Table S1

Characterization and integration site mapping of LV-based TrapSil HeLa and MEF clones.

(0.08 MB DOC)

Click here for additional data file.^{(77K, doc)}

Table S2

Integration site mapping of the repressible and the irrepressible LV TrapSil populations. After proviral integration site mapping, the position of TrapSil integrants was determined relative to the trapped promoter for both the repressible and the irrepressible populations. More specifically, we computed two gene lists (LV repressible and LV irrepressible), which included all the trapped promoters with a single annotated transcriptional start site (TSS), which were further used to build Figure 1C. Each of these gene lists are represented by a table containing gene name (column 1), gene length (column 2), and integrant coordinate relative to the TSS of the gene (column 3). The gene name column can contain more than one identifier due to UCSC known gene annotation, which merges different gene annotations.

(0.06 MB XLS)

Click here for additional data file.^{(59K, xls)}

Table S3

List of primers used in this study.

(0.05 MB XLS)

Click here for additional data file.^{(46K, xls)}

Table S4

Fold change values calculated for different TrapSil clones by comparing ChIP enrichments in the presence or absence of tTRKRAB binding.

(0.04 MB XLS)

Click here for additional data file.^{(37K, xls)}

Acknowledgments

We thank Alexey Ivanov, Frank Rauscher, David Schultz, Francesca Santoni de Sio, and Luigi Naldini for the gift of reagents; Andrea Corsinotti for the KAP1 knockout MEF cell line; and Audrey Letourneau, Sujana Nylakonda, Sandra Offner, and Sonia Verp for technical assistance. In addition, we are grateful to Rick Bushman and Nirav Malani for help with analyzing TrapSil integrant sequences, Natalia Soshnikova for advice on chromatin IP experiments, Patrick Tschopp for help with the graphical representations, Cristina Marquez for support in the statistical analyses, and Miguel Garcia for cell sorting. We acknowledge Dirk Schübeler, Gwen MacDonald, Nahid Iglesias, and members of the Trono lab, especially Helen Rowe and Johan Jakobsson, for helpful comments on the manuscript.

Footnotes

The authors have declared that no competing interests exist.

This work was supported by grants from the Swiss National Science Foundation, the Infectigen Association, and the European Union. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

References

1. Ding G, Lorenz P, Kreutzer M, Li Y, Thiesen HJ. SysZNF: the C2H2 zinc finger gene database. Nucleic Acids Res. 2009;37:D267–273. [PMC free article] [PubMed]

2. Emerson RO, Thomas JH. Adaptive evolution in zinc finger transcription factors. PLoS Genet. 2009;5:e1000325. doi: 10.1371/journal.pgen.1000325. [PMC free article] [PubMed]

3. Huntley S, Baggott DM, Hamilton AT, Tran-Gyamfi M, Yang S, et al. A comprehensive catalog of human KRAB-associated zinc finger genes: insights into the evolutionary history of a large family of transcriptional repressors. Genome Res. 2006;16:669–677. [PubMed]

4. Tadepally HD, Burger G, Aubry M. Evolution of C2H2-zinc finger genes and subfamilies in mammals: species-specific duplication and loss of clusters, genes and effector domains. BMC Evol Biol. 2008;8:176. [PMC free article] [PubMed]

5. Urrutia R. KRAB-containing zinc-finger repressor proteins. Genome Biol. 2003;4:231. [PMC free article] [PubMed]

6. Vaquerizas JM, Kummerfeld SK, Teichmann SA, Luscombe NM. A census of human transcription factors: function, expression and evolution. Nat Rev Genet. 2009;10:252–263. [PubMed]

7. Hamilton AT, Huntley S, Tran-Gyamfi M, Baggott DM, Gordon L, et al. Evolutionary expansion and divergence in the ZNF91 subfamily of primate-specific zinc finger genes. Genome Res. 2006;16:584–594. [PubMed]

8. Cammas F, Herzog M, Lerouge T, Chambon P, Losson R. Association of the transcriptional corepressor TIF1beta with heterochromatin protein 1 (HP1): an essential role for progression through differentiation. Genes Dev. 2004;18:2147–2160. [PubMed]

9. Cammas F, Mark M, Dolle P, Dierich A, Chambon P, et al. Mice lacking the transcriptional corepressor TIF1beta are defective in early postimplantation development. Development. 2000;127:2955–2963. [PubMed]

10. Hu G, Kim J, Xu Q, Leng Y, Orkin SH, et al. A genome-wide RNAi screen identifies a new transcriptional module required for self-renewal. Genes Dev. 2009;23:837–848. [PubMed]

11. Jakobsson J, Cordero MI, Bisaz R, Groner AC, Busskamp V, et al. KAP1-mediated epigenetic repression in the forebrain modulates behavioral vulnerability to stress. Neuron. 2008;60:818–831. [PubMed]

12. Li X, Ito M, Zhou F, Youngson N, Zuo X, et al. A maternal-zygotic effect gene, Zfp57, maintains both maternal and paternal imprints. Dev Cell. 2008;15:547–557. [PMC free article] [PubMed]

13. Tian C, Xing G, Xie P, Lu K, Nie J, et al. KRAB-type zinc-finger protein Apak specifically regulates p53-dependent apoptosis. Nat Cell Biol. 2009;11:580–591. [PubMed]

14. Zheng L, Pan H, Li S, Flesken-Nikitin A, Chen PL, et al. Sequence-specific transcriptional corepressor function for BRCA1 through a novel zinc finger protein, ZBRK1. Mol Cell. 2000;6:757–768. [PubMed]

15. Rowe HM, Jakobsson J, Mesnard D, Rougemont J, Reynard S, et al. KAP1 controls endogenous retroviruses in embryonic stem cells. Nature. 463:237–240. [PubMed]

16. Abrink M, Ortiz JA, Mark C, Sanchez C, Looman C, et al. Conserved interaction between distinct Kruppel-associated box domains and the transcriptional intermediary factor 1 beta. Proc Natl Acad Sci U S A. 2001;98:1422–1426. [PubMed]

17. Friedman JR, Fredericks WJ, Jensen DE, Speicher DW, Huang XP, et al. KAP-1, a novel corepressor for the highly conserved KRAB repression domain. Genes Dev. 1996;10:2067–2078. [PubMed]

18. Le Douarin B, Nielsen AL, Garnier JM, Ichinose H, Jeanmougin F, et al. A possible involvement of TIF1 alpha and TIF1 beta in the epigenetic control of transcription by nuclear receptors. Embo J. 1996;15:6701–6715. [PubMed]

19. Moosmann P, Georgiev O, Le Douarin B, Bourquin JP, Schaffner W. Transcriptional repression by RING finger protein TIF1 beta that interacts with the KRAB repressor domain of KOX1. Nucleic Acids Res. 1996;24:4859–4867. [PMC free article] [PubMed]

20. Lechner MS, Begg GE, Speicher DW, Rauscher FJ., 3rd Molecular determinants for targeting heterochromatin protein 1-mediated gene silencing: direct chromoshadow domain-KAP-1 corepressor interaction is essential. Mol Cell Biol. 2000;20:6449–6465. [PMC free article] [PubMed]

21. Nielsen AL, Ortiz JA, You J, Oulad-Abdelghani M, Khechumian R, et al. Interaction with members of the heterochromatin protein 1 (HP1) family and histone deacetylation are differentially involved in transcriptional silencing by members of the TIF1 family. Embo J. 1999;18:6385–6395. [PubMed]

22. Ryan RF, Schultz DC, Ayyanathan K, Singh PB, Friedman JR, et al. KAP-1 corepressor protein interacts and colocalizes with heterochromatic and euchromatic HP1 proteins: a potential role for Kruppel-associated box-zinc finger proteins in heterochromatin-mediated gene silencing. Mol Cell Biol. 1999;19:4366–4378. [PMC free article] [PubMed]

23. Schultz DC, Friedman JR, Rauscher FJ., 3rd Targeting histone deacetylase complexes via KRAB-zinc finger proteins: the PHD and bromodomains of KAP-1 form a cooperative unit that recruits a novel isoform of the Mi-2alpha subunit of NuRD. Genes Dev. 2001;15:428–443. [PubMed]

24. Schultz DC, Ayyanathan K, Negorev D, Maul GG, Rauscher FJ., 3rd SETDB1: a novel KAP-1-associated histone H3, lysine 9-specific methyltransferase that contributes to HP1-mediated silencing of euchromatic genes by KRAB zinc-finger proteins. Genes Dev. 2002;16:919–932. [PubMed]

25. Sripathy SP, Stevens J, Schultz DC. The KAP1 corepressor functions to coordinate the assembly of de novo HP1-demarcated microenvironments of heterochromatin required for KRAB zinc finger protein-mediated transcriptional repression. Mol Cell Biol. 2006;26:8623–8638. [PMC free article] [PubMed]

26. Underhill C, Qutob MS, Yee SP, Torchia J. A novel nuclear receptor corepressor complex, N-CoR, contains components of the mammalian SWI/SNF complex and the corepressor KAP-1. J Biol Chem. 2000;275:40463–40470. [PubMed]

27. Wiznerowicz M, Jakobsson J, Szulc J, Liao S, Quazzola A, et al. The Kruppel-associated box repressor domain can trigger de novo promoter methylation during mouse early embryogenesis. J Biol Chem. 2007;282:34535–34541. [PubMed]

28. Dialynas GK, Vitalini MW, Wallrath LL. Linking Heterochromatin Protein 1 (HP1) to cancer progression. Mutat Res. 2008;647:13–20. [PMC free article] [PubMed]

29. Nielsen AL, Oulad-Abdelghani M, Ortiz JA, Remboutsika E, Chambon P, et al. Heterochromatin formation in mammalian cells: interaction between histones and HP1 proteins. Mol Cell. 2001;7:729–739. [PubMed]

30. Platero JS, Hartnett T, Eissenberg JC. Functional analysis of the chromo domain of HP1. Embo J. 1995;14:3977–3986. [PubMed]

31. Thiru A, Nietlispach D, Mott HR, Okuwaki M, Lyon D, et al. Structural basis of HP1/PXVXL motif peptide interactions and HP1 localisation to heterochromatin. Embo J. 2004;23:489–499. [PubMed]

32. Kwon SH, Workman JL. The heterochromatin protein 1 (HP1) family: put away a bias toward HP1. Mol Cells. 2008;26:217–227. [PubMed]

33. Aagaard L, Laible G, Selenko P, Schmid M, Dorn R, et al. Functional mammalian homologues of the Drosophila PEV-modifier Su(var)3-9 encode centromere-associated proteins which complex with the heterochromatin component M31. Embo J. 1999;18:1923–1938. [PubMed]

34. Bannister AJ, Zegerman P, Partridge JF, Miska EA, Thomas JO, et al. Selective recognition of methylated lysine 9 on histone H3 by the HP1 chromo domain. Nature. 2001;410:120–124. [PubMed]

35. Lachner M, O'Carroll D, Rea S, Mechtler K, Jenuwein T. Methylation of histone H3 lysine 9 creates a binding site for HP1 proteins. Nature. 2001;410:116–120. [PubMed]

36. Nakayama J, Rice JC, Strahl BD, Allis CD, Grewal SI. Role of histone H3 lysine 9 methylation in epigenetic control of heterochromatin assembly. Science. 2001;292:110–113. [PubMed]

37. Ayyanathan K, Lechner MS, Bell P, Maul GG, Schultz DC, et al. Regulated recruitment of HP1 to a euchromatic gene induces mitotically heritable, epigenetic gene silencing: a mammalian cell culture model of gene variegation. Genes Dev. 2003;17:1855–1869. [PubMed]

38. O'Geen H, Squazzo SL, Iyengar S, Blahnik K, Rinn JL, et al. Genome-wide analysis of KAP1 binding suggests autoregulation of KRAB-ZNFs. PLoS Genet. 2007;3:e89. doi: 10.1371/journal.pgen.0030089. [PubMed]

39. Vogel MJ, Guelen L, de Wit E, Peric-Hupkes D, Loden M, et al. Human heterochromatin proteins form large domains containing KRAB-ZNF genes. Genome Res. 2006;16:1493–1504. [PubMed]

40. Deuschle U, Meyer WK, Thiesen HJ. Tetracycline-reversible silencing of eukaryotic promoters. Mol Cell Biol. 1995;15:1907–1914. [PMC free article] [PubMed]

41. Moosmann P, Georgiev O, Thiesen HJ, Hagmann M, Schaffner W. Silencing of RNA polymerases II and III-dependent transcription by the KRAB protein domain of KOX1, a Kruppel-type zinc finger factor. Biol Chem. 1997;378:669–677. [PubMed]

42. Mitchell RS, Beitzel BF, Schroder AR, Shinn P, Chen H, et al. Retroviral DNA integration: ASLV, HIV, and MLV show distinct target site preferences. PLoS Biol. 2004;2:e234. doi: 10.1371/journal.pbio.0020234. [PMC free article] [PubMed]

43. Wu X, Li Y, Crise B, Burgess SM. Transcription start regions in the human genome are favored targets for MLV integration. Science. 2003;300:1749–1751. [PubMed]

44. Schroder AR, Shinn P, Chen H, Berry C, Ecker JR, et al. HIV-1 integration in the human genome favors active genes and local hotspots. Cell. 2002;110:521–529. [PubMed]

45. Wang GP, Ciuffi A, Leipzig J, Berry CC, Bushman FD. HIV integration site selection: analysis by massively parallel pyrosequencing reveals association with epigenetic modifications. Genome Res. 2007;17:1186–1194. [PubMed]

46. Brasher SV, Smith BO, Fogh RH, Nietlispach D, Thiru A, et al. The structure of mouse HP1 suggests a unique mode of single peptide recognition by the shadow chromo domain dimer. Embo J. 2000;19:1587–1597. [PubMed]

47. Murzina N, Verreault A, Laue E, Stillman B. Heterochromatin dynamics in mouse cells: interaction between chromatin assembly factor 1 and HP1 proteins. Mol Cell. 1999;4:529–540. [PubMed]

48. Cammas F, Janoshazi A, Lerouge T, Losson R. Dynamic and selective interactions of the transcriptional corepressor TIF1 beta with the heterochromatin protein HP1 isotypes during cell differentiation. Differentiation. 2007;75:627–637. [PubMed]

49. Frietze S, Lan X, Jin VX, Farnham PJ. Genomic targets of the KRAB and scan domain-containing zinc finger protein 263 (ZNF263). J Biol Chem 2009 [PubMed]

50. Grewal SI, Jia S. Heterochromatin revisited. Nat Rev Genet. 2007;8:35–46. [PubMed]

51. Stewart MD, Li J, Wong J. Relationship between histone H3 lysine 9 methylation, transcription repression, and heterochromatin protein 1 recruitment. Mol Cell Biol. 2005;25:2525–2538. [PMC free article] [PubMed]

52. Riclet R, Chendeb M, Vonesch JL, Koczan D, Thiesen HJ, et al. Disruption of the interaction between transcriptional intermediary factor 1{beta} and heterochromatin protein 1 leads to a switch from DNA hyper- to hypomethylation and H3K9 to H3K27 trimethylation on the MEST promoter correlating with gene reactivation. Mol Biol Cell. 2009;20:296–305. [PMC free article] [PubMed]

53. Talbert PB, Henikoff S. Spreading of silent chromatin: inaction at a distance. Nat Rev Genet. 2006;7:793–803. [PubMed]

54. de Wit E, Greil F, van Steensel B. High-resolution mapping reveals links of HP1 with active and inactive chromatin components. PLoS Genet. 2007;3:e38. doi: 10.1371/journal.pgen.0030038. [PubMed]

55. Johansson AM, Stenberg P, Pettersson F, Larsson J. POF and HP1 bind expressed exons, suggesting a balancing mechanism for gene regulation. PLoS Genet. 2007;3:e209. doi: 10.1371/journal.pgen.0030209. [PubMed]

56. Burke SJ, Collier JJ, Scott DK. cAMP prevents glucose-mediated modifications of histone H3 and recruitment of the RNA polymerase II holoenzyme to the L-PK gene promoter. J Mol Biol. 2009;392:578–588. [PMC free article] [PubMed]

57. Cloos PA, Christensen J, Agger K, Helin K. Erasing the methyl mark: histone demethylases at the center of cellular differentiation and disease. Genes Dev. 2008;22:1115–1140. [PubMed]

58. Loh YH, Zhang W, Chen X, George J, Ng HH. Jmjd1a and Jmjd2c histone H3 Lys 9 demethylases regulate self-renewal in embryonic stem cells. Genes Dev. 2007;21:2545–2557. [PubMed]

59. Lomvardas S, Thanos D. Nucleosome sliding via TBP DNA binding in vivo. Cell. 2001;106:685–696. [PubMed]

60. Antoniou M, Harland L, Mustoe T, Williams S, Holdstock J, et al. Transgenes encompassing dual-promoter CpG islands from the human TBP and HNRPA2B1 loci are resistant to heterochromatin-mediated silencing. Genomics. 2003;82:269–279. [PubMed]

61. Lindahl Allen M, Antoniou M. Correlation of DNA methylation with histone modifications across the HNRPA2B1-CBX3 ubiquitously-acting chromatin open element (UCOE). Epigenetics. 2007;2:227–236. [PubMed]

62. Clark SJ. Action at a distance: epigenetic silencing of large chromosomal regions in carcinogenesis. Hum Mol Genet. 2007;16 Spec No 1:R88–95. [PubMed]

63. Beer DG, Kardia SL, Huang CC, Giordano TJ, Levin AM, et al. Gene-expression profiles predict survival of patients with lung adenocarcinoma. Nat Med. 2002;8:816–824. [PubMed]

64. Bhattacharjee A, Richards WG, Staunton J, Li C, Monti S, et al. Classification of human lung carcinomas by mRNA expression profiling reveals distinct adenocarcinoma subclasses. Proc Natl Acad Sci U S A. 2001;98:13790–13795. [PubMed]

65. De Palma M, Montini E, Santoni de Sio FR, Benedicenti F, Gentile A, et al. Promoter trapping reveals significant differences in integration site selection between MLV and HIV vectors in primary hematopoietic cells. Blood. 2005;105:2307–2315. [PubMed]

66. Wiznerowicz M, Trono D. Conditional suppression of cellular genes: lentivirus vector-mediated drug-inducible RNA interference. J Virol. 2003;77:8957–8961. [PMC free article] [PubMed]

67. Zufferey R, Nagy D, Mandel RJ, Naldini L, Trono D. Multiply attenuated lentiviral vector achieves efficient gene delivery in vivo. Nat Biotechnol. 1997;15:871–875. [PubMed]

68. Chen WV, Soriano P. Gene trap mutagenesis in embryonic stem cells. Methods Enzymol. 2003;365:367–386. [PubMed]

69. Ewing B, Hillier L, Wendl MC, Green P. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998;8:175–185. [PubMed]

70. Kent WJ. BLAT–the BLAST-like alignment tool. Genome Res. 2002;12:656–664. [PubMed]

71. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:RESEARCH0034. [PMC free article] [PubMed]

72. Iseli C, Ambrosini G, Bucher P, Jongeneel CV. Indexing strategies for rapid searches of short words in genome sequences. PLoS One. 2007;2:e579. doi: 10.1371/journal.pone.0000579. [PMC free article] [PubMed]

73. Li H, Rauch T, Chen ZX, Szabo PE, Riggs AD, et al. The histone methyltransferase SETDB1 and the DNA methyltransferase DNMT3A interact directly and localize to promoters silenced in cancer cells. J Biol Chem. 2006;281:19489–19500. [PubMed]

Articles from PLoS Genetics are provided here courtesy of
Public Library of Science