Open Access

H3K9me3-binding proteins are dispensable for SETDB1/H3K9me3-dependent retroviral silencing

  • Irina A Maksakova1,
  • Preeti Goyal1,
  • Jörn Bullwinkel3,
  • Jeremy P Brown3,
  • Misha Bilenky4,
  • Dixie L Mager1, 2,
  • Prim B Singh3 and
  • Matthew C Lorincz1Email author
Epigenetics & Chromatin20114:12

DOI: 10.1186/1756-8935-4-12

Received: 21 April 2011

Accepted: 20 July 2011

Published: 20 July 2011



Endogenous retroviruses (ERVs) are parasitic sequences whose derepression is associated with cancer and genomic instability. Many ERV families are silenced in mouse embryonic stem cells (mESCs) via SETDB1-deposited trimethylated lysine 9 of histone 3 (H3K9me3), but the mechanism of H3K9me3-dependent repression remains unknown. Multiple proteins, including members of the heterochromatin protein 1 (HP1) family, bind H3K9me2/3 and are involved in transcriptional silencing in model organisms. In this work, we address the role of such H3K9me2/3 "readers" in the silencing of ERVs in mESCs.


We demonstrate that despite the reported function of HP1 proteins in H3K9me-dependent gene repression and the critical role of H3K9me3 in transcriptional silencing of class I and class II ERVs, the depletion of HP1α, HP1β and HP1γ, alone or in combination, is not sufficient for derepression of these elements in mESCs. While loss of HP1α or HP1β leads to modest defects in DNA methylation of ERVs or spreading of H4K20me3 into flanking genomic sequence, respectively, neither protein affects H3K9me3 or H4K20me3 in ERV bodies. Furthermore, using novel ERV reporter constructs targeted to a specific genomic site, we demonstrate that, relative to Setdb1, knockdown of the remaining known H3K9me3 readers expressed in mESCs, including Cdyl, Cdyl2, Cbx2, Cbx7, Mpp8, Uhrf1 and Jarid1a-c, leads to only modest proviral reactivation.


Taken together, these results reveal that each of the known H3K9me3-binding proteins is dispensable for SETDB1-mediated ERV silencing. We speculate that H3K9me3 might maintain ERVs in a silent state in mESCs by directly inhibiting deposition of active covalent histone marks.


endogenous retrovirus ERV heterochromatin protein 1 HP1 Cbx1 Cbx3 Cbx5 H3K9me3 retroviral repression transcriptional silencing mouse embryonic stem cells


Endogenous retroviral sequences (ERVs) are relics of ancient retroviral integration into the germline. These parasitic elements are abundant in mammals, occupying approximately 8% of the mouse genome and 10% of the human genome [1, 2]. ERVs are subdivided into three diverse classes based on the similarity of their reverse transcriptase genes or their relationship to different genera of exogenous retroviruses. In the mouse, class I ERVs, similar to gammaretroviruses, include active families such as murine leukaemia viruses (MLVs) and murine retroviruses that use tRNAGln (GLN). Class II ERVs are similar to alpha- and betaretroviruses and include Mus musculus ERV using tRNALys type 10C (MMERVK10C), the highly retrotranspositionally active intracisternal A-type particles (IAPEz) and early transposon/Mus musculus type D retrovirus (ETn/MusD) families. Class III ERVs, the oldest and most abundant ERVs, are most similar to spumaviruses and are represented by mouse endogenous retrovirus type L (MERV-L) and mouse apparent LTR retrotransposons (MaLR) [3, 4]. Numerous regulatory motifs in the ERV long terminal repeats (LTRs) can initiate high levels of transcription in tissues and cell lines [5], and there is extensive evidence of aberrant ERV-driven gene expression in cancers [611] and tissues of aging mice [12, 13]. In an effort to counteract the potentially detrimental effects of ERVs, eukaryotic genomes have evolved multiple lines of defence against active exogenous and endogenous retroviruses [14], including DNA methylation and repressive histone modifications.

DNA methylation was the first epigenetic mark recognized to contribute to ERV silencing, with dramatic upregulation of ERVs observed in DNA methylation-deficient somatic cells [15, 16]. However, genome-wide chromatin immunoprecipitation (ChIP) followed by ChIP sequencing (ChIP-seq) [1719] or ChIP followed by quantitative PCR (qPCR) [20] revealed that in mouse embryonic stem cells (mESCs), class I and class II ERVs are enriched for the repressive histone H3 lysine 9 trimethylation (H3K9me3) deposited by lysine methyltransferase (KMTase) SETDB1/ESET/KMT1E [20]. SETDB1 is in turn thought to be recruited to ERVs via the obligatory corepressor KRAB-associated protein 1 (KAP-1) [21], presumably through sequence-specific KAP-1-binding zinc finger proteins such as ZFP809 in the case of MLVs [22]. Moreover, we and others have recently shown that in mESCs, H3K9me3 and SETDB1 play a greater role than DNA methylation in the silencing of class I and class II ERVs [20, 23]. IAP and ETn/MusD retrotransposons, the two most active class II mouse ERV families and the source of numerous recent germline mutations [24], are among the families with the highest H3K9me3 enrichment levels. Intriguingly, these families are dramatically upregulated in SETDB1 knockout (SETDB1 KO) mESCs [19, 20], confirming that they have a high potential for activation in the absence of H3K9me3. In contrast, the class III MERV-L and MaLR families, which are devoid of the H3K9me3 mark in mESCs, are repressed by the histone lysine-specific demethylase 1 (LSD1/KDM1A) [25], revealing that different ERV classes are regulated by distinct epigenetic modifications in these pluripotent cells.

Acetylation of lysine residues on the N-terminal tails of histones, including H3K9, directly influences the state of chromatin compaction by reducing the affinity of histones for DNA [26, 27]. In contrast, methylation per se of such lysine residues is less likely to directly affect chromatin structure, as this modification does not alter their charge. Rather, the prevailing view is that specific proteins, the so-called "readers," bind to methylated lysines and coordinate the biological outcome associated with such covalent histone marks. H3K9me3, for example, which is essential for the establishment and maintenance of the silent chromatin state [2831], is bound by three isoforms of heterochromatin protein 1 (HP1) in the mouse genome: HP1α (encoded by Cbx5), HP1β (encoded by Cbx1) and HP1γ (encoded by Cbx3) [32]. HP1 is a highly conserved family; its members are frequently present in several copies in eukaryotic genomes and play both structural and gene regulatory roles [3335]. The chromodomain of HP1 is responsible for binding H3K9me2/3 [36, 37], and a chromoshadow domain is required for HP1 homo- and heterodimerization and the recruitment of other proteins [38, 39].

Although their exact function in transcriptional regulation and cross-talk with histone and DNA methylation varies between species, the ability of HP1s to modulate gene expression via H3K9me2/3 binding has been reported in multiple systems [33, 4042]. In fission yeast, for example, two HP1 homologues, Swi6 and Chp2, are both required for assembly of repressive chromatin [43]. In mammalian cells, targeting of HP1α, HP1β and HP1γ to heterologous loci is sufficient to induce recruitment of SETDB1 and deposition of H3K9me3 [44], and HP1 has been implicated in SUV39H1-mediated silencing of euchromatic genes [45].

A role for HP1 proteins in silencing of repetitive and/or transposable elements has been well documented in several model organisms. In Drosophila, two families of transposons are derepressed in larvae with mutant HP1a and, to a lesser extent, mutant HP1c [46]. HP1d/Rhino is required for transposon silencing in the female germline of Drosophila, but this silencing seems to stem from Rhino's role in Piwi-interacting RNA (piRNA) production rather than establishment of repressive chromatin [47]. At transposable elements in Neurospora, DNA methylation is dependent on methylated H3K9 bound by HP1 [48, 49]. In Arabidopsis, however, H3K9me3-directed DNA methylation applies only to CpNpG methylation, not to CpG methylation, of transposons [50, 51]. HP1γ is a negative regulator of HIV in human cell lines [52] and of non-LTR LINE1 retrotransposons in male mouse germ cells [53]. On the contrary, HP1γ has also been implicated in activating gene expression through its association with elongating RNA polymerase II [54, 55]. The latter example notwithstanding, HP1 proteins are excellent candidates for the role of downstream effectors of H3K9me3-dependent silencing affecting ERVs in mESCs. Indeed, an intact HP1-binding domain of KAP-1 is essential for complete restriction of MLV in mouse embryonic carcinoma cells [56]. Furthermore, direct interaction of HP1 and KAP-1, as well as binding of HP1 to H3K9me3, is necessary for the full extent of silencing mediated by these factors [5761]. Moreover, we have recently demonstrated by ChIP-qPCR that HP1α, HP1β and HP1γ are enriched on IAPEz, MusD and MLV ERV sequences in mESCs, albeit at modest levels, and that this binding is partially dependent on SETDB1-deposited H3K9me3 [20]. On the basis of these observations, we hypothesized that HP1s might play a role in H3K9me3-mediated ERV silencing in mESCs and possibly in early embryos.

In addition to their reported roles in transcriptional silencing, HP1 proteins are required for heterochromatin spreading in specific genomic contexts in Drosophila[62, 63], yeast [64] and mammals [42, 57]. The presence of both chromodomains and chromoshadow domains suggests that HP1 proteins may bind H3K9me3 and recruit additional proteins, such as SUV39H1/2 or SETDB1-bound KAP-1 [61, 65, 66], to facilitate the spreading of the repressive H3K9me3 mark [67, 68]. Intriguingly, repetitive elements may act as foci of de novo heterochromatin formation and spreading, as H3K9me3 is enriched at sequences flanking ERVs [18, 19]. Conversely, in Neurospora, HP1 is a component of a histone demethylase-containing complex that prevents spreading of heterochromatin [69].

In addition to HP1s, many other mouse chromodomain proteins [70] are reported to bind H3K9me3 in vitro, including CDYL, CDYL2, CBX2, CBX4, CBX7 and M-phase phosphoprotein 8 (MPP8) [7178]. Furthermore, nonchromodomain proteins with affinity for H3K9me3 have also been identified [79]. Although MPP8 and CBX7 have been shown to negatively influence transcription of specific genes [71, 80], the functional and biological significance of the interaction of most of these H3K9me3 readers with H3K9me3 remains poorly understood.

To determine what role, if any, H3K9me3 readers play in silencing of ERVs and spreading of repressive chromatin from these repetitive elements, we first generated Cbx1 (HP1β) knock-out (KO) and Cbx5 (HP1α) KO mESCs [40, 81]. Surprisingly, we observed no upregulation of ERVs in Cbx5-/- mESCs and only modest upregulation of several ERV families in Cbx1-/- mESCs compared to that seen in Setdb1 KO mESCs. We found that both HP1α and HP1β are dispensable for DNA methylation of the ETnII/MusD family of ERVs, although HP1α has a modest influence on DNA methylation of IAP elements. Furthermore, we demonstrate that while deposition of H4K20me3 at major satellite repeats is dependent in part on HP1α, as reported previously [82], HP1α and HP1β are dispensable for deposition of H4K20me3 at ERVs and play only a modest role in spreading of H4K20me3 into sequences flanking these elements. Finally, employing RNAi and newly derived mESC lines harbouring silenced IAP, MusD and exogenous MLV-based reporters, we show that depletion of all of the HP1 proteins, alone or in combination, or each of the remaining known H3K9me3-binding proteins, has only a modest effect on ERV derepression, indicating that at classes I and II ERVs, H3K9me3 inhibits transcription independently of HP1 and other known H3K9me3 readers.


Catalytic activity of SETDB1 is largely required for ERV silencing

We recently showed by ChIP-qPCR [20] and ChIP-seq [19] analyses that numerous class I and class II ERV families are marked by H3K9me3. Furthermore, we demonstrated the critical role of SETDB1, the KMTase that deposits this mark, in transcriptional repression of these ERVs. Mapping all H3K9me3 ChIP-seq reads along the span of the consensus sequences of class I and class II ERVs, including IAPEz, MusD, MMERVK10C, MLV and GLN, confirms a high but nonuniform level of H3K9me3 along these elements in wt mESCs and a significantly lower level of H3K9me3 in Setdb1 KO mESCs (Figure 1A and Figure S3 in Additional file 1). Consistent with these data and those published in a previous report [18], analysis of the uniquely mapped ChIP-seq reads reveals a high level of H3K9me3 in the regions flanking IAPEz, MusD and MLV ERVs (Figure 1B and Figure S4 in Additional file 1).
Figure 1

