Stem Cell Reports Ar ticle
DNA Methyltransferases Modulate Hepatogenic Lineage Plasticity of Mesenchymal Stromal Cells Chien-Wei Lee,1,14 Wei-Chih Huang,3 Hsien-Da Huang,3,4,5,6 Yi-Hsiang Huang,7 Jennifer H. Ho,8 Muh-Hwa Yang,9,10,11,12,13 Vincent W. Yang,17 and Oscar K. Lee2,14,15,16,* 1Program
in Molecular Medicine, National Yang-Ming University and Academia Sinica, Taipei 11221, Taiwan City Hospital, Taipei 10341, Taiwan 3Institute of Bioinformatics and Systems Biology 4Department of Biological Science and Technology 5Center for Bioinformatics Research National Chiao Tung University, HsinChu 30010, Taiwan 6Department of Biomedical Science and Environmental Biology, Kaohsiung Medical University, Kaohsiung 80708, Taiwan 7Division of Gastroenterology, Department of Medicine, Taipei Veterans General Hospital, Taipei 11221, Taiwan 8Center for Stem Cell Research, Taipei Medical University-Wan Fang Hospital, Taipei 11031, Taiwan 9Institute of Biotechnology in Medicine 10Genome Research Center 11Immunity and Inflammation Research Center National Yang-Ming University, Taipei 11221, Taiwan 12Division of Hematology and Oncology, Department of Medicine, Taipei Veterans General Hospital, Taipei 11217, Taiwan 13Genomics Research Center, Academia Sinica, Taipei 11529, Taiwan 14Institute of Clinical Medicine, National Yang-Ming University 15Stem Cell Research Center, National Yang-Ming University 16Department of Medical Research, Taipei Veterans General Hospital Taipei 11221, Taiwan 17Department of Medicine, Stony Brook University, Stony Brook, NY 11794, USA *Correspondence:
[email protected] http://dx.doi.org/10.1016/j.stemcr.2017.05.008 2Taipei
SUMMARY The irreversibility of developmental processes in mammalian cells has been challenged by rising evidence that de-differentiation of hepatocytes occurs in adult liver. However, whether reversibility exists in mesenchymal stromal cell (MSC)-derived hepatocytes (dHeps) remains elusive. In this study, we find that hepatogenic differentiation (HD) of MSCs is a reversible process and is modulated by DNA methyltransferases (DNMTs). DNMTs are regulated by transforming growth factor b1 (TGFb1), which in turn controls hepatogenic differentiation and de-differentiation. In addition, a stepwise reduction in TGFb1 concentrations in culture media increases DNMT1 and decreases DNMT3 in primary hepatocytes (Heps) and confers Heps with multi-differentiation potentials similarly to MSCs. Hepatic lineage reversibility of MSCs and lineage conversion of Heps are regulated by DNMTs in response to TGFb1. This previously unrecognized TGFb1-DNMTs-MSC-HD axis may further increase the understanding the normal and pathological processes in the liver, as well as functions of MSCs after transplantation to treat liver diseases.
INTRODUCTION Hepatocytes possess de-differentiation and re-differentiation potential to achieve liver regeneration (Tarlow et al., 2014). Such findings have challenged the long-standing belief that differentiation is an irreversible process in mammals (Jopling et al., 2011; Metcalf, 2007). However, to date detailed mechanisms of how somatic cells undergo de-differentiation and subsequent re-differentiation remain undefined, and whether lineage reversibility exists in adult stem cells remains controversial (Zhau et al., 2011). One useful model for studying lineage reversibility is the multi-potent mesenchymal stromal cells (MSCs) (Makino et al., 1999; Pittenger et al., 1999; Wislet-Gendebien et al., 2005), which can differentiate into hepatocyte-like cells (dHeps) (Jiang et al., 2014; Kuo et al., 2008; Lee et al., 2004) that possess therapeutic potentials for liver dis-
eases. Currently, whether and how differentiated hepatic progenies from MSCs can be de-differentiated is unknown. Studies do show, however, DNA methylation plays an important role in controlling MSC differentiation (Tsai et al., 2012). DNA methyltransferases (DNMTs) control gene transcription and cellular phenotypic changes during liver organogenesis (Snykers et al., 2009). It has been reported that inhibition of DNMTs increases liver-specific gene expression to maintain a hepatic fate; a concomitant decrease of DNMT1 and increase of DNMT3 expression is associated with hepatic maturation in mouse hepatoblasts (Gailhouste et al., 2013). DNMT1 is the key to methylation maintenance, while DNMT3 is responsible for the initiation of de novo methylation (Li, 2002). However, the function roles of each DNMT in hepatic lineage commitment of MSCs are unknown. The purpose of this study is to investigate the hepatic lineage reversibility in MSCs as well as the
Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017 j ª 2017 247 This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
(legend on next page) 248 Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017
underlying molecular mechanisms. We hypothesize that alteration of culture conditions may achieve hepatic lineage reversibility between MSCs and their differentiated hepatic progenies.
RESULTS Mouse MSC-Derived Differentiated Hepatocyte-like Cells Possess De-differentiation Potential Mouse MSCs were induced into dHeps via our previously reported two-step protocol (Jiang et al., 2014; Kuo et al., 2008; Lee et al., 2004). During 4 weeks of hepatogenic differentiation (HD), the bipolar fibroblast-like morphology of the MSCs changed into cobblestone-like shapes from day 14 and then progressively changed into hexagonal shapes, a typical epithelial morphology at day 28 (Figure 1A). Expression of hepatocyte-specific genes was significantly upregulated after differentiation induction in a time-dependent manner (Figures 1B–1D). dHeps possessed hepatocyte surface markers similar to those of mouse Heps (Figure 1E). The in vitro functional characteristics were further investigated by periodic acid-Schiff (PAS) staining (Figure 1F), measurements of glycogen, albumin, and urea production (Figure 1G), as well as low-density lipoprotein (LDL) uptake assays (Figure 1H). All the hepatic function markers significantly and progressively increased during differentiation. In general, although slightly lower, dHeps expressed function markers similar to those of Heps (Figures 1B–1H). Expression levels of hepatic-specific genes of dHeps are close to those of embryonic day 12.5 (E12.5) fetal livers and FL83B cells, and differ relatively from those of Heps (Figure S2). Compared with Heps, dHeps also expressed higher levels of Sca-1 (Figure 1E), Dlk1, and Cd24a (Figure S2), which are hepatic progenitor markers (Miyajima et al., 2014; Petersen et al., 2003; Qiu et al., 2011). These
findings indicated that dHeps are similar to hepatic progenitors. Although dHeps are not identical to terminally differentiated mature Heps, they are committed hepatic cells. After 28 days of HD, we de-differentiated dHeps by replacing differentiation medium with mouse MSC-maintenance medium for 28 days. Following hepatogenic de-differentiation (dHD), cellular morphology reverted back to fibroblast-like shapes (Figure 1A). Furthermore, liver-specific gene expression, surface marker presentation, and hepatic functions all switched back to the original levels observed in original MSCs (Figures 1B–1H). We next examined cell-cycle profiles during HD and dHD. After differentiation for 28 days, the cells arrested at the G0/G1 phase but re-entered the cell cycle after de-differentiation for 28 days (Figure 1I). The reversibility of the cell cycle was accompanied by changes in expression of genes encoding cell-cycle checkpoint regulators as well as proliferation markers Ki67 and Pcna (Figures 1J and 1K). The results indicated that hepatic lineage reversibility in MSCs was accompanied by cell-cycle exit and re-entry. To explore the transcriptomic and proteomic changes during HD and dHD, we performed microarray analysis and high-density antibody arrays on MSCs (HD 0), dHeps (HD 28), and de-differentiated cells from dHeps (ddHeps, dHD 28). Hierarchical clustering (Figures 2A and 2B) showed that ddHeps clustered closely together with MSCs and were separated from dHeps and Heps. The profiles of genes involved in fibroblast markers, MSC markers (Figure 2C), lipid, glycolysis, cholesterol, and drug metabolism (Figure 2D) were very similar between dHeps and Heps, whereas ddHeps and MSCs displayed similar gene profiles. Collectively, these results indicated that the ddHep transcriptome and proteome had indeed reverted back to the MSC state. Using the KEGG database, gene ontology was analyzed to understand the functional significance of differential gene expression during HD and dHD. Pathways involved in
Figure 1. Hepatogenic Differentiation and De-differentiation of Mouse Mesenchymal Stromal Cells (A) Representative morphological changes in MSCs during hepatogenic differentiation (HD) and hepatogenic de-differentiation (dHD). Scale bar, 100 mm. (B–D) Results of (B) qRT-PCR (n = 4 independent experiments), (C) a representative immunoblot, and (D) immunocytochemical staining for hepatogenic-specific gene expression. b-Actin was used as an internal control. Scale bar, 50 mm. (E) Flow-cytometry analysis for surface phenotypic characterization. (F) PAS (periodic acid-Schiff) staining for glycogen storage. (G and H) Hepatic functions were measured via (G) albumin production (n = 7 independent experiments), glycogen storage (n = 4 independent experiments), urea production (n = 3 independent experiments), and (H) LDL-uptake assays. LDLR, LDL receptor. (I) Results of flow-cytometry analysis showing the relative percentage of cells in various cell-cycle phases (n = 5 independent experiments). (J and K) Results of (J) qRT-PCR (n = 4 independent experiments) and (K) a representative immunoblot for cell-cycle-related genes. ‘‘HD 0’’ refers to MSCs; ‘‘HD 28’’ refers to MSC-derived dHeps (hepatocyte-like cells) after 28 days of HD; ‘‘dHD 28’’ refers to dHep-derived ddHeps (de-differentiated cells) after 28 days of dHD; ‘‘Heps’’ refers to mouse primary hepatocytes. All quantitative data are presented as means ± SD. Statistical analyses were performed using Student’s paired t test, with significance set at **p < 0.001. NS, not significant. See also Figures S1–S5 and S7. Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017 249
(legend on next page) 250 Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017
both general (e.g., cell cycle) and specific functions (e.g., lysosome and metabolism) were demonstrated (Table S1). A substantial overlap of KEGG pathways between HD and dHD was also observed. Moreover, microarray results revealed that upregulated pathways during HD substantially overlapped with downregulated pathways during dHD (Figure 2E, left panel), and vice versa (Figure 2E, right panel). These patterns indicated that dHD involved the pathways associated with HD and is essentially a reversal of the latter process. ddHeps Exhibit Potential for Multi-lineage Differentiation To investigate whether ddHeps regained multi-lineage differentiation potential, an essential function of MSCs, we examined their hepatogenic, osteogenic, and adipogenic differentiation capabilities. Hepatic induction caused ddHeps to become functional dHeps, evidenced by morphological changes, the capacity for glycogen storage and albumin production (Figure 3A), as well as the reexpression of hepatogenic-specific genes. Positive function assays, including glycogen storage, albumin production, and generation of urea, and downregulation of Ki67 also supported the manifestation of hepatic re-differentiation (Figures 3B and 3C). Next, osteogenic differentiation of ddHeps was evidenced by morphologic change, alkaline phosphatase activity, and calcium mineralization of the extracellular matrix (Figure 3D), visible through von Kossa and alizarin red S staining. Osteogenic-specific gene expression was also increased significantly (Figure 3E). Finally, ddHeps differentiation into adipocytes was shown by the accumulation of neutral lipid vacuoles (Figure 3F), positive oil red O staining, and increased adipogenic-specific gene expression (Figure 3G). To further confirm dHD and to rule out the possibility of incomplete HD, we sorted dHeps by staining for CD26, which is only expressed on mature hepatocytes but not in MSCs or immature hepatic progenitor cells. After dHD for 4 weeks, ddHeps derived from CD26high dHeps expressed a typical MSC phenotype (Figures S3A–S3D) and regained multi-lineage differentiation potentials (Figures S3E–S3G). Expression of lineage-specific genes in ddHeps upon various lineage induction is similar to that of primary
MSCs (Figures S3E–S3G). ddHeps lacked cell-surface antigens CD45 and CD90 (Figure S3D) and did not express markers of liver progenitors reported in the literature (Chen et al., 2012; Oikawa et al., 2009). Thus, ddHeps are indeed closer to MSCs rather than committed liver progenitor cells, oval cells, or hepatoblasts. Hepatic Lineage Reversibility Is Associated with MET and EMT Mesenchymal-to-epithelial transition (MET) and epithelial-to-mesenchymal transition (EMT) are reversible processes that dictate cell morphology, migration, and differentiation (Thiery and Sleeman, 2006). To examine whether MET and EMT were involved during HD and dHD, we surveyed the expression of mesenchymal and epithelial marker genes. During HD of MSCs, epithelial markers such as Cdh1, Epcam, Gja1, and Tjp1 were upregulated, whereas mesenchymal markers such as Vim, Fn1, Acta1, and Acta2, as well as the master EMT regulators Tgfb1 and Snai1, were downregulated. A reverse pattern was observed during dHD (Figures S4A–S4C). In addition, migratory capacity was decreased after HD and recovered upon dHD (Figure S4D). We concluded that HD is associated with MET, whereas dHD is associated with EMT. Reversible HD Also Occurs in Human MSCs We applied the same approach to induce HD and dHD in human MSCs. Reversible changes of morphology, hepatogenic-specific gene profile, hepatic functions, MET markers, and migration capacity (Figures S5A–S5F) were also observed in human MSCs. Moreover, human ddHeps retained multi-lineage differentiation potential (Figures S5G–S5I). Overall, the results from mouse and human MSCs were consistent; the fate of human MSCs could also be reversed after differentiation into the hepatic lineage. DNA Methyltransferases Regulate HD of MSCs and dHD To investigate the role of epigenetic effects during HD and dHD of MSCs, we used real-time qRT-PCR to screen the expression of 19 epigenetic modulators in DNA and histone methylation pathways (Figure 4A). We found that
Figure 2. Microarrays and High-Density Antibody Arrays of Transcriptomes and Proteomes Reveals Hepatogenic Differentiation and De-differentiation (A–D) Hierarchical clustering via correlations in the transcriptome (A) and the proteome (B). Expression patterns of genes involved in fibroblast, MSC, and EMT markers (C). Expression patterns of genes involved in lipid, glucose, cholesterol, and drug metabolism (D). Each column represents a single array sample. (E) Venn diagrams showing the overlap of up- and downregulated KEGG pathways within the transcriptome during HD and dHD. ‘‘HD 0’’ refers to MSCs; ‘‘HD 28’’ refers to MSC-derived dHeps after 28 days of HD; ‘‘dHD 28’’ refers to dHep-derived ddHeps after 28 days of dHD; ‘‘Heps’’ refers to mouse primary hepatocytes; ‘‘MLCs’’ refers to MSC-like cells derived from Heps. See also Tables S1 and S2. Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017 251
Figure 3. dHeps Can Also De-differentiate and Regain Multi-lineage Differentiation Potentials (A) Representative morphological changes and staining (PAS and albumin [ALB]) of ddHeps during hepatogenic re-differentiation (reHD). (B) qRT-PCR of hepatogenic-specific genes inddHep-derived dHeps (n = 3 independent experiments). (C) Hepatic functions of ddHep-derived dHeps were evaluated with albumin production (n = 3 independent experiments), glycogen storage (n = 6 independent experiments), and urea production (n = 3 independent experiments). (legend continued on next page) 252 Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017
DNMT1 expression was decreased during HD but was strongly increased after dHD (Figures 4A and 4B), to levels higher than the nascent MSCs (HD 0). In contrast, DNMT3a and DNMT3b expression was increased during HD and returned to the original levels similar to MSCs after dHD. In addition, pretreatment with 5-azacytidine (a panDNMT inhibitor) reduced global DNA methylation in MSCs (Figure 4C). Over time, increasing doses of 5-azacytidine enhanced hepatogenic-specific gene expression and hepatic functions (Figures 4D and 4E). These results indicated that DNMTs are involved in the HD of MSCs. We then performed loss-of-function assays to determine the roles of individual DNMTs in HD. Small interfering RNAs (siRNAs) against DNMT1, DNMT3a, and DNMT3b were added individually during HD (Figures 4F and 4G). The results revealed that in HD 14, DNMT1 knockdown significantly enhanced the expression of hepatogenicspecific genes, MET/EMT markers (Figure 4H), and the metabolically important cytochrome P450 (CYP) family members (Figure 4I). Furthermore, hepatic functions were enhanced (Figure 4J). An acceleration of HD via DNMT1 knockdown was observed: HD 14 cells with DNMT1 knockdown are similar to HD 28 cells in terms of gene expression and hepatic functions (Figures 4H–4J). Knockdown of DNMT3a only downregulated several hepatogenic markers (Figure 4H), but significantly reduced the expression of all CYP enzymes examined (Figure 4I). Knockdown of DNMT3b had no influence on HD (Figures 4H–4J). Simultaneous silencing of DNMT3a and DNMT3b did not further repress HD at day 14 (Figures 4H–4J). In terms of morphological changes during dHD, DNMT1 knockdown maintained epithelial morphology and increased albumin expression compared with parental and control groups. In contrast, DNMT3a knockdown promoted a fibroblast-like morphology and decreased albumin expression (Figures 4K and 4L). Collectively, these data showed that DNMT1 extensively regulated hepatic lineage commitment but DNMT3a specifically regulated the CYP family during the HD of MSCs. Moreover, DNMT1 silencing accelerated HD.
15 candidate genes, four (Tgfb1, Thbs1, Thbs2, and Thbs3) were involved in the TGFb signaling pathways. Among these four, Tgfb1 was highly expressed in both arrays (Figure 5A). Increased TGFb1 expression was strongly associated with reversible differentiation (Figures 5B and S5E). To test how DNMTs responded to TGFb1, we altered the concentrations of TGFb1 during HD. As TGFb1 dosage increased, DNMT1 expression also increased in dHeps, whereas DNMT3a and DNMT3b expression decreased (Figure 5C). These effects were eliminated by LY-364947, an Alk-5 (TGFbRI/II) inhibitor, indicating that extracellular TGFb1 signaling is Alk-5 dependent. To identify the signaling pathways involved in TGFb1-regulated DNMT expression during HD, we tested both Smad and nonSmad signaling pathway inhibitors (Figure 5D). We showed that PD0325901 and SIS3 treatment reversed DNMT1 upregulation, SB203580 and LY294002 treatment reversed DNMT3a/3b downregulation, and SP600125 treatment did both. Therefore, TGFb1 may regulate DNMT1 and DNMT3a/3b together via the JNK pathway, whereas it separately regulates DNMT1 via the ERK1/2 and canonical Smad2/3 pathways, as well as DNMT3a/3b via the p38 and PI3K-Akt pathways. Of note, neither DNMT knockdowns (Figures 5E and 5F) nor 5-azacytidine pretreatment (Figure 5G) affected TGFb1 expression. These results demonstrated that DNMTs regulate HD in response to TGFb1 signaling. Moreover, a feedback loop does not appear to exist in the TGFb1-DNMT axis during the HD of MSCs. To determine whether the downregulation of TGFb1 signaling is essential for HD, we treated MSCs with TGFb1 during HD. Our results showed that TGFb1 repressed morphologic changes (Figure 5H), hepatogenicspecific gene expression (Figure 5I), and hepatic functions (Figure 5J), whereas it increased cell migratory capacity (Figure 5K). All effects from TGFb1 treatment were reduced or eliminated by LY-364947. Our results confirmed that TGFb1 signaling must be downregulated for the HD of MSCs.
Identification of Signaling Pathways Involved in DNMT-Modulated HD and dHD To identify the signaling pathways involved in DNMTdependent HD, we integrated the microarray and protein array data associated with DNMT1 expression. Among
Murine Hepatocytes Can Be Converted to Functional MSC-like Cells Since dHeps are similar to hepatic progenitors and are different from Heps, we further investigate whether similar phenomenon of de-differentiation from dHeps to MSCs can
(D) Representative morphological changes and staining (alkaline phosphatase [ALP], von Kossa, and alizarin red S) of ddHeps during osteogenic re-differentiation (reOD). (E) qRT-PCR of osteogenic-specific genes in ddHep-derived osteoblasts (n = 3 independent experiments). (F) Representative morphological changes and staining (oil red O) of ddHeps during adipogenic re-differentiation (reAD). (G) qRT-PCR of adipogenic-specific genes in ddHep-derived adipocytes (n = 3 independent experiments). All quantitative data are presented as means ± SD. Statistical analyses were performed using Student’s paired t test, with significance set at *p < 0.05, **p < 0.005. NS, not significant. Scale bars represent 100 mm for phase contrast and staining. See also Figures S5 and S7. Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017 253
(legend on next page) 254 Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017
also occur in Heps. To investigate whether TGFb1 can induce lineage conversion in hepatocytes, we treated primary hepatocytes with TGFb1 or TGFb1 receptor inhibitor (LY-364947) for 3–7 days. Treatment of Heps with TGFb1 for 7 days resulted in a fibroblast-like morphology and a loss of hepatogenic-specific genes. These effects were eliminated by LY364947 administration (Figure S7). Expression of hepatogenic-specific genes was increased in LY-364947-treated Heps compared with medium-only control on day 7 (Figure S7). Blocking of the endogenous TGFb1 signaling repressed the loss of hepatocyte characteristics, supporting that TGFb1 is necessary and sufficient to trigger hepaticmesenchymal lineage conversion. Currently known actions of TGFb1 include inducing apoptosis, de-differentiating fetal hepatocytes to liver progenitors, and mediating the acquisition of mesenchymal morphology (Caja et al., 2011; Tu et al., 2014). Additionally TGFb1 is involved in the myofibroblastic (Kakudo et al., 2012), chondrogenic (Toh et al., 2005), and cardiomyogenic differentiation (Mohanty et al., 2013) of MSCs. Therefore, we attempted to convert purified CD26high Heps to MSC-like cells (MLCs) by treating them with stepwise reduction of TGFb1 (Figure 6A). This stepwise protocol converted CD26high Heps to an MSC-like status in terms of morphology, hepatogenic-specific gene expression, cell-surface phenotype, hepatic functions (Figures 6B–6G), and EMT markers (Figures 6H and 6I). MLCs re-entered into cell cycle and increased the cell number in culture, and re-thawed MLCs retained proliferative capacity (Figures 6J and 6K). The transcriptomes and proteomes of MLCs and MSCs were similar (Figures 2A and 2B). The two cell types also exhibited similarities in the expression of metabolismrelated genes (Figure 2D) and other gene sets (Figure 2C). Gene ontology analyses revealed substantial overlap of enriched pathways during Hep lineage conversion to MLCs and during reversible HD of MSCs (Figure 2E and
Table S2). After treatment with TGFb1 we noted the induction of endogenous TGFb1 expression in MLCs, which was maintained after exogenous TGFb1 removal (Figures 6H and 6L). In addition, DNMT1 knockdown prevented TGFb1-induced loss of hepatogenic-specific genes in Heps (Figure 6M), and DNMT expression patterns were consistent between Hep-derived MLCs and MSCs (Figures 6N and 6O), providing evidence that DNMTs mediated the TGFb1-induced lineage conversion of Heps. Finally, MLCs possessed multi-lineage potentials comparable with those of MSCs (Figures 7A–7G); while they exhibited higher hepatic potential, particularly in the expression of CYP family enzymes (Figure 7B). Together, these results indicated that Hep-converted MLCs are very similar to primary MSCs. Reversible Differentiation Does Not Induce Tumorigenicity Long-term culture expansion, changes to methylation states, and DNMT deficiency may cause genomic instability and lead to malignant transformation (Foudah et al., 2009; Josse et al., 2010; Karpf and Matsui, 2005; Miura et al., 2006; Rosland et al., 2009). Moreover, DNMT1 expression has been associated with hepatocellular carcinomas (Raggi et al., 2014). We found a deletion of Y chromosome and a duplication of chromosome 16 after the long-term differentiation and de-differentiation process revealed by comparative genomic hybridization arrays (n = 1, data not shown). To test whether differentiation and de-differentiation process induces tumorigenicity, we subcutaneously inoculated nude mice with mouse MSCs, ddHeps, dHeps with siRNA against DNMT1, and MLCs. We observed no tumor formation after 13 weeks, indicating that tumorigenicity is not induced by long-term reversible differentiation, genetic manipulation of DNMTs, or lineage conversion (Figures 7H and 7I).
Figure 4. DNA Methyltransferases Regulate Hepatogenic Differentiation and De-differentiation of Mouse MSCs (A) A heatmap summarizing relative mRNA levels of 19 epigenetic modulators (n = 4). (B) Representative immunoblot for DNA methyltransferases (DNMTs) expression in MSCs during HD and dHD. (C–I) DNMT-inhibitor (5-azacytidine [5-aza]) treatment inhibited global DNA methylation (n = 3 independent experiments). After 14 days of HD, enhancement from pan-DNMT inhibitor pretreatment occurred in hepatogenic-specific gene expression, as visualized with qRT-PCR (n = 4 independent experiments) (D), and in hepatic functions (E), as represented by albumin production (n = 2 independent experiments, four technical replicates), glycogen storage (n = 2 independent experiments, 12 technical replicates), and urea production (n = 2 independent experiments, five technical replicates). Symbols indicate treatments with different 5-aza concentrations (+, 0.5 mM; ++, 1 mM; +++, 2 mM). Efficiency of DNMT knockdown determined with qRT-PCR (F) (n = 6 independent experiments) and immunoblot (G) at HD 14. Knockdown of DNMTs altered hepatogenic-specific genes, EMT/MET markers (H), and CYP enzymes (I) at HD 14 (n = 6 independent experiments). (J–L) Knockdown of DNMTs altered hepatic functions at HD 14, as represented by albumin production (n = 6 independent experiments), glycogen storage (n = 6 independent experiments), and urea production (n = 6 independent experiments). Knockdown of DNMTs altered cellular morphology at dHD 3 (K) and modulated albumin expression, as seen in the representative images of immunocytochemical staining (L). All quantitative data are presented as means ± SD. Statistical analyses were performed using Student’s paired t test, with significance set at *p < 0.05, **p < 0.005. NS, not significant. Scale bars represent 100 mm for phase contrast and staining. See also Figure S7. Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017 255
Figure 5. TGFb1 Regulates DNMTs and Represses Hepatogenic Differentiation of Mouse MSCs (A) Schematic illustration of signaling pathways. PCC, Pearson’s correlation coefficient. (B) ELISA results for TGFb1 expression of MSCs during HD and dHD (n = 5 independent experiments). (C) qRT-PCR of DNMTs from MSCs treated with TGFb1 during HD (n = 6 independent experiments). Symbols indicate treatments of different TGFb1 concentrations (+, 1 ng/mL; ++, 10 ng/mL; +++, 100 ng/mL). (D–F) qRT-PCR of DNMTs from MSCs treated with TGFb1 and various inhibitors during HD (n = 4 independent experiments) (D). LY-364947, Alk-5 inhibitor, 400 nM; SIS3, Smad2/3 inhibitor, 10 mM; SP600125, JNK inhibitor, 30 mM; SB203580, p38 inhibitor, 30 mM; PD0325901, MEK-ERK1/2 inhibitor, 1 mM; and LY294002, PI3K-Akt inhibitor, 25 mM. The mRNA expression from each treatment was statistically (legend continued on next page) 256 Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017
DISCUSSION In this study, we show that hepatic lineages are reversible between adult MSCs and hepatocytes. Such reversibility is regulated by the TGFb1-DNMT axis. To the best of our knowledge, this is also the demonstration that hepatocytes can be converted into multi-potent MLCs without ectopic gene delivery. Our findings support the occurrence of lineage reversibility in mammalian cells, indicating that such a phenomenon is not restricted to plants and non-mammalian vertebrates (Brockes, 1997; Brockes and Kumar, 2002; Weigel and Jurgens, 2002). Our findings demonstrate a role of DNMTs in hepatic lineage reversibility of MSCs, with individual DNMTs exerting opposite effects; DNMT1 repressed HD and promoted dHD; while DNMT3a promoted HD and repressed dHD. The two enzymes seem to have both unique and overlapping target cytosine-phosphate-guanine (CpG) sites, as shown by the promotion of all hepatogenic-specific genes and hepatic functions via DNMT1 silencing, but a more specific repression of the CYP family via DNMT3a silencing (Figure 4). In addition, the inhibitory effect of DNMT1 appears to be more predominant than the activating effect of DNMT3a, because a pan-DNMT inhibitor resulted in the promotion of HD (Figure 4). Although DNMT3b expression is also strongly associated with HD and responded to TGFb1 stimulation, DNMT3b knockdown fails to affect hepatic gene expression and function. While more data are necessary to clarify the function of DNMT3b in HD, we speculate that the enzyme may cooperate with other chromatin modifiers. Epigenetic effects exerted by these modifiers may have been sufficient for HD regulation, explaining why we do not observe significant changes resulting from DNMT3b knockdown (Baubec et al., 2015; Morselli et al., 2015). In addition to demonstrating the importance of DNMTs, we show that TGFb1 differentially regulated DNMTs, and that TGFb1 is both necessary and sufficient to induce lineage conversion of hepatocytes into MLCs with phenotypes and multipotency resembling MSCs. Upon withdrawal of TGFb1 on day 28, we observed that Hep-converted MLCs
possessed MSC-like DNMTs expression patterns (Figures 6N and 6O) due to the induction of endogenous TGFb1 expression (Figures 6H and 6L). This outcome supports the existence of an MSC-specific methylome that maintains MLC phenotypes and multipotency. a-Fetoprotein (AFP) is a marker indicating early hepatic fate. Although AFP is not detectable in hepatocytes in vivo, expression of AFP may be detected in the in vitro culture conditions owing to a loss of contact (Gleiberman et al., 1989). Moreover, TGFb1 downregulated albumin and AFP expression and induced fibroblast morphology and mesenchymal marker expression (mesenchymal fate) in Hep3B cells (Caja et al., 2011). Therefore, AFP may serve as a marker to distinguish hepatic fate from mesenchymal fate and allows us to track cellular conversion from the former to the latter. Fibrosis is often a consequence of various chronic liver diseases and leads to hepatic failure as well as cancer. TGFb1, a master inducer of EMT, orchestrates fibrogenesis. Persistent TGFb1 stimulation converts both hepatocytes and MSCs into myofibroblasts and participates in liver fibrosis, and even contributes to the development of tumor (Sekiya et al., 2016). While blocking TGFb signaling in hepatocytes alone attenuates liver fibrosis (Dooley et al., 2008) and facilitates liver regeneration (Karkampouna et al., 2016), MLCs may be at an intermediate state during the transition to myofibroblasts. Therefore, inhibition of the TGFb1-DNMT axis in vivo may convert Hepderived myofibroblasts back into MLCs or hepatocytes, and may attenuate liver fibrosis and facilitate liver regeneration. Our previous studies established a protocol for functional dHeps with therapeutic potentials when transplanted in vivo (Jiang et al., 2014; Kuo et al., 2008; Lee et al., 2004). However, dHeps required 28 days of in vitro differentiation, a lengthy duration that hampers their clinical application. In this study, we successfully accelerate differentiation to 14 days by silencing DNMT1 (Figure 4). Transplantation of dHeps with TGFb receptor inhibitor or DNMT1 antagonist may stabilize hepatic characteristics of transplanted dHeps. On the other hand, although transplanted dHeps may de-differentiate
compared with a TGFb1-treated group. Control refers to MSC-derived dHeps after 14 days. Results from qRT-PCR (E) (n = 6 independent experiments) and ELISA (F) (n = 3 independent experiments) of TGFb1 expression at HD 14 after DNMT knockdown. (G) DNMT-inhibitor pretreatment did not affect TGFb1 expression (n = 3 independent experiments). Symbols indicate treatments with different 5-azacytidine concentrations (+, 0.5 mM; ++, 1 mM; +++, 2 mM). (H) Representative morphological changes of MSCs treated with TGFb1 during HD. (I) qRT-PCR of hepatogenic-specific genes in dHeps under TGFb1 treatment (n = 5 independent experiments). (J) Hepatic functions of dHeps under TGFb1 treatment, as represented by albumin production (n = 7 independent experiments), glycogen storage (n = 5 independent experiments), and urea production (n = 6 independent experiments). (K) Quantification of transwell migration assay (n = 3 independent experiments) of MSCs treated with TGFb1 during HD. All quantitative data are presented as means ± SD. Statistical analyses were performed using Student’s paired t test, with significance set at *p < 0.05, **p < 0.005. NS, not significant. Scale bars represent 100 mm for phase contrast and staining. See also Figures S6 and S7. Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017 257
(legend on next page) 258 Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017
into ddHeps, ddHeps may promote liver regeneration like MSCs through immunomodulation and trophic support (Christ et al., 2015). Hep-converted MLCs have proliferation ability and higher hepatic lineage potential (Figure 7B). Therefore, these cells may also serve as an alternative cell source for liver regeneration (Du et al., 2014; Takeuchi et al., 2015). In summary, hepatic lineage reversibility exists and is regulated by the previously unrecognized TGFb1-DNMT axis. The ability to convert Heps into MLCs with altered extracellular concentrations of TGFb1 suggests that terminally differentiated hepatocytes may exhibit considerably greater plasticity than previously believed. Our findings have not only offered therapeutic insights for use of MSCs or differentiated hepatic progenies for transplantation to treat liver diseases, but have also contributed to the development of approaches for achieving tissue regeneration through promoting de-differentiation and lineage conversion in vivo.
sion of surface markers and ability of multi-potent differentiation into osteoblasts, adipocytes, and chondroblasts in vitro (Figure S1). Mouse MSCs (mMSCs) were cultured in an mMSC-maintenance medium, consisting of low-glucose DMEM supplemented with 10% fetal bovine serum (FBS) (Gibco, Thermo Fisher Scientific); human MSCs (hMSCs) were cultured in hMSC-maintenance medium, consisting of Iscove’s modified Dulbecco’s medium (IMDM) supplemented with 10% FBS, 10 ng/mL epidermal growth factor, and 10 ng/mL basic fibroblast growth factor (bFGF) (R&D Systems). Human primary hepatocytes (hHeps) were purchased from ScienCell and cultured in Hepatocyte Medium (ScienCell). Mouse hepatocytes were isolated from 10- to 12-week-old BALB/c mice with a protocol in the literature (Klaunig et al., 1981a, 1981b), and were cultured in Hep-maintenance medium consisting of DMEM/F12 supplemented with 10% FBS. FL83B cells were cultured in F-12K medium (Gibco, Thermo Fisher Scientific) supplemented with 10% FBS. Low-glucose DMEM, IMDM, DMEM/ F12, 5-azacytidine, TGFb1, and LY-364947 were purchased from Sigma-Aldrich; SP600125, SB203580, PD325901, LY294002, and SIS3 were purchased from Abcam.
Hepatogenic Differentiation and De-differentiation
EXPERIMENTAL PROCEDURES All animal studies were approved by the Taipei Veterans General Hospital Institutional Animal Care and Use Committee (IACUC 2014-043).
Cell Culture and Reagents MSCs were obtained from bone marrow and were maintained in culture as previously described (Jiang et al., 2014; Lee et al., 2004). The characteristics of MSCs were confirmed by their expres-
HD was induced following our previously reported two-step protocol (Lee et al., 2004). In brief, 1.0 3 104/cm2 MSCs were seeded and cultured for 2 days in the maintenance medium. Two days post seeding, cells were treated with a step-1 differentiation medium, consisting of IMDM supplemented with 20 ng/mL hepatocyte growth factor and 10 ng/mL bFGF (R&D Systems), and 0.61 g/L nicotinamide (Sigma-Aldrich). Seven days later, cells were treated with a step-2 differentiation medium, consisting of IMDM supplemented with 20 ng/mL oncostatin M (Invitrogen), 1 mmol/L dexamethasone (Sigma-Aldrich), and 50 mg/mL ITS+ premix (BD Biosciences) for 3 weeks. The differentiation medium was changed
Figure 6. Stepwise TGFb1 Treatment Converts Mouse Primary Hepatocytes into MSC-like Cells though Initiating the Secretion of Endogenous TGFb1 (A) Schematic illustration of the experimental approach. (B) Representative flow-cytometry plots showing CD26 expression and the gating strategies used to isolate CD26low and CD26high hepatocytes. (C) Representative morphology and PAS staining for MLCs. (D and E) Loss of hepatogenic-specific gene expression in MLCs, analyzed with qRT-PCR (D) (n = 4 independent experiments) and immunoblot (E). (F) Loss of hepatic function in MLCs, as represented by albumin production (n = 4 independent experiments), glycogen storage (n = 3 independent experiments), and urea production (n = 4 independent experiments). (G) Surface phenotypic characterization analyzed with flow cytometry. (H and I) Gene expression of MET/EMT markers during de-differentiation of mouse Heps, analyzed with qRT-PCR (H) (n = 4 independent experiments) and immunoblot (I). (J and K) Cell numbers counted at each time point (n = 3 independent experiments) (J) and relative percentage of cells in various cell-cycle phases (n = 5 independent experiments) (K). ‘‘Re-thawed MLCs’’ refers to MLCs re-thawed from cryopreservation and cultured for another week at least. (L) ELISA of TGFb1 expression in Heps, MLCs, and MSCs (n = 4 independent experiments); the amount of TGFb1 from serum has been deducted. (M) qRT-PCR of hepatogenic-specific genes in mouse Heps under TGFb1 and siDnmt1 treatment for 3 days (n = 3 independent experiments). The plus and minus symbols respectively indicate treatments with and without TGFb1 (10 ng/mL), siControl (50 nM), and siDnmt1 (50 nM). (N and O) Expression of DNMTs in hepatocytes, MLCs, and MSCs, analyzed by qRT-PCR (N) (n = 4 independent experiments) and immunoblot (O). *p < 0.05, **p < 0.005. NS, not significant. Scale bars represent 100 mm for phase contrast and staining. See also Figures S6 and S7. Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017 259
(legend on next page) 260 Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017
once per week. To induce dHD in MSC-derived dHeps, we replaced the step-2 medium with maintenance medium (low-glucose DMEM supplemented with 10% FBS). A protocol with stepwise reduction of TGFb1 was used to induce the lineage conversion of Heps. Specifically, Heps were cultured in an mMSC-maintenance medium with 10 ng/mL TGFb1 for the first 7 days, then with 5 ng/mL TGFb1 for another 7 days, and finally with 1 ng/mL TGFb1 for the following 14 days. After 28 days, cells were cultured in maintenance medium without TGFb1. Other experimental procedures and statistical analysis are described in Supplemental Experimental Procedures.
from Ministry of Education. The authors also acknowledge the technical support provided by Flow Cytometry Core Facility of National Yang Ming University and Imaging Core Facility of Nanotechnology of the UST-YMU. We acknowledge Prof. Ting-Fen Tsai from National Yang-Ming University, Taiwan for the kind gift of E12.5 fetal livers and thank Prof. Chi-Chang Juan from National Yang-Ming University, Taiwan for the kind gift of FL83B cells. Received: June 13, 2016 Revised: May 4, 2017 Accepted: May 5, 2017 Published: June 8, 2017
ACCESSION NUMBERS With regard to transcript profiling the accession numbers are: microarray, GEO: GSE73617; protein array, GEO: GSE73762.
SUPPLEMENTAL INFORMATION Supplemental Information includes Supplemental Experimental Procedures, seven figures, and three tables and can be found with this article online at http://dx.doi.org/10.1016/j.stemcr.2017. 05.008.
AUTHOR CONTRIBUTIONS C.-W.L. designed the study, performed experiments, analyzed data, and wrote the manuscript; W.-C.H. and H.-D.H. performed microarray and protein array analysis; Y.-H.H., J.H.H., and M.-H.Y. provided feedback for the manuscript; V.W.Y. and O.K.L. edited the manuscript and supervised the study.
ACKNOWLEDGMENTS This work was supported in part by the Novel Bioengineering and Technological Approaches to Solve Two Major Health Problems in Taiwan sponsored by the Taiwan Ministry of Science and Technology Academic Excellence Program under grant number MOST 105-2633-B-009-003. The authors acknowledge financial support from the Ministry of Science and Technology, Taiwan (MOST1032314-B-010-053-MY3, MOST 106-2321-B-010-008, MOST 1062911-I-010-502, and MOST 105-2633-B-010-002). This study was also supported by Aiming for the Top University Plan, a grant
REFERENCES Baubec, T., Colombo, D.F., Wirbelauer, C., Schmidt, J., Burger, L., Krebs, A.R., Akalin, A., and Schubeler, D. (2015). Genomic profiling of DNA methyltransferases reveals a role for DNMT3B in genic methylation. Nature 520, 243–247. Brockes, J.P. (1997). Amphibian limb regeneration: rebuilding a complex structure. Science 276, 81–87. Brockes, J.P., and Kumar, A. (2002). Plasticity and reprogramming of differentiated cells in amphibian regeneration. Nat. Rev. Mol. Cell Biol. 3, 566–574. Caja, L., Bertran, E., Campbell, J., Fausto, N., and Fabregat, I. (2011). The transforming growth factor-beta (TGF-beta) mediates acquisition of a mesenchymal stem cell-like phenotype in human liver cells. J. Cell. Physiol. 226, 1214–1223. Chen, Y., Wong, P.P., Sjeklocha, L., Steer, C.J., and Sahin, M.B. (2012). Mature hepatocytes exhibit unexpected plasticity by direct dedifferentiation into liver progenitor cells in culture. Hepatology 55, 563–574. Christ, B., Bruckner, S., and Winkler, S. (2015). The therapeutic promise of mesenchymal stem cells for liver restoration. Trends Mol. Med. 21, 673–686. Dooley, S., Hamzavi, J., Ciuclan, L., Godoy, P., Ilkavets, I., Ehnert, S., Ueberham, E., Gebhardt, R., Kanzler, S., Geier, A., et al. (2008). Hepatocyte-specific Smad7 expression attenuates TGF-beta-mediated fibrogenesis and protects against liver damage. Gastroenterology 135, 642–659.
Figure 7. Mouse Hepatocyte-Derived MSC-like Cells Retain Multi-lineage Differentiation Potential In Vitro (A) Representative morphology and PAS staining for MLC-derived dHeps. (B) qRT-PCR analysis of hepatogenic-specific gene expression (n = 3 independent experiments). (C) Cellular hepatic functions (n = 4 independent experiments). (D) Representative morphology and staining (alkaline phosphatase [ALP], von Kossa, and alizarin red S [n = 3 independent experiments]) of MLC-derived osteocytes. (E) qRT-PCR analysis of osteogenic-specific gene expression (n = 3 independent experiments). (F) Representative morphology and staining (oil red O) of MLC-derived adipocytes. (G) qRT-PCR analysis of adipogenic-specific gene expression (n = 3 independent experiments). (H) Representative photograph comparing tumorigenicity from mouse MSCs, ddHeps, dHeps with siDnmt1, and Hep3B cells in BALB/c nude mice after inoculation for 13 weeks (n = 6 mice per group). Black arrows indicate injection sites. (I) Representative photograph comparing tumorigenicity of MLCs and Hep3B after inoculation for 8 weeks. Black arrows indicate injection sites. All quantitative data are presented as means ± SD. Statistical analyses were performed using Student’s paired t test, with significance set at *p < 0.05, **p < 0.005. NS, not significant. Scale bars represent 100 mm for phase contrast and staining. Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017 261
Du, Y., Wang, J., Jia, J., Song, N., Xiang, C., Xu, J., Hou, Z., Su, X., Liu, B., Jiang, T., et al. (2014). Human hepatocytes with drug metabolic function induced from fibroblasts by lineage reprogramming. Cell Stem Cell 14, 394–403.
differentiation of human mesenchymal stem cells. Hepatology 40, 1275–1284.
Foudah, D., Redaelli, S., Donzelli, E., Bentivegna, A., Miloso, M., Dalpra, L., and Tredici, G. (2009). Monitoring the genomic stability of in vitro cultured rat bone-marrow-derived mesenchymal stem cells. Chromosome Res. 17, 1025–1039.
Makino, S., Fukuda, K., Miyoshi, S., Konishi, F., Kodama, H., Pan, J., Sano, M., Takahashi, T., Hori, S., Abe, H., et al. (1999). Cardiomyocytes can be generated from marrow stromal cells in vitro. J. Clin. Invest. 103, 697–705.
Gailhouste, L., Gomez-Santos, L., Hagiwara, K., Hatada, I., Kitagawa, N., Kawaharada, K., Thirion, M., Kosaka, N., Takahashi, R.U., Shibata, T., et al. (2013). miR-148a plays a pivotal role in the liver by promoting the hepatospecific phenotype and suppressing the invasiveness of transformed cells. Hepatology 58, 1153–1165.
Metcalf, D. (2007). Concise review: hematopoietic stem cells and tissue stem cells: current concepts and unanswered questions. Stem Cells 25, 2390–2395.
Gleiberman, A.S., Sharovskaya, Y., and Chailakhjan, L.M. (1989). ‘‘Contact inhibition’’ of alpha-fetoprotein synthesis and junctional communication in adult mouse hepatocyte culture. Exp. Cell Res. 184, 228–234. Jiang, W.C., Cheng, Y.H., Yen, M.H., Chang, Y., Yang, V.W., and Lee, O.K. (2014). Cryo-chemical decellularization of the whole liver for mesenchymal stem cells-based functional hepatic tissue engineering. Biomaterials 35, 3607–3617. Jopling, C., Boue, S., and Izpisua Belmonte, J.C. (2011). Dedifferentiation, transdifferentiation and reprogramming: three routes to regeneration. Nat. Rev. Mol. Cell Biol. 12, 79–89. Josse, C., Schoemans, R., Niessen, N.A., Delgaudine, M., Hellin, A.C., Herens, C., Delvenne, P., and Bours, V. (2010). Systematic chromosomal aberrations found in murine bone marrow-derived mesenchymal stem cells. Stem Cells Dev. 19, 1167–1173. Kakudo, N., Kushida, S., Suzuki, K., Ogura, T., Notodihardjo, P.V., Hara, T., and Kusumoto, K. (2012). Effects of transforming growth factor-beta1 on cell motility, collagen gel contraction, myofibroblastic differentiation, and extracellular matrix expression of human adipose-derived stem cell. Hum. Cell 25, 87–95. Karkampouna, S., Goumans, M.J., Ten Dijke, P., Dooley, S., and Kruithof-de Julio, M. (2016). Inhibition of TGFbeta type I receptor activity facilitates liver regeneration upon acute CCl4 intoxication in mice. Arch. Toxicol. 90, 347–357. Karpf, A.R., and Matsui, S. (2005). Genetic disruption of cytosine DNA methyltransferase enzymes induces chromosomal instability in human cancer cells. Cancer Res. 65, 8635–8639. Klaunig, J.E., Goldblatt, P.J., Hinton, D.E., Lipsky, M.M., Chacko, J., and Trump, B.F. (1981a). Mouse liver cell culture. I. Hepatocyte isolation. In Vitro 17, 913–925. Klaunig, J.E., Goldblatt, P.J., Hinton, D.E., Lipsky, M.M., and Trump, B.F. (1981b). Mouse liver cell culture. II. Primary culture. In Vitro 17, 926–934. Kuo, T.K., Hung, S.P., Chuang, C.H., Chen, C.T., Shih, Y.R., Fang, S.C., Yang, V.W., and Lee, O.K. (2008). Stem cell therapy for liver disease: parameters governing the success of using bone marrow mesenchymal stem cells. Gastroenterology 134, 2111– 2121, 2121.e1–3. Lee, K.D., Kuo, T.K., Whang-Peng, J., Chung, Y.F., Lin, C.T., Chou, S.H., Chen, J.R., Chen, Y.P., and Lee, O.K. (2004). In vitro hepatic
262 Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017
Li, E. (2002). Chromatin modification and epigenetic reprogramming in mammalian development. Nat. Rev. Genet. 3, 662–673.
Miura, M., Miura, Y., Padilla-Nash, H.M., Molinolo, A.A., Fu, B., Patel, V., Seo, B.M., Sonoyama, W., Zheng, J.J., Baker, C.C., et al. (2006). Accumulated chromosomal instability in murine bone marrow mesenchymal stem cells leads to malignant transformation. Stem Cells 24, 1095–1103. Miyajima, A., Tanaka, M., and Itoh, T. (2014). Stem/progenitor cells in liver development, homeostasis, regeneration, and reprogramming. Cell Stem Cell 14, 561–574. Mohanty, S., Bose, S., Jain, K.G., Bhargava, B., and Airan, B. (2013). TGFbeta1 contributes to cardiomyogenic-like differentiation of human bone marrow mesenchymal stem cells. Int. J. Cardiol. 163, 93–99. Morselli, M., Pastor, W.A., Montanini, B., Nee, K., Ferrari, R., Fu, K., Bonora, G., Rubbi, L., Clark, A.T., Ottonello, S., et al. (2015). In vivo targeting of de novo DNA methylation by histone modifications in yeast and mouse. Elife 4, e06205. Oikawa, T., Kamiya, A., Kakinuma, S., Zeniya, M., Nishinakamura, R., Tajiri, H., and Nakauchi, H. (2009). Sall4 regulates cell fate decision in fetal hepatic stem/progenitor cells. Gastroenterology 136, 1000–1011. Petersen, B.E., Grossbard, B., Hatch, H., Pi, L., Deng, J., and Scott, E.W. (2003). Mouse A6-positive hepatic oval cells also express several hematopoietic stem cell markers. Hepatology 37, 632–640. Pittenger, M.F., Mackay, A.M., Beck, S.C., Jaiswal, R.K., Douglas, R., Mosca, J.D., Moorman, M.A., Simonetti, D.W., Craig, S., and Marshak, D.R. (1999). Multilineage potential of adult human mesenchymal stem cells. Science 284, 143–147. Qiu, Q., Hernandez, J.C., Dean, A.M., Rao, P.H., and Darlington, G.J. (2011). CD24-positive cells from normal adult mouse liver are hepatocyte progenitor cells. Stem Cells Dev. 20, 2177–2188. Raggi, C., Factor, V.M., Seo, D., Holczbauer, A., Gillen, M.C., Marquardt, J.U., Andersen, J.B., Durkin, M., and Thorgeirsson, S.S. (2014). Epigenetic reprogramming modulates malignant properties of human liver cancer. Hepatology 59, 2251–2262. Rosland, G.V., Svendsen, A., Torsvik, A., Sobala, E., McCormack, E., Immervoll, H., Mysliwietz, J., Tonn, J.C., Goldbrunner, R., Lonning, P.E., et al. (2009). Long-term cultures of bone marrowderived human mesenchymal stem cells frequently undergo spontaneous malignant transformation. Cancer Res. 69, 5331–5339. Sekiya, S., Miura, S., Matsuda-Ito, K., and Suzuki, A. (2016). Myofibroblasts derived from hepatic progenitor cells create the tumor microenvironment. Stem Cell Reports 7, 1130–1139. Snykers, S., Henkens, T., De Rop, E., Vinken, M., Fraczek, J., De Kock, J., De Prins, E., Geerts, A., Rogiers, V., and Vanhaecke, T.
(2009). Role of epigenetics in liver-specific gene transcription, hepatocyte differentiation and stem cell reprogrammation. J. Hepatol. 51, 187–211. Takeuchi, K., Goto, H., Ito, Y., Sato, M., Matsumoto, S., Senba, T., Yamada, H., and Umehara, K. (2015). Dehydroepiandrosterone sulfate and cytochrome P450 inducers alleviate fatty liver in male rats fed an orotic acid-supplemented diet. J. Toxicol. Sci. 40, 181–191. Tarlow, B.D., Pelz, C., Naugler, W.E., Wakefield, L., Wilson, E.M., Finegold, M.J., and Grompe, M. (2014). Bipotential adult liver progenitors are derived from chronically injured mature hepatocytes. Cell Stem Cell 15, 605–618.
Tsai, C.C., Su, P.F., Huang, Y.F., Yew, T.L., and Hung, S.C. (2012). Oct4 and Nanog directly regulate Dnmt1 to maintain self-renewal and undifferentiated state in mesenchymal stem cells. Mol. Cell 47, 169–182. Tu, X., Zhang, H., Zhang, J., Zhao, S., Zheng, X., Zhang, Z., Zhu, J., Chen, J., Dong, L., and Zang, Y. (2014). MicroRNA-101 suppresses liver fibrosis by targeting the TGFbeta signalling pathway. J. Pathol. 234, 46–59. Weigel, D., and Jurgens, G. (2002). Stem cells that make stems. Nature 415, 751–754.
Thiery, J.P., and Sleeman, J.P. (2006). Complex networks orchestrate epithelial-mesenchymal transitions. Nat. Rev. Mol. Cell Biol. 7, 131–142.
Wislet-Gendebien, S., Hans, G., Leprince, P., Rigo, J.M., Moonen, G., and Rogister, B. (2005). Plasticity of cultured mesenchymal stem cells: switch from nestin-positive to excitable neuron-like phenotype. Stem Cells 23, 392–402.
Toh, W.S., Liu, H., Heng, B.C., Rufaihah, A.J., Ye, C.P., and Cao, T. (2005). Combined effects of TGFbeta1 and BMP2 in serum-free chondrogenic differentiation of mesenchymal stem cells induced hyaline-like cartilage formation. Growth Factors 23, 313–321.
Zhau, H.E., He, H., Wang, C.Y., Zayzafoon, M., Morrissey, C., Vessella, R.L., Marshall, F.F., Chung, L.W., and Wang, R. (2011). Human prostate cancer harbors the stem cell properties of bone marrow mesenchymal stem cells. Clin. Cancer Res. 17, 2159–2169.
Stem Cell Reports j Vol. 9 j 247–263 j July 11, 2017 263
Stem Cell Reports, Volume 9
Supplemental Information
DNA Methyltransferases Modulate Hepatogenic Lineage Plasticity of Mesenchymal Stromal Cells Chien-Wei Lee, Wei-Chih Huang, Hsien-Da Huang, Yi-Hsiang Huang, Jennifer H. Ho, MuhHwa Yang, Vincent W. Yang, and Oscar K. Lee
Figure S1
A
B
CD29
CD34
CD44
CD105
CD90
CD117
C
CD73
Sca-1
CD45
CD11b
D
1 μm
Figure S1 Characterization of Mouse Mesenchymal Stromal Cells, Related to Figure 1.
Figure S2
A HD 28
FL83B
Heps
B HD 28
FL83B
E12.5 liver
Heps
Alb Tat Krt18 Dpp4 Hnf1β Hnf3β Hnf4α Asgr1 G6pc Cyp2b10 Αfp Dlk1 Cd24a
C
High
Low
D HD 28
FL83B
Heps
αFP **
800
*
αFP
X1000
600 400 200 0
s
3B
H ep
FL 8
H
D
28
ALB
X100000
NS
5
s H ep
FL 8
H
D
28
3B
0
**
**
H ep s
CK18 600 400 200 30
FL 83 B
β-actin
**
10
X100000
DLK1
CK18
CK18
Integrated Intensity per cell
αFP
ALB
ALB 200 150 100 15
20 10
H D
28
0
Figure S2 Comparison of Mouse MSC-Derived Hepatocyte-Like Cells, FL83B Cell Line and Primary Hepatocytes with Regard to Their Hepatic-Specific Genes Expression, Related to Figure 1.
Figure S3 D mesoderm
CD26high dHeps
osteoblasts
Isolation by FASC with hepatogenic surface marker
endoderm
MSCs
dHeps
MSCs*
dHeps
CD26high dHeps
88.4%
100%
4%
56.3%
5.5%
100%
1.1%
75.5%
0.4%
5.8%
0.8%
41.5%
CD29
dHeps CD26high
ddHeps CD90.2
A
mesoderm
ddHeps
0.2%
4.9%
0.7%
9.1%
0.3 %
2.7%
0.3%
2.9%
92.3%
4%
1.5%
0.1%
CD117 CD45
Hepatocyte markers VE marker
CD106
CD105
Relative mRNA expression
C
CD44
82%
CD73
mBM-MSC surface markers
0%
CD90
4.9%
84%
CD26
CD26high
CD31
CD26low
CD34
B
Hematopoietic lineage surface markers
adipocytes
E reHD 14
reHD 21
reHD 28
PAS
Relative mRNA expression
Morphology
reHD 0
Alb ** **
16 14 12 10 8 6 4 2 0
30 25 20 15 10 5 0
αFP ** **
20
Hnf1β ** **
15 10 5 0
14 12 10 8 6 4 2 0
Dpp4 ** **
HD markers
Tat * **
14 12 10 8 6 4 2 0
αFP
30
HNF3β
Hnf3β * *
20
Hnf4α * **
20
15
15
15
10
10
10
5
5
5
0
0
0
G6pc * **
60 50 40 30 20 10 0
Cyp2b10 * **
Krt18 * **
16 14 12 10 8 6 4 2 0
reOD 0
reOD 7
Tat
**
0.6
ALP
O.D. 405nm
0.5 0.4 0.3 0.2 0.1 0 reOD 0
Cyp2b10 Krt18 Dpp4
0
reOD 7
HD 28
OD markers
reOD 14
Von kossa staining
Alizarin red S
G6pc
reHD 28
Relative mRNA expression
reOD 14
Morphology
reOD 7
Hnf4α
10
F reOD 0
Hnf1β
20 20
Alb
8 7 6 5 4 3 2 1 0
Cdh11 * **
4 3.5 3 2.5 2 1.5 1 0.5 0
Runx2 ** **
8 7 6 5 4 3 2 1 0
Sp7 * **
Bglap1 400 350 300 250 200 150 100 50 0
*
**
8
Cdh11 Sparc Col1α1
4
Col1α2 Runx2
16 14 12 10 8 6 4 2 0
Sparc * *
10 8 6 4 2 0
Col1α1 **
12 10 8 6 4 2 0
Col1α2 * **
20
Spp1 * **
Bglap1
2
Spp1 Sp7
15 10 5 0
1
reOD 14
reOD 14
OD 14
G reAD 7
reAD 14 Relative mRNA expression
Oil Red O
Morphology
reAD 0
AD markers 25 5 4
Fabp4 ** *
10 8
3
6
2 1
4 2
0
0
Fasn * **
12 10 8 6 4 2 0
PPARα ** *
10 8 6 4 2 0
Prdm16 ** **
12 10 8 6 4 2 0
Ucp-1 * **
Fabp4
20
Fasn
15
Prdm16
Pparα Ucp1
10 5 0
reAD 14
Figure S3 CD26high dHeps Derived from mBM-MSCs Can De-Differentiate into Functional MSCs In Vitro, Related to Figure 1 to 3.
AD 14
0hr
D HD 0
HD 0
dHD 7 HD 14
HD 14
dHD 14 HD 21 HD 28
HD 21
dHD 21 dHD 7
dHD 28 dHD 14 dHD 21
300
50
20
ep s
20
ep s
80
H
28 ep s
D H
1500
H
dH
Heps
28
dHD 28
D
HD 28
dH
0
28
1
D
0
0.0
28
2
D
3
H
0.5
1.0
0
1.0 1.5
28
5
D
1.0
H
1.5
Tgfb1
dH
0
HD 0
Snai1
28
X10000
Vim
D
Fn1
0.5
H
0.0
2.5 2.0 1.5 1.0 0.5 0.0
X10000
0
Tjp1
D
0
10
H D H 0 D H 14 D H 21 D d H 28 dH D D 7 dH 1 D 4 dH 2 D 1 2 H 8 ep s
20
20
H
0.5
Epcam
H D H 0 D H 14 D H 21 D dH 28 dH D D 7 dH 1 D 4 dH 2 D 1 2 H 8 ep s
5
40
H D H 0 D H 14 D H 21 D d H 28 dH D D 7 dH 1 D 4 dH 2 D 1 2 H 8 ep s
60
H D H 0 D H 14 D H 21 D dH 28 dH D D 7 dH 1 D 4 dH 2 D 1 2 H 8 ep s
1.5
H D H 0 D H 14 D H 21 D d H 28 dH D D 7 dH 1 D 4 dH 2 D 1 2 H 8 ep s
H D H 0 D H 14 D H 21 D d H 28 dH D D 7 dH 1 D 4 dH 2 D 1 2 H 8 ep s
Relative mRNA expression 2.0
H D H 0 D H 14 D H 21 D dH 28 dH D D 7 dH 1 D 4 dH 2 D 1 2 H 8 ep s
H D H 0 D H 14 D H 21 D dH 28 dH D D 7 dH 1 D 4 dH 2 D 1 2 H 8 ep s
Cdh1
D
D
Integrated Intensity per cell
CDH1
A
H
H
X10000
FN1
1700 900 100 10
Migrated Cells/ field
18hr
SNAI1
Figure S4 B
2.0
Acta2
1.5
1.0
0.0
CDH1 TJP1
4
FN1
0.5
SNAI1
0.0
β-ACTIN
C CDH1
1000
**
** **
500 0
FN1
60
40
*
0
** **
SNAI1
40
30
**
10 0
** **
dHD 28
HD 28 500
400
200
** ** ** ** **
100
0
Figure S4 Mesenchymal-to-Epithelial Transition Was Involved in Reversible Hepatogenic Differentiation of Mouse MSCs, Related to Figure 1.
Oil Red O
Morphology
reAD 0
reAD 0
reAD 7
reAD 7
reAD 14
reAD 14
200
100
0
COL1A2
25
ADIPOQ
4
6
2
0
8
20
6
15
15
4
10
10
2
5
5
0
0
0
**
**
**
**
3
2 1.5
SPARC
25
FABP4
100
1
20
0
0
14
**
14
**
D
BGLAP
14
8
7
0
D
5
*
D
** 2.5
SNAI1
2.0
1.5
1.0
0.5
Migrated Cells/ field
CYP3A4 activity assay
re O D
RUNX2 **
re O
CYP2B6
re A
** *
re O
10 120 80 40 8
7
UPP2
D
** **
re O
FN1
7
TGFβ2
**
D
** 50 40 30 15
**
re H D 0 re HD 14 re HD 28 re HD 42 hH ep s
ALB
**
re A
0
200
0
**
D
UPP2
40
30
50 45 40 35 20 15 10 5 0
** **
0.5
0.0
350
**
250
200
HD 0
**
**
**
HD 42
H H D0 D H 14 D H 21 D H 28 D H 35 D d H 42 dH D d HD 1 7 d H D 24 D 1 hH 2 ep 8 s
H H D0 D H 14 D H 21 D H 28 D H 35 D d 4 dHHD 2 d HD 1 7 d HD 24 D 1 hH 2 ep 8 s
0
re O
10
0
0
20
5
D
10
re O
50
0
0.5
H H D0 D H 14 D H 21 D H 28 D H 35 D d 4 dHHD 2 D dH 17 d HD 24 D 1 hH 2 ep 8 s
H H D0 D H 14 D H 21 D H 28 D H 35 D d 4 dHHD 2 d HD 1 7 d HD 24 D 1 hH 2 ep 8 s
ALB
D
150
Fold of CYP3A4 activity
H H D0 D H 14 D H 21 D H 28 D H 35 D d 4 dHHD 2 d HD 1 7 d HD 24 D 1 hH 2 ep 8 s
Relative mRNA expression 10000 7500 5000 15
re A
1.0
14
1.0
14
1.5
D
0.0
14
5
H D 0 H D 14 H D 28 H D 4 dH 2 D 1 dH 4 D 2 hH 8 ep s
100
D
0.5
7
1.0
D
1.0
re O D
** H H D0 D H 14 D H 21 D H 28 D H 35 D d 4 dHHD 2 D dH 17 dHD 24 D 1 hH 2 ep8 s
1.5
re O
10
re O
0.0
H H D0 D H 14 D H 21 D H 28 D H 35 D d 4 dHHD 2 D dH 17 d HD 24 D 1 hH 2 ep 8 s
0
7
15
re HD 0 re H D 14 re H D 28 re H D 42 hH ep s
42 14
**
D
** D 2 hH 8 ep s
dH
D
28
Albumin production
re A
300
dH
H D
0
14
20
re O
10
*
0
*
D
H H D
**
7
reOD 14 60
re AD
reOD 14 *
re O
reOD 14 **
0
reOD 7 reHD 42 25000 20000 15000 10000 10 8 6 4 2 0
D
0.5
re O
1.5 80
0
reOD 7 HD 42
re AD
TGFβ1
14
CDH2
14
2.0 **
14
0.0
re O D
0.5 **
D
1.5 40
re O
2.0 H D
C
D
reOD 0 reOD 7 dHD28
re A
reOD 0 400
μg/ 106 cells
150
7
reOD 0 reHD 28 reHD 42 **
D
reHD 0 reHD 28 **
7
0.5
Urea production
D
1.0 **
0
ACTA2 NS
re O
CDH1 dHD28
7
0
reHD 0 600
D
0
HD 42
re O
1.5 800
H
0
28 D 4 dH 2 D 1 dH 4 D 2 hH 8 ep s
50 200
D
*
H
Glycogen production
H H D0 D H 14 D H 21 D H 28 D H 35 D dH 42 dH D dHD 17 dHD 24 D 1 hH 2 ep 8 s
800 600 400 20 15 10 5 0 0
**
14
100 **
dHD21
0
0.0
H H D0 D H 14 D H 21 D H 28 D H 35 D d H 42 dH D d HD 1 7 d HD 2 4 D 1 hH 2 ep 8 s
D **
D
**
H
NS
H D
dHD 14
H
250
HD 28
re HD 0 re H D 14 re H D 28 re H D 42 hH ep s
0
14
dHD 7
H H D0 D H 14 D H 21 D H 28 D H 35 D d 4 dHHD 2 D dH 17 dHD 24 D 1 hH 2 ep 8 s
D
200
28 H D 4 dH 2 D 1 dH 4 D 2 hH 8 ep s
H
D
H D H
HD 14
D
0.0
H H D0 D H 14 D H 21 D H 28 D H 35 D d 4 dHHD 2 D dH 1 7 d H D 24 D 1 hH 2 ep 8 s
H H D0 D H 14 D H 21 D H 28 D H 35 D d 4 dHHD 2 dHD 17 dHD 24 D 1 hH 2 ep 8 s
HD 0 dHD21
re O
0.0
H H D0 D H 14 D H 21 D H 28 D H 35 D d H 42 dH D d HD 1 7 d HD 24 D 1 hH 2 ep 8 s
A
re AD
D
I
re O
Morphology
Morphology
dHD 14
Relative mRNA expression
PAS
PAS dHD 7 HD 28
re AD
Morphology
G HD 14
Relative mRNA expression
ALP
μg/ 106 cells
HD 0
Relative mRNA expression
Von kossa
Relative mRNA expression
Figure S5 B HNF4α
1500 1000 500 20 15 10 5 0
TAT
50
E
40
30
20
10 0
F
**
50 0
HD 0 HD 42
dHD 28 hHeps
2.0
VIM
1.5
1.0
Migration assay
300
** NS
**
150
100 50
0 dHD 28
CYP2B6 **
6
4
**
2
0
SP7
1.0
**
4
0.5
0.0
SPP1
20
**
PPARG
80
**
60
40
Figure S5 Hepatogenic Differentiation of Human MSCs Is a Reversible Process, Related to Figure 1 to 3. hHeps
Figure S6
A
Control+LY364947
TGFβ1
TGFβ1+LY364947
7 days
3 days
Control
B Alb 4 3.5
Relative mRNA expression
Hnf1β 5
**
3
**
1.5
*
1
TGFβ1 Alk-5 inhibitor
-
-
+
+
-
-
+
+
-
+
-
+
-
+
-
+
3d
0
**
1.2
-
-
+
+
-
-
+
+
-
+
-
+
-
+
-
+
0.6
3
0.4
2
0.2
1
0
TGFβ1 Alk-5 inhibitor
-
-
+
+
-
-
+
+
-
+
-
+
-
+
-
+
3d
7d
0
-
-
+
+
-
-
+
+
-
+
-
+
-
+
-
+
-
-
+
+
-
-
+
+
-
+
-
+
-
+
-
+
3d
Krt18
**
Dpp4 **
**
**
-
+
+
-
-
+
+
-
+
-
+
-
+
-
+
7d
0
**
4
**
3
**
2
0.5
-
*
5
** 1
*
7d
6
1.5
3d
0
7d
2
5 4
0.5
3d
**
6
*
0
7d
Tat
1
**
1
1
7
0.8
**
1.5
1.5
3d
*
2
*
0.5
7d
** *
2
**
Cyp2b10 1.4
**
2.5
1
0.5 0
3.5
*
2
*
Hnf4α 2.5
**
3
3
**
2
*
4
**
2.5
Hnf3β **
4
*
1
-
-
+
+
-
-
+
+
-
+
-
+
-
+
-
+
3d
7d
0
-
-
+
+
-
-
+
+
-
+
-
+
-
+
-
+
3d
Figure S6 TGFβ1 Induced De-Differentiation in Mouse Hepatocytes, Related to Figure 6.
7d
Figure S7 A
Figure 1C
Figure 1K
Figure 4B
CCND2
250
DNMT1
130
ALB
Figure 4D
250 130 100
DNMT3b
PCNA
DNMT1
130 250
250 β-actin αFP CK18
DNMT1
130
DNMT1
130 100
Longer exposure
130 100
DNMT3a
130 100
DNMT3a
DNMT3b
130 100
Shorter exposure
130
PCNA
55 43 34 26
β-ACTIN
DNMT3b
100
Longer exposure
Short exposure
β-actin
130 100
DNMT3a
55
β-ACTIN
35
Figure 6E
Figure 6I
170 130 100 72 55 43 34 26
ALB
Figure 6O 250
CDH1
130 130
Figure S4B
Figure S2C
DNMT1
DNMT3a
ALB
CDH1
100 CK18
250
CK18
FN1
130
DNMT3b
100 αFP
72 55 43 34 26
αFP
250
SNAI1 TJP1
72 55 43 34 26
β-ACTIN
DLK1
TJP1
FN1
SNAI1 β-ACTIN
130 β-ACTIN
72 55 43 34
β-ACTIN β-ACTIN
B Rabbit IgG
Mouse IgG2a, k
Mouse IgG1, k
MSCs
mHeps
Figure S7 Uncropped Scans of Western Blots and Isotype Control for Immunocytochemical staining, Related to Figure 1, 3, 4, 5, 6 and 7.
Supplemental Figure Legends
Figure S1 Characterization of Mouse Mesenchymal Stromal Cells, Related to Figure 1. (A) Surface phenotype of MSCs. (B) Alkaline-phosphatase staining at 14 days of osteogenic induction in MSCs. Scale bar at 50 μm. (C) Oil Red O staining at 28 days of adipogenic induction in MSCs. Scale bar at 50 μm. (D) Safranin O staining at 28 days of chrondrogenic induction in MSCs. Scale bar at 1 μm. Figure S2 Comparison of Mouse MSC-Derived Hepatocyte-Like Cells, FL83B Cell Line and Primary Hepatocytes with Regard to Their Hepatic-Specific Genes Expression, Related to Figure 1. (A) Cellular morphology of MSC-derived hepatocyte-like cells, FL83B and Heps. Scale bar, 100 μm. (B) A heat-map summarizing relative mRNA levels of hepatogenic-specific genes (n = 3 independent experiments). Result of qRT-PCR is normalized relative to MSCs. (C-D) Results of a representative immunoblot, and (D) immunocytochemical staining for hepatogenic-specific gene expression. Scale bar, 50 μm. “HD 28” refers to MSC-derived hepatocyte-like cells (dHeps). “FL83B” is a hepatocyte cell line from a 15-17 day old fetal mouse. “Heps” refers to mouse primary hepatocytes. All quantitative data are presented as means ± SD. Statistical analyses were performed using Student’s paired t-test, with significance set at P < 0.05. (*P < 0.05; **P < 0.005). Figure S3 CD26high dHeps Derived from mBM-MSCs Can De-Differentiate into Functional MSCs In Vitro, Related to Figure 1 to 3. (A) Schematic illustration of the experimental approach. (B) Representative flow-cytometry plots showing CD26 expression and the gating strategies used to isolate CD26low and CD26high dHeps. (C) qRT-PCR of hepatogenic-specific gene expression in CD26low and CD26high dHeps (n = 4 independent experiments). (D) Surface phenotypic characterization analyzed with flow cytometry during dedifferentiation of CD26high dHeps. CD26high dHeps-derived ddHeps retained hepatogenic (E), osteogenic (F), and adipogenic (G) differentiation potential in vitro (n = 3 independent experiments). Results from Alizarin red S staining are shown (n = 3 independent experiments). Scale bar at 100 μm for phase contrast and staining. All quantitative data are presented as means ± SD. Statistical analyses were performed using Student’s paired t-test, with significance set at P < 0.05. (*P < 0.05; **P < 0.005). “HD 28,” “OD 14,” and “AD 14” represent primary MSCs after lineage commitments to hepatocyte-like cells, osteocytes, and adipocytes, respectively. “reHD 0,” “reOD 0,” and “reAD 0” represent CD26high dHepsderived ddHeps before lineage commitments to hepatocyte-like cells, osteocytes, and adipocytes, respectively. All quantitative data are presented as means ± SD. Statistical analyses were performed using Student’s paired t-test, with significance set at P < 0.05. (*P < 0.05; **P < 0.005). Figure S4 Mesenchymal-to-Epithelial Transition Was Involved in Reversible Hepatogenic
Differentiation of Mouse MSCs, Related to Figure 1. Gene expression of MET-markers analyzed with qRT-PCR (n = 4 independent experiments) (A), immunoblot (B), and immunocytochemical staining (n = 3 independent experiments) (C). Scale bar, 50 μm. (D) Scratch wound-healing assay (upper panel) and trans-well migration assay (lower panels) to measure cell migration capacity (n = 3 independent experiments). Quantification of three independent experiments (lower right panel) from trans-well migration assay (lower left panel). Scale bar, 100 μm. All quantitative data are presented as means ± SD. Statistical analyses were performed using Student’s paired t-test, with significance set at P < 0.05. (*P < 0.05; **P < 0.005). Figure S5 Hepatogenic Differentiation of Human MSCs Is a Reversible Process, Related to Figure 1 to 3. (A) Representative morphological changes and PAS staining of human MSCs during hepatogenic differentiation and de-differentiation. (B) Characterization of hepatogenic-specific genes with qRT-PCR (n = 3 independent experiments). (C) Hepatic functions were measured via albumin production (n = 4 independent experiments), CYP3A4 activity (n = 2, six technical replicates for each point), glycogen storage (n = 2 independent experiments, six technical replicates for each point), and urea production (n = 3 independent experiments). (D) Gene expression of MET/EMT markers analyzed with qRT-PCR (n = 3 independent experiments). (E) Trans-well migration assays of cell migration. The lower left graph shows the quantification of three independent experiments. (F) TGFβ1 expression was measured with ELISA (n = 3 independent experiments); the amount of TGFβ1 from serum has been deducted from results. De-differentiated human MSCs retained hepatogenic (G), osteogenic (H), and adipogenic differentiation (I) potential in vitro (n = 3 independent experiments). Scale bar at 100 μm for phase contrast and staining. All quantitative data are represented as means ± SD. Statistical analyses were performed using Student’s paired t-test, with significance set at P < 0.05. (*P < 0.05; **P < 0.005). Figure S6 TGFβ1 Induced De-Differentiation in Mouse Hepatocytes, Related to Figure 6. (A) Representative morphology of mouse hepatocytes under TGFβ1 treatment. (B) qRT-PCR of hepatogenic-specific genes in mouse Heps under TGFβ1 treatment for 3 and 7 days (n = 4 independent experiments). Scale bar at 100 μm for phase contrast. The + and - symbols respectively indicate treatments with and without TGFβ1 (10 ng/mL) and Alk-5 inhibitor (400 nM). All quantitative data are presented as means ± SD. Statistical analyses were performed using Student’s paired t-test, with significance set at P < 0.05. (*P < 0.05; **P < 0.005). Figure S7 Uncropped Scans of Western Blots and Isotype Control for Immunocytochemical staining, Related to Figure 1, 3, 4, 5, 6 and 7. (A) Uncropped scans of western blots. (B) Rabbit IgG for αFP, CK18 and SNAI1 staining; mouse IgG2a, k for ALB and CDH1 staining; mouse IgG1, k for FN1 staining. Scale bar, 50 μm.
Supplemental Tables Table S1 Functional Enrichment of KEGG Pathways in Hepatogenic Differentiation and DeDifferentiation of Mouse MSCs, Related to Figure 2. KEGG pathway KEGG ID
Up-regulation in HD; down-regulation in dHD
P Value HD 0 to HD 28
HD 28 to dHD 28
mmu00982
Drug metabolism
1.40E-15
1.69E-15
mmu04610
Complement and coagulation cascades
1.91E-15
3.36E-17
4.76E-14
4.08E-13
mmu00980
Metabolism of xenobiotics by cytochrome P450
mmu04142
Lysosome
1.39E-10
3.27E-12
mmu01100
Metabolic pathways
8.93E-10
2.09E-09
mmu00140
Steroid hormone biosynthesis
1.24E-09
1.67E-07
mmu00830
Retinol metabolism
1.07E-08
1.38E-08
mmu03320
PPAR signaling pathway
2.43E-06
2.40E-06
mmu00380
Tryptophan metabolism
9.03E-06
4.65E-04
mmu00480
Glutathione metabolism
1.43E-05
1.91E-06
mmu00590
Arachidonic acid metabolism
1.59E-04
1.58E-04
mmu04146
Peroxisome
6.39E-04
5.84E-05
mmu00350
Tyrosine metabolism
6.66E-04
4.92E-03
mmu00120
Primary bile acid biosynthesis
9.39E-04
1.49E-03
mmu00591
Linoleic acid metabolism
1.45E-03
2.92E-03
mmu00360
Phenylalanine metabolism
1.46E-03
2.43E-03
KEGG pathway KEGG ID
Down-regulation in HD; Up-regulation in dHD
PValue HD 0 to HD 28
HD 28 to dHD 28
mmu03030
DNA replication
2.05E-21
4.62E-23
mmu04110
Cell cycle
2.08E-18
1.39E-17
mmu00240
Pyrimidine metabolism
4.30E-09
2.21E-09
mmu03440
Homologous recombination
1.01E-07
1.88E-07
mmu03430
Mismatch repair
1.01E-06
5.99E-09
mmu00230
Purine metabolism
3.73E-06
1.08E-05
mmu03013
RNA transport
1.45E-05
4.78E-11
mmu03420
Nucleotide excision repair
1.63E-05
2.79E-08
mmu04914
Progesterone-mediated oocyte maturation
5.99E-05
1.33E-04
mmu03410
Base excision repair
7.62E-05
5.07E-08
mmu00670
One carbon pool by folate
3.70E-04
2.32E-04
mmu03040
Spliceosome
2.23E-03
1.32E-08
mmu03008
Ribosome biogenesis in eukaryotes
5.06E-03
7.77E-04
mmu01130
Biosynthesis of antibiotics
6.14E-03
9.28E-02
mmu04114
Oocyte meiosis
6.51E-03
5.38E-03
mmu03018
RNA degradation
3.33E-02
2.03E-03
Table S2 Top 10 Functional Enrichment of KEGG Pathway during MSC-Like Cells Generation from Mouse Hepatocytes, Related To Figure 2.
KEGG ID
KEGG pathway Up-regulation
P Value
mmu04110
Cell cycle
8.55E-25
mmu03040
Spliceosome
3.22E-17
mmu04510
Focal adhesion
4.44E-15
mmu03013
RNA transport
9.02E-12
mmu03030
DNA replication
9.47E-10
mmu04151
PI3K-Akt signaling pathway
3.64E-09
mmu04114
Oocyte meiosis
8.99E-09
mmu04390
Hippo signaling pathway
1.44E-08
mmu00240
Pyrimidine metabolism
9.42E-08
mmu04810
Regulation of actin cytoskeleton
1.25E-07
KEGG ID
KEGG pathway Down-regulation
P Value
mmu01100
Metabolic pathways
9.49E-50
mmu04146
Peroxisome
4.63E-22
mmu00280
Valine, leucine and isoleucine degradation
2.87E-21
mmu04610
Complement and coagulation cascades
1.27E-18
mmu00140
Steroid hormone biosynthesis
1.83E-16
mmu00071
Fatty acid degradation
4.77E-16
mmu00830
Retinol metabolism
6.69E-16
mmu01130
Biosynthesis of antibiotics
5.85E-15
mmu00380
Tryptophan metabolism
1.04E-12
mmu00980
Metabolism of xenobiotics by cytochrome P450
2.56E-12
Table S3 Primer and Primary Antibodies List, Related to Figure 1, 3, 4, 5, 6 and 7. ICC, immunocytochemical assay; IB, immunoblot; FL, flow cytometry; M, mouse; H, Human; NA, not accessed. * The clone numbers are provided for all monoclonal antibodies used in this study. NA means that the antibodies are polyclonal.
Supplemental Experimental Procedures Osteogenic, adipogenic and chondrogenic differentiation We induced osteogenic and adipogenic re-differentiation via treating with osteogenic medium or adipogenic medium for two weeks. The differentiation medium was changed once per week. Mouse osteogenic medium consists of DMDM supplemented with 100 nM dexamethasone, 10 mM β-glycerol phosphate, and 50 μg/mL ascorbate-2-phosphate. Human osteogenic medium consists of IMDM supplemented with 100 nM dexamethasone, 10 mM β-glycerol phosphate, and 50 μg/mL ascorbate-2phosphate. Osteogenesis was assessed by qRT-PCR, alkaline phosphatase, and von Kossa and Alizarin red S stainings. Mouse adipogenic medium consists of DMDM supplemented with 0.5 mM 3-isobutyl1-methylxanthine, 1 μM dexamethasone, 50 μM indomethacin, 5 μg/mL Insulin, and 10% FBS (Gibco®, Thermo Fisher Scientific, Waltham, MA). Human adipogenic medium consists of IMDM supplemented with 0.5 mM 3-isobutyl-1-methylxanthine, 1 μM dexamethasone, 50 μM indomethacin. Adipogenesis was assessed by qRT-PCR and Oil Red O staining. To induce chrondrogenic differentiation, MSCs were transferred into a 15 mL polypropylene tube and centrifuged at 1000 rpm for 5 min, to form a pelleted micromass at the bottom of the tube, and then treated with chondrogenic medium for 4 weeks. Chondrogenic medium consists of high-glucose DMEM supplemented with 0.1 μM dexamethasone, 50 μg/mL ascorbic acid, 100 μg/mL sodium pyruvate, 40 μg/mL proline, 10 ng/mL hTGFβ1, and 50 mg/mL ITS+ premix (BD Bioscience, San Jose, CA). Chondrogensis was assessed by Safranin O staining. SigmaAldrich (St. Louis, MO) supplied dexamethasone, β-glycerol phosphate, ascorbate-2-phosphate, 3isobutyl-1-methylxanthine, indomethacin, Insulin, ascorbic acid, sodium pyruvate, proline and hTGFβ1. Quantitative real-time quantitative polymerase chain reaction RNA was extracted using an RNeasy Mini kit (Qiagen, Valencia, CA) or Trizol reagent (Sigma-Aldrich, St. Louis, MO) and reverse transcribed by Moloney murine leukemia virus reverse transcriptase (Epicentre, Madison, WI). Quantitative real-time quantitative polymerase chain reaction (qRT-PCR) was performed using a ABI StepOnePlus (Applied Biosystems, Waltham, MA). All quantifications were normalized to an endogenous control, Gapdh or RPS27a. Primers for mouse and human genes are listed in Supplemental Table S3. Hepatic functional assays Glycogen storage and urea concentrations were measured by colorimetric assays (Biovision, Milpitas, CA), cytochrome P450 3A4 (CYP3A4) activity was measured by cell-based P450-Glo™ assays (Promega, Madison, WI) through a luminescent method, and albumin secretion was measured by an ELISA assay (Immunology Consultants Laboratory, Inc., Portland, OR) per the manufacturers’ recommendations and analyzed with microplate readers (TECAN, Männedorf, Switzerland). Lowdensity lipoprotein (LDL) uptake by cells was assessed by LDL Uptake Cell-Based Assay kit (Cayman, Ann Arbor, MI) according to the manufacturer’s instructions. MSC-maintenance medium and Step-2 differentiation medium served as background controls for urea assay, CYP3A4 activity, and albumin
secretion assay. Cytochemical, Immunoblot and Immunocytochemical Analysis For glycogen storage determination (PAS staining), cells were fixed in 4% formaldehyde for 10 min, permeabilized with 0.1% Triton X-100 for 10 min, oxidized in 1% periodic acid for 5 min, rinsed 3 times with ddH2O, treated with Schiff’s reagent for 30 min, and rinsed 3 times with ddH2O. For alkaline phosphatase detection (ALP staining), cells were fixed with 4% formaldehyde and treated with NBT/BCIP RT in the dark for 1 hour, and rinsed 3 times in ddH2O. For evaluation of mineralization, cells were fixed with 4% formaldehyde and assayed by von Kossa staining using 1% silver nitrate under ultraviolet light for 45 min, and assayed by Alizarin Red S (ARS) staining using 40 mM Alizarin-red S solution pH 4.1 for 30 min. ARS staining extent was extracted by acetic acid and semi-quantified as previously described (Gregory et al., 2004). For Oil Red O staining, cells were fixed with 4% formaldehyde, treated with 60% isopropanol for 5 min, and stained with Oil Red O for 10 min. Samples were assessed under an optical microscope. For Safranin O staining, paraffin section were deparaffinized by xylene and gradient alcohol treatment and were hydrated by distilled water. Samples were stained with hematoxylin for 3 min and washed in running tap water for 10 min. Samples were rinsed with 1% acetic acid for 10 s and washed in running tap water. Samples were stained in 0.1% Safranin O solution for 5 min. Samples were dehydrated with 95% ethyl alcohol, absolute ethyl alcohol, and xylene, using 2 changes each, 2 min each time. For immunoblotting, cells were lysed in cell lysis buffer containing protease inhibitors and phosphatase inhibitors. Quantified samples were separated by polyacrylamide gels and transferred to PVDF membranes. For immunocytochemical staining, cells were fixed in 4% formaldehyde for 10 min, and permeabilized with 0.1% Triton X-100 for 10 min. After blocking, samples were incubated with primary antibodies against target proteins overnight at 4⁰C, followed by a secondary antibody at room temperature for 1 hour. Sigma-Aldrich (St. Louis, MO) supplied Triton X-100, periodic acid, Schiff’s reagent, formaldehyde, NBT/BCIP, silver nitrate, Alizarin-red S, Oil Red O, Safranin O, cell lysis buffer, protease inhibitor, phosphatase inhibitors, and secondary antibody. Primary antibodies for immunoblotting and immunocytochemical staining were as follows: anti-Fn1 (BD Bioscience, San Jose, CA); anti-CK18 (Epitomic, Burlingame, CA); anti-Cdh1 (Cell Signaling, Danvers, MA). Abcam (Cambridge, UK) supplied anti-Albumin (western blotting), anti-αFP, anti-Cd26, anti-Pcna, anti-Snai1, anti-Tjp1, anti-DNMT1, anti-DNMT3a, and anti-DNMT3b. Sigma-Aldrich (St. Louis, MO) supplied anti-Ccnd2 and anti-β-actin. R&D Systems (Minneapolis, MN) supplied anti-Albumin (for immunocytochemical staining and flow cytomerty). For cell surface phenotyping, cells were detached and stained with fluorescein- or phycoerythrin-coupled antibodies; for cell cycle analysis, cells were detached and stained with propidium iodide per the manufacturer’s recommendations (BD Biosciences, San Jose, CA) and analyzed with FACS Calibur. Primary antibodies for flow cytometry are as follows: BD Biosciences (San Jose, CA) supplied anti-Cd117, anti-Cd44, anti-Cd106, anti-Cd26, anti-Sca-1, antiCd31, anti-Cd90.1, anti-Cd34, anti-Cd45, anti-Cd90.2). eBioscience (San Diego, CA) supplied antiCd29, anti-Cd73, anti-Cd105. Primary antibodies list are listed in Supplemental Table S3 in detail.
Cell growth curve MSCs, MLCs and re-thawed MLCs were seeded at density of 1.8 x 104 and Heps were seeded at a density of 1 x 105 in 6-well plate. The cells were detached and the numbers of cells were counted at each time point. TGFβ1 ELISA assay Human and mouse TGFβ1 secretion were measured by ELISA assay (R&D system, Minneapolis, MN) per the manufacturers’ recommendations and analyzed with microplate readers (TECAN, Männedorf, Switzerland). Scratch wound-healing assay Cell migration was analyzed by a scratch wound-healing assay. Cells were grown to full confluence, and a scratch wound in the monolayer was performed by dragging a 1 mL pipette tip across the layer. Detached cells were washed away with PBS. The remaining cells were cultured in medium, and closure of the wound was observed under a microscope post inflicting the wound for 18 hours. Trans-well migration assay The assay was performed in 24-well transwells (Corning, NY). A total of 2.5 x 104 cells in 0.2 mL medium were seeded to a membrane in the upper well of a trans-well apparatus and allowed to migrate for 18 hours. Cells were fixed with 4% formaldehyde for 10 min and permeablized by 0.1% Triton X100 for 10 min. After removing the non-migrated cells, hematoxylin staining was performed at room temperature for 20 min to visualize the cells passing through the membrane. The migration activity of cells was determined by counting the cells under a microscope, at ×100 magnification. siRNA transfection siRNA oligos were purchased from Invitrogen (Grand Island, NY). Briefly, 50 nM of target siRNA and control siRNA (Stealth RNAi™ siRNA Negative Control, Med GC) were transfected with RNAimax transfection reagent (Invitrogen, Grand Island, NY) at 37⁰C. Cells were washed with PBS and the medium was replaced with differentiation or de-differentiation medium after 18 hours. For Fig. 4B-D, the siRNA against the DNMTs was transfected at day 7 and day 11 of differentiation, and the samples were collected at day 14 and day 28. For Figure 4E, the siRNA against the DNMTs was transfected at day 28 of differentiation with MSC-maintenance medium. Target sequences: mouse DNMT1 siRNA, CAGAGCAGAUUGAGAAGUAUGAUAA, GAGGACAUCCUGUGGUGCACUGAAA, CAGCACAGAGUAUCAGGAUGAUAAA. Microarray and antibody array analysis
mouse and
DNMT3a mouse
DNMT3b
siRNA, siRNA,
MSCs, un-isolated dHep, ddHeps, isolated CD26high Heps and un-isolated MLCs were used to perform microarray and antibody array. For transcriptome analysis, total RNA were hybridized to the Mouse Genome 430 2.0 array (Affymetrix, Santa Clara, CA) at the National Research Progress for Genomic Medicine Microarray and Gene Expression Analysis Core Facility, National Yang-Ming University VYM Genome Research Center, Taiwan. For proteomic analysis, total proteins were analyzed on Lseries mouse antibody array L-308 (RayBiotech, Norcross, GA). All data have been uploaded to the Gene Expression Omnibus (GEO) repository (GSE73617 for microarray and GSE73762 for protein array). One-step Tukey’s BiWeight (TBW) method (Hubbell et al., 2002) were applied to protein expression data to adjust the negative signals. Data were preprocessed and analyzed using packages of Bioconductor (Gentleman et al., 2004) on the R programming language (Ihaka and Gentleman, 1996). The ‘affy’ package (version 1.42.3) was used for data preprocessing and normalization; the ‘affyQCReport’ package (1.42.0) was used to assess quality of microarray data; correlations between microarray samples was shown in heatmap created via the ‘gplots’ package (version 2.16.0). Heatmaps and Hierarchical clustering of gene expression profiles were created by TM4:MeV software (Saeed et al., 2003) (http://www.tm4.org/mev.html). Enrichment analysis was carried out via the Database for Annotation, Visualization and Integrated Discovery (DAVID) (Huang da et al., 2009a, b) and Kyoto Encyclopedia of Genes and Genomes (KEGG). For signaling pathway identification, first 299 intersected genes/proteins were selected and the Pearson’s correlation coefficient (PCC) between the microarray data and protein data for each gene/protein was calculated. Next, the candidates were chosen according to 3 criteria: (1) PCC between DNMT1 and each candidate should be greater than 0.6; (2) gene expression and protein expression of each candidate in MSC should be greater than Heps, as well as DNMT1; (3) PCC between gene expression and protein expression of each candidate should be greater than 0.7. Final, all candidates were arranged by PCC between gene expression and protein expression in descending order. 4 out of 9 selected candidates with high PCC (Tgfb1, Thbs2, Thbs1 and Thbs3) were associated with the TGF-beta signaling pathway. Tumorigenicity assays For tumorigenicity assays, 1 x 106 cells suspended in 100 μL culture medium were injected subcutaneously into 6 week-old male BALB/c nude mice with 6 mice in each group. Animals were monitored regularly for tumor occurrence. An approval prior to the commencement of the experiment from IACUC has been obtained regarding the use of animal. Animal Experimentation In the study, we assured that all animals received humane care according to the criteria outlined in the “Guide for the Care and Use of Laboratory Animals” prepared by the National Academy of Sciences and published by the National Institutes of Health (NIH publication 86-23 revised 1985). All animal studies were approved by the Taipei Veterans General Hospital Institutional Animal Care and Use Committee (IACUC 2014-043).
Statistical analysis All of the quantitative results are shown as the means ± SD from at least three independent measurements. Statistical analyses were performed using Student’s paired t test. Value of P