Transcriptome Analysis of Dnmt3l Knock-Out Mice Derived Multipotent Mesenchymal Stem/Stromal Cells During Osteogenic Differentiation

Abstract

Multipotent mesenchymal stem/stromal cells (MSCs) exhibit great potential for cell-based therapy. Proper epigenomic signatures in MSCs are important for the maintenance and the subsequent differentiation potential. The DNA methyltransferase 3-like (DNMT3L) that was mainly expressed in the embryonic stem (ES) cells and the developing germ cells plays an important role in shaping the epigenetic landscape. Here, we report the reduced colony forming ability and impaired in vitro osteogenesis in Dnmt3l-knockout-mice-derived MSCs (Dnmt3l KO MSCs). By comparing the transcriptome between undifferentiated Dnmt3l KO MSCs and the MSCs from the wild-type littermates, some of the differentially regulated genes (DEGs) were found to be associated with bone-morphology-related phenotypes. On the third day of osteogenic induction, differentiating Dnmt3l KO MSCs were enriched for genes associated with nucleosome structure, peptide binding and extracellular matrix modulation. Differentially expressed transposable elements in many subfamilies reflected the change of corresponding regional epigenomic signatures. Interestingly, DNMT3L protein is not expressed in cultured MSCs. Therefore, the observed defects in Dnmt3l KO MSCs are unlikely a direct effect from missing DNMT3L in this cell type; instead, we hypothesized them as an outcome of the pre-deposited epigenetic signatures from the DNMT3L-expressing progenitors. We observed that 24 out of the 107 upregulated DEGs in Dnmt3l KO MSCs were hypermethylated in their gene bodies of DNMT3L knock-down ES cells. Among these 24 genes, some were associated with skeletal development or homeostasis. However, we did not observe reduced bone development, or reduced bone density through aging in vivo. The stronger phenotype in vitro suggested the involvement of potential spreading and amplification of the pre-deposited epigenetic defects over passages, and the contribution of oxidative stress during in vitro culture. We demonstrated that transient deficiency of epigenetic co-factor in ES cells or progenitor cells caused compromised property in differentiating cells much later. In order to facilitate safer practice in cell-based therapy, we suggest more in-depth examination shall be implemented for cells before transplantation, even on the epigenetic level, to avoid long-term risk afterward.

Introduction

Stem cells that are capable of self-renewal and differentiation, play very important roles to maintain dynamic but sustainable tissue homeostasis (Collins et al., 2005). This delicate balance relies on cooperative interactions between transcription networks and epigenetic regulations, including DNA methylation, histone modifications and non-coding RNA-mediated modulation of gene expression (Reik, 2007; Lunyak and Rosenfeld, 2008; Bloushtain-Qimron et al., 2009; Atlasi and Stunnenberg, 2017). Multipotent mesenchymal stem/stromal cells (MSCs) have a trilineage differentiation ability toward osteoblasts, chondrocytes and adipocytes (Pittenger et al., 1999). In accordance with their ability to differentiate into multiple lineages, and having immunomodulatory function(s), MSCs can be isolated from various origins, such as bone marrow, adipose tissues (Zuk et al., 2002), umbilical cord blood (Erices et al., 2000), skeletal muscle (Young et al., 2001), amniotic fluid (In t Anker et al., 2003), placenta (Miao et al., 2006), and dental pulp (Lin et al., 2011; Berebichez-Fridman and Montero-Olvera, 2018; Spagnuolo et al., 2018). With the long history of successful MSC-containing bone marrow transplantations, as well as the availability from other accessible origins like tooth, placenta, umbilical cord blood, amniotic fluid and adipose tissues, MSCs are one of the top candidate cell types available for allogenic and autologous cell therapy in regenerative medicine (Ballini et al., 2017; Han et al., 2019). The therapeutic potential of MSCs is mainly based on their ability to modulate immune response (Abdi et al., 2008; Gao et al., 2016). The strong anti-inflammatory effect and paracrine activity from their secretory proteins and functional RNAs can be delivered into remote organs or tissues via extracellular vesicles (Chen et al., 2019). Preclinical and clinical trials on MSCs based therapy have been carried out in many countries, in order to improve the treatment for the injured bone, cartilage, and the autoimmune diseases.

However, one of the limitations of using MSCs in regenerative medical applications is their gradual loss of multipotency and therapeutic properties under culture conditions. MSCs cultured for prolonged periods or isolated from elder individual are less efficient in their trilineage differentiation ability, with the impairment in osteogenesis being particularly significant (Digirolamo et al., 1999; Bonab et al., 2006; Zhou et al., 2008; Yang, 2018). Optimization of culture conditions can sustain MSC properties better in vitro (Sotiropoulou et al., 2006), and make them safer and more suitable for clinical applications. These include the introduction of better culture surface (Engler et al., 2006; Lee et al., 2017), hypoxia condition (Wang et al., 2020), providing scaffold and other biomaterials (Meinel et al., 2004; Marrelli et al., 2016), to maintain better multipotency or differentiation outcome for the cultured MSCs. In addition, the replacement of FBS by chemically defined or standardized supplements (Bieback et al., 2009; Marrazzo et al., 2016) can facilitate the clinical-grade production of MSCs.

While the culture condition can be optimized to certain extent, the intrinsic defects from the isolated MSCs cannot be easily fixed. Here we report an unexpected observation of compromised osteogenesis differentiation ability of MSCs isolated from DNMT3L deficient mutant mice. DNMT3L is a germ and ES cell enriched epigenetic cofactor (Bourc’his et al., 2001; Hata et al., 2002; Liao et al., 2014). We and others have demonstrated that DNMT3L maintains the quiescence of spermatogonial progenitor cells and prevents exhaustion of stem cell populations to maintain male germ line homeostasis (Bourc’his and Bestor, 2004; Liao et al., 2014). Dnmt3l knock-out mice are infertile (Bourc’his and Bestor, 2004; Webster et al., 2005; Hata et al., 2006), but otherwise develop normally into adulthood without reported somatic phenotypes. DNMT3L does not have enzymatic activity but interacts with DNMT3A and DNMT3B to facilitate de novo DNA methylation and thus influences gene expression (Chedin et al., 2002; Guenatri et al., 2013). DNMT3L binds to histone H3 tails in a H3K4methylation sensitive manner, and recruits other histone modifiers through its PHD domain (Aapola et al., 2000; Ooi et al., 2007; Otani et al., 2009; Hashimoto et al., 2010; Zhang et al., 2010). We further demonstrated that ectopic DNMT3L expression can promote the assembly of the HDAC1/TRIM28/SETDB1/DNMT3A/DNMT3L complex and repress transcription of newly infected retroviral sequence independent of DNA methylation (Kao et al., 2014). There has been very limited description of potential DNMT3L functions beyond germ lines and ES cells, partly due to the difficulties in demonstrating its expression in specific progenitor cell types in somatic lineages.

Recently we demonstrated that transient expression of ectopic DNMT3L in later passaged MEFs were sufficient to cause long term epigenomic landscape changes and halt senescence progression (Yu et al., 2020). The transiently expressed DNMT3L facilitated the short-term formation of DNMT3L-DNMT3A-KAP1-SETDB-HDAC1 complex as well as guiding them to certain endogenous retroviruses and retrotransposons to introduce H3K9me3 and reduce histone acetylation in aging fibroblasts (Kao et al., 2014; Yu et al., 2020). DNMT3L also interacted with polycomb group members to facilitate repressive H3K27me3 modifications on certain aging associated derepressed genes. The long-term repressive histone modifications and dramatically prolonged cell proliferation still maintained long after the ectopic DNMT3L is silenced (Yu et al., 2020). In the current study, we tackled a potential long-term effect of transient endogenous DNMT3L expression in ES cells and progenitor cells, associated with cellular property changes of the Dnmt3l KO mice derived MSCs (Dnmt3l KO MSCs) under culture condition.

We compared the in vitro CFU-F formation rate, proliferation and differentiation abilities between MSCs derived from Dnmt3l KO mice and their WT littermates. Significantly reduced CFU-Fs and impaired osteogenic ability were observed in Dnmt3l KO mice derived BM-MSCs as well as in another epiphysis-derived MSCs (EMSCs) population (Cheng et al., 2012). The molecular mechanisms associated with this impaired osteogenic ability were investigated by comparing and contrasting the transcriptomes of undifferentiated and differentiating MSCs derived from Dnmt3l KO mice and their WT littermates. The injury repair ability of the endogenous MSCs between wild type and Dnmt3l KO mice were also compared. As the canonical DNMT3L protein was not detectable in cultured MSCs, various hypotheses were proposed and tested to determine how transiently expressed DNMT3L in ES cells or subsequent progenitor cells contribute to the long-term effect on the subsequent differentiating cell lineages, and their adaptation toward in vitro culture conditions.