Catalytically active SETDB1 is required for endogenous retrovirus silencing. (A) Profiling of trimethylated lysine 9 of histone 3 (H3K9me3) along the length of IAPEz endogenous retroviruses (ERVs) in the TT2 wild type (TT2 wt) and Setdb1 knockout (Setdb1 KO) mouse embryonic stem cells (mESCs) (see Figure S3 in Additional file 1 for profiles of murine leukaemia virus (MLV), MusD, MMERVK10C and GLN ERVs). The profile was generated by aligning chromatin immunoprecipitation assay sequencing (ChIP-seq) reads from TT2 wt and Setdb1 KO mESCs [19] to the consensus sequence of IAPEz. H3K9me3 enrichment levels are presented as reads per kilobase per million mapped reads values (RPKM). (B) Profiling of H3K9me3 and H4K20me3 in the genomic regions flanking 599 IAPEz elements in TT2 wt and Setdb1 KO mESCs (see Figure S4 in Additional file 1 for MusD and MLV profiles). H3K9me3 ChIP-seq reads from TT2 wt (C57BL/6 ± CBA) and Setdb1 KO mESCs [19] were used, along with H4K20me3 ChIP-seq from the wt V6.5 mESCs (129SvJae ± C57BL/6) [18]. Reads were aligned to the mouse genome (mm9), and the density of reads mapping to the 7-kb regions flanking intact IAPEz ERV families was plotted for H3K9me3 in TT2 wt and Setdb1 KO mESCs and for H4K20me3 in V6.5 wt mESCs. Vertical lines indicate the 5' and 3' boundaries of the ERV. The average mappability for 50-bp reads was confirmed to be, on average, uniform in the assayed 7 kb region (data not shown), ruling out the possibility of mapping bias. (C) Setdb1 deletion was induced with 4-hydroxytamoxifen (4-OHT) in mESCs containing no transgene (KO), a wt transgene (KO TG+) or a transgene with a mutation rendering SETDB1 catalytically inactive (KO C1243A) [20]. Expression is normalized to β-actin relative to wt. Data are presented as means ± standard deviations (SD) for three technical replicates. (D) To establish the expression levels of Setdb1 in the KO and transgenic lines, quantitative RT-PCR (qRT-PCR) was performed with Setdb1-specific primers, and expression was normalized to β-actin relative to wt. Data are presented as means ± SD for three technical replicates.

To confirm that the KMTase activity of SETDB1 is critical for ERV silencing in mESCs [20], we analyzed the Setdb1 conditional KO mESC line, either unmodified (SETDB1 KO) or stably expressing wild-type (wt) (SETDB1 KO TG+) or KMTase-defective (SETDB1 KO C1243A) SETDB1 transgenes, the latter harbouring a single amino acid change in the catalytic domain [20]. As expected, robust derepression of ERVs is observed in the SETDB1 KO line (Figure 1C). Despite the fact that the SETDB1 C1243A line expresses an approximately threefold higher level of Setdb1 than wt cells (Figure 1D and [20]), derepression of several of these ERVs is also observed in this transgenic line, confirming that SETDB1 KMTase activity is essential for ERV silencing. Interestingly, the extent of derepression was dependent on the ERV family. The level of upregulation of MusD and IAPEz elements was equivalent in the SETDB1 KO and catalytic mutant lines, suggesting that silencing of these elements depends on the KMTase activity of SETDB1. MMERVK10C and GLN show a lower level of derepression in the SETDB1 C1243A line than the SETDB1 KO line, and MLV remains completely restricted in the SETDB1 C1243A line. Similar results were noted previously in Northern blot analyses [20]. Taken together, these results indicate that different ERV families are subject to SETDB1-mediated silencing generally dependent on SETDB1 catalytic activity.

Depletion of HP1β but not HP1α leads to modest upregulation of select ERV families

Having confirmed that the KMTase activity of SETDB1 is required for efficient silencing of MMERVK10C, MusD and IAPEz, we next sought to determine whether the archetypal heterochromatic H3K9me2/3 readers HP1α and HP1β [40], both of which are enriched on IAPEz, MusD and MLV ERVs [20], are the effectors of transcriptional suppression of these elements. We generated Cbx5 (HP1α) KO mESCs (Figure S1 in Additional file 1) and Cbx1 (HP1β) KO mESCs (Figure S2 in Additional file 1) and confirmed downregulation of the corresponding genes at the mRNA level by qRT-PCR and at the protein level by Western blot analysis. Equivalent levels of expression of the pluripotency factor Nanog were detected in these lines, indicating that deletion of HP1 proteins does not stimulate differentiation (Figure 2A). Interestingly, while compensatory upregulation of the Cbx1 and Cbx3 genes was observed at the mRNA level in the Cbx5-/- line, upregulation of these genes was not observed at the protein level (Figure 2B).
Figure 2

Expression of heterochromatin protein 1 genes and ERVs in the Cbx5-/- and Cbx1-/- mESCs. (A) qRT-PCR with primers specific for Cbx5 (encoding HP1α) Cbx1 (encoding HP1β), Cbx3 (encoding HP1γ) and the pluripotency factor Nanog in the Cbx5-/- and Cbx1-/- mESC lines confirms the KOs and reveals compensatory upregulation of the genes encoding the remaining HP1 proteins in the Cbx5 -/- line. Expression levels were normalized to β-actin relative to wt, and the data are presented as means ± SD for three technical replicates. (B) Western blot analysis of whole-cell lysates confirms the lack of expression of HP1α and HP1β in the Cbx5-/- and Cbx1-/- mESC lines, respectively. (C) Expression of representative ERV families in the HM1 (wt), Cbx5-/- and Cbx1-/- mESCs was determined by qRT-PCR. Expression levels were normalized to β-actin relative to wt. Data are presented as means ± SD of four independent experiments, each of which was performed in triplicate. (D) Northern blot analysis of RNA isolated from the parental HM1,Cbx5-/- and Cbx1-/- mESC lines using probes specific for ETnII, MusD and IAP ERVs are shown. RNA from Dnmt1 -/- mESCs, in which IAP elements are upregulated approximately fourfold [20, 86] and MusD elements are upregulated approximately 1.5-fold [86], was used as a control.

Surprisingly, unlike deletion of Setdb1, deletion of Cbx5 does not lead to upregulation of any members of the ERV families analyzed, as determined by qRT-PCR (Figure 2C) or Northern blot analysis (Figure 2D). Similarly, deletion of Cbx1 has no effect on MusD, MLV or GLN elements. Although Cbx1 deletion does result in modest derepression of MMERVK10C (approximately 3-fold) and IAPEz (approximately 1.5-fold) relative to the parental HM1 line, these ERVs show approximately 47-fold and approximately 3-fold upregulation respectively, in the Setdb1 KO line, relative to the parental TT2 line, (see Figure 1). Taken together, these results indicate that in contrast to SETDB1, HP1α and HP1β play no role or a relatively minor role, respectively, in class II ERV silencing in mESCs.

Depletion of HP1α results in a modest reduction of DNA methylation at IAPEz ERVs

We recently demonstrated that while G9a is dispensable for silencing of ERVs, this H3K9 KMTase is required for efficient DNA methylation of these elements in mESCs [83]. Similarly, DNA methylation of major satellite repeats is dependent upon the H3K9 KMTase SUV39H1/2 in mESCs [84]. Intriguingly, HP1 proteins are required for DNA methylation of repetitive elements in Neurospora[48, 85], but the role of HP1 proteins in DNA methylation of ERVs in mESCs has not been explored. To address this question, ETnII/MusD and IAPEz families, shown previously to be densely DNA methylated in mESCs [20, 83, 86], were analyzed by bisulphite sequencing using genomic DNA isolated from wt, Cbx1-/- and Cbx5 -/- mESCs. In wt cells, several copies of ETnII and MusD were either completely unmethylated or hypomethylated specifically at the 5' end of the LTR (Figures 3A and 3C) as observed previously [86]. The number of methylated CpG sites per element of this family remained similar in either of the Cbx KO lines. In contrast, while the level of DNA methylation was very high at IAPEz elements in wt cells, several IAP molecules showed reduced levels of DNA methylation in the Cbx5-/- cell line (Figures 3B and 3C), indicating that HP1α plays a role in DNA methylation of a subset of IAP elements, presumably dependent upon their genomic location. Nevertheless, as discussed above, this modest decrease in DNA methylation did not result in derepression of IAP elements in these cells.
Figure 3

