PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of bloodadvBlood Advances
 
Blood Adv. 2017 September 12; 1(20): 1699–1711.
Published online 2017 September 6. doi:  10.1182/bloodadvances.2017005710
PMCID: PMC5728331

RUNX1 regulates site specificity of DNA demethylation by recruitment of DNA demethylation machineries in hematopoietic cells

Key Points

  • Ectopic expression of RUNX1 induces binding site–directed DNA demethylation, in which hematopoietic gene promoters are included.
  • RUNX1 binding sites are enriched in demethylated regions during hematopoietic development.

Abstract

RUNX1 is an essential master transcription factor in hematopoietic development and plays important roles in immune functions. Although the gene regulatory mechanism of RUNX1 has been characterized extensively, the epigenetic role of RUNX1 remains unclear. Here, we demonstrate that RUNX1 contributes DNA demethylation in a binding site–directed manner in human hematopoietic cells. Overexpression analysis of RUNX1 showed the RUNX1-binding site–directed DNA demethylation. The RUNX1-mediated DNA demethylation was also observed in DNA replication–arrested cells, suggesting an involvement of active demethylation mechanism. Coimmunoprecipitation in hematopoietic cells showed physical interactions between RUNX1 and DNA demethylation machinery enzymes TET2, TET3, TDG, and GADD45. Further chromatin immunoprecipitation sequencing revealed colocalization of RUNX1 and TET2 in the same genomic regions, indicating recruitment of DNA demethylation machinery by RUNX1. Finally, methylome analysis revealed significant overrepresentation of RUNX1-binding sites at demethylated regions during hematopoietic development. Collectively, the present data provide evidence that RUNX1 contributes site specificity of DNA demethylation by recruitment of TET and other demethylation-related enzymes to its binding sites in hematopoietic cells.

Visual Abstract

figure advances005710absf1

Introduction

RUNX1 is an essential master transcription factor (TF) during hematopoietic development. All blood cell types differentiate from pluripotent stem cells, such as embryonic stem cells and induced pluripotent stem cells (iPSCs), through hematopoietic stem cells in vitro and in vivo.1 RUNX1-deficient mouse embryonic stem cells do not commit to the hematopoietic lineage.2 In addition, RUNX1 mutations and chromosomal truncations are frequently observed in acute myeloid leukemia.3,4 Thus, RUNX1 plays a pivotal role in normal hematopoietic development.

RUNX1 was identified as a transcriptional activator that regulates expression of downstream hematopoietic genes, such as interleukin-3, granulocyte-macrophage colony-stimulating factor, and PU.1 (SPI1).5-7 In addition, RUNX1 interacts with either activation cofactors, such as p300 and CBP (Creb Binding Protein), histone acetyltransferases,8,9 or repression cofactors such as SUV39H1, a histone methyltransferase.10 RUNX1 has also been shown to enhance EBF-derived DNA demethylation at the mb-1 promoter region.11 Accordingly, RUNX1 appears to contribute to not only direct gene activation but also epigenetic regulation by interacting with functional epigenetic enzymes.

DNA methylation is a crucial epigenetic modification because methylation of cytosine residues in cytosine guanine dinucleotide (CpG) dinucleotides of gene promoters represses transcription by prevention of TF binding, and therefore, DNA demethylation is fundamental for cell type–specific gene expression.12 DNA methylation is precisely regulated during cell differentiations, and abnormal DNA methylation is frequently observed in cancers. Hence, DNA methylation inhibitors such as 5-aza-2′-deoxycytidine are expected to be anticancer drugs.

De novo methylation of unmethylated cytosines is performed by DNA methyltransferase (DNMT) 3a/b, and the methylated status is inherited to daughter cells by DNMT1.13-15 In contrast, DNA demethylation occurs through passive16 and active17 mechanisms. Specifically, passive demethylation occurs after dilution of hemimethylated cytosines during DNA replication in the absence of DNMT1 activity, whereas multiple enzymes perform active DNA demethylation in a DNA replication-independent manner. Among these mechanisms, ten-eleven translocation methylcytosine dioxygenase (TET) oxidizes methylated cytosine nucleotides to hydroxyethyl cytosine, formyl cytosine, and then carboxyl cytosine.18-20 Subsequently, formyl and carboxyl cytosines are recognized as abnormal bases that are replaced by following base excision repair (BER) by thymine–DNA glycosylase (TDG), GADD45A, and various other proteins.21

During hematopoietic development, spatiotemporal changes in DNA methylation occur at the genome-wide level.22 Recently, several TFs, which play important roles in cellular differentiations, were reported to induce DNA demethylation by interacting with TET proteins.23-26 However, it has not been elucidated whether RUNX1 also contributes to regulating DNA methylation. Here, we demonstrate that RUNX1 contributes site-directed DNA demethylation by recruitment of DNA demethylation enzymes. Our findings provide evidence of a novel epigenetic function of RUNX1 and its potential contribution to hematopoietic development.

Methods

Cells and cell culture

Experiments were performed with human iPSCs (RIKEN BRC, HPS0002), CD34+ hematopoietic cells (Lonza, Basel, Switzerland; 2M-101c), CD14+ monocytes (MON; Lonza, 2W-400C), CD3+ T cells (HemaCare BioResearch, Van Nuys, CA; PB03C-2), CD19+ B cells (HemaCare BioResearch; PB19C-2), and peripheral blood mononuclear cells (Lonza, CC-2702). 293T cells were acquired from the BRC and cultured in Dulbecco’s modified Eagle medium (Wako, Osaka, Japan) supplemented with 10% fetal bovine serum, 100 U/mL penicillin, and 100 µg/mL streptomycin. Jurkat cells were purchased from the BRC and cultured in RPMI-1640 (Wako) supplemented with 10% fetal bovine serum, 100 U/mL penicillin, and 100 µg/mL streptomycin.

Methylation array analysis