Results

Multipotent MSCs Derived From Dnmt3l KO Mice Have Reduced Colony Forming Capacity and Osteogenic Potential Suggesting Compromised Potency in vitro

Multipotent MSCs are capable of self-renewal and can differentiate into multiple lineages in vitro. We compared and contrasted BM-MSC properties between Dnmt3l KO mice and their WT littermate. We did not detect differences in the frequency of MSC cell surface markers tested between genotypes (Supplementary Figure 1). A significantly reduced CFU-F was observed from BM-MSCs and EMSCs (another subtype of MSCs derived from epiphysis from Dnmt3l KO mice than those obtained from WT littermates (Figures 1A,B and Supplementary Figure 2A). To clarify whether Dnmt3l KO MSCs and EMSCs have an overall lower proliferation activity, or just proliferate poorly in low density, we checked the proliferation curves of BM-MSCs and EMSCs from paired Dnmt3l KO mice and wild type littermates. We did not observe major proliferation differences from paired BM-MSCs at passage 4 (Figure 1C), where we used for the differentiation assay. On the other hand, at passage 6, Dnmt3l KO EMSCs had reduced proliferation activity compared to those from wild type animals (Supplementary Figure 2B). This is consistent with our previous observations of premature senescence phenotype in cultured Dnmt3l KO MEFs (Liao et al., 2015), as well as slightly shorter lifespan for Dnmt3l KO compared to wild type littermates (data not shown).

FIGURE 1

In addition to self-renewal property, MSCs are able to differentiate into osteoblasts, chondrocytes and adipocytes in vitro upon induction (Caplan, 1991). We quantified calcium precipitation by ARS staining after 14 days of osteogenic induction and showed that Dnmt3l-KO-mice-bone-marrow-derived MSCs (Dnmt3l KO MSCs), and Dnmt3l KO EMSCs dramatically lost their osteogenic ability (Figure 1D and Supplementary Figure 2C). The compromised osteogenic abilities of Dnmt3l KO BM-MSCs could not be attributed merely to reduced proliferation activity of cells in the passage used for differentiation analysis (Figure 1C), which is also indirectly supported by the reasonable adipogenic ability observed on these cells (Supplementary Figure 3). It is reported that aging MSCs favor adipogenic differentiation over osteogenic differentiation (Infante and Rodriguez, 2018). The severely compromised osteogenesis and relatively mild adipogenesis defects are reminiscent of a trend for premature aging event.

Lack of Canonical DNMT3L Expression in MSCs

We observed very strikingly reduced CFU-F and osteogenesis deficit of MSCs isolated from Dnmt3l KO mice (Figure 1B). The intuitive explanation could be that DNMT3L is expressed and plays a direct role in cultured MSCs when the cells are out of their endogenous niches. However, the expression of DNMT3L in MSCs has not been documented. We performed strand-specific RNA sequencing analysis for undifferentiated WT and Dnmt3l KO MSCs at passage 4, as well as for the differentiating cells of the two genotypes 3 days after osteogenic induction. We confirmed that Dnmt3l was not expressed in Dnmt3l KO MSCs. However, low levels of Dnmt3l associated sequences could be identified from undifferentiated and differentiating wild type BM-MSCs (Figure 2A). Interestingly, all of the Dnmt3l RNA sequences from WT MSCs mapped to exon 9a onward, missing exons 1–8 (Figure 2A). The absence of exons 1–8 and the presence of exon 9a–10–11–12–13 of the Dnmt3l transcript in WT BM-MSCs were supported by our RT-qPCR analysis with various primer pairs (data not shown). Three alternative transcripts of Dnmt3l starting from exon 9b were described during sperm development (Shovlin et al., 2007). However, a potential isoform starting from exon 9a instead of exon 9b has not been previously described.

FIGURE 2

On the protein level, we could not detect the canonical ∼50 kDa DNMT3L protein (Figure 2B), consistent with the lack of RNA expression from exons 1–8. Neither could we detect distinct immunoblotting signals of the estimated 126-AA-long protein from the possible Dnmt3l isoform starting from exon 9a (Figure 2C). We have not excluded the possibility that this exon 9a initiated Dnmt3l RNA isoform may have a function. Our preliminary RT-qPCR data demonstrated higher Dnmt3l expression from the 3′-end exons in EMSCs than in BM-MSCs, associated with more severe potency defects in Dnmt3l EMSCs.

Transcriptome Analysis of Undifferentiated and Differentiating MSCs Derived From Dnmt3l KO Mice Compared to Those of Wild-Type Littermates

In order to gain insight into the mechanistic significance of the in vitro osteogenesis defects in Dnmt3l KO MSC, we performed WGCNA, an unbiased and unsupervised analysis that identifies modules corresponding to clusters of co-expressed transcripts among multiple transcriptomes. The results showed two major gene modules (Figure 3A, marked by turquoise and blue). The first module contains 928 genes enriched in “cell, morphogenesis and stem cell differentiation” by gene ontology analysis. The second module has 898 genes that are specifically associated with “the maintenance of DNA methylation.” In addition, PCA revealed four distinct groups signifying that cell differentiation and Dnmt3l genotype contribute to significant transcriptome changes (Figure 3B).

FIGURE 3

In addition to global gene expression patterns observed by PCA, DEG (at least log₂ fold change > 1 with FDR < 0.05) between undifferentiated and differentiating MSCs, or between different Dnmt3l genotypes of the same differentiation stage, were indicated in the dendrogram (Figure 3C). In Figure 3C, each row presents a DEG. Approximately 300 genes (Figure 3C, DEG_wt_d3lKO_un group) were differentially expressed in Dnmt3l KO MSCs in the undifferentiated state, of which 106 were upregulated and 202 were downregulated. These DEGs were highly correlated to receptor binding, cytokine activity, cargo receptor activity, and ECM structural constituent, which may directly or indirectly influence the MSC differentiation potential upon induction of differentiation (Figure 3D). In addition, some of these DEGs were overlapped with signature genes associated with bone related phenotypes as summarized in the MGI mammalian phenotypes database (Table 1). For example, Postn encodes the matricellular protein POSTN, which is downregulated in Dnmt3l KO MSCs. Postn is transcriptionally activated by RUNX2 in the MC3T3-E1 osteoblast precursor cell line, suggesting that POSTN, also known as OSF-2, is involved in early osteoblast development (Stock et al., 2004). This gene is highlighted in the categories of abnormal long bone morphology (MP:0003723), short femur (MP:0003109), and short tibia (MP:0002764). Previous studies also demonstrated that POSTN may play a crucial role in bone formation and homeostasis (Bonnet et al., 2012; Merle and Garnero, 2012; Zhang et al., 2017), which is highly correlated to the jeopardized osteogenic differentiation ability observed in Dnmt3l KO MSCs.

Upon comparing transcriptomes of MSCs derived from Dnmt3l KO and WT mice after 3 days of induction toward osteogenesis, a more significant change was observed (n = 567 DEGs, 239 genes were upregulated and 328 genes were downregulated) (Figure 3C, DEG_wt_d3lKO_D group). We noticed a clear enrichment of these DEGs in molecular binding related properties, including nucleosome and nucleosomal DNA, peptide, kinetochore, odorant, and collagen binding, as well as genes associated with general molecular function regulation. In addition, with enrichment of DEGs in GO terms associated with collagen binding, hydrolase activity and ECM structural constituent, we suggest that MSCs derived from Dnmt3l KO mice may have suboptimal modulatory activities and responsiveness toward ECMs (Figure 3E).

With candidate approach, we also have observed the expression profile of several other bone related biomarkers, including Sp7, Col1a1 and Col1a2, from our sequencing datasets. Dnmt3l genotype associated, statistically significant differences were observed in all of these genes at day 3 of in vitro MSC differentiation toward osteoblast, while significant differences can be found in Col1a1 and Col1a2, as well as the aforementioned Postn genes in MSCs of different Dnmt3l genotype even before differentiation (Supplementary Figure 4).

Searching for Potential DNMT3L-Mediated Epigenetic Memory From ES Cells

Given DNMT3L protein was not detected in MSCs, we explored the possibility of DNMT3L-mediated epigenetic modulation in pluripotent stem cells; we hypothesized that the DNMTs may exhibit a long term effect that could alter the differentiation potential of the cultured MSCs in vitro. We re-analyzed publicly available DNMT3L ChIP-Seq (GSE49178), DNMT3L KD ES cells gene expression array (GSE44643) and MeDIP sequencing (GSE44642) data (Neri et al., 2013).

To understand which genes are directly regulated by DNMT3L, we categorized DNMT3L binding sites in gene bodies, promoters or repeat sequences from DNMT3L ChIP-Seq data from ES cells. We observed that DNMT3L binding sites were enriched in 5′-UTRs, coding regions and TEs (Figure 4A). To verify the lasting impact of DNMT3L from ES cells to MSCs, we first sought to identify the overlap between DNMT3L binding sites in ES cells and DEGs between the WT and Dnmt3l KO MSCs. However, in both undifferentiated and differentiating MSCs, only 11 Dnmt3l genotype associated DEGs were also associated with DNMT3L binding in ES cells (Figure 4B and Supplementary Figure 5); suggesting most DEGs in MSCs were not the direct targets of DNMT3L in ESCs.

FIGURE 4

To investigate whether DNMT3L-dependent gene expression patterns carry over from ES cells to MSCs, DEGs from an expression array of ES cells were compared to MSC RNA-Sequencing data. Only 4 out of 309 genes were differentially expressed in the Dnmt3l KD ES cells and undifferentiated Dnmt3l KO MSCs (Figure 4C), and 7 out of 567 DEGs were identified in differentiating MSCs (Figure 4D), indicating that the role of DNMT3L legacy from ES cells on MSC properties is not a primary effect in gene expression.

DNMT3L is a well-known epigenetic cofactor, functioning within the de novo DNA methylation machinery. With the faithful memory of DNA methylation after cell division, alterations in Dnmt3l genotype-dependent DNA methylation level in ES cells could therefore be one of the most promising long-lasting epigenomic marks that are potentially responsible for the compromised Dnmt3l KO MSCs properties. To understand the influence of DNMT3L on DNA methylation, we analyzed the DNA methylation patterns of WT and Dnmt3l KD ES cells. According to the metagene plots, the methylation of Dnmt3l KD ES cells was generally higher than WT cells around genes, particularly at the promoter region (Figures 4E,F), and the same result was observed in the enrichment analysis (Figure 4G). We further found that these changes in DNA methylation occurred in genes from ES are associated with increased expression (between Dnmt3l KO and WT) in MSC (Figure 4H), indicating a positive correlation between changes in DNMT3L-induced methylation in ES cells and gene expression in MSCs, as a potential signal of DNMT3L legacy.

We further clarified whether the differential gene expressions in in vitro cultured Dnmt3l KO and WT MSCs, a cell type no longer expressing DNMT3L, may have association with DNMT3L-dependent DNA methylation mark in ES cells. Apart from the epigenomic modulation of regulatory sequences, DNA methylation on actively transcribed genes, had also been demonstrated (Bender et al., 1999).

We observed that 24 genes were upregulated in Dnmt3l KO MSCs and hypermethylated at the gene bodies in Dnmt3l KD ES cells (Figure 4I and more comprehensive analysis summarized in Supplementary Figure 6), with a very low p-value < 1.387e-04, indicating the significant overlap based on the hypergeometric distribution test. Of the 24 genes, some gene products have been reported to tie with skeletal development or homeostasis, with a few of them acting directly in MSCs to regulate osteogenic differentiation, namely, HIEVP3 and GATA2 inhibiting BMP-induced osteogenesis (Tolkachov et al., 2018; Li et al., 2020), DLX3 being a negative modulator (Levi and Gitton, 2014), while FGF18 and Wnt10a providing stimulatory signal (Charoenlarp et al., 2017; Jing et al., 2018). Those results suggest that the occurrence of DNMT3L-mediated epigenetic legacy, if any, may be observed in limited genes in cultured MSCs.

Altered Transposable Element Expression Profile Suggesting Changes of Global Epigenomic Landscape

Because of the known transposon element (TE) silencing activity of DNMT3L (Bourc’his and Bestor, 2004), we assessed the possibility of correlating specific families of deregulated TEs to the impaired properties of Dnmt3l KO MSCs. The DETEs with four-fold changes were subjected to unsupervised heat map analysis (Supplementary Figure 7A). Genotype-dependent differences (four-fold changes) were manifested in an undifferentiated but not at differentiating state based on PCA analysis (Supplementary Figure 7B). Even without considering vast copy numbers of each TE sequences, there are more DETEs between WT and Dnmt3l KO genotypes in undifferentiated MSCs than the number of DEGs (Supplementary Figure 7C). DETEs populated in specific subfamilies such as DNA transposon, LINE, SINE and LTR, suggesting the accumulative effects from DNMT3L targeting these four subfamilies of TEs (Figure 5A). The significant reduction of upregulated TEs from undifferentiated Dnmt3l KO MSCs across many TE families (blue bars from Figure 5A) indicated possible epigenomic signature changes across many genomic loci from this cell type. In the differentiating MSCs, DETEs were enriched in the subfamilies of DNA transposon (hAT-Tip 100), LINEs (L2 and CR1), SINE (B4) and ERVK of the LTR family (green and purple bars in Figure 5A).

FIGURE 5

With the known property of DNMT3L in TE modulation, we also analyzed the profile of DNMT3L associated DNA methylation marks from ES cells in relation to TE expressions in cultured MSCs of different Dnmt3l genotype. The DNA methylation is relatively low at most TE regions in ES cells. Dnmt3l KD induced further reduction was particularly enriched in LINE and LTR retrotransposons in ES cells (Figures 5B,C). In addition, we investigated the potential role of Dnmt3l in undifferentiated and differentiating MSCs using differentially expressed analyses of TEs. The direct correlation between changes in DNA methylation and gene expression of TEs (Figure 5D) did not provide clear evidence of an association in a specific direction. However, our results displayed that the change of global epigenomic landscapes are associated with the altered expression of TEs. These affected TEs or their neighboring genes in Dnmt3l KO MSCs may contribute to the observed osteogenesis defect phenotype in vitro.

Bone Morphology Analysis in vivo

In contrast to the strong osteogenesis defects in vitro, we did not observe obvious differences in trabecular bone phenotype in 5-months-old Dnmt3l KO females (Supplementary Figure 8). In 23-months-old wild type female mice, we can indeed observe the aging induced reductions in trabecular number, connectivity density, percent bone volume, which give higher risk for vertebral fracture, as well as increase in SMI and reductions in bone thickness and bone surface density expected for aging mice (Qiu et al., 2006; Ostertag et al., 2009; Chen et al., 2013). However, we did not observe the expected osteoporosis related symptoms in Dnmt3l KO aging female mice. In contrast, a protective effect was observed in the parameters including connectivity density and trabecular thickness (Supplementary Figure 8).

In an attempt to induce accelerated aging, we introduced high fat diet and time-zone shifting model in 6-month-old Dnmt3l KO male mice and their littermates. The light-dark cycle was shifted 6 h forward each week, for six times in total. Circadian disruption has been associated with metabolism defect and aging-like phenotype (Sato et al., 2017). With the combined jet-lag paradigm and high fat diet treatment, we no longer observed Dnmt3l genotype specific differences in bone properties (Supplementary Figure 9). We concluded that Dnmt3l KO male mice have stronger response to lifestyle intervention-induced trabecular bone mass reduction.

Discussion

Mesenchymal stem/stromal cells hold promises for regenerative medical applications either by cell transplantation directly or by applying their secretory products (Madrigal et al., 2014; Han et al., 2019). In vitro amplification is needed for most of these practices. Various advancements have made to maintain and amplify MSCs in culture condition (Sotiropoulou et al., 2006). However, genotype effect cannot be easily fixed. Apart from infertility, Dnmt3l KO mice does not cause any other obvious defects. It is therefore a surprise for us to see obvious osteogenesis impairment of MSCs derived from Dnmt3l KO in this current study, especially when DNMT3L protein expression cannot be detected in these cell type. Although we cannot exclude the possibility that a small population of MSCs may express Dnmt3l, it is notable that the impaired CFU-F forming and osteogenic activities occurred in global population of MSCs and not the small population of Dnmt3l-expressing MSCs. The current study also identified a potential novel isoform of Dnmt3l transcript starting from exon 9a. Our preliminary RT-qPCR data suggested higher expression from the 3′-end of Dnmt3l in EMSCs than in BM-MSCs, associated with more severe potency defects in Dnmt3l EMSCs. The possibility for this partial transcript of Dnmt3l exon 9a onward being functional could not be excluded at this stage.

However, it is still more plausible that the differences in the cultured MSCs was due to an accumulative cellular memory pre-deposited in DNMT3L expressing progenitor cells. DNMT3L is predominantly expressed in ES cells. We have recently demonstrated that transient ectopic DNMT3L expression in fibroblast is sufficient to introduce long-term effect on epigenome and cell proliferation property, over 40 passages after the ectopic DNMT3L is silenced (Yu et al., 2020), speculating that transiently expressed DNMT3L in ES cells or other progenitor cells, may also have long lasting effects that protect MSCs from premature senescence. The hypothesis of “DNMT3L-mediated epigenetic legacy” (Figure 6) assumes that DNMT3L sustains a normal epigenome or gene expression profile in ES cells. Without DNMT3L expression in the pluripotent stem cells of Dnmt3l KO preimplantation embryos, some of the aberrant epigenetic signatures may be carried over to somatic progenitor cells, and therefore compromising the property of Dnmt3l KO MSCs under in vitro culture and differentiation condition. When comparing DNMT3L binding sites and DNMT3L associated gene expression pattern and methylation status in ES cells, with the DEGs between MSCs from Dnmt3l KO and wild type littermates. There is no strong correlation between DNMT3L bound regions or DNMT3L associated DEGs in ES cells, with Dnmt3l genotype related DEGs in MSCs. Interestingly, DNA methylation at the gene body and promoter regions is most affected in Dnmt3l KD ES cells. Our results showed that gene bodies contain the most overlapped targets, with the effect of DNA methylation potentially carried over from ES cells into MSCs (Figure 4I). Among the 107 DEGs between undifferentiated Dnmt3l KO and wild type MSCs, 24 bears DNMT3L associated DNA methylation on their gene bodies in ES cells. Several of these 24 genes have direct or indirect implications in bone morphogenesis or homeostasis. For example, HIVEP3 (also called Schnurri-3, Shn3) inhibits BMP9-induced osteogenic differentiation from human amniotic MSCs (Tolkachov et al., 2018; Li et al., 2020), and GATA2 is implicated in interfering BMP2-induced SMAD signaling and thus inhibiting osteoblast differentiation demonstrated by an ectopic expression experiment (Tolkachov et al., 2018; Li et al., 2020). FGF18 stimulates rat and mouse MSCs for their osteogenic differentiation (Charoenlarp et al., 2017; Jing et al., 2018). DLX3 serves as a negative regulator of osteogenesis and contributes to bone homeostasis with other transcription factors, such as DLX5 (Levi and Gitton, 2014). Also, Cyp26b1 has effect in osteoblast and Cxcl5 functions for mobilizing HSC away from BM. Although these last two factors may not directly contribute to the phenotype observed in our MSC in vitro osteogenesis impairment phenotype, they provide implications for the miscommunication between the niche MSCs and the DNMT3L expressing HSCs, the progenitor of osteoclast cells. Although there are only limited numbers of genes having correlation between their DNMT3L associated DNA methylation in ESCs and their Dnmt3l genotype associated expression patterns in MSCs, the several highlighted genes may be important enough to contribute in part to the osteogenesis defects observed. In addition, potential DNMT3L associated legacy from progenitor cells or niche interacting defects may go beyond DNA methylation level, and beyond coding genes. The fact that more TEs related sequences are deregulated in cultured Dnmt3l KO MSCs, associated to the deregulation of methylation landscape in Dnmt3l KD ESCs, indirectly support this notion.

FIGURE 6

The serious loss of osteogenesis ability and retained adipogenesis ability of the cultured Dnmt3l KO MSCs out of their endogenous niche resembled the tilted osteogenesis to adipogenesis balance of aging MSCs (Yang et al., 2018). Several mechanisms may simultaneously contribute to MSC aging and senescence related phenomena, including telomere shortening, mitochondrial dysfunction resulting from oxidative stress, the accumulation of DNA damage, and genome wide shifting of the epigenetic modification landscape (Li et al., 2017). Whereas the physiological oxygen level in healthy human bone marrow is approximately 7.5% (Harrison et al., 2002; Jagannathan et al., 2016), conventional cell culture systems use roughly 19% oxygen, which may pose oxidative stress to MSCs, leading to mitochondrial dysfunction and cellular senescence. Indeed, recent studies indicated that hypoxic conditions (1% oxygen) prevent MSCs from entering senescence via the HIF-1α pathway (Tsai et al., 2011). However, osteogenesis may be inhibited under this extreme hypoxic condition (Yang et al., 2011). Interestingly, MSCs consume more oxygen for increased mitochondrial oxidative phosphorylation and exhibit downregulated levels of HIF-1α during osteogenesis, whereas glycolysis activity remains consistent (Shum et al., 2016). It has been shown that HIF-1α requires pyruvate for cancer cells under aerobic conditions (Lu et al., 2002). Accordingly, we observed that the levels of glyceraldehyde 3-phpsphate dehydrogenase (GAPDH), a key player in pyruvate metabolism of glycolysis, drastically fluctuate between undifferentiated and differentiating MSCs (data not shown) despite culturing MSCs under 19% oxygen and with pyruvate. Although the detailed mechanisms for the failure of osteogenic induction in Dnmt3l KO MSCs require much further study, it is possible that metabolic defect could be one of the legacies pre-deposited from DNMT3L-expressing progenitors. Our ongoing study concerning dietary alteration-associated effects in paired Dnmt3l KO and WT littermates also shows similar implications.

Our previous study demonstrated that MEFs derived from Dnmt3l KO embryos enter senescence prematurely. This phenomenon is associated with a reduction of suppressive H3K9me3 and H3K27me3 marks and the premature extrusion of HDAC1 to the cytoplasm (Liao et al., 2015). On the other hand, ectopic expression of DNMT3L in aging MEFs significantly increases the nuclear localization of HDAC1. In addition, exogenous DNMT3L is sufficient to recruit DNMT3A/KAP1/SETDB1/HDAC1 to both newly infected retroviruses and some endogenous retroviruses by inducing H3K9me3 in aging MEFs (Kao et al., 2014). As non-canonical bivalent H3K4me3/H3K9me3 signature was observed on adipogenic determining genes, Cebp-alpha and Ppar-gamma (Matsumura et al., 2015) in MSC cell line, we looked into whether the SETDB1-dependent repressive H3K9me3 marks on those adipogenic master regulator genes accounted for a possible DNMT3L associated epigenetic legacy in MSCs that prevented premature determination of adipogenic fate. Although we have not obtained any conclusion from H3K9me3 ChIP-qPCR experiments, lacks of Cebp-alpha and Ppar-gamma upregulations in undifferentiated Dnmt3l KO MSCs or their osteogenic differentiating derivatives suggest the impaired non-canonical bivalent mark on adiogenic genes is unlikely the primary reason behind the titled osteogenesis to adipogenesis balance observed in Dnmt3l KO MSCs in vitro.

Despite the serious osteogenesis impairment of Dnmt3l KO MSCs in vitro, counterintuitively, does not correlate to strong bone formation defects or osteoporosis phenotype in vivo. One may consider the balance between activities of osteoblasts and osteoclasts in vivo for the bone density phenotype. Since DNMT3L has been shown to be expressed in hematopoietic stem cells (HSCs), the progenitor cells of osteoclasts (Fierro et al., 2017), DNMT3L deficiency in HSCs may potentially contribute to reduced osteoclast function and hence the increased bone density observed in Dnmt3l KO mice. MSCs are also determined as one of the cell types building HSC niche (Ehninger and Trumpp, 2011), suggesting potential intercommunication between these two cell lineages. Interestingly, expression of the Cxcl5 gene that encodes a chemokine to induce mobilization of HSC from BM (Yoon et al., 2012) was elevated in Dnmt3l KO MSCs, whereas expression of the factors for HSC maintenance, Angpt1 and Cxcl12 were significantly decreased (Supplementary Table 3). Altogether, in Dnmt3l KO mice, the defective osteogenesis of MSCs as disclosed in our in vitro study, combined with the proposed inefficiency of bone breakdown, would give rise to only subtle, if any, phenotypes in vivo concerning the bone homeostasis. The in-depth study of DNMT3L function in HSCs and their derivatives would be another important episode for understanding the function and properties of this important epigenomic modulator DNMT3L beyond its highly expressing pluripotent stem cells and germ cells.

We demonstrated in this study that long term effects on cell fate and property can be induced by transient expression or deficiency of epigenetic modifiers and cofactors. These effects may even be amplified from spreading of the epigenetic alternations through excessive cell proliferations in vitro or in vivo. This is particularly important for the cell-based regenerative medicine when transient exposure of stem cells under aberrant condition would have long lasting effect that may not be perceived before transplantation, but conveying long-term risk thereafter.

Materials and Methods

Animals

Due to the infertility phenotype for both male and female Dnmt3l KO mice, the Dnmt3l^+/tmEnls mutant allele was maintained as heterozygous colony. Dnmt3l KO mice and Dnmt3l^+/+ littermate controls were generated from intercrossing mice with the heterozygous Dnmt3l^+/tmEnls mutant allele (Hata et al., 2002) in a C57BL/6 genetic background. In average, one Dnmt3l^–/– KO mice and one Dnmt3l^+/+ littermate control from one gender were generated from one litter. Female mice (8–16 weeks old) were euthanized to harvest MSCs for experiments. The care of mice and the experimental procedures involving animals were approved by the Institutional Animal Care and Use Committee (IACUC) of National Taiwan University (approval number NTU-104-EL-00031).

Isolation of Mouse Bone Marrow MSCs

For each experiments, we collected bone marrow MSCs from three Dnmt3l KO mice of the same gender and age, as well as the counterparts from the three Dnmt3l^+/+ littermate controls, to prepare mixed MSCs from each genotype. The procedure used to isolate mouse bone marrow MSCs was previously described (Cheng et al., 2012, 2014; Lee et al., 2017). Briefly, the femur and tibia were separated from the mouse after euthanasia by quickly removing the muscle and connective tissue with forceps and scissors. A 23-gauge needle containing 12 mL of culture medium [αMEM (Gibco, Grand Island, NY, United States) supplemented with 20% FBS (HyClone, Logan, UT, United States) and antibiotics (100 U/mL penicillin and 100 μg/mL streptomycin; Gibco)] was used to flush out femur bone marrow. A similar procedure was performed with a 26-gauge needle and 6 mL of culture medium to acquire tibia bone marrow. The bone marrow and minced femur and tibia were cultured in 100-mm Petri dishes (TPP, Trasadingen, CH, Switzerland). After 9 days of culturing with medium changes every 3 days to deplete non-adherent cells, adherent cells were detached with 0.25% trypsin/EDTA (Gibco) for 5 min in an incubator. Cells were then transferred to a 70-μm nylon cell strainer (Falcon, Corning, NY, United States) to remove cell clumps. MSCs were isolated using a MACS system (Miltenyi Biotec, Teterow, DE, Germany) with CD11b (Miltenyi Biotec) and CD45 (Miltenyi Biotec) microbead antibody double-negative selection according to the manufacturer’s protocol to avoid hematopoietic lineage contamination. Briefly, 10⁷ cells were resuspended in 10 μL of CD11b antibody, 10 μL of CD45 antibody and 180 μL of MACS^® running buffer (Miltenyi Biotec), and the cells were left to react for 15 min at 4°C. One microliter of running buffer was added to stop the reaction, and the cells were centrifuged for 5 min at 1200 rpm. Subsequently, the supernatant was removed, and another 1 mL of running buffer was added to resuspend the cells. A MACS^® LD column (Miltenyi Biotec) was inserted into a MidiMACS^TM separator (Miltenyi Biotec), and a collection tube was placed under the column to harvest MSCs. The column was rinsed with 2 mL of running buffer. Afterward, the cells were transferred into the column. When the column was about to dry out, 1 mL of running buffer was added to wash out unlabeled cells (i.e., non-hematopoietic cells). The procedure was repeated, and enriched MSCs were harvested from the effluent. After 5 min of centrifugation at 1200 rpm, running buffer was removed and replaced with culture medium. MSCs were seeded at a density of 5 × 10⁴ cells per cm².

Isolation of an MSC Subtype: Epiphysis-Derived Mesenchymal Stem/Stromal Cells (EMSCs)

The procedure used to harvest the EMSC subtype was based on our previous report (Cheng et al., 2012). We isolated the femurs and tibias of mice and thoroughly removed the adherent soft tissues. The epiphyses were then removed using shears and rinsing out of any possible bone marrow fluids with culture medium. The cleaned intact epiphyses (without mincing) were then transferred to 22.1 cm² culture dishes (TPP) and cultured in complete medium composed of MEM alpha (Sigma-Aldrich) supplemented with 20% FBS (HyClone), 2 mM L-glutamine (100 U/ml penicillin, and 100 μg/ml streptomycin; Invitrogen). The cells were incubated under a humidified atmosphere containing 95% air and 5% CO₂ at 37°C. The non-adherent cells were removed by changing the medium every 3 days. After 11 days, the primary culture reached approximately 70% confluence, and the cells were detached with 0.25% trypsin/0.1 mM EDTA (trypsin/EDTA; Invitrogen) for 1 min at 37°C. Subsequently, the reaction was stopped by adding complete medium. The cells that could not be lifted within 1 min were discarded. The detached cells were centrifuged at 400 × g for 5 min, resuspended in the complete medium, and subjected to negative selection using the CD11b and CD45 antibodies as described in the BM-MSC isolation section above.

Colony-Forming Unit Assay

One or two thousand cells from passage 3 were seeded into 6-well plates (TPP). After 9 days of culture, cells were fixed with methanol (J. T. Baker, Center Valley, PA, United States) and visualized by Giemsa stain (Gibco). Colonies greater than 2 mm in diameter were defined as valid.

Cell Growth Assay

The PDT was determined to assess the growth of MSCs. In brief, Dnmt3l KO MSCs and wild type MSCs from passage 3 were seeded for 33 wells with a density of 2 × 10⁴ cells per well in 12-well plate. These cells were allowed to adhere for 24 h at 37°C under an atmosphere with 95% air and 5% CO₂. Cells from three wells of each genotype were isolated and evaluated by trypan blue staining each day for 11 days to calculate the number of living cells.

Osteogenic Differentiation

Osteogenic induction medium was composed of 10% FBS-αMEM supplemented with 0.1 μM dexamethasone (Sigma-Aldrich, St. Louis, MO, United States), 10 mM β-glycerophosphate (Sigma-Aldrich) and 50 μM L-ascorbic acid 2-phosphate (Sigma-Aldrich) (Jaiswal et al., 1997). For each sample, 5 × 10⁴ cells per cm² were seeded into a 6-well plate. After reaching confluence in 2 days, the cells were incubated with osteogenic induction medium. The medium was changed every 3 days, and analyses were performed after 14 days of induction. To assess calcium precipitation, the induction medium was removed, and the cells were washed once with PBS. Next, the cells were fixed with 10% formaldehyde (Sigma-Aldrich) for 10 min with gentle shaking, and after another wash with PBS, the cells were stained with 500 μL of 2% ARS (pH 4.1–4.3) (Sigma-Aldrich) for 15 min. After extensive PBS irrigation, the precipitate was dissolved in 500 μL of 10% cetylpyridinium chloride (Sigma-Aldrich) in 8 mM Na₂PO₄ (Sigma-Aldrich) and 1.5 mM KH₂PO₄ (Sigma-Aldrich). The absorbance at 550 nm was recorded, and quantification of calcium was determined using an optimal ARS standard curve. The optimal dilution was performed to ensure that the sample absorbance fell within the standard curve.

RNA Extraction, Library Preparation and Sequencing

Total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, United States) and a RNeasy mini kit (Qiagen, Hilden, DE, Germany) following the manufacturer’s instructions. RNA was treated with DNase 1 (Qiagen) to prevent DNA contamination. The quality of total RNA was checked using an Agilent Bioanalyzer 2100 with an Agilent RNA 6000 Nano Kit, and the RIN of the total RNA samples were all greater than 8. We used 2.5 μg of total RNA as input for strand-specific RNA library preparation following the removal of rRNA using a Ribosomal Illumina Ribo-Zero Gold rRNA Removal Kit (Human/Mouse/Rat) (Illumina, San Diego, CA, United States). The rRNA-depleted RNA libraries were constructed using an Illumina TruSeq Stranded mRNA Sample Prep kit following the user guide at the fragment RNA step. The libraries were sequenced using an Illumina HiSeq 2500 system with Rapid run mode (50-bp single-end reads). Each sample was sequenced to obtain approximately 70.7 M raw reads.