DNA methylation and chromatin marks at ERVs in Cbx5-/- and Cbx1-/-mESCs. (A) An approximately 600-bp fragment of the LTR and downstream region of ETnII/MusD ERVs was analyzed by bisulphite sequencing using primers that detect 105 ETnII/MusD elements, according to in silico PCR analysis (UCSC Genome Browser, http://​genome.​ucsc.​edu/​). (B) An approximately 400-bp fragment of the LTR and downstream region of IAP ERVs was analyzed by bisulphite sequencing using primers that detect 1,461 IAP elements. (C) Bisulphite-sequenced molecules were binned into four categories, depending on the number of methylated CpG sites detected, and the data are presented as the percentage of all clones for each cell line in each bin. While the ETnII/MusD family shows no difference in DNA methylation, several IAP molecules in the Cbx5 -/- cell line exhibit partial demethylation. (D) Native ChIP (N-ChIP) with antibodies against H3K9me3, H4K20me3 and pan-H3 was followed by qPCR using primers specific for major satellite repeats and IAP, MLV and MusD ERVs, and the data are presented as means ± SD for three technical replicates. The level of H4K20me3 was reduced by more than 50% at major satellite repeats in the Cbx5 -/- mESCs but remained at the same level at ERVs.

Neither HP1α nor HP1β are essential for H4K20me3 deposition at ERVs

Although H4K20me3 is dispensable for proviral silencing in mESCs, its deposition by SUV4-20H at ERVs requires SETDB1-deposited H3K9me3 [20]. On the basis of the fact that in mouse embryonic fibroblasts (MEFs), H4K20me3 at satellite repeats is dependent on SUV39H1/2-deposited H3K9me3 and subsequent binding of HP1 to this mark [82], we investigated whether H4K20me3 at ERVs is also dependent upon HP1 proteins in mESCs. Native ChIP (N-ChIP) followed by qPCR revealed that H4K20me3 enrichment was reduced by more than 50% at major satellite repeats in the Cbx5 -/- line (Figure 3D), demonstrating that as in MEFs [82], HP1 proteins are required for efficient H4K20me3 deposition at pericentric heterochromatin in mESCs. However, this mark is not entirely lost in either of the KO lines, presumably due to partial redundancy of HP1 proteins at major satellites. In line with these findings, it was recently shown that HP1β is dispensable for H4K20me3 and H3K9me3 deposition and localization in heterochromatin of mouse neurons [81]. Similarly, H4K20me3 levels at IAPEz, ETnII/MusD and MLV ERVs in the Cbx1-/- and Cbx5-/- lines remained at levels similar to the wt parent line, demonstrating either that H4K20me3 is deposited independently of HP1 binding or that these proteins act redundantly to promote deposition of H4K20me3 at these elements. As expected, H3K9me3 also remained unaltered in the absence of HP1α or HP1β (Figure 3D).

HP1β plays a role in the spreading of H4K20me3 but not H3K9me3 from ERVs into flanking genomic regions

Intriguingly, while HP1 homologs play a positive role in heterochromatin spreading in Drosophila[62, 63] and mammals [42, 57], HP1 plays a critical role in inhibiting aberrant spreading of heterochromatin in Neurospora[69]. Genome-wide analysis of H3K9me3 in wt mESCs reveals high levels of H3K9me3 in the immediate flanks of ERVs, including IAPEz, MusD and MLV elements, with progressively lower levels of this mark at distances farther from the ERV integration site [18, 19]. As expected, deposition of H3K9me3 in these regions is SETDB1-dependent [19] (see Figure 1B and Figure S4 in Additional file 1). Notably, the profile of H4K20me3 in the genomic regions flanking IAP, MusD and MLV elements is similar to that of H3K9me3 and the relative levels of both marks are consistent with their abundance in each ERV family (that is, IAP > MusD > MLV [20]). To determine whether spreading of H3K9me3 and/or H4K20me3 is affected in Cbx1-/- and/or Cbx5-/- mESCs, we examined these marks at the flanks of two randomly chosen full-length IAPEz elements and three genomic locations distal to the integration sites of these ERVs by ChIP-qPCR (Figure 4). IAPEz elements were chosen because, among the ERVs analyzed, on average, this family showed the highest mean H3K9me3 density in flanking genomic regions (Figure 1B and Figure S4 in Additional file 1). As expected, in wt cells, the levels of both H3K9me3 and H4K20me3 generally declined as the distance from the IAP increased, dropping substantially at approximately 3.5 kb. Depletion of either HP1 protein did not show a consistent effect on the spreading of H3K9me3 into the flanks of the selected IAP elements, since at the majority of regions surveyed enrichment was not statistically significantly different in each of the KO lines from the wt control. HP1β may be involved in propagation of H3K9me3 beyond 2 kb from the IAP assayed on chromosome 2, however, suggesting that at least at some loci, HP1β may facilitate the spreading of H3K9me3.
Figure 4

HP1β plays a role in H4K20me3 but not H3K9me3 spreading from ERVs into flanking genomic DNA. N-ChIP was performed with H3K9me3-, H4K20me3- and pan-H3-specific antibodies using chromatin isolated from HM1, Cbx5 -/- and Cbx1 -/- mESCs. The level of enrichment of these modifications at the flanks of two full-length IAP elements on chromosomes 2 and 5 as well as at positions approximately 1 kb, 2 kb and 3.5 kb distal to these flanking regions, was determined by qPCR. Data are presented as means ± SD of three technical replicates, and pairs of control and experimental samples with *P < 0.05 and **P < 0.01 (two-tailed Student's t-test) are shown. H3K9me3 enrichment levels across these genomic regions as determined using our previously published ChIP-seq data sets [19] are also shown for wt and Setdb1 KO mESCs.

Analysis of H4K20me3 in the same regions revealed no decrease in this mark in the Cbx5-/- line. In contrast, relative to the HM1 parent line, the Cbx1-/- line showed a consistent, approximately 1.5- to 2-fold decrease (P < 0.05, two-tailed Student's t- test) in H4K20me3 at both loci in distal regions 2 and 3 (Figure 4). Thus, while neither HP1 protein is required for deposition of H4K20me3 at the ERVs themselves (see Figure 3D), HP1β may generally be involved in the spreading of this covalent mark into the genomic regions flanking these repetitive elements.

Application of novel ERV reporter lines in a siRNA-based screen of H3K9me3-binding proteins

In addition to HP1 proteins, a number of other chromodomain proteins, including CDYL2, CBX2, CBX4, CBX7 and MPP8, as well as the Tudor domain-containing protein TDRD7, were recently shown to bind H3K9me3 in vitro[7173, 75, 76]. To address whether any of the H3K9me3 readers expressed in mESCs (all of those mentioned above with the exception of Cbx4 and Tdrd7) play a role in SETDB1-dependent silencing, we used recombinase-mediated cassette exchange (RMCE) [87, 88] (Figure 5A) to derive novel mESC lines with single-copy proviral reporters integrated at a specific genomic site. Specifically, constructs harbouring the green fluorescent protein (GFP) gene downstream of the MusD or IAP LTR promoters were generated and introduced into the same genomic site in the mESC line HA36 (a gift from F Lienert and D Schübeler) via RMCE. In parallel, the MFG-GFP construct [89] derived from the Moloney murine leukaemia virus (MMLV) and efficiently silenced in mESCs and embryonic carcinoma cells [20, 9093], and a cytomegalovirus (CMV)-GFP cassette were introduced into the same site. Following Cre-mediated recombination and a five-day negative selection with ganciclovir to exclude cells harbouring the original hygromycin B-herpes simplex virus thymidine kinase fusion (HyTK) cassette, each of the LTR reporters became silenced, while the CMV promoter maintained expression (data not shown). To select cells that contain the ERV-driven GFP gene silenced via the SETDB1 pathway, we transiently transfected the GFP-negative ganciclovir-resistant pools with siRNA specific for Setdb1. Depending on the cassette, GFP expression was induced in 20% to 65% of viable cells and these GFP-positive cells were isolated by fluorescence-activated cell sorting (FACS). The LTR reporter cassettes were progressively resilenced over approximately three weeks in culture (Figure 5B), and the negative populations were sorted at day 12 to be used as reporters. Importantly, GFP was efficiently reactivated in each population upon subsequent treatment of these pools with Setdb1 siRNA (Figure 5C), confirming that silencing of these LTR reporters is SETDB1-dependent at this integration site.
Figure 5

Silencing kinetics and reactivation of ERV reporters integrated in a specific genomic site. (A) Scheme for targeting of ERV reporters into a specific genomic site in mESCs via recombinase-mediated cassette exchange (RMCE). The mESC line HA36 contains a hygromycin B and herpes simplex virus thymidine kinase (HyTK) cassette between inverted Lox sites (L1 and 1L). MFG-green fluorescent protein (GFP), MusD-GFP and IAP-GFP proviral reporter cassettes, which contain the Moloney murine leukaemia virus, MusD (approximately +130 bp of downstream sequence) and IAP (approximately +450 bp of downstream sequence) LTRs, respectively, flanked by L1 and 1L sites, were cotransfected into the HA36 line with a Cre recombinase expression vector. Negative selection with ganciclovir eliminated cells with the original HyTK cassette, yielding pools of cells harbouring the proviral reporter cassettes predominantly integrated at the same site. (B) The kinetics of silencing of the MFG, MusD and IAP cassettes after reactivation of the RMCE pool with siRNA against Setdb1 are shown. (C) GFP-negative cells were sorted at day 12 postinduction with Setdb1 siRNA. Robust reactivation of GFP } expression from each of these pools of cells was observed upon secondary Setdb1 knockdown (KD). Flow cytometry data are presented as contour plots and histograms of 10,000 viable (propidium iodide (PI)-negative) cells.

To determine whether any of the remaining chromodomain-containing H3K9me3 readers expressed in mESCs are required for SETDB1-mediated silencing, we knocked down Cbx3 (HP1γ) as well as Cdyl2, Cbx1, Cbx2, Cbx5, Cbx7 and Mpp8 in the above-described reporter lines and a previously described pool of mESCs harbouring the silent murine stem cell virus (MSCV) provirus [20]. As expected, treatment of the MFG, MSCV, IAP and MusD reporter lines with Setdb1 and Kap1 specific siRNAs induced GFP expression in approximately 45% and approximately 25% of cells, respectively (Figure 6A, upper panel). In contrast, knockdown (KD) of each of the H3K9me3-binding proteins failed to induce GFP expression to the levels seen upon KD of Setdb1 or Kap1, despite efficient depletion of the target mRNAs (Figure 6A, lower panel). KD of genes encoding other chromodomain-encoding proteins with H3K9me-binding properties, such as Cdyl[72] and Chd4[94, 95], also did not result in reporter reactivation (Figure S5 in Additional file 1).
Figure 6

Reactivation of ERV reporters and ERVs upon siRNA-mediated KD of H3K9me3-binding proteins. (A) The percentage of enhanced green fluorescent protein-positive mESCs with reactivated ERV reporters was determined by flow cytometry (upper panel) on day 5 after the second transfection with siRNA against specified H3K9me3 readers. At least 10,000 cells were collected for each sample. Data are presented as means ± SD of three biological replicates. KD efficiency was determined by qRT-PCR at 30 hours after the second siRNA transfection (lower panel). Data are presented as means ± SD of three technical replicates. (B) Relative expression of ERVs at day 5 after the second KD with the indicated siRNA pool (upper panel), along with the efficiency of each KD, as determined by qRT-PCR at 30 hours after the second siRNA transfection (lower panel) is presented as means ± SD of three technical replicates. For each amplicon, expression was normalized to β-actin relative to scramble siRNA KD. (C) Western blot analysis of HP1 proteins in single- and triple-KD cells at day 5 after the second siRNA transfection is shown. H3 was used as a loading control.

KD of Cbx3 and Mpp8 did induce GFP expression in about 10% of treated cells, raising the possibility that these H3K9me3 readers act in a redundant manner to maintain these ERV reporters in a silent state. However, simultaneous KD of Cbx3 in combination with Mpp8 (Figure S6 in Additional file 1) or of Cbx3 in combination with Cbx1 and Cbx5 (Cbx1/3/5) (Figure 6A) did not significantly increase the percentage of GFP-positive cells over that observed with individual KD, despite efficient depletion of each mRNA. KD of Cbx1 and Cbx3 in the Cbx5 -/- mESCs and KD of Cbx3 and Cbx5 in the Cbx1 -/- mESCs showed similar results (data not shown). Thus, none of the assayed chromodomain-encoding proteins with H3K9-binding activity are essential for proviral silencing.

The H3K4me2/3 demethylase JARID1C (SMCX, which does not harbour a chromodomain, is capable of binding H3K9me3 via its plant homeodomain (PHD) [96], and its yeast homologue, Lid2, interacts directly with the H3K9 HMTase Clr4 [97]. These interactions suggest that JARID1C may direct H3K4 demethylation to loci marked by H3K9me3, promoting silencing. However, KD of the Jarid1 genes expressed in mESCs, including Jarid1a, Jarid1b and Jarid1c, either alone or in combination, leads to only modest reactivation of the proviral reporters, indicating that H3K9me3-recognizing H3K4 demethylases are not critical for maintenance of ERV silencing (Figure S7 in Additional file 1). Similarly, KD of Uhrf1 (NP95 in mouse and ICBP90 in human), which was recently shown to bind H3K9me3 via its PHD or SRA (SET- and RING-associated) domain [79, 98100], and/or KD of a related gene, Uhrf2 yields minimal upregulation of the four ERV reporters (Figure S8 in Additional file 1). Based on its pericentric localization [98], the main function of ICBP90 may lie in replication of heterochromatin and transcriptional regulation of major satellites [101], which show SUV39H1/2-dependent H3K9me3. Taken together, these data reveal that none of the known H3K9me3 readers are essential for silencing of ERVs that are repressed by the SETDB1 pathway.

To determine whether Cbx3 and Mpp8, the H3K9me3 readers which showed the highest reactivation of the LTR reporters, are required for silencing of ERVs, we performed qRT-PCR on cDNA isolated from wt TT2 and siRNA-treated mESCs. In Setdb1 KD cells, the MMERVK10C and MusD families showed the highest level of derepression, as expected. The same families, however, are only modestly upregulated upon KD of Cbx3 or Mpp8 (Figure 6B, upper panel), despite reduction of the target mRNA to 9% to 30% of wt levels (Figure 6B, lower panel) and dramatic downregulation of Cbx3 at the protein level, as determined by Western blot analysis (Figure 6C).

Finally, to determine whether HP1 proteins act redundantly to silence ERVs, we performed simultaneous KD of Cbx1, Cbx3 and Cbx5. Strikingly we observed only modest reactivation of each of the ERVs analyzed (Figure 6B, upper panel). Similar levels of upregulation of MMERVK10C and IAPEz ERVs in the Cbx1-/- mESCs (3.0- and 1.5-fold, respectively) (see Figure 2C), the Cbx3 KD (2.3- and 1.8-fold, respectively) and the triple Cbx1/3/5 KD (2.6- and 2.1-fold, respectively) suggest that Cbx1 and Cbx3 account for most of the HP1-mediated silencing of these ERV families. However, MusD elements, which are upregulated approximately fourfold in the triple KD, were not upregulated in any of the KOs, suggesting that all three HP1 proteins must be depleted to generate the relatively modest level of derepression observed for this family. Although we cannot rule out the possibility that an insufficient level or duration of HP1 depletion upon KD is responsible for these negative results, Western blot analysis revealed almost complete loss of all HP1 proteins in cells simultaneously depleted of Cbx1, Cbx3 and Cbx5 (Figure 6C). On the basis of the lack of ERV upregulation upon simultaneous KD of all three HP1 isoforms, we postulate that redundancy in HP1 function might not be the major factor preventing broad ERV reactivation. Similarly, as the maximum level of ERV reactivation upon KD of the remaining H3K9me3 readers is considerably lower than that observed in Setdb1 KO or Setdb1 KD cells, we conclude that none of these H3K9me3 readers play a major role in SETDB1-mediated ERV silencing in mESCs.


Role of H3K9me3 and HP1 in silencing of ERVs

We and others have recently shown that the H3K9me3 KMTase SETDB1 is critical for silencing of class I and class II ERVs in mESCs [20, 21]. However, the mechanism by which H3K9me3 modification leads to their transcriptional repression is currently unclear. In the present study, we have shown that HP1α plays a modest role in maintaining DNA methylation of IAPEz ERVs, while HP1β plays a modest role in promoting the spreading of H4K20me3 into the regions flanking these elements. HP1β also contributes to silencing of select IAPEz and MMERVK10C elements, but has no effect on DNA methylation of the ERVs analyzed. However, individual depletion of HP1α, HP1β and all other candidate H3K9me3 readers does not result in upregulation of ERVs or ERV reporters to a level observed in Setdb1 KO or Setdb1 KD mESCs, indicating that these factors either play only a modest role in silencing or act redundantly in this process. Strikingly, robust proviral derepression was not observed, even after simultaneous depletion of all three HP1 proteins, ruling out the latter, at least for these readers. Nevertheless, we cannot exclude the possibility that an as yet unidentified H3K9me3-binding protein and/or functional redundancy between H3K9me readers other than HP1 proteins may be required for H3K9me3-mediated ERV repression.

H3K9me3-dependent, H3K9me3 reader-independent proviral silencing?

Consistent with our observation that HP1s do not play a major role in transcriptional silencing of ERVs in mESCs, tethering of HP1 proteins in Drosophila inactivates only a limited number of reporter lines [102, 103]. Indeed, H3K9me3-dependent silencing may occur through mechanisms independent of H3K9me3 readers, such as by preventing the binding of transcription factors essential for transcription and/or the recruitment of the RNA polymerase II complex itself.

Specifically, H3K9me3 may directly or indirectly inhibit the deposition of active covalent histone marks. Acetylation of H3 at lysine 9 (H3K9Ac), for example, which is incompatible with methylation at this residue, promotes recruitment of chromatin remodelers and binding of RNA polymerase II in promoter regions [104108], and the histone acetyltransferases GCN5 and PCAF, which acetylate H3K9 [109, 110], are required for expression of specific genes [110] and retroviral elements [111]. Furthermore, hyperacetylation of H3 and H4 occurs concomitantly with IAP upregulation in MEFs and early embryos deficient in lymphoid-specific helicase (LSH) [112], implicating histone acetylation in ERV transcriptional activity. Intriguingly, in Xenopus oocytes expressing human H3K9 KMTases and HP1, H3K9me3 mediates transcriptional repression independently of HP1 recruitment through a mechanism that involves histone deacetylation [59].

Alternatively (or in addition), H3K9me3 may block transcription by indirectly inhibiting phosphorylation at serine 10 (H3S10ph) in the proviral promoter region. Intriguingly, transcriptional activation of the mouse mammary tumour retrovirus is dependent on H3S10ph and hyperacetylation of H3, mediated by binding of the nuclear factor 1 (NF-1) transcription factor to the proviral LTR [113115]. Predicted NF-1-binding sites are also found in the LTRs of other ERVs, including ETn elements [116], implicating a broad role for H3S10ph in transcription of these elements. While experiments directly addressing whether H3K9me3 blocks phosphorylation of H3S10ph have not been conducted in mammalian cells, H3K9me3 severely inhibits H3S10ph mediated by the Ipl1/aurora kinase in yeast [117].

Finally, H3K4 di- and trimethylation, marks also associated with the promoter regions of transcriptionally active genes, may also be inhibited by the presence of H3K9me3. Indeed, the H3K4 methyltransferases ASH1L [118] and SET7 [119] are less efficient in depositing H3K4me on histones marked with H3K9me in human cell lines, and H3K9me3 and H3K4me3 are mutually exclusive marks in mESCs [18]. While H3K4me2 is detected in the promoter region of IAP elements in Lsh -/- MEFs concomitant with their upregulation [120], the appearance of such active marks may be a consequence of, rather than a prerequisite for, transcriptional activation. To directly address the role of H3K4 methylation in retroviral expression, we sought to determine whether KD of Wdr5, a subunit of MLL/SET1 H3K4 methyltransferase complexes, inhibits Setdb1 KD-induced activation of the ERV reporters. We found that simultaneous KD of Setdb1 and Wdr5, reduced the level of reactivation of all ERV reporters, especially MusD and IAP (Figure S9 in Additional file 1), indicating that H3K4me3, the catalytic product of WDR5-containing complexes [121], is indeed required for optimal transcription of ERVs. Thus, the presence of H3K9me3 may effectively block transcription by inhibiting deposition of H3K4me3 and/or the other active marks mentioned above.

Heterochromatin spreading into sequences flanking ERV

Heterochromatin spreading is thought to involve a reiterative process of HP1 proteins binding to H3K9me2/3 [36, 37] followed by the recruitment of protein complexes with H3K9me2/3 catalytic activity, such as SUV39H1/2 [65] and SETDB1 [61]. Consistent with this model, HP1 proteins have been implicated in heterochromatin spreading in Drosophila[62, 63], yeast [64] and mammals [42, 57]. Moreover, H3K9me3, a hallmark of silent chromatin, is abundant in the vicinity of ERVs [18, 19]. However, our results indicate that HP1α and HP1β play only a modest role, if any, in the spreading of H3K9me3 into the sequences flanking ERVs. In contrast, HP1β is required for efficient spreading of H4K20me3 at the IAP ERVs analyzed. Although the biological role of H4K20me3 spreading is still unclear, recent studies have indicated that this covalent mark is involved in the maintenance of genomic stability [122124]. Intriguingly, a role for HP1 in the DNA damage response independent of H3K9me3 has also been reported [125, 126]. The availability of HP1β KO embryos will allow for studies aimed at addressing whether the distribution of H4K20me3 is dependent upon this protein and whether DNA damage repair pathways are perturbed in vivo.


In this work, we demonstrate the surprising finding that despite the accepted function of HP1 proteins in H3K9me-dependent gene silencing and the critical role of H3K9me3 in transcriptional repression of class I and class II ERVs, HP1α and HP1β are not required for silencing of these repetitive elements. Furthermore, while neither HP1α nor HP1β is essential for DNA methylation or the deposition of H3K9me3 or H4K20me3 within ERVs, HP1β plays a role in the spreading of the latter into sequences flanking these elements. Using a RNAi-based screen with newly derived mESCs harbouring novel ERV reporters, we have shown that the remaining proteins reported to bind H3K9me3 in vitro, including HP1γ, CDYL, CDYL2, CBX2, CBX7, MPP8, UHRF1 and JARID1A-C, are also dispensable for ERV silencing. The lack of proviral derepression in these experiments may be explained by functional redundancy of these or as yet unidentified H3K9me3 readers. Alternatively, H3K9me3 may repress ERV transcription via inhibiting deposition of covalent histone modifications required for transcription. The ERV reporter cell lines generated here should be useful in future screens of factors predicted to play a role in proviral expression. Regardless, additional studies aimed at delineating the functional significance of H3K9 readers, including nuclear processes not directly related to transcription, are clearly warranted.

Materials and methods

Cell culture, constructs and recombinase-mediated cassette exchange

To produce the Cbx5 -/- and Cbx1 mESC lines, each Cbx allele was targeted sequentially using two different targeting vectors. See Figures S1 and S2 in Additional file 1 for a detailed description of the ESC-KO derivation. mESCs were cultured in DMEM supplemented with 15% fetal bovine serum (HyClone Laboratories, Logan, UT, USA), 20 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid, 0.1 mM nonessential amino acids, 0.1 mM 2-mercaptoethanol, 100 U/mL penicillin, 0.05 mM streptomycin, leukaemia-inhibitory factor and 2 mM L-glutamine on gelatinized plates. For RMCE into HA36 cells, CMV was cut out of the L1-CMV-GFP-1L vector [127] by restriction with Cla I and Nhe I restriction enzymes. IAP and MusD LTR, together with the downstream sequence, were cloned into the resulting Cla I-Nhe I site upstream of the enhanced green fluorescent protein (EGFP) gene. MusD from the C57BL/6 genomic DNA on chr8:131270355-131277831 (mm9) was cloned, an element similar in sequence to those commonly expressed in wt cells [20]. The Nhe I site at nt 444 prevented us from including a longer fragment. However, this sequence still included the 319 bp 5'-LTR and 125 bp immediately downstream of it. A fragment containing a LTR and a downstream sequence, approximately 800 bp in total, was cloned for an IAP reporter. The element chosen was the one at the site of a novel insertion into the A/WySn mouse strain [128, 129] and was cloned from the DNA of the respective strain. All inserts were confirmed by sequencing. The primer sequences are given in Table S1 in Additional file 1.

Recombinase-mediated cassette exchange, transfection and transgene selection

For targeting of the ERV reporter constructs into the genome, Cre RMCE was used [87, 130]. The HA36 mESC line contains a cassette with the HyTK fusion gene at the random integration site which allows CMV-GFP expression for multiple passages (cell line a gift from F Lienert and D Schübeler). This selectable marker allows for positive selection through resistance to hygromycin B and for negative selection through sensitivity to ganciclovir. HA36 mESCs were cultured in 25 μg/mL hygromycin B for 14 days before transfection to select for cells expressing the fusion gene. Cells were transfected with Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) in a 24-well plate according to the manufacturer's recommendations. Briefly, 1.5 μg of a cassette with a MusD, IAP or MFG insert was cotransfected with 0.5 μg of CMV-Cre plasmid using 2 μL of Lipofectamine 2000 per well. After three days, cells were transferred to medium containing 3 μM ganciclovir to select against cells still expressing the HyTK fusion gene. Cells were grown in ganciclovir-containing medium for five or more days, with subculturing performed when necessary.

siRNA-mediated knockdown

For reporter assays, 10,000 mESCs per well of a 96-well plate were seeded into antibiotic-free mESC medium the day before transfection. Transfection was performed according to the manufacturer's protocol using 100 nM concentrations of each siRNA (siGENOME SMARTpool reagent Dharmacon, Lafayette, CO, USA) and 0.4 μL of DharmaFECT 1 siRNA transfection reagent (Dharmacon) per well. On the first day after transfection, approximately 1/5th of the cells were transferred into another 96-well plate containing antibiotic-free mESC medium, and the KD was repeated on the third day. The next day, approximately 1/2 of the cells were transferred into a 24-well or 12-well plate, and flow cytometry was performed on day 4 or 5 after the second KD. For RNA or protein collection, the first KD was performed in a 12-well plate and the cells were transferred to two 6-cm dishes the next day. The day after the second KD in 6-cm dishes, three-fourths of the cells were collected for RNA for confirmation of KD efficiency, and the rest were plated onto two 10-cm dishes for expansion and collection for RNA or protein on day 4 after the second KD.

Preparation of genomic DNA, bisulphite treatment, T/A cloning and sequencing

Genomic DNA was extracted using DNAzol reagent (Invitrogen), and bisulphite conversion of DNA was performed using the EZ DNA Methylation Kit (Zymo Research, Orange, CA, USA) according to the manufacturer's protocol. The approximately 370 bp of IAP and approximately 590 bp of ETnII/MusD element sequence containing the LTR and the downstream region were amplified from converted DNA by PCR using Platinum Taq (Invitrogen). The primer sequences are given in Additional file 1, Table S1. PCR products from three separate PCRs for each sample were cloned using the pGEM-T Easy Vector System kit (Promega, Madison, WI, USA). All sequences had a conversion rate of > 98%. QUMA http://​quma.​cdb.​riken.​jp/​top/​index.​html, with some follow-up processing, was used for analysis of bisulphite data [131].

Native chromatin immunoprecipitation assay and quantitative PCR

Briefly, 1 × 107 mESCs } for each cell line were resuspended in douncing buffer and homogenized through a 25-gauge 5/8-inchneedle syringe for 20 repetitions. A quantity of 1.875 μL of 20 U/μL micrococcal nuclease (MNase; Worthington Biochemical Corp., Lakewood, NJ, USA) was added and incubated at 37°C for 7 minutes. The reaction was quenched with 0.5 M ethylenediaminetetraacetic acid and incubated on ice for 5 minutes; then 1 mL of hypotonic buffer was added and incubated on ice for 1 hour. Cellular debris was pelleted, and the supernatant was recovered. Protein A/G Sepharose beads were blocked with single-stranded salmon sperm DNA and BSA, washed and resuspended in immunoprecipitation buffer. Blocked protein A/G Sepharose beads were added to the digested chromatin fractions and rotated at 4°C for 2 hours to preclear chromatin. A quantity of 100 μL of the precleared chromatin was purified by phenol-chloroform extraction, and DNA fragment sizes were analyzed and confirmed to correspond to one to three nucleosome fragments. Chromatin was subdivided into aliquots for each immunoprecipitated sample. Antibodies specific for unmodified H3 (H9289; Sigma-Aldrich, St Louis, MO, USA), H3K9me3 (Active Motif 39161, Carlsbad, CA, USA), H4K20me3 (Active Motif 39180) and control immunoglobulin G (I8140; Sigma-Aldrich, St Louis, MO, USA) were added to each tube and rotated at 4°C for 1 hour. The antibody-protein-DNA complex was precipitated by adding 20 μL of the blocked protein A/G Sepharose beads and rotated at 4°C overnight. The complex was washed and eluted, and immunoprecipitated material was purified using the QIAquick PCR Purification Kit (Qiagen, Germantown, MD, USA). The purified DNA was analyzed by qPCR with respect to input using EvaGreen dye (Biotium, Hayward, CA, USA) and Maxima Hot Start Taq DNA Polymerase (Fermentas, Vilnius, Lithuania). Primers are listed in Table S1 in Additional file 1.

RNA isolation, reverse transcription and quantitative RT-PCR

RNA was isolated using GenElute™ Mammalian Total RNA Miniprep Kit (Sigma-Aldrich) and reverse-transcribed using SuperScript III Reverse Transcriptase (Invitrogen) as per the manufacturers' instructions. Quantitative RT-PCR was carried out using SsoFAST™ EvaGreen Supermix (Bio-Rad Laboratories, Hercules, CA, USA) on StepOne™ version 2.1 software (Applied Biosystems, Foster City, CA, USA) in a total volume of 20 μL. Data are presented as means ± standard deviations of three technical replicates. Primer efficiencies were around 100%. Dissociation curve analysis was performed after the end of the PCR to confirm the presence of a single and specific product.

Whole-cell protein extracts and Western blot analysis

Briefly, cells were resuspended in 2 ± Laemmli buffer and incubated at 100°C for 10 minutes. Cells were then homogenized through a 25-gauge needle syringe 10 to 15 times. Extracts were run on SDS-PAGE gels and transferred onto a membrane. Primary antibodies used were α-HP1α (05-684, 1:200 dilution; Upstate Biotechnology, Lake Placid, NY, USA), α-HP1β (MCA 1946, 1:100 dilution; AbD Serotec, Burlington, ON, Canada) and α-H3 (Active Motif 39163, 1:200 dilution, Carlsbad, CA, USA). Secondary antibodies used at 1:10,000 dilutions were IRdye 800CW (926-32210) and IRdye 680 (926-32221), both from LI-COR Biosciences (Lincoln, NE, USA). Membranes were analyzed using the Odyssey Infrared Imaging System (LI-COR Biosciences).

Northern blot analysis

For each lane, 6 mg of RNA were denatured, electrophoresed in 1% agarose/3.7% formaldehyde gel in 1 ± 3-(N-morpholino)propanesulfonic acid buffer, transferred overnight onto a Zeta-Probe nylon membrane (Bio-Rad Laboratories, Hercules, CA, USA) and baked at 80°C. ETnII/MusD-, IAP- and Gapdh-specific probes were synthesized by PCR. Primer sequences are given in Additional file 1, Table S1. An amplified DNA fragment was α-32P-labeled using the Random Primers DNA Labeling System (Invitrogen). Membranes were prehybridized in ExpressHyb hybridization solution (Clontech, Mountain View, CA, USA) for two to four hours at 68°C, hybridized overnight at the same temperature in fresh ExpressHyb solution, washed according to the manufacturer's instructions and exposed to film.

Fluorescence-activated cell sorting and analysis of cassette integration

FACS analysis was performed using BD FACSAria III cell sorter with BD FACSDiva software (BD Biosciences), and flow analysis was performed using a BD LSR II flow cytometer. Viable cells were gated on the basis of propidium iodide exclusion. At least 10,000 propidium iodide-negative events were analyzed. Untransfected cells were used as a control for baseline EGFP fluorescence.

H3K9me3 profiling of endogenous retroviruses

To determine the H3K9me3 status of ERVs in TT2 wt versus Setdb1 KO mESC lines, we generated average H3K9me3 profiles for representative ERVs upregulated in the latter [19], including MusD, MMERVK10C, IAPEz, MLV and GLN elements. For each ERV family, all sequenced 50-bp reads from our previously published TT2 and Setdb1 KO H3K9me3 native ChIP-seq data sets [19] were aligned to the corresponding consensus sequences (including internal regions and corresponding LTRs) from Repbase http://​www.​girinst.​org/​repbase/​[132] for all ERVs except MusD. For MusD, a representative active element was used (139824 to 132348 nt of AC084696, reverse strand). The Burrows-Wheeler Aligner http://​bio-bwa.​sourceforge.​net/​[133] was employed with default parameters (allowing up to two mismatches in the 32-bp seed and one gap). Reads were directionally extended by 150 bp, and extended reads were profiled along the element. All mapped reads were taken into account, and the profiles for each library were normalized by the total number of reads uniquely mapped to the mm9 genome. For reads that were aligned into multiple locations (LTRs), we considered only one randomly selected alignment location. The irregular nature of the profile is most likely attributable to SNPs and insertions and/or deletions in the consensus vs. genomic reads.

H3K9me3 profiling in the sequences flanking endogenous retroviruses

To compare the average density of H3K9me3 in the genomic regions flanking ERVs, H3K9me3 N-ChIP-seq data sets for TT2 wt and Setdb1 KO mESCs [19] were used. Intact elements were selected for three ERV families: MusD, IAPEz, and MLV. For MusD, IAPEz and MLV, 195, 599 and 51 elements, respectively, satisfied the length and sequence similarity criteria that we applied [19]. All H3K9me3 reads that passed the quality score threshold above 7 were aligned to the mouse genome (mm9) using the Burrows-Wheeler Aligner [133] and directionally extended by 150 bp [19]. Only reads uniquely aligned to the regions within 7 kb on either side of intact elements were considered. If multiple reads were mapped to the same location, only one copy of the read was counted. To generate the profiles shown for the TT2 wt and Setdb1 KO cell lines, extended reads were first agglomerated for 5'- and 3'-flanks. Subsequently, the data were normalized to the total number of included elements and weighted by the total number of aligned reads to the genome for each sample.

Note added in proof

While this manuscript was under review, an article published by Shang and colleagues (PNAS 2011, 108(18):7541-7546) revealed that the H3K9me3 demethylase JMJD2B greatly facilitates H3K4 methylation by purified MLL2 in vitro (demonstrating that H3K9 demethylation is required for efficient H3K4 methylation) and is required for transcription of MLL2 targets in vivo.



bovine serum albumin


Dulbecco's modified Eagle's medium


quantitative polymerase chain reaction


reverse transcriptase


small interfering RNA


single-nucleotide polymorphism.



We thank Danny Leung, Sandra Lee, Lucia Lam and the University of British Columbia flow cytometry facility for technical support. We also thank Florian Lienert and Dirk Schübeler for providing the HA36 ES cell line, En Li for providing the Dnmt1-KO line, Yoichi Shinkai for providing the Setdb1 KO line and Mark Bedford for helpful suggestions. This work was supported by CIHR grant 77805 (to ML) and CIHR grant 92090 (to ML and DM). This work was also supported by Biotechnology and Biological Sciences Research Council core strategic grants and Deutsche Forschungsgemeinschaft grant SI 1209/2-1 (to PS). ML is a Scholar of the MSFHR and a CIHR New Investigator.

Authors’ Affiliations

Department of Medical Genetics, Life Sciences Institute, University of British Columbia
Terry Fox Laboratory, BC Cancer Agency
Division of Immunoepigenetics, Department of Immunology and Cell Biology, Research Center Borstel
Canada's Michael Smith Genome Sciences Centre, BC Cancer Agency


  1. International Mouse Genome Sequencing Consortium: Initial sequencing and comparative analysis of the mouse genome. Nature. 2002, 420: 520-562. 10.1038/nature01262.View Article
  2. International Human Genome Sequencing Consortium: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.View Article
  3. Stocking C, Kozak C: Endogenous retroviruses. Cell Mol Life Sci. 2008, 65: 3383-3398. 10.1007/s00018-008-8497-0.PubMedView Article
  4. Gifford R, Kabat P, Martin J, Lynch C, Tristem M: Evolution and distribution of class II-related endogenous retroviruses. J Virol. 2005, 79: 6478-6486. 10.1128/JVI.79.10.6478-6486.2005.PubMed CentralPubMedView Article
  5. Gimenez J, Montgiraud C, Pichon JP, Bonnaud B, Arsac M, Ruel K, Bouton O, Mallet F: Custom human endogenous retroviruses dedicated microarray identifies self-induced HERV-W family elements reactivated in testicular cancer upon methylation control. Nucleic Acids Res. 2010, 38: 2229-2246. 10.1093/nar/gkp1214.PubMed CentralPubMedView Article
  6. Lamprecht B, Walter K, Kreher S, Kumar R, Hummel M, Lenze D, Köchert K, Bouhlel MA, Richter J, Soler E, Stadhouders R, Jöhrens K, Wurster KD, Callen DF, Harte MF, Giefing M, Barlow R, Stein H, Anagnostopoulos I, Janz M, Cockerill PN, Siebert R, Dörken B, Bonifer C, Mathas S: Derepression of an endogenous long terminal repeat activates the CSF1R proto-oncogene in human lymphoma. Nat Med. 2010, 16: 571-579. 10.1038/nm.2129.PubMedView Article
  7. Moyes D, Griffiths DJ, Venables PJ: Insertional polymorphisms: a new lease of life for endogenous retroviruses in human disease. Trends Genet. 2007, 23: 326-333. 10.1016/j.tig.2007.05.004.PubMedView Article
  8. McLaughlin-Drubin ME, Munger K: Viruses associated with human cancer. Biochim Biophys Acta. 2008, 1782: 127-150.PubMed CentralPubMedView Article
  9. Howard G, Eiges R, Gaudet F, Jaenisch R, Eden A: Activation and transposition of endogenous retroviral elements in hypomethylation induced tumors in mice. Oncogene. 2008, 27: 404-408. 10.1038/sj.onc.1210631.PubMedView Article
  10. Lee JS, Haruna T, Ishimoto A, Honjo T, Yanagawa S: Intracisternal type A particle-mediated activation of the Notch4/int3 gene in a mouse mammary tumor: generation of truncated Notch4/int3 mRNAs by retroviral splicing events. J Virol. 1999, 73: 5166-5171.PubMed CentralPubMed
  11. Romanish MT, Cohen CJ, Mager DL: Potential mechanisms of endogenous retroviral-mediated genomic instability in human cancer. Semin Cancer Biol. 2010, 20: 246-253. 10.1016/j.semcancer.2010.05.005.PubMedView Article
  12. Puech A, Dupressoir A, Loireau MP, Mattei MG, Heidmann T: Characterization of two age-induced intracisternal A-particle-related transcripts in the mouse liver: transcriptional read-through into an open reading frame with similarities to the yeast ccr4 transcription factor. J Biol Chem. 1997, 272: 5995-6003. 10.1074/jbc.272.9.5995.PubMedView Article
  13. Barbot W, Dupressoir A, Lazar V, Heidmann T: Epigenetic regulation of an IAP retrotransposon in the aging mouse: progressive demethylation and de-silencing of the element by its repetitive induction. Nucleic Acids Res. 2002, 30: 2365-2373. 10.1093/nar/30.11.2365.PubMed CentralPubMedView Article
  14. Goff SP: Retrovirus restriction factors. Mol Cell. 2004, 16: 849-859.PubMedView Article
  15. Yoder JA, Walsh CP, Bestor TH: Cytosine methylation and the ecology of intragenomic parasites. Trends Genet. 1997, 13: 335-340. 10.1016/S0168-9525(97)01181-5.PubMedView Article
  16. Walsh CP, Chaillet JR, Bestor TH: Transcription of IAP endogenous retroviruses is constrained by cytosine methylation. Nat Genet. 1998, 20: 116-117. 10.1038/2413.PubMedView Article
  17. Martens JH, O'Sullivan RJ, Braunschweig U, Opravil S, Radolf M, Steinlein P, Jenuwein T: The profile of repeat-associated histone lysine methylation states in the mouse epigenome. EMBO J. 2005, 24: 800-812. 10.1038/sj.emboj.7600545.PubMed CentralPubMedView Article
  18. Mikkelsen TS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, Alvarez P, Brockman W, Kim TK, Koche RP, Lee W, Mendenhall E, O'Donovan A, Presser A, Russ C, Xie X, Meissner A, Wernig M, Jaenisch R, Nusbaum C, Lander ES, Bernstein BE: Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature. 2007, 448: 553-560. 10.1038/nature06008.PubMed CentralPubMedView Article
  19. Karimi MM, Goyal P, Maksakova IA, Bilenky M, Leung D, Tang JX, Shinkai Y, Mager DL, Jones S, Hirst M, Lorincz MC: DNA methylation and SETDB1/H3K9me3 regulate predominantly distinct sets of genes, retroelements, and chimeric transcripts in mESCs. Cell Stem Cell. 2011, 8: 676-687. 10.1016/j.stem.2011.04.004.PubMedView Article
  20. Matsui T, Leung D, Miyashita H, Maksakova IA, Miyachi H, Kimura H, Tachibana M, Lorincz MC, Shinkai Y: Proviral silencing in embryonic stem cells requires the histone methyltransferase ESET. Nature. 2010, 464: 927-931. 10.1038/nature08858.PubMedView Article
  21. Rowe HM, Jakobsson J, Mesnard D, Rougemont J, Reynard S, Aktas T, Maillard PV, Layard-Liesching H, Verp S, Marquis J, Spitz F, Constam DB, Trono D: KAP1 controls endogenous retroviruses in embryonic stem cells. Nature. 2010, 463: 237-240. 10.1038/nature08674.PubMedView Article
  22. Wolf D, Goff SP: Embryonic stem cells use ZFP809 to silence retroviral DNAs. Nature. 2009, 458: 1201-1204. 10.1038/nature07844.PubMed CentralPubMedView Article
  23. Hutnick LK, Huang X, Loo TC, Ma Z, Fan G: Repression of retrotransposal elements in mouse embryonic stem cells is primarily mediated by a DNA methylation-independent mechanism. J Biol Chem. 2010, 285: 21082-21091. 10.1074/jbc.M110.125674.PubMed CentralPubMedView Article
  24. Maksakova IA, Romanish MT, Gagnier L, Dunn CA, van de Lagemaat LN, Mager DL: Retroviral elements and their hosts: insertional mutagenesis in the mouse germ line. PLoS Genet. 2006, 2: e2-10.1371/journal.pgen.0020002.PubMed CentralPubMedView Article
  25. Macfarlan TS, Gifford WD, Agarwal S, Driscoll S, Lettieri K, Wang J, Andrews SE, Franco L, Rosenfeld MG, Ren B, Pfaff SL: Endogenous retroviruses and neighboring genes are coordinately repressed by LSD1/KDM1A. Genes Dev. 2011, 25: 594-607. 10.1101/gad.2008511.PubMed CentralPubMedView Article
  26. Tse C, Sera T, Wolffe AP, Hansen JC: Disruption of higher-order folding by core histone acetylation dramatically enhances transcription of nucleosomal arrays by RNA polymerase III. Mol Cell Biol. 1998, 18: 4629-4638.PubMed CentralPubMedView Article
  27. Krajewski WA: Histone hyperacetylation facilitates chromatin remodelling in a Drosophila embryo cell-free system. Mol Gen Genet. 2000, 263: 38-47. 10.1007/PL00008674.PubMedView Article
  28. Strahl BD, Allis CD: The language of covalent histone modifications. Nature. 2000, 403: 41-45. 10.1038/47412.PubMedView Article
  29. Jenuwein T, Allis CD: Translating the histone code. Science. 2001, 293: 1074-1080. 10.1126/science.1063127.PubMedView Article
  30. Wysocka J: Identifying novel proteins recognizing histone modifications using peptide pull-down assay. Methods. 2006, 40: 339-343. 10.1016/j.ymeth.2006.05.028.PubMed CentralPubMedView Article
  31. Daniel JA, Pray-Grant MG, Grant PA: Effector proteins for methylated histones: an expanding family. Cell Cycle. 2005, 4: 919-926. 10.4161/cc.4.7.1824.PubMedView Article
  32. Zeng W, Ball AR, Yokomori K: HP1: heterochromatin binding proteins working the genome. Epigenetics. 2010, 5: 287-292. 10.4161/epi.5.4.11683.PubMed CentralPubMedView Article
  33. Vermaak D, Malik HS: Multiple roles for heterochromatin protein 1 genes in Drosophila. Annu Rev Genet. 2009, 43: 467-492. 10.1146/annurev-genet-102108-134802.PubMedView Article
  34. Kwon SH, Workman JL: The changing faces of HP1: from heterochromatin formation and gene silencing to euchromatic gene expression. HP1 acts as a positive regulator of transcription. Bioessays. 2011, 33: 280-289. 10.1002/bies.201000138.PubMedView Article
  35. Singh PB, Georgatos SD: HP1: facts, open questions, and speculation. J Struct Biol. 2002, 140: 10-16. 10.1016/S1047-8477(02)00536-1.PubMedView Article
  36. Bannister AJ, Zegerman P, Partridge JF, Miska EA, Thomas JO, Allshire RC, Kouzarides T: Selective recognition of methylated lysine 9 on histone H3 by the HP1 chromo domain. Nature. 2001, 410: 120-124. 10.1038/35065138.PubMedView Article
  37. 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. 10.1038/35065132.PubMedView Article
  38. Nielsen AL, Oulad-Abdelghani M, Ortiz JA, Remboutsika E, Chambon P, Losson R: Heterochromatin formation in mammalian cells: interaction between histones and HP1 proteins. Mol Cell. 2001, 7: 729-739. 10.1016/S1097-2765(01)00218-0.PubMedView Article
  39. Thiru A, Nietlispach D, Mott HR, Okuwaki M, Lyon D, Nielsen PR, Hirshberg M, Verreault A, Murzina NV, Laue ED: Structural basis of HP1/PXVXL motif peptide interactions and HP1 localisation to heterochromatin. EMBO J. 2004, 23: 489-499. 10.1038/sj.emboj.7600088.PubMed CentralPubMedView Article
  40. Singh PB: HP1 proteins: what is the essential interaction?. Genetika. 2010, 46: 1424-1429.PubMed
  41. Rountree MR, Selker EU: DNA methylation and the formation of heterochromatin in Neurospora crassa. Heredity. 2010, 105: 38-44. 10.1038/hdy.2010.44.PubMedView Article
  42. Groner AC, Meylan S, Ciuffi A, Zangger N, Ambrosini G, Dénervaud N, Bucher P, Trono D: KRAB-zinc finger proteins and KAP1 can mediate long-range transcriptional repression through heterochromatin spreading. PLoS Genet. 2010, 6: e1000869-10.1371/journal.pgen.1000869.PubMed CentralPubMedView Article
  43. Sadaie M, Kawaguchi R, Ohtani Y, Arisaka F, Tanaka K, Shirahige K, Nakayama J: Balance between distinct HP1 family proteins controls heterochromatin assembly in fission yeast. Mol Cell Biol. 2008, 28: 6973-6988. 10.1128/MCB.00791-08.PubMed CentralPubMedView Article
  44. Kourmouli N, Sun YM, van der Sar S, Singh PB, Brown JP: Epigenetic regulation of mammalian pericentric heterochromatin in vivo by HP1. Biochem Biophys Res Commun. 2005, 337: 901-907. 10.1016/j.bbrc.2005.09.132.PubMedView Article
  45. Nielsen SJ, Schneider R, Bauer UM, Bannister AJ, Morrison A, O'Carroll D, Firestein R, Cleary M, Jenuwein T, Herrera RE, Kouzarides T: Rb targets histone H3 methylation and HP1 to promoters. Nature. 2001, 412: 561-565. 10.1038/35087620.PubMedView Article
  46. Kwon SH, Florens L, Swanson SK, Washburn MP, Abmayr SM, Workman JL: Heterochromatin protein 1 (HP1) connects the FACT histone chaperone complex to the phosphorylated CTD of RNA polymerase II. Genes Dev. 2010, 24: 2133-2145. 10.1101/gad.1959110.PubMed CentralPubMedView Article
  47. Klattenhoff C, Xi H, Li C, Lee S, Xu J, Khurana JS, Zhang F, Schultz N, Koppetsch BS, Nowosielska A, Seitz H, Zamore PD, Weng Z, Theurkauf WE: The Drosophila HP1 homologue Rhino is required for transposon silencing and piRNA production by dual-strand clusters. Cell. 2009, 138: 1137-1149. 10.1016/j.cell.2009.07.014.PubMed CentralPubMedView Article
  48. Freitag M, Hickey PC, Khlafallah TK, Read ND, Selker EU: HP1 is essential for DNA methylation in Neurospora. Mol Cell. 2004, 13: 427-434. 10.1016/S1097-2765(04)00024-3.PubMedView Article
  49. Tamaru H, Selker EU: A histone H3 methyltransferase controls DNA methylation in Neurospora crassa. Nature. 2001, 414: 277-283. 10.1038/35104508.PubMedView Article
  50. Johnson LM, Cao X, Jacobsen SE: Interplay between two epigenetic marks: DNA methylation and histone H3 lysine 9 methylation. Curr Biol. 2002, 12: 1360-1367. 10.1016/S0960-9822(02)00976-4.PubMedView Article
  51. Jackson JP, Lindroth AM, Cao X, Jacobsen SE: Control of CpNpG DNA methylation by the KRYPTONITE histone H3 methyltransferase. Nature. 2002, 416: 556-560. 10.1038/nature731.PubMedView Article
  52. du Chéné I, Basyuk E, Lin YL, Triboulet R, Knezevich A, Chable-Bessia C, Mettling C, Baillat V, Reynes J, Corbeau P, Bertrand E, Marcello A, Emiliani S, Kiernan R, Benkirane M: Suv39H1 and HP1γ are responsible for chromatin-mediated HIV-1 transcriptional silencing and post-integration latency. EMBO J. 2007, 26: 424-435. 10.1038/sj.emboj.7601517.PubMedView Article
  53. Brown JP, Bullwinkel J, Baron-Lühr B, Billur M, Schneider P, Winking H, Singh PB: HP1γ function is required for male germ cell survival and spermatogenesis. Epigenetics Chromatin. 2010, 3: 9-10.1186/1756-8935-3-9.PubMed CentralPubMedView Article
  54. Vakoc CR, Mandat SA, Olenchock BA, Blobel GA: Histone H3 lysine 9 methylation and HP1γ are associated with transcription elongation through mammalian chromatin. Mol Cell. 2005, 19: 381-391. 10.1016/j.molcel.2005.06.011.PubMedView Article
  55. Mateescu B, Bourachot B, Rachez C, Ogryzko V, Muchardt C: Regulation of an inducible promoter by an HP1β-HP1γ switch. EMBO Rep. 2008, 9: 267-272. 10.1038/embor.2008.1.PubMed CentralPubMedView Article
  56. Wolf D, Cammas F, Losson R, Goff SP: Primer binding site-dependent restriction of murine leukemia virus requires HP1 binding by TRIM28. J Virol. 2008, 82: 4675-4679. 10.1128/JVI.02445-07.PubMed CentralPubMedView Article
  57. 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. 10.1128/MCB.00487-06.PubMed CentralPubMedView Article
  58. Schultz DC, Ayyanathan K, Negorev D, Maul GG, Rauscher FJ: 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. 10.1101/gad.973302.PubMed CentralPubMedView Article
  59. 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. 10.1128/MCB.25.7.2525-2538.2005.PubMed CentralPubMedView Article
  60. Nielsen AL, Ortiz JA, You J, Oulad-Abdelghani M, Khechumian R, Gansmuller A, Chambon P, Losson R: 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. 10.1093/emboj/18.22.6385.PubMed CentralPubMedView Article
  61. Ryan RF, Schultz DC, Ayyanathan K, Singh PB, Friedman JR, Fredericks WJ, Rauscher FJ: KAP-1 corepressor protein interacts and colocalizes with heterochromatic and euchromatic HP1 proteins: a potential role for Krüppel-associated box-zinc finger proteins in heterochromatin-mediated gene silencing. Mol Cell Biol. 1999, 19: 4366-4378.PubMed CentralPubMedView Article
  62. Danzer JR, Wallrath LL: Mechanisms of HP1-mediated gene silencing in Drosophila. Development. 2004, 131: 3571-3580. 10.1242/dev.01223.PubMedView Article
  63. Hines KA, Cryderman DE, Flannery KM, Yang H, Vitalini MW, Hazelrigg T, Mizzen CA, Wallrath LL: Domains of heterochromatin protein 1 required for Drosophila melanogaster heterochromatin spreading. Genetics. 2009, 182: 967-977. 10.1534/genetics.109.105338.PubMed CentralPubMedView Article
  64. Partridge JF, Borgstrom B, Allshire RC: Distinct protein interaction domains and protein spreading in a complex centromere. Genes Dev. 2000, 14: 783-791.PubMed CentralPubMed
  65. Aagaard L, Laible G, Selenko P, Schmid M, Dorn R, Schotta G, Kuhfittig S, Wolf A, Lebersorger A, Singh PB, Reuter G, Jenuwein T: 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. 10.1093/emboj/18.7.1923.PubMed CentralPubMedView Article
  66. Li Y, Kirschmann DA, Wallrath LL: Does heterochromatin protein 1 always follow code?. Proc Natl Acad Sci USA. 2002, 99 (Suppl 4): 16462-16469.PubMed CentralPubMedView Article
  67. Locke SM, Martienssen RA: Slicing and spreading of heterochromatic silencing by RNA interference. Cold Spring Harb Symp Quant Biol. 2006, 71: 497-503. 10.1101/sqb.2006.71.062.PubMedView Article
  68. Talbert PB, Henikoff S: Spreading of silent chromatin: inaction at a distance. Nat Rev Genet. 2006, 7: 793-803.PubMedView Article
  69. Honda S, Lewis ZA, Huarte M, Cho LY, David LL, Shi Y, Selker EU: The DMM complex prevents spreading of DNA methylation from transposons to nearby genes in Neurospora crassa. Genes Dev. 2010, 24: 443-454. 10.1101/gad.1893210.PubMed CentralPubMedView Article
  70. Tajul-Arifin K, Teasdale R, Ravasi T, Hume DA, Mattick JS: Identification and analysis of chromodomain-containing proteins encoded in the mouse transcriptome. Genome Res. 2003, 13: 1416-1429. 10.1101/gr.1015703.PubMed CentralPubMedView Article
  71. Kokura K, Sun L, Bedford MT, Fang J: Methyl-H3K9-binding protein MPP8 mediates E-cadherin gene silencing and promotes tumour cell motility and invasion. EMBO J. 2010, 29: 3673-3687. 10.1038/emboj.2010.239.PubMed CentralPubMedView Article
  72. Bua DJ, Kuo AJ, Cheung P, Liu CL, Migliori V, Espejo A, Casadio F, Bassi C, Amati B, Bedford MT, Guccione E, Gozani O: Epigenome microarray platform for proteome-wide dissection of chromatin-signaling networks. PLoS One. 2009, 4: e6789-10.1371/journal.pone.0006789.PubMed CentralPubMedView Article
  73. Liu H, Galka M, Iberg A, Wang Z, Li L, Voss C, Jiang X, Lajoie G, Huang Z, Bedford MT, Li SS: Systematic identification of methyllysine-driven interactions for histone and nonhistone targets. J Proteome Res. 2010, 9: 5827-5836. 10.1021/pr100597b.PubMedView Article
  74. Fischle W, Franz H, Jacobs SA, Allis CD, Khorasanizadeh S: Specificity of the chromodomain Y chromosome family of chromodomains for lysine-methylated ARK(S/T) motifs. J Biol Chem. 2008, 283: 19626-19635. 10.1074/jbc.M802655200.PubMed CentralPubMedView Article
  75. Kim J, Daniel J, Espejo A, Lake A, Krishna M, Xia L, Zhang Y, Bedford MT: Tudor, MBT and chromo domains gauge the degree of lysine methylation. EMBO Rep. 2006, 7: 397-403.PubMed CentralPubMed
  76. Mulligan P, Westbrook TF, Ottinger M, Pavlova N, Chang B, Macia E, Shi YJ, Barretina J, Liu J, Howley PM, Elledge SJ, Shi Y: CDYL bridges REST and histone methyltransferases for gene repression and suppression of cellular transformation. Mol Cell. 2008, 32: 718-726. 10.1016/j.molcel.2008.10.025.PubMedView Article
  77. Bernstein E, Duncan EM, Masui O, Gil J, Heard E, Allis CD: Mouse polycomb proteins bind differentially to methylated histone H3 and RNA and are enriched in facultative heterochromatin. Mol Cell Biol. 2006, 26: 2560-2569. 10.1128/MCB.26.7.2560-2569.2006.PubMed CentralPubMedView Article
  78. Quinn AM, Bedford MT, Espejo A, Spannhoff A, Austin CP, Oppermann U, Simeonov A: A homogeneous method for investigation of methylation-dependent protein-protein interactions in epigenetics. Nucleic Acids Res. 2009, 38: e11.PubMed CentralPubMedView Article
  79. Rottach A, Frauer C, Pichler G, Bonapace IM, Spada F, Leonhardt H: The multi-domain protein Np95 connects DNA methylation and histone modification. Nucleic Acids Res. 2010, 38: 1796-1804. 10.1093/nar/gkp1152.PubMed CentralPubMedView Article
  80. Bernard D, Martinez-Leal JF, Rizzo S, Martinez D, Hudson D, Visakorpi T, Peters G, Carnero A, Beach D, Gil J: CBX7 controls the growth of normal and tumor-derived prostate cells by repressing the Ink4a/Arf locus. Oncogene. 2005, 24: 5543-5551. 10.1038/sj.onc.1208735.PubMedView Article
  81. Aucott R, Bullwinkel J, Yu Y, Shi W, Billur M, Brown JP, Menzel U, Kioussis D, Wang G, Reisert I, Weimer J, Pandita RK, Sharma GG, Pandita TK, Fundele R, Singh PB: HP1-β is required for development of the cerebral neocortex and neuromuscular junctions. J Cell Biol. 2008, 183: 597-606. 10.1083/jcb.200804041.PubMed CentralPubMedView Article
  82. Schotta G, Lachner M, Sarma K, Ebert A, Sengupta R, Reuter G, Reinberg D, Jenuwein T: A silencing pathway to induce H3-K9 and H4-K20 trimethylation at constitutive heterochromatin. Genes Dev. 2004, 18: 1251-1262. 10.1101/gad.300704.PubMed CentralPubMedView Article
  83. Dong KB, Maksakova IA, Mohn F, Leung D, Appanah R, Lee S, Yang HW, Lam LL, Mager DL, Schübeler D, Tachibana M, Shinkai Y, Lorincz MC: DNA methylation in ES cells requires the lysine methyltransferase G9a but not its catalytic activity. EMBO J. 2008, 27: 2691-2701. 10.1038/emboj.2008.193.PubMed CentralPubMedView Article
  84. Lehnertz B, Ueda Y, Derijck AA, Braunschweig U, Perez-Burgos L, Kubicek S, Chen T, Li E, Jenuwein T, Peters AH: Suv39h-mediated histone H3 lysine 9 methylation directs DNA methylation to major satellite repeats at pericentric heterochromatin. Curr Biol. 2003, 13: 1192-1200. 10.1016/S0960-9822(03)00432-9.PubMedView Article
  85. Honda S, Selker EU: Direct interaction between DNA methyltransferase DIM-2 and HP1 is required for DNA methylation in Neurospora crassa. Mol Cell Biol. 2008, 28: 6044-6055. 10.1128/MCB.00823-08.PubMed CentralPubMedView Article
  86. Maksakova IA, Zhang Y, Mager DL: Preferential epigenetic suppression of the autonomous MusD over the nonautonomous ETn mouse retrotransposons. Mol Cell Biol. 2009, 29: 2456-2468. 10.1128/MCB.01383-08.PubMed CentralPubMedView Article
  87. Feng YQ, Seibler J, Alami R, Eisen A, Westerman KA, Leboulch P, Fiering S, Bouhassira EE: Site-specific chromosomal integration in mammalian cells: highly efficient CRE recombinase-mediated cassette exchange. J Mol Biol. 1999, 292: 779-785. 10.1006/jmbi.1999.3113.PubMedView Article
  88. Schübeler D, Lorincz MC, Groudine M: Targeting silence: the use of site-specific recombination to introduce in vitro methylated DNA into the genome. Sci STKE. 2001, 2001 (83): l1.
  89. Lorincz MC, Schübeler D, Groudine M: Methylation-mediated proviral silencing is associated with MeCP2 recruitment and localized histone H3 deacetylation. Mol Cell Biol. 2001, 21: 7913-7922. 10.1128/MCB.21.23.7913-7922.2001.PubMed CentralPubMedView Article
  90. Teich NM, Weiss RA, Martin GR, Lowy DR: Virus infection of murine teratocarcinoma stem cell lines. Cell. 1977, 12: 973-982. 10.1016/0092-8674(77)90162-3.PubMedView Article
  91. Pannell D, Osborne CS, Yao S, Sukonnik T, Pasceri P, Karaiskakis A, Okano M, Li E, Lipshitz HD, Ellis J: Retrovirus vector silencing is de novo methylase independent and marked by a repressive histone code. EMBO J. 2000, 19: 5884-5894. 10.1093/emboj/19.21.5884.PubMed CentralPubMedView Article
  92. Niwa O, Yokota Y, Ishida H, Sugahara T: Independent mechanisms involved in suppression of the Moloney leukemia virus genome during differentiation of murine teratocarcinoma cells. Cell. 1983, 32: 1105-1113. 10.1016/0092-8674(83)90294-5.PubMedView Article
  93. Cheng L, Du C, Murray D, Tong X, Zhang YA, Chen BP, Hawley RG: A GFP reporter system to assess gene transfer and expression in human hematopoietic progenitor cells. Gene Ther. 1997, 4: 1013-1022. 10.1038/ Article
  94. Musselman CA, Mansfield RE, Garske AL, Davrazou F, Kwan AH, Oliver SS, O'Leary H, Denu JM, Mackay JP, Kutateladze TG: Binding of the CHD4 PHD2 finger to histone H3 is modulated by covalent modifications. Biochem J. 2009, 423: 179-187. 10.1042/BJ20090870.PubMed CentralPubMedView Article
  95. Mansfield RE, Musselman CA, Kwan AH, Oliver SS, Garske AL, Davrazou F, Denu JM, Kutateladze TG, Mackay JP: The plant homeodomain (PHD) fingers of CHD4 are histone H3-binding modules with preference for unmodified H3K4 and methylated H3K9. J Biol Chem. 2011, 286: 11779-11791. 10.1074/jbc.M110.208207.PubMed CentralPubMedView Article
  96. Iwase S, Lan F, Bayliss P, de la Torre-Ubieta L, Huarte M, Qi HH, Whetstine JR, Bonni A, Roberts TM, Shi Y: The X-linked mental retardation gene SMCX/JARID1C defines a family of histone H3 lysine 4 demethylases. Cell. 2007, 128: 1077-1088. 10.1016/j.cell.2007.02.017.PubMedView Article
  97. Li F, Huarte M, Zaratiegui M, Vaughn MW, Shi Y, Martienssen R, Cande WZ: Lid2 is required for coordinating H3K4 and H3K9 methylation of heterochromatin and euchromatin. Cell. 2008, 135: 272-283. 10.1016/j.cell.2008.08.036.PubMed CentralPubMedView Article
  98. Karagianni P, Amazit L, Qin J, Wong J: ICBP90, a novel methyl K9 H3 binding protein linking protein ubiquitination with heterochromatin formation. Mol Cell Biol. 2008, 28: 705-717. 10.1128/MCB.01598-07.PubMed CentralPubMedView Article
  99. Bartke T, Vermeulen M, Xhemalce B, Robson SC, Mann M, Kouzarides T: Nucleosome-interacting proteins regulated by DNA and histone methylation. Cell. 2010, 143: 470-484. 10.1016/j.cell.2010.10.012.PubMed CentralPubMedView Article
  100. Nady N, Lemak A, Walker JR, Avvakumov GV, Kareta MS, Achour M, Xue S, Duan S, Allali-Hassani A, Zuo X, Wang YX, Bronner C, Chédin F, Arrowsmith CH, Dhe-Paganon S: Recognition of multivalent histone states associated with heterochromatin by UHRF1 protein. J Biol Chem. 2011, 286: 24300-24311. 10.1074/jbc.M111.234104.PubMed CentralPubMedView Article
  101. Papait R, Pistore C, Negri D, Pecoraro D, Cantarini L, Bonapace IM: Np95 is implicated in pericentromeric heterochromatin replication and in major satellite silencing. Mol Biol Cell. 2007, 18: 1098-1106. 10.1091/mbc.E06-09-0874.PubMed CentralPubMedView Article
  102. Seum C, Delattre M, Spierer A, Spierer P: Ectopic HP1 promotes chromosome loops and variegated silencing in Drosophila. EMBO J. 2001, 20: 812-818. 10.1093/emboj/20.4.812.PubMed CentralPubMedView Article
  103. Seum C, Spierer A, Delattre M, Pauli D, Spierer P: A GAL4-HP1 fusion protein targeted near heterochromatin promotes gene silencing. Chromosoma. 2000, 109: 453-459. 10.1007/s004120000101.PubMedView Article
  104. Agalioti T, Lomvardas S, Parekh B, Yie J, Maniatis T, Thanos D: Ordered recruitment of chromatin modifying and general transcription factors to the IFN-β promoter. Cell. 2000, 103: 667-678. 10.1016/S0092-8674(00)00169-0.PubMedView Article
  105. Kasten M, Szerlong H, Erdjument-Bromage H, Tempst P, Werner M, Cairns BR: Tandem bromodomains in the chromatin remodeler RSC recognize acetylated histone H3 Lys14. EMBO J. 2004, 23: 1348-1359. 10.1038/sj.emboj.7600143.PubMed CentralPubMedView Article
  106. Vicent GP, Zaurin R, Nacht AS, Li A, Font-Mateu J, Le Dily F, Vermeulen M, Mann M, Beato M: Two chromatin remodeling activities cooperate during activation of hormone responsive promoters. PLoS Genet. 2009, 5: e1000567-10.1371/journal.pgen.1000567.PubMed CentralPubMedView Article
  107. Hassan AH, Neely KE, Workman JL: Histone acetyltransferase complexes stabilize SWI/SNF binding to promoter nucleosomes. Cell. 2001, 104: 817-827. 10.1016/S0092-8674(01)00279-3.PubMedView Article
  108. Wang Z, Zang C, Cui K, Schones DE, Barski A, Peng W, Zhao K: Genome-wide mapping of HATs and HDACs reveals distinct functions in active and inactive genes. Cell. 2009, 138: 1019-1031. 10.1016/j.cell.2009.06.049.PubMed CentralPubMedView Article
  109. Jin Q, Yu LR, Wang L, Zhang Z, Kasper LH, Lee JE, Wang C, Brindle PK, Dent SY, Ge K: Distinct roles of GCN5/PCAF-mediated H3K9ac and CBP/p300-mediated H3K18/27ac in nuclear receptor transactivation. EMBO J. 2011, 30: 249-262. 10.1038/emboj.2010.318.PubMed CentralPubMedView Article
  110. Nagy Z, Riss A, Fujiyama S, Krebs A, Orpinell M, Jansen P, Cohen A, Stunnenberg HG, Kato S, Tora L: The metazoan ATAC and SAGA coactivator HAT complexes regulate different sets of inducible target genes. Cell Mol Life Sci. 2010, 67: 611-628. 10.1007/s00018-009-0199-8.PubMedView Article
  111. Nagy Z, Tora L: Distinct GCN5/PCAF-containing complexes function as co-activators and are involved in transcription factor and global histone acetylation. Oncogene. 2007, 26: 5341-5357. 10.1038/sj.onc.1210604.PubMedView Article
  112. Huang J, Fan T, Yan Q, Zhu H, Fox S, Issaq HJ, Best L, Gangi L, Munroe D, Muegge K: Lsh, an epigenetic guardian of repetitive elements. Nucleic Acids Res. 2004, 32: 5019-5028. 10.1093/nar/gkh821.PubMed CentralPubMedView Article
  113. Vicent GP, Ballaré C, Nacht AS, Clausell J, Subtil-Rodríguez A, Quiles I, Jordan A, Beato M: Induction of progesterone target genes requires activation of Erk and Msk kinases and phosphorylation of histone H3. Mol Cell. 2006, 24: 367-381. 10.1016/j.molcel.2006.10.011.PubMedView Article
  114. Hebbar PB, Archer TK: Nuclear factor 1 is required for both hormone-dependent chromatin remodeling and transcriptional activation of the mouse mammary tumor virus promoter. Mol Cell Biol. 2003, 23: 887-898. 10.1128/MCB.23.3.887-898.2003.PubMed CentralPubMedView Article
  115. Vicent GP, Zaurin R, Nacht AS, Font-Mateu J, Le Dily F, Beato M: Nuclear factor 1 synergizes with progesterone receptor on the mouse mammary tumor virus promoter wrapped around a histone H3/H4 tetramer by facilitating access to the central hormone-responsive elements. J Biol Chem. 2010, 285: 2622-2631. 10.1074/jbc.M109.060848.PubMed CentralPubMedView Article
  116. Baust C, Gagnier L, Baillie GJ, Harris MJ, Juriloff DM, Mager DL: Structure and expression of mobile ETnII retroelements and their coding-competent MusD relatives in the mouse. J Virol. 2003, 77: 11448-11458. 10.1128/JVI.77.21.11448-11458.2003.PubMed CentralPubMedView Article
  117. Rea S, Eisenhaber F, O'Carroll D, Strahl BD, Sun ZW, Schmid M, Opravil S, Mechtler K, Ponting CP, Allis CD, Jenuwein T: Regulation of chromatin structure by site-specific histone H3 methyltransferases. Nature. 2000, 406: 593-599. 10.1038/35020506.PubMedView Article
  118. Gregory GD, Vakoc CR, Rozovskaia T, Zheng X, Patel S, Nakamura T, Canaani E, Blobel GA: Mammalian ASH1L is a histone methyltransferase that occupies the transcribed region of active genes. Mol Cellular Biol. 2007, 27: 8466-8479. 10.1128/MCB.00993-07.View Article
  119. Wang H, Cao R, Xia L, Erdjument-Bromage H, Borchers C, Tempst P, Zhang Y: Purification and functional characterization of a histone H3-lysine 4-specific methyltransferase. Mol Cell. 2001, 8: 1207-1217. 10.1016/S1097-2765(01)00405-1.PubMedView Article
  120. Yan Q, Huang J, Fan T, Zhu H, Muegge K: Lsh, a modulator of CpG methylation, is crucial for normal histone methylation. EMBO J. 2003, 22: 5154-5162. 10.1093/emboj/cdg493.PubMed CentralPubMedView Article
  121. Wysocka J, Swigut T, Milne TA, Dou Y, Zhang X, Burlingame AL, Roeder RG, Brivanlou AH, Allis CD: WDR5 associates with histone H3 methylated at K4 and is essential for H3 K4 methylation and vertebrate development. Cell. 2005, 121: 859-872. 10.1016/j.cell.2005.03.036.PubMedView Article
  122. Balakrishnan L, Milavetz B: Decoding the histone H4 lysine 20 methylation mark. Crit Rev Biochem Mol Biol. 2010, 45: 440-452. 10.3109/10409238.2010.504700.PubMedView Article
  123. Gonzalo S, García-Cao M, Fraga MF, Schotta G, Peters AH, Cotter SE, Eguía R, Dean DC, Esteller M, Jenuwein T, Blasco MA: Role of the RB1 family in stabilizing histone methylation at constitutive heterochromatin. Nat Cell Biol. 2005, 7: 420-428. 10.1038/ncb1235.PubMedView Article
  124. Schotta G, Sengupta R, Kubicek S, Malin S, Kauer M, Callén E, Celeste A, Pagani M, Opravil S, De La Rosa-Velazquez IA, Espejo A, Bedford MT, Nussenzweig A, Busslinger M, Jenuwein T: A chromatin-wide transition to H4K20 monomethylation impairs genome integrity and programmed DNA rearrangements in the mouse. Genes Dev. 2008, 22: 2048-2061. 10.1101/gad.476008.PubMed CentralPubMedView Article
  125. Ayoub N, Jeyasekharan AD, Venkitaraman AR: Mobilization and recruitment of HP1β: a bimodal response to DNA breakage. Cell Cycle. 2009, 8: 2945-2950.PubMed
  126. Luijsterburg MS, Dinant C, Lans H, Stap J, Wiernasz E, Lagerwerf S, Warmerdam DO, Lindh M, Brink MC, Dobrucki JW, Aten JA, Fousteri MI, Jansen G, Dantuma NP, Vermeulen W, Mullenders LH, Houtsmuller AB, Verschure PJ, van Driel R: Heterochromatin protein 1 is recruited to various types of DNA damage. J Cell Biol. 2009, 185: 577-586. 10.1083/jcb.200810035.PubMed CentralPubMedView Article
  127. Feng YQ, Lorincz MC, Fiering S, Greally JM, Bouhassira EE: Position effects are influenced by the orientation of a transgene with respect to flanking chromatin. Mol Cell Biol. 2001, 21: 298-309. 10.1128/MCB.21.1.298-309.2001.PubMed CentralPubMedView Article
  128. Juriloff DM, Harris MJ, Dewell SL, Brown CJ, Mager DL, Gagnier L, Mah DG: Investigations of the genomic region that contains the clf1 mutation, a causal gene in multifactorial cleft lip and palate in mice. Birth Defects Res A Clin Mol Teratol. 2005, 73: 103-113. 10.1002/bdra.20106.PubMedView Article
  129. Plamondon JA, Harris MJ, Mager DL, Gagnier L, Juriloff DM: The clf2 gene has an epigenetic role in the multifactorial etiology of cleft lip and palate in the A/WySn mouse strain. Birth Defects Res A Clin Mol Teratol.
  130. Bouhassira EE, Westerman K, Leboulch P: Transcriptional behavior of LCR enhancer elements integrated at the same chromosomal locus by recombinase-mediated cassette exchange. Blood. 1997, 90: 3332-3344.PubMed
  131. Kumaki Y, Oda M, Okano M: QUMA: quantification tool for methylation analysis. Nucleic Acids Res. 2008, W170-W175. 36 Web Server
  132. Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110: 462-467. 10.1159/000084979.PubMedView Article
  133. Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralPubMedView Article


© Maksakova et al; licensee BioMed Central Ltd. 2011

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.