Title: Autocatalytic differentiation of epigenetic modifications within the Arabidopsis genome
Abstract: Article10 September 2010free access Autocatalytic differentiation of epigenetic modifications within the Arabidopsis genome Soichi Inagaki Soichi Inagaki Department of Integrated Genetics, National Institute of Genetics, Mishima, Shizuoka, JapanPresent address: Graduate School of Biological Sciences, Nara Institute of Science and Technology, 8916-5 Takayama, Ikoma, Nara 630-0192, Japan Search for more papers by this author Asuka Miura-Kamio Asuka Miura-Kamio Department of Integrated Genetics, National Institute of Genetics, Mishima, Shizuoka, Japan Search for more papers by this author Yasukazu Nakamura Yasukazu Nakamura Center for Information Biology, National Institute of Genetics, Mishima, Shizuoka, Japan Search for more papers by this author Falong Lu Falong Lu State Key Laboratory of Plant Genomics and National Center for Plant Gene Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, PR, China Search for more papers by this author Xia Cui Xia Cui State Key Laboratory of Plant Genomics and National Center for Plant Gene Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, PR, China Search for more papers by this author Xiaofeng Cao Corresponding Author Xiaofeng Cao State Key Laboratory of Plant Genomics and National Center for Plant Gene Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, PR, China Search for more papers by this author Hiroshi Kimura Hiroshi Kimura Graduate School of Frontier Biosciences, Osaka University, Suita, Osaka, Japan Search for more papers by this author Hidetoshi Saze Hidetoshi Saze Department of Integrated Genetics, National Institute of Genetics, Mishima, Shizuoka, Japan Search for more papers by this author Tetsuji Kakutani Corresponding Author Tetsuji Kakutani Department of Integrated Genetics, National Institute of Genetics, Mishima, Shizuoka, Japan Search for more papers by this author Soichi Inagaki Soichi Inagaki Department of Integrated Genetics, National Institute of Genetics, Mishima, Shizuoka, JapanPresent address: Graduate School of Biological Sciences, Nara Institute of Science and Technology, 8916-5 Takayama, Ikoma, Nara 630-0192, Japan Search for more papers by this author Asuka Miura-Kamio Asuka Miura-Kamio Department of Integrated Genetics, National Institute of Genetics, Mishima, Shizuoka, Japan Search for more papers by this author Yasukazu Nakamura Yasukazu Nakamura Center for Information Biology, National Institute of Genetics, Mishima, Shizuoka, Japan Search for more papers by this author Falong Lu Falong Lu State Key Laboratory of Plant Genomics and National Center for Plant Gene Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, PR, China Search for more papers by this author Xia Cui Xia Cui State Key Laboratory of Plant Genomics and National Center for Plant Gene Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, PR, China Search for more papers by this author Xiaofeng Cao Corresponding Author Xiaofeng Cao State Key Laboratory of Plant Genomics and National Center for Plant Gene Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, PR, China Search for more papers by this author Hiroshi Kimura Hiroshi Kimura Graduate School of Frontier Biosciences, Osaka University, Suita, Osaka, Japan Search for more papers by this author Hidetoshi Saze Hidetoshi Saze Department of Integrated Genetics, National Institute of Genetics, Mishima, Shizuoka, Japan Search for more papers by this author Tetsuji Kakutani Corresponding Author Tetsuji Kakutani Department of Integrated Genetics, National Institute of Genetics, Mishima, Shizuoka, Japan Search for more papers by this author Author Information Soichi Inagaki1, Asuka Miura-Kamio1, Yasukazu Nakamura2, Falong Lu3, Xia Cui3, Xiaofeng Cao 3, Hiroshi Kimura4, Hidetoshi Saze1 and Tetsuji Kakutani 1 1Department of Integrated Genetics, National Institute of Genetics, Mishima, Shizuoka, Japan 2Center for Information Biology, National Institute of Genetics, Mishima, Shizuoka, Japan 3State Key Laboratory of Plant Genomics and National Center for Plant Gene Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, PR, China 4Graduate School of Frontier Biosciences, Osaka University, Suita, Osaka, Japan *Corresponding authors: Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China, Tel.: +86 10 64869203; Fax.: +86 10 64873428; E-mail: [email protected] of Integrated Genetics, National Institute of Genetics, Mishima, Shizuoka 411-8540, Japan. Tel.: +81 55 981 6801; Fax: +81 55 981 6804; E-mail: [email protected] The EMBO Journal (2010)29:3496-3506https://doi.org/10.1038/emboj.2010.227 PDFDownload PDF of article text and main figures. Peer ReviewDownload a summary of the editorial decision process including editorial decision letters, reviewer comments and author responses to feedback. ToolsAdd to favoritesDownload CitationsTrack CitationsPermissions ShareFacebookTwitterLinked InMendeleyWechatReddit Figures & Info In diverse eukaryotes, constitutively silent sequences, such as transposons and repeats, are marked by methylation at histone H3 lysine 9 (H3K9me). Although selective H3K9me is critical for maintaining genome integrity, mechanisms to exclude H3K9me from active genes remain largely unexplored. Here, we show in Arabidopsis that the exclusion depends on a histone demethylase gene, IBM1 (increase in BONSAI methylation). Loss-of-function ibm1 mutation results in ectopic H3K9me and non-CG methylation in thousands of genes. The ibm1-induced genic H3K9me depends on both histone methylase KYP/SUVH4 and DNA methylase CMT3, suggesting interdependence of two epigenetic marks—H3K9me and non-CG methylation. Notably, IBM1 enhances loss of H3K9me in transcriptionally de-repressed sequences. Furthermore, disruption of transcription in genes induces ectopic non-CG methylation, which mimics the loss of IBM1 function. We propose that active chromatin is stabilized by an autocatalytic loop of transcription and H3K9 demethylation. This process counteracts a similarly autocatalytic accumulation of silent epigenetic marks, H3K9me and non-CG methylation. Introduction Chromatin modifications, such as methylation of histone H3 at lysine 9 (H3K9me) and DNA methylation at cytosine residues (Bender, 2004; Chan et al, 2005; Grewal and Elgin, 2007), are important for constitutively silencing transposons (Yoder et al, 1997; Walsh et al, 1998; Miura et al, 2001; Lippman et al, 2004; Ding et al, 2007; Mirouze et al, 2009; Tsukahara et al, 2009). On the other hand, these modifications are normally excluded from active genes to ensure their appropriate expression. Mechanisms for the exclusion, however, remain largely unexplored. In Arabidopsis, methylation in genes is found only at CG sites, although transposons are methylated at all contexts of cytosines, CG, CHG, and CHH (H can be A, C, or T) (Zhang et al, 2006; Zilberman et al, 2007; Cokus et al, 2008; Lister et al, 2008). The CHG methylation in transposons is associated with dimethylation of H3K9 (H3K9me2) (Bernatavichute et al, 2008). H3K9me2 is catalysed by multiple histone methyltransferases, KRYPTONITE/SUVH4 (hereafter KYP), SUVH5, and SUVH6 (Jackson et al, 2002; Malagnac et al, 2002; Ebbs et al, 2005; Ebbs and Bender, 2006). These histone methyltransferases act upstream of DNA methyltransferase CHROMOMETHYLASE3 (CMT3), which catalyses DNA methylation at CHG sites (Bartee et al, 2001; Lindroth et al, 2001; Chan et al, 2005). In addition to the positive regulators of DNA methylation, negative regulators have been recently identified (Agius et al, 2006; Gehring et al, 2006; Saze et al, 2008). One of the negative regulators of DNA methylation is IBM1 (INCREASE IN BONSAI METHYLATION 1); loss-of-function mutation in the IBM1 gene induces ectopic DNA methylation at non-CG sites (mainly at CHG sites) in thousands of genes, whereas transposons are unaffected (Saze et al, 2008; Miura et al, 2009). Interestingly, the ibm1-induced genic DNA methylation occurs across each individual transcriptional unit (Miura et al, 2009). It is not known how a transcriptional unit can be recognized by the machinery that excludes the transposon-specific modification. The IBM1 protein has a jumonji-C (jmjC) domain, and it is a member of the JHDM2/KDM3 family, which contains H3K9 demethylases (Klose et al, 2006; Lu et al, 2008; Sun and Zhou, 2008; Liu et al, 2010). To understand the mechanisms that establish differential epigenetic modifications between active genes and silent transposons, we examined genome-wide pattern of H3K9me2 using mutants for IBM1, KYP, and CMT3. The ibm1 mutation induced accumulation of H3K9me2 in thousands of genes, which parallels the increases in DNA methylation. Although IBM1 generally affected H3K9me2 of genes, transposons could also be affected, once they were activated by kyp or cmt3 mutation, suggesting that transcription is an important prerequisite for the IBM1 function. Furthermore, blockage of transcription by T-DNA insertion in the target genes of IBM1 caused ectopic CHG methylation in the body of these genes, which mimics the effect of the loss of IBM1 function. These results lead us to propose that transcription is critical for IBM1 to remove H3K9me2, which mediates stabilization of active chromatin by the autocatalytic loop. Results IBM1 is necessary for excluding H3K9me2 from genes To understand how genomic H3K9me2 patterns are generated, we examined the function of the putative H3K9 demethylase, IBM1, as well as the H3K9 methylase KYP and the CHG methylase CMT3. First, we examined demethylase activity of IBM1 in vivo by overexpressing FLAG-tagged IBM1 in Nicotiana benthamiana leaf. Double immunostaining of nuclei with antibodies against FLAG and methylated histone at different residues shows the histone methylation status at each residue in FLAG-IBM1-overexpressing cells. In the FLAG-stained cells, H3K9me2 was hardly detectable and the H3K9me1 signal was also much reduced compared with that of the cells without FLAG labelling (Figure 1A; Supplementary Figure S1A). The effects on H3K9me2 and H3K9me1 levels were abolished when two amino-acid residues are substituted (H679A and D681A; referred to as FLAG-IBM1m) in the conserved Fe(II) binding site of the IBM1 protein (Figure 1B; Supplementary Figure S1B). Figure 1.IBM1 downregulates H3K9me2 at genes. (A, B) Overexpression of FLAG-IBM1 in tobacco leaf cells induced reduction in H3K9me2 and H3K9me1. Nuclei transfected with FLAG-IBM1 (A) or a mutated construct, FLAG-IBM1m (B), were mixed with control nuclei without the transfection (see Materials and methods for details). The transfected nuclei and control nuclei could be distinguished by immunostaining with anti-FLAG (shown in red). Histone methylation was visualized by immunostaining with rabbit polyclonal methylation-specific antibodies (green). Nuclei were visualized by DAPI staining (blue). Arrows indicate nuclei transfected with FLAG-IBM1 or FLAG-IBM1m. Quantification of these signals, together with those from other samples, is shown in Supplementary Figure S1A and B. (C, D) Representative examples of H3K9me2 analysed by tilling array for a gene-rich region (C, position 20 829 000–20 934 000 on chromosome 1) and a transposon-rich region (D, positions 16 494 000–16 592 000 on chromosome 1). Profile of DNA methylation (Miura et al, 2009) is also shown. Blue and red boxes indicate genes and transposons, respectively. They are oriented 5′–3′ on the top and 3′–5′ on the bottom. Each vertical bar represents the log2 ratio of the immunoprecipitated DNA to input control for the DNA methylation (DNAme) and H3K9me2, respectively. Each brown bar represents difference of signals for each comparison. (E) H3K9me2 level in each gene (blue dot) and transposon (red dot), compared between wild type and ibm1. (F) Comparison between the effect of ibm1 mutation on DNA methylation (DNAme) and H3K9me2 in each gene (blue dot) and transposon (red dot). The changes in DNA methylation and H3K9me2 levels were not proportional, most likely because each gene differs in the CG methylation level, which is not affected by the ibm1 mutation (Miura et al, 2009). Download figure Download PowerPoint H3K9me3 and the three types of methylation at each of other sites, H3K4, H3K27, and H3K36, were indistinguishable between cells with or without the FLAG-IBM1 signal (Supplementary Figure S1C–F). Among the 12 modifications examined, H3K9me2 and H3K9me1 were the only modifications affected by the FLAG-IBM1. These results support the prediction that IBM1 is a histone demethylase specific for H3K9me2 and H3K9me1. We then examined H3K9me2 genome-wide in wild-type Columbia-0 (Col-0) and ibm1 mutant samples by chromatin immunoprecipitation (ChIP) coupled with DNA microarray analysis (ChIP-chip). In order to minimize cross-reacting signals from other modifications, we used a monoclonal antibody specific for H3K9me2 (Kimura et al, 2008). In wild type, H3K9me2 level was high in transposons and repeats, but it was very low in annotated expressed genes (Figure 1C and D; Supplementary Figure S2). The ibm1 mutation caused ectopic H3K9me2 accumulation in transcribed regions of a large number of genes (Figure 1C–E). Essentially the same set of genes was affected in two biological replicates (Supplementary Figures S3). On the other hand, transposons were unaffected by the ibm1 mutation (Figure 1D and E). The ibm1-induced genic H3K9me2 is accompanied by DNA methylation changes (Figure 1F). As is the case in the ibm1-induced DNA methylation in genes, the ibm1-induced H3K9me2 showed specific characteristics in relation to the transcription units (Miura et al, 2009); long genes and constitutively transcribed genes were affected severely (Supplementary Figure S4A–E), and central region of a transcription unit showed more signal than the 5′ and 3′ regions (Supplementary Figure S5A–C). The kyp and cmt3 mutations partially abolish H3K9me2 in transposons In addition to the negative regulator of H3K9me2, we examined the contribution of positive regulators, KYP and CMT3. Compared with wild-type, kyp mutant samples showed reduced H3K9me2 levels in many transposons both in pericentromeric regions and arm regions, although significant levels of H3K9me2 still remained (Figure 2A–C; Supplementary Figure S5D). Other methyltransferase(s) such as SUVH5 and SUVH6 may contribute to the remaining H3K9me2 in transposons and repeats, as has been suggested previously (Ebbs and Bender, 2006). Consistent with the proposed functional redundancy, the H3K9me2 levels were reduced in the triple mutant for KYP, SUVH5, and SUVH6 (suvh456) compared with the single mutant kyp (Figure 2D). Figure 2.Effects of kyp and cmt3 mutations on H3K9me2 level. (A, B) Global view of H3K9me2 patterns along the chromosome 3, compared between wild type and kyp. Average signals in each 100 kb regions (A) or the differences of the signals between wild type and kyp (B) were plotted. (C) H3K9me2 levels in genes (blue dots) and transposons (red dots), compared between wild type and kyp. (D) H3K9me2 levels in wild type (Col), wild-type sibling of suvh456 (SUVH456), kyp, and suvh456 triple mutant examined with ChIP-PCR. Amplifications from input control and ChIP DNA are shown. (E, F) Global view of H3K9me2 patterns in wild type and cmt3 (E) or the difference in H3K9me2 between wild type and cmt3 (F). (G, H) H3K9me2 levels in genes (blue dots) and transposons (red dots), compared between wild type and cmt3 (G) or cmt3 and kyp (H). Download figure Download PowerPoint The cmt3 mutant also showed some reduction of H3K9me2 in heterochromatic regions, although the effect was generally less than that seen in kyp samples (Figure 2E–H; Supplementary Figure S5D). The dependence of H3K9me2 on CMT3 function suggests a function for DNA methylation in maintaining H3K9me2 in transposons. We then examined the relationships between the levels of DNA methylation at different contexts and the H3K9me2 levels in wild-type, kyp, and cmt3 samples. Consistent with previous report (Bernatavichute et al, 2008), the frequency of methylated CHG per length ([mCHG]) was positively correlated with H3K9me2 level in transposons (Supplementary Figure S6A). The kyp and cmt3 mutant samples showed a greater decrease in H3K9me2 at transposons with higher [mCHG] (Supplementary Figure S6C, E, G, and H), suggesting that CHG methylation is indeed required for a high level of H3K9me2 at transposons. Interestingly, the decreases of H3K9me2 in kyp and cmt3 mutants were more pronounced at transposons with lower [mCG] than those with higher [mCG] (Supplementary Figure S6B, D, F, I, and J). The results suggest that an additional pathway, which depends on CG methylation, may also be involved in the H3K9me2 of transposons (Soppe et al, 2002; Tariq et al, 2003). The kyp and cmt3 mutations fully suppress the ibm1-induced H3K9me2 Although the kyp mutation induced only a partial reduction of H3K9me2 in transposons (Figure 2A–C), this mutation suppressed the developmental phenotypes induced by the ibm1 mutation (Saze et al, 2008), which are probably due to ectopic H3K9me2 and/or CHG methylation in genic regions. We therefore examined the effects of the ibm1 mutation in the kyp mutant background. Notably, in the kyp mutant background, ibm1-induced genic H3K9me2 was lost almost completely (Figure 3A and B; Supplementary Figure S7). This result is in striking contrast to the effect of the kyp mutation on H3K9me2 in transposons, which is only partial (Figure 2A–C). Taken together, the results suggest that different H3K9 methyltransferases function in genes and transposons: only KYP in genes and KYP and others in transposons. Figure 3.The ibm1-induced genic H3K9me2 depends on KYP and CMT3. (A) and (D) are from the data shown in Figure 1E. Effect of the ibm1 mutation on H3K9me2 levels was also examined in the background of kyp (B, E) or cmt3 (C, F) mutant. Genes (A–C) and transposons (D–F) are separately shown. Download figure Download PowerPoint Surprisingly, the cmt3 mutation also suppressed the ibm1-induced genic H3K9me2 almost completely (Figure 3C; Supplementary Figure S7). This result suggests that CMT3-mediated CHG methylation is required for the increase of H3K9me2 at genes in the ibm1 mutant. On the other hand, the impact of the cmt3 mutation on H3K9me2 in transposons was only subtle (Figure 2E–G). To test whether CHG methylation, rather than an unidentified function of CMT3, mediates its effect on H3K9me2 induction in ibm1, we calculated the proportion of CHG sites (total number of CHG sites in transcribed region per length of annotation unit; hereafter [CHG]) for all annotated transcription units and examined the relationship of the [CHG] with the effect of the ibm1 mutation on DNA methylation and H3K9me2. [CHG] strongly correlated with the degree of response to the ibm1 mutation not only in DNA methylation, but also in H3K9me2 (Figure 4A and B). Together with the result of ibm1 cmt3 double mutant, these results suggest that CMT3-mediated CHG methylation is required for the induction of genic H3K9me2 in ibm1. Figure 4.[CHG] and [CG] correlate with the effect of ibm1 on H3K9me2. (A, B) Relationships between [CHG] and the ibm1-induced changes in H3K9me2 (A) or DNA methylation (B) level in each gene. (C, D) Relationships between [CG] and the ibm1-induced changes in H3K9me2 (C) or DNA methylation (D) level in each gene. The effects of [CHG] and [CG] are independent of the effect of gene length (Supplementary Figure S8A–D). (E) Positive correlation between [CHG] and [CG]. (F) The proportion of genes with increased H3K9me2 in ibm1 (Supplementary Figure S3A), within four gene groups with high and low [CHG] and [CG]. Actual gene numbers (ibm1 affected/all) in each group are also indicated. Download figure Download PowerPoint Unexpectedly, we also found that [CG] of genes negatively correlated with the induction of H3K9me2 and DNA methylation in the ibm1 mutant (Figure 4C and D). Despite the opposite effects, [CHG] and [CG] show a weak positive correlation (Figure 4E), suggesting that [CHG] and [CG] independently affect H3K9me2 and DNA methylation. Genes with higher [CHG] and lower [CG] were most strongly affected by the ibm1 mutation (Figure 4F). IBM1 affects transposons de-repressed in kyp or cmt3 mutant backgrounds In kyp or cmt3 mutant backgrounds, the ibm1 mutation did not affect H3K9me2 in genes (Figure 3A–C). However, a significant number of transposons, which were not affected by the ibm1 single mutation (Figure 3D), responded to the ibm1 mutation in those backgrounds (Figure 3E and F). These transposons showed a decrease of H3K9me2 in kyp and cmt3 single mutants that was recovered in the ibm1 kyp and ibm1 cmt3 double mutants (Figure 5A–C). Not only H3K9me2, but also CHG methylation was recovered in the ibm1 kyp double mutant for these transposons (Figure 5D). These results suggest that a subset of transposons can be targets of IBM1 in a kyp or cmt3 background. Figure 5.De-repressed transposons can be targets of IBM1. (A) Scatter plot showing the relationship between changes in H3K9me2 from wild type to kyp (x axis) and from kyp to ibm1 kyp (y axis). (B) Scatter plot showing the relationship between changes in H3K9me2 from wild type to cmt3 (x axis) and from cmt3 to ibm1 cmt3 (y axis). Genes are shown in blue and transposons in red. (C) ChIP-chip result of four transposons showing decrease in H3K9me2 in kyp and cmt3, and restoration in ibm1 kyp and ibm1 cmt3 double mutants. (D) DNA methylation at CHG sites analysed by digestion with methylation-sensitive restriction enzymes MspI and BglII and subsequent PCR (details in Materials and methods). Homozygous plants segregating in self-pollinated progeny of a double heterozygote (IBM1/ibm1 KYP/kyp) were used for analyses. Compared with the kyp IBM1 plants with the hypomethylated transposons, the kyp ibm1 plants showed increased signal reflecting robust methylation. The sibling wild-type KYP plants showed high-level methylation irrespective of the genotypes of the IBM1 gene. (E) RT–PCR analysis for the expression of the four transposons in the wild type, ibm1, kyp, and kyp ibm1 double mutant. Total RNA isolated from individual plants was used for each RT–PCR analysis. GAPC gene was used to normalize RNA template amounts. Negative control (no RT) lacked reverse transcriptase. We also examined expression pattern of transposons in kyp and cmt3 mutants genome wide (Supplementary Figure S9). Other transposons de-repressed by kyp or cmt3 mutation also showed increased H3K9me2 in the ibm1 mutant. Download figure Download PowerPoint One possible mechanism for this observation is that even transposons can be targets of IBM1 when they are transcribed. Indeed, these transposons were transcribed in the kyp or cmt3 single mutant (Figure 5E; Supplementary Figure S9). Even though the transcriptional activation tended to be weak in the absence of IBM1 function (ibm1 kyp in Figure 5E), robust activation could be achieved when IBM1 functions (IBM1 kyp in Figure 5E). IBM1 apparently stabilizes the transcriptionally active state, most likely by active H3K9 demethylation. T-DNA insertion induces non-CG methylation within IBM1 target genes These observations suggest that transcription of the target gene is necessary for IBM1 function. In order to test this prediction directly, we examined the effect of blocking transcription by T-DNA insertion. DNA methylation was examined in T-DNA insertion alleles for three of the IBM1 target genes (At4g11130, At1g63020, and At4g32700). In all of them, non-CG (mostly CHG) methylation was induced, a situation which is analogous to the effects of the ibm1 mutation (Figure 6). This observation also supports the prediction that IBM1 function depends on transcription of the target gene. Figure 6.Bisulphite sequencing analysis of cytosine methylation at T-DNA-inserted genes. (A–C) Schematic representations of At4g11130 (A), At1g63020 (B), and At4g32700 (C) loci and cytosine methylation level analysed by bisulphite sequencing. The percentage of methylated cytosines is indicated by vertical bar for each site (red, CG; blue, CHG; green, CHH). T-DNA insertion mutants showed decreases of expressions (Supplementary Figure S10E), and increases of non-CG methylation, as in ibm1 mutant, at downstream regions of T-DNA inserted genes. Mutants for the genes in same pathway (RDR2 and NRPD1a) did not affect the DNA methylation at reciprocal gene. (D) Proportions of methylated cytosines at different contexts in regions shown in (A–C). Download figure Download PowerPoint In addition to the region downstream of the T-DNA insertion (Figure 6), the upstream regions sometimes, but not always, showed increased non-CG methylation (Supplementary Figure S10), suggesting that function of IBM1 might also be affected by the insertion downstream. Discussion Differentiation of active and inactive chromatin by transcription and histone demethylation In the absence of IBM1 function, thousands of genes accumulated epigenetic marks of silent chromatin, such as H3K9me2 and cytosine methylation at non-CG sites. The IBM1 protein is necessary for excluding these marks from transcribed genes. Although transposons were generally not affected by the ibm1 mutation (Figure 1; Miura et al, 2009), de-repressed transposons behaved like genes in regard to their response to ibm1 (Figure 5; Supplementary Figure S9). Conversely, loss of transcription by T-DNA insertion in the target gene of IBM1 caused induction of non-CG methylation (Figure 6). These results lead us to propose that IBM1 demethylates H3K9me2 from transcribed genes. This would stabilize the expressed state of active genes by preventing accumulation of H3K9me2 (Figure 7A). Figure 7.Proposing model for differential epigenetic modifications between genes and transposons. (A) IBM1 demethylates H3K9 of transcribed genes, which in turn ensure robust transcription, generating the autocatalytic loop. In the absence of IBM1 function, H3K9me2 and mCHG accumulate even in transcribed genes by the self-reinforcing mechanism of KYP and CMT3 (grey arrows). (B) Even in wild type, the self-reinforcing loop of KYP and CMT3 operates in transposons. In transposons, additional backup of H3K9me2 machinery exists that depends on SUVH5 and/or SUVH6. When de-repressed, the autocatalytic loop of transcription and H3K9 demethylation runs even in transposons (grey arrows). The genes and transposons differ in that (i) IBM1 functions only in transcribed genes and that (ii) the KYP alone can function in genes for the methylation of H3K9, but additional SUVH proteins operate in transposons. Download figure Download PowerPoint How could transcription affect IBM1 function? The IBM1 protein may be recruited to transcriptional machinery, as H3K4 and H3K36 methyltransferases are (Li et al, 2002; Krogan et al, 2003; Ng et al, 2003). If so, this recruitment might account for the long-standing enigmatic observations in the mutant of CG methylase gene MET1. The met1 mutation induces not only global loss of CG methylation and de-repression of transposons, but also ectopic methylation at non-CG sites in a subset of genes (Jacobsen and Meyerowitz, 1997; Cokus et al, 2008; Lister et al, 2008). This might be due to depletion of the IBM1 protein in genes, which is caused by de-repression of transposons and recruitment of IBM1 protein there. Consistent with this interpretation, genes showing ectopic CHG methylation in ibm1 and met1 largely overlap (Miura et al, 2009). Interplay between DNA methylation and H3K9me2 The ibm1-induced genic H3K9me2 required CMT3, the DNA methylase for CHG sites (Figure 3C). How does CHG methylation affect H3K9me2? One proposed mechanism posits the direct recognition of DNA methylation by the histone methyltransferase; SRA domains within KYP and other SUVH proteins were recently shown to have methylcytosine-binding activity (Unoki et al, 2004; Johnson et al, 2007; Woo et al, 2007). Interestingly, the SRA domain of KYP binds preferentially to methylated CHG and CHH than methylated CG in vitro (Johnson et al, 2007). If KYP recognizes only CHG and CHH methylation in vivo, it could explain why the cmt3 mutation alone fully suppresses the ibm1-induced H3K9me2 despite the presence of CG methylation. All observations are consistent with the model that genic H3K9me2 and CHG methylation are established by the self-reinforcing mechanism (Figure 7). Unlike the ibm1-induced genic H3K9me2, the H3K9me2 in transposons was not much affected by the cmt3 mutation (Figure 2E–G). Cytosine methylation at context(s) other than CHG may also function upstream of H3K9me2. Interestingly, cmt3 mutation affects H3K9me2 at transposons with low CG methylation (Supplementary Figure S6G–J). This observation is consistent with the interpretation that the remaining H3K9me2 in transposons depends on CG methylation. In fact, mutation in the CG methylase gene MET1 induces a reduction of H3K9me in transposon-rich regions (Soppe et al, 2002; Tariq et al, 2003). Furthermore, siRNA-guided DNA methylation also occurs in transposon sequences (Matzke et al, 2009), which can be excluded from genes by DNA demethylases (Penterman et al, 2007; Zhu et al, 2007). Taken together, H3K9me of transposons seems to be targeted by multiple pathways, resulting in robust