RNA-Seq Analysis

The mRNA reads were strand-specifically aligned to the mouse reference genome (mm10) using Tophat2 (Kim et al., 2013; Trapnell et al., 2013), with the gene annotation downloaded from iGenome. The reads per gene counts were generated using Cuffdiff (Trapnell et al., 2013), which was also used to identify DEGs between samples using a ≥2-fold change in expression and a FDR less than 5%. The transposon annotation was downloaded from UCSC (Karolchik et al., 2004). The expression counts of transposons were generated using multiBamCov from BEDtools (Quinlan and Hall, 2010), and a transposon was considered to be expressed if it was expressed in at least one or more samples. The normalization of sequencing depth across samples and differential transposon expression analysis was performed using edgeR (Robinson et al., 2010) with a ≥2-fold change in expression and FDR 5%. The minimum expressed TE replicate had to have a minimum of 10 reads among all samples after normalizing for sequencing depth across samples.

Weighted Gene Co-expression Network Analysis (WGCNA)

Weighted gene co-expression network analysis was performed using the R package for transcriptome data following the standard method described on the authors’ website (Langfelder and Horvath, 2008). This unsupervised and unbiased analysis identifies distinct co-expression gene modules by clustering transcripts with similar expression patterns across samples. Transcriptionally variable genes with an RPKM of ≥0.5 were selected for the analysis. These genes were hierarchically clustered based on a dissimilarity measure of topological overlap matrix. The resulting dendrogram was used for module detection (minimum size 100; height cutoff 0.23). Gene modules were labeled in unique colors, with unassigned genes labeled in gray.

