Structural Basis of a Histone H3 Lysine 4 Demethylase Required for Stem Elongation in Rice

Histone lysine methylation is an important epigenetic modification in regulating chromatin structure and gene expression. Histone H3 lysine 4 methylation (H3K4me), which can be in a mono-, di-, or trimethylated state, has been shown to play an important role in gene expression involved in plant developmental control and stress adaptation. However, the resetting mechanism of this epigenetic modification is not yet fully understood. In this work, we identified a JmjC domain-containing protein, JMJ703, as a histone lysine demethylase that specifically reverses all three forms of H3K4me in rice. Loss-of-function mutation of the gene affected stem elongation and plant growth, which may be related to increased expression of cytokinin oxidase genes in the mutant. Analysis of crystal structure of the catalytic core domain (c-JMJ703) of the protein revealed a general structural similarity with mammalian and yeast JMJD2 proteins that are H3K9 and H3K36 demethylases. However, several specific features were observed in the structure of c-JMJ703. Key residues that interact with cofactors Fe(II) and N-oxalylglycine and the methylated H3K4 substrate peptide were identified and were shown to be essential for the demethylase activity in vivo. Several key residues are specifically conserved in known H3K4 demethylases, suggesting that they may be involved in the specificity for H3K4 demethylation.


Introduction
Histone methylation is an important epigenetic modification for chromatin structure, genome function, and gene expression in eukaryotic cells [1]. Histone methylation can be reversed by histone demethylases [1]. Lysine Specific Demethylase 1 (LSD1) is the first identified histone demethylase characterized as a member of the flavin-dependent amine oxidase family [2]. The second class of histone demethylases featured with the jumonji C (JmjC) domain has been shown to catalyze histone lysine demethylation through ferrous ion (Fe(II)) and a-ketoglutaric acid (a-KG)dependent oxidative reactions [3].
Structurally related JmjC domain-containing proteins are classified into 7 subgroups based on phylogenetic analysis of members from yeast and animal cells [4]. Among them, the JmjC domain-containing histone demethylase 1 (JHDM1) subgroup has been demonstrated to reverse mono-and dimethylated histone H3 lysine 36 (H3K36me1/me2) [3]. The JmjC domain-containing histone demethylase 2 (JHDM2) subgroup has been identified to reverse dimethylated histone H3 lysine 9 (H3K9me2) [5], whereas the JMJD2 (also called JHDM3) group specifically reverses methylated H3K9 and/or H3K36 [6]. UTX/UTY is involved in the reversal of histone H3 lysine 27 methylation (H3K27me) [7], whereas PHF8 which belongs to PHF2/PHF8 subfamily demethylates H3K9me1/2, H3K27me2 and monomethylated histone H4 lysine 20 (H4K20me1) [8][9][10]. Aside from these members, the JARID subgroup is responsible for H3K4 demethylation [11]. The JmjC domain-only subgroup has the smallest molecular architecture; members of which have been shown to catalyze divergent reactions, such as histone arginine demethylation [12] and asparagine protein hydroxylation [13]. More recently, JMJD6 is shown to catalyze protein lysyl-hydroxylation [14]. Plant JmjC proteins have been shown to play important roles in the regulation of epigenetic processes, growth and development [15]. Although conserved with yeast and animal homologues, plant JmjC proteins display several distinct features. For instance, the UTX subgroup proteins that exhibit H3K27 demethylase activity are not found in plants. Recent data show that RELATIVE OF EARLY FLOWERING 6 (REF6), a member of the JMJD2 subgroup, can demethylate H3K27 in Arabidopsis [16]. Conversely, there exists a subgroup of plant JmjC proteins that include additional protein modules that are missing from animal or yeast homologues [17]. The amino acid sequence of the JmjC domains of this subgroup is closely related to JARID, while the overall domain organization of the core protein is similar to that of the JMJD2 group [17,18]. To study this plant specific group of JmjC proteins, we analyzed the developmental function, enzymatic activity, and crystal structure of JMJ703, a member of this subgroup in rice. Our results demonstrate that JMJ703 is essential for plant cell division and stem elongation and specifically demethylates mono-, di-, and trimethylated H3K4 in vivo and in vitro. The high resolution structure of the catalytic core of JMJ703 (c-JMJ703) in complex with cofactors and substrate peptide reveals that the overall folding of c-JMJ703 is similar to those defined in animal and yeast JMJD2 that are H3K9 and H3K36 demethylases. However, the crystal structure of c-JMJ703 displays a number of specific features in cofactor interaction and in substrate peptide binding. Substitution mutation analysis indicated that residues implicated in the specific structures are essential for the enzymatic activity, some of which are conserved within JARID proteins and plant H3K4 demethylases and may be involved in the specificity of this class of enzymes.

