Abstract
Poorly differentiated tumors usually exhibit phenotypes similar to that of their developmental precursor cells. Tumor cells that acquire the lineage progenitor cells feature usually exploit developmental signaling to potentiate cancer progression. However, the underlying molecular events remain elusive. In this study, based on analysis of an in vitro hepatocyte differentiation model, the maternal factor PGC7 (also known as DPPA3, STELLA) was found closely associated with liver development and tumor differentiation in hepatocellular carcinoma (HCC). Expression of PGC7 decreased during hepatocyte maturation and increased progressively from well-differentiated HCCs to poorly differentiated HCCs. Whole-genome methylation sequencing found that PGC7 could induce promoter demethylation of genes related to development. Pathway-based network analysis indicated that downstream targets of PGC7 might form networks associated with developmental transcription factor activation. Overexpression of PGC7 conferred progenitor-like features of HCC cells both in vitro and in vivo. Mechanism studies revealed that PGC7 could impede nuclear translocation of UHRF1, and thus facilitate promoter demethylation of GLI1 and MYCN, both of which are important regulators of HCC self-renewal and differentiation. Depletion or inhibition of GLI1 effectively downregulated MYCN, abolished the effect of PGC7, and sensitized HCC cells to sorafenib treatment. In addition, we found a significant correlation of PGC7 with GLI1/MYCN and lineage differentiation markers in clinical HCC patients. PGC7 expression might drive HCC toward a “dedifferentiated” progenitor lineage through facilitating promoter demethylation of key developmental transcription factors; further inhibition of PGC7/GLI1/MYCN might reverse poorly differentiated HCCs and provide novel therapeutic strategies.
Similar content being viewed by others
Introduction
It has been well-studied that cancer cells usually exhibited phenotypic similarity with their lineage precursor cells during development [1]. The gene expression signature during normal organ development was reactivated in tumor cells, which substantially contribute to tumor malignant transformation [2]. Hepatocellular carcinoma (HCC) is one of the frequently diagnosed cancers worldwide with an inferior prognosis [3]. Poorly differentiated HCCs retaining the developing signaling of their liver precursors are usually highly aggressive with poor clinical outcome [4]. For example, HCCs which regain the expression of bipotential hepatic progenitor markers such as cytokeratin 7 (CK7) or CK19 is predicted to have an extremely poor prognosis [5]. This evidence convinced us that the mechanisms governing tumor cells lineage reversion should be of critical importance during HCC malignant transformation.
Recently we established an in vitro hepatocyte differentiation model to mimic liver development and HCC progression [6]. The differentiation process comprises four stages—embryonic stem cell, endoderm, liver progenitor cell, and premature hepatocyte stages. Combining clinical liver cancer transcriptomic data, we analyzed different gene expression patterns and mainly focus on actively expressed genes derived from liver progenitors and tumor tissues. Based on the selection criteria, the maternal factor PGC7 arouse our interest since it ranks at the top of the list, and has a unique expression pattern strongly associated with tumor differentiation and poor clinical outcomes.
PGC7/DPPA3 belongs to the developmental pluripotency-associated protein (DPPA) family and was found frequently expressed in pluripotent cells [7]. PGC7 protects DNA methylation pattern in early embryogenesis [8] by binding to either dimethylated histone H3K9 (H3K9me2) [9] or TET2/TET3 [10]. In addition, PGC7 was found to safeguard oocyte methylome by interacting with UHRF1 and DNMT1 [11]. Overexpression of PGC7 enhanced the metastatic ability of melanoma cells coupled with global DNA demethylation [12]. Expression of PGC7 was also found indispensable for the generation of induced pluripotent stem cells (iPSCs) [13]. All of the evidence indicated that PGC7 might be critical in lineage-related tumorigenesis through remodeling DNA methylation patterns.
In the present study, we established a close association of PGC7 expression with liver development and HCC lineage reversion. Overexpression of PGC7 conferred lineage progenitor-like features of tumor cells both in vitro and in vivo. Mechanism studies revealed that PGC7 could facilitate a developmental pluripotency-related transcriptional program through remodeling DNA methylation patterns at key transcription factors. We proposed a novel mechanism in which ectopic expression of PGC7 could promote HCC oncogenic dedifferentiation and maintain an epigenetic status suitable for liver progenitor cells, which further contributed to metastasis and poor prognosis of HCC.
Methods
Cell Lines and HCC clinical specimens
107 pairs of primary HCC specimens and their adjacent non-tumor tissues used for TMA were collected from patients who underwent hepatectomy from HCC at Sun Yat-Sen University Cancer Center (Guangzhou, China). All HCC patients gave written informed consent on the use of clinical specimens for medical research. The samples used in this study were approved by the Committees for the Ethical Review of Research Involving Human Subjects at the Sun Yat-Sen University Cancer Center. Histological examination was carried out by pathologists. Tumors grades were determined according to the American Joint Committee on Cancer (AJCC)/Union for International Cancer Control (UICC) grading system. Investigators were blinded for patient histopathological diagnosis when counting PGC7 positive staining. Human immortalized liver cell lines MIHA and HCC cell lines MHCC-97H, SNU378, H2P, Huh7, Hep3B, SNU449, H2M, SNU475, PLC-8024 were used in this study and tested for mycoplasma contamination. STR DNA profiling analysis was conducted for cell lines authentication. Cells were maintained in high-glucose Dulbecco’s modified Eagle medium (Gibco BRL, Grand Island, NY) supplemented with 10% fetal bovine serum (Gibco BRL) and 1% penicillin/streptomycin. Cells were incubated at 37 °C in a humidified incubator containing 5% CO2.
Establishment of an in vitro hepatocyte differentiation model
The derivation of hESCs and their use for research was approved by the ethical committee of the CITIC-Xiangya Reproductive & Genetic Hospital. In brief, hESCs were cultured at a density of about 2500 cells/cm2 on a feeder layer of mitotically inactivated human embryonic fibroblasts. The culturing medium of hESCs was KO/DMEM medium (Life Technologies) supplemented with 15% knockout serum replacement (SR), 2mM l-glutamine, 2 mM nonessential amino acids, 0.1 mM β-mercaptoethanol, and 4 ng/ml of basic fibroblast growth factor. hESCs were passaged when the confluence was between 50 and 70%, then cultured in RPMI-1640 supplemented with 100 ng/ml activin A (R&D Systems) and 25 ng/ml Wnt3a (R&D Systems) for 3 days until definitive endoderm was specified. To induce hepatic endoderm, cells were growing in KO/DMEM medium supplemented with 25 nm/ml keratinocyte growth factor (R&D Systems) and 2% FBS (Gibco) for two days, followed by culturing in the KO/DMEM medium containing 20% SR, 1% nonessential amino acids, 1 mM glutamine, 0.1 mM 2-mercaptoethanol and 1% DMSO for 4–7 days. To generate mature hepatocytes, cells were cultured in the mature medium consisting of 10% FBS, 10 ng/ml hepatocyte growth factor (R&D Systems), 20 ng/ml Oncostatin M (R&D Systems), and 0.5 μM dexamethasone (R&D Systems) for 7 more days.
In vitro and in vivo functional assays
The in vitro self-renewal capability was assessed by sphere formation assay, drug resistance assay, and ALDH activity assay in cell lines and patient-derived organoids. In vivo tumorigenic ability was investigated in a xenograft mouse model. A detailed description of the HCC organoid establishment and mice model could be found in the Supplementary Materials and Methods section.
Methylation sequencing analysis
The platform and normalization methods of whole-genome methylation sequencing were introduced in the supplementary materials and methods section. Sequences of primers used in Qrt-PCR and bisulfite genomic sequencing were listed in Supplementary Table 1.
Statistics
Statistical analysis was performed using SPSS (version 16.0; SPSS, Chicago, IL). Pearson χ2 test was used to analyze the association of PGC7 expression with clinicopathologic parameters. Kaplan–Meier plot and log-rank test were used for survival analysis. Independent Student t-test was used for most studies as indicated in the figure legends. The Levene or Brown–Forsynth test was used to compare the variance between the two groups. For limiting dilution assay, the tumor-initiating frequency and statistics were calculated using ELDA online platform (http://bioinf.wehi.edu.au/software/elda/). The data are presented as the mean ± standard deviation of three independent experiments. The P-values were denoted as *P < 0.05, **P < 0.01, ***P < 0.001 in all figures.
Results
The distinct gene expression pattern is identified in an in vitro hepatocyte differentiation model
Recently we established a novel in vitro hepatocyte differentiation model, in which human embryonic stem cells (hESCs) were differentiated into hepatocytes with defined culture medium and growth factors. Transcriptomic data were collected for the four developmental stages—embryonic stem cell, endoderm, liver progenitor, and premature hepatocyte, together with two paired HCC clinical tissues. STEM-based trend analysis was performed and as shown in Supplementary Fig. 1a, 13 gene expression patterns were significantly enriched. The genes in profile 31 and 28 shared a similar pattern, in which they have a maximal expression in liver progenitor stage (profile 31) or premature hepatocytes (profile 28), and then markedly decrease or are entirely silenced in adult hepatocytes. Interestingly, hepatic stem cell/progenitor or hepatoblast markers reported by previous studies including E-cadherin [14], Claudin3 [15], TACSTD2 [16] were identified in profile 31, while CD133 [17], AFP [18], DLK1 [19], and NOPE [20] were found in profile 28 (Fig. 1a). Bioinformatics-aided pathway and gene ontology analysis revealed that this cluster of genes was closely associated with cellular development, cell differentiation, and signaling regulating pluripotency of stem cells (Supplementary Fig. 1b), among which PGC7 arouse our interest as it ranks at the top of the list.
PGC7 is highly expressed in pluripotent cells and poorly differentiated tumor tissues
The expression pattern of PGC7 in the hepatocyte differentiation model was further confirmed by qRT-PCR. PGC7 reached its peak expression in the liver progenitor stage and then rapidly dropped after differentiation (Fig. 1b). We analyzed PGC7 expression in Christodoulou’s cohort from the GEO database (GSE27087) and found similar to GLI1, SOX2, NANOG, and OCT3/4, PGC7 expression decreased upon differentiation from iPSCs or ES cells to endoderm (Supplementary Fig. 1c), indicating its role in pluripotency maintenance. The expression of PGC7 was detected by qRT-PCR in 60 pairs of HCC specimens, and it was upregulated significantly in advanced-stage tumors, compared with early-stage tumors and normal liver tissues (Fig. 1c, d). As poorly differentiated tumors always preserve the molecular signature of their developmental precursor cells, we firstly examined PGC7 expression in our HCC cohort divided by the differentiation stage. The result revealed PGC7 was evidently upregulated in poorly differentiated tumor tissues, compared with well-differentiated tumor and normal liver tissues (Fig. 1e). Molecular markers that are activated during liver development or regeneration such as CD133, CK19, AFP, SOX9, EpCAM, CD90, CD105 as well as hepatocyte terminal differentiation markers G6PC and ALB were detected by qRT-PCR in clinical samples. The association study exhibited that PGC7 expression was positively correlated with the developmental markers and negatively correlated with the differentiation markers mentioned above (Fig. 1f and Supplementary Fig. 1d). All of this evidence has linked PGC7 with liver development and HCC lineage reversion.
PGC7 positive cells are markedly associated with poor clinical outcome
To examine the expression and distribution of PGC7 in clinical specimens, tissue microarray (TMA) comprising 107 pairs of HCC tissues and their non-tumor counterparts were used to perform IHC staining. PGC7 expression was detected in about 67.2% (72/107) HCC specimens but almost absent in the majority of normal liver tissues, with only 14 (13.1%) cases showed positive staining. PGC7 was strongly expressed in 0.08 to 4.55% of cancer cells scattered among the tumor tissue. Nevertheless, in normal liver tissue, the percentage of positive cells was around 0.025 to 0.3% (Fig. 1g, h). The clinical pathologic study revealed PGC7 positive staining correlates with poor survival outcome in HCC patients and was significantly associated with advanced tumor stage (Pearson χ2 test, P < 0.05) and relatively high AFP level (Pearson χ2 test, P < 0.05) (Fig. 1h). Co-staining of PGC7 and hepatic progenitor marker CK19 in HCC samples revealed that the majority of PGC7 positive cells were co-expressed with CK19, indicating the marking of undifferentiated cells by PGC7 (Supplementary Fig. 1e).
Furthermore, the frequency of PGC7 positive cells was detected by flow cytometry in 9 HCC cell lines and 1 immortalized liver cell lines. The proportion of PGC7 positive cells ranged from 0.26 to 18.7% (Fig. 1i). Also, we found PGC7 was strongly expressed in liver cancer stem cells (CSCs) enriched from CD133+ populations of Huh7 and Hep3B cells (Supplementary Fig. 1f) as well as tumor spheroids (Supplementary Fig. 1g). Moreover, cryosectioned spheroids produced from PLC-8024 or SNU475 were double-stained with antibodies specific to PGC7 and CD133. We noticed that most of the PGC7 positive cells were co-expressed with CD133 (Fig. 1j and Supplementary Fig. 1h). Collectively, PGC7 may represent and mark a specific population of cancer cells with lineage reversed progenitor features and poor differentiation status, resulting in cancer progression.
PGC7 promotes HCC oncogenic dedifferentiation both in vitro and in vivo
To investigate the role of PGC7 in HCC tumor differentiation and lineage reversion, we established several cell lines with PGC7 overexpressing or silencing. Western blot and immunofluorescence staining confirmed the exogenous expression of PGC7 and its nuclear localization (Fig. 2a and Supplementary Fig. 2a). Lineage progenitor cells usually exhibited enhanced self-renewal ability. We found that overexpression of PGC7 markedly improved both primary and secondary sphere formation ability in MIHA and 97H cell lines (Fig. 2b). Conversely, knockdown of PGC7 significantly decreased the size and number of spheres formed in PLC-8024 cells (Supplementary Fig. 2b, c). Pluripotent markers were found upregulated in 3D spheroids (Supplementary Fig. 2d), among which NANOG, SOX2, and OCT4 exhibited increased expression in PGC7-transfected cells (Supplementary Fig. 2e). Low sensitivity to chemotherapeutic drugs is one of the most important hallmarks of tumor progenitor cells and maintains the major challenge in HCC treatment. To test whether PGC7 affects the sensitivity of HCC cells to chemotherapeutic or targeted drugs, PGC7-transfected cells were treated with cisplatin, 5-fluorouracil (5-Fu), and sorafenib. Cell viability and flow cytometry analysis showed that the frequencies of apoptotic cells were significantly decreased in PGC7-transfected MIHA and 97H cells (Fig. 2c, d), as well as HCC patient-derived organoids (PDOs) (Fig. 2e, f), indicating PGC7 confers chemoresistance property to HCC cells. Conversely, silencing of PGC7 markedly inhibited viability and increased the proportions of apoptotic cells (Supplementary Fig. 2f, g). The effect on drug-sensitivity was also confirmed in vivo using the xenograft mice model. Tumors induced by PGC7-overexpressed cells showed retarded response to sorafenib treatment, compared with Vec-transfected cells (Supplementary Fig. 2h). As high ALDH (aldehyde dehydrogenase) activity serves as a universal marker of lineage precursor cells [21], we then tested ALDH activity in PDOs transfected with PGC7. Compared with Vec-transfected PDOs, PGC7-transfected PDOs demonstrated enhanced ALDH activity (Fig. 2g).
To further confirm whether PGC7 could promote HCC lineage reversion in vivo, empty vector- and PGC7-transfected cells were subcutaneously injected into the left and right dorsal flank of nude mice (n = 6). Results showed the PGC7-transfected cells overtly increased xenograft tumor growth (Fig. 2h). In contrast, no tumor formation was observed for PLC-8024 cells with PGC7 silencing (Supplementary Fig. 2i). We next subcutaneously implanted 500×, 100×, 10×, and 5 × 104 cells into NOD/SCID mice. PGC7-transfected cells demonstrated enhanced tumor-initiating capacity and increased tumor progenitor ratio (Fig. 2i). Moreover, IHC staining was performed in the xenograft tumors and stronger intensity of AFP and CK19 were found in tumors induced by PGC7-overexpressed cells compared with Vec-transfected cells (Supplementary Fig. 2j). To summarize, these data suggest that PGC7 is involved in regulating HCC lineage reversion both in vitro and in vivo.
PGC7 remodels DNA methylation patterns for key developmental transcription factors related to pluripotency and differentiation
Because PGC7 has been reported to preserve DNA methylation pattern in either hyper- or hypomethylated status of early embryos or mice oocyte [9, 11], we propose that PGC7 may exert its function on HCC cells through regulating DNA methylation. Whole-genome bisulfite genomic sequencing (WGBS) was performed to compare the methylation status between Vec- and PGC7-transfected MIHA cells. The total sequence yield was 89.67 and 96.11 gigabases (Gb), with an average read depth of 21.63 and 23.42 per strand (Supplementary Fig. 3a). In each cell type, over 94% of the cytosines in the genome were covered by at least one read depth (Supplementary Fig. 3b). Approximately 24 million and 28 million methylcytosines were detected in Vec- and PGC7-transfected cells respectively, and over 96% of the methylcytosines were in the CG context (Supplementary Fig. 3c).
To examine the difference of methylation distribution between Vec- and PGC7-transfected cells, we firstly compared the methylation level and density throughout various genomic features. The results revealed PGC7 introduction decreased the methylation level of mCG, while increased the relative methylation density of non-CG (mCHG and mCHH, where H = A, C, or T) across all the genomic features (Fig. 3a). The result is consistent with the theory that non-CG methylation primarily exists in pluripotent stem cells and disappeared after differentiation [22]. By looking at the promoter regions (±2 kb to the TSS), lower methylation densities of mCG and mC were observed in PGC7-overexpressed cells, indicating the transcription activity was elevated (Fig. 3b). To identify differentially methylated promoters (DMPs), a sliding window approach was used to scan the genome and 537 genes with promoter demethylation were discovered (Supplementary Table 2). Bioinformatics-aided analysis using Metascape [23] revealed that the potential targets of PGC7 formed clusters associated with cell cycle, CDC42 signaling, and PLK1 pathway (Supplementary Fig. 3d, e). Interestingly, by analyzing the 537 genes with promoter demethylation, we found 44 of them were transcription factors and formed a regulatory network of embryonic development, cell fate commitment, and transcriptional misregulation in cancer (Fig. 3c, d). This evidence implied that PGC7 might activate the oncogenic developmental network through remodeling DNA methylation patterns in HCC.
PGC7 activates key transcription factors to regulate the self-renewal and differentiation of HCC
Expression of transcription factors showing promoter demethylation in PGC7-overexpressed cells was analyzed by qRT-PCR. Among the potential PGC7 targets, GLI1 and MYCN are the most significantly upregulated transcription factors (Supplementary Table 3). GLI1 belongs to the GLI family proteins and has been shown to be involved in cell fate determination [24]. The MYCN gene is a member of the MYC family and was reported to be important in regulating the developmental process and oncogenesis [25]. In our hepatocyte differentiation model, MYCN shared a similar expression pattern with PGC7, and GLI1 was highly expressed in embryonic stem cells and decreased with differentiation, indicating their functions in pluripotency maintenance lineage reversion (Fig. 3e, f).
To explore the role of GLI1/MYCN in HCC, we examined their expression levels in CD133+ CSCs and oncospheres derived from Huh7 cells. Both of them were significantly enriched in liver CSCs (Fig. 4a, b). Their upregulation induced by PGC7 was also confirmed in multiple cell lines and HCC organoids in both mRNA (Fig. 4c, d) and protein level (Fig. 4e), while decreased expression was observed in PGC7-silenced HCC cells (Supplementary Fig. 4a). Notably, MYCN was discovered as a downstream target of GLI1 according to the TRANSFAC database based on CHIP-seq analysis and literature mining. To determine whether PGC7 promotes lineage reversion in a GLI1-dependent manner, shRNA was used to knockdown GLI1 expression in HCC cells. Silencing GLI1 in both Vec- and PGC7-transfected 97H and MIHA cells or in parental PLC-8024 cells remarkably suppressed MYCN expression (Fig. 4f, g and Supplementary Fig. 4b), verifying MYCN as a downstream target gene of GLI1 in HCC. Functional assays including sphere formation (Fig. 4h), ALDH activity (Fig. 4i), and sorafenib-induced apoptosis (Fig. 4j) were then applied and results showed that GLI1 knockdown significantly impaired the progenitor-like features of HCC cells induced by PGC7. We next tested whether the inhibitor of GLI1 (GANT61) exerts a similar effect as GLI1 knockdown. After treatment of 5 μM GANT61 for 3 days, HCC cells demonstrated a reduction of GLI1 and MYCN expression (Fig. 5a and Supplementary Fig. 5a). Fewer spheres (Fig. 5b) and lower ALDH activity (Fig. 5c, d) were observed after GANT61 treatment in HCC cell lines and organoids with PGC7 ectopic expression. To investigate the effect of MYCN on HCC self-renewal, the inhibitor JQ1 was applied in PLC-8024 cells. Treatment of 6 μM JQ1 for 3 days evidently reduced the expression of MYCN but not GLI1 (Supplementary Fig. 5b) and suppressed the sphere formation activity (Supplementary Fig. 5c). These results indicate that PGC7 promotes lineage reversion of HCC via GLI1/MYCN activation.
Combined treatment of sorafenib with GANT61 results in maximal tumor suppression both in vitro and in vivo
To examine the role of GLI1 in regulating sorafenib resistance induced by PGC7, we treated PGC7-overexpressed HCC cells and organoids with both sorafenib and GANT61. Results showed GLI1 inhibition sensitized HCC cells (Fig. 5e) and organoids (Fig. 5f) to sorafenib treatment. Given the crucial role of GLI1 in regulating sorafenib resistance, we evaluated the therapeutic effect of GANT61 treatment and its combined effect with sorafenib in vivo using PGC7-transfected HCC cells. Mice with established tumors were given either vehicle control, GANT61, sorafenib, or the combination of both for 16 days. Compared with the vehicle control group, GANT61, sorafenib, and combined treatment groups all demonstrated reduced tumor size and weight. Besides, the combined treatment group showed maximal tumor suppression (Fig. 5g, h), while no significant difference in mice body weight was observed (Supplementary Fig. 5d). The therapeutic potential of JQ1 was also investigated. Treatment of 10 μM JQ1 for 2 days sensitized HCC cells to sorafenib-induced apoptosis (Supplementary Fig. 5e). Combined treatment of JQ1 with sorafenib in xenograft mice achieved the maximum tumor suppression (Supplementary Fig. 5f, g). These findings indicated that inhibition of PGC7/GLI1/MYCN might reverse poorly differentiated HCCs and provide new therapeutic strategies.
PGC7 interacts with UHRF1 to induce promoter demethylation of key transcription factor
To decipher the mechanism of PGC7 in regulating DNA methylation, we analyzed the interacting proteins of PGC7 using the BioGrid platform (https://thebiogrid.org/). Cytoscape was then used to depict the protein interactions according to evidence value (Fig. 6a). Among the proteins that interacted with PGC7, UHRF1 and UHRF2 arouse our interest with their association with DNA methylation. Immunoprecipitation analysis demonstrated PGC7 binds to UHRF1 rather than UHRF2 or DNMT1 in PGC7-overexpressed 97H and MIHA cells (Fig. 6b and Supplementary Fig. 6a). UHRF1 has been reported to maintain DNA methylation via DNMT1 recruitment [26]. To assess the molecular mechanism more specifically, we investigated the nuclear localization of UHRF1 and DNMT1 in cells synchronized to the S phase (Supplementary Fig. 6b) by dual immunofluorescence staining. We found ectopic expression of PGC7 sequestered UHRF1 from the nuclei to the cytoplasm (Fig. 6c and Supplementary Fig. 6c), while the loss of PGC7-induced UHRF1 nuclear enrichment (Supplementary Fig. 6d). The quantification result exhibited that the percentage of UHRF1 and DNMT1 colocalization in control cells (over 80%) was significantly higher than that in PGC7 expressing cells (<40%) in both cell lines (Fig. 6c, right panel). The nuclear sub-localization of DNMT1 was also investigated and the results showed PGC7 disturbed the localization of DNMT1 on the nuclear speckles (SC35 staining) in both cell lines (Fig. 6d), indicating UHRF1–DNMT1 function was disturbed by PGC7 in HCC cells.
Then, we studied the effect of UHRF1 on PGC7-induced promoter demethylation. Bisulfite genomic sequencing in the promoter region of GLI1 was conducted. The result showed that overexpression of UHRF1 could abolish the effect of PGC7-induced GLI1 promoter demethylation but had no effect on Vec-transfected HCC cells (Fig. 6e and Supplementary Fig. 6e). Conversely, silencing of PGC7 increased the methylation level of the GLI1 promoter (Supplementary Fig. 6f, g). To investigate the role of UHRF1 in regulating PGC7-mediated HCC lineage reversion, sphere formation, and xenograft mice experiments were performed. Fewer spheroids and smaller tumors were generated by 97H-PGC7 cells with UHRF1 ectopic expression, compared with cells without UHRF1 expression (Supplementary Fig. 6h, i). Collectively, these data strongly suggested that binding of PGC7 to UHRF1 disrupts the localization of maintenance DNA methyltransferase DNMT1 in the nucleus, thus results in GLI1 promoter demethylation and HCC lineage reversion (Supplementary Fig. 6j).
GLI1/MYCN activation is correlated with clinical severity and prognosis of HCC patients
As GLI1 and MYCN are two important downstream targets of PGC7 in HCC cells, we proceeded to investigate their co-expression pattern in HCC. An association study showed that PGC7 expression was positively correlated with GLI1 (Pearson R = 0.72; P = 0.000) and MYCN (Pearson R = 0.72; P = 0.000) in our own HCC cohort (Fig. 7a). As expected, the expression of GLI1 was also positively correlated with MYCN both in our HCC cohort (Pearson R = 0.64; P = 0.000, Fig. 7a) and TCGA database (Pearson R = 0.173; P < 0.001, Fig. 7b). GLI1 and MYCN expression was further examined in xenograft tumors. IHC staining (Fig. 7c) and qRT-PCR (Fig. 7d) results demonstrated that both GLI1 and MYCN were markedly increased in xenograft tumors induced by PGC7-transfected HCC cells. Double immunofluorescence staining revealed the colocalization pattern of GLI1 and MYCN in xenograft tumors with PGC7 overexpression (tM1 = 0.52, Supplementary Fig. 7).
Since GLI1/MYCN had a crucial role in HCC lineage reversion, the relationship of GLI1/MYCN activation with HCC progression was further explored. We observed GLI1 was highly expressed in the tumor tissues of HCC patients in the TCGA database (Fig. 7e), and the expression levels were consistent with clinicopathological stages (Fig. 7f). Importantly, GLI1 expression increased progressively from well-differentiated tumors to poor or undifferentiated ones (Fig. 7g) and was associated with poor prognosis (Fig. 7h). Similar observations were achieved when analyzing MYCN expression in the TCGA database. The expression of MYCN was increased gradually with tumor stage development and upregulated in poor or undifferentiated tumors compared with normal tissues (Fig. 7i–k), which was significantly associated with poor clinical outcome (Fig. 7l). In summary, we found that PGC7 is dynamically expressed during liver development and HCC progression. PGC7 was highly expressed in the fetal liver, decreased during hepatocyte maturation, then reactivated in HCCs, and increased progressively from well-differentiated tumors to poorly differentiated tumors. Ectopic expression of PGC7 in HCC induces promoter demethylation of key developmental transcriptional factors and activates their transcription. This further enables HCC cells to acquire lineage progenitor-like features including enhanced self-renewal and drug resistance, which have significantly contributed to the aggressiveness of HCC (Fig. 7m).
Discussion
Accumulating evidence indicated that cancer malignant transformation is associated with the activation of tumor precursor cells [27]. In the process of organ development, lineage-specific regulators remain active and gradually decrease until the terminal differentiation. However, in some cases they were found reactivated in tumor cells, conferring precursor-like features and contributing to cancer progression [28]. Epithelial-to-mesenchymal transition (EMT) is one of the evidence indicating the association between developmental signaling and cancer progression [29, 30]. To investigate the reactivation of key developmental regulators in HCC, we established an in vitro hepatocyte differentiation model, in which hESCs were differentiated into human hepatocytes via several stages. With this model, we could easily identify genes most active in the liver progenitor stage, gradually silenced in mature hepatocytes and reactivated in HCC, indicating their roles in regulating lineage reversion of HCC.
PGC7, also known as DPPA3 and STELLA, was initially identified to be preferentially expressed in primordial germ cells (PGCs) [31]. The chromosomal location of PGC7 is 12p13.31, a region consisting of conserved gene clusters including NANOG and GDF3 [7], which are highly expressed in fetal germ cells and hESCs [32]. Notably, PGC7 expression emerged earlier than endogenous pluripotency genes including Oct3/4, Sox2, Nanog, and Klf4 during somatic cell reprogramming, and is indispensable for the generation of fully reprogrammed iPSCs [13, 33]. All the evidence suggests PGC7 has a critical role in pluripotency maintenance. Here we demonstrate that PGC7 expression reached its peak value in the liver progenitor stage, then gradually disappeared in mature hepatocytes and reactivated in HCC tumor. Our histological and functional studies reveal that PGC7 expression is enriched in liver CSCs and conferred progenitor-like features in HCC.
Based on previous literature, the mechanism of how PGC7 regulates DNA methylation status remains controversial. Some studies conclude that PGC7 protects DNA methylation of imprinted loci and is essential for early embryo development [9, 10]. A recent study demonstrated PGC7 maintains the mice oocyte genome at a hypomethylated status via preventing de novo DNA methylation mediated by DNMT1 and UHRF1 [11]. Here we demonstrated a single-base-resolution map of DNA methylation in liver cells after ectopic expression of PGC7. We found that the average mCG level was decreased by PGC7 across various genomic features including promoter regions. Mechanistically, we observed PGC7 overexpression sequestered UHRF1 from the nucleus, disturbing the correct localization and function of DNMT1. However, several challenges still remain regarding the detailed interactive mode between PGC7 and UHRF1. For example, how PGC7 sequestered UHRF1 from the nucleus when it was predominantly expressed in the nucleus? One possible explanation is PGC7 shuttles between nuclei and cytosol. Upon overexpression, the amount of PGC7 protein is far more abundant than UHRF1, leading to UHRF1 export regardless of the pattern of PGC7 distribution. In this situation, the uncoupling mechanism in the cytosol that releases UHRF1 from PGC7 is yet to be determined.
GLI family proteins are the ultimate effectors of Hedgehog signaling. As zinc-finger transcription factors, GLI1 was reported to induce Hh target genes, having essential roles in embryonic development, tumorigenesis, and maintenance of CSCs [34, 35]. In the present study, we found MYCN is one of the downstream targets of GLI1, and both of them could be activated via promoter demethylation induced by PGC7. Small-molecule inhibitor GANT61 [36] was used to block GLI-mediated transcription, and the self-renewal features of HCC cells were suppressed. It is important to notice that the application of GANT61 increased the sensitivity of Sorafenib both in vitro and in vivo, suggesting that GANT61 might be used as a novel inhibitor to augment the efficacy of sorafenib in HCC treatment.
Above all, evidence from in vitro, in vivo, patient-derived organoids and clinical data analysis suggested an important role of PGC7 in HCC lineage reversion and malignant progression. Further inhibition of PGC7 and the related developmental signaling networks might help to reverse poorly differentiated HCCs and sensitize them to chemotherapeutic or targeted drugs. This might provide a novel therapeutic strategy in HCC treatment.
Data availability
The data sets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
References
Berman JJ. Tumor taxonomy for the developmental lineage classification of neoplasms. BMC Cancer. 2004;4:88.
Kho AT, Zhao Q, Cai Z, Butte AJ, Kim JY, Pomeroy SL, et al. Conserved mechanisms across development and tumorigenesis revealed by a mouse development perspective of human cancers. Genes Dev. 2004;18:629–40.
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. Cancer J Clin. 2011;61:69–90.
Lee JS, Heo J, Libbrecht L, Chu IS, Kaposi-Novak P, Calvisi DF, et al. A novel prognostic subtype of human hepatocellular carcinoma derived from hepatic progenitor cells. Nat Med. 2006;12:410–6.
Mishra L, Banker T, Murray J, Byers S, Thenappan A, He AR, et al. Liver stem cells and hepatocellular carcinoma. Hepatology. 2009;49:318–29.
Liu M, Yan Q, Sun Y, Nam Y, Hu L, Loong JH, et al. A hepatocyte differentiation model reveals two subtypes of liver cancer with different oncofetal properties and therapeutic targets. Proc Natl Acad Sci USA. 2020;117:6103–13.
Payer B, Saitou M, Barton SC, Thresher R, Dixon JP, Zahn D, et al. Stella is a maternal effect gene required for normal early development in mice. Curr Biol. 2003;13:2110–7.
Nakamura T, Arai Y, Umehara H, Masuhara M, Kimura T, Taniguchi H, et al. PGC7/Stella protects against DNA demethylation in early embryogenesis. Nat Cell Biol. 2007;9:64–71.
Nakamura T, Liu YJ, Nakashima H, Umehara H, Inoue K, Matoba S, et al. PGC7 binds histone H3K9me2 to protect against conversion of 5mC to 5hmC in early embryos. Nature. 2012;486:415–9.
Bian C, Yu X. PGC7 suppresses TET3 for protecting DNA methylation. Nucleic Acids Res. 2014;42:2893–905.
Li Y, Zhang Z, Chen J, Liu W, Lai W, Liu B, et al. Stella safeguards the oocyte methylome by preventing de novo methylation mediated by DNMT1. Nature. 2018;564:136–40.
Funaki S, Nakamura T, Nakatani T, Umehara H, Nakashima H, Okumura M, et al. Global DNA hypomethylation coupled to cellular transformation and metastatic ability. FEBS Lett. 2015;589:4053–60.
Xu X, Smorag L, Nakamura T, Kimura T, Dressel R, Fitzner A, et al. Dppa3 expression is critical for generation of fully reprogrammed iPS cells and maintenance of Dlk1-Dio3 imprinting. Nat Commun. 2015;6:6008.
Nitou M, Ishikawa K, Shiojiri N. Immunohistochemical analysis of development of desmin-positive hepatic stellate cells in mouse liver. J Anat. 2000;197:635–46.
Turner R, Lozoya O, Wang Y, Cardinale V, Gaudio E, Alpini G, et al. Human hepatic stem cell and maturational liver lineage biology. Hepatology. 2011;53:1035–45.
Okabe M, Tsukahara Y, Tanaka M, Suzuki K, Saito S, Kamiya Y, et al. Potential hepatic stem cells reside in EpCAM+ cells of normal and injured mouse liver. Development. 2009;136:1951–60.
Rountree CB, Barsky L, Ge S, Zhu J, Senadheera S, Crooks GM. A CD133-expressing murine liver oval cell population with bilineage potential. Stem Cells. 2007;25:2419–29.
Dolle L, Best J, Mei J, Al Battah F, Reynaert H, van Grunsven LA, et al. The quest for liver progenitor cells: a practical point of view. J Hepatol. 2010;52:117–29.
Yanai H, Nakamura K, Hijioka S, Kamei A, Ikari T, Ishikawa Y, et al. Dlk-1, a cell surface antigen on foetal hepatic stem/progenitor cells, is expressed in hepatocellular, colon, pancreas and breast carcinomas at a high frequency. J Biochem. 2010;148:85–92.
Nierhoff D, Levoci L, Schulte S, Goeser T, Rogler LE, Shafritz DA. New cell surface markers for murine fetal hepatic stem cells identified through high density complementary DNA microarrays. Hepatology. 2007;46:535–47.
Minn I, Wang H, Mease RC, Byun Y, Yang X, Wang J, et al. A red-shifted fluorescent substrate for aldehyde dehydrogenase. Nat Commun. 2014;5:3662.
Lister R, Pelizzola M, Dowen RH, Hawkins RD, Hon G, Tonti-Filippini J, et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 2009;462:315–22.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.
Ruiz i Altaba A. The works of GLI and the power of Hedgehog. Nat Cell Biol. 1999;1:E147–8.
Cliff TS, Wu T, Boward BR, Yin A, Yin H, Glushka JN, et al. MYC controls human pluripotent stem cell fate decisions through regulation of metabolic flux. Cell Stem Cell. 2017;21:502–9.
Sharif J, Muto M, Takebayashi S, Suetake I, Iwamatsu A, Endo TA, et al. The SRA protein Np95 mediates epigenetic inheritance by recruiting Dnmt1 to methylated DNA. Nature. 2007;450:908–12.
Garraway LA, Sellers WR. Lineage dependency and lineage-survival oncogenes in human cancer. Nat Rev Cancer. 2006;6:593–602.
Johnston RN, Pai SB, Pai RB. The origin of the cancer cell: oncogeny reverses phylogeny. Biochem Cell Biol. 1992;70:831–4.
Thiery JP, Acloque H, Huang RY, Nieto MA. Epithelial-mesenchymal transitions in development and disease. Cell. 2009;139:871–90.
Singh A, Settleman J. EMT, cancer stem cells and drug resistance: an emerging axis of evil in the war on cancer. Oncogene. 2010;29:4741–51.
Sato M, Kimura T, Kurokawa K, Fujita Y, Abe K, Masuhara M, et al. Identification of PGC7, a new gene expressed specifically in preimplantation embryos and germ cells. Mechanisms Dev. 2002;113:91–4.
Clark AT, Rodriguez RT, Bodnar MS, Abeyta MJ, Cedars MI, Turek PJ, et al. Human STELLAR, NANOG, and GDF3 genes are expressed in pluripotent cells and map to chromosome 12p13, a hotspot for teratocarcinoma. Stem Cells. 2004;22:169–79.
Xu X, Pantakani DV, Luhrig S, Tan X, Khromov T, Nolte J, et al. Stage-specific germ-cell marker genes are expressed in all mouse pluripotent cell types and emerge early during induced pluripotency. PLoS ONE. 2011;6:e22413.
Varnat F, Duquet A, Malerba M, Zbinden M, Mas C, Gervaz P, et al. Human colon cancer epithelial cells harbour active HEDGEHOG-GLI signalling that is essential for tumour growth, recurrence, metastasis and stem cell survival and expansion. EMBO Mol Med. 2009;1:338–51.
Syn WK, Jung Y, Omenetti A, Abdelmalek M, Guy CD, Yang L, et al. Hedgehog-mediated epithelial-to-mesenchymal transition and fibrogenic repair in nonalcoholic fatty liver disease. Gastroenterology. 2009;137:1478.
Lauth M, Bergstrom A, Shimokawa T, Toftgard R. Inhibition of GLI-mediated transcription and tumor cell growth by small-molecule antagonists. Proc Natl Acad Sci USA. 2007;104:8455–60.
Funding
This work was supported by the National Natural Science Foundation of China (Grant No. 81903049, 81772554, and 81802316), High Level-Hospital Program, Health Commission of Guangdong Province, China (Grant No. HKUSZH201902032), Shenzhen Fundamental Research Program (Grant No. JCYJ20180508153249223), Shenzhen Science and Technology Program (Grant No. KQTD20180411185028798), Hong Kong Research Grant Council (RGC) including GRF (17143716 and 767313), Collaborative Research Funds (C7065-18GF and C7026-18GF), Theme-based Research Scheme (T12-704/16-R), and National Key Sci-Tech Special Project of Infectious Diseases (2013ZX10002-011-005). Professor X-Y Guan is Sophie YM Chan Professor in Cancer Research.
Author information
Authors and Affiliations
Contributions
Conceptualization QY and X-YG; methodology QY, YZ, BLL, and XNF; investigation QY, YZ, BLL, XNF, DDY, SL, XW, TLW, LQG, ML, LXJ, TW, YXJ, LL, LZS, YT, NNZ, YL, Y-FY, and SM; writing—original drafts, QY; writing—review and editing, QY and X-YG; funding acquisition, QY and X-YG; resources, Y-FY; supervision, X-YG.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Ethics approval and consent to participate
All HCC patients gave written informed consent on the use of clinical specimens for medical research. The samples used in this study were approved by the Committees for the Ethical Review of Research Involving Human Subjects at the Sun Yat-Sen University Cancer Center.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Edited by M. Sibilia
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yan, Q., Zhang, Y., Fang, X. et al. PGC7 promotes tumor oncogenic dedifferentiation through remodeling DNA methylation pattern for key developmental transcription factors. Cell Death Differ 28, 1955–1970 (2021). https://doi.org/10.1038/s41418-020-00726-3
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41418-020-00726-3
This article is cited by
-
Genomic and transcriptomic profiling of hepatocellular carcinoma reveals a rare molecular subtype
Discover Oncology (2024)
-
LOXL1 promotes tumor cell malignancy and restricts CD8 + T cell infiltration in colorectal cancer
Cell Biology and Toxicology (2024)
-
UHRF1 inhibition epigenetically reprograms cancer stem cells to suppress the tumorigenic phenotype of hepatocellular carcinoma
Cell Death & Disease (2023)
-
SPINK1-induced tumor plasticity provides a therapeutic window for chemotherapy in hepatocellular carcinoma
Nature Communications (2023)
-
PITX2C increases the stemness features of hepatocellular carcinoma cells by up-regulating key developmental factors in liver progenitor
Journal of Experimental & Clinical Cancer Research (2022)