Bioinformatic Analysis on ChIP-Seq and MeDIP-Seq Datasets

The ChIP-Seq of DNMT3L from mouse ES cells was from GSE49178 (Neri et al., 2013), and the DNA methylation data (MeDIP-Seq) from mouse Dnmt3l KD ES cells was from GSE44642 (Neri et al., 2013). All reads were aligned to the mouse reference genome (mm10) using Bowtie2 (Langmead and Salzberg, 2012), with peaks called using MACS2 (Feng et al., 2012).

Western Blot

Cells were homogenized and lysed in RIPA lysis buffer (EMD Millipore, Burlington, MA, United States) supplemented with 1 mM phenylmethylsulfonyl fluoride (Sigma-Aldrich) and protease inhibitor cocktail (Sigma-Aldrich). Lysate (5–10 μg per lane) were separated by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS–PAGE) and transferred to PVDF membranes (EMD Millipore). The membranes were blocked with 5% BSA (Sigma-Aldrich) in 1x PBST for 1 h at room temperature or 4°C overnight. The incubation with primary antibody anti-DNMT3L (E1Y7Q, Cell signaling technology, 1:1000 dilution), and anti-GAPDH (ab181602, Abcam, 1:1000 dilution) for 1 h at room temperature with gentle shaking. Membranes were incubated with a horseradish peroxidase - conjugated secondary antibody (1:5000 dilution) for 1 h under room temperature. The proteins were detected using a chemiluminescent reagent (Millipore), and the images were performed by GeneGnome XRQ-chemiluminescence image system (SYNGENE).