Genomic DNA was isolated using a NucleoSpin Tissue Kit (Macherey-Nagel, Düren, Germany) and examined using an Infinium Humanmethylation450 BeadChip (Illumina, San Diego, CA) according to the manufacturer’s instructions. Data were quantile-normalized, and M values were calculated using the lumi Bioconductor package. CpGs that were differentially methylated between samples (DMCs) were selected using a cutoff M-value difference (ΔM) of >2.

Overrepresentation analyses of TF-binding motif

The TFBMs were searched in ±5-kbp regions from DMCs using the Bioconductor Biostrings package based on the hg19 genomic sequence. Motifs were counted in a 100-bp window sliding by 50-bp steps in ±200-bp regions from DMCs. Comparisons with the same numbers of randomly selected probes were performed using an exact test for the Poisson distribution model (P < .00001). Enrichment scores were then computed for significantly differing motifs using the following formula:

equation image

Here, i and j represent identified motifs in DMC regions and randomly selected controls, respectively, and xt and xc indicate motif positions in DMC regions and randomly selected controls, respectively. The smoothing parameter h was set at 50, and K was calculated as a Gaussian kernel function, as follows:

equation image

The significance of the enrichment is further assessed by positional profile: whether the peak of enrichment score is greater than barrier of Q3 + interquartile range × 3 in ±5-kbp regions.

ChIP-sequencing

Chromatin immunoprecipitation (ChIP) assays were performed as described previously27 using antibodies shown in supplemental Table 5. ChIPed DNA was subjected to library construction using a Library Preparation Kit (KAPA Biosystems, Wilmington, MA) and sequenced using the HiSequation 2500 system (Illumina). FASTQ data were mapped to human hg19 using Bowtie2. Peak calling and visualization were performed using MACS228 (cutoff P = .0001) and Ngsplot.

Bisulfite sequencing

Bisulfite conversion was performed using an EZ DNA Methylation-Gold Kit (Zymo Research, Irvine, CA), and the target genomic regions were amplified using Epi-Taq HS (Takara Bio, Shiga, Japan) with the primers shown in supplemental Table 6. Polymerase chain reaction (PCR) products were cloned into pTA2 plasmids (Toyobo, Osaka, Japan), and 24 clones per target were sequenced. Visualization and statistical analysis were performed using QUMA.

Lentivirus preparation and transduction

The lentiviral vector CSII-EF-RfA was kindly provided by H. Miyoshi (RIKEN BRC), and puromycin-resistance gene fragment was inserted into the CSII-EF-RfA (CSII-EF-RfA-IRES2-Puro). Entry clone of RUNX1 was synthesized by Codon Devices and of TET2 and TET3 were acquired from the ORFeome collection (HsCD00377074) and GeneCopoeia (GC-H2292), respectively. The RUNX1 gene was subcloned into CSII-EF-RfA-Bsd,29 and TET2 and TET3 genes were subcloned into CSII-EF-RfA-Puro. Lentiviral vector preparation and transduction were then performed as described previously.29 Following transfection, antibiotic selection was performed using 2 μg/mL puromycin or 20 μg/mL blasticidin for 1 week.

qRT-PCR

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) analyses were performed as described previously29 with the primers shown in supplemental Table 6. Gene expression was evaluated by ΔΔCt method compared with mock control sample.

Mitomycin C treatment

Confluent 293T cells were treated with mitomycin C at a final concentration of 5 μg/mL for 3 hours at 37°C, washed in phosphate-buffered saline twice, and then transduced with the RUNX1 lentiviral vector and cultured for 7 days.

3-Aminobenzamide treatment

Cells were plated at 1 × 104cells per well to 6-well plates. One day after the cell plating, the medium was replaced with fresh medium supplemented with 10 mM 3-aminobenzamide and RUNX1 lentivirus vector or mock control. After the medium replacement, the cells were cultured for 7 days.

qMSP

DNA was extracted using the NucleoSpin Tissue Kit, and the DNA methylation status was analyzed by quantitative methylation-specific PCR (qMSP) as described previously,30 with the primers shown in supplemental Table 6.

Co-IP and western blot

Cells were lysed in lysis buffer (50 mM Tris-HCl [pH 7.4], 1% NP-40, 0.25% Na-deoxycholate, 150 mM NaCl, 1 mM EDTA, Protease Inhibitor Cocktail). Antibodies were incubated with protein G Dynabeads (Thermo Fisher Scientific, Waltham, MA) for 1 hour, incubated with cell lysates at 4°C overnight, and washed twice in Tris-buffered saline washing buffer (0.05% IGEPA CA-630). Bead-bound proteins were eluted in elution buffer. Western blot was performed as described previously31 using the primary antibodies listed in supplemental Table 5 and a VeriBlot (Abcam, Cambridge, UK) IP secondary antibody.

HT immunoprecipitation analyses

HaloTag-fused expression vectors pFC14K-RUNX1, pFC14K-TDG, pFN21A-TET2, pFN21A-TET3, and pFN21A-GADD45A, expression vectors pCDNA3.1/V5-RUNX1, pCDNA3.1/V5-TDG, and pCDNA3.1/V5-GADD45A, and the V5-tag fused expression vector pCDNA3.2/nV5-TET2 and pCDNA3.2/nV5-TET3 were transfected into 293T cells using Fugne HD (Promega, Madison, WI) according to the manufacturer’s instructions. Cells were collected at 48 hours after transfection and high-throughput (HT)-based coimmunoprecipitation (co-IP) assays were performed using the HaloTag Mammalian Pull-Down System (Promega) according to the manufacturer’s instructions, followed by western blot.

Construction of deletion mutants

Entry clones of RUNX1 deletion mutants were constructed from a full-length RUNX1 by inverse PCR with the primers shown in supplemental Table 6 and were then subcloned into pCDNA3.2/nV5-DEST (Thermo Fisher Scientific).

RUNX1 knockdown

