- Research article
- Open Access
DNA methylation analysis identifies key transcription factors involved in mesenchymal stem cell osteogenic differentiation
Biological Research volume 56, Article number: 9 (2023)
Knowledge about regulating transcription factors (TFs) for osteoblastogenesis from mesenchymal stem cells (MSCs) is limited. Therefore, we investigated the relationship between genomic regions subject to DNA-methylation changes during osteoblastogenesis and the TFs known to directly interact with these regulatory regions.
The genome-wide DNA-methylation signature of MSCs differentiated to osteoblasts and adipocytes was determined using the Illumina HumanMethylation450 BeadChip array. During adipogenesis no CpGs passed our test for significant methylation changes. Oppositely, during osteoblastogenesis we identified 2462 differently significantly methylated CpGs (adj. p < 0.05). These resided outside of CpGs islands and were significantly enriched in enhancer regions. We confirmed the correlation between DNA-methylation and gene expression. Accordingly, we developed a bioinformatic tool to analyse differentially methylated regions and the TFs interacting with them. By overlaying our osteoblastogenesis differentially methylated regions with ENCODE TF ChIP-seq data we obtained a set of candidate TFs associated to DNA-methylation changes. Among them, ZEB1 TF was highly related with DNA-methylation. Using RNA interference, we confirmed that ZEB1, and ZEB2, played a key role in adipogenesis and osteoblastogenesis processes. For clinical relevance, ZEB1 mRNA expression in human bone samples was evaluated. This expression positively correlated with weight, body mass index, and PPARγ expression.
In this work we describe an osteoblastogenesis-associated DNA-methylation profile and, using these data, validate a novel computational tool to identify key TFs associated to age-related disease processes. By means of this tool we identified and confirmed ZEB TFs as mediators involved in the MSCs differentiation to osteoblasts and adipocytes, and obesity-related bone adiposity.
Musculoskeletal disorders (MSD) are a set of pathologies that affect the locomotor system (bones, joints, peri-articular structures, and muscles). Research interest in these diseases has increased due to their elevated prevalence because of an ageing demographic resulting in a burgeoning economic cost to national health systems [1, 2]. In fact, 20% of Europeans undergo treatment or take supplements to help them cope with their MSD . Besides ageing, the increase in the prevalence of major MSD conditions is also associated with an increase in obesity [4, 5] and a sedentary lifestyle [6, 7]. A common link amongst major MSDs are alterations to bone, such as osteoporosis and osteopenia, which have been associated with an increase in patient fragility  and in turn to a reduced life expectancy . Bone alterations are also present in other MSDs, such as osteoarthritis where osteoporotic and sclerotic bone regions co-exist and both are linked to disease progression . Likewise, in rheumatoid arthritis, as well as in other inflammatory arthropathies, there are local and systemic bone alterations .
Approximately 15% of the human skeleton is renewed every year . This means that maintenance of an optimal bone mass depends on the precise balance between bone resorption and bone formation. When this balance is altered pathological conditions such as osteoporosis or osteopetrosis occur. Bone formation is carried out by osteoblasts, which share with adipocytes, their mesenchymal origin . Therefore, the process of bone formation in turn depends on the balance between the differentiation of mesenchymal stem cells (MSCs) towards osteoblast or adipocyte cell fates . It is considered that both processes, osteoblastogenesis and adipogenesis, are competing and reciprocal [12, 13]. In fact, bone marrow adiposity has been associated with bone loss in several MSDs as well as in other associated pathologies or conditions such as obesity and ageing [12, 14].
The commitment of the MSC towards osteoblast or adipocyte cell fate is controlled by certain transcription factors (TFs) such as runt related TF 2 (Runx2) and SP7 for osteoblastogenesis, and peroxisome proliferator activated receptor gamma (PPARγ) and C/EBPs  for adipogenesis. These TFs integrate the cell environment and the signalling of diverse pathways helping to initiate or block the differentiation process . This can potentially have clinical implications for example, long-term pharmacological activation of PPARγ, which promotes adipogenesis and inhibits osteoblastogenesis, increases fracture rates among patients with diabetes . Advances have been made investigating the role of TFs and their associated signalling networks in osteoblastogenesis and adipogenesis, however, despite the well-established transcriptional cascade associated with adipocyte differentiation  a better understanding of the TF network involved in osteoblast differentiation is required.
TF activity can be inhibited by DNA methylation through the blockade of their interaction with the DNA . This inhibition has been associated with repression of gene expression, which supports the key role of DNA methylation in multiple processes, including development and tumorigenesis [17, 18]. DNA methylation is considered a heritable repressive epigenetic mark that consists of the covalent addition of a methyl group in the 5′ position of a cytosine in a CpG dinucleotide . The abundance of these CpGs in the genome is not homogeneous, with more CpG-rich regions being present at TF-binding sites . As a result DNA methylation has been proposed as a major regulator of TF activity .
There is limited information about the TFs that regulate the osteoblastogenesis process. Therefore, by means of a combined genome-wide methylation analysis and bioinformatic approach and using the link between DNA methylation and TF activity, we identified the TFs potentially affected by the changes in the DNA methylation during osteoblastogenesis. Our analysis suggested the DNA binding sites for the TF ZEB1/ZEB2 (Zinc Finger E-Box Binding Homeobox 1/2) were hyper-methylated due to the process of osteoblast differentiation. Functional data confirmed the role of these TFs on both osteoblastogenesis and adipogenesis processes, emphasising the relevance of DNA methylation on the activity of key TFs during cellular differentiation. Consistent with this, ZEB1 expression in human pathological bone samples revealed a potential link between this TF and metabolic-mediated bone alterations.
Material and methods
IGF-1 was purchased from Peprotech (Rocky Hill, NJ, USA). DMEM, Foetal bovine serum, β-Glycerol Phosphate, dexamethasone, ascorbic acid 2-phosphate, insulin, 3-Isobutyl-1-methylxanthine, indomethacin, rosiglitazone, Cetylpyridinium, Oil Red, Alizarin Red were purchased from Sigma-Aldrich (St. Louise, MO, USA). Other products were also purchased from Sigma-Aldrich unless otherwise indicated.
MSCs were purchased from Lonza (BAS, Switzerland). All the cells came from women donors, MSCs were cultured, characterized, and their trilineage potential was determined as previously described [19, 20]. The procedure for their culture and differentiation to osteoblast or adipocytes was performed as we previously described .
Cells undergoing osteoblast differentiation were fixed in 70% cold ethanol (5 min, − 20 ºC). After drying the wells, to reveal calcium-rich mineralisation deposits the cells were incubated at room temperature with a solution of Alizarin Red (40 mM, pH 4.2) for 20–30 min. Prior to acquiring the images, the cells were gently washed with distilled water to avoid unspecific staining. For quantitation, the staining was eluted with 10% (w/v) Cetylpyridinium solubilized in 10 mM sodium phosphate buffer (pH 7.0), and the absorbance measured at 570 nm.
Cells undergoing adipogenesis were fixed with formalin for 1 h. After washing the wells with distilled water and 60% isopropanol the wells were dried. To reveal the presence of lipid droplets the cells were stained with a 21% (w/v) solution of Oil Red O for 10 min. Prior to acquiring the images, the wells were gently washed with distilled water to avoid unspecific staining. Staining images were quantified using the image analysis software ImageJ.
Bone samples from patients were obtained after total knee/hip replacement surgery for osteoarthritic and osteoporotic conditions. Healthy bone was from cadavers. Both healthy and pathological bone samples were comparable in terms of age and sex. The Ethics Committee for Research at Santiago-Lugo Area approved the protocol. Informed consent was obtained from all patients or patients’ families. Clinical data regarding weight and height was obtained from the clinical records, where available. BMI was calculated as weight (kg)/height (m)2.
To isolate bone RNA, bone explants were obtained assuring only trabecular bone was processed, thus without fat, cartilage, nor other fibroblastic or stromal tissues. Bone was repeatedly washed with Phosphate Buffered Saline until clean, frozen to − 80 ºC, and pulverised using a CellCrusher (Cellcrusher, Co. Cork, Ireland) following the manufacturer instructions. Per each 500 µl of bone powder, 1 ml of TriReagent was added to perform RNA extraction.
RNA extraction and real-time reverse transcription PCR
Cell cultures were disrupted, and RNA extraction was performed using Qiagen RNeasy mini kit (Qiagen, Crawley, UK) following manufacturer instructions. Alternatively, when the experiments were performed in 96-well plates cell cultures were disrupted in Ambion Cells-to-cDNA II Cell Lysis buffer (Life Technologies, Carlsbad, CA, USA). Total RNA was then extracted and converted to cDNA using M-MLV reverse transcriptase (Invitrogen, Waltham, MA, USA) and RT-PCR was performed using TaqMan® probes. Changes in gene expression levels were calculated as described previously .
DNA methylation and RNA expression arrays
Global DNA methylation analysis was performed using the Illumina Infinium HumanMethylation450K BeadChip array using DNA from three donor samples for each condition analysed (undifferentiated MSCs, MSCs differentiated to adipocytes, and MSCs differentiated to osteoblasts). All the samples were derived from the same donor to reduce variability. DNA, isolated using DNeasy Blood & Tissue Kit (Qiagen GmbH, Hilden, Germany), was bisulphite converted using the EpiTect® 96 Bisulphite Kit (Qiagen GmbH, Hilden, Germany), and 200 ng of bisulfite converted DNA was analysed using the array by the service provider Edinburgh Clinical Research Facility. The raw data were extracted using GenomeStudio (Illumina, San Diego, CA, USA) which provides the methylation data as β values: β = M/(M + U), where M and U represent the fluorescent signal of the methylation and unmethylated probes respectively. β values range from 0 (no methylation) to 1 (100% methylation).
For gene expression, an Illumina HumanHT-12 v4 Expression BeadChip array was performed using 200 ng of RNA [with an RNA integrity score > 7 (Agilent bioanalyzer 2100)] for each sample. Three samples from different donors and for each condition (undifferentiated MSCs, and MSCs differentiated to osteoblasts) were analysed. The RNA samples were processed according to the manufacturer’s protocol and the array was performed by Central Biotechnology Services, Cardiff University. Gene expression data were analysed essentially as previously described . All data are available on request.
HumanMethylation450K BeadChip data were normalised and analysed using the R language and the Tost analysis pipeline . Further analyses to identify differentially methylated CpGs during differentiation were performed essentially using our own scripts developed from the limma package . Data were transformed from β values to the more statistically valid M-values . To explore the relationship between the distribution of differentially methylated CpGs and TF-binding and thus activity during the differentiation process we created an R-script, Regulatory Element Interrogation Script (REINS), available as an R Markdown document (Additional file 1). This script allows the download, management and overlap of the information provided by the ENCODE database , essentially ChIP-seq information about the genome-wide binding sites for 161 TFs in 91 different cell types, and differently methylated CpGs from any source including HumanMethylation450K BeadChip arrays. The script normalizes as a percentage the number of CpGs associated to a TF with the total number of CpGs. This allows the comparison of the overlap among different TFs with both hypo- or hyper-methylated CpGs. This normalization, “TF Relevance” (TFR), was performed using the following equations:
Finally, for the same TF, an unbalance in its ratio between TFRs for the hypo- and hyper-methylation (RRT, Relative Relevance of a TF) was calculated.
A computational model of BMP2 signalling was constructed to evaluate ALPL induction in the context of ZEB TFs activity. The design and construction of this model is described in Additional file 2.
RNA-mediated interference (ZEBs expression inhibition)
For siRNA transfection 50 nM siRNA was transfected into 50% confluent MSCs using Dharmafect™ 1 lipid reagent (Thermo Fisher, Waltham, MA, USA). MSC were plated at a density of 5000 cells/well in 96-well plates for 24 h then media were replaced with osteogenic or adipogenic medium. After day 7 of differentiation the expression of key marker genes of differentiation was evaluated by RT-PCR. Cell cultures were cultured for 14 or 21 days to evaluate lipid accumulation or mineral deposition, respectively. Dharmacon siRNA SMARTpools® (Thermo Fisher, Waltham, MA, USA) of 4 specific siRNA duplexes (total of 50 nM siRNA) were used to target ZEB1 and ZEB2 TFs. Depletion of gene-specific mRNA levels was calculated by comparison of expression levels with cells transfected with 50 nM siCONTROL (non-targeting siRNA 2, cat. 001210–02; Dharmacon, Lafayette, CO, USA).
Data are expressed as mean ± standard error of the mean (SEM) for at least 3 independent experiments. Statistical differences were determined using a one-way analysis of variance or Kruskal–Wallis test, followed by a Bonferoni or Dunn’s post-hoc test, respectively, or Student’s t test or Mann–Whitney when appropriate. Contingency table statistical analysis was performed using Fisher exact tests. Correlations were assessed with Spearman tests. Statistical analysis was performed using Prism software (GraphPad Software Inc), p < 0.05 was considered significant. The statistical analysis of the array was performed using R software and Bioconductor R-packages .
Phenotypical characterization of MSC differentiation to osteoblasts and adipocytes
We examined the DNA methylation profile associated with the differentiation of human MSC to osteoblasts and adipocytes.
MSCs were differentiated in adipogenic or osteoblastogenic media for 14 or 21 days, respectively. Differentiation was confirmed by cytological staining with alizarin red for osteoblast and oil red for adipocytes (Fig. 1A). To further confirm the differentiation well-characterised osteoblast and adipocyte differentiation marker genes were measured by RT-PCR with expression consistent with the expected differentiation status and staining (Fig. 1B and C).
Global DNA methylation of MSC differentiation
DNA was extracted from MSC and differentiated cells and their global CpG methylation profile was assayed. For osteoblastogenesis, 2462 CpGs (1984 hypo-methylated and 478 hyper-methylated) showed a significant (adjusted p ≤ 0.05) change in their methylation (Fig. 2A) (Additional file 3). The overlapping genes of these CpGs (Additional file 3) included many in osteoblast metabolism such as RUNX2, GPNMB, CTSK, WNT5A, COL1A1 and PTH1R. However, surprisingly, only significant changes in CpG methylation were observed for osteoblastogenic differentiation (Fig. 2B), which involves a limited role for DNA methylation in establishing and maintaining the adipogenic phenotype (Additional file 4) (Fig. 2B).
Genomic topology of the differentially methylated CpGs during osteoblastogenesis
We examined the distribution of the 2462 differentially methylated CpGs within the genome and in relation to CpG island features. We found that osteoblastogenesis-dependent methylation remodelling was significantly associated with “open sea regions”. In contrast, methylation remodelling was less frequent than expected at CpG islands and their shores (Fig. 2C). Interestingly, only at “open sea regions” were the proportion of hypo-methylated CpGs higher than those hyper-methylated (Fig. 2C).
Chromatin profiling is a powerful tool to detect regulatory features within the underlying DNA . Therefore, we overlapped the location of the osteoblastogenesis differentially hypo- or hyper-methylated CpGs with the location of the fifteen chromatin states used to segment the genome by ENCODE for the mesoderm cell-line GM12878 , normalising for CpG location frequency bias from the array. A significantly higher frequency of differentially methylated CpGs were located at enhancer, active transcription, and heterochromatin chromatin states (Additional file 6). CpGs defined to be located within insulator, poised promoter or polycomb-repressed states were more frequently hyper-methylated, while the reciprocal was evident at transcription elongation state where CpGs were more frequently hypo-methylated (Fig. 2D).
Next, we determined the methylation levels of all CpGs for each chromatin state in either the initial MSCs or following differentiation into osteoblasts. As predicted, promoter regions were very hypo-methylated while regions associated with transcription elongation or heterochromatin were hyper-methylated. When comparing CpG methylation levels for the two cell states, MSC and osteoblasts, data revealed that the methylation of the CpGs that overlapped with the insulator state and with the heterochromatin state were significantly higher in the osteoblasts than MSCs whilst those associated with chromatin states for transcription transition and elongation, and weak enhancer activity were lower (cut off 0.1% value; p < 0.001) (Fig. 2E).
Correlation of DNA methylation and gene expression
To correlate DNA methylation with gene expression, we performed whole genome transcriptome analysis of RNA from the donor MSCs and differentiated osteoblasts. This analysis identified 2315 genes/transcripts as significantly differentially expressed (adjusted p < 0.05) during the differentiation process (Fig. 3A) (Additional file 5). The analysis of this gene list using Ingenuity Pathway Analysis (IPA) identified that the major bone related pathways were significantly enriched (osteoblast differentiation, WNT pathway activation) (Table 1), (Additional file 7). Interestingly, there was also a significant enrichment in genes associated with adipogenesis inhibition (Table 1), supporting the fact that osteoblastogenesis and adipogenesis are opposing cell fates.
Next, we explored whether gene expression and DNA methylation levels correlated in both MSCs and differentiated osteoblast. Using the methylation level of the significantly differentially methylated CpGs we observed a significant negative correlation only in osteoblasts (MSCs: r = − 0.0242, p = 0.6055 vs Osteoblasts: r = -− 0.377, p = 0.0031). Furthermore, when focussing on the methylation level of significantly differentially methylated CpGs that occur in 5’UTR regions, the negative correlation with gene expression improved for osteoblasts (r = − 0.2500; p = 0.0499), in line with data for other cell-types . This correlation further improved when analysing genes above a robust expression level threshold (> 1/3 of the maximum expression) (r = − 0.6289; p = 0.0013) (Fig. 3B).
Osteoblastogenesis DNA methylation changes and POL II-chromatin interaction loci
DNA methylation changes do not correlate with the gene expression and biochemical changes that occur during adipogenesis. Similarly, though DNA methylation did associate with osteoblastogenesis the correlation with gene expression was low, except for the described 5’UTR gene region CpGs and genes with a robust expression level. Using Chromatin Interaction Analysis by Paired-End Tag Sequencing (ChIA-PET), Li et al.  defined regions of the genome, essentially enhancers, topologically-associated with RNA polymerase II (RNAPII). Following the same approach, we determined, using ENCODE data for the cell-line MCF7, if the genomic loci of our significantly differentially methylated CpGs during osteoblastogenesis were RNAPII bound regions -indicative of active transcription. We observed that there was a significantly different (p = 0.0436) association of hypo- versus hyper-methylated CpGs with RNAPII regions, with a greater overlap for hypo-methylated CpGs. (Fig. 3C)—implying that hypo-methylation of CpGs during osteoblastogenesis correlates with active gene expression.
Identifying a link between DNA methylation and TF-binding sites
We investigated whether our DNA methylation signature could identify key TFs that drive or elicit MSC differentiation towards osteoblasts. To do this we designed an R script named REINS that overlaps and normalizes changes in DNA methylation with the TF-binding sites determined by ChIP-seq. For our analysis we used ChIP-seq data derived from the ENCODE project (161 TFs studied across 91 different cell-types). To validate the script we used publicly available data of DNA methylation from different tissues and cells, namely; adipose tissue, muscle, pancreas, thymus, and human pluripotent stem cells (hPSCs) . The polycomb repressive complex 2 (PRC2) related TFs (EZH2, SUZ12, CTBP2) have been tightly associated with the epigenetic repression of stem cell genes during cell differentiation . Thus, the identification of the differential activation state of these TFs [30, 31] on hPSCs (hypo-methylated > hyper-methylated) in comparison with other tissues (hyper-methylated > hypo-methylated) was considered a positive validation (Additional file 8A–C). Likewise, the prediction of the potential critical role of the transcription factor ZEB1 on the thymus (hypo-methylated > hyper-methylated) was also considered a positive validation (Additional file 8D). After the validation, we used REINS with the DNA methylation data of MSCs differentiated to osteoblasts. TFs were sorted according to their normalized ratio of enrichment for binding sites for hypo-/hyper-differentially methylated CpGs, with values > 0 more enriched for de-methylation and < 0 more enriched for methylation (Fig. 3D). The top significant (p < 0.001) TFs enriched, whose binding sites were relatively more hypo-methylated, included the bone anabolism related TFs SMARCC1, MAFK, JUNB, FOS, JUN and STAT3 . Factors that promote osteoblast differentiation are frequently inhibitors of adipocyte differentiation, and vice versa . Accordingly, we noted that a significant proportion of the top significant (p < 0.001) TFs whose binding sites were hyper-methylated were associated with the promotion of adipogenesis (EZH2, ZEB1, NFY, ZNF143, SIN3A, SREBP1, SUZ12, SAP30 and CTBP2) [33,34,35,36,37,38,39].
Expression of ZEB TFs during MSC differentiation to osteoblasts
To validate the data obtained with REINS, we picked the TF most associated to DNA hyper-methylation during osteoblastogenesis, ZEB1 (log2RRT = -− 3.053). Thus, we evaluated the expression of ZEB1 during our in vitro differentiation model. Counterintuitively, based on the methylation data, ZEB1 mRNA expression increased during osteoblastogenesis (Fig. 4A), which was also positively correlated with the induction of ALPL expression, a differentiation marker (Fig. 4B). Interestingly, ZEB1 consensus binding sites can also be occupied by a related TF with repressor activity, ZEB2. ZEB2 mRNA (ΔCt value) expression (data not shown) was relatively higher than ZEB1 in MSC but showed a similar increase during osteoblastogenesis (Fig. 4C). However, early in the osteoblastogenesis the kinetics of ZEB1 and ZEB2 induction differed, as evidenced by the significantly different induction ratio ZEB1/ZEB2 at day 3 (Fig. 4D). To examine how the induction of both ZEB factors, and their different induction-kinetics, could be involved in osteoblastogenesis we constructed a computational model, which included known functions of ZEB1/ZEB2 in bone morphogenic protein (BMP) signalling  (Additional file 9). Since we only wanted to investigate the coexistence and interaction of both ZEB TFs in this model, we did not include other signalling pathways relevant for differentiation. Confirming the experimental data, in the model ALPL expression was induced during differentiation despite the presence of the repressor activity of ZEB2 (Fig. 4E, F), mainly due to the change in the induction kinetics of ZEB1 and ZEB2.
Role of ZEBs TFs on MSC differentiation
Since ZEB1 is involved in both adipose and bone metabolism and the information on the function of ZEB2 is limited we studied the contribution of both TFs to the differentiation of MSC to adipocytes and osteoblasts. We depleted MSC of either ZEB1 or ZEB2 by siRNA prior to their differentiation to osteoblasts or adipocytes (Fig. 5A, B). After 7 days of differentiation the depletion of ZEB1 expression did not alter the expression of osteoblast differentiation markers genes but did significantly reduce the expression of adipogenic markers (Fig. 5C–F). ZEB2 depletion enhanced ALPL expression (Fig. 5C) but not RUNX2 (Fig. 5D) and though adipogenic marker genes were increased this was not significant (Fig. 5E, F) unless compared to their expression level following ZEB1 loss.
Although the transient depletion of either ZEB1 or ZEB2 did not significantly alter the mineralization in fully differentiated (day 21) osteoblasts (Fig. 5G), significant differences were observed in the lipid accumulation between fully differentiated adipocytes with an inhibited expression of ZEB2 compared with ZEB1 (Fig. 5H). These data again pointed that the alteration of ZEB1/ZEB2 relative levels at the start of differentiation affected the outcome of the adipogenic process.
Bone ZEB1 expression positively correlates with body weight and BMI
Regarding the observed implication of ZEB1 in adipogenic metabolism, we aimed to determine whether the expression of this TF was modulated in certain pathologies or physiological situations where bone integrity is affected. Neither osteoarthritic nor osteoporotic bone exhibited variations of ZEB1 expression in comparison to healthy bone (Fig. 6A, B).
As previously described, ZEB1 is involved in obesity development. Considering that obesity is defined as increased fat accumulation, we used body weight to investigate whether bone ZEB1 expression was linked to this pathology, identifying a significant correlation between the expression of this TF and body weight (Fig. 6C). Moreover, to discriminate among healthy, overweight, and obese patients we explored the same correlation using BMI instead of weight, which further improved the correlation coefficient (Fig. 6D).
Obesity has been associated with bone adiposity . Likewise, bone adiposity has been linked to the activation of PPARγ . Thus, we aimed to determine whether PPARG expression correlated with ZEB1 expression in bone samples. Data obtained showed a positive correlation between both genes (Fig. 6E), however PPARG expression did not correlate with weight nor BMI (Fig. 6F, G), indicating a relevant role for ZEB1 in obesity-associated bone alterations.
Here we have described the methylation profile of MSCs undergoing osteoblastogenesis and shown that changes in DNA-methylation occur outside of CpG islands. Likewise, we have associated these methylation changes to established chromatin states and linked the changes in DNA-methylation and gene expression. We have also developed a software tool to investigate the relationship between TFs and DNA-methylation profiles, which we used to define a set of candidate TFs potentially linked to the activation or repression of osteoblastogenesis. Accordingly, we validated the role of one of these inactivated TFs, ZEB1, and its opposing counterpart, ZEB2, on the differentiation of MSCs. We observed that modulating the ratio of these opposing factors affected both osteo- and adipogenesis, experimentally validating our bioinformatic tool.
Bone alterations are a common link amongst many musculoskeletal pathologies and have been associated with increased patient fragility . Further, bone-marrow adiposity increases with age in part because of an increase in MSC differentiation to adipocytes rather than osteoblasts [4, 42, 43]. Epigenomic state has been studied and no significant methylation changes were found in adipogenesis [44, 45], however, similar osteoblastogenic information is limited . As a result, we determined DNA-methylation signature of both processes with an aim of identifying regulatory elements and TFs involved.
In this work we were able to determine the changes in the methylation profile associated with MSC osteoblastogenesis, which revealed that the methylation remodelling occurred outside of the CpG islands, consistent with previous reports describing methylation differences between tissues , and suggested that DNA methylation could contribute to, or reinforce, the differentiation process. Additionally, by overlaying chromatin profiling data , we observed that osteoblastogenic DNA methylation remodelling (both hypo- and hyper-methylation) occurred at enhancer chromatin states. Interestingly, an increase in methylation was evident in the CpGs of insulator, poised promoter, and polycomb-repressed states, with the opposite true for transcription elongation state. This not only added a new layer of evidence to the involvement of DNA methylation on the regulation of osteoblastogenesis but also described the insulators as important points of regulation. DNA-binding of the well-characterized insulator CCCTC-binding factor is highly sensitive to DNA methylation , which suggests that osteoblastogenesis-associated methylation changes could alter CCCTC-binding factor binding, and therefore alter 3D genome architecture and thus influence function [47,48,49].
Unexpectedly, we did not identify any significant changes in the DNA methylation profile during adipogenesis. Supporting this observation, Noer et al.  described that DNA methylation of adipogenic promoters did not reflect transcriptional status, nor potential for gene expression, in adipocytes differentiated from MSCs. Likewise, it was determined that DNA methylation remained stable during adipocyte differentiation, implying that DNA methylation may not be a determinant of the adipogenic differentiation process . Conversely, crude pharmacological inhibition of DNA methylation (with 5-aza-2′-deoxycytidine) in 3T3-L1 preadipocytes inhibited adipogenesis but promoted osteoblastogenesis . MSCs are clearly delicately balanced for their differentiation commitment with regards adipo-osteogenic differentiation of MSCs [12, 13]. Our data suggests that adipogenesis can be achieved without a substantial modification of the default methylation profile of MSCs and that this phenotype is somewhat plastic.
When MSC undergo osteogenic differentiation, the phenotype may be more stable due to its maintenance through a defined DNA methylation programme. Accordingly, during osteoblastogenesis we observed a negative correlation between DNA methylation and medium to highly expressed genes, and a significant link between DNA hypo-methylation and DNA regions related to active gene expression (RNAPOLII regions). Supporting this hypothesis we observed through our novel bioinformatic tool “REINS” that DNA hyper-methylation during osteoblastogenesis was enriched at binding sites of several pro-adipogenic transcriptional regulators [33,34,35,36,37,38,39, 52,53,54]—thus having the potential to inhibit their DNA-binding and therefore activity. Also enriched in osteoblastogenic hyper-methylated regions were binding sites for the transcriptional repressor SETDB1, a H3–K9 histone methyltransferase. SETDB1 binds to methylated DNA via a methyl-CpG-binding domain , thus enrichment of this factor is consistent with its described role as an inhibitor of adipogenesis . Although REINS was validated because it identified the inactive epigenetic repression state of PRC2 TFs on hPSCs, the script is currently limited to the data of the 161 TFs present in the ChIP-seq data from the ENCODE project. As a result, the analysis is biased against TFs that regulate osteoblast metabolism and differentiation for which ChIP-seq data is lacking.
Among the TF potentially inactive or repressed during the osteoblastogenesis, we studied ZEB1, which had the largest enrichment in the ratio DNA hyper-methylated/ hypo-methylated predicted by REINS. ZEB1 is a zinc-finger protein involved in adipocyte differentiation in mouse, as well as in obesity development in humans [33, 57]. It is also involved in thymus and skeletal development [58, 59] with its deletion in mouse being related to craniofacial abnormalities (shortened jaw), limb defects, and fusion of ribs . Consistent with these data, REINS also predicted the potential activation state of this TF on thymus and adipose tissues (Additional file 8D). However, the potential inactive state predicted for ZEB1 in osteoblastogenesis by the tool was contrary to its reported activities promoting BMP2-mediated ALPL expression , skeletal abnormalities presented in ZEB1-null mice , and our observed increase in ZEB1 mRNA expression during osteoblastogenesis, which correlated with ALPL mRNA expression. In fact, ZEB1, along with orthologous ZEB2, are members of the ZEB family of transcription factors, and both can bind to E-box sequences CACCT(G). ZEB2 deficiency in humans has also been associated with craniofacial abnormalities, such as the excessive growth of the jaw (Mowat-Wilson syndrome) . Interestingly, the functionality of ZEB proteins appear cellular-context dependent and, of relevance here, ZEB1 and ZEB2 have reported opposing activities (enhancer and repressor, respectively) upon activation of osteoblastic TGFβ or BMP signalling pathways [40, 61, 62]. Bone marrow has a higher level of ZEB2 expression compared to ZEB1 expression , which is consistent with their expression levels in our MSCs. Although both factors increased in mRNA expression during osteoblastogenesis, the ratio of these inductions (ZEB1 vs ZEB2) significantly increased at day 3 of differentiation. These data may suggest ZEB1 could be more active in early differentiation. At the late stages of differentiation, however, when ZEB2 would be predicted to dominate, we found a potential link of ZEB binding sites with hyper-methylation, which we speculate is to reduce the impact of ZEB2 acting repressively on osteoblastogenesis or to regulate the potential ZEB1-mediated adipogenic behaviour described previously [33, 57]. In line with this, our computational model revealed that the presence of the repressive ZEB2 to be compatible with the induction of ALPL along the differentiation process.
To examine this hypothesis, we depleted ZEB1 and ZEB2 in MSCs prior to inducing their differentiation to osteoblasts and adipocytes. Removal of ZEB2 enhanced expression of the osteoblast differentiation marker ALPL, again supporting its role as an osteoblastogenic repressor. Interestingly, ZEB1 depletion did not affect osteoblast differentiation or mineralization suggesting that the role of ZEB1 during the process could be limited or only evident upon removal of the repressive ZEB2.
To the contrary, ZEB1 depletion in MSCs significantly reduced the expression of the adipogenic differentiation markers, consistent with its role as an adipogenic promoter . ZEB1 depletion partially reduced lipid accumulation, though this did not reach significance possibly due to the transient nature of siRNA-mediated depletion after the 14 days of differentiation required for lipid generation. Meanwhile, ZEB2 depletion increased adipogenic gene expression and oil staining when compared with the levels in ZEB1 depleted cells, again suggesting that the modulation of the ratio of the expression of these TFs in the MSCs regulates the adipogenic process.
Regarding the essential role of ZEB1 in adipogenesis, we observed a positive correlation between bone ZEB1 expression and body weight or BMI. These data are consistent with the described relation between ZEB1 and obesity  and evidence how this metabolic pathology not only affects fat tissue, but also bone. Underpinning these correlations, we have found that ZEB1 also correlated with PPARG in bone samples. In agreement with this, it has been reported that ZEB1 knockdown diminished PPARG expression in 3T3-L1 cells . Since PPARG is a major adipocytic inducer, and its activation has been associated with bone adiposity and fragility , the positive correlation of ZEB1 with this TF could support the role for ZEB1 in bone adiposity. Interestingly, PPARG expression did not correlate with weight or BMI, which suggests that the obesity effects on bone adiposity might imply ZEB1 modulation rather than major adipogenesis inducer PPARG.
In conclusion, in this work we have characterised the DNA-methylation changes that occur during MSC-osteoblastogenesis. We propose that DNA methylation may be important for imprinting and maintaining the osteoblastic phenotype but not that of adipocytes. We present a script that could help to identify key TFs-associated with a given process through the analysis of DNA-methylation. Moreover, we have validated the algorithm on tissue-specific and osteoblastogenic DNA-methylation signatures, identifying ZEB TF ratios as modulators of differentiation of MSCs to osteoblasts and adipocytes, and ZEB1 as a critical TF in obesity-related bone adiposity.
Availability of data and materials
The datasets generated during and/or analysed during the current study are not publicly available but are available from the corresponding author on reasonable request.
Loeser RF, Gandhi U, Long DL, Yin W, Chubinskaya S. Aging and oxidative stress reduce the response of human articular chondrocytes to insulin-like growth factor 1 and osteogenic protein 1. Arthritis Rheumatol. 2014;66(8):2201–9. https://doi.org/10.1002/art.38641.
Shum LC, White NS, Nadtochiy SM, de Bentley KLM, Brookes PS, Jonason JH, et al. Cyclophilin D knock-out mice show enhanced resistance to osteoporosis and to metabolic changes observed in aging bone. PLOS ONE. 2016;11(5):e0155709. https://doi.org/10.1371/journal.pone.0155709.
Woolf AD, Zeidler H, Haglund U, Carr AJ, Chaussade S, Cucinotta D, et al. Musculoskeletal pain in Europe: its impact and a comparison of population and medical perceptions of treatment in eight European countries. Ann Rheumatic Dis. 2004;63(4):342–7.
Rosen CJ, Bouxsein ML. Mechanisms of disease: is osteoporosis the obesity of bone? Nat Clin Practice Rheumatol. 2006;2(1):35–43. https://doi.org/10.1038/ncprheum0070.
Zhuo Q, Yang W, Chen J, Wang Y. Metabolic syndrome meets osteoarthritis . Nature Reviews Rheumatology. Nature Publishing Group; 2012 . 8: 729–37. Accessed 07 Oct 2021.https://www.nature.com/articles/nrrheum.2012.135.
Goodship AE, Cunningham JL, Oganov V, Darling J, Miles AW, Owen GW. Bone loss during long term space flight is prevented by the application of a short term impulsive mechanical stimulus. Acta Astronaut. 1998;43(3–6):65–75.
Lee DS, Choi H, Han BS, Kim WK, Lee SC, Oh KJ, et al. C-Jun regulates adipocyte differentiation via the KLF15-mediated mode. Biochem Biophys Res Commun. 2016;469(3):552–8. https://doi.org/10.1016/j.bbrc.2015.12.035.
Yamamoto Y, Turkiewicz A, Wingstrand H, Englund M, Englund M. Fragility fractures in patients with rheumatoid arthritis and osteoarthritis compared with the general population. J Rheumatol. 2015;42(14):2055–8.
Benet-Travé J, Domínguez-García Á, Sales-Pérez JM, Orozco-Delclós R, Salleras-Sanmartí L, et al. In-Hospital Case-Fatality of Aged Patients with Hip Fracture in Catalonia. Spain Eur J Epidemiol. 1997;13(6):681–6.
Chiba K, Uetani M, Kido Y, Ito M, Okazaki N, Taguchi K, et al. Osteoporotic changes of subchondral trabecular bone in osteoarthritis of the knee: a 3-T MRI study. Osteoporos Int. 2012;23(2):589–97.
Favero M, Giusti A, Geusens P, Goldring SR, Lems W, Schett G, et al. OsteoRheumatology: a new discipline? RMD Open. 2015;1(Suppl 1):e000083.
Chen Q, Shou P, Zheng C, Jiang M, Cao G, Yang Q, et al. Fate decision of mesenchymal stem cells: adipocytes or osteoblasts? Cell Death and Differ. 2016;23(7):1128–39.
Takada I, Kouzmenko AP, Kato S. Molecular switching of osteoblastogenesis versus adipogenesis: implications for targeted therapies. Expert Opin Ther Targets. 2009;13(5):593–603.
Meunier P, Aaron J, Edouard C, Vignon G, et al. Osteoporosis and the replacement of cell populations of the marrow by adipose tissue. A quantitative study of 84 iliac bone biopsies. Clin Orthop Related Res. 1971;80:147–54.
Wan Y. PPARγ in bone homeostasis. Trends Endocrinol Metab. 2010;21(12):722–8.
Rosen ED, MacDougald OA. Adipocyte differentiation from the inside out. Nat Rev Mol Cell Biol. 2006;7(12):885–96.
Deaton AM, Bird A. CpG islands and the regulation of transcription. 2011;25(10).
Jin B, Li Y, Robertson KD. DNA methylation: superior or subordinate in the epigenetic hierarchy? Gene Cancer. 2011;2(6):607–17. https://doi.org/10.1177/1947601910393957.
Barter MJ, Tselepi M, Gómez R, Woods S, Hui W, Smith GR, et al. Genome-wide microRNA and gene analysis of mesenchymal stem cell chondrogenesis identifies an essential role and multiple targets for miR-140–5p. Stem Cells. 2015;33(11):3266–80. https://doi.org/10.1002/stem.2093.
Zeybel M, Hardy T, Wong YK, Mathers JC, Fox CR, Gackowska A, et al. Multigenerational epigenetic adaptation of the hepatic wound-healing response. Nat Med. 2012;18(9):1369–77.
Barter MJ, Gomez R, Hyatt S, Cheung K, Skelton AJ, Xu Y, et al. The long non-coding RNA ROCR contributes to SOX9 expression and chondrogenic differentiation of human mesenchymal stem cells. Development. 2017;144(24):4510–21.
Touleimat N, Tost J. Complete pipeline for infinium® human methylation 450K beadchip data processing using subset quantile normalization for accurate DNA methylation estimation. Epigenomics. 2012;4(3):325–41. https://doi.org/10.2217/epi.12.21.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47–e47.
Du P, Zhang X, Huang C-C, Jafari N, Kibbe WA, Hou L, et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinform. 2010;11(1):1–9. https://doi.org/10.1186/1471-2105-11-587.
Gerstein MB, Kundaje A, Hariharan M, Landt SG,Yan K-KK, Cheng C, et al. Architecture of the human regulatory network derived from ENCODE data. Nature. 2012;489(7414):91–100.
Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, et al. Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011;473(7345):43–9.
Martino D, Saffery R. Characteristics of DNA methylation and gene expression in regulatory features on the infinium 450k Beadchip. bioRxiv. 2015. https://doi.org/10.1101/032862v1.
Li G, Ruan X, Auerbach RK, Sandhu KS, Zheng M, Wang P, et al. extensive promoter-centered chromatin interactions provide a topological basis for transcription regulation. Cell. 2012;148(1):84–98.
Nazor KL, Altun G, Lynch C, Tran H, Harness JV, Slavin I, et al. Recurrent variations in DNA methylation in human pluripotent stem cells and their differentiated derivatives. Cell Stem Cell. 2012;10(5):620–34.
Kim TW, Kang B-H, Jang H, Kwak S, Shin J, Kim H, et al. Ctbp2 modulates NuRD-mediated deacetylation of H3K27 and facilitates PRC2-mediated H3K27me3 in active embryonic stem cell genes during exit from pluripotency. Stem Cells. 2015;33(8):2442–55. https://doi.org/10.1002/stem.2046.
Burdach S, Plehm S, Unland R, Borkhardt A, Staege MS, Müller-Tidow C, et al. Epigenetic maintenance of stemness and malignancy in peripheral neuroectodermal tumors by EZH2. Cell Cycle. 2009;8(13):1991–6. https://doi.org/10.4161/cc.8.13.8929.
McCabe LR, Kockx M, Lian J, Stein J, Stein G. Selective expression of fos- and jun-related genes during osteoblast proliferation and differentiation. Exp Cell Res. 1995;218(1):255–62.
Gubelmann C, Schwalie PC, Raghav SK, Röder E, Delessa T, Kiehlmann E, et al. Identification of the transcription factor ZEB1 as a central component of the adipogenic gene regulatory network. eLife. 2014;3(August 2014):1–30.
Hemming S, Cakouros D, Isenmann S, Cooper L, Menicanin D, Zannettino A, et al. EZH2 and KDM6A act as an epigenetic switch to regulate mesenchymal stem cell lineage specification. Stem Cells. 2014;32(3):802–15. https://doi.org/10.1002/stem.1573.
Kim JB, Spiegelman BM. ADD1/SREBP1 promotes adipocyte differentiation and gene expression linked to fatty acid metabolism. Genes Dev. 1996;10(9):1096–107.
Laudes M, Bilkovski R, Oberhauser F, Droste A, Gomolka M, Leeser U, et al. Transcription factor FBI-1 acts as a dual regulator in adipogenesis by coordinated regulation of cyclin-A and E2F–4. J Mol Med. 2008;86(5):597–608. https://doi.org/10.1007/s00109-008-0326-2.
Lu YH, Dallner OS, Birsoy K, Fayzikhodjaeva G, Friedman JM. Nuclear Factor-Y is an adipogenic factor that regulates leptin gene expression. Mol Metab. 2015;4(5):392–405.
Raghav SK, Waszak SM, Krier I, Gubelmann C, Isakova A, Mikkelsen TS, et al. Integrative genomics identifies the corepressor SMRT as a gatekeeper of adipogenesis through the transcription factors C/EBPβ and KAISO. Mol Cell. 2012;46(3):335–50.
Wang L, Jin Q, Lee J-E, Su I, Ge K. Histone H3K27 methyltransferase Ezh2 represses Wnt genes to facilitate adipogenesis. Proc Nat Academ Sci. 2010;107(16):7317–22.
Postigo AA. Opposing functions of ZEB proteins in the regulation of the TGFβ/BMP signaling pathway. The EMBO J. 2003;22(10):2443–52. https://doi.org/10.1093/emboj/cdg225.
Liu L, Aronson J, Huang S, Lu Y, Czernik P, Rahman S, et al. Rosiglitazone inhibits bone regeneration and causes significant accumulation of fat at sites of new bone formation. Calcified Tissue Int. 2012;91(2):139.
Astudillo P, Ríos S, Pastenes L, Pino AM, Rodríguez JP. Increased adipogenesis of osteoporotic human-mesenchymal stem cells (MSCs) characterizes by impaired leptin action. J Cell Biochem. 2008;103(4):1054–65. https://doi.org/10.1002/jcb.21516.
Rodríguez JP, Garat S, Gajardo H, Pino AM, Seitz G. Abnormal osteogenesis in osteoporotic patients is reflected by altered mesenchymal stem cells dynamics. J Cell Biochem. 1999;75(3):414–23.
Tang QQ, Lane MD. Adipogenesis: from stem cell to adipocyte. Annu Rev Biochem. 2012;81(1):715–36. https://doi.org/10.1146/annurev-biochem-052110-115718.
van den Dungen MW, Murk AJ, Kok DE, Steegenga WT. Comprehensive DNA methylation and gene expression profiling in differentiating human adipocytes. J Cell Biochem. 2016;117(12):2707–18. https://doi.org/10.1002/jcb.25568.
Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, et al. The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009;41(2):178–86.
Wang H, Maurano MT, Qu H, Varley KE, Gertz J, Pauli F, et al. Widespread plasticity in CTCF occupancy linked to DNA methylation. Genom Res. 2012;22(9):1680–8.
Guo Y, Xu Q, Canzio D, Shou J, Li J, Gorkin DU, et al. CRISPR Inversion of CTCF sites alters genome topology and enhancer/promoter function. Cell. 2015;162(4):900–10.
Nichols MH, Corces VG. A CTCF code for 3D genome architecture. Cell. 2015;162(4):703–5.
Noer A, Sørensen AL, Boquest AC, Collas P. Stable CpG hypomethylation of adipogenic promoters in freshly isolated cultured, and differentiated mesenchymal stem cells from adipose tissue. Mol Biol Cell. 2006;17(8):3543–56. https://doi.org/10.1091/mbc.e06-04-0322.
Chen Y-SS, Wu R, Yang X, Kou S, MacDougald OA, Yu L, et al. Inhibiting DNA methylation switches adipogenesis to osteoblastogenesis by activating Wnt10a. Sci Rep. 2016;6(1):1–12.
Fromm-Dornieden C, von der Heyde S, Lytovchenko O, Salinas-Riester G, Brenig B, Beissbarth T, et al. Novel polysome messages and changes in translational activity appear after induction of adipogenesis in 3T3-L1 cells. BMC Mol Biol. 2012;13(1):1–12. https://doi.org/10.1186/1471-2199-13-9.
Jimenez MA, Åkerblad P, Sigvardsson M, Rosen ED. Critical role for Ebf1 and Ebf2 in the adipogenic transcriptional cascade. Mol Cell Biol. 2007;27(2):743–57.
Liu TM, Wu YN, Guo XM, Hui JHP, Lee EH, Lim B. Effects of ectopic nanog and Oct4 overexpression on mesenchymal stem cells. Stem Cells Dev. 2009;18(7):1013–21. https://doi.org/10.1089/scd.2008.0335.
Roloff TC, Ropers HH, Nuber UA. Comparative study of methyl-CpG-binding domain proteins. BMC Genom. 2003. https://doi.org/10.1186/1471-2164-4-1.
Okamura M, Inagaki T, Tanaka T, Sakai J. Role of histone methylation and demethylation in adipogenesis and obesity [Internet]. Vol. 6, Organogenesis. Organogenesis; 2010;p. 24–32
Hinney A, Ziegler A, Oeffner F, Wedewardt C, Vogel M, Wulftange H, et al. Independent confirmation of a major locus for obesity on chromosome 10. J Clin Endocrinol Metab. 2000;85(8):2962–5.
Higashi Y, Moribe H, Takagi T, Sekido R, Kawakami K, Kikutani H, et al. Impairment of T Cell Development in δEF1 Mutant Mice. J Exp Med. 1997;185(8):1467–80.
Takagi T, Moribe H, Kondoh H, Higashi Y. DeltaEF1, a zinc finger and homeodomain transcription factor, is required for skeleton patterning in multiple lineages. Development. 1998;125:21–31.
Mowat DR, Wilson MJ, Goossens M. Mowat-Wilson syndrome. J Med Genet. 2003;40(5):305–10.
Gheldof A, Hulpiau P, van Roy F, De Craene B, Berx G. Evolutionary functional analysis and molecular regulation of the ZEB transcription factors. Cell Mol Life Sci. 2012;69(15):2527–41. https://doi.org/10.1007/s00018-012-0935-3.
Postigo AA, Depp JL, Taylor JJ, Kroll KL. Regulation of Smad signaling through a differential recruitment of coactivators and corepressors by ZEB proteins. EMBO J. 2003;22(10):2453–62. https://doi.org/10.1093/emboj/cdg226.
Uhlén M, Fagerberg L, Hallström BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Tissue-based map of the human proteome. Science. 2015. https://doi.org/10.1126/science.1260419.
The authors acknowledge Dr. Louise Reynard for advice.
This research was supported by the JGW Patterson Foundation and Arthritis Research UK (programme Grant No. 19424), National Institute for Health Research Newcastle Biomedical Research Centre and the Northumberland, Tyne, and Wear Comprehensive Local Research Network. R.G. work was funded by the Instituto de Salud Carlos III through a Miguel Servet II programme. A.A-P. is supported by Fundación IDIS. R.G. and A.A-P. work was funded by the European Union via “Fondo de Investigación Sanitaria” from Fondo Europeo de Desarrollo Regional (FEDER) (PI16/01870, PI19/01446), and Fundación Mutua Madrileña (MMA2018_2).
Ethics approval and consent to participate
The Ethics Committee for Research at Santiago-Lugo Area approved the protocol for bone samples obtention. Informed consent was obtained from all patients or patients’ families (Registration No 2016/258).
Consent for publication
All the authors have given their consent for the publication of this work.
The authors have no relevant financial or non-financial interests to disclose.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
REINS NCL 2020. Script, Regulatory Element Interrogation Script (REINS), available as an R Markdown document, to explore the relationship between the distribution of differentially methylated CpGs and TF-binding and thus activity during the differentiation process.
Additional file 2:
Model design. A computational model of Bmp2 signalling was constructed to evaluate ALP induction in the context of ZEB TFs activity.
Additional file 3:
CPG list and overlapping genes. List of differently methylated CPGs associated to osteoblast differentiation and their overlapping genes.
Additional file 4:
CpG_OB_AP_MSC. Methylation status (M-values) of CPGs differently methylated across the MSC differentiation to osteoblast and adipocyte.
Additional file 5:
RNA array. Differently expressed genes between stem cells differentiated (21 days) or not to osteoblast.
Additional file 6: Figure S1.
Observed vs expected ratio of CpGs enrichment among the different chromatin states and their cumulative methylation profile in MSCs. A) The bars represent the ratio of the percentage of differently methylated CpGs vs the percentage of the total CpGs studied for each chromatin state. B) Cumulative plot of the methylation of the CpGs overlapping top enriched (Strong_ Enhancer, Weak_ Enhancer, and Txn_Transition) and depleted chromatin states (Active_promoter and Poised_ Promoter). Y axis: cumulative percentage of CpGs overlapping to each chromatin state. X axis: fold change methylation MSCs vs Osteoblasts.
Additional file 7: Figure S2.
Output of ingenuity pathway analysis. Representation of the significant enrichment on ‘inhibition of MMPs’ as well as ‘WNT pathway activation’.
Additional file 8: Figure S3.
Validation of REINS algorithm. Transcription factor relevance (TFR) metric for hypo (Ho)- and hypermethylated (Hr) CpGs associated to MSCs differentiation to osteoblasts (OB), to different tissues (muscle, adipose tissue, pancreas, thymus) and to human pluripotent stem cells (hPSC). A) TFR for the repressor TF CTBP2. B) TFR for the TF EZH2 a member of the polycomb repressive complex 2 (PRC2). C) TFR for the TF SUZ12 a member of the polycomb repressive complex 2 (PRC2). D) TFR for the TF ZEB1.
Additional file 9: Figure S4.
Computational simulation of Zeb1 and Zeb2 interaction in the context of BMP2-induced target gene expression.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Gómez, R., Barter, M.J., Alonso-Pérez, A. et al. DNA methylation analysis identifies key transcription factors involved in mesenchymal stem cell osteogenic differentiation. Biol Res 56, 9 (2023). https://doi.org/10.1186/s40659-023-00417-6
- Stem cell differentiation
- Bone adiposity