Research | Open | Published:
PU.1 target genes undergo Tet2-coupled demethylation and DNMT3b-mediated methylation in monocyte-to-osteoclast differentiation
Genome Biologyvolume 14, Article number: R99 (2013)
DNA methylation is a key epigenetic mechanism for driving and stabilizing cell-fate decisions. Local deposition and removal of DNA methylation are tightly coupled with transcription factor binding, although the relationship varies with the specific differentiation process. Conversion of monocytes to osteoclasts is a unique terminal differentiation process within the hematopoietic system. This differentiation model is relevant to autoimmune disease and cancer, and there is abundant knowledge on the sets of transcription factors involved.
Here we focused on DNA methylation changes during osteoclastogenesis. Hypermethylation and hypomethylation changes took place in several thousand genes, including all relevant osteoclast differentiation and function categories. Hypomethylation occurred in association with changes in 5-hydroxymethylcytosine, a proposed intermediate toward demethylation. Transcription factor binding motif analysis revealed an over-representation of PU.1, NF-κB, and AP-1 (Jun/Fos) binding motifs in genes undergoing DNA methylation changes. Among these, only PU.1 motifs were significantly enriched in both hypermethylated and hypomethylated genes; ChIP-seq data analysis confirmed its association to both gene sets. Moreover, PU.1 interacts with both DNMT3b and TET2, suggesting its participation in driving hypermethylation and hydroxymethylation-mediated hypomethylation. Consistent with this, siRNA-mediated PU.1 knockdown in primary monocytes impaired the acquisition of DNA methylation and expression changes, and reduced the association of TET2 and DNMT3b at PU.1 targets during osteoclast differentiation.
The work described here identifies key changes in DNA methylation during monocyte-to-osteoclast differentiation and reveals novel roles for PU.1 in this process.
DNA methylation plays a fundamental role in differentiation as it drives and stabilizes gene activity states during cell-fate decisions. Recent reports have shown a close relationship between the participation of transcription factors during differentiation and the generation of cell type-specific epigenetic signatures [1–3]. Several mechanisms explain the co-occurrence of DNA methylation changes and transcription factor binding, including the active recruitment of enzymes involved in DNA methylation deposition, interference, or alternative use of the same genomic regions. One of the best models for investigating these mechanisms is the hematopoietic differentiation system given the profound knowledge on the transcription factors implicated at different stages. Many studies have focused on hematopoiesis in order to learn about the type, distribution, and role of epigenetic changes, particularly DNA methylation during differentiation. However, the role of DNA methylation changes and the mechanisms participating in their acquisition in terminal differentiation processes remain elusive, even though these are among the most important since they produce functional cell types with very specific roles.
A singular differentiation process within the hematopoietic system is represented by differentiation from monocytes (MOs) to osteoclasts (OCs), which are giant, multinucleated cells that are specialized in degrading bone . OCs differentiate from monocyte/macrophage progenitors following M-CSF  and RANKL  stimulation. Osteoclastogenesis requires cell fusion, cytoskeleton re-organization  and the activation of the specific gene sets necessary for bone catabolism. The signaling pathways activated after M-CSF and RANKL induction have been extensively described, and act through TRAF-6 [8, 9], immunoreceptor tyrosine-based activation motif (ITAM)  adaptors DAP12  and FcRγ  associated with their respective receptors, TREM-2  and OSCAR, as well as calcium oscillations . Signals end in the activation of NF-kB, MAPK, and c-Jun, leading to the activation of NFATc1 , the master transcription factor of osteoclastogenesis, together with PU.1 and MITF , which is already present in the progenitors. These transcription factors bind to the promoter and help upregulating OC markers such as dendritic cell-specific transmembrane protein (DC-STAMP/TM7SF4) , tartrate-resistant acid phosphatase (TRACP/ACP5) , cathepsin K (CTSK) , matrix metalloproteinase 9 (MMP9) , and carbonic anhydrase 2 (CA2).
OC deregulation is involved in several pathological contexts, either in the form of deficient function, as in the case in osteopetrosis , or aberrant hyperactivation, as in osteoporosis . These cells are also involved in autoimmune rheumatic disease. For instance, in rheumatoid arthritis aberrantly activated OCs are major effectors of joint destruction . Moreover, OCs cause bone complications in several diseases, such as multiple myeloma , prostate cancer, and breast cancer , and there is also a specific tumor with OC origin, the giant cell tumor of bone .
In vitro generation of OCs allows this cell type to be investigated, whereas isolating primary bone OCs for this purpose is very difficult. MOs stimulated with RANKL and M-CSF generate functional OCs , which degrade bone and express OC markers . As indicated, the involvement of transcription factors in this model has been well studied, however very few reports have analyzed the role of epigenetic changes during osteoclastogenesis, and these focus mainly on histone modifications [29, 30]. Given the relationship between transcription factors and DNA methylation, we hypothesized that examining DNA methylation changes would provide clues about the involvement of specific factors in the dynamics and hierarchy of these changes in terminal differentiation.
In this study, we compared the DNA methylation profiles of MOs and derived OCs following M-CSF and RANKL stimulation. We found that osteoclastogenesis was associated with the drastic reshaping of the DNA methylation landscape. Hypermethylation and hypomethylation occur in many relevant functional categories and key genes, including those whose functions are crucial to OC biology, like CTSK, ACP5, and DC-STAMP. Hypomethylation occurred early, concomitantly with transcription changes, was DNA replication-independent and associated with a change in 5-hydroxymethylcytosine, which has been proposed as an intermediate in the process of demethylation. Inspection of transcription factor binding motif over-representation in genes undergoing DNA methylation changes revealed the enrichment of the PU.1 binding motif in hypermethylated genes and AP-1, NF-kB, and also PU.1 motifs among hypomethylated genes. In fact, analysis of PU.1 ChIPseq data showed its general association to a high number of both hypo- and hypermethylated sites. Chromatin immunoprecipitation assays and immunoprecipitation experiments suggested a potential novel role for PU.1 recruiting DNMT3B to hypermethylated promoters, and TET2, which converts 5-methylcytosine to 5-hydroxymethylcytosine, to genes that become demethylated. This has been demonstrated by performing siRNA-mediated downregulation of PU.1 which partially impaired DNA methylation, expression, and recruitment of TET2 and DNMT3B to PU.1 targets, supporting the participation of PU.1 in the acquisition of DNA methylation changes at their target sites.
Cell differentiation and fusion in osteoclastogenesis are accompanied by hypomethylation and hypermethylation of key functional pathways and genes
To investigate the acquisition of DNA methylation changes during monocyte-to-osteoclast differentiation we first obtained three sets of matching samples corresponding to MOs (CD14+ cells) from peripheral blood and OCs derived from the same CD14+ cells, 21 days after the addition of M-CSF and RANKL. The quality of mature, bone-resorbing OCs obtained under these conditions was confirmed by several methods, including the presence of more than three nuclei in TRAP-positive cells (in some cases, up to 40 nuclei per cell were counted), the upregulation of OC markers, such as CA2, CTSK, ACP5/TRACP, and MMP9, and downregulation of the monocytic gene CX3CR1 (Additional file 1). At 21 days, over 84% of the nuclei detected in these preparations could be considered to be osteoclastic nuclei (in polykaryons, nuclei and not cells were counted) (Additional file 1). We then performed DNA methylation profiling using bead arrays that interrogate the DNA methylation status of >450,000 CpG sites across the entire genome covering 99% of RefSeq genes. Statistical analysis of the combined data from the three pairs of samples revealed that 3,515 genes (8,028 CpGs) displayed differential methylation (FC ≥2 or FC ≤0.5; false discovery rate (FDR) ≤0.05). Specifically, we identified 1,895 hypomethylated genes (3,597 CpG sites) and 2,054 hypermethylated genes (4,429 CpGs) (Figure 1A and Additional file 2). Changes corresponding to the average three pairs of monocytes/osteoclasts (Figure 1B) were almost identical to the pattern obtained for each individual pair of samples (Additional file 3), highlighting the specificity of the differences observed.
Over one-third of the differentially methylated CpG-containing probes (33% for hypomethylated CpGs, 45% for hypermethylated CpGs) mapped to gene promoters, the best-described regulatory region for DNA methylation, although DNA methylation changes also occurred at a similar scale in gene bodies (51% for hypomethylated CpGs, 40% for hypermethylated CpGs) (Figure 1C). Gene ontology analysis of hypomethylated CpGs revealed significant enrichment (FDR ≤0.05) for a variety of functional categories of relevance in OC differentiation and function (Figure 1D). We observed very high significance for terms like immune response (FDR = 4.25E-25) and signal transduction (FDR = 1.45E-21), but also more specific categories such as ruffle organization (FDR = 9.91E-2), calcium ion transport (FDR = 4.6E-2), and OC differentiation (FDR = 1.94E-1). In the case of hypermethylated genes, we also found highly significant enrichment of signal transduction (FDR = 4.09E-17), and enrichment of categories related to other hematopoietic cell types, suggesting that hypermethylation and associated silencing take place in gene sets that become silent in differentiated OCs (Figure 1D). Together, these data indicate that DNA hypomethylation is targeted to genomic regions that are activated during osteoclastogenesis, and hypermethylation silences alternative lineage genes that are not expressed in OCs.
Remarkably, among the group of hypomethylated genes (Additional file 2), we identified changes in several of the archetypal OC genes near their transcription start sites. For example, CTSK, the lysosomal cysteine proteinase involved in bone remodeling and resorption, is hypomethylated more than 60%. The ACP5/TRACP gene is hypomethylated around 47%. Finally, TM7SF4, which encodes for DC-STAMP, a seven-pass transmembrane protein involved in signal transduction in OCs and dendritic cells, undergoes 59% hypomethylation. We also observed significant hypomethylation at the osteoclast-specific transcription factor gene NFATC1, although in this case hypomethylation occurred at CpG sites located in its gene body region. Conversely, CX3CR1, an important factor for MO adhesion to blood vessels that is downregulated during osteoclastogenesis, displayed an increase in methylation of over 28% (Additional file 2).
To confirm that differences in DNA methylation identified between MOs and OCs were robust, we carried out bisulfite genomic pyrosequencing of the aforementioned selection of genes, looking at CpG sites corresponding to the oligonucleotide probe represented in the methylation array. In all cases, bisulfite pyrosequencing confirmed the results of the beadchip array (Figure 1E and Additional file 4). This analysis showed a very close correlation between the array and the pyrosequencing data (R2 = 0.9707) (Figure 1E).
We also investigated the coordinated hypomethylation or hypermethylation of adjacent CpGs by analyzing the different sequence window lengths (from 500 bp to 1,000,000 bp). With the largest sequence windows we were able to observe the coordinated hypermethylation of multiple CpGs across several genes, like those in the HOXA gene cluster. However, the majority of CpGs undergoing coordinated methylation changes were identified within the single gene level. By analyzing CpGs that are concomitantly deregulated within a 500-bp window, we identified several genes displaying coordinated hypomethylation or hypermethylation of many CpG sites (Additional file 5). Among these, we identified several CpG clusters in genes potentially involved in OC function and/or differentiation, including 10 CpGs at the promoter of the TM4SF19 gene, also known as OC maturation-associated gene 4 protein, and nine CpGs in the gene body of ARID5B, the AT-rich interactive domain 5B (MRF1-like) (Figure 1F).
To examine the specificity of the DNA methylation changes further we performed bisulfite sequencing of repetitive elements (Sat2, D4Z4, and NBL2 repeats) and ribosomal RNA genes (Figure 1G and Additional file 3). We also performed genome-wide amplification of unmethylated DNA Alu repeats (AUMA), the most common family of repetitive elements that are present in tandem or interspersed in the genome . These experiments showed no significant DNA methylation changes in any of these repetitive elements (Additional file 3), reinforcing the notion of the high specificity of hypomethylation and hypermethylation of the identified gene sets.
Hypomethylation is replication-independent and involves changes in 5-hydroxymethylcytosine
To investigate the dynamics of DNA methylation in relation to gene expression changes we first examined how DNA methylation changes are associated with expression changes and then compared the dynamics of DNA methylation and expression changes.
We used osteoclastogenesis expression data (available from the ArrayExpress database under accession number E-MEXP-2019) on 0, 5, and 20 days . Our analysis showed that most changes occurred within the first 5 days, since the expression changes between 0 and 5 days were very similar to those observed between 0 and 20 days, and very few genes changed between 5 and 20 days (Figure 2A and Additional file 6). The 0-to-20-day comparison showed that 2,895 genes were upregulated (FC >2; FDR <0.05) and 1,858 were downregulated (FC <0.5; FDR <0.05). We found different relationships between DNA methylation changes and gene expression (Figure 2B). An inverse relationship between DNA methylation and gene expression was mainly observed for changes occurring in CpGs in the proximity of the TSS and within the first exon (Figure 2C) and it was less frequent in those at gene bodies and 3′UTR (Figure 2C). Comparing DNA methylation and expression data revealed that 452 genes were both hypomethylated and overexpressed and 280 genes were both hypermethylated and repressed at the selected thresholds (Additional file 7). We selected a panel of 10 genes from those undergoing hypomethylation and hypermethylation to investigate the dynamics of DNA methylation and expression changes, and performed bisulfite pyrosequencing and quantitative RT-PCR over the entire osteoclastogenesis for three sets of samples (Figure 2D). We found that the promoters of genes like ACP5, CTSK, TM7SF4, and TM4SF19 rapidly became hypomethylated following RANKL and M-CSF stimulation (Figure 2D, top). In fact, around 60% of the entire range of hypomethylation occurred between days 0 and 4. Changes in mRNA levels occurred at a similar pace or, in some cases, in an even more gradual manner and were slightly delayed with respect to changes in DNA methylation. In contrast, hypermethylated genes like PPP1R16B, CD6, and NR4A2 (Figure 2D, bottom) displayed loss of expression before experiencing an increase in DNA methylation, highlighting the different dynamics and mechanisms involved in hypomethylation and hypermethylation events.
It is well established that osteoclastogenesis occurs in the absence of cell division. We tested the levels of cell division in our monocyte-to-osteoclast differentiation experiments by treating cells with BrdU pulses. Consistent with previous observations, fewer than 9.8% were found to be BrdU-positive between 1 and 4 days, confirming the virtual absence of replication (Figure 2E and Additional file 8). This implies that the large DNA methylation changes observed in this time period are independent of DNA replication. This conclusion is also supported by the fact that treatment with 5-Aza-2′-deoxycytidine (5azadC), a pharmacological compound that results in replication-coupled DNA demethylation , had no significant effect on osteoclastogenesis (Figure 2F).
The existence of DNA methylation changes in the absence of replication is particularly significant for genes undergoing demethylation, given the controversy around active DNA demethylation mechanisms. In this context, recent studies have drawn attention towards a family of enzymes, the Tet proteins, which convert 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) [34, 35] and other modified forms of cytosine, 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC) . 5hmC, 5fC, and 5caC may represent intermediates in an active demethylation pathway that ultimately replaces 5mC with cytosine in non-dividing cells [37, 38]. To establish the potential involvement of these mechanisms, we here focused on the 5hmC levels at early time points in several of the genes that are hypomethylated during osteoclastogenesis, using a method that cleaves DNA that has C, 5mC, but not the glucosyl-5hmC produced as a result of treatment with the 5-hydroxymethylcytosine specific glucosyltransferase enzyme (Figure 2G). For several genes that become hypomethylated, like ACP5 and TM4SF19, we observed an initial increase in 5hmC levels followed by a slight but significant decrease (Figure 2H). In other genes, like TM7SF4 and CD59, of there were high levels of 5hmC before the addition of RANKL/M-CSF as if these genes were already primed for demethylation. In any case, our results suggested the participation of hydroxymethylation, and therefore the activity of Tet proteins, in genes that undergo a reduction in DNA methylation.
Sequences undergoing DNA methylation changes are enriched for binding motifs for AP-1, NF-kB, and PU.1, key transcription factors in osteoclastogenesis
Different studies have recently shown that transcription factor binding events are associated with changes in the DNA methylation profiles and the response to different situations [2, 39, 40]. To address this further, we first investigated the potential over-representation of transcription factor binding motifs among the sequences undergoing DNA methylation changes during OC differentiation using the TRANSFAC database and focusing on a region of 500 bp around the CpG sites identified as undergoing hypomethylation or hypermethylation. We noted highly significant overrepresentation of a small selection of transcription factor binding motifs for genes that undergo hypomethylation or hypermethylation (Figure 3A). We observed that the over-representation of binding motifs was very specific to the direction of the DNA methylation change (hypomethylation or hypermethylation).
In the case of hypomethylated genes, we observed highly significant enrichment of binding motifs of the AP-1 family and NF-kB subunits (Figure 3A). We also observed enrichment of PU.1 (FDR 1.07E-12). In fact, 39% of all hypomethylated genes had binding motifs for AP-1, 15% genes had NF-kB binding motifs, and another 15% genes had binding motifs for PU.1 or other ETS-related factors (PU.1 alone, 9%) (Figure 3B). As aforementioned, these three groups of TFs play critical roles in osteoclastogenesis . For instance, c-Fos, a component of the dimeric TF AP-1, regulate the switch between monocytes/macrophages and OC differentiation. Fra-1 is downstream to c-Fos, whereas PU.1 and NF-kB are upstream. NF-kB is critical in the expression of a variety of cytokines involved in OC differentiation. In the case of hypermethylated genes, we identified even greater enrichment of the binding motifs of ETS-related transcription factors, especially PU.1 (Figure 3A). In fact, the PU.1 binding motif is present in 15% of all hypermethylated genes (Figure 3B). Other motifs of ETS-related transcription factors from our list of hypermethylated genes included SPIB, ESE1, ETS1, ETS2, and others (Figure 3A). Much lower or insignificant levels of enrichment were obtained for AP-1 family members and NF-kB subunits among the hypermethylated genes. Previous studies have shown that genes that become methylated during hematopoietic differentiation are characterized by ETS transcription factors . This appears to be particularly relevant in monocytic differentiation . Interestingly, most of the reports about the role of PU.1 in osteoclastogenesis are associated with the activation of osteoclast-specific genes. However, in relation with methylation changes, PU.1 appears to be better correlated with those changes in the direction of repression. Overall, the analysis of transcription factor motifs showed that several of the factors associated with osteoclastogenesis had a significant over-representation of their binding motifs among the sets of hypo-and hypermethylated genes (Additional file 9).
To confirm the association of some of these factors with genes that become hypo- and hypermethylated we performed chromatin immunoprecipitation (ChIP) assays with a selection of transcription factors including PU.1, the NFkB subunit p65 and c-Fos, given their known role in osteoclastogenesis as well as the presence of binding motifs for them around hypo- and hypermethylated genes (in the case of c-Fos, it was chosen as a component of the dimeric TF AP-1). To select candidate genes we considered genes with motifs for these three factors among the list of hypo- and hypermethylated genes. For instance, genes that become hypomethylated and have binding sites for p65 include CCL5, IL1R, and TNFR5SF. In the case of transcription factor c-Fos, we looked at genes that become hypomethylated like IL7R, CD59, and IL1R. In the case of PU.1, we chose key genes with PU.1 binding near the differentially methylated CpG, including ACP5 and TM4SF7 (hypomethylated) and CX3CR1 (hypermethylated). ChIP assays demonstrated the interaction of these factors with most of the aforementioned promoters, even before the stimulation with M-CSF and RANKL (Figure 3C), as if these genes were primed by these factors in monocytes. No binding was observed in control sequences like Sat2 repeats and the MYOD1 and TDRD1 promoters. Interestingly, in the case of PU.1, with both hypo- and hypermethylated genes displaying binding at 0 days, genes becoming demethylated showed a slight increase in PU.1 binding, whereas hypermethylated genes showed a slight decrease in PU.1 association (Figure 3C).
PU.1 recruits DNMT3b and TET2 to hypermethylated and hypomethylated genes
To investigate the potential role of the aforementioned transcription factors in the acquisition of DNA methylation changes we chose PU.1 and NF-kB (p65 subunit) as two representative examples. We first checked their expression levels during osteoclastogenesis, by carrying out qRT-PCR and western blot assays. mRNA and protein analysis (Figure 4A and 4B) both confirmed the expression of these factors. PU.1 revealed an increase at the mRNA levels, although there was no change at the protein level. In the case of p65 NF-kB we only observed a clear increase at the protein level (Figure 4B). In parallel, we also confirmed the presence of DNMT3b, a de novo DNA methyltransferase, and the ten eleven translocation (TET) protein TET2, as enzymatic activities potentially related with DNA demethylation (Figure 4A and 4B). TET proteins are responsible for conversion of 5mC in 5hmC , 5fC, and 5cac . Recent evidences support a role for TET-dependent active DNA demethylation process [42, 43]. We focused on TET2 given their high levels in hematopoietic cells of myeloid origin [44, 45]. Also, we have recently reported that TET2 plays a role in derepressing genes in pre-B cell to macrophage differentiation , and recent data shows that TET2 is required for active DNA demethylation in primary human MOs . In fact TET1 and TET3 were undetectable in western blot (not shown) and qRT-PCR evidenced their low levels in this cell type (Figure 4A, only shown for TET1).
The confirmed binding of factors like PU.1 and the p65 subunit of NF-kB to hypo- and hypermethylated genes (Figure 3C) raised the possibility of their potential direct interaction with factors involved in maintaining the DNA methylation homeostasis. Some of these interactions have already been explored. For instance, PU.1 physically interacts with the de novo DNA methyltransferases DNMT3A and DNMT3B . Such an interaction, if it occurred in osteoclastogenesis, could provide a potential mechanism to explain how PU.1 target genes become hypermethylated. One would expect that these transcription factors could also interact with factors participating in demethylation processes. Our previous results suggested the existence of 5hmC enrichment in genes that become hypomethylated, and therefore it is reasonable to test whether NF-kB p65 and PU.1 interact with Tet proteins, the enzymes catalyzing hydroxylation of 5mC.
We therefore tested the recruitment by NF-kB p65 and PU.1 of both DNMT3b and TET2 by carrying out immunoprecipitation assays with osteoclastogenesis samples 0, 2, and 4 days after stimulation with M-CSF and RANKL. Our results showed that PU.1 directly interacted with both DNMT3b and TET2 (Figure 4C). It is plausible that these two interactions may involve different subpopulations of PU.1, for instance with specific post-translational modifications like Ser phosphorylation. However, we did not address this aspect at this point. In the case of NF-kB, we did not observe binding with either of these factors (Figure 4C). This could perhaps be explained by the fact that p65 is shuttling back to the cytoplasm much of the time.
To confirm the interaction between PU.1 and DNMT3b and TET2, we performed reciprocal immunoprepitation experiments with anti-DNMT3b and anti-TET2. These confirmed the direct interaction with PU.1 (Figure 4C). Our results suggested that PU.1 may play a dual coupling transcription factor that can interact with the DNA methyltransferases and enzymes perhaps participating in, or leading to, demethylation. It is likely that other factors participate in the recruitment of these enzymes, however at this stage we focused on PU.1 because of its ability to bind both DNMT3b and Tet2 and its association with both hyper- and hypomethylated sequences.
We then investigated the dual role of PU.1 in recruiting TET2 and DNMT3b to promoters. To this end we performed chromatin immunoprecipitation assays with PU.1, TET2, and DNMT3b in MOs at 0 and 2 days following stimulation with M-CSF/RANKL. We amplified gene promoters with predicted binding sites for PU.1 that become both demethylated (ACP5, TMS7SF4, and TM4SF19) as well as hypermethylated (CX3CR1 and NR4A2) and used a non-target of PU.1 (MYOD1) as negative control. Our analysis showed binding of PU.1 at both 0 and 2 days (Figure 4D). For genes that become hypomethylated, we observed an increased recruitment of TET2 at these promoters after 2 days, whereas DNMT3b was initially enriched but its association with these promoters was lost after M-CSF and RANKL stimulation (Figure 4D). In genes that become hypermethylated (CX3CR1 and NR4A2), we also observed association of PU.1 at both 0 and 2 days. However we again observed a slight decrease at 2 days together. We also observed increased recruitment of DNMT3b at 2 days after M-CSF/ RANKL stimulation. We did not observe association of PU.1, DNMT3b, and TET2 in the negative control for PU.1 binding, the MyoD promoter.
To evaluate the extent to which hypo- and hypermethylated genes correlate with PU.1 occupancy, we used our DNA methylation data and PU.1 ChIPseq data (GSE31621) obtained in MOs . Most of the individual example genes previously analyzed displayed PU.1 binding overlapping or in the proximity of the CpG sites undergoing a methylation change (Figure 4E). To systematize the analysis we used a window of 500 bp centered around the CpG displaying DNA methylation changes. Under these conditions we found that 10.7% of all hypomethylated CpGs located in promoter regions genes and 25.1% of all hypermethylated CpGs located in promoter regions had PU.1 peaks within this 500-bp window (Figure 4F). These numbers were similar when focusing on CpGs located in distal regions (Figure 4F). We also compared the ChIPseq data to the prediction by TRANSFAC analysis and observed that the overlap between the two sets of list was around 20.6% for hypomethylated genes and 46.9% for hypermethylated genes, again using the same 500 bp for both datasets (Figure 4F). These analyses reinforced the notion of PU.1 associated with a high number of genes undergoing DNA methylation changes; however, it also reveals the weakness in the predictive power of TRANSFAC motif searches and the need of experimental validation of its results.
Dowregulation of PU.1 in MOs impairs activation of OC markers, hypomethylation, and recruitment of DNMT3b and TET2
To investigate a potential causal relationship between PU.1 and DNA methylation changes in monocyte-to-osteoclast differentiation we investigated the effects of ablating PU.1 expression in MOs. We therefore downregulated PU.1 levels in MOs using transient transfection experiments with a mix of two siRNAs targeting exon2 and the 3′UTR of PU.1 (Figure 5A). In parallel, we used a control siRNA. Following transfection we stimulated differentiation using RANKL/M-CSF. In these conditions, we checked by qRT-PCR and western blot the effects on PU.1 levels at 1, 2, 4, and 6 days following RANKL/M-CSF stimulation of MOs and confirmed the PU.1 downregulation close to 60% (Figure 5B and 5C). We then observed that the upregulation of genes like ACP5 and CTSK (both PU.1-direct targets) was partially impaired (Figure 5D). In the case of genes like CX3CR1 and NR4A2 we determined that downregulation was also impaired in PU.1-siRNA treated MOs. Interestingly, we also analyzed two genes that are not direct PU.1 targets, one upregulated and hypomethylated during osteoclastogenesis (PLA2G4E) and a second one, highly methylated, that does not experience DNA methylation changes during OC differentiation (FSCN3). PU.1-siRNA treatment had only small effects on gene expression changes during osteoclast differentiation (perhaps due to indirect effects) when compared to control siRNA, confirming the specificity of the changes observed for the other genes (Figure 5D, bottom).
We then tested the effects of PU.1 downregulation in DNA methylation changes. We looked at both PU.1-target genes that become hypomethylated and hypermethylated. In both cases, we observed that downregulation of PU.1 impaired the acquisition of DNA methylation changes, in contrast with the changes observed for control siRNA-treated MOs (Figure 5D). In the case of TM7SF4, one of the key genes undergoing hypomethylation, we did not detect an effect of PU.1 downregulation on its DNA methylation dynamics (Additional file 10). However, this could perhaps be explained because this gene undergoes changes before downregulation of PU.1 by siRNA is effective, within day 1 (Additional file 10) and suggests the participation of other factors in this process. At any rate, the observed effects only occurred in PU.1 targets. It did not affect genes that are not targeted by PU.1 (PLA2G4E, FSCN3). In the case of PLA2G4E, PU.1-siRNA treatment did not impair the loss of methylation that occurred in the control experiment. For FSCN3, we observed no loss of methylation in any case (Figure 5D).
Finally, we compared the effect of PU.1 downregulation in the recruitment of DNMT3b and TET2 to hyper- and hypomethylated promoters (Figure 5E and Additional file 10). As expected, we observed that PU.1 dowregulation resulted in a decrease of the levels of PU.1 associated with the promoters of both hypo- and hypermethylated genes. Most importantly, it also reduced the association of DNMT3b and TET2 reinforcing the notion of the role of these factors and their association with PU.1 in the DNA methylation changes occurring at these CpG sites (Figure 5E). The time course analysis (at 2 and 6 days) of these results also revealed a complex dynamics for the PU.1, TET2, and DNMT3b interactions with their target genes, particularly in the case of hypermethylated genes. It is possible that perturbation of PU.1 levels could be compensated by additional factors that participate in the acquisition of DNA methylation changes of these genes. These aspects will need to be further investigated.
Our results provide evidence of the participation of transcription factors, focusing on PU.1, in determining changes in DNA methylation during monocyte-to-osteoclast differentiation. First, a detailed analysis of the sequences undergoing DNA methylation changes produced evidences of the participation of several transcription factors, given the specific over-representation of certain motifs in hypo- and hypermethylated genes. This initial analysis was validated in several candidate genes and using ChIPseq data for human primary monocytes . Second, further analyses on one these candidate transcription factors, PU.1, and manipulation of its levels revealed a novel role for this factor in mediating DNA methylation changes during osteoclastogenesis, by direct binding of both DNMT3B and TET2.
In general, DNA methylation changes in differentiation or any other dynamic process are of interest for two reasons: (1) these changes are generally associated with gene expression changes, particularly when associated with promoters or gene bodies, and reveal aspects intrinsic to identity and function of the corresponding cell types; (2) they can be considered as epigenetic footprints that, despite not necessarily being associated with an expression or organizational change, reveal a change in the milieu of a particular CpG and therefore can be used to trace the participation of specific transcription factors or other nuclear elements in that environment/neighborhood. This information can then be used to reconstruct cell signaling events, transcription factors involved and mechanisms participating in differentiation. In this sense, our data show that DNA methylation changes are involved in the differentiation dynamics and stabilization of the OC phenotype since they are concomitant with, or even precede, expression changes. These data are closely correlated with gene expression changes, and a majority of genes that undergo hypomethylation or hypermethylation at their promoters or gene bodies also experience a change in expression, although the relationship varies between different gene sets. Finally, gene ontology analysis reveals that all relevant functional categories and the majority of key genes for differentiation or the activity of functional OCs undergo DNA methylation changes and that genes within all relevant functional categories undergo DNA methylation changes.
Our study suggests that both hypomethylation and hypermethylation events are equally important. Hypomethylation events, in many cases associated with gene activation, affect genes that are specific to this differentiation process or are related with the function of differentiated OCs. In contrast, the identity of genes affected by hypermethylation events is less closely correlated with OC function, given that most of them are related with gene repression. In fact, we found that hypermethylation affects genes that are specific to other cellular types. Given that osteoclastogenesis involves cell fusion and the generation of highly polyploid cells, we had speculated whether the existence of redundant copies of genetic material could lead to massive gene repression, and the silencing of extra copies. However, hypermethylation does not seem to be predominant over hypomethylation. The two activities are very specific to particular gene sets and there are no indications of changes in repetitive elements.
A number of transcription factors are essential for OC formation. Some of these factors are involved in various differentiation processes. Among these, PU.1, c-Fos, NF-kB, and other factors are essential for osteoclastogenesis. In fact, NF-kB- and PU.1-deficient mice show a macrophage differentiation failure, and osteoclastogenesis is inhibited at an early stage of differentiation. c-Fos is a component of the dimeric TF AP-1, which also includes FosB, Fra-1, Fra-2, and Jun proteins such as c-Jun, JunB, and JunD. Other key factors involved in OC differentiation include C/EBPα  and Bach1 . Osteoclastogenesis also depends on the activity of more specific transcription factors like NFATc1 and MITF. Interestingly, the analysis of the presence of transcription factor binding sites in sequences that undergo DNA methylation changes shows a significant enrichment in binding motifs of transcription factors that are key in OC differentiation, some of which we have validated for a selection of putative target genes.
One of the most interesting factors in this process is the ETS factor PU.1. In fact, PU.1 is the earliest molecule known to influence the differentiation and commitment of precursor myeloid cells to the OC lineage. PU.1 functions in concert with other transcription factors, including c-Myb, C/EBPα, cJun, and others, to activate osteoclast-specific genes.
Our results reveal two hitherto undescribed roles for PU.1 in the context of monocyte-to-OC differentiation. First, we have identified the association of PU.1 with genes that become repressed through hypermethylation and describe its direct interaction with DNMT3b in the context of osteoclastogenesis. Second, we identify a novel interaction between PU.1 and TET2 and their association with genes that become demethylated. Our study shows that PU.1 may act as a dual adaptor during osteoclastogenesis, in the directions of hypomethylation and hypermethylation. This is compatible with previous data on genome wide DNA methylation profiling comparing cell types across the hematopietic differentiation system where an over-representation of ETS transcription factor binding sites was found . In monocyte-to-osteoclast differentiation, PU.1 is best known for its role in the activation of osteoclast-specific genes. However, studies in other models have previously shown that PU.1 can participate in the repression of genes in concert with elements of the epigenetic machinery. For instance, PU.1 is known to generate a repressive chromatin structure characterized by H3K9me3 in myeloid and erythroid differentiation . Also, PU.1 has been shown to act in concert with MITF to recruit co-repressors to osteoclast-specific in committed myeloid precursors capable of forming either macrophages or OC . Moreover, previous studies have shown that PU.1 can form a complex with DNMT3a and DNMT3b . However, this is the first report where the association between PU.1 and DNMTs in association with gene repression is shown in this context.
Moreover, our findings constitute the first report where the binding of PU.1 to TET2 has been described. Several recent reports have pointed at TET2-mediated hydroxylation of 5-methylcytosine as an intermediate step towards demethylation  and our data show changes in 5hmC at genes that become demethylated in osteoclastogenesis, reinforcing the possibility that PU.1-mediated recruitment of TET2 is leading to 5hmC-mediated demethylation. However the detailed mechanisms that couple hydroxylation of 5mC and demethylation are still objects of debate.
The manipulation of PU.1 levels by using siRNAs has shown that PU.1 has a direct role in recruiting DNMT3b and TET2 to its target promoters, as well as showing how impaired association of PU.1 results in defective acquisition of DNA methylation changes in both directions as well as reduced effect on gene expression changes. Therefore, our data reveal a novel role of PU.1 as a dual adaptor with the ability to bind both epigenetically repressive and epigenetically activating events and targeting DNA methylation changes in both directions (Figure 6). The incomplete impairment of DNA methylation and expression changes, as well as partial loss of Tet2 and DNMT3b following PU.1 knock-down indicates that additional transcription factors are also participating in this process. In future studies, it will also be interesting to identify the mechanisms that operate in the specific recruitment of PU.1-TET2 to genes that become demethylated, and to determine how PU.1-DNMT3b is recruited to genes that become hypermethylated. It is likely that specific transcription factors play a role, and specific post-translational modifications in PU.1 may participate in the coupling of its associated complexes to specific factors. In this context, Ser phosphorylation of PU.1 has already been shown to play a role in its recruitment to promoters  and could also participate in discriminating interaction with epigenetic modifiers.
Our study has allowed us to identify key DNA methylation changes during OC differentiation and has revealed an implication of PU.1 in the acquisition of DNA methylation and expression changes as well as identifying novel interactions with DNMT3b and TET2.
Our study of the DNA methylation changes in monocyte-to-osteoclast differentiation reveals the occurrence of both hypomethylation and hypermethylation changes. These changes occur in the virtual absence of DNA replication suggesting the participation of active mechanisms, particularly relevant for hypomethylation events, for which the mechanisms are still subject of debate. Also, when comparing the dynamics of DNA methylation and expression changes, hypomethylation occurs concomitant or even earlier than expression changes. In contrast, for the majority of genes becoming hypermethylated, hypermethylation follows expression changes. Hypomethylation takes place in relevant functional categories related with OC differentiation and most of the genes that are necessary for OC function undergo hypomethylation including ACP5, CTSK, and TM7SF4 among others. The analysis of over-representation of transcription factor binding motifs reveals the enrichment of specific motifs for hypomethylated and hypermethylated genes. Among these, PU.1 and other ETS-related binding motifs are highly enriched in both hypomethylated and hypermethylated genes. We have demonstrated that PU.1 is bound to both hypo- and hypermethylated promoters and that it is able to recruit both DNMT3b and TET2. Most importantly, downregulation of PU.1 with siRNAs not only shows a reduction in the recruitment of these two enzymes to PU.1 target genes but also results in a specific reduction in the acquisition of DNA methylation and expression changes at those targets. Our results demonstrate a key role of PU.1 in driving DNA methylation changes during OC differentiation.
Materials and methods
Differentiation of OCs from peripheral blood mononuclear cells
Human samples (blood) used in this study came from anonymous blood donors and were obtained from the Catalan Blood and Tissue Bank (Banc de Sang i Teixits) in Barcelona as thrombocyte concentrates (buffy coats). The anonymous blood donors received oral and written information about the possibility that their blood would be used for research purposes, and any questions that arose were then answered. Prior to obtaining the first blood sample the donors signed a consent form at the Banc de Teixits. The Banc de Teixits follows the principles set out in the WMA Declaration of Helsinki. The blood was carefully layered on a Ficoll-Paque gradient (Amersham, Buckinghamshire, UK) and centrifuged at 2,000 rpm for 30 min without braking. After centrifugation, peripheral blood mononuclear cells (PBMCs), in the interface between the plasma and the Ficoll-Paque gradient, were collected and washed twice with ice-cold PBS, followed by centrifugation at 2,000 rpm for 5 min. Pure CD14+ cells were isolated from PBMCs using positive selection with MACS magnetic CD14 antibody (Miltenyi Biotec). Cells were then resuspended in α-minimal essential medium (α-MEM, Glutamax no nucleosides) (Invitrogen, Carlsbad, CA, USA) containing 10% fetal bovine serum, 100 units/mL penicillin, 100 μg/mL streptomycin and antimycotic and supplemented with 25 ng/mL human M-CSF and 50 ng/mL hRANKL soluble (PeproTech EC, London, UK). Depending on the amount needed, cells were seeded at a density of 3 · 105 cells/well in 96-well plates, 5 · 106 cells/well in 6-well plates or 40 · 106 cells in 10 mm plates and cultured for 21 days (unless otherwise noted); medium and cytokines were changed twice a week. The presence of OCs was checked by tartrate-resistant acid phosphatase (TRAP) staining using the Leukocyte Acid Phosphatase Assay Kit (Sigma-Aldrich) according to the manufacturer’s instructions. A phalloidin/DAPI stain allowed us to confirm that the populations were highly enriched in multinuclear cells, some of them containing more than 40 nuclei. We used several methods to determine that on day 21 almost 85% of the nuclei detected were ‘osteoclastic nuclei’ (in polykaryons, nuclei and not cells were quantified). OCs (TRAP-positive cells with more than three nuclei) were also analyzed at the mRNA level: upregulation of key OC markers (TRAP/ACP5, CA2, MMP9, and CTSK) and the downregulation of the MO marker CX3CR1 were confirmed.
Treatment of MOs with 5-aza-2-deoxycytidine
In some cases we performed monocyte-to-osteoclast differentiation experiments in the presence of different subtoxic concentrations of the DNA replication-coupled demethylating drug 5-aza-2-deoxycytidine (at 50 nM, 500 mM) for 72 h.
Visualization of OCs with phalloidin and DAPI staining
PBMCs or pure isolated CD14+ cells were seeded and cultured in glass Lab-Tek Chamber Slides (Thermo Fisher Scientific) for 21 days in the presence of hM-CSF and hRANKL. OCs were then washed twice with PBS and fixed (3.7% paraformaldehyde, 15 min). Cells were permeabilized with 0.1% (V/V) Triton X-100 for 5 min and stained for F-actin with 5 U/mL Alexa Fluor® 647-Phalloidin (Invitrogen). Cells were then mounted in Mowiol-DAPI mounting medium. Cultures were visualized by CLSM (Leica TCP SP2 AOBS confocal microscope).
DNA methylation profiling using universal bead arrays
Infinium HumanMethylation450 BeadChips (Illumina, Inc.) were used to analyze DNA methylation. This array allows interrogating >485,000 methylation sites per sample at single-nucleotide resolution, covering 99% of RefSeq genes, with an average of 17 CpG sites per gene region distributed across the promoter, 5′UTR, first exon, gene body and 3′UTR. It covers 96% of CpG islands, with additional coverage in CpG island shores and the regions flanking them. DNA samples were bisulfite converted using the EZ DNA methylation kit (Zymo Research, Orange, CA, USA). After bisulfite treatment, the remaining assay steps were performed following the specifications and using the reagents supplied and recommended by the manufacturer. The array was hybridized using a temperature gradient program, and arrays were imaged using a BeadArray Reader (Illumina, Inc.). The image processing and intensity data extraction software and procedures were those previously described . Each methylation data point is obtained from a combination of the Cy3 and Cy5 fluorescent intensities from the M (methylated) and U (unmethylated) alleles. Background intensity computed from a set of negative controls was subtracted from each data point. For representation and further analysis we used both Beta values and M values. The Beta-value is the ratio of the methylated probe intensity and the overall intensity (sum of methylated and unmethylated probe intensities). The M value is calculated as the log2 ratio of the intensities of methylated probe versus unmethylated probe. The Beta value ranges from 0 to 1 and is more intuitive and was used in heatmaps and in comparisons with DNA methylation percentages from bisulfite pyrosequencing experiments, however for statistics purposes it is more adequate the use of M values .
Detection of differentially methylated CpGs
The approach to select differentially methylated CpGs was implemented in R , a well-known language in statistical computing. In order to process Illumina Infinium HumanMethylation450 methylation data, we used the methods supplied in limma , genefilter, and lumi  packages from Bioconductor repository. Previous to statistical analysis, a pre-process stage is applied, the main steps are: (1) color balance adjustment, that is, normalization between two color channels; (2) performing quantile normalization based on color balance adjusted data; and (3) variance filtering by IQR (interquartile range) using 0.50 for threshold value. Subsequently, for statistical analysis, eBayes moderated t-statistics test was carried out from limma package . Specifically, a paired limma was performed as designed in IMA package . To choose significant differences in methylated CpGs several criteria have been proposed. In this study, we considered a probe as differentially methylated if: (1) it has a fold-change >2 for hypermethylated and <0.5 hypomethylated; and (2) the statistical test was significant (P value <0.01 and FDR <0.05).
Identification of genomic clusters of differentially methylated CpGs
A clustering method was applied to the differenced methylated CpGs from charm package . We re-implemented the code to invoke the main clustering function using genomic CpG localization: identify differentially methylated regions (DMRs) by grouping differentially methylated probes (DMPs). The maximum allowable gap between probe positions for probes to be clustered into the same region was set to 500 bp. It has been shown that in many cases methylation changes are observed over a range of CpGs, which may be identified for instance at shores close to transcription starting sites. We considered that DMR are more robust signals than DMPs. In this analysis, the considered list of CpGs attains a P value <0.01 and FDR <0.05.
Bisulfite sequencing and pyrosequencing
We used bisulfite pyrosequencing to validate CpG methylation changes resulting from the analysis with the Infinium HumanMethylation450 BeadChips. Bisulfite modification of genomic DNA isolated from MOs, OCs, and samples from time course or PU.1-knockdown experiments was carried out as described by Herman et al. . A total of 2 μL of the converted DNA (corresponding to approximately 20–30 ng) were then used as a template in each subsequent PCR. Primers for PCR amplification and sequencing were designed with the PyroMark® Assay Design 2.0 software (Qiagen). PCRs were performed with the HotStart Taq DNA polymerase PCR kit (Qiagen), and the success of amplification was assessed by agarose gel electrophoresis. PCR products were pyrosequenced with the Pyromark™ Q24 system (Qiagen). In the case of repetitive elements (Sat2, D4Z4, NBL2, 18S rRNA, and 28 rRNA) we performed standard bisulfite sequencing of a minimum of 10 clones. Results from bisulfite pyrosequencing and sequencing of multiple clones are presented as a percentage of methylation. All primer sequences are listed in Additional file 11. Raw data for bisulfite sequencing of all samples is presented in Additional file 4.
Gene expression data analysis and comparison of DNA expression data versusDNA methylation data
In order to compare expression data versus methylation data, we used CD14+ and OC expression data from ArrayExpress database  under the accession name (E-MEXP-2019) from a previous publication . Affymetrix GeneChip Human Genome U133 Plus 2.0 expression data was processed using limma  and affy  packages from bioconductor. The preprocessing stage is divided into three major steps: (1) background correction; (2) normalization; and (3) reporter summarization. Here, the expresso function in the affy package was chosen for preprocessing. Thus, the RMA method  was applied for background correction. Then, a quantile normalization was performed. In addition, we introduced a specific step for PM (perfect match probes) adjustment, utilizing the PM-only model based expression index (option ‘pmonly’). And finally, for summarization step, the median polish method was taken. Next, as previously in the methylation analysis, a variance filtering by IQR using 0.50 for threshold value was executed. After preprocessing, a statistical analysis was applied, using eBayes moderated t-statistics test from limma package. Subsequently, expression genes matching to methylated genes were studied. Genes differentially expressed between MOs and Mo-OCs groups were selected with a criteria of P value <0.01 and FDR <0.05 as calculated by Benjamini-Hochberg and a fold-change of expression >2 or <0.5. Validation of expression data was performed by quantitative RT-PCR. All primer sequences are listed in Additional file 11.
Gene ontology analysis
Gene ontology (GO) was analyzed with the FatiGO tool , which uses Fisher’s exact test to detect significant over-representation of GO terms in one of the sets (list of selected genes) with respect to the other (the rest of the genome). Multiple test correction to account for the multiple hypotheses tested (one for each GO term) was applied to reduce false-positive results. GO terms with adjusted values of P <0.05 were considered significant.
Analysis of transcription factor binding
We used the STORM algorithm  to identify potential over-representation of transcription factor motifs in the 500-bp region around the center of the hypomethylated/hypermethylated CpG sites (as well as for all other CpGs-containing probes contained in the array) assuming cutoff values of P = 0.00002 (for hypo-/hypermethylated probes) and 0.00001 (for all other probes), using position frequency matrices (PFMs) from the TRANSFAC database (Professional version, release 2009.4) . Enrichment analysis of predicted TF in the probes of significant hypomethylated probes (n = 421) was conducted using GiTools [67, 68]. We calculated two-tailed probabilities, and a final adjusted FDR P value (with 0.25 cutoff) was considered statistically significant.
We downloaded PU.1 ChIPseq data for CD14+ MOs generated by Michael Rehli’s laboratory  from the Gene Expression Omnibus (GSE31621). The genomic locations of the calculated peaks were mapped to GRCh37.p10 human alignment obtained from Biomart , by using bedtools (intersect function) in order to obtain the PU.1 occupied genes. To determine whether a given CpG (from the Illumina bead array) was positive for PU.1 binding, we used the same 500-bp window used for TRANSFAC analysis.
Graphics and heatmaps
All graphs were created using Prism5 Graphpad. Heatmaps were generated from the expression or methylation data using the Genesis program (Graz University of Technology) .
BrdU proliferation assays
BrdU was used at a final concentration of 300 μm, as previously described. On the days specified, BrdU pulsing solution was added to each well for 2 to 4 days. For confocal microscopy of monocyte-to-osteoclast differentiation samples, CD14+ cells were seeded on Millicell EZ 8-well glass slides (Millipore) and cultured in differentiation media. At different times BrdU was added to the medium and after 2 to 4 days cells were fixed (4% paraformaldehyde, 30 min, RT), permeabilized (PBS-BSA-Triton X-100 0.8% (PBT), 10 min, RT) and treated with HCl 2 N for 30 min. After DNA opening, HCl was neutralized by two 5-min washes with NaBo (0.1 M, pH 8.5) and two 5-min washes with PBT. Cells were incubated with anti-BrdU antibody (18 h at 4°C, 1:1,000 dilution) and an anti-mouse Alexa-568 conjugated antibody was added to visualize the BrdU-positive nuclei. A phalloidin incubation step and Mowiol-DAPI mounting media were used.
Transfection of primary human MOs
We used two different Silencer® select pre-designed siRNAs against human PU.1 (one targeting exon 2 and another targeting the 3′UTR) and a Silencer® select negative control to perform PU.1 knockdown experiments in peripheral blood MOs. We used Lipofectamine RNAiMAX Transfection Reagent (Invitrogen) for efficient siRNA transfection. mRNA and protein levels were examined by quantitative RT-PCR and western blot at 1, 2, 4, and 6 days after siRNA transfection. In this case MO samples were prepared by incubating PBMCs in plates in α-MEM without serum for 30 min and washing out the unattached cells. Under these conditions over 80% are MOs. This alternative protocol was used for increased viability following transfection. These experiments were performed with three biological replicates.
Chromatin immunoprecipitation (ChIP) assays and immunoprecipitation experiments
Immunoprecipitation was performed by standard procedures in CD14+ cells at 0, 2, and 4 days after treatment with M-CSF and RANKL. Cell extracts were prepared in 50 mM Tris–HCl, pH 7.5, 1 mM EDTA, 150 mM NaCl, 1% Triton-X-100, and protease cocktail inhibitors (Complete, Roche Molecular Biochemicals). Cellular extracts and samples from immunoprecipitation experiments were electrophoresed and western blotted following standard procedures.
For chromatin immunoprecipitation (ChIP) assays, CD14+ at 0, 2, and 4 days after treatment with M-CSF and RANKL were cross-linked with 1% formaldehyde and subjected to immunoprecipitation after sonication. ChIP experiments were performed as described . Analysis was performed by real-time quantitative PCR. Data are represented as the ratio of the bound fraction over the input for each specific factor. We used a mouse monoclonal antibody against the TET2 N-t for ChIPs and a rabbit polyclonal antibody against TET2 for western blot. For DNMT3b we used a rabbit polyclonal against amino acids 1–230 of human DNMT3b (sc-20704, Santa Cruz Biotechnology). We also used a rabbit polyclonal against the C-t of PU.1 (sc-352, Santa Cruz Biotechnology), a rabbit polyclonal against the N-t of c-Fos (sc-52, Santa Cruz Biotechnology) and a rabbit polyclonal against the C-t of NF-kB p65 (sc-372, Santa Cruz Biotechnology). IgG was used as a negative control. Primer sequences were designed to contain either predicted or known TF binding (from TRANSFAC or ChIPseq data) as close as possible from the CpG undergoing methylation changes. Primer sequences are shown in Additional file 11. These experiments were performed with three biological replicates.
5hmC was analyzed using the Quest 5hmC Detection system (Zymo). Genomic DNA was treated with a specific 5hmC glucosyltransferase (GT) or left untreated (No GT, 0% 5hmC). DNA was then digested with MspI (100U) at 37°C overnight, followed by column purification. The MspI-resistant fraction (bearing the glucosile group, and therefore the original 5hmC) was quantitated by qPCR using primers designed around at least one MspI site (CCGG), and normalized to the amplification of the same region in the original DNA input. The amplification obtained in the untreated (no GT, MspI sensible) was then substracted to the samples in order to calculate the level of 0% 5hmC. The resulting values were the percentage of 5hmC present in each of the samples. Primer sequences are shown in Additional file 11.
Amplification of unmethylated Alus
This method, aiming at the amplification of unmethylated Alus (AUMA), was performed as described [31, 39]. Products were resolved on denaturing sequencing gels. Bands were visualized by silver staining the gels. AUMA fingerprints were visually checked for methylation differences between bands in different samples.
Amplification of unmethylated Alu repeats.
Pham TH, Benner C, Lichtinger M, Schwarzfischer L, Hu Y, Andreesen R, Chen W, Rehli M: Dynamic epigenetic enhancer signatures reveal key transcription factors associated with monocytic differentiation states. Blood. 2012, 119: e161-e171. 10.1182/blood-2012-01-402453.
Hogart A, Lichtenberg J, Ajay SS, Anderson S, Margulies EH, Bodine DM: Genome-wide DNA methylation profiles in hematopoietic stem and progenitor cells reveal overrepresentation of ETS transcription factor binding sites. Genome Res. 2012, 22: 1407-1418. 10.1101/gr.132878.111.
Zhang JA, Mortazavi A, Williams BA, Wold BJ, Rothenberg EV: Dynamic transformations of genome-wide epigenetic marking and transcriptional control establish T cell identity. Cell. 2012, 149: 467-482. 10.1016/j.cell.2012.01.056.
Blair HC, Teitelbaum SL, Ghiselli R, Gluck S: Osteoclastic bone resorption by a polarized vacuolar proton pump. Science. 1989, 245: 855-857. 10.1126/science.2528207.
Wiktor-Jedrzejczak W, Bartocci A, Ferrante AW, Ahmed-Ansari A, Sell KW, Pollard JW, Stanley ER: Total absence of colony-stimulating factor 1 in the macrophage-deficient osteopetrotic (op/op) mouse. Proc Natl Acad Sci USA. 1990, 87: 4828-4832. 10.1073/pnas.87.12.4828.
Lacey DL, Timms E, Tan HL, Kelley MJ, Dunstan CR, Burgess T, Elliott R, Colombero A, Elliott G, Scully S, et al: Osteoprotegerin ligand is a cytokine that regulates osteoclast differentiation and activation. Cell. 1998, 93: 165-176. 10.1016/S0092-8674(00)81569-X.
Saltel F, Chabadel A, Bonnelye E, Jurdic P: Actin cytoskeletal organisation in osteoclasts: a model to decipher transmigration and matrix degradation. Eur J Cell Biol. 2008, 87: 459-468. 10.1016/j.ejcb.2008.01.001.
Wong BR, Besser D, Kim N, Arron JR, Vologodskaia M, Hanafusa H, Choi Y: TRANCE, a TNF family member, activates Akt/PKB through a signaling complex involving TRAF6 and c-Src. Mol Cell. 1999, 4: 1041-1049. 10.1016/S1097-2765(00)80232-4.
Kobayashi N, Kadono Y, Naito A, Matsumoto K, Yamamoto T, Tanaka S, Inoue J: Segregation of TRAF6-mediated signaling pathways clarifies its role in osteoclastogenesis. Embo J. 2001, 20: 1271-1280. 10.1093/emboj/20.6.1271.
Blank U, Launay P, Benhamou M, Monteiro RC: Inhibitory ITAMs as novel regulators of immunity. Immunol Rev. 2009, 232: 59-71. 10.1111/j.1600-065X.2009.00832.x.
Humphrey MB, Ogasawara K, Yao W, Spusta SC, Daws MR, Lane NE, Lanier LL, Nakamura MC: The signaling adapter protein DAP12 regulates multinucleation during osteoclast development. J Bone Miner Res. 2004, 19: 224-234.
Koga T, Inui M, Inoue K, Kim S, Suematsu A, Kobayashi E, Iwata T, Ohnishi H, Matozaki T, Kodama T, et al: Costimulatory signals mediated by the ITAM motif cooperate with RANKL for bone homeostasis. Nature. 2004, 428: 758-763. 10.1038/nature02444.
Humphrey MB, Daws MR, Spusta SC, Niemi EC, Torchia JA, Lanier LL, Seaman WE, Nakamura MC: TREM2, a DAP12-associated receptor, regulates osteoclast differentiation and function. J Bone Miner Res. 2006, 21: 237-245.
Negishi-Koga T, Takayanagi H: Ca2 + -NFATc1 signaling is an essential axis of osteoclast differentiation. Immunol Rev. 2009, 231: 241-256. 10.1111/j.1600-065X.2009.00821.x.
Takayanagi H, Kim S, Koga T, Nishina H, Isshiki M, Yoshida H, Saiura A, Isobe M, Yokochi T, Inoue J, Wagner EF, Mak TW, Kodama T, Taniguchi T: Induction and activation of the transcription factor NFATc1 (NFAT2) integrate RANKL signaling in terminal differentiation of osteoclasts. Dev Cell. 2002, 3: 889-901. 10.1016/S1534-5807(02)00369-6.
Sharma SM, Bronisz A, Hu R, Patel K, Mansky KC, Sif S, Ostrowski MC: MITF and PU.1 recruit p38 MAPK and NFATc1 to target genes during osteoclast differentiation. J Biol Chem. 2007, 282: 15921-15929. 10.1074/jbc.M609723200.
Kim K, Lee SH, Ha Kim J, Choi Y, Kim N: NFATc1 induces osteoclast fusion via up-regulation of Atp6v0d2 and the dendritic cell-specific transmembrane protein (DC-STAMP). Mol Endocrinol. 2008, 22: 176-185.
Yu M, Moreno JL, Stains JP, Keegan AD: Complex regulation of tartrate-resistant acid phosphatase (TRAP) expression by interleukin 4 (IL-4): IL-4 indirectly suppresses receptor activator of NF-kappaB ligand (RANKL)-mediated TRAP expression but modestly induces its expression directly. J Biol Chem. 2009, 284: 32968-32979. 10.1074/jbc.M109.001016.
Matsumoto M, Kogawa M, Wada S, Takayanagi H, Tsujimoto M, Katayama S, Hisatake K, Nogi Y: Essential role of p38 mitogen-activated protein kinase in cathepsin K gene expression during osteoclastogenesis through association of NFATc1 and PU.1. J Biol Chem. 2004, 279: 45969-45979. 10.1074/jbc.M408795200.
Sundaram K, Nishimura R, Senn J, Youssef RF, London SD, Reddy SV: RANK ligand signaling modulates the matrix metalloproteinase-9 gene expression during osteoclast differentiation. Exp Cell Res. 2007, 313: 168-178. 10.1016/j.yexcr.2006.10.001.
Tolar J, Teitelbaum SL, Orchard PJ: Osteopetrosis. N Engl J Med. 2004, 351: 2839-2849. 10.1056/NEJMra040952.
Rachner TD, Khosla S, Hofbauer LC: Osteoporosis: now and the future. Lancet. 2011, 377: 1276-1287. 10.1016/S0140-6736(10)62349-5.
Scott DL, Wolfe F, Huizinga TW: Rheumatoid arthritis. Lancet. 2010, 376: 1094-1108. 10.1016/S0140-6736(10)60826-4.
Mundy GR, Raisz LG, Cooper RA, Schechter GP, Salmon SE: Evidence for the secretion of an osteoclast stimulating factor in myeloma. N Engl J Med. 1974, 291: 1041-1046. 10.1056/NEJM197411142912001.
Yoneda T: Cellular and molecular mechanisms of breast and prostate cancer metastasis to bone. Eur J Cancer. 1998, 34: 240-245. 10.1016/S0959-8049(97)10132-0.
Mii Y, Miyauchi Y, Morishita T, Miura S, Honoki K, Aoki M, Tamai S: Osteoclast origin of giant cells in giant cell tumors of bone: ultrastructural and cytochemical study of six cases. Ultrastruct Pathol. 1991, 15: 623-629. 10.3109/01913129109023192.
Nicholson GC, Malakellis M, Collier FM, Cameron PU, Holloway WR, Gough TJ, Gregorio-King C, Kirkland MA, Myers DE: Induction of osteoclasts from CD14-positive human peripheral blood mononuclear cells by receptor activator of nuclear factor kappaB ligand (RANKL). Clin Sci (Lond). 2000, 99: 133-140. 10.1042/CS19990355.
Sorensen MG, Henriksen K, Schaller S, Henriksen DB, Nielsen FC, Dziegiel MH, Karsdal MA: Characterization of osteoclasts derived from CD14+ monocytes isolated from peripheral blood. J Bone Miner Metab. 2007, 25: 36-45.
Youn MY, Takada I, Imai Y, Yasuda H, Kato S: Transcriptionally active nuclei are selective in mature multinucleated osteoclasts. Genes Cells. 2010, 15: 1025-1035. 10.1111/j.1365-2443.2010.01441.x.
Saltman LH, Javed A, Ribadeneyra J, Hussain S, Young DW, Osdoby P, Amcheslavsky A, van Wijnen AJ, Stein JL, Stein GS, Lian JB, Bar-Shavit Z: Organization of transcriptional regulatory machinery in osteoclast nuclei: compartmentalization of Runx1. J Cell Physiol. 2005, 204: 871-880. 10.1002/jcp.20329.
Rodriguez J, Vives L, Jorda M, Morales C, Munoz M, Vendrell E, Peinado MA: Genome-wide tracking of unmethylated DNA Alu repeats in normal and cancer cells. Nucleic Acids Res. 2008, 36: 770-784.
Gallois A, Lachuer J, Yvert G, Wierinckx A, Brunet F, Rabourdin-Combe C, Delprat C, Jurdic P, Mazzorana M: Genome-wide expression analyses establish dendritic cells as a new osteoclast precursor able to generate bone-resorbing cells more efficiently than monocytes. J Bone Miner Res. 2010, 25: 661-672. 10.1359/jbmr.090829.
Momparler RL: Pharmacology of 5-Aza-2′-deoxycytidine (decitabine). Semin Hematol. 2005, 42: S9-S16.
Tahiliani M, Koh KP, Shen Y, Pastor WA, Bandukwala H, Brudno Y, Agarwal S, Iyer LM, Liu DR, Aravind L, Rao A: Conversion of 5-methylcytosine to 5-hydroxymethylcytosine in mammalian DNA by MLL partner TET1. Science. 2009, 324: 930-935. 10.1126/science.1170116.
Kriaucionis S, Heintz N: The nuclear DNA base 5-hydroxymethylcytosine is present in Purkinje neurons and the brain. Science. 2009, 324: 929-930. 10.1126/science.1169786.
Ito S, Shen L, Dai Q, Wu SC, Collins LB, Swenberg JA, He C, Zhang Y: Tet proteins can convert 5-methylcytosine to 5-formylcytosine and 5-carboxylcytosine. Science. 2011, 333: 1300-1303. 10.1126/science.1210597.
Wu H, Zhang Y: Mechanisms and functions of Tet protein-mediated 5-methylcytosine oxidation. Genes Dev. 2011, 25: 2436-2452. 10.1101/gad.179184.111.
Williams K, Christensen J, Helin K: DNA methylation: TET proteins-guardians of CpG islands?. EMBO Rep. 2012, 13: 28-35.
Hernando H, Shannon-Lowe C, Islam AB, Al-Shahrour F, Rodriguez-Ubreva J, Rodriguez-Cortez VC, Javierre BM, Mangas C, Fernandez AF, Parra M, Delecluse HJ, Esteller M, Lopez-Granados E, Fraga MF, Lopez-Bigas N, Ballestar E: The B cell transcription program mediates hypomethylation and overexpression of key genes in Epstein-Barr virus-associated proliferative conversion. Genome Biol. 2013, 14: R3-10.1186/gb-2013-14-1-r3.
Stadler MB, Murr R, Burger L, Ivanek R, Lienert F, Scholer A, van Nimwegen E, Wirbelauer C, Oakeley EJ, Gaidatzis D, Tiwari VK, Schubeler D: DNA-binding factors shape the mouse methylome at distal regulatory regions. Nature. 2012, 480: 490-495.
Edwards JR, Mundy GR: Advances in osteoclast biology: old findings and new insights from mouse models. Nat Rev Rheumatol. 2011, 7: 235-243. 10.1038/nrrheum.2011.23.
Shen L, Wu H, Diep D, Yamaguchi S, D’Alessio AC, Fung HL, Zhang K, Zhang Y: Genome-wide analysis reveals TET- and TDG-dependent 5-methylcytosine oxidation dynamics. Cell. 2013, 153: 692-706. 10.1016/j.cell.2013.04.002.
Cortellino S, Xu J, Sannai M, Moore R, Caretti E, Cigliano A, Le Coz M, Devarajan K, Wessels A, Soprano D, Abramowitz LK, Bartolomei MS, Rambow F, Bassi MR, Bruno T, Fanciulli M, Renner C, Klein-Szanto AJ, Matsumoto Y, Kobi D, Davidson I, Alberti C, Larue L, Bellacosa A: Thymine DNA glycosylase is essential for active DNA demethylation by linked deamination-base excision repair. Cell. 2011, 146: 67-79. 10.1016/j.cell.2011.06.020.
Kallin EM, Rodriguez-Ubreva J, Christensen J, Cimmino L, Aifantis I, Helin K, Ballestar E, Graf T: Tet2 facilitates the derepression of myeloid target genes during CEBPalpha-induced transdifferentiation of pre-B cells. Mol Cell. 2012, 48: 266-276. 10.1016/j.molcel.2012.08.007.
Klug M, Schmidhofer S, Gebhard C, Andreesen R, Rehli M: 5-Hydroxymethylcytosine is an essential intermediate of active DNA demethylation processes in primary human monocytes. Genome Biol. 2013, 14: R46-10.1186/gb-2013-14-5-r46.
Suzuki M, Yamada T, Kihara-Negishi F, Sakurai T, Hara E, Tenen DG, Hozumi N, Oikawa T: Site-specific DNA methylation by a complex of PU.1 and Dnmt3a/b. Oncogene. 2006, 25: 2477-2488. 10.1038/sj.onc.1209272.
Chen W, Zhu G, Hao L, Wu M, Ci H, Li YP: C/EBPalpha regulates osteoclast lineage commitment. Proc Natl Acad Sci USA. 2013, 110: 7294-7299. 10.1073/pnas.1211383110.
Hama M, Kirino Y, Takeno M, Takase K, Miyazaki T, Yoshimi R, Ueda A, Itoh-Nakadai A, Muto A, Igarashi K, Ishigatsubo Y: Bach1 regulates osteoclastogenesis in a mouse model via both heme oxygenase 1-dependent and heme oxygenase 1-independent pathways. Arthritis Rheum. 2012, 64: 1518-1528. 10.1002/art.33497.
Stopka T, Amanatullah DF, Papetti M, Skoultchi AI: PU.1 inhibits the erythroid program by binding to GATA-1 on DNA and creating a repressive chromatin structure. EMBO J. 2005, 24: 3712-3723. 10.1038/sj.emboj.7600834.
Hu R, Sharma SM, Bronisz A, Srinivasan R, Sankar U, Ostrowski MC: Eos, MITF, and PU.1 recruit corepressors to osteoclast-specific genes in committed myeloid progenitors. Mol Cell Biol. 2007, 27: 4018-4027. 10.1128/MCB.01839-06.
Vincent JJ, Huang Y, Chen PY, Feng S, Calvopina JH, Nee K, Lee SA, Le T, Yoon AJ, Faull K, Fan G, Rao A, Jacobsen SE, Pellegrini M, Clark AT: Stage-specific roles for tet1 and tet2 in DNA demethylation in primordial germ cells. Cell Stem Cell. 2013, 12: 470-478. 10.1016/j.stem.2013.01.016.
Azim AC, Wang X, Park GY, Sadikot RT, Cao H, Mathew B, Atchison M, van Breemen RB, Joo M, Christman JW: NF-kappaB-inducing kinase regulates cyclooxygenase 2 gene expression in macrophages by phosphorylation of PU.1. J Immunol. 2007, 179: 7868-7875.
Bibikova M, Lin Z, Zhou L, Chudin E, Garcia EW, Wu B, Doucet D, Thomas NJ, Wang Y, Vollmer E, Goldmann T, Seifart C, Jiang W, Barker DL, Chee MS, Floros J, Fan JB: High-throughput DNA methylation profiling using universal bead arrays. Genome Res. 2006, 16: 383-393. 10.1101/gr.4410706.
Du P, Zhang X, Huang CC, Jafari N, Kibbe WA, Hou L, Lin SM: Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinforma. 2010, 11: 587-10.1186/1471-2105-11-587.
The-R-Development-Core-Team: R: A language and environment for statistical computing. 2012, R Foundation for Statistical Computing
Smyth GK: Limma: linear models for microarray data. ‘Bioinformatics and Computational Biology Solutions using R and Bioconductor’. Edited by: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W. 2005, New York: Springer, 397-420.
Du P, Kibbe WA, Lin SM: Lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008, 24: 1547-1548. 10.1093/bioinformatics/btn224.
Wang D, Yan L, Hu Q, Sucheston LE, Higgins MJ, Ambrosone CB, Johnson CS, Smiraglia DJ, Liu S: IMA: an R package for high-throughput analysis of Illumina’s 450 K Infinium methylation data. Bioinformatics. 2012, 28: 729-730. 10.1093/bioinformatics/bts013.
Aryee MJ, Wu Z, Ladd-Acosta C, Herb B, Feinberg AP, Yegnasubramanian S, Irizarry RA: Accurate genome-scale percentage DNA methylation estimates from microarray data. Biostatistics. 2011, 12: 197-210. 10.1093/biostatistics/kxq055.
Herman JG, Graff JR, Myohanen S, Nelkin BD, Baylin SB: Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands. Proc Natl Acad Sci USA. 1996, 93: 9821-9826. 10.1073/pnas.93.18.9821.
Gautier L, Cope L, Irizarry RA, Bolstad BM: Affy-analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20: 307-315. 10.1093/bioinformatics/btg405.
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31: e15-10.1093/nar/gng015.
Al-Shahrour F, Diaz-Uriarte R, Dopazo J: FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics. 2004, 20: 578-580. 10.1093/bioinformatics/btg455.
Schones DE, Smith AD, Zhang MQ: Statistical significance of cis-regulatory modules. BMC Bioinforma. 2007, 8: 19-10.1186/1471-2105-8-19.
Matys V, Fricke E, Geffers R, Gossling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel AE, Kel-Margoulis OV, Kloos DU, Land S, Lewicki-Potapov B, Michael H, Munch R, Reuter I, Rotert S, Saxel H, Scheer M, Thiele S, Wingender E: TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res. 2003, 31: 374-378. 10.1093/nar/gkg108.
Perez-Llamas C, Lopez-Bigas N: Gitools: analysis and visualisation of genomic data using interactive heat-maps. PLoS One. 2011, 6: e19541-10.1371/journal.pone.0019541.
Flicek P, Amode MR, Barrell D, Beal K, Brent S, Carvalho-Silva D, Clapham P, Coates G, Fairley S, Fitzgerald S, Gil L, Hendrix M, Hourlier T, Johnson N, Kahari AK, Keefe D, Keenan S, Kinsella R, Komorowska M, Koscielny G, Kulesha E, Larsson P, Longden I, McLaren W, Muffato M, Overduin B, Pignatelli M, Pritchard B, et al: Ensembl 2012. Nucleic Acids Res. 2012, 40: D84-D90. 10.1093/nar/gkr991.
Smink JJ, Begay V, Schoenmaker T, Sterneck E, de Vries TJ, Leutz A: Transcription factor C/EBPbeta isoform ratio regulates osteoclastogenesis through MafB. EMBO J. 2009, 28: 1769-1781. 10.1038/emboj.2009.127.
We thank Dr Fátima Al-Shahrour and Dr Núria López-Bigas for helpful suggestions about the bioinformatic analyses, and Dr Mercedes Garayoa and Dr Antonio Garcia-Gomez for their helpful suggestions about protocols and providing antibodies. This work was supported by grant SAF2011-29635 from the Spanish Ministry of Science and Innovation, grant CIVP16A1834 from Fundación Ramón Areces and grant 2009SGR184 from AGAUR (Catalan Government). LR is supported by a PFIS predoctoral fellowship.
Methylation array data for this publication has been deposited in NCBI’s Gene Expression Omnibus and is accessible through GEO Series accession number GSE46648.
The authors declare that they have no competing interests.
LR and EB conceived experiments; LR, MG, JR-U, and HH performed experiments; AI and JMU performed biocomputing analysis; LR, JMU, JC, KH, CGV, and EB analyzed the data; EB wrote the paper. All authors read and approved the final manuscript.