Correction: An AP Endonuclease Functions in Active DNA Demethylation and Gene Imprinting in Arabidopsis

[This corrects the article DOI: 10.1371/journal.pgen.1004905.].


Introduction
DNA methylation is a stable epigenetic mark that regulates numerous aspects of the genome, including transposon silencing and gene expression [1][2][3][4][5][6][7]. In plants, DNA methylation can occur within CG, CHG, and CHH motifs (H represents A, T, or C). Genome-wide mapping of DNA methylation in Arabidopsis has revealed that methylation in gene bodies is predominantly at CG context whereas methylation in transposon-and other repeatenriched heterochromatin regions can be within all three motifs [8]. Although the function of abundant CG methylation within genic regions remains unclear, DNA methylation generally correlates with histone modifications that repress transcription activities [1,9,10]. DNA methylation patterns are coordinately controlled by methylation and demethylation reactions. In Arabidopsis, symmetric CG and CHG methylation can be maintained by DNA METHYLTRANSFERASE 1 (MET1) and CHROMOMETHYLASE 3 (CMT3), respectively, during DNA replication. In contrast, asymmetric CHH methylation cannot be maintained and is established de novo by DOMAINS REAR-RANGED METHYLASE 2 (DRM2), which can be targeted to specific sequences by the RNA-directed DNA methylation (RdDM) pathway [1,10,11]. DNA methylation is antagonized by an active DNA demethylation pathway that includes the DNA glycosylases REPRESSOR OF SILENCING1 (ROS1), DEME-TER (DME), DEMETER-LIKE2 (DML2) and DEMETER-LIKE3 (DML3) [12][13][14].
ROS1, DME, DML2 and DML3 are all bifunctional DNA glycosylases that initiate active DNA demethylation by removing the 5-methylcytosine (5-meC) base and subsequently cleaving the phosphodiester backbone by either bor b, d-elimination [12,[14][15][16]. When b, d-elimination occurs, a gap with a 39phosphate group is generated. Our previous work demonstrated that the 39 DNA phosphatase ZDP catalyzes the conversion of 39-phosphate group to a 39-hydroxyl (39-OH), enabling DNA polymerase and ligase activities to fill in the gap [17]. The belimination product, a gap with a blocking 39-phosphor-a, bunsaturated aldehyde , also must be converted to a 39-OH to allow completion of the demethylation process through single-nucleotide insertion or long patch DNA synthesis by DNA polymerase and ligase [18]. However, the enzymes that may function downstream of ROS1 and DME in the b-elimination pathway have not been identified.
The mutation of ROS1 leads to hypermethylation and transcriptional silencing of a luciferase reporter gene driven by the RD29A promoter, as well as of the endogenous RD29A gene [13]. ROS1 dysfunction also causes DNA hypermethylation in thousands of endogenous genomic regions [19]. zdp mutants also show hypermethylation in the RD29A promoter and many endogenous loci. However, the hypermethylation in the RD29A promoter caused by zdp mutations is not as high as that caused by ros1 mutations, and there are many ROS1 targets that are not hypermethylated in zdp mutants [17]. These observations suggest that there may be an alternative, ZDP-independent branch of the DNA demethylation pathway downstream of ROS1 and other DNA glycosylases/lyases. Although ROS1 functions in almost all plant tissues [13], DME is preferentially expressed in the central cell of the female gametophyte and is important for the regulation of gene imprinting in the endosperm [20][21][22]. In Arabidopsis, the imprinted protein-coding genes include FWA (Flowering Wageningen), MEA (MEDEA) and FIS2 (Fertilization-Independent Seed 2) and the list is expanding [21][22][23][24][25]. The loss-of-function mutation of DME results in aberrant endosperm and embryo development because of DNA hypermethylation and down-regulation of the maternal alleles of imprinted genes [26]. DME is also necessary for DNA demethylation in the companion cells in the male gametophyte [27][28][29]. SSRP1, a chromatin remodeling protein, was identified as another factor required for gene imprinting and the mutation of SSRP1 gives rise to a maternal lethality phenotype similar to that caused by DME mutations [30]. Therefore, it is possible that ZDP and other protein(s) acting downstream of the 5-meC DNA glycosylases/lyases may also affect gene imprinting in Arabidopsis. Intriguingly though, neither ZDP mutants nor mutants in other DNA repair enzymes that may be downstream of DNA glycosylases/lyases show developmental phenotypes associated with defective gene imprinting.
In this study, we characterized the functions of Arabidopsis APE-like proteins in the processing of 39-blocking ends generated by ROS1 and examined methylome changes induced by ape mutations. We found that purified APE1L can process 39-PUA termini to generate 39-OH ends. APE1L also displays a weak activity in converting 39-phosphate termini to 39-OH ends. ape1l-1 mutants show altered methylation patterns in thousands of genomic regions. Interestingly, we found that the ape1l +/2 zdp 2/2 mutant is maternally lethal, giving rise to a seed abortion phenotype resembling that of dme mutants. The maternal alleles of the imprinted genes FWA and MEA are hypermethylated, and their expression levels are reduced in the endosperm of such abnormal seeds of the double mutant. Thus, APE1L functions downstream of the ROS1/DME subfamily of DNA glycosylases/ lyases in active DNA demethylation and genomic imprinting in Arabidopsis.

Results
APE1L possesses a potent activity against 39-PUA termini generated by ROS1 The Arabidopsis genome encodes three AP endonuclease-like proteins: APE1L, APE2 and ARP [31]. We purified recombinant full-length APE1L, APE2 and ARP proteins, and found that all three enzymes exhibit AP endonuclease activity in vitro. APE1L, but not APE2 or ARP, also displayed a 39-phosphatase activity (Fig. 1A). We next wanted to determine if these proteins can process the 39-PUA termini generated by ROS1 after the belimination reaction. We first incubated ROS1 with a 51-mer duplex DNA substrate containing a 5-meC residue at position 29 in the 59-end labeled strand (Fig. 1B). As expected, the DNA glycosylase/lyase activity of ROS1 generated a mixture of band b, d -elimination products, with either 39-PUA or 39-phosphate ends, respectively (Fig. 1B, lane 1). These products were then purified and combined with either APE1L, APE2 or ARP proteins. We found that APE1L efficiently processed the 39-PUA to generate a 39-OH terminus. In comparison, 10-fold higher amounts of APE2 or ARP proteins displayed either weak [32] or undetectable (APE2) activity against 39-PUA ends (Fig. 1B).
To confirm that APE1L is responsible for the detected enzymatic activity we generated an APE1L mutant, N224D. Residue N224 corresponds to N212 of human APE1, which is essential for the enzymatic activity of the mammalian protein [33]. Substitution of N224 by aspartic acid almost completely abolished the activity of APE1L on the 39-PUA termini (Fig. 1C). The mutation also greatly reduced the AP endonuclease activity on a synthetic AP site and the 39 phosphatase activity on 39-phosphate ends. Altogether these results indicate that, in addition to its AP endonuclease activity, APE1L possesses a potent 39-phosphodiesterase activity that can efficiently process the 39-PUA blocking ends generated by ROS1.
APE1L is able to process 39-PUA and 39-phosphate termini in the presence of ROS1, but does not increase the enzymatic turnover of ROS1 ROS1 remains bound to its reaction products, which contributes at least partially to the highly distributive behavior of the enzyme in vitro [34]. To determine whether APE1L is able to process 39-PUA and/or 39-phosphate termini in the presence of ROS1, we incubated ROS1 and a duplex DNA substrate containing a single 5-meC residue, with WT or N224D APE1L ( Fig. 2A). We found that a 39-OH terminus is efficiently generated in the presence of WT but not mutant APE1L ( Fig. 2A). The emergence of the 39-OH terminus is concomitant with the loss of Author Summary DNA cytosine methylation (5-methylcytosine, 5-meC) is an important epigenetic mark, and methylation patterns are coordinately controlled by methylation and demethylation reactions during development and reproduction. In plants, REPRESSOR OF SILENCING (ROS1) is one of the well characterized 5-meC DNA glycosylases that initiate active DNA demethylation by 5-meC excision. Our previous work showed that a 39-DNA phosphatase, ZDP, functions downstream of ROS1 during active DNA demethylation in Arabidopsis. Here we found that the apurinic/apyrimidinic endonuclease APE1L functions downstream of ROS1 in a ZDP-independent branch of the active DNA demethylation pathway in Arabidopsis. In plants, gene imprinting requires the 5-meC DNA glycosylase Demeter (DME) that has been proposed to initiate a base excision repair pathway for active DNA demethylation in the central cell in female gametophyte. However, besides DME, no other base excision repair enzymes have been found to be important for gene imprinting. Our results show that APE1L and ZDP act jointly downstream of DME to regulate gene imprinting in plants, and suggest that DME-initiated active DNA demethylation in the central cell and endosperm uses both APE-and ZDP-dependent mechanisms.
both 39-PUA and 39-phosphate ends, suggesting that the 39-OH terminus is produced by the 39-phosphodiesterase activity of APE1L on the 39-blocking ends generated by ROS1. Quantification of the reaction products revealed that the total amount of strand incision is not increased in the presence of APE1L (Fig. 2B).
To assess whether APE1L modulates the DNA glycosylase/lyase activity of ROS1, we performed the reaction in the absence of Mg 2+ , which is required for APE1L but not ROS1 activity. We found that the enzymatic activity of ROS1 is not increased in the presence of APE1L (S1 Fig.). Thus, APE1L is able to access the Reactions were stopped at 2 or 4 h, as indicated, and products were separated using a 12% denaturing polyacrylamide gel and detected by fluorescence scanning. doi:10.1371/journal.pgen.1004905.g001 39-blocked termini generated by ROS1 but does not increase the turnover of this DNA glycosylase. These results suggest that APE1L does not displace ROS1 from DNA.
APE1L and ROS1 interact in vitro and in vivo and form a ternary complex with a gapped DNA substrate We next used in vitro pull-down assays to test whether ROS1 and APE1L can physically interact (Fig. 3A). His-tagged ROS1 (His-ROS1) was incubated with either Maltose Binding Protein (MBP) or MBP-APE1L bound to an amylose column. We found that MBP-APE1L, but not MBP, associates with His-ROS1, suggesting that APE1L and ROS1 directly interact in vitro.
To gain insights into the transfer of DNA demethylation intermediates between ROS1 and APE1L, we performed electrophoretic mobility shift assays with a gapped DNA substrate ( Fig. 3B-C). MBP-APE1L alone is not able to form a stable complex with the substrate, judging by the smeared band next to the position of the free probe (Fig. 3B, lanes 2 and 1). A mobility shift was observed when the DNA substrate was incubated with His-ROS1, consistent with complex formation (Fig. 3B, lane 3). As we have previously reported [17], part of the labeled probe remained trapped in the wells, hinting at the formation of insoluble His-ROS1-DNA complexes. Next, we incubated the gapped DNA substrate and His-ROS1 with increasing concentrations of MBP-APE1L to assess complex formation. With increasing MBP-APE1L, the band corresponding to the ROS1-DNA complex and the labeled material in the well gradually disappeared, concomitant with the appearance of a discrete, new band (Fig. 3B, lanes 4-7). Importantly, this band was only detected when both ROS1 and APE1L were present in the binding reaction. These results suggest that ROS1, APE1L, and gapped DNA form a ternary complex and that ROS1 is required for APE1L to stably associate with the DNA substrate.
To further examine complex formation we performed supershift experiments using antibodies against MBP-APE1L and His-ROS1 (Fig. 3C). We found that adding anti-MBP to a binding reaction containing MBP-APE1L, His-ROS1 and DNA generated an additional shift, thus confirming the presence of APE1L in the complex. However, a supershift was not observed in the presence of the anti-His antibody. We reasoned that access to the His epitope on His-ROS1 might be restricted in the complex. Therefore, as an alternative approach, we compared the mobility shifts generated from binding reactions containing the DNA gapped substrate, MBP-APE1L and either His-ROS1 or MPB-ROS1 (Fig. 3D). We found that MBP-ROS1, which has a higher molecular weight than His-ROS1, gave rise to a higher molecular weight gel shift, thus confirming that ROS1 is also present in the complex. The most likely interpretation for these results is that ROS1, APE1L, and the gapped DNA substrate form a ternary complex. To further confirm the interaction between APE1L and ROS1, we performed a firefly luciferase complementation imaging assay [35] in tobacco leaves. We found that APE1L can interact with ROS1 in the tobacco leaves (Fig. 3E).

APE1L co-localizes with ROS1 in vivo
Our previous data show that ZDP, a component of the active DNA demethylation pathway, co-localizes with ROS1 in subnuclear foci [17]. To determine the subnuclear localization of APE1L protein, we generated antibodies specific to APE1L and used them for immunolocalization of APE1L in Arabidopsis leaf nuclei. As shown in Fig. 4A, APE1L is broadly distributed throughout the nucleus. In 62% of the cells examined, APE1L is enriched in the nucleolus whereas in 38% of the cells, APE1L localizes to small nucleoplasmic foci. Only very weak signals were observed when the antibodies were applied to nuclei preparations of ape1l-1 mutant plants, indicating that the staining patterns in wild type plants reflect APE1L localization rather than non-specific binding of the antibody (Fig. 4A).
To test whether APE1L co-localizes with ROS1 or ZDP, we performed co-immunofluorescence. In our experiments, FLAGtagged ROS1 was expressed from its native promoter in ros1-1 mutants and visualized with anti-FLAG antibodies. We observed APE1L co-localization with ROS1 within both nucleoplasmic foci and also the nucleolus in about 10% of cells, as shown by the strong yellow signals (Fig. 4B). In 54% of the cells, APE1L colocalizes with ROS1 in the nucleolus but not in nucleoplasmic foci, whereas in 36% of the cells, APE1L and ROS1 do not substantially co-localize (Fig. 4B). APE1L and ZDP also  co-localize in nucleoplasmic foci in approximately 28% of cells (Fig. 4C). Thus, APE1L co-localizes with components of the DNA demethylase machinery in distinct subnuclear structures in a subset of cells.

APE1L dysfunction causes genome-wide alterations in DNA methylation
To evaluate possible roles of APE1L in active DNA demethylation initiated by the ROS1 subfamily of DNA glycosylases/ lyases, two T-DNA insertion lines were isolated for APE1L (S2A Fig.). RT-PCR analysis with APE1L-specific primers corresponding to the full-length open reading frame of the gene detected the expected product in wild-type plants in both the Ws and Col backgrounds, but not in ape1l-1, which is in Ws. In contrast, ape1l-2 shows almost the same expression level as wild-type plants (S2B Fig.). Since the endonuclease ARP shows weak activity against the 39-PUA blocking ends generated by ROS1 in vitro, we also isolated two T-DNA insertion lines for ARP (S2A Fig.) and confirmed by RT-PCR that they have a complete loss of mRNA expression (S2B Fig.). One of the mutants, arp-1, was used for further experiments. To examine the general DNA methylation status in the ape1l-1 and arp-1 mutants, we compared the susceptibility of 5S rDNA and 180-bp centromeric repeat regions to the restriction enzymes HpaII and MspI. These enzymes recognize the same site (CCGG), but HpaII cleavage is methylation-inhibited whereas methylation does not affect cleavage by MspI. DNA cleavage was assessed by Southern analysis. Similar to the zdp-1 and ros1-4 mutations, the ape1l-1 or arp-1 mutation does not affect the DNA methylation levels at the 5S rDNA or 180-bp centromeric repeats (S3 Fig.), suggesting that the ape1l-1 and arp-1 mutants do not have changes in their overall DNA methylation patterns.
We performed whole genome bisulfite sequencing using DNA from 14-day-old ape1l-1, arp-1, zdp-1 and their corresponding wild-type control plants. The CG methylation levels in wild type (Col-0) and zdp-1 mutant are similar, but the CHG and CHH levels are mildly elevated in zdp-1 (S4A Fig.). For ape1l-1, its overall genome methylation level in CG, CHG and CHH contexts is slightly higher than that in Ws (S4B Fig.). In total, we identified 6389 DMRs (differentially methylated regions) in ape1l-1 mutant plants, including 3497 hyper-DMRs that have a significant increase in methylation and 2892 hypo-DMRs that have a significant reduction in methylation (S4C Fig.; S3 Table). In contrast, arp-1 only affects methylation levels at 403 genomic regions, including 162 hyper-DMRs and 241 hypo-DMRs (S4C Fig.). 1559 hyper-DMRs and 612 hypo-DMRs were identified from zdp-1 (S4C Fig.; S4 Table). The hyper-DMRs and hypo-DMRs identified in ape1l-1, zdp-1 and ros1-4 are evenly distributed along the five chromosomes (S4D Fig.). To determine whether APE1L and ZDP mutations affected DNA demethylation in specific genomic regions, we analyzed intergenic regions, transposable elements (TEs) outside of genes, TEs overlapping with genes and genic regions. Unlike zdp-1, ros-1 mutants and ros1-3;dml2-1;dml3-1 (rdd) triple mutants, which have less than 43% of hypermethylated (hyper-) DMRs distributed in gene regions, in ape1l-1 and arp-1 more than 60% of the hyper-DMRs are distributed in gene regions (S4E Fig.). In contrast, the percentages of hyper-DMRs distributed in TEs in ape1l-1 and arp-1 are lower than those in zdp-1, ros-1 and rdd mutants (S4E Fig.). These data indicate that the APE1L and ARP mutations preferentially impact DNA demethylation of gene regions while the ZDP and ROS1 mutations have a greater impact on TE regions. The distribution patterns of classified hypo-DMRs are different from those of hyper-DMRs. The percentages of hypo-DMRs in gene regions are higher than 70% in rdd, zdp-1 and ape1l-1. The arp-1 has a low percentage of hypo-DMRs in gene regions but high percentage of hypo-DMRs in intergenic regions (S4E Fig.). The ape1l-1 mutation affects CHG and CHH demethylation more profoundly than CG demethylation, both in gene regions or in TEs (S5A Fig.). We also examined the effect of APE1L mutation on TEs of different lengths and found that the ape1l-1 mutation has a bigger impact on shorter genes but longer TEs (S5B and S5C Figs.). Unlike ape1l-1, zdp-1 shows almost the same DNA methylation pattern for both gene regions and TEs (S5 Fig.).
Compared to the high level of overlap (70.9%) between zdp-1 and rdd hyper-DMRs, less than 50% of the hyper-DMRs in ape1l overlap with those in rdd (S4C and S5D-S5E Figs.). One reason for this relatively low level of overlap may be the difference in genetic backgrounds; the ape1l-1 mutant is in the Ws background whereas the other mutants are in Col. When the hyper-DMRs in ros1-1 (C24 background) and ros1-4 (Col-0 background) were compared, the overlap was also quite low (52%). For the hyper-DMRs, the level of overlap between ape1l-1 and zdp-1 is also very low (14%) (S5F Fig.) even though some loci do show hypermethylation in ape1l-1 as well as zdp-1 (S5G-S5J Fig.). These results are consistent with the notion that APE1L and ZDP largely represent two different mechanisms (AP endonuclease vs 39-phosphatase) downstream of the DNA glycosylases/lyases, despite their redundant functions (as 39-phosphatases).
ROS1 and ZDP mRNA levels are decreased in RdDM pathway mutants [17]. We examined the expression of APE1L and ARP in the RdDM mutants nrpd1-3 and nrpe1-11, and found no substantial decreases in the mRNA levels in the mutants compared to Col (S6A Fig.). We also measured the expression levels of ROS1 and ZDP in ape1l-1 and arp-1 mutants, and found that the expression levels are similar in the mutants compared to those in the Col or Ws wild type control plants (S6B Fig.). Also, unlike the zdp-1 mutant, which is hypersensitive to MMS induced DNA damage, the ape1l-1 and arp-1 mutants show a sensitivity level similar to that of wild-type plants (S7A Fig.).

The double mutant of APE1L and ZDP are embryonic lethal
To study the potential genetic interactions between APE1L and ZDP, we crossed ape1l-1 and zdp-1 mutant plants. Interestingly, we found that ape1l +/2 zdp 2/2 and ape1l 2/2 zdp +/2 plants produce many aborted seeds, suggesting that the double mutations of APE1L and ZDP are lethal (Fig. 5A). We grew the viable seeds, genotyped the seedlings, and found no ape1l 2/2 zdp 2/2 plants ( Table 1). The ratio of aborted seeds is 48.7% in self-pollinated ape1l +/2 zdp 2/2 plants and 26.5% in self-pollinated ape1l 2/2 zdp +/2 plants (S5 Table). Approximately seven days after pollination, the seeds fated to abortion show white color and plump phenotypes (S8A Fig.). The endosperm in those seeds fails to undergo cellularization and the growth of their embryos is arrested (S8B Fig.). Later, those seeds accumulate brown pigments and collapse.
We examined the morphology of aborting seeds from ape1l +/2 zdp 2/2 and ape1l 2/2 zdp +/2 mutants using differential interference contrast microscopy. The major defects of aborting seeds are arrested embryo growth at the heart stage or earlier (S8B Fig.) and abnormal sizes of endosperm nuclei (S8C Fig.). In some aborting seeds, the embryos are invisible, indicating that the embryos are arrested very early in development. The aborting seeds of ape1l +/2 zdp 2/2 and dme +/2 both display arrested embryo growth. Unlike ape1l +/2 zdp 2/2 mutant seeds, dme +/2 mutant seeds display clumps of unknown structures but there were no aberrant endosperm nuclei (S8B-C Fig.).
We noticed that the ape1l +/2 zdp 2/2 mutant has abnormal segregation ratio (4.07:1:0), which does not fit the expected segregation ratio of maternally lethal plant (1:1:0). Alexander staining and in vitro germination assay were carried out to examine the pollen development in different mutants. The ape1l +/2 zdp 2/2 mutant showed defects in pollen development and germination (S9 Fig.), suggesting that the ape1l +/2 zdp 2/2 mutation not only leads to maternal lethality but also gives rise to paternal defects.
The double mutations of APE1L and ZDP cause DNA hypermethylation and down-regulation of imprinted genes in the endosperm Maternal lethality phenotypes can be caused by aberrant expression of maternally imprinted genes and defects in the central cell or the endosperm [12,26,30]. FWA and FIS2 are two wellstudied maternally imprinted genes, and their maternal expression in the endosperm relies on active DNA demethylation initiated by DME [12,23]. We investigated whether the methylation of the FWA and FIS2 promoters in endosperm tissues is affected by APE1L and ZDP double mutations (Fig. 6A). The ape1l +/2 zdp 2/2 plants were backcrossed to zdp 2/2 plants three times to minimize the Ws background. To examine the methylation levels of DME target genes in our mutants, we employed the method of Buzas et al. [36] where the DNA methylation specific restriction enzyme McrBC is used to digest DNA before doing q-PCR in seeds at 3 days post manual pollination. We found that after digestion with McrBC, the amount of DNA recovered from FWA and FIS2 promoter regions (where is methylated in wild type leaf) was reduced in both dme and ape1l 2/2 zdp 2/2 mutants compared with wild type, but there was no difference in the unmethylated FWA gene body region (Fig. 6A). These results indicate that the ape1l 2/2 zdp 2/2 endosperm has hypermethylation in FWA and FIS2 promoter regions.
In order to measure the mRNA levels of FWA and MEA in Col and ape1l 2/2 zdp 2/2 endosperms, we carried out real-time PCR and found that the expression levels of FWA and MEA but not the DME and FIE mRNAs are down-regulated in the ape1l 2/2 zdp 2/2 mutant endosperm (Fig. 6B). To confirm and further analyze the FWA expression change, we introduced a pFWA::DFWA-GFP reporter into the ape1l +/2 zdp 2/2 and ape1l 2/2 zdp +/2 mutants by crossing the mutants with a transgenic line expressing the reporter [23]. Both ape1l +/2 zdp 2/2 and ape1l 2/2 zdp +/2 plants produce about 50% seeds defective in pFWA::DFWA-GFP expression (Fig. 6C-6D and S6 Table). To our surprise, ape1l 2/2 zdp +/2 mutant also produced 50% GFP-off seeds even though it is not maternally lethal and it produces about 75% viable seeds (Table 1). It turns out that hypermethylation of pFWA::DFWA-GFP promoter and silencing of FWA-GFP can occur in mutants which do not show maternal lethality. In addition, it seems that GFP-off seeds can be viable, so 75% viable seeds may be comprised of 50% GFP-on seeds and 25% GFP-off seeds. Taken together, our data suggest that DNA hypermethylation and down-regulation of imprinted genes occur and may be the cause of defects in the ape1l 2/2 zdp 2/2 endosperm.

Discussion
Active DNA demethylation in plants is initiated by the ROS1 subfamily of 5-meC DNA glycosylases/lyases and presumably completed through a base excision repair pathway [2,37]. Previous work has reported that the 39-phosphatase ZDP and the scaffold DNA repair protein XRCC1 also function in active DNA demethylation in Arabidopsis [17,38]. AP endonucleases are known to catalyze post-excision events during base excision repair. Our study here demonstrates that APE1L, one of the Arabidopsis AP endonucleases, functions in active DNA demethylation by processing b-elimination products of the bifunctional 5-meC DNA glycosylases/lyases and generating a 39-OH group. APE1L-mediated reaction comprises a new branch of the DNA demethylation pathway downstream of ROS1, DME, DML2 and DML3 (Fig. 7). Our biochemical data show that APE1L has an additional, weak 39-phosphatase activity, and thus may also function in the other branch, perhaps redundantly with ZDP, to process b, d-elimination products. Interestingly, it has been recently reported that the wheat homolog of APE1L also possesses 39-phosphatase and 39-phosphodiesterase activities [39]. Our results suggest that APE1L not only functions downstream of ROS1, DML2 and DML3 in vegetative tissues to prevent DNA hypermethylation but also functions together with ZDP downstream of DME to control DNA demethylation and gene imprinting in the central cell and endosperm and is thus important for seed development.
Active DNA demethylation in mammals can be initiated through the deamination of 5meC by AID to generate thymine, or the oxidation of 5meC to generate 5-hydroxymethylcytosine (5hmC), and further to 5-formylcytosine (5fC) and 5-carboxycytosine (5caC) by the TET family of DNA dioxygenases [2,[40][41][42]. 5fC and 5caC can be excised by the monofunctional DNA glycosylase TDG, whereas thymine can be removed by the monofunctional DNA glycosylase MBD4. Thus, a base excision repair pathway is required for completing the DNA cleavage and cytosine insertion steps during active DNA demethylation in mammals. Little is known about the DNA repair factors involved in active DNA demethylation in mammals, but it is likely that mammalian APE functions in active DNA demethylation downstream of the DNA glycosylases.
The ape1l-1 mutation leads to DNA hypermethylation in thousands of genomic regions, indicating that APE1L is required for DNA demethylation in these regions in Arabidopsis. Like mutations in 5-methylcytosine DNA glycosylases/lyases such as ROS1, mutations in DNA repair enzymes downstream of these enzymes are expected to preclude active DNA demethylation and cause hypermethylation. Coordinating the DNA glycosylase/lyase and repair activity would be predicted to prevent an otherwise fatal accumulation of strand breaks throughout the genome [17]. APE1L and ROS1 physically interact in vitro and co-localize in vivo, strongly suggesting that these proteins form a complex which coordinates their activities. One may ask why the DNA demethylation pathway includes both lyase activity of ROS1 and AP endonuclease activity of APE1L. In a recent study, it was reported that Wheat APE1L has weak endonuclease activity but robust 39-repair phosphodiesterase and 39-phosphatase activities [43]. Even though we detected the endonuclease activity of Arabidopsis APE1L in vitro, it is possible that, like Wheat APE1L, Arabidopsis APE1L is weak in cleaving DNA backbone at AP sites when involved in DNA demethylation. In this case, the lysase activity of ROS1 is required for generating the DNA gap.
The Arabidopsis ARP endonuclease also processes the 39-PUA generated by ROS1 in vitro, although this activity is much weaker than APE1L. Our whole genome bisulfite sequencing data identified only a small number of DMRs in the arp-1 mutant. Therefore, ARP is unlikely to play a major role in DNA demethylation, at least under normal growth conditions. Interestingly, we detected many genomic regions that are hypomethylated in the ape1l-1 mutant. APE1L is a multifunctional enzyme; its APE and 39-phosphatase activities may contribute to other DNA repair pathways in addition to active DNA demethylation, Thus, APE1L dysfunction may affect many DNA-related processes that directly or indirectly cause DNA hypomethylation. Compared to the ros1-4 and rdd mutations, ape1l and arp mutations induce higher percentages of hypermethylation in genic regions, whereas the zdp mutation induces a higher percentage of hypermethylation in TEs. The mechanisms underlying this genomic specificity are unclear, but it is possible that APE1L and ARP function redundantly in the demethylation of TEs, such that mutating either one individually does not cause hypermethylation. Unlike ZDP, which processes 39-phosphate blocking ends and promotes the release of ROS1 from its products, APE1L converts both 39-phosphate and 39-PUA to 39-OH, but does not increase the turnover of ROS1. Although both ZDP and APE1L interact with ROS1 in vitro and co-localize with ROS1 in vivo, ZDP and APE1L do not show extensive co-localization. It is possible that ZDP and APE1L exist mostly in two different protein complexes (Fig. 7). ZDP dysfunction caused DNA hypermethylation and transcriptional silencing of a luciferase reporter driven by the RD29A promoter, although the mutant phenotype was less severe than ros1 mutants [17]. We hypothesized that at some DNA demethylation target regions, such as the RD29A promoter, the DNA glycosylases/lyases may use both band b, d-elimination activities and thus require both APE and ZDP to process the intermediates and prevent transcriptional silencing. However, we found that the ape1l-1 and arp mutations did not affect expression of the reporter gene (S7B Fig.). It is possible that APE1L may function redundantly with ARP and/or ZDP in demethylation of the RD29A promoter. zdp mutant showed sensitivity to MMS but ape1l and arp mutants are not sensitive to MMS probably because they carry out different reactions. In addition, APE1L, APE2 and ARP may play redundant roles in repairing MMS-induced DNA damage, such that the single mutation or double mutations are not sufficient to induce sensitivity to MMS.
The choice between the APE branch and the ZDP branch of the active DNA demethylation pathway depends on the elimination mechanism used by the DNA glycosylases/lyase enzymes. It is unclear when and where a DNA glycosylases/lyase employs belimination, b, d-elimination, or both. Knowing which genomic regions depend on APE1L and which depend on ZDP for demethylation would be helpful. However, because the zdp-1 and ape1l-1 mutants are in different ecotypes, it is not ideal to compare the genomic regions targeted by the two different branches of the demethylation pathway.
The double mutations of APE1L and APE2 are embryonic lethal, but not paternally or maternally lethal based on our results and the segregation ratio of selfed ape1l +/2 ape2 2/2 reported previously [31]. It is possible that the lethal phenotype caused by APE1L and APE2 double mutations reflect deficiencies of DNA repair. Interestingly, we found that the ape1l +/2 zdp 2/2 mutant shows a maternal lethality phenotype, which has been shown to occur in other mutants that are defective in DNA demethylation, such as the dme and ssrp1 mutants [26,30]. Unexpectedly, only the ape1l +/2 zdp 2/2 mutant shows maternal lethality but the ape1l 2/2 zdp +/2 mutant is not maternally lethal. As a result of maternal lethality, about 50% of seeds abort in dme +/2 and ape1l +/2 zdp 2/2 mutants. In contrast, about 25% seeds abort in the ape1l 2/2 zdp +/2 mutant. All of the aborting seeds display embryos arrested at early growth stages presumably because an abnormal endosperm cannot support normal growth of the embryo. The morphology of aborted seeds in the ape1l +/2 zdp 2/2 and ape1l 2/2 zdp +/2 mutants is almost the same as that in the ape1l +/2 ape2 2/2 mutant, which is not maternally lethal and gives about 25% aborted seeds [31]. It is likely that the ape1l +/2 ape2 2/2 mutant is also defective Fig. 7. A working model of the active DNA demethylation pathway in Arabidopsis. ROS1 (or DME, DML2 or DML3) is a bifunctional DNA glycosylase/lyase that removes the 5-methylcytosine base and then cleaves the DNA backbone at the abasic site using bor b, d -elimination, resulting in a gap with PUA or phosphate at the 39 terminus which needs to be removed by APE1L or ZDP. Then the gap will be filled with an unmethylated cytosine nucleotide by an as yet unknown DNA polymerase and the processed strand will be sealed by a DNA ligase. The question mark represents unknown DNA polymerase and ligase enzymes. doi:10.1371/journal.pgen.1004905.g007 in DNA demethylation. Alternatively, this type of morphology (arrested embryo and aberrant endosperm) may reflect deficiencies of base excision repair.
It is likely that APE1L and ZDP function downstream of DME in the active DNA demethylation pathway that controls seed development. However, the aborting seeds in ape1l +/2 zdp 2/2 mutants have varied sizes of endosperm nuclei but the aborting seeds in dme +/2 mutants have endosperm nuclei of uniform sizes. This phenotypic difference may arise because APE1L and ZDP have multiple functions in DNA demethylation and repair, whereas DME only participates in DNA demethylation. As in dme +/2 mutants, the seed abortion phenotype in ape1l +/2 zdp 2/2 mutants is associated with the hypermethylation of the FWA promoter and the MEA ISR, and reduced FWA and MEA expression. Similar to the ape1l +/2 zdp 2/2 mutant, the ape1l 2/2 zdp +/2 mutant also produces about 50% GFP-off seeds, suggesting that these two types of mutants are similarly defective in DNA demethylation of imprinted genes. The phenotype of the ape1l 2/2 zdp +/2 mutant in pFWA-GFP silencing (50% GFP-off) and seeds viability (25% aborted and 75% viable) resembles that of the recently discovered atdre2 mutant [44]. Some other factors beyond DNA demethylation or some dosage effects must be differentially involved in different types of mutants, leading to maternally lethality in some mutants but not in others, even though they are all defective in the expression of imprinted genes. In summary, our results show that APE1L and ZDP are important regulators of gene imprinting in plants, and suggest that DMEinitiated active DNA demethylation in the central cell and endosperm employs both APE-and ZDP-dependent mechanisms.

Protein expression and purification
Full-length APE1L and APE2 cDNAs were subloned into pMAL-c2X (New England Biolabs) to generate MBP-APE1L and MBP-APE2 fusion proteins. The full-length ARP cDNA was subcloned into pET28a (Novagen) to generate a His-ARP fusion protein. Expression was induced in Escherichia coli BL21 (DE3) dcm 2 Codon Plus cells (Stratagene). MBP-APE1L and MBP-APE2 were purified by amylose affinity chromatography (New England Biolabs) and His-ARP was purified by affinity chromatography on a Ni 2+ -nitrilotriacetic acid column (Amersham Biosciences). His-ROS1 and MBP-ROS1 were expressed and purified as previously described [15,34].

Site directed mutagenesis
Site directed mutagenesis of APE1L was performed using the Quick-Change II XL kit (Stratagene) according to the manufacturer's instructions. The N224D mutation was introduced into pMal-APE1L by using the oligonucleotides APE1LN212D_F4 and APE1LN212D_R4 (see S1 Table). The mutant sequence was confirmed by DNA sequencing, and the construct was used to transform E. coli strain BL21 (DE3) dcm 2 Codon Plus cells (Stratagene). Mutant protein was expressed and purified as described above for APE1L.

DNA substrates
Oligonucleotides used to prepare DNA substrates (see S2 Table) were synthesized by Integrated DNA Technologies [45] and purified by PAGE before use. Double-stranded DNA substrates were prepared by mixing a 5 mM solution of a 59-fluoresceinlabeled or 59-Alexa Fluor-labeled oligonucleotide (upper strand) with a 10 mM solution of an unlabeled oligomer (lower strand). For preparation of 1-nt gapped DNA, a 5 mM solution of the corresponding 59-labelled oligonucleotide was mixed with 10 mM solutions of unlabelled 59-phosphorylated oligonucleotides P30_51 and CGR. Annealing reactions were performed at 95uC for 5 min, followed by slow cooling to room temperature.

Enzyme assays
To detect 5-meC DNA glycosylase/lyase activity, purified His-ROS1 (35 nM) was incubated at 30uC for 4 h with a Alexa Fluorlabeled DNA duplex (20 nM), containing a single 5-meC, in a reaction mixture containing 50 mM Tris-HCl pH 8.0, 1 mM DTT, 0.1 mg/ml BSA. In reactions containing APE1L, the mixture also included 200 mM NaCl and 1 mM MgCl 2 . Reactions were stopped by adding 20 mM EDTA, 0.6% sodium dodecyl sulfate, and 0.5 mg/ml proteinase K, and the mixtures were incubated at 37uC for 30 min. DNA was extracted with phenol/chloroform/isoamyl alcohol (25:24:1) and ethanol precipitated at 220uC in the presence of 0.3 mM NaCl and 16 mg/ml glycogen. When the ROS1 reaction products were used as purified substrates for AP endonucleases (see below), samples were resuspended in 5 ml of distilled water. Otherwise, they were resuspended in 10 ml of 90% formamide, heated at 95uC for 5 min, and separated in a 12% denaturing polyacrylamide gel containing 7 M urea. Alexa Fluor-labeled DNA was visualized using the blue fluorescence mode of the FLA-5100 imager and analyzed using Multigauge software (Fujifilm).
The AP endonuclease activity was detected using a DNA substrate containing a synthetic AP site (tetrahydrofuran, THF) opposite G. The 39-phosphatase activity was assayed on a 1-nt gapped substrate containing 39-phosphate and 59-phosphate ends. The 39-phosphodiesterase activity was tested on purified ROS1 products, which contain a mixture of fragments with 39-PUA and 39-phosphate termini. In all assays, purified AP endonucleases were incubated with DNA substrates (20 or 40 nM) at 30uC for the indicated times in a reaction mixture containing 50 mM Tris-HCl pH 8.0, 200 mM NaCl, 1 mM DTT, 0.1 mg/ml BSA and 1 mM MgCl 2 . Reactions were stopped and products analyzed as indicated above.

Pull down assays
Purified MBP alone or MBP-APE1L (200 pmol) in 100 ml of Column Buffer (20 mM Tris, pH 7.4, 1 mM EDTA, 1 mM DTT, 0.5% Triton X-100) was added to 100 ml of amylose resin (New England Biolabs) and incubated for 1 h at 4uC. The resin was washed twice with 600 ml of Binding Buffer (10 mM Tris, pH 8.0, 1 mM DTT, 0.01 mg/ml BSA). Purified His-ROS1 (15 pmol) was incubated at 25uC for 1 h with either MBP or MBP-APE1L bound to resin. The resin was washed twice with Binding Buffer. Bound proteins were analyzed by Western blot using antibodies against His 6 tag (Novagen).

Electrophoretic Mobility Shift Assay (EMSA)
EMSAs were performed using an Alexa Fluor-labeled duplex containing a gap flanked by 39-phosphate and 59-phosphate termini prepared as described above. The labeled duplex substrate (10 nM) was incubated with MBP-APE1L and/or His-ROS1 at the indicated concentrations in DNA-binding reaction mixtures (10 ml) containing 10 mM Tris HCl, pH 8.0, 1 mM DTT, 10 mg/ ml BSA. After 15 min incubation at 25uC, reactions were immediately loaded onto 0.2% agarose gels in 16 Tris acetate/ EDTA. Electrophoresis was carried out in 16 Tris acetate/EDTA for 40 min at 80 V at room temperature. Alexa Fluor-labeled DNA was visualized in a FLA-5100 imager and analyzed using MultiGauge software (Fujifilm).

Firefly luciferase complementation imaging assay
To investigate the interaction between APE1L and ROS1, two constructs was generated: APE1L-Cluc and ROS1-Nluc. The BamHI and SalI sites were used for cloning APE1L genomic DNA into pCAMBIA1300-CLUC vector. ROS1 was introduced to NLUC by In-Fusion HD Cloning Kit (Clontech). For protein interaction analysis, two combinatory constructs were transformed simultaneously into Nicotiana benthamiana leaves. To prevent the silencing of those genes, a virus p19 protein gene containing construct was transformed at the same time. After 3 d, 1 mM luciferin was sprayed onto the lower epidermis and kept in the dark for 5 min, then a CCD camera (1300B; Roper) was used to capture the fluorescence signal at 21uC.

Plant materials
Two T-DNA insertion mutants of the APE1L gene (At3g48425), INRA Flag240B06 and Salk_024194C, were used and they were referred to as ape1l-1 and ape1l-2 respectively. T-DNA insertions are present in the fifth exon and fourth intron of ARP in arp-1 (SALK_021478) and arp-2 (SAIL_866_H10) respectively. For all plants, seeds were sown on 1/2 MS plates containing 2% sucrose and 0.7% agar, stratified for 48 hours at 4uC and grown under long day conditions at 22uC. They were collected at 14 days or transplanted to soil.

Immunofluorescence
Immunofluorescence was performed in 2-to 3-week-old leaves as described by Pontes et al., [46]. Nuclei preparations were incubated overnight at room temperature with rabbit anti-APE1L (anti-APE1L antibodies were generated by injecting rabbits with a recombinant full length APE1L protein that was purified by affinity chromatography), anti-ZDP [17] and mouse anti-Flag (F3165, Sigma). Primary antibodies were visualized using mouse Alexa 488-conjugated and rabbit Alexa-594 secondary antibody at 1:200 dilution (Molecular Probes) for 2 h at 37uC. DNA was counterstained using DAPI in Prolong Gold (Invitrogen). Nuclei were examined with a Nikon Eclipse E800i epifluorescence microscope equipped with a Photometrics Coolsnap ES.

RNA purification and real-time PCR
Total RNA was extracted from 2-week-old seedlings using the RNeasy Plant Mini Kit (QIAGEN). 2-mg RNA was used for the first-strand cDNA synthesis with the Super script III First-Strand Synthesis System (Invitrogen) for RT-PCR following the manufacturer's instructions. The cDNA synthesis reaction was then diluted five times, and 1 ml was used as template in a 20-ml PCR reaction with iQ SYBR Green Supermix (Bio-Rad). All reactions were carried out on the iQ5 Multicolour Real-Time PCR Detection System (Bio-Rad). The comparative threshold cycle (Ct) method was used for determining relative transcript levels (Bulletin 5279, Real-Time PCR Applications Guide, Bio-Rad), with TUB8 as an internal control.
Whole genome bisulfite sequencing and data analysis DNA was extracted from 2 g of 12-day-old seedlings grown in a growth chamber and sent to BGI (Shenzhen, China) for bisulfite treatment, library preparation, and sequencing.

Microscopy
Images of seed phenotypes were captured using an Olympus SZX7 microscope equipped with a Canon Powershot A640 camera. For cleared whole-mount observation, immature seeds, that are 8 days after pollination, were cleared using chloral hydrate, glycerol, and water (8 g: 1 ml: 2 ml) and photographed using a Leica DM6000 B differential interference contrast microscope equipped with a Leica DFC 425 camera. Fluorescence was detected with an Olympus BX53 fluorescence microscope equipped with an Olympus DP80 digital camera.

McrBC assay
The McrBC assay was performed according to Buzas et al [36]. Briefly, wild type and the apel1 +/2 zdp 2/2 mutant were pollinated with Ler pollen. 3 days after pollination, pools of GFP-on and GFPoff seeds were selected under a dissecting fluorescence microscope and more than 300 seeds were used for DNA extraction. Genomic DNA concentration was measured by Nanodrop. Approximately 1 mg of DNA was digested with 1 mL of McrBC overnight at 37uC. After digestion, DNA methylation levels at the specific loci were determined by real-time PCR using absolute quantification against a 1:1 mixture of genomic DNA extracted from Col-0 and Ler leaves. Primers are listed in S1 Table. RNA purification and real-time PCR using endosperm tissues Female ape1l +/2 zdp 2/2 plants (Col-0) were crossed with male wild type plants. The endosperm plus seed coat fraction was collected for RNA purification using the Trizol method. DNAase treatment and LiCl precipitation were applied to remove DNA and polysaccharide contaminations, respectively. RNA was reverse transcribed into cDNA by the SuperScript III First-Strand Synthesis System (Invitrogen) with an oligo dT primer. Real-time PCR analysis was performed using SYBR Premix Ex Taq (TaKaRa) and CFX96 real-time system (Bio-Rad). ACT11 was used as the internal control.

Accession numbers
We used whole-genome bisulfite sequencing to analyze the methylomes of Ws, ape1l-1, arp-1 and zdp-1 mutant plants. The data set was deposited at NCBI (GSE52983).

S1 Fig
The DNA glycosylase/lyase activity of ROS1 is not increased in the presence of APE1L. Purified His-ROS1 (2 or 10 nM) was incubated in a reaction buffer lacking Mg 2+ with a DNA substrate containing a 5-meC (20 nM) either in the absence or presence of purified MBP-APE1L (2, 10, 20 or 100 nM). Reactions were stopped after 16 hours and products were separated in a 12% denaturing polyacrylamide gel and detected by fluorescence scanning. S8 Fig Effects of ape1l and zdp double mutations on seed development. (A) A wild type seed, and an aborting seed showing endosperm over-proliferation (B) Cleared seed samples showing arrested embryo or invisible embryo in ape1l +/2 zdp 2/2 , ape1l 2/2 zdp +/2 and dme +/2 mutants as compared to seeds in Col. Images were captured at 8 DAP. (C) Images of cleared seed samples showing unequal sizes of endosperm nuclei in ape1l +/2 zdp 2/2 and ape1l 2/2 zdp +/2 mutants but not in Col and dme +/2 . Arrows indicate abnormally large endosperm nuclei. Images were captured at 8 DAP under a differential interference contrast microscope. (PDF)   Table The GFP phenotype in self crosses of ape1l 2/2 zdp +/2 , ape1l +/2 zdp 2/2 and pFWA-GFP. (DOCX)