Assessment of Trabecular Bone Parameters

The entire legs of mice with different age, gender and Dnmt3l genotypes were collected. After carefully removing the skin and most of the muscle, the intact femur, tibia, fibula with intact patella and joint were stored in 70% EtOH, and subject to 3D micro-CT imaging. The area of distal femur 0.4 to 1.9 mm from growth plate were used to analyze bone density and other trabecular parameters, with source voltage of 70 kV, source current of 57 μA, exposure for 1669 ms with AI 0.5 mm filter and image pixel size of 9 μm. Among the parameters, the elevated trabecular separation, reduced trabecular number, connectivity density and percent bone volume are associated with vertebral fracture (Qiu et al., 2006; Ostertag et al., 2009). On the other hand, increased SMI, decreased trabecular thickness and bone surface density are linked to with aging associated decline (Stauber and Muller, 2006; Chen et al., 2013).

References

• 1

  AapolaU.KawasakiK.ScottH. S.OllilaJ.VihinenM.HeinoM.et al (2000). Isolation and initial characterization of a novel zinc finger gene, DNMT3L, on 21q22.3, related to the cytosine-5-methyltransferase 3 gene family.Genomics65293–298. 10.1006/geno.2000.6168

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AbdiR.FiorinaP.AdraC. N.AtkinsonM.SayeghM. H. (2008). Immunomodulation by mesenchymal stem cells: a potential therapeutic strategy for type 1 diabetes.Diabetes571759–1767. 10.2337/db08-0180

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  AtlasiY.StunnenbergH. G. (2017). The interplay of epigenetic marks during stem cell differentiation and development.Nat. Rev. Genet.18643–658. 10.1038/nrg.2017.57

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BalliniA.ScaccoS.ColettiD.PluchinoS.TatulloM. (2017). Mesenchymal Stem Cells as Promoters, Enhancers, and Playmakers of the Translational Regenerative Medicine.Stem Cells Int.20173292810. 10.1155/2017/3292810

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BenderC. M.GonzalgoM. L.GonzalesF. A.NguyenC. T.RobertsonK. D.JonesP. A. (1999). Roles of cell division and gene transcription in the methylation of CpG islands.Mol. Cell Biol.196690–6698.

  • Google Scholar
  • View reference in article

• 6

  Berebichez-FridmanR.Montero-OlveraP. R. (2018). Sources and Clinical Applications of Mesenchymal Stem Cells: State-of-the-art review.Sultan Qaboos Univ. Med. J.18e264–e277. 10.18295/squmj.2018.18.03.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  BiebackK.HeckerA.KocaomerA.LannertH.SchallmoserK.StrunkD.et al (2009). Human alternatives to fetal bovine serum for the expansion of mesenchymal stromal cells from bone marrow.Stem Cells272331–2341. 10.1002/stem.139

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  Bloushtain-QimronN.YaoJ.ShipitsinM.MaruyamaR.PolyakK. (2009). Epigenetic patterns of embryonic and adult stem cells.Cell Cycle8809–817. 10.4161/cc.8.6.7938

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  BonabM. M.AlimoghaddamK.TalebianF.GhaffariS. H.GhavamzadehA.NikbinB. (2006). Aging of mesenchymal stem cell in vitro.BMC Cell Biol.7:14. 10.1186/1471-2121-7-14

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  BonnetN.ConwayS. J.FerrariS. L. (2012). Regulation of beta catenin signaling and parathyroid hormone anabolic effects in bone by the matricellular protein periostin.Proc. Natl. Acad. Sci. U S A10915048–15053. 10.1073/pnas.1203085109

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  Bourc’hisD.BestorT. H. (2004). Meiotic catastrophe and retrotransposon reactivation in male germ cells lacking Dnmt3L.Nature43196–99. 10.1038/nature02886

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  Bourc’hisD.XuG. L.LinC. S.BollmanB.BestorT. H. (2001). Dnmt3L and the establishment of maternal genomic imprints.Science2942536–2539. 10.1126/science.1065848

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  CaplanA. I. (1991). Mesenchymal stem cells.J. Orthop. Res.9641–650. 10.1002/jor.1100090504

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  CharoenlarpP.RajendranA. K.IsekiS. (2017). Role of fibroblast growth factors in bone regeneration.Inflamm. Regen.37:10. 10.1186/s41232-017-0043-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  ChedinF.LieberM. R.HsiehC. L. (2002). The DNA methyltransferase-like protein DNMT3L stimulates de novo methylation by Dnmt3a.Proc. Natl. Acad. Sci. U S A9916916–16921. 10.1073/pnas.262443999

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  ChenH.ZhouX.FujitaH.OnozukaM.KuboK. Y. (2013). Age-related changes in trabecular and cortical bone microstructure.Int. J. Endocrinol.2013:213234. 10.1155/2013/213234

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  ChenS. Y.LinM. C.TsaiJ. S.HeP. L.LuoW. T.HerschmanH.et al (2019). EP4 Antagonist-Elicited Extracellular Vesicles from Mesenchymal Stem Cells Rescue Cognition/Learning Deficiencies by Restoring Brain Cellular Functions.Stem Cells Transl. Med.8707–723. 10.1002/sctm.18-0284

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  ChengC. C.LeeY. H.LinS. P.HuangfuW. C.LiuI. H. (2014). Cell-autonomous heparanase modulates self-renewal and migration in bone marrow-derived mesenchymal stem cells.J. Biomed. Sci.21:21. 10.1186/1423-0127-21-21

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  ChengC. C.LianW. S.HsiaoF. S.ILiuH.LinS. P.LeeY. H.et al (2012). Isolation and characterization of novel murine epiphysis derived mesenchymal stem cells.PLoS One7:e36085. 10.1371/journal.pone.0036085

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  CollinsC. A.IOlsenP. S.ZammitL.HeslopA.PetrieT. A.Partridgeet al (2005). Stem cell function, self-renewal, and behavioral heterogeneity of cells from the adult muscle satellite cell niche.Cell122289–301. 10.1016/j.cell.2005.05.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  DigirolamoC. M.StokesD.ColterD.PhinneyD. G.ClassR.ProckopD. J. (1999). Propagation and senescence of human marrow stromal cells in culture: a simple colony-forming assay identifies samples with the greatest potential to propagate and differentiate.Br. J. Haematol.107275–281. 10.1046/j.1365-2141.1999.01715.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  EhningerA.TrumppA. (2011). The bone marrow stem cell niche grows up: mesenchymal stem cells and macrophages move in.J. Exp. Med.208421–428. 10.1084/jem.20110132

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  EnglerA. J.SenS.SweeneyH. L.DischerD. E. (2006). Matrix elasticity directs stem cell lineage specification.Cell126677–689. 10.1016/j.cell.2006.06.044

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  EricesA.CongetP.MinguellJ. J. (2000). Mesenchymal progenitor cells in human umbilical cord blood.Br. J. Haematol.109235–242.

  • Google Scholar
  • View reference in article

• 25

  FengJ.LiuT.QinB.ZhangY.LiuX. S. (2012). Identifying ChIP-seq enrichment using MACS.Nat. Protoc.71728–1740. 10.1038/nprot.2012.101

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  FierroF. A.NoltaJ. A.AdamopoulosI. E. (2017). Concise Review: Stem Cells in Osteoimmunology.Stem Cells351461–1467. 10.1002/stem.2625

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  GaoF.ChiuS. M.MotanD. A.ZhangZ.ChenL.JiH. L.et al (2016). Mesenchymal stem cells and immunomodulation: current status and future prospects.Cell Death Dis.7e2062. 10.1038/cddis.2015.327

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  GuenatriM.DuffieR.IranzoJ.FauqueP.Bourc’hisD. (2013). Plasticity in Dnmt3L-dependent and -independent modes of de novo methylation in the developing mouse embryo.Development140562–572. 10.1242/dev.089268

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  HanY.LiX.ZhangY.HanY.ChangF.DingJ. (2019). Mesenchymal Stem Cells for Regenerative Medicine.Cells8:8080886. 10.3390/cells8080886

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  HarrisonJ. S.RameshwarP.ChangV.BandariP. (2002). Oxygen saturation in the bone marrow of healthy volunteers.Blood99:394.

  • Google Scholar
  • View reference in article

• 31

  HashimotoH.VertinoP. M.ChengX. (2010). Molecular coupling of DNA methylation and histone methylation.Epigenomics2657–669. 10.2217/epi.10.44

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  HataK.KusumiM.YokomineT.LiE.SasakiH. (2006). Meiotic and epigenetic aberrations in Dnmt3L-deficient male germ cells.Mol. Reprod Dev.73116–122. 10.1002/mrd.20387

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  HataK.OkanoM.LeiH.LiE. (2002). Dnmt3L cooperates with the Dnmt3 family of de novo DNA methyltransferases to establish maternal imprints in mice.Development1291983–1993.

  • Google Scholar
  • View reference in article

• 34

  In t AnkerP. S.ScherjonS. A.NoortW. A.ClaasF. H.WillemzeR.et al (2003). Amniotic fluid as a novel source of mesenchymal stem cells for therapeutic transplantation.Blood1021548–1549. 10.1182/blood-2003-04-1291

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  InfanteA.RodriguezC. I. (2018). Osteogenesis and aging: lessons from mesenchymal stem cells.Stem Cell Res. Ther.9:244. 10.1186/s13287-018-0995-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  JagannathanL.CuddapahS.CostaM. (2016). Oxidative stress under ambient and physiological oxygen tension in tissue culture.Curr. Pharmacol. Rep.264–72. 10.1007/s40495-016-0050-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  JaiswalN.HaynesworthS. E.CaplanA. I.BruderS. P. (1997). Osteogenic differentiation of purified, culture-expanded human mesenchymal stem cells in vitro.J. Cell Biochem.64295–312.

  • Google Scholar
  • View reference in article

• 38

  JingH.SuX.GaoB.ShuaiY.ChenJ.DengZ.et al (2018). Epigenetic inhibition of Wnt pathway suppresses osteogenic differentiation of BMSCs during osteoporosis.Cell Death Dis.9:176. 10.1038/s41419-017-0231-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  KaoT. H.LiaoH. F.WolfD.TaiK. Y.ChuangC. Y.LeeH. S.et al (2014). Ectopic DNMT3L triggers assembly of a repressive complex for retroviral silencing in somatic cells.J. Virol.8810680–10695. 10.1128/JVI.01176-14

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  KarolchikD.HinrichsA. S.FureyT. S.RoskinK. M.SugnetC. W.HausslerD.et al (2004). The UCSC Table Browser data retrieval tool.Nucleic Acids Res.32D493–D496. 10.1093/nar/gkh103

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  KimD.PerteaG.TrapnellC.PimentelH.KelleyR.SalzbergS. L. (2013). TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions.Genome Biol.14:R36. 10.1186/gb-2013-14-4-r36

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  LangfelderP.HorvathS. (2008). WGCNA: an R package for weighted correlation network analysis.BMC Bioinformatics9:559. 10.1186/1471-2105-9-559

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  LangmeadB.SalzbergS. L. (2012). Fast gapped-read alignment with Bowtie 2.Nat. Methods9357–359. 10.1038/nmeth.1923

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  LeeM. K.LinS. P.HuangFuW. C.YangD. S.LiuI. H. (2017). Endothelial-derived extracellular matrix ameliorate the stemness deprivation during ex vivo expansion of mouse bone marrow-derived mesenchymal stem cells.PLoS One12:e0184111. 10.1371/journal.pone.0184111

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  LeviG.GittonY. (2014). Dlx genes and the maintenance of bone homeostasis and skeletal integrity.Cell Death Differ.211345–1346. 10.1038/cdd.2014.94

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  LiY.LiuZ.TangY.FengW.ZhaoC.LiaoJ.et al (2020). Schnurri-3 regulates BMP9-induced osteogenic differentiation and angiogenesis of human amniotic mesenchymal stem cells through Runx2 and VEGF.Cell Death Dis.11:72. 10.1038/s41419-020-2279-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  LiY.WuQ.WangY.LiL.BuH.BaoJ. (2017). Senescence of mesenchymal stem cells (Review).Int. J. Mol. Med.39775–782. 10.3892/ijmm.2017.2912

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  LiaoH. F.ChenW. S.ChenY. H.KaoT. H.TsengY. T.LeeC. Y.et al (2014). DNMT3L promotes quiescence in postnatal spermatogonial progenitor cells.Development1412402–2413. 10.1242/dev.105130

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  LiaoH. F.MoC. F.WuS. C.ChengD. H.YuC. Y.ChangK. W.et al (2015). Dnmt3l-knockout donor cells improve somatic cell nuclear transfer reprogramming efficiency.Reproduction150245–256. 10.1530/REP-15-0031

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  LinC. Y.LinH. H.TsaiM. H.LinS. P.ChenM. H. (2011). Zinc chloride for odontogenesis of dental pulp stem cells via metallothionein up-regulation.J. Endod.37211–216. 10.1016/j.joen.2010.11.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  LuH.ForbesR. A.VermaA. (2002). Hypoxia-inducible factor 1 activation by aerobic glycolysis implicates the Warburg effect in carcinogenesis.J. Biol. Chem.27723111–23115. 10.1074/jbc.M202487200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  LunyakV. V.RosenfeldM. G. (2008). Epigenetic regulation of stem cell fate.Hum. Mol. Genet.17R28–R36. 10.1093/hmg/ddn149

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  MadrigalM.RaoK. S.RiordanN. H. (2014). A review of therapeutic effects of mesenchymal stem cell secretions and induction of secretory modification by different culture methods.J. Transl. Med.12:260. 10.1186/s12967-014-0260-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  MarrazzoP.PaduanoF.PalmieriF.MarrelliM.TatulloM. (2016). Highly Efficient In Vitro Reparative Behaviour of Dental Pulp Stem Cells Cultured with Standardised Platelet Lysate Supplementation.Stem Cells Int.2016:7230987. 10.1155/2016/7230987

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  MarrelliM.PujiaA.PalmieriF.GattoR.FalisiG.GargariM.et al (2016). Innovative approach for the in vitro research on biomedical scaffolds designed and customized with CAD-CAM technology.Int. J. Immunopathol. Pharmacol.29778–783. 10.1177/0394632016646121

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  MatsumuraY.NakakiR.InagakiT.YoshidaA.KanoY.KimuraH.et al (2015). H3K4/H3K9me3 Bivalent Chromatin Domains Targeted by Lineage-Specific DNA Methylation Pauses Adipocyte Differentiation.Mol. Cell60584–596. 10.1016/j.molcel.2015.10.025

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  MeinelL.KarageorgiouV.FajardoR.SnyderB.Shinde-PatilV.ZichnerL.et al (2004). Bone tissue engineering using human mesenchymal stem cells: effects of scaffold material and medium flow.Ann. Biomed. Eng.32112–122. 10.1023/b:abme.0000007796.48329.b4

  • CrossRef
  • Google Scholar

• 58

  MerleB.GarneroP. (2012). The multiple facets of periostin in bone metabolism.Osteoporos Int.231199–1212. 10.1007/s00198-011-1892-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  MiaoZ.JinJ.ChenL.ZhuJ.HuangW.ZhaoJ.et al (2006). Isolation of mesenchymal stem cells from human placenta: comparison with human bone marrow mesenchymal stem cells.Cell Biol. Int.30681–687. 10.1016/j.cellbi.2006.03.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  NeriF.KrepelovaA.IncarnatoD.MaldottiM.ParlatoC.GalvagniF.et al (2013). Dnmt3L antagonizes DNA methylation at bivalent promoters and favors DNA methylation at gene bodies in ESCs.Cell155121–134. 10.1016/j.cell.2013.08.056

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  OoiS. K.QiuC.BernsteinE.LiK.JiaD.YangZ.et al (2007). DNMT3L connects unmethylated lysine 4 of histone H3 to de novo methylation of DNA.Nature448714–717. 10.1038/nature05987

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  OstertagA.Cohen-SolalM.AudranM.LegrandE.MartyC.ChappardD.et al (2009). Vertebral fractures are associated with increased cortical porosity in iliac crest bone biopsy of men with idiopathic osteoporosis.Bone44413–417. 10.1016/j.bone.2008.11.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  OtaniJ.NankumoT.AritaK.InamotoS.AriyoshiM.ShirakawaM. (2009). Structural basis for recognition of H3K4 methylation status by the DNA methyltransferase 3A ATRX-DNMT3-DNMT3L domain.EMBO Rep.101235–1241. 10.1038/embor.2009.218

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  PittengerM. F.MackayA. M.BeckS. C.JaiswalR. K.DouglasR.MoscaJ. D.et al (1999). Multilineage potential of adult human mesenchymal stem cells.Science284143–147.

  • Google Scholar

• 65

  QiuS.RaoD. S.PalnitkarS.ParfittA. M. (2006). Independent and combined contributions of cancellous and cortical bone deficits to vertebral fracture risk in postmenopausal women.J. Bone Miner Res.211791–1796. 10.1359/jbmr.060801

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  QuinlanA. R.HallI. M. (2010). BEDTools: a flexible suite of utilities for comparing genomic features.Bioinformatics26841–842. 10.1093/bioinformatics/btq033

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67

  ReikW. (2007). Stability and flexibility of epigenetic gene regulation in mammalian development.Nature447425–432. 10.1038/nature05918

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  RobinsonM. D.McCarthyD. J.SmythG. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.Bioinformatics26139–140. 10.1093/bioinformatics/btp616

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  SatoS.SolanasG.PeixotoF. O.BeeL.SymeonidiA.SchmidtM. S.et al (2017). Circadian Reprogramming in the Liver Identifies Metabolic Pathways of Aging.Cell66:e611. 10.1016/j.cell.2017.07.042

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  ShovlinT. C.Bourc’hisD.La SalleS.O’DohertyA.TraslerJ. M.BestorT. H.et al (2007). Sex-specific promoters regulate Dnmt3L expression in mouse germ cells.Hum. Reprod22457–467. 10.1093/humrep/del379

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71

  ShumL. C.WhiteN. S.MillsB. N.BentleyK. L.EliseevR. A. (2016). Energy Metabolism in Mesenchymal Stem Cells During Osteogenic Differentiation.Stem Cells Dev.25114–122. 10.1089/scd.2015.0193

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  SotiropoulouP. A.PerezS. A.SalagianniM.BaxevanisC. N.PapamichailM. (2006). Characterization of the optimal culture conditions for clinical scale production of human mesenchymal stem cells.Stem Cells24462–471. 10.1634/stemcells.2004-0331

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  SpagnuoloG.CodispotiB.MarrelliM.RengoC.RengoS.TatulloM. (2018). Commitment of Oral-Derived Stem Cells in Dental and Maxillofacial Applications.Dent J.6:6040072. 10.3390/dj6040072

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  StauberM.MullerR. (2006). Age-related changes in trabecular bone microstructures: global and local morphometry.Osteoporos Int.17616–626. 10.1007/s00198-005-0025-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75

  StockM.SchaferH.FliegaufM.OttoF. (2004). Identification of novel genes of the bone-specific transcription factor Runx2.J. Bone Miner Res.19959–972. 10.1359/jbmr.2004.19.6.959

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  TolkachovA.FischerC.AmbrosiT. H.BotheM.HanC. T.MuenznerM.et al (2018). Loss of the Hematopoietic Stem Cell Factor GATA2 in the Osteogenic Lineage Impairs Trabecularization and Mechanical Strength of Bone.Mol. Cell Biol.38599–517. 10.1128/MCB.00599-17

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 77

  TrapnellC.HendricksonD. G.SauvageauM.GoffL.RinnJ. L.PachterL. (2013). Differential analysis of gene regulation at transcript resolution with RNA-seq.Nat. Biotechnol.3146–53. 10.1038/nbt.2450

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78

  TsaiC. C.ChenY. J.YewT. L.ChenL. L.WangJ. Y.ChiuC. H.et al (2011). Hypoxia inhibits senescence and maintains mesenchymal stem cell properties through down-regulation of E2A-p21 by HIF-TWIST.Blood117459–469. 10.1182/blood-2010-05-287508

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79

  WangJ. P.LiaoY. T.WuS. H.ChiangE. R.HsuS. H.TsengT. C.et al (2020). Mesenchymal stem cells from a hypoxic culture improve nerve regeneration.J. Tissue Eng. Regen. Med.141804–1814. 10.1002/term.3136

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 80

  WebsterK. E.O’BryanM. K.FletcherS.CrewtherP. E.AapolaU.CraigJ.et al (2005). Meiotic and epigenetic defects in Dnmt3L-knockout mouse spermatogenesis.Proc. Natl. Acad. Sci. U S A1024068–4073. 10.1073/pnas.0500702102

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 81

  YangD. C.YangM. H.TsaiC. C.HuangT. F.ChenY. H.HungS. C. (2011). Hypoxia inhibits osteogenesis in human mesenchymal stem cells through direct regulation of RUNX2 by TWIST.PLoS One6:e23965. 10.1371/journal.pone.0023965

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 82

  YangY. K. (2018). Aging of mesenchymal stem cells: Implication in regenerative medicine.Regen Ther.9120–122. 10.1016/j.reth.2018.09.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83

  YangY. K.OgandoC. R.Wang SeeC.ChangT. Y.BarabinoG. A. (2018). Changes in phenotype and differentiation potential of human mesenchymal stem cells aging in vitro.Stem Cell Res. Ther.9:131. 10.1186/s13287-018-0876-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 84

  YoonK. A.ChoH. S.ShinH. I.ChoJ. Y. (2012). Differential regulation of CXCL5 by FGF2 in osteoblastic and endothelial niche cells supports hematopoietic stem cell migration.Stem Cells Dev.213391–3402. 10.1089/scd.2012.0128

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85

  YoungH. E.SteeleT. A.BrayR. A.HudsonJ.FloydJ. A.HawkinsK.et al (2001). Human reserve pluripotent mesenchymal stem cells are present in the connective tissues of skeletal muscle and dermis derived from fetal, adult, and geriatric donors.Anat. Rec.26451–62.

  • Google Scholar

• 86

  YuY. C.HuiT. Z.KaoT. H.LiaoH. F.YangC. Y.HouC. C.et al (2020). Transient DNMT3L Expression Reinforces Chromatin Surveillance to Halt Senescence Progression in Mouse Embryonic Fibroblast.Front. Cell Dev. Biol.8:103. 10.3389/fcell.2020.00103

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 87

  ZhangF.LuoK.RongZ.WangZ.LuoF.ZhangZ.et al (2017). Periostin Upregulates Wnt/beta-Catenin Signaling to Promote the Osteogenesis of CTLA4-Modified Human Bone Marrow-Mesenchymal Stem Cells.Sci. Rep.7:41634. 10.1038/srep41634

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 88

  ZhangY.JurkowskaR.SoeroesS.RajaveluA.DhayalanA.BockI.et al (2010). Chromatin methylation activity of Dnmt3a and Dnmt3a/3L is guided by interaction of the ADD domain with the histone H3 tail.Nucleic Acids Res.384246–4253. 10.1093/nar/gkq147

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 89

  ZhouS.GreenbergerJ. S.EpperlyM. W.GoffJ. P.AdlerC.LeboffM. S.et al (2008). Age-related intrinsic changes in human bone-marrow-derived mesenchymal stem cells and their differentiation to osteoblasts.Aging Cell7335–343. 10.1111/j.1474-9726.2008.00377.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 90

  ZukP. A.ZhuM.AshjianP.De UgarteD. A.HuangJ. I.MizunoH.et al (2002). Human adipose tissue is a source of multipotent stem cells.Mol. Biol. Cell134279–4295. 10.1091/mbc.e02-02-0105

  • Pubmed Abstract
  • CrossRef
  • Google Scholar