pLKO.1 shRUNX1 puro (#45816) and empty pLKO.1 puro (#8453) plasmids were purchased from Addgene. Lentiviral vector preparation and transduction were performed as described above. Knockdown efficiency was confirmed by qRT-PCR (supplemental Table 6).

Results

Ectopic expression of RUNX1 induces genome-wide binding site-directed DNA demethylation

To investigate whether RUNX1 is involved in DNA methylation regulation, we performed overexpression analyses in 293T cells that do not express RUNX1. Following ectopic expression of RUNX1 or the mock control by lentiviral vectors in 293T cells (293T-RUNX1 and 293T-mock, respectively), single-base resolution DNA methylation arrays were used to characterize global DNA methylation changes between 293T-RUNX1 and 293T-mock cells (Figure 1A). Comparing M-values, a statistical metric for log-scale methylation levels32 of 293T-RUNX1 cells with that of 293T-mock cells, we identified 145 and 59 demethylated and methylated CpGs (Figure 1B; supplemental Table 1). Interestingly, the RUNX1 promoter region was demethylated by the RUNX1 overexpression, suggesting an autofeedback loop, which consists of transcription and epigenetic alterations.

Figure 1.
Induction of DNA demethylation by RUNX1 overexpression. (A) Confirmation of RUNX1 overexpression by western blotting. (B) Scatter plot showing DMCs caused by RUNX1 overexpression in 293T cells. x- and y-axes show M-values for 293T-mock and RUNX1-overexpressing ...

RUNX1 binds to a specific consensus motif with a core sequence motif [TGTGG(TTT/TCA)].33 Therefore, we hypothesized that, if RUNX1 is involved in regulation of DNA methylation in a binding site-directed manner, the RUNX1-binding motif will be overrepresented at differentially methylated regions in RUNX1-overexpressing cells. To test this hypothesis, we performed TF-binding motif (TFBM) overrepresentation analyses by using an algorithm that was developed in our laboratory (supplemental Figure 1A, see “Methods”) in RUNX1-overexpressing cells. The result revealed significant overrepresentation of the RUNX1-binding motif in the immediate vicinity of the demethylated regions (P = 1 × 10−30, exact test for Poisson distribution; Figure 1C) but not in that of the methylated regions (supplemental Figure 1B). This RUNX1-binding motif overrepresentation at demethylated regions by RUNX1 overexpression was validated in an independent analysis by using a different methylation array platform (Human MethylationEpic BeadChip [Illumina]) (P < 1 × 10−100, exact test for Poisson distribution) (supplemental Figure 1C-D).

To further confirm link between DNA demethylation and RUNX1 binding, we also performed a chromatin immunoprecipitation-sequencing (ChIP-seq) analysis of RUNX1 in RUNX1-overexpressing 293T cells. The ChIP-seq analysis identified 96 158 peaks. Of these, 391 and 265 peaks were harbored in ±5 kb from demethylated and randomly selected probes, respectively, showing significant overrepresentation of RUNX1 ChIP peaks around the demethylated CpGs (P = 1.93 × 10−13, exact test for Poisson distribution) (Figure 1D-E).

In subsequent bisulfite sequencing analyses of 2 demethylated regions, cg00041401 (PTPN22) and cg07236781 (RUNX3), in RUNX1-overexpressing cells, both demethylated CpGs and neighboring CpGs were demethylated (Figure 1F).

Overall, these results indicate that RUNX1 induces DNA demethylation at its binding regions.

RUNX1-mediated DNA demethylation occurs via a combination of DNA replication-dependent and -independent mechanisms

The 2 mechanisms of DNA demethylation (active and passive) are DNA replication-independent and -dependent, respectively. To determine which mechanism is responsible for RUNX1-mediated DNA demethylation, we examined RUNX1-mediated DNA demethylation in DNA replication-arrested 293T cells using qMSP. Proliferation of 293T cells was arrested by treatment with 5 μg/mL mitomycin C (Figure 2A), followed by RUNX1 overexpression and qMSP analyses of 3 RUNX1-mediated DNA demethylation target sites that were randomly selected from the above methylation array results (Figure 1B). qMSP revealed that, although the demethylation was weaker than in mitomycin C untreated 293T cells, 2 RUNX1 targets were significantly demethylated and 1 was slightly demethylated, albeit not statistically significant, by RUNX1 overexpression compared with mock control overexpression in DNA replication-arrested 293T cells (Figure 2B). Conversely, 2 randomly selected negative control regions were not demethylated (Figure 2B). Furthermore, treatment of an inhibitor of active DNA demethylation, 3-aminobenzamide,34 significantly reduced the RUNX1-mediated DNA demethylation (supplemental Figure 2). These results suggest that both active and passive mechanisms contribute to RUNX1-mediated DNA demethylation.

Figure 2.
RUNX1-mediated DNA demethylation in mitomycin C–treated 293T cells. (A) 293T cell growth after mitomycin C (MMC) treatment. 293T cells were treated with the indicated concentrations of MMC, and then cell densities were measured after 1, 3, and ...

TET enhances RUNX1-mediated DNA demethylation

Although TET proteins (TET1, TET2, and TET3) play key roles in active DNA demethylation,20 expressions of these genes are low to middle in HEK293 cells (supplemental Figure 3A-B). Therefore, we hypothesized that overexpression of TETs enhance the RUNX1-mediated DNA demethylation. DNA methylation array analyses were performed in TET2- or TET3- and RUNX1-co-overexpressing cells (293T-T2R, 293T-T3R, respectively) (Figure 3A-B), because TET2 and TET3 are expressed in hematopoietic lineage cells as well as RUNX1. Because TET2 or TET3 single overexpression induced DNA methylation alterations in the absence of RUNX1, which can be considered nonspecific effects (supplemental Figure 3C-D), we compared 293T-T2R and 293T-T3R with TET2- and TET3-overexpressing 293T cells (293T-T2 and 293T-T3, respectively) to clarify the RUNX1 overexpression effect. These comparative analyses showed that RUNX1 overexpression together with TET2 or TET3 expression led to more demethylated CpGs (565 and 357, respectively) compared with RUNX1 single overexpression (Figure 3C). Of these demethylated CpGs, 97 (66.9%, P = 2.2 × 10−16; χ2 test) and 83 (57.2%, P = 2.2 × 10−16; χ2 test) in 293T-T2R and 293T-T3R cells overlapped with those in RUNX1-overexpressing cells, respectively, and 251 (44.4% and 70.3%) demethylated CpGs in 293T-T2R and 293T-T3R cells overlapped each other. In the proximal regions of these demethylated CpGs, RUNX1-binding motifs were significantly overrepresented (TET2 vs TET2 + RUNX1; P = 1.23 × 10−161, TET3 vs TET3 + RUNX1; P < 1 × 10−500), ruling out the possibility of nonspecific demethylation (Figure 3D). Taken together, these data indicate that TET2 and TET3 enhance RUNX1-mediated DNA demethylation.

Figure 3.
TETs enhance RUNX1-mediated DNA demethylation. (A) Confirmation of RUNX1, TET2, and TET3 overexpression in mock vector–overexpressing (293T-mock), RUNX1-overexpressing (293T-RUNX1), TET2-overexpressing (293T-TET2), TET3-overexpressing (293T-TET3), ...

RUNX1 physically interacts with TET, TDG, and GADD45A

We next examined physical interactions between RUNX1 and TET2. We performed co-IP assays followed by western blot analyses in the Jurkat T-lymphocyte line. Jurkat expressed RUNX1 and TET2 (supplemental Figure 4A), and the RUNX1-binding motif was overrepresented in demethylated regions compared with iPSCs (P < 1 × 10−500, exact test for Poisson distribution) (supplemental Figure 4B). Therefore, we used Jurkat as a model hematopoietic cell line for molecular analyses. Co-IP followed by western blot showed that RUNX1 and TET2 proteins coprecipitated each other (Figure 4A), indicating an endogenous physical interaction between RUNX1 and TET2. HaloTag (HT)-based co-IP assays also confirmed the interaction between RUNX1 and TET2 in 293T cells (Figure 4B; supplemental Figure 4C). The HT-based co-IP assay further revealed interactions between RUNX1 and BER pathway enzymes TDG and GADD45A35 that are involved in active DNA demethylation (Figure 4C-D; supplemental Figure 4D-E). We also detected interaction between RUNX1 and TET3 by using the HT-based co-IP assay in 293T cells (supplemental Figure 4F-G). Cumulatively, these results suggest that RUNX1 forms a DNA demethylation machinery complex by interacting with enzymes involved in active DNA demethylation.

Figure 4.
Physical interactions between RUNX1 and DNA demethylation proteins. (A) Co-IP of endogenous RUNX1 or TET2 followed by western blotting in Jurkat. Input indicates 0.05%, 0.1%, and 1% (from left) of nonimmunoprecipitated cell lysate; IgG, control IP with ...

To identify the TET2-binding domain of RUNX1, we performed HT-based co-IP analyses using RUNX1 deletion mutants. Because the DNA-binding domain of RUNX1 (RUNT domain) is necessary for sequence-specific binding to DNA, we designed 4 deletion mutants with intact RUNT domains (Figure 4E). V5-tagged RUNX1 deletion mutants and HT-fused TET2 were overexpressed in 293T cells followed by HT co-IP analyses. Three of the 4 RUNX1 mutants clearly coprecipitated with TET2, whereas the Δ300-480 RUNX1 mutant showed a very weak interaction (Figure 4F), suggesting that RUNX1 binds to TET via a region from 300 to 480 aa. To confirm disruption of the RUNX1-mediated DNA demethylation activity in the Δ300-480 RUNX1 mutant, we performed qMSP analysis. The results showed a lack of demethylation activity in the Δ300-480 RUNX1 mutant (Figure 4G), supporting the interaction of RUNX1 with TET2 through the 300-480 aa region. Overall, these results indicate that the interaction with TET2 is necessary for RUNX1-mediated DNA demethylation.

Genome-wide colocalization of RUNX1 and TET2

To determine whether RUNX1 and TET2 colocalize in the same chromatin regions, we performed ChIP-seq analyses of RUNX1 and TET2 in Jurkat. A total of 25 881 625 and 9 466 629 reads were obtained from RUNX1 and TET2 ChIP-seqs, respectively, which were mapped onto the hg19 human genome (Figure 5A). We identified 50 334 and 126 184 peaks in RUNX1 and TET2 ChIP-seqs, respectively. Of these peaks, 40 982 of RUNX1 peaks (81.4%) and 38 345 of TET2 peaks (30.4%) overlapped each other, indicating that a substantial number of RUNX1-binding sites overlap with TET2-binding sites (Figure 5B). Genome-wide analysis also showed enrichment of TET2 reads around the RUNX1 ChIP peaks, which was significantly reduced by RUNX1 knockdown (P = .028, Student t test) (Figure 5C-E; supplemental Figure 5A). The enrichment of TET2 reads around RUNX1 ChIP peaks was validated in human peripheral blood mononuclear cells (supplemental Figure 5B). These results suggest that TET2 colocalizes with RUNX1 in a RUNX1-dependent manner. Furthermore, heatmap analysis revealed that both RUNX1 and TET2 ChIP-seq reads were highly enriched at transcription start sites (TSSs) (Figure 5F) and CpG islands (Figure 5G) as described previously.36,37 Collectively, these results suggest RUNX1:TET2 colocalization at gene regulatory regions such as promoters and CpG islands.

Figure 5.
ChIP-seq analyses of RUNX1 and TET2. (A) A ZENBU browser screenshot showing a typical relationship between RUNX1 and TET2 ChIP-seq data, and the methylation level. The DNA methylation track is the M-value measured by the methylation array. Green and purple ...

To investigate differences in the DNA methylation levels between regions with RUNX1 and with TET2, we superimposed ChIP-seq data and methylation array data. RUNX1 tended to strongly bind to hypomethylated regions and weakly bind to hypermethylated regions (Figure 5H top). Interestingly, regions with TET2 tended not to harbor hypomethylation (Figure 5H middle). Therefore, regions with both RUNX1 and TET2 biased toward hypermethylation with weak RUNX1 binding (Figure 5H bottom). Taken together, the results suggest that colocalization of RUNX1 and TET2 is affected by DNA methylation level of their binding region.

RUNX1-binding sites are significantly overrepresented in demethylated regions during hematopoietic development

To investigate whether the RUNX1-mediated DNA demethylation is involved in hematopoietic development, we performed TFBM overrepresentation analysis of differentially methylated regions during hematopoietic development. DNA methylation arrays were used to characterize global DNA methylation changes during hematopoietic differentiation of iPSCs, CD34+ HPCs, CD14+ MONs, CD19+ B cells, and CD3+ TCs (Figure 6A). To identify differentially methylated CpGs, we compared M-values between iPSCs and HPCs (iPSCs-HPCs) and between HPCs and 3 terminally differentiated hematopoietic cell types (HPCs-MONs, HPCs-BCs, and HPCs-TCs). These analyses identified 24 289, 1814, 2360, and 11 733 methylated (Figure 6B green dots) and 32 434, 12 410, 10 278, and 26 376 demethylated (Figure 6B purple dots) CpGs in iPSC-HPCs, HPCs-MONs, HPCs-BCs, and HPCs-TCs, respectively (supplemental Table 2). Importantly, out of demethylated CpGs by the RUNX1 overexpression, 82 (56.6%) overlapped with demethylated CpGs in any of the hematopoietic differentiations (P = 1.58 × 10−106, χ2 test) (supplemental Table 1, “Overlap” column), including CpGs harbored in the promoter (from 1.5 kbp upstream of the TSS to the first exon) of hematopoietic genes, such as PTPN22, RUNX1, RUNX3, IL17RE, LSP1, and BST2 (Tetherin).38-44 Further characterization of the 82 overlapped demethylated regions by using Roadmap epigenetics datasets45 revealed that H3K4me1, an enhancer mark, was highly enriched at the demethylated regions in hematopoietic cell types, suggesting that RUNX1-mediated DNA demethylation predominantly targets enhancer regions (supplemental Figure 6A).

Figure 6.
RUNX1-binding motif overrepresentation analysis of differentially methylated regions during hematopoietic development. (A) Differentiation hierarchy of hematopoietic development. (B) Scatter plots of M-values between induced pluripotent stem cells (iPSCs) ...

Subsequently, to assess TFBM overrepresentation in an unbiased way, we screened overrepresented TFBMs. Out of 113 human and mouse TFBMs, which are registered in the JASPAR CORE,46 29, 10, 7, and 12 were overrepresented at demethylated regions in iPSCs-HPCs, HPCs-MONs, HPCs-BCs, and HPCs-TCs, respectively, which include hematopoietic linage-specific TF responsible motifs, such as SPI1, RUNX1, SPIB, NR4A2, IRF1, IRF2, and AP-147-51 (supplemental Table 3). Notably, the RUNX1-binding motif was overrepresented at the demethylated regions of all differentiations (Figure 6C top). On the other hand, we identified 7 and 4 TFBMs at the methylated regions in iPSCs-HPCs, and HPCs-TCs, respectively, which did not include the RUNX1-binding motif (Figure 6C bottom; supplemental Table 4). The RUNX1-binding motif overrepresentation was positively correlated with RUNX1 expression, which is computed from the FANTOM5 database52 (R = 0.79) (supplemental Figure 6B). Furthermore, the RUNX1-binding motif overrepresentation in demethylated regions was reproduced in different donor samples of published methylation array data (GSE40790): between iPSCs (iGP91-02) and 2 lines of human peripheral blood CD34+ cells (CD34_1 and CD34_2), or human peripheral blood CD14+ cells (CD14) (supplemental Figure 6C). Thus, these results suggest that RUNX1 is 1 of the TFs that contribute to the DNA demethylation regulation in hematopoietic differentiation.

Because 1 TFBM sometimes corresponds to multiple TFs, we examined whether the RUNX1-binding motif overrepresentation at demethylated regions reflect actual RUNX1 binding. We compared the demethylated regions of iPSC-HPCs with RUNX1 ChIP-seq peaks of HPCs. We identified 8768 RUNX1-binding peaks in publically available RUNX1-ChIP-seq data of HPCs (GSE45144). Of these, 3130 peaks were located at ±5-kb regions from demethylated CpGs, which showed statistically significant overrepresentation compared with ±5-kb regions from randomly selected probes (P = 3.85 × 10−10, exact test for Poisson distribution), indicating the DNA demethylation links to RUNX1 binding (Figure 6D).

Collectively, these results suggest an involvement of RUNX1 in DNA demethylation at a subset of its binding site proximal regions at multiple hematopoietic development stages.

Discussion

In this study, we demonstrated that RUNX1 contributes site specificity of DNA demethylation by both active and passive mechanisms. In the active DNA methylation, we showed recruitment of DNA demethylation machinery by RUNX1. A possible mechanism of TF-mediated passive DNA demethylation is that DNMT1 recruitment is physically inhibited by TF binding,53,54 which may also occur in RUNX1-mediated DNA demethylation .In addition, hydroxymethyl-, formyl-, and carboxyl-cytosines, which are oxidized by RUNX1-mediated recruitment of TET, may block maintenance methylation during cell division, although it is still controversial.55-57

We revealed enhancement of RUNX1-mediated DNA demethylation by TET2 or TET3 co-overexpression. In the context of RUNX1 overexpression, RUNX1 proteins may be in excess compared with endogenous TET proteins because TET expression level is not high. Therefore, additional overexpression of TET may compensate for the imbalance of RUNX1 and TET proteins, resulting in enhancement of RUNX1-mediated DNA demethylation. Another possible mechanism of demethylation enhancement by TET overexpression may be increased accessibility of DNA. However, because we compared TET and RUNX1 co-overexpression with TET single overexpression to remove the effect of TET alone, this possibility can also be eliminated.

In zebrafish, Tet2 and Tet3 have overlapping functions in promotion of normal development.58 Furthermore, both TET2 and TET3 regulate GlcNAcylation by interacting with O-GlcNAc transferase.59 Our data indicate that both TET2 and TET3 enhance RUNX1-mediated DNA demethylation, and each effect is highly overlapped, supporting the functional overlap between TET2 and TET3.

Co-IP analyses showed interactions between RUNX1 and active DNA demethylation machinery enzymes. BER is a pivotal step in active DNA demethylation. Interestingly, RUNX1 regulates the expression of BER-related genes, including OGG1 and OLE.60,61 Thus, RUNX1-mediated DNA demethylation may be enhanced by a transcription-based auto-feedback loop.

ChIP-seq analyses revealed that TET2 was not associated with RUNX1 at hypomethylated genomic regions, although RUNX1 and TET2 were generally colocalized. This result suggests that TET2-RUNX1 complexes dissociate rapidly following to the completion of demethylation. A recent study showed that IDAX protein binds to both unmethylated CpG and TET2, and then downregulates TET2 protein expression, although the related mechanisms have not been described.62 This finding suggests the involvement of IDAX in the reduction of TET2 from demethylated regions.

The proposed mechanism of RUNX1-mediated demethylation is consistent with that of other reported TFs that induce DNA demethylation in a TET-dependent manner.23-26 Thus, TFs with a TET-dependent demethylation function may comprise a functional subclass. However, because a limited number of TFs have been identified in this subclass, it is unclear how major this subclass is. Therefore, whether this mechanism is a general function of TFs should be investigated further.

Our data showed that RUNX1-mediated DNA demethylation may be involved in hematopoietic development. Although it is widely established that RUNX1 is essential for hematopoietic commitment, and dysfunction of RUNX1 is linked to hematological diseases such as leukemia, the contribution of the demethylation function of RUNX1 to normal development and diseases is completely unknown. Therefore, physiological and pathological analyses of RUNX1-mediated DNA demethylation are of great interest.

DNA demethylation is an essential step for gene activation, changing “silenced state” to “poised state,”63 enabling additional TF binding.64 Thus, RUNX1-mediated DNA demethylation may unlock the silenced hematopoietic genes, and other factor(s) subsequently regulates the expression.

DNA methylation inhibitors such as 5-azacytidine are used to treat precancerous lesions in myelodysplastic syndrome, suggesting that aberrant DNA methylation is a major factor in oncogenesis. However, nonspecificity of DNA methylation inhibitors may lead to significant side effects. Because RUNX1 induces binding site-specific DNA demethylation, our findings may contribute to the development of highly specific epigenome drugs for cancer therapeutics that induce fewer adverse effects.

Acknowledgments

The authors thank Dave Tang and Marina Lizio for helpful advice on bioinformatics analysis and the members of RIKEN Center for Life Science Technologies Division of Genomic Technologies.

This work was supported by a research grant from the Ministry of Education, Culture, Sport, Science, and Technology of Japan for the RIKEN Center for Life Science Technologies.

Authorship

Contribution: T.S., Y.H., and H.S. were instrumental in the conceptualization; T.S. oversaw the methodology and the software and wrote the original draft; T.S., Y.S., E.F., S.M., M.K., H.N., and S.E. oversaw the investigation; H.S. wrote, reviewed, and edited the paper; and Y.H. and H.S. acquired the funding and supervised the study.

Conflict-of-interest disclosure: The authors declare no competing financial interests.

Correspondence: Harukazu Suzuki, RIKEN Center for Life Science Technologies, 1-7-22, Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan; e-mail: pj.nekir@ikuzus.uzakurah.

References

1. Murry CE, Keller G Differentiation of embryonic stem cells to clinically relevant populations: lessons from embryonic development. Cell. 2008;132(4):661-680. [PubMed]
2. Lacaud G, Gore L, Kennedy M, et al. Runx1 is essential for hematopoietic commitment at the hemangioblast stage of development in vitro. Blood. 2002;100(2):458-466. [PubMed]
3. Harada H, Harada Y, Niimi H, Kyo T, Kimura A, Inaba T High incidence of somatic mutations in the AML1/RUNX1 gene in myelodysplastic syndrome and low blast percentage myeloid leukemia with myelodysplasia. Blood. 2004;103(6):2316-2324. [PubMed]
4. Schnittger S, Dicker F, Kern W, et al. RUNX1 mutations are frequent in de novo AML with noncomplex karyotype and confer an unfavorable prognosis. Blood. 2011;117(8):2348-2357. [PubMed]
5. Cameron S, Taylor DS, TePas EC, Speck NA, Mathey-Prevot B Identification of a critical regulatory site in the human interleukin-3 promoter by in vivo footprinting. Blood. 1994;83(10):2851-2859. [PubMed]
6. Takahashi A, Satake M, Yamaguchi-Iwai Y, et al. Positive and negative regulation of granulocyte-macrophage colony-stimulating factor promoter activity by AML1-related transcription factor, PEBP2. Blood. 1995;86(2):607-616. [PubMed]
7. Huang G, Zhang P, Hirai H, et al. PU.1 is a major downstream target of AML1 (RUNX1) in adult mouse hematopoiesis. Nat Genet. 2008;40(1):51-60. [PubMed]
8. Kitabayashi I, Yokoyama A, Shimizu K, Ohki M Interaction and functional cooperation of the leukemia-associated factors AML1 and p300 in myeloid cell differentiation. EMBO J. 1998;17(11):2994-3004. [PubMed]
9. Oakford PC, James SR, Qadi A, et al. Transcriptional and epigenetic regulation of the GM-CSF promoter by RUNX1. Leuk Res. 2010;34(9):1203-1213. [PubMed]
10. Reed-Inderbitzin E, Moreno-Miralles I, Vanden-Eynden SK, et al. RUNX1 associates with histone deacetylases and SUV39H1 to repress transcription. Oncogene. 2006;25(42):5777-5786. [PubMed]
11. Maier H, Ostraat R, Gao H, et al. Early B cell factor cooperates with Runx1 and mediates epigenetic changes associated with mb-1 transcription. Nat Immunol. 2004;5(10):1069-1077. [PubMed]
12. Watt F, Molloy PL Cytosine methylation prevents binding to DNA of a HeLa cell transcription factor required for optimal expression of the adenovirus major late promoter. Genes Dev. 1988;2(9):1136-1143. [PubMed]
13. Bird A. DNA methylation patterns and epigenetic memory. Genes Dev. 2002;16(1):6-21. [PubMed]
14. Goll MG, Bestor TH Eukaryotic cytosine methyltransferases. Annu Rev Biochem. 2005;74:481-514. [PubMed]
15. Cheng X, Blumenthal RM Mammalian DNA methyltransferases: a structural perspective. Structure. 2008;16(3):341-350. [PMC free article] [PubMed]
16. Ohno R, Nakayama M, Naruse C, et al. A replication-dependent passive mechanism modulates DNA demethylation in mouse primordial germ cells. Development. 2013;140(14):2892-2903. [PubMed]
17. Kohli RM, Zhang Y TET enzymes, TDG and the dynamics of DNA demethylation. Nature. 2013;502(7472):472-479. [PMC free article] [PubMed]
18. Tahiliani M, Koh KP, Shen Y, et al. Conversion of 5-methylcytosine to 5-hydroxymethylcytosine in mammalian DNA by MLL partner TET1. Science. 2009;324(5929):930-935. [PMC free article] [PubMed]
19. Ito S, D’Alessio AC, Taranova OV, Hong K, Sowers LC, Zhang Y Role of Tet proteins in 5mC to 5hmC conversion, ES-cell self-renewal and inner cell mass specification. Nature. 2010;466(7310):1129-1133. [PMC free article] [PubMed]
20. Ito S, Shen L, Dai Q, et al. Tet proteins can convert 5-methylcytosine to 5-formylcytosine and 5-carboxylcytosine. Science. 2011;333(6047):1300-1303. [PMC free article] [PubMed]
21. Maiti A, Drohat AC Thymine DNA glycosylase can rapidly excise 5-formylcytosine and 5-carboxylcytosine: potential implications for active demethylation of CpG sites. J Biol Chem. 2011;286(41):35334-35338. [PMC free article] [PubMed]
22. Hodges E, Molaro A, Dos Santos CO, et al. Directional DNA methylation changes and complex intermediate states accompany lineage specificity in the adult hematopoietic compartment. Mol Cell. 2011;44(1):17-28. [PMC free article] [PubMed]
23. Costa Y, Ding J, Theunissen TW, et al. NANOG-dependent function of TET1 and TET2 in establishment of pluripotency. Nature. 2013;495(7441):370-374. [PMC free article] [PubMed]
24. de la Rica L, Rodríguez-Ubreva J, García M, et al. PU.1 target genes undergo Tet2-coupled demethylation and DNMT3b-mediated methylation in monocyte-to-osteoclast differentiation. Genome Biol. 2013;14(9):R99. [PMC free article] [PubMed]
25. Guilhamon P, Eskandarpour M, Halai D, et al. Meta-analysis of IDH-mutant cancers identifies EBF1 as an interaction partner for TET2. Nat Commun. 2013;4:2166. [PMC free article] [PubMed]
26. Fujiki K, Shinoda A, Kano F, Sato R, Shirahige K, Murata M PPARγ-induced PARylation promotes local DNA demethylation by production of 5-hydroxymethylcytosine. Nat Commun. 2013;4:2262. [PubMed]
27. Nishida H, Suzuki T, Ookawa H, Tomaru Y, Hayashizaki Y Comparative analysis of expression of histone H2a genes in mouse. BMC Genomics. 2005;6:108. [PMC free article] [PubMed]
28. Zhang Y, Liu T, Meyer CA, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137. [PMC free article] [PubMed]
29. Suzuki T, Nakano-Ikegaya M, Yabukami-Okuda H, et al. Reconstruction of monocyte transcriptional regulatory network accompanies monocytic functions in human fibroblasts. PLoS One. 2012;7(3):e33474. [PMC free article] [PubMed]
30. Nikolaidis G, Raji OY, Markopoulou S, et al. DNA methylation biomarkers offer improved diagnostic efficiency in lung cancer. Cancer Res. 2012;72(22):5692-5701. [PMC free article] [PubMed]
31. Ando Y, Tomaru Y, Morinaga A, et al. Nuclear pore complex protein mediated nuclear localization of dicer protein in human cells. PLoS One. 2011;6(8):e23385. [PMC free article] [PubMed]
32. Du P, Zhang X, Huang CC, et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010;11:587. [PMC free article] [PubMed]
33. Melnikova IN, Crute BE, Wang S, Speck NA Sequence specificity of the core-binding factor. J Virol. 1993;67(4):2408-2411. [PMC free article] [PubMed]
34. Jost JP, Jost YC Transient DNA demethylation in differentiating mouse myoblasts correlates with higher activity of 5-methyldeoxycytidine excision repair. J Biol Chem. 1994;269(13):10040-10043. [PubMed]
35. Cortellino S, Xu J, Sannai M, et al. Thymine DNA glycosylase is essential for active DNA demethylation by linked deamination-base excision repair. Cell. 2011;146(1):67-79. [PMC free article] [PubMed]
36. Regha K, Assi SA, Tsoulaki O, Gilmour J, Lacaud G, Bonifer C Developmental-stage-dependent transcriptional response to leukaemic oncogene expression. Nat Commun. 2015;6:7203. [PMC free article] [PubMed]
37. Chen Q, Chen Y, Bian C, Fujiki R, Yu X TET2 promotes histone O-GlcNAcylation during gene transcription. Nature. 2013;493(7433):561-564. [PMC free article] [PubMed]
38. Bottini N, Musumeci L, Alonso A, et al. A functional variant of lymphoid tyrosine phosphatase is associated with type I diabetes. Nat Genet. 2004;36(4):337-338. [PubMed]
39. Taniuchi I, Osato M, Egawa T, et al. Differential requirements for Runx proteins in CD4 repression and epigenetic silencing during T lymphocyte development. Cell. 2002;111(5):621-633. [PubMed]
40. Schmidt T, Karsunky H, Rödel B, Zevnik B, Elsässer HP, Möröy T Evidence implicating Gfi-1 and Pim-1 in pre-T-cell differentiation steps associated with beta-selection. EMBO J. 1998;17(18):5349-5359. [PubMed]
41. Song X, Zhu S, Shi P, et al. IL-17RE is the functional receptor for IL-17C and mediates mucosal immunity to infection with intestinal pathogens. Nat Immunol. 2011;12(12):1151-1158. [PubMed]
42. Higuchi Y, Yasui A, Matsuura K, Yamamoto S CD156 transgenic mice. Different responses between inflammatory types. Pathobiology. 2002;70(1):47-54. [PubMed]
43. Jongstra-Bilen J, Jongstra J Leukocyte-specific protein 1 (LSP1): a regulator of leukocyte emigration in inflammation. Immunol Res. 2006;35(1-2):65-74. [PubMed]
44. Perez-Caballero D, Zang T, Ebrahimi A, et al. Tetherin inhibits HIV-1 release by directly tethering virions to cells. Cell. 2009;139(3):499-511. [PMC free article] [PubMed]
45. Kundaje A, Meuleman W, Ernst J, et al. ; Roadmap Epigenomics Consortium. Integrative analysis of 111 reference human epigenomes. Nature. 2015;518(7539):317-330. [PMC free article] [PubMed]
46. Sandelin A, Alkema W, Engström P, Wasserman WW, Lenhard B JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004;32(Database issue):D91-D94. [PMC free article] [PubMed]
47. Okuda T, Takeda K, Fujita Y, et al. Biological characteristics of the leukemia-associated transcriptional factor AML1 disclosed by hematopoietic rescue of AML1-deficient embryonic stem cells by using a knock-in strategy. Mol Cell Biol. 2000;20(1):319-328. [PMC free article] [PubMed]
48. Schotte R, Rissoan MC, Bendriss-Vermare N, et al. The transcription factor Spi-B is expressed in plasmacytoid DC precursors and inhibits T-, B-, and NK-cell development. Blood. 2003;101(3):1015-1023. [PubMed]
49. Sekiya T, Kashiwagi I, Inoue N, et al. The nuclear orphan receptor Nr4a2 induces Foxp3 and regulates differentiation of CD4+ T cells. Nat Commun. 2011;2:269. [PMC free article] [PubMed]
50. Lohoff M, Mak TW Roles of interferon-regulatory factors in T-helper-cell differentiation. Nat Rev Immunol. 2005;5(2):125-135. [PubMed]
51. Liebermann DA, Gregory B, Hoffman B AP-1 (Fos/Jun) transcription factors in hematopoietic differentiation and apoptosis. Int J Oncol. 1998;12(3):685-700. [PubMed]
52. Forrest AR, Kawaji H, Rehli M, et al. ; FANTOM Consortium and the RIKEN PMI and CLST (DGT). A promoter-level mammalian expression atlas. Nature. 2014;507(7493):462-470. [PMC free article] [PubMed]
53. Schoenherr CJ, Levorse JM, Tilghman SM CTCF maintains differential methylation at the Igf2/H19 locus. Nat Genet. 2003;33(1):66-69. [PubMed]
54. Wang H, Maurano MT, Qu H, et al. Widespread plasticity in CTCF occupancy linked to DNA methylation. Genome Res. 2012;22(9):1680-1688. [PubMed]
55. Hashimoto H, Liu Y, Upadhyay AK, et al. Recognition and potential mechanisms for replication and erasure of cytosine hydroxymethylation. Nucleic Acids Res. 2012;40(11):4841-4849. [PMC free article] [PubMed]
56. Valinluck V, Sowers LC Endogenous cytosine damage products alter the site selectivity of human DNA maintenance methyltransferase DNMT1. Cancer Res. 2007;67(3):946-950. [PubMed]
57. Kubosaki A, Tomaru Y, Furuhata E, et al. CpG site-specific alteration of hydroxymethylcytosine to methylcytosine beyond DNA replication. Biochem Biophys Res Commun. 2012;426(1):141-147. [PubMed]
58. Li C, Lan Y, Schwartz-Orbach L, et al. Overlapping requirements for Tet2 and Tet3 in normal development and hematopoietic stem cell emergence. Cell Reports. 2015;12(7):1133-1143. [PMC free article] [PubMed]
59. Deplus R, Delatte B, Schwinn MK, et al. TET2 and TET3 regulate GlcNAcylation and H3K4 methylation through OGT and SET1/COMPASS. EMBO J. 2013;32(5):645-655. [PMC free article] [PubMed]
60. Liddiard K, Hills R, Burnett AK, Darley RL, Tonks A OGG1 is a novel prognostic indicator in acute myeloid leukaemia. Oncogene. 2010;29(13):2005-2012. [PubMed]
61. Krejci O, Wunderlich M, Geiger H, et al. p53 signaling in response to increased DNA damage sensitizes AML1-ETO cells to stress-induced death. Blood. 2008;111(4):2190-2199. [PubMed]
62. Ko M, An J, Bandukwala HS, et al. Modulation of TET2 expression and 5-methylcytosine oxidation by the CXXC domain protein IDAX. Nature. 2013;497(7447):122-126. [PMC free article] [PubMed]
63. Meissner A, Mikkelsen TS, Gu H, et al. Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature. 2008;454(7205):766-770. [PMC free article] [PubMed]
64. Domcke S, Bardet AF, Adrian Ginno P, Hartl D, Burger L, Schübeler D Competition between DNA methylation and transcription factors determines binding of NRF1. Nature. 2015;528(7583):575-579. [PubMed]

Articles from Blood Advances are provided here courtesy of American Society of Hematology