JMJ703 is essential for rice stem elongation
To investigate the developmental function of JMJ703, we characterized a T-DNA insertion mutant and generated several RNAi lines of the gene ( Figure 1A). Both knockout and knockdown plants were semi-dwarf and produced smaller seeds ( Figure 1B, Figure S1). The phenotypes co-segregated with the T-DNA insertion or the transgene (Table S1). Histological study of stem epidermal tissues revealed no difference in cell length between wild type and the mutant ( Figure 1C), suggesting that the shorter stem phenotype may be caused by a slower cell division rate in the jmj703 mutant plants. Quantitative RT-PCR analysis revealed that the expression of several cell cycle-related genes was unaffected by the mutation (not shown). However, analysis of the cytokinin oxidase (CKX) gene family that reduces active cytokinin levels revealed that several members were highly induced in the young stem of mutant plants ( Figure 1D). Chromatin immunoprecipitation assays revealed that H3K4me3 was clearly increased over the promoter region of the CKX genes in the mutants ( Figure 1D). These data suggested that JMJ703 might regulate H3K4me3 on CKX genes and that the mutant phenotype might be due to cytokinin deficiency caused by increased H3K4me3 and increased expression of CKX genes.
JMJ703 is a histone H3K4 demethylase JMJ703 protein contains several modules, including JmjN, JmjC, C5HC2 zinc finger, FYRN, and FYRC ( Figure 2A). The JmjC domain is critical for histone demethylase activity [3]. To investigate the substrate specificity of JMJ703, the FLAG:HAtagged JmjN-JmjC-zinc finger region (between amino acids 113 and 700, called FA-J3NCZ) of the protein was transiently overexpressed in tobacco leaves for in vivo histone demethylase assays. As shown in Figure 2B, nuclei expressing FA-J3NCZ showed a clear decrease of H3K4me1/2/3 compared to non-transfected nuclei, whereas no difference was observed for H3K27me3. FA-J3NCZ was then affinity-purified from transfected tobacco cells for in vitro histone demethylase assays. Consistent with the in vivo results, FA-J3NCZ could demethylate H3K4me1/2/3 in vitro ( Figure 2C). By contrast, no activity of FA-J3NCZ to demethylate H3K9me1/2/3, or H3K36me1/2/3 was detected either in tobacco cells or in vitro ( Figure S2A-S2B), indicating that JMJ703 is an H3K4-specific demethylase.
Although the structures of a number of animal and yeast JmjC proteins have been defined [19][20][21], the structure of an H3K4 demethylase has not yet been reported. To study the structure of JMJ703, we made many attempts, including N-and C-terminal truncations and limited proteolysis, to obtain JMJ703 crystals. Finally, the fragment spanning amino acids 139 to 498 (termed c-JMJ703), covering the JmjN and JmjC domains, was deemed to be suitable for crystallographic investigation. We first performed surface plasmon resonance (SPR) experiments to measure the binding affinities of c-JMJ703 (J3NC) to H3K4 peptides with mono-, di-, or trimethylation. SPR results indicated that c-JMJ703 (J3NC) bound to H3K4me1, H3K4me2, and H3K4me3 peptides with dissociation constant K d values of 28.9 mM, 19.3 mM, and 30.1 mM, respectively ( Figure 2D). In addition, c-JMJ703 (J3NC) could also bind to H3K4me0 peptide, but with a much lower K d value (76.6 mM), while no binding activity to H3K9me0 and H3K36me0 peptide was observed ( Figure 2D). Meanwhile, the fragment corresponding to the JmjN-JmjC-zinc finger region (J3NCZ) showed a higher affinity to H3K4me1/2/3 peptides with a K d value of 12.6 mM, 15.1 mM and 15.9 mM, respectively ( Figure  S2C). However, neither c-JMJ703 nor J3NCZ bound to H3K9me3 peptide, suggesting that the zinc finger enhances the substrate binding affinity of JMJ703 but is not essential for the binding specificity ( Figure S2D).

Crystal structure features of c-JMJ703
The crystal structure of c-JMJ703 alone or in complex with a-KG (termed as c-JMJ703-a-KG) or with NOG (N-oxalylglycine, a non-catalytic analog of a-KG) and H3K4me3 peptide (termed as c-JMJ703-NOG-H3K4me3) was determined via the molecular replacement method with the modified crystal structure of the core of human JMJD2A (PDB code 2OQ6) [22] as the initial searching model. The final model was best at a resolution of 2.35 Å with a final R work value of 19.3% (R free = 22.5%). The crystals belonged to the P6 3 space group with a slight change in unit cell parameters among the three crystals. One c-JMJ703 molecule was identified in the asymmetric unit with a Matthews coefficient of 2.3 Å 3 /Da (corresponding to 46% solvent content) [23]. There are five solvent-exposed regions in c-JMJ703 structure, including P195-K199, S224-R261, R288-S295, T329-Y349, and Q363-V377, which could not be built due to lack of interpretable electron

Author Summary
Genomic DNA is associated with histone proteins to form the basic structure of chromatin. Lysine residues within the N-terminal end of histones H3 and H4 can be methylated, which may have a positive or a negative effect on the activity of associated DNA or genes, depending on the position of the lysines in the histones. Histone lysine methylation can be reversed by histone demethylases. However, it is not very clear how the specificity of histone demethylases to different histone lysines is determined. In this work we have identified a rice histone demethylase, namely JMJ703, which specifically demethylates methylated histone H3 lysine 4. We found that loss of the enzyme reduces cell division rate of the stem and the size of plant stature, indicating the importance of the protein in plant growth. The crystal structure of the catalytic domain of the protein shares a general similarity with that of mammalian and yeast proteins that demethylate methylated histone H3 lysine 9 and lysine 36, but displays several distinct structural features that are important for substrate and cofactor binding and enzymatic activity of the protein. We found that key amino acids involved in the specific structures are conserved within known H3 lysine 4 demethylases, which may be involved in the specificity to histone H3 lysine 4.
density, suggesting their intrinsic structural flexibility. Moreover, although inter-molecular interactions were found in the c-JMJ703 crystal, gel filtration revealed that c-JMJ703 existed as a monomer in solution ( Figure S3), suggesting that the monomer may be the biological unit.
The c-JMJ703 molecule presented a canonical overall folding of JMJD2 proteins and contained four of the five domains defined in the structure of c-JMJD2A: the JmjN domain (A139-K199), the long b-hairpin (D200-T271), the mixed domain (L272-V377) and the JmjC domain (L378-A498) ( Figure 3A, Figure S4). The JmjC domain, sandwiched by the JmjN and the long b-hairpin with the mixed domain, adopted a jellyroll-like structure with two fourstranded b-sheets as a cupin fold ( Figure 3A) [24].
Although c-JMJ703 shared low primary sequence similarity (less than 25% of the sequence identities) with the reported structural homologues [20,22], the core portion, especially the JmjC domain and the catalytic center of c-JMJ703, presented a topology similar to that of JMJD2 proteins with root mean square derivations of 1.75 Å and 1.81 Å relative to c-JMJD2A and c-Rph1, respectively ( Figure 3B-3E). Nevertheless, c-JMJ703 displayed several significant structural differences. First, the JmjN domain presented a number of distinct features ( Figure 3B). The orientation of the Nterminus of c-JMJ703 was opposite to that of c-JMJD2A and c-Rph1 ( Figure 3B). Second, the long b-hairpin domain contained two b-strands (b2 and b3, aa 209-215 and aa 267-270) and a short a-helix (a4, aa 216-219). The two b-strands of c-JMJ703 were shorter than that of c-JMJD2A and c-Rph1, but had an extra long insertion between them. However, except for the short ahelix (a4), this insertion was mostly invisible in the structure ( Figure 3C, Figure S4A). The short a-helix (a4) represented a sharp difference between c-JMJ703 and c-JMJD2A/c-Rph1. Moreover, the mixed domain of c-JMJ703, composed of several different structural elements, also showed clear structural differences. The first a-helix (a5, aa 272-286) of the mixed domain was two turns longer than that of c-JMJD2A and c-Rph1 ( Figure 3D). The residues Q363 to V377 in the mixed domain of c-JMJ703 were structurally disordered and presented an uninterpretable electron density, whereas the corresponding region of c-JMJD2A is an a-helix. An additional short helix (a7, aa 356-361) was observed in c-JMJ703, whereas a loop was present in c-JMJD2A ( Figure 3D). Several distinct differences could also be observed in the JmjC domain of JMJ703 compared to mammalian and yeast homologues, particularly in the helix-rich region (a8 to a11, aa 419-459) ( Figure 3E, Figure S4A). Residues W381, C392, F437, Q440, L443, H445, L447 and V448 are conserved in and specific to JARID and plant H3K4 demethylases ( Figure 4A). Most of these residues were located in the helix-rich region (a8 to a11) except W381 and C392 ( Figure 4B). A potential extension path of H3K4me peptide toward its C-terminal was proposed based on the locations of these H3K4 demethylase-specific residues ( Figure 4B). Substitution mutations of W381 and L447 abolished the H3K4 demethylase activities (Table 1, Figure S6), indicating that these residues are essential for the demethylation function of the protein.

Catalytic machinery of c-JMJ703
The unambiguous electron density denoted Fe(II) (determined by inductively coupled plasma mass spectrometry) bound to c-JMJ703 and identified the active site within the JmjC domain ( Figure 5). Three key residues, H394, E396, and H482, are perfectly conserved in JMJD2 proteins. They chelated Fe(II) in the active site through their hydrophilic side chains. Fe(II) also interacted with the C-1 carboxyl and C-2 oxo groups of NOG ( Figure 5). Substitution mutations H394A, E396A, and H482A abolished the activity of JMJ703 to demethylate H3K4me1/2/3 ( Table 1, Figure S6), confirming their critical role in the enzymatic activity of the protein.
The electron density indicated the presence of substrate peptide in the c-JMJ703-NOG-H3K4me3 complex (Figure 3, Figure 6). Three of the ten residues in the H3K4me3 peptide (ARTK-me3QTARKS) were visible in the c-JMJ703-NOG-H3K4me3 complex structure (Figure 3, Figure 5, Figure 6). This complex structure provided a view of the substrate peptide-binding mode of an H3K4 demethylase. Compared with the structure of c-JMJD2A in complex with methylated H3K9/H3K36 peptide, the methylated H3K4 and its two flanking residues adopted a different binding conformation to c-JMJ703. H3K9me3/H3K36me3 and their two flanking residues stretched along the long axis of a-KG in c-JMJD2A-a-KG-H3K9me3/c-JMJD2A-NOG-H3K36me3 structure (PDB code 2Q8C/2Q8E, Figure 7B-7C), whereas H3K4me3 and its two flanking residues stretched along the short axis of NOG and was almost perpendicular to H3K9me3/ H3K36me3 and their flanking residues ( Figure 7A).
The methyl group binding pocket of JmjC domain is unique among methylated peptide-binding proteins due to the polar rather than hydrophobic environment [25]. The methyl groups could not be defined properly in the c-JMJ703-NOG-H3K4me3 complex structure because of low occupancy as indicated by obviously higher B-factor for atoms in H3K4me3 peptide. Based on the crystal structure of c-JMJD2A in complex with H3K9me3 (PDB code 2Q8C) [26], four of the five residues that have been shown to be important in methyl group binding in c-JMJD2A [25] were conserved in c-JMJ703, namely, G376, Y383, E396, and N496. S288 in c-JMJD2A (S335 in Rph1) is substituted by an alanine (A494) in c-JMJ703. G376 had no electron density and was invisible in structure, whereas Y383, E396, A494, and N496 were well-defined and adopted very similar conformations as in c-JMJD2A [26]. Substitution mutations of these potential methyl group-binding residues (G376A, Y383A, E396A, N496A) generally impaired the H3K4 demethylase activity of JMJ703 in tobacco cells, with the exceptions of Y383A, which retained a residual activity to demethylate H3K4me2, and N496A, which was still active to demethylate H3K4me2/3 (Table 1, Figure S6).

Distinct conformations of Y321 and N404 in stabilizing a-KG
In the crystal structure of c-JMJ703-NOG-H3K4me3, NOG formed three hydrogen bonds with the side chains of N404 and K412 ( Figure 5) similar to the JMJD2 homologues. K412A mutation also abolished demethylation activity of H3K4 in all three methylation states (Table 1, Figure S6). Unexpectedly, the N404A mutation produced no clear effect on the demethylase activity of JMJ703, whereas the mutation of the counterpart in Rph1 abolished the H3K36me3 demethylase activity [20]. Moreover, NOG interacted indirectly with N490 through two water molecules ( Figure 5), a finding that is not observed in mammalian and yeast homologues. Although most of the residues in a-KG/NOG stabilization showed high structural similarities with its homologues, several features, particularly those of Y321 and N404, distinguished JMJ703 from the JMJD2 members. Y321 and N404 showed significant conformational shifts compared with their counterparts in c-Rph1 and c-JMJD2A ( Figure 8A-8C and Figure S5). Alternative conformation of the Tyr corresponding to Y321 in other JmjC proteins has been observed in JMJD2D (PDB code 3DXU). Interestingly, substitution of Y321 by alanine decreased H3K4me1 demethylase activity but did not affect that of H3K4me2 and H3K4me3 (Table 1, Figure S6), whereas substitution by alanine of the corresponding residue in Rph1 leads to loss of the H3K36me3 demethylase activity of the protein [20]. These data suggested a different NOG-binding mechanism between JMJ703 and Rph1.

Discussion
In this work, we have provided evidence that JMJ703 is a histone demethylase specifically reversing the three forms of methylated H3K4. Arabidopsis JMJ14 that belongs to the same group as JMJ703 also has an H3K4 demethylase activity [27]. About two-thirds of plant genes contain at least one type of H3K4me, among which H3K4me3 is associated with gene activation in Arabidopsis and rice [28,29]. Demethylation of the three forms of H3K4me implies that JMJ703 is important for rice gene expression. The phenotype caused by T-DNA insertion and RNAi and increased expression of CKX genes suggest that JMJ703 may regulate growth hormone metabolism in rice, which is different from JMJ14 that has functions in floral transition, DNA methylation, and RNA silencing in Arabidopsis [27,[30][31][32][33]. H3K4me3 on specific plant genes can be induced by external signals, such as light, drought, and submergence, among others [34][35][36]. H3K4me3 is suggested to be a mark of active genes and may play a role in plant adaptation to environmental cues. Whether JMJ703 is involved in resetting stress-induced gene expression requires further analysis.
Although JmjC domain of JMJ703 is closely related to that of the JARID proteins [17], the domain organization of the catalytic core of the enzyme (JmjN-long b-hairpin-mixed domain-JmjC) resembles that of JMJD2 members (c-JMJD2A and c-Rph1). However, our functional analysis has revealed that JMJ703 (i.e. FA-J3NCZ) was unable to demethylate H3K9me1/2/3 or H3K36me1/2/3 ( Figure S2A-S2B). This suggests that the enzymatic specificity may be determined essentially in the JmjC domain. Our data indicate that the helix-rich region (a8-a11, aa 419-459) of the JmjC domain of JMJ703 is substantially divergent from that of JMJD2A and Rph1. For instance, K241 of c-JMJD2A, which is conserved in c-Rph1, is considered to be a key catalytic residue that recruits and positions O 2 between Fe(II) and the methyl group to participate in the reaction [25]. Substitution mutation of K241L abolishes the activity of c-JMJD2A [25]. However, this position is replaced by a leucine in JMJ703 (L447) and in other plant homologues and JARID proteins ( Figure 4). L447 was found to be essential for the H3K4 demethylase activity of JMJ703 (Table 1, Figure S6). This leucine together with the other residues that are specifically conserved in JARID/H3K4 demethylases and mostly found in the helix-rich region (a8-a11, aa 419-459) of c-JMJ703 (Table 1, Figure 4 and Figure S6), may be involved in the specificity of H3K4 demethylases.
Our data indicate that although most of the JMJ703 residues in a-KG/NOG stabilization showed a high structural similarity with JMJD2 proteins, the residues Y321 and N404 displayed specific features (Figure 8, Figure S5). The conformation of Y321 in c-JMJ703 indicates that the residue may not be in direct contact with a-KG or NOG. Substitution mutation of Y321 only affected H3K4me1 demethylase activity. In addition, N404A substitution did not clearly alter the H3K4 demethylase activity.
The methyl group binding pocket in c-JMJ703 differs also from that of c-JMJD2A and Rph1 in one key residue (A494 in JMJ703 compared to S288 in JMJD2A or S355 in Rph1). S288 of JMJD2A has been shown to be important to reverse di-and trimethylated H3K9 [25,26,37]. However, JMJD2D, which harbors an alanine in this position, is 7-fold and 60-fold more efficient than JMJD2A in the demethylation of H3K9me3 and H3K9me2, respectively. JMJD2D has residual activity to H3K9me1 that is absent from JMJD2A [26]. The presence of alanine at this position in JMJ703 (A494) supports the demethylation function of the protein of all three forms of methylated H3K4 (Figure 2). Substitution of this alanine by a serine abolished the demethylase activity of H3K4me1/2, but retained the activity toward H3K4me3 (Table 1, Figure S6). Interestingly, this alanine residue is conserved in most of the plant JMJD2 members ( Figure S4), suggesting that these plant proteins may be able to demethylate the three methylated forms of relevant histone lysines.
In addition, the zinc finger that overlaps partially with the JmjC domain of JMJD2A is not found in JMJ703. Although the residues involved in the zinc finger formation in JMJD2A are conserved in Rph1, no zinc ion is observed in its structure and this domain is disordered and invisible in the c-Rph1 structure [20]. By contrast, JMJ703 contains a distal zinc finger motif relative to the JmjC domain (about 100 amino acids downstream the JmjC domain), which is present in Rph1, but absent from JMJD2A. The distal zinc finger is found in a number of other JmjC proteins and is shown to be involved in nuclear localization of JMJ706 [17]. These observations together with our data showing that the inclusion of the motif increased the substrate peptide (H3K4me1/ 2/3) binding activity to two-folds but not alter the binding specificity, suggest that this zinc finger may not be a component of the catalytic core of JmjC proteins, but is important for their enzymatic activity. Our data also revealed a number of other clear structural differences in the plant JmjC protein. The N-terminal end of c-JMJ703 adopted a completely different conformation from that of c-JMJD2A and c-Rph1. The long b-hairpin domain and mixed domain of c-JMJ703 also showed significant divergence at the amino acid sequence level compared to c-JMJD2A and c-Rph1 ( Figure S4). However, the divergent regions are mostly invisible in c-JMJ703 and c-Rph1, indicating intrinsic disorder.
In summary, we have shown that JMJ703 is an H3K4 demethylase that is important for rice plant development. Although JMJ703 displays an overall structure similar to that of JMJD2 H3K9 and H3K36 demethylases, it exhibits a number of important structural features which are conserved within JARID and plant H3K4 demethylases. In particular, our data suggest that the helix-rich region (a8-a11, aa 419-459) of the catalytic JmjC domain of JMJ703 may be involved in the determination of substrate specificity of the H3K4 demethylase.

RNAi vector construction and rice transformation
A cDNA fragment between nucleotides 2282 bp and 2834 bp relative to the translation start site was amplified using primers RiJ3-F and RiJ3-R and then cloned into T vector (Promega). The fragment was cloned into double strand RNAi vector pDS1301 [38]. Rice transformation with callus generated from the rice cultivar Zhonghua11 (ZH11) was performed as previously described [38].

Histone demethylation assay
The cDNA fragment of the J3NCZ region was amplified using primers FAJ3NCZ-F and FAJ3NCZ-R and then cloned into pFA121 vector modified from pBI121 (GenBank: AF485783.1), in which GUS was replaced by 26FLAG:26HA tag. The pFA121-J3NCZ plasmid was transferred into Agrobacterium tumefaciens strain EHA105 cells for tobacco infection with the help of Agrobacterium cells harboring P19 plasmids [39]. For in vivo histone demethylation assay, tobacco infection and nuclei isolation were performed as previously described [27]. The immunostaining protocol was modified from http://sites.bio.indiana.edu/,pikaardlab. Briefly, the nucleus solution was placed on a poly-lysine coated slide, air dried, and then refixed with 4% formaldehyde in KPBS (1.28 M NaCl, 20 mM KCl, 80 mM Na 2 HPO 4 , 20 mM KH 2 PO 4 , pH7.2) containing 1% Triton X-100 for 20 min. After 3 washes with KPBS-1% Triton, slides were blocked with blocking solution (1% BSA in KPBS-1% Triton) and incubated at 37uC for 30 min. The blocking solution was washed off with KPBS-1% Triton and the slides were incubated with both anti-HA and histone méthylation specific antibodies diluted in blocking solution to 1:300 at 4uC overnight. After washing with KPBS, slides were blocked and incubated with fluorescent-labeled antibodies at 37uC for 2 h. After washing with KPBS thrice, nuclei were stained with 5 mg/ mL 4,6-diamidino-2-phenylindole (DAPI), coated with a drop of VECTASHIELD Mounting Medium (H-1000, Vector Laboratories), and then covered with a coverslip followed by confocal microscope (Leica) detection. In vitro histone demethylation assays were performed as previously described [6]. FA-J3NCZ was purified after overexpression with Anti-FLAG M2 magnetic beads according the manufacturer's instruction (Sigma).

Chromatin immunoprecipitation
Chromatin immunoprecipitation assays were performed as previously described [38]. Briefly, 2 g of 10 day-old rice seedlings were fixed with 1% formaldehyde, then chromatin was extracted and immunoprecipitated with anti-H3K4me3 (ab8580, Abcam) and anti-H3K4me2 (04-790, Millipore) antibodies. After reversing crosslink and protease K treatment, DNA was recovered for realtime PCR analysis with primers listed in Table S2.
Protein production for crystallization DNA fragment encoding amino acids 139 to 498 of wild type JMJ703 was amplified by PCR, ligated into pET-28a expression vector using BamHI and XhoI and transformed into Escherichia coli Rosetta (DE3) cells. Transformed E. coli cells were cultured at 37uC in LB medium containing 50 mg/L kanamycin until the OD 600 nm reached 0.8. The culture was cooled to 16uC and subsequently induced with 0.5 mM isopropyl b-D-1-thiogalactopyranoside. Cells were harvested after overnight induction by centrifugation at 5,000 g for 10 min at 4uC. The cell pellets were resuspended in lysis buffer containing 50 mM Tris (pH 8.0), 500 mM NaCl, and 10% (v/v) glycerol and then disrupted by an ultra-high pressure cell disrupter (JNBIO, Guangzhou, China) at low temperature. Cell debris was removed by centrifugation at 25,000 g for 30 min at 4uC. The 66His tagged protein was purified by Ni-NTA affinity chromatography, cleaved with thrombin (Sigma) overnight at 4uC, and eluted with lysis buffer. The purified c-JMJ703 protein was concentrated and further purified by Hitrap Q (GE Healthcare) anion-exchange chromatography using a 0.05 M to 1 M NaCl gradient in 25 mM HEPES (pH 7.5). The target protein was confirmed to have purity over 95% by SDS-PAGE and concentrated to 10 mg/mL before crystal growth or storage.

Crystallization
Crystallization was performed at 18uC by the hanging-drop vapor-diffusion technique. Crystals were obtained by mixing 1 mL of the protein solution with an equal volume of a reservoir solution; the mixture drop was equilibrated against 500 mL of the reservoir solution. Crystals were obtained with a reservoir solution containing 0.05 M potassium phosphate monobasic, 25% (w/v) polyethylene glycol 8,000 and reached final dimensions of  506506100 mm3 within 2 days. Crystals of apo c-JMJ703 were initially obtained but showed poor diffraction at 3.0 Å resolution and poor reproducibility. Crystals of c-JMJ703 in complex with NOG displayed much better diffraction quality at 2.3 Å , indicating that NOG helps stabilize the conformation of the c-JMJ703 polypeptide. Crystals of the c-JMJ703 substrate complex were obtained by co-crystalizing c-JMJ703 with 10 mM peptide. Crystals were cryo-protected by soaking in a cryo-protectant consisting of the reservoir solution with additional 15% (v/v) glycol. Cryo-protected crystals were then flash-cooled in liquid nitrogen and transferred into a dry nitrogen stream at 100 K for X-ray data collection.
X-Ray data collection, processing, and structure determination Diffraction data for apo c-JMJ703 crystals were collected at a resolution of 3.0 Å at 100 K using a MAResearch M165 CCD detector in beamline 1W2A at the Beijing Synchrotron Radiation Facility. The data set for c-JMJ703-a-KG and c-JMJ703-NOG-H3K4me3 were collected at 2.3 Å and 2.4 Å in beamline BL17A (Photon Factory, Japan) and BL17U1 (SSRF, Shanghai, China), respectively, with an ADSC Q315 CCD detector at the wavelength of 1.0000 Å . All data sets were indexed, integrated, and scaled using the HKL2000 package [40]. The crystals belonged to space group P6 3 with cell parameters a = 53.0 Å , b = 86.1 Å , c = 87.3 Å , and b = 90.7u. All crystals showed very small shifts in their cell parameters.
The molecular replacement method was used to calculate the phases using the PHASER program [41] and the modified crystal structure of a catalytic core of the human c-JMJD2A (PDB code 2OQ6) as the initial searching model. Manual model building and refinement were performed with the programs COOT [42] and PHENIX [43]. Solvent molecules were located from stereochemically reasonable peaks in the sA-weighted 2Fo-Fc difference Fourier electron density map (1.2 s). Model geometry was verified using the program PROCHECK (Table S3) [44]. Coordinates have been deposited in PDB with accession 4IGP, 4IGO and 4IGQ for apo c-JMJ703, c-JMJ703-a-KG and c-JMJ703-NOG-H3K4me3, respectively.

Surface plasmon resonance binding analysis
Analyses were carried out at 25uC with the BIAcore 3000 system. Ten mg/mL c-JMJ703 in 10 mM sodium acetate buffer (pH 5.0) was covalently coupled to a CM5 chip (Biacore) using an Amine Coupling Kit (Biacore) according to the manufacturer's instructions. Using HEPES buffer (10 mM HEPES pH 7.5, 150 mM NaCl, 0.005% Tween20), various concentrations of the H3K4me1/2/3 peptide were injected through a flow cell that was not activated and then through another flow cell containing c-JMJ703 at a rate of 10 mL/min for 2 min. The c-JMJ703 surface was regenerated between two injections by running 15 mL of 5 mM NaOH twice through the flow cell at 30 mL/min. The sensorgram obtained for the inactivated flow cell was subtracted to correct for nonspecific binding and the bulk signal from the peptide in solution. Data were analyzed using BIAevaluation 4.1 software.   Figure S4 Key catalytic core sequence comparison between JMJ703, structurally studied JmjC proteins and other plant JmjC proteins. (A). Sequence alignment of c-JMJ703, c-JMJ16, c-JMJD2A and c-Rph1. The secondary structure of c-JMJ703 is shown and labeled above the alignment. The color scheme is the same as Figure 3. Residues involved in Fe(II), a-KG and methyl group binding and previously proposed O 2 -recruiting lysine are shown beneath the alignment in brown triangle, green square, light blue circle and red circle, respectively. (B). Sequence alignment of the rice (Jmj701-720) and Arabidopsis JmjC (Jmj11-32) proteins with representative animal/yeast proteins. JmjC proteins for which the crystal structure has been determined are highlighted by blue. Only key residue regions are shown. Key residues identified in the structure of JMJ703 are indicated by stars at the bottom, with residues involved in Fe(II) binding in purple, a-KG binding in brown, methyl group binding in blue and the previously proposed O 2 -recruiting leucine in yellow. (TIF) Figure S5 Electron density of bound a-KG/NOG and interacting residues of c-JMJ703. The residues interact with bound a-KG (A) and NOG (B) in the complex structure of c-JMJ703-a-KG and c-JMJ703-NOG-H3K4me3 are shown as blue and green sticks, respectively. Fe(II) and solvent molecules, which mediate interaction between polypeptide and compounds, are presented by colored spheres. All components are covered by electron density (2Fo-Fc map at 1.1s). (TIF) Figure S6 Demethylation assays for JMJ703 substitution mutants that have been produced based on the structural data. The mutant names are labeled at the top left corner of each panel. Image panels from left to right are staining by DAPI, anti-HA and anti-methylated histones, respectively. Bar = 10 mm. At least 30 nuclei that expressed JMJ703 per transfection were observed and imaged. (TIF)