CD73 sustained cancer-stem-cell traits by promoting SOX9 expression and stability in hepatocellular carcinoma
Journal of Hematology & Oncology volume 13, Article number: 11 (2020)
Aberrant AKT activation contributes to cancer stem cell (CSC) traits in hepatocellular carcinoma (HCC). We previously reported that CD73 activated AKT signaling via the Rap1/P110β cascade. Here, we further explored the roles of CD73 in regulating CSC characteristics of HCC.
CD73 expression modulations were conducted by lentiviral transfections. CD73+ fractions were purified by magnetic-based sorting, and fluorescent-activated cell sorting was used to assess differentiation potentials. A sphere-forming assay was performed to evaluate CSC traits in vitro, subcutaneous NOD/SCID mice models were generated to assess in vivo CSC features, and colony formation assays assessed drug resistance capacities. Stemness-associated gene expression was also determined, and underlying mechanisms were investigated by evaluating immunoprecipitation and ubiquitylation.
We found CD73 expression was positively associated with sphere-forming capacity and elevated in HCC spheroids. CD73 knockdown hindered sphere formation, Lenvatinib resistance, and stemness-associated gene expression, while CD73 overexpression achieved the opposite effects. Moreover, CD73 knockdown significantly inhibited the in vivo tumor propagation capacity. Notably, we found that CD73+ cells exhibited substantially stronger CSC traits than their CD73– counterparts. Mechanistically, CD73 exerted its pro-stemness activity through dual AKT-dependent mechanisms: activating SOX9 transcription via c-Myc, and preventing SOX9 degradation by inhibiting glycogen synthase kinase 3β. Clinically, the combined analysis of CD73 and SOX9 achieved a more accurate prediction of prognosis.
Collectively, CD73 plays a critical role in sustaining CSCs traits by upregulating SOX9 expression and enhancing its protein stability. Targeting CD73 might be a promising strategy to eradicate CSCs and reverse Lenvatinib resistance in HCC.
Hepatocellular carcinoma (HCC) ranks as the leading lethal malignant tumor worldwide [1,2,3]. Though great improvements have been made recently, the prognosis of HCC patients remains unfavorable, as evidenced by 25–30% overall 5-year survival and 50–70% recurrence/metastasis rates within 5 years after radical resection [4,5,6,7]. However, the detailed mechanisms underlying liver cancer carcinogenesis and progression remain largely unknown. Hence, the identification of key mediators of HCC development and progression, and the clarification of molecular mechanisms could help improve HCC patient prognosis.
Accumulating evidence has revealed that cellular heterogeneity is maintained in most solid tumors including HCC [8,9,10,11,12]. Cancer stem cells (CSCs) are a small proportion of tumor cells that exhibit capacities of self-renewal and differentiation [13, 14]. In HCC, CSCs initiate tumor development, induce tumor progression and modulate chemotherapy resistance [15,16,17]. Therefore, targeting therapy to eradicate CSCs has the potential to hinder HCC progression. Aberrant signaling activation is involved in maintaining CSC traits in distinct cancers such as gastric and colorectal carcinoma, and HCC . The activation of AKT signaling is implicated in HCC initiation and considered a hallmark to reflect the acquisition of CSC traits [19,20,21,22]. However, it is unclear how key upstream regulators of AKT signaling confer CSC traits to HCC cells.
We previously identified CD73 as a novel indicator of poor prognosis in HCC . Importantly, CD73 promoted HCC progression and metastasis by activating AKT signaling , leading us to hypothesize that CD73 could serve as a novel biomarker for CSCs in HCC. Here, we report that CD73 expression is essential for maintaining CSC traits in HCC cells and that CD73+ HCC cells exhibit substantially greater CSC potential than their CD73– counterparts. Further investigation demonstrated that CD73 not only promotes the transcription of SOX9, but also increases its protein stability in an AKT/glycogen synthase kinase (GSK)3β-dependent manner, which sustained the stemness of HCC cells.
Materials and methods
Patients and clinical samples
Primary HCC samples were obtained from cohort 1 (n = 25, collected from April to June 2018; fresh cancerous tissues were collected and used for sphere-forming assays); cohort 2 (n = 7, collected from July to August 2018; fresh cancerous tissues were collected and used for CD73+ sorting); cohort 3 (n = 212, collected from January to December 2011; used for immunohistochemistry [IHC] staining). All enrolled patients had undergone curative resection at Zhongshan Hospital.
The present study was approved by the Zhongshan Hospital Research Ethics Committee, and all individuals provided their informed consent. HCC diagnosis was based on histopathology examination according to the American Association for Study of Liver Disease guidelines. Follow-up was conducted as previously described  and ended in December 2018. Time to recurrence (TTR) and overall survival (OS) were defined according to previous reports [23,24,25].
Cell lines and animals
HCC cell lines HCCLM3, Hep3B, MHCC97L, and HepG2 were from Yang Xu (Zhongshan Hospital, Fudan University). All cell lines were cultured in DMEM medium supplemented with 10% FBS, 100 unit/ml streptomycin, 100 μg/ml penicillin. Male 4- to 6-week-old non-obese diabetic severe combined immunodeficiency (NOD-SCID) mice were obtained from the Chinese Academy of Science. All model mice were maintained in specific pathogen-free conditions. Humane care of animals was objected to the “Guide for the Care and Use of Laboratory Animals” criteria of the National Academy of Science (National Institute of Health publication 86-23, revised 1985) .
Preparation of primary HCC cells
Primary single HCC cell suspensions were obtained according to our previous study .
Sphere-forming assays were conducted according to our previous study . For HCC cell line culture, cells were seeded at a density of 2000 per well in a 6-well plate. For primary HCC cell culture, tumor cells were seeded at a density of 20,000 per well in a 6-well plate.
RT-PCR and western blot assays
Total RNA extraction was performed using a RNeasy mini kit (Qiagen, Germany), while cDNA synthesis was performed using the Quantitect Reverse Transcription Kit (Qiagen, Germany) as we did previously [23, 27]. mRNA expressions of Target genes as well as internal control were determined with FastStart Universal SYBR Green Master (Roche diagnostic, Germany) and Lightcycle 480 (Roche diagnostic, Germany). ΔCq method (Cqtarget–Cqcontorl) was applied to quantify the relative expression of target genes. Primers used in the present study were listed in Additional file 1: Table S1.
western blot (WB) assays were performed according to our previous studies [23, 28]. Antibodies and corresponding dilutions were listed in Additional file 1: Table S2. All experiments were conducted in triplicate.
The hU6-MCS-CMV-Puromycin lentiviral vectors were purchased from GeneChem Co Ltd. (Shanghai, China). Two distinct shRNAs targeting CD73 and three shRNAs targeting Sox9 were purchased from Merdiobio Co Ltd. (Shanghai, China). shRNA for c-Myc was purchased from Santa Cruz (USA). shRNAs were cloned into lentiviral vectors according to the manufacturer’s instructions. Lentivirus were further transfected into indicated HCC cell lines, and the knockdown efficiencies were validated by WB assays. Both two shRNAs targeting CD73 exerted satisfactory effects, which was in accordance with our previous study . #1 and #2 shRNAs targeting Sox9 were selected for further experiments due to their satisfactory knockdown efficiencies.
For overexpression, expression plasmid of CD73, wild-type SOX9, and mutant SOX9 (T236A) were obtained from Merdiobio Co Ltd. (Shanghai, China). Stable overexpression was verified by WB assays.
Evaluations of self-renewal and differentiated capacities
Differentiated capacity was assessed via culturing sphere cells in medium supplemented with 10% FBS according to our previous study and expression of stemness-associated genes were determined by RT-PCR . EpCAM was set as an internal control reflecting CSC traits due to its universal expression pattern in HCC stem cells . For self-renewal capacity evaluation, serial sphere formation was conducted in continuous 6 passages according to previous studies [27, 30].
In vivo serial dilution xenograft tumor formation
NOD-SCID mice were randomly divided into groups (six per group) and maintained in SPF environment . After quantification of cell number, indicated HCC cells were suspended in a DMEM/MatriGel (Corning, USA) mixture (volume ratio, 1:1), and then injected subcutaneously into the flanks of recipient mice. Tumor formation was monitored as previously described [26, 27], and the incidence was recorded. In vivo tumorigenicity assessment was terminated 12 weeks after injection, at which point mice with no obvious tumor nodules observed at the injection site were considered as negative results.
Isolation of CD73+ cells by magnetic bead cell sorting
For magnetic cell sorting, HCC cells were labeled with primary CD73 antibody (Abcam, USA) followed by incubation with anti-mouse IgG microbeads according to the manufacturer’s instructions. Then, MACS was carried out with miniMACS Starting Kit according to the operation handbook of the manufacturer. To achieve a high purity (> 95%), we conducted positive selection with LS column three times. Aliquots of CD73+ and CD73− populations were evaluated for purity with a FACS Aria II (BD Biosciences). For CD73+ proportion assessment, indicated cells were also labeled with primary CD73 antibody followed by incubation with FITC-conjugated anti-mouse IgG (Abcam, USA). Then FACS analysis was performed with a FACS Aria II. All experiments were conducted in triplicate.
Colony formation assays
For Lenvatinib resistance evaluations, colony formation assays were performed according to our previous studies with tiny modification . Indicated cells were seeded at a density of 4000 cells per well in a 6-well plate.
Drugs and reagents
Lenvatinib was purchased from Selleck and applied as indicated concentrations. Giemsa staining solution was obtained from Sigma. MK-2206 (10 μM) and SC-79 (5 μM) were purchased from Selleck. 10058-F4 (5μM) was purchased from MCE.
Luciferase reporter assays
The luciferase reporter assays were performed as the previous study did . The luciferase activities were assessed with Dual-Luciferase Reporter Assay System (Promega, USA) according to the manufacturer’s instructions at 48 h after transfection. All experiments were performed in triplicate.
Immunohistochemistry IHC staining was conducted with a tissue microarray according to our previous study . The detailed dilution ratio of antibodies used could be seen in Additional file 1: Table S2. Negative controls were included in all assays by omitting the primary antibody. Quantification criteria of IHC staining were the same as the previous study by two independent pathologists .
TCGA database mining
mRNA expression levels of CD73, SOX9, and c-Myc from TCGA database were obtained via StarBase 3.0, and correlations among these three genes were also evaluated based on StarBase 3.0 database .
Protein complexes were precipitated from whole-cell lysates according to our previous study with specific SOX9 antibody, followed by precipitation with protein G beads (Thermo, USA) . To obtain the immunoprecipitates, beads were further boiled in loading buffer followed by centrifugation to obtain supernatants. To assess the expression levels of certain protein within immuoprecipitates, SDS-PAGE and WB assays were performed as previously report described .
Protein stability with cycloheximide
SOX9 protein stability in different modulated HCC cells was assessed with the treatment of cycloheximide (CHX) (10 μg/ml) for 1, 2, and 8 h as the previous study did .
For in vivo ubiquitylation assay, indicated plasmids were transfected into cells along with MG132 treatment (20 μM) for 6 h prior to harvest. Subsequently, protein from whole-cell lysates was extracted followed by sonicated and diluted 10 times with NP-40 lysis buffer according to a previous study . Then, immunoprecipitation was conducted with SOX9 antibody followed by WB assays with indicated antibodies.
Statistical analyses were performed using SPSS 21.0 software (IBM, Chicago, IL, USA). Experimental values for continuous variables were expressed as the mean ± standard error of the mean. The chi-squared test, Fisher’s exact probability tests, and the Student’s t test were used as appropriate to evaluate the significance of differences in data between groups. If variances within groups were not homogeneous, a non-parametric Mann–Whitney test was used. Prognostic value was evaluated by Kaplan–Meier survival curves, log-rank tests, and Cox proportional hazards models. A P value less than 0.05 was considered significant (Additional file 2).
CD73 expression was associated with sphere-forming capacity and was elevated in HCC spheroids
We first evaluated the association between CD73 expression and sphere-forming capacity in 25 fresh resection HCC samples, of which 12 formed spheres within 2 weeks. CD73 protein expression levels were significantly positively associated with the number of spheres formed (R2 = 0.736, P < 0.001, Fig. 1a upper). Similarly, we detected a significant positive correlation between CD73 protein expression and sphere numbers in six HCC cell lines (R2 = 0.316, P = 0.008, Fig. 1a lower). Moreover, HCC spheroid cells derived from both clinical samples and cell lines exhibited higher CD73 protein expression compared with their parental cells respectively according to western blotting assays (Additional file 3: Figure S1A and B).
CD73 expression conferred CSC traits to HCC cells
We knocked down CD73 expression in two CD73-high expression HCC cell lines, Hep3B, and HCCLM3, and overexpressed CD73 in two CD73-low expression cell lines, HepG2, and MHCC97L. After carrying out sphere-forming assays, we found that CD73 knockdown greatly hindered sphere formation (Fig. 1b), whereas CD73 overexpression remarkably increased sphere numbers (Fig. 1c). To validate these results, Hep3B and HCCLM3 spheres were transfected with CD73 short hairpin (sh)RNAs. We observed a significant decrease in sphere number 72 h after transfection in both cell lines (Fig. 1d), and similar results were observed in spheres derived from two clinical samples (Fig. 1e). Limiting dilution xenograft assays showed that CD73 knockdown significantly reduced tumor initiation and tumorigenic cell frequency compared with control cells (Fig. 1f).
Serial sphere formation assays revealed that CD73 knockdown also greatly reduced the ability of cells to self-renew (Fig. 1g), whereas CD73 overexpression achieved the opposite effect (Fig. 1h). Three rounds of serial passaging were performed to investigate dynamic changes in CD73 mRNA expression, and the expression of EpCAM as a universal CSC marker was measured as an internal control to reflect CSC traits . CD73 mRNA expression in Hep3B and HCCLM3 cells was significantly upregulated in sphere cells and showed a notable decrease following 10% FBS-induced differentiation (Additional file 4: Figure S2A). Consistently, similar dynamic change patterns in CD73 mRNA expression were detected in cells derived from two clinical samples (Additional file 4: Figure S2B).
Additionally, CD73 knockdown remarkably sensitized HCC cells to Lenvatinib treatment (Fig. 1i), while CD73 overexpression induced Lenvatinib resistance (Fig. 1j). Collectively, these data suggest that CD73 promoted the self-renewal of HCC cells and in vivo tumor propagation.
CD73 is essential for the HCC stemness-associated phenotype
RT-PCR assays indicated that CD73 knockdown significantly reduced the mRNA expression of stemness-associated genes such as EpCAM, Nanog, SOX2, Oct4, SOX9, and c-Myc, while increasing the expression of albumin and cytokeratin 8 (CK8), which were considered mature liver cell markers  (Fig. 2a). Conversely, CD73 overexpression achieved the opposite effects (Fig. 2b). Western blotting assays further confirmed these findings (Fig. 2c, d). Notably, after transfection with CD73 shRNAs, the expression of EpCAM and SOX9 decreased in a time-dependent manner in Hep3B and HCCLM3 sphere cells, while CK8 showed increased expression (Fig. 2e, f). Consistently, the transfection of CD73 shRNAs in spheres derived from two clinical HCC samples exerted similar dynamic change patterns of EpCAM, SOX9, and CK8 to cell lines (Fig. 2g, h). These data confirmed that CD73 is essential for maintaining stemness-associated molecular phenotypes in HCC.
CD73+ HCC cells possess CSC characteristics
We next examined CD73 expression by FACS in six HCC cell lines. Consistent with western blotting findings, HCCLM3, MHCC97H, and Hep3B cells possessed a high percentage of CD73+ cells, while HepG2, and MHCC97L cells had relatively low CD73+ cell numbers (Additional file 5: Figure S3). MACS-purified CD73– and CD73+ HCC cells underwent sphere formation assays, and CD73+ cells derived from both HCC cell lines and fresh clinical samples were observed to form more spheres than their CD73– counterparts (Fig. 3a, b). Moreover, CD73+ cells exhibited a stemness molecular phenotype, evidenced by the higher expression of stemness-associated genes, while CD73– fractions demonstrated higher expression of albumin and CK8, indicating a mature hepatocyte-like phenotype (Fig. 3c, d). We further purified CD73+ and CD73– cells via MACS, then cultured them in DMEM medium supplemented with 10% FBS for 7 days. CD73 expression was then determined in each subpopulation by FACS. CCK8 assays revealed that both CD73+ and CD73– fractions could proliferate in the given medium. However, CD73+ cells dramatically decreased in number after 3 days, then almost reverted to the presorting level, while CD73– cells maintained their low number after 7 days; this suggested that CD73– cells arose from CD73+ cells and not vice versa (Fig. 3e). CD73+ cells also exhibited greater resistance capacities to Lenvatinib than their CD73– counterparts (Fig. 3f). After purified cells were subcutaneously inoculated into NOD/SCID mice, we observed a difference in tumor incidence between CD73+ and CD73– cells (Fig. 3g).
We further knocked down CD73 expression in purified CD73+ cells, which significantly reduced the expression of stemness-associated genes, while promoting CK8 expression (Fig. 3h). Moreover, CD73 knockdown also greatly hindered the sphere formation capacities of CD73+ cells (Fig. 3i), and induced sensitization to Lenvatnib (Fig. 3j). Taken together, these data suggested that CD73+ cells possess CSC characteristics, and that CD73 expression is essential for sustaining stemness in the CD73+ subpopulation.
CD73 promotes CSC traits by up-regulating SOX9
The above data not only revealed that SOX9 mRNA expression was reduced more than that of other stemness-associated genes by CD73 knockdown, but also showed the highest increase after CD73 overexpression in HCC cells (Fig. 2a, b). Interestingly, SOX9 was reported to be regulated by the AKT signaling pathway , which is activated by CD73 according to our previous study . Therefore, we speculated that SOX9 could be a key downstream regulator of CD73. To validate this, we first depleted SOX9 expression in CD73-high expression Hep3B cells. EpCAM mRNA expression was significantly decreased, whereas CK8 mRNA expression was greatly upregulated following SOX9 knockdown, which mimicked the effects of CD73 knockdown (Fig. 4a). More importantly, SOX9 knockdown effectively abolished the effects of CD73 on EpCAM and CK8 mRNA expression in MHCC97L cells (Fig. 4b). These findings were further confirmed by western blotting (Fig. 4c). Similarly, SOX9 knockdown mimicked the inhibition effects of CD73 knockdown on sphere formation and Lenvatinib resistance capacities in Hep3B cells and abolished the promotional effects of CD73 overexpression on sphere formation and Lenvatinib resistance capacities in MHCC97L cells (Fig. 4d, e), indicating the potential role of SOX9 in CD73-mediated CSC traits.
We further purified CD73+ cells from two different HCC cell lines, followed by SOX9 shRNA transfection. Western blotting assays indicated that SOX9 knockdown reduced EpCAM expression in CD73+ cells almost to the levels of CD73– cells (Fig. 4f). Further investigation demonstrated that SOX9 knockdown restrained sphere formation and Lenvatinb resistance capacities of CD73+ cells (Fig. 4g, h). Collectively, our data demonstrated that CD73 controlled CSC traits by upregulating SOX9 expression.
CD73 upregulates SOX9 by AKT signaling
We next explored whether CD73 upregulates SOX9 by activating the AKT signaling pathway. First, we treated CD73-high expression cells with the AKT antagonist MK-2206, which greatly reduced pAKT (Ser473) and pGSK3 (Ser9) levels. Importantly, SOX9 expression was also dramatically reduced by AKT inhibition (Fig. 5a). By contrast, treatment with the AKT agonist SC-79 significantly increased SOX9 expression in CD73-low expression cells (Fig. 5a). Additionally, as expected, AKT antagonist treatment significantly restrained sphere formation and Lenvatinib resistance capacities of CD73-high expression cells, whereas the AKT agonist achieved the opposite effects in CD73-low expression cells (Fig. 5b, c). Further study showed that re-activation of AKT by SC-79 successfully rescued SOX9 expression after inhibition by CD73 knockdown, along with restored sphere formation and Lenvatinib resistance capacities (Fig. 5d–f). Conversely, the inactivation of AKT signaling effectively abolished SOX9 upregulation caused by CD73 overexpression, together with abolished effects on sphere formation and Lenvatinib resistance capacities (Fig. 5g–i). These data imply that CD73 modulated SOX9 expression in an AKT-dependent manner.
We then purified CD73+ and CD73– cells, and treated CD73+ cells with MK-2206 and CD73– cells with Sc-79. We found that AKT inhibition reduced the expression of SOX9 in CD73+ cells, whereas AKT activation greatly increased SOX9 expression in CD73– cells (Fig. 5j). Further, as expected, AKT inhibition attenuated sphere formation and Lenvatinib resistance capacities in CD73+ cells, while AKT activation achieved the opposite effects in CD73– cells (Fig. 5k, l). Finally, to confirm the critical role of SOX9 in the regulatory process, we knocked down SOX9 expression in SC-79-treated CD73-low expression cells. This substantially attenuated the effects of AKT activation on CSC traits, as evidenced by decreased EpCAM expression, sphere formation capacity, and Lenvatinib resistance potentials (Fig. 5m–o). Hence, our data clearly demonstrated that CD73 modulates SOX9 expression by activating AKT signaling.
CD73 activates SOX9 transcription through c-Myc
An increase in SOX9 expression might reflect activated SOX9 transcription and/or enhanced SOX9 protein stability. Indeed, SOX9 mRNA levels were positively correlated with CD73 expression (Fig. 3a). c-Myc serves as a key downstream target of AKT signaling, inducing the transcription of various genes upon AKT activation . Importantly, the results from our cohort and the TCGA database showed that CD73 mRNA expression was positively correlated with both SOX9 and c-Myc, and that SOX9 expression was significantly correlated with that of c-Myc (Fig. 6a). Thus, we hypothesized that CD73 triggered SOX9 transcription via c-Myc.
In support of this, we first observed notable SOX9 downregulation at both the mRNA and protein level following the silencing of c-Myc in CD73-high expression cells (Fig. 6b upper). Moreover, the treatment of cells with the specific c-Myc inhibitor 10058-F4 significantly reduced SOX9 expression (Fig. 6b lower). Interestingly, we found that either c-Myc silencing or inhibitor treatment abolished the effect of CD73 overexpression on promoting SOX9 expression in CD73-low expression cells (Fig. 6c). Similarly, c-Myc silencing or inhibition significantly reduced SOX9 mRNA and protein levels in CD73+ fractions (Additional file 6: Figure S4A and B). Luciferase reporter assays demonstrated significantly higher SOX9 promoter activity when SOX9 and c-Myc were co-transfected into 293 T cells (Additional file 6: Figure S4C). Moreover, CD73 overexpression markedly activated the SOX9 promoter, while this effect was almost completely abolished by c-Myc silencing in CD73-low expression cells (Fig. 6d). Similar results were observed in sorted CD73– cells (Additional file 6: Figure S4D). Together, these data imply that CD73 mainly promotes SOX9 transcription through c-Myc.
CD73 prevents SOX9 ubiquitination and proteasome degradation by inhibiting GSK3β
Active GSK3β phosphorylates SOX9, resulting in SOX9 degradation via the ubiquitin–proteasome pathway . Because GSK3β activity was greatly hindered by CD73 through activating AKT , we next explored whether CD73 used this strategy to enhance SOX9 protein stability. We observed that the proteasome inhibitor MG132 partially restored SOX9 expression in CD73-silenced cells (Fig. 6e), suggesting that CD73 also contributes to the stability of SOX9 protein. Next, we detected a significant reduction in SOX9 protein stability when CD73 was silenced or AKT signaling was inhibited, as evidenced by a shorter SOX9 half-life. However, either an AKT agonist or GSK3β inhibitor effectively restored SOX9 stability caused by CD73 knockdown. By contrast, CD73 overexpression or the use of a GSK3β inhibitor greatly prolonged the SOX9 half-life in CD73-low expression cells, whereas an AKT inhibitor attenuated the effect of CD73 overexpression (Fig. 6f, Additional file 6: Figure S4E).
A previous study reported that residue T236 was critical for GSK3β-induced phosphorylation, resulting in SOX9 degradation . To confirm this, we co-transfected wild-type SOX9 (WT-SOX9) or mutant SOX9 (SOX9-T236A) with GSK3β into 293 T cells. We found that GSK3β markedly promoted WT-SOX9 ubiquitination levels, resulting in a reduction in SOX9 expression. However, GSK3β failed to induce SOX9-T236A ubiquitination, also reducing SOX9 expression (Additional file 6: Figure S4F). Immunoprecipitation followed by western blotting was performed to determine SOX9 ubiquitination levels under different treatments. CD73 knockdown or AKT inhibition greatly enhanced the interaction between SOX9 and GSK3β/FBW7, resulting in increased SOX9 ubiquitination levels. However, AKT activation or GSK3β inhibition abolished the effects of CD73 knockdown on SOX9 ubiquitination (Fig. 6g). By contrast, CD73 overexpression and GSK3β inhibition markedly reduced the interaction between SOX9 and GSK3β/FBW7, leading to lower SOX9 ubiquitination levels, while AKT inhibition attenuated the effects of CD73 overexpression (Fig. 6h). Collectively, our data demonstrated that CD73 both promoted SOX9 transcription and enhanced SOX9 protein stability by activating AKT signaling (Fig. 6i).
The combined analysis of CD73 and SOX9 is a promising approach to evaluate prognosis in HCC
We observed positive correlations between CD73 and SOX9 (r = 0.583, Fig. 7a, b), as well as between c-Myc and SOX9 and CD73 and c-Myc, which confirmed our earlier results based on mRNA expression. High expression of CD73 or SOX9 indicated a significantly worse HCC prognosis, as evidenced by shorter TTR and OS (all P < 0.001, Fig. 7c, d). Importantly, multivariate Cox regression identified both CD73 and SOX9 as independent predictors for shorter TTR [CD73: HR 2.251, 95% confidence interval (CI) 1.533–3.304, P < 0.001; SOX9: HR 1.697, 95% CI 1.138–2.478, P = 0.009] as well as OS [CD73: HR 1.981, 95% CI 1.208–2.998, P = 0.001; SOX9: HR 1.734, 95% CI 1.139–2.642, P = 0.010]. Therefore, we speculated that the combined analysis of CD73 and SOX9 would provide a more powerful tool to predict the prognosis of HCC patients. To validate this, we stratified HCC patients into four subgroups: (i) CD73-high and SOX9-high (n = 62); (ii) CD73-high and SOX9-low (n = 37); (iii) CD73-low and SOX9-high (n = 25); and (iv) CD73-low and SOX9-low (n = 88). TTR was significantly shorter in group i than groups ii (P = 0.017), iii (P = 0.003), and iv (P < 0.001). Patients in group ii (P = 0.005) or group iii (P = 0.013) also experienced a shorter TTR than those in group iv. However, TTR showed no significant difference between groups ii and iii (P = 0.819, Fig. 7e). Similarly, OS was significantly shorter in group i than groups ii (P = 0.017), iii (P = 0.015), and iv (P < 0.001). Patients in group ii (P = 0.025) or group iii (P = 0.016) experienced a shorter OS than those in group iv. However, OS did not differ significantly between groups ii and iii (P = 0.925, Fig. 7e). ROC curve analysis demonstrated that combining the analyses of CD73 and SOX9 exerted a greater power in distinguishing patients with worse prognosis (AUC-ROC; TTR: CD73 0.675, SOX9 0.695, combined 0.741; OS: CD73 0.642, SOX9 0.668, combined 0.701; Fig. 7f). Collectively, our data suggest that this combination is a promising tool for predicting prognosis in HCC patients.
Tumor heterogeneity is widely accepted, and the CSC model has been verified in various types of solid tumors including HCC . Here, we report CD73 as a novel surface marker of CSC in HCC. We found that CD73 expression was essential for the capacities of self-renewal, differentiation, and the generation of new tumors. Importantly, we showed that the CD73+ cell fraction, but not CD73– cells, exerted typical CSC traits. Further investigations revealed that CD73 promoted CSC traits through dual AKT-dependent mechanisms: activating SOX9 transcription via c-Myc, and preventing SOX9 degradation by inhibiting GSK3β. Moreover, CD73 expression positively correlated with SOX9, suggesting that a combination of these two biomarkers has the potential to precisely predict the prognosis of HCC patients.
CD73 is a well-known surface marker for identifying mesenchymal stem cells . Recent investigations revealed the potential for CD73+ photoreceptor precursors in future cell replacement therapy , and showed that CD73 could mark a multipotent stromal population . CD73 was also reported to mark a fraction of mammary cells endowed with lineage plasticity . These findings indicate the importance of CD73 in regulating stemness while increasing evidence suggests a role for CD73 in regulating CSC traits in several types of solid tumors [42, 43]. However, there was no specific data to connect CD73 with CSCs in HCC.
We previously revealed CD73 to be a critical oncogene for HCC progression that could trigger epithelial–mesenchymal transition , which is considered to confer stemness traits to cancer cells . This strongly suggested that CD73 is a potential marker for identifying CSCs in HCC. Here, we found that CD73 was not only enriched in HCC spheres, but that it was also essential for the sphere formation capacity of HCC cells. Interestingly, we found CD73 also mainly depended on its enzyme activity to promote CSC traits in HCC (Additional file 7: Figure S5). To further determine whether CD73 was an ideal marker for liver CSC, we used MACS sorting and employed CD73 expression modulation to show that CD73 expression was essential for self-renewal, differentiation, and in vivo tumor propagation. Consistently, CD73+ HCC cells exhibited a significantly greater stemness potential than their CD73– counterparts. Notably, CD73 was also closely associated with an undifferentiated phenotype, as evidenced by the increased expression of stemness genes. Moreover, single-cell sorting showed that CD73+ HCC cells could differentiate into CD73– cells, but the reverse did not occur. Collectively, our data clearly show that CD73 is a novel, ideal, and promising surface marker for HCC.
Lenvatinib was recently approved as a novel molecular target regimen for first-line therapy in HCC ; this has resulted in a meaningful improvement in several second endpoints including disease-free progression . However, resistance to Lenvatinib is observed in clinical practice , which greatly hinders the therapeutic effects and contributes to poor prognosis. Here, we found that overexpression of CD73 conferred HCC cells with significant resistance to Lenvatinib. Moreover, purified CD73+ cells exhibited outstanding drug resistance compared with their CD73– counterparts. Colony formation assays of the effects of CD73 on Lenvatinib resistance reflected the survival potential of HCC cells under Lenvatinib treatment. Thus, our data demonstrated that CD73 is a critical regulator contributing to Lenvatinib resistance; therefore, indicating that targeting of CD73 or elimination of CD73+ cells is a promising strategy for overcoming Lenvatinib resistance and prolonging OS. Intriguingly, CD73 is also critical for promoting other multiple kinase inhibitor resistance, such as sorafenib or Cabozantinib (Additional file 8: Figure S6), which suggests the potential role of CD73 in sustaining multiple kinase resistance of HCC, and more investigations are needed to confirm our findings.
Accumulating evidence has indicated that aberrant activation of AKT signaling is a vital process for maintaining CSC traits in HCC [11, 20, 46]. We previously found that CD73 effectively activated AKT signaling via the Rap1/P110β cascade . Because CD73 appears to be a novel biomarker for CSCs in HCC, we speculated that CD73 also regulated CSC traits by activating AKT. In support of this, our data showed that AKT inactivation successfully hampered CSC traits in CD73-high expression cell lines or sorted CD73+ cells, whereas AKT activation triggered by SC-79 achieved the opposite effects in CD73-low expression cell lines or sorted CD73– cells. These results clearly confirmed our hypothesis and demonstrated the involvement of AKT signaling in CD73-regulated CSC traits.
SOX9 is an important transcription factor for regulating high mobility group box DNA binding and transactivation domains . Previous studies implicated the contributions of SOX9 to stemness characteristics in HCC [47, 48]. However, the regulation of SOX9 transcription and protein stability was not fully investigated in HCC. Here, we showed that SOX9 is a critical downstream regulator of CD73 in HCC. Intriguingly, CD73 upregulated SOX9 expression via two distinct mechanisms: promoting SOX9 transcription by c-Myc and promoting SOX9 protein stability by inhibiting GSK3β activity, which resulted in an enhancement effect to further facilitate stemness characteristics. Although SOX9 was previously shown to activate Wnt signaling  which triggered Myc expression, the regulatory relationship between SOX9 and c-Myc remained vague in HCC. Here, we observed that c-Myc knockdown dramatically inhibited SOX9 mRNA expression, whereas SOX9 silencing had no significant effect on c-Myc expression, indicating the vital role of c-Myc in regulating SOX9 transcription in HCC. However, we failed to identify specific DNA-binding seeds of c-Myc on the SOX9 promoter (data not shown). Therefore, we speculate that c-Myc serves as a coregulator for other factors to indirectly enable SOX9 transcription. We plan to elucidate the mechanism underlying c-Myc-mediated SOX9 transcription in future work. It is conceivable that SOX9 is phosphorylated at the A236 residue by GSK3β, resulting in proteasomal degradation . Because GSK3β is mainly inactivated by AKT , we further inferred that CD73 also enhances SOX9 protein stability via AKT activation. Here, we observed that CD73 could act as a “switch” to control SOX9 ubiquitination following proteasome degradation through inhibiting GSK3β by activating AKT signaling. Importantly, our findings also implicate a new regulatory mechanism of GSK3β in restricting CSC traits in HCC.
In summary, our present study revealed CD73 as a novel surface marker for identifying CSC in HCC. CD73 sustained CSC traits both by upregulating SOX9 expression and maintaining its protein stability. Our findings suggest a potential target for overcoming Lenvatinib resistance and a more comprehensive understanding of the regulatory mechanism involved in liver CSC.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Cancer stem cell
Glycogen synthase kinase 3-beta
Magnetic bead cell sorting
Non-obese diabetic severe combined immunodeficiency
The cancer genome atlas
Time to recurrence
Yin Z, Dong C, Jiang K, Xu Z, Li R, Guo K, et al. Heterogeneity of cancer-associated fibroblasts and roles in the progression, prognosis, and therapy of hepatocellular carcinoma. J Hematol Oncol. 2019;12:101.
Jayachandran A, Dhungel B, Steel JC. Epithelial-to-mesenchymal plasticity of cancer stem cells: therapeutic targets in hepatocellular carcinoma. J Hematol Oncol. 2016;9:74.
Lu C, Rong D, Zhang B, Zheng W, Wang X, Chen Z, et al. Current perspectives on the immunosuppressive tumor microenvironment in hepatocellular carcinoma: challenges and opportunities. Mol Cancer. 2019;18:130.
Nault JC, Ningarhari M, Rebouissou S, Zucman-Rossi J. The role of telomeres and telomerase in cirrhosis and liver cancer. Nat Rev Gastroenterol Hepatol. 2019;16:544–58.
Calderaro J, Ziol M, Paradis V, Zucman-Rossi J. Molecular and histological correlations in liver cancer. J Hepatol. 2019;71:616–30.
Wu Q, Zhou L, Lv D, Zhu X, Tang H. Exosome-mediated communication in the tumor microenvironment contributes to hepatocellular carcinoma development and progression. J Hematol Oncol. 2019;12:53.
Kanwal F, Singal AG. Surveillance for hepatocellular carcinoma: current best practice and future direction. Gastroenterology. 2019;157:54–64.
Yin X, Zhang BH, Zheng SS, Gao DM, Qiu SJ, Wu WZ, et al. Coexpression of gene Oct4 and Nanog initiates stem cell characteristics in hepatocellular carcinoma and promotes epithelial-mesenchymal transition through activation of Stat3/snail signaling. J Hematol Oncol. 2015;8:23.
Fan Z, Duan J, Wang L, Xiao S, Li L, Yan X, et al. PTK2 promotes cancer stem cell traits in hepatocellular carcinoma by activating Wnt/beta-catenin signaling. Cancer Lett. 2019;450:132–43.
Cao J, Zhao M, Liu J, Zhang X, Pei Y, Wang J, et al. RACK1 promotes self-renewal and Chemoresistance of Cancer stem cells in human hepatocellular carcinoma through stabilizing Nanog. Theranostics. 2019;9:811–28.
Gu Y, Wei X, Sun Y, Gao H, Zheng X, Wong LL, et al. miR-192-5p silencing by genetic aberrations is a key event in hepatocellular carcinomas with Cancer stem cell features. Cancer Res. 2019;79:941–53.
Majumdar A, Curley SA, Wu X, Brown P, Hwang JP, Shetty K, et al. Hepatic stem cells and transforming growth factor beta in hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol. 2012;9:530–8.
Lv H, Lv G, Han Q, Yang W, Wang H. Noncoding RNAs in liver cancer stem cells: the big impact of little things. Cancer Lett. 2018;418:51–63.
Gordeeva O. Cancer-testis antigens: unique cancer stem cell biomarkers and targets for cancer therapy. Semin Cancer Biol. 2018;53:75–89.
Ji J, Wang XW. Clinical implications of cancer stem cell biology in hepatocellular carcinoma. Semin Oncol. 2012;39:461–72.
Huo X, Han S, Wu G, Latchoumanin O, Zhou G, Hebbard L, et al. Dysregulated long noncoding RNAs (lncRNAs) in hepatocellular carcinoma: implications for tumorigenesis, disease progression, and liver cancer stem cells. Mol Cancer. 2017;16:165.
Klingenberg M, Matsuda A, Diederichs S, Patel T. Non-coding RNA in hepatocellular carcinoma: mechanisms, biomarkers and therapeutic targets. J Hepatol. 2017;67:603–18.
Bharti R, Dey G, Mandal M. Cancer development, chemoresistance, epithelial to mesenchymal transition and stem cells: a snapshot of IL-6 mediated involvement. Cancer Lett. 2016;375:51–61.
Chen HA, Kuo TC, Tseng CF, Ma JT, Yang ST, Yen CJ, et al. Angiopoietin-like protein 1 antagonizes MET receptor activity to repress sorafenib resistance and cancer stemness in hepatocellular carcinoma. Hepatology. 2016;64:1637–51.
Wei X, You X, Zhang J, Zhou C. MicroRNA-1305 inhibits the Stemness of LCSCs and tumorigenesis by repressing the UBE2T-dependent Akt-signaling pathway. Mol Ther Nucleic Acids. 2019;16:721–32.
Huan HB, Yang DP, Wen XD, Chen XJ, Zhang L, Wu LL, et al. HOXB7 accelerates the malignant progression of hepatocellular carcinoma by promoting stemness and epithelial-mesenchymal transition. J Exp Clin Cancer Res. 2017;36:86.
Zhang PP, Wang PQ, Qiao CP, Zhang Q, Zhang JP, Chen F, et al. Differentiation therapy of hepatocellular carcinoma by inhibiting the activity of AKT/GSK-3beta/beta-catenin axis and TGF-beta induced EMT with sophocarpine. Cancer Lett. 2016;376:95–103.
Ma XL, Shen MN, Hu B, Wang BL, Yang WJ, Lv LH, et al. CD73 promotes hepatocellular carcinoma progression and metastasis via activating PI3K/AKT signaling by inducing Rap1-mediated membrane localization of P110beta and predicts poor prognosis. J Hematol Oncol. 2019;12:37.
Ma XL, Jiang M, Zhao Y, Wang BL, Shen MN, Zhou Y, et al. Application of serum Annexin A3 in diagnosis, outcome prediction and therapeutic response evaluation for patients with hepatocellular carcinoma. Ann Surg Oncol. 2018;25:1686–94.
Ma XL, Zhou JY, Gao XH, Tian L, Wu J, Zhang CY, et al. Application of the albumin-bilirubin grade for predicting prognosis after curative resection of patients with early-stage hepatocellular carcinoma. Clin Chim Acta. 2016;462:15–22.
Zhou SL, Yin D, Hu ZQ, Luo CB, Zhou ZJ, Xin HY, et al. A Positive Feedback Loop Between Cancer Stem-Like Cells and Tumor-Associated Neutrophils Controls Hepatocellular Carcinoma Progression. Hepatology. 2019. https://doi.org/10.1002/hep.30630 Epub ahead of print.
Ma XL, Sun YF, Wang BL, Shen MN, Zhou Y, Chen JW, et al. Sphere-forming culture enriches liver cancer stem cells and reveals Stearoyl-CoA desaturase 1 as a potential therapeutic target. BMC Cancer. 2019;19:760.
Fu PY, Hu B, Ma XL, Yang ZF, Yu MC, Sun HX, et al. New insight into BIRC3: a novel prognostic indicator and a potential therapeutic target for liver cancer. J Cell Biochem. 2019;120:6035–45.
Ho DW, Tsui YM, Sze KM, Chan LK, Cheung TT, Lee E, et al. Single-cell transcriptomics reveals the landscape of intra-tumoral heterogeneity and stemness-related subpopulations in liver cancer. Cancer Lett. 2019;459:176–85.
Wu J, Zhu P, Lu T, Du Y, Wang Y, He L, et al. The long non-coding RNA LncHDAC2 drives the self-renewal of liver cancer stem cells via activation of hedgehog signaling. J Hepatol. 2019;70:918–29.
Liu Y, Zhang JB, Qin Y, Wang W, Wei L, Teng Y, et al. PROX1 promotes hepatocellular carcinoma metastasis by way of up-regulating hypoxia-inducible factor 1alpha expression and protein stability. Hepatology. 2013;58:692–705.
Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42:D92–7.
Zhang L, Chen J, Ning D, Liu Q, Wang C, Zhang Z, et al. FBXO22 promotes the development of hepatocellular carcinoma by regulating the ubiquitination and degradation of p21. J Exp Clin Cancer Res. 2019;38:101.
Suryo RA, Savov V, Brunner A, Bolin S, Weishaupt H, Malyukova A, et al. FBW7 suppression leads to SOX9 stabilization and increased malignancy in medulloblastoma. EMBO J. 2016;35:2192–212.
Yu Q, Li W, Xie D, Zheng X, Huang T, Xue P, et al. PI3Kgamma promotes vascular smooth muscle cell phenotypic modulation and transplant arteriosclerosis via a SOX9-dependent mechanism. EBioMedicine. 2018;36:39–53.
Zhang F, Li K, Yao X, Wang H, Li W, Wu J, et al. A miR-567-PIK3AP1-PI3K/AKT-c-Myc feedback loop regulates tumour growth and chemoresistance in gastric cancer. EBioMedicine. 2019;44:311–21.
Chen S, Cao W, Yue P, Hao C, Khuri FR, Sun SY. Celecoxib promotes c-FLIP degradation through Akt-independent inhibition of GSK3. Cancer Res. 2011;71:6270–81.
Chen C, Zhao S, Karnad A, Freeman JW. The biology and role of CD44 in cancer progression: therapeutic implications. J Hematol Oncol. 2018;11:64.
Gagliardi G, Ben MK, Chaffiol A, Slembrouck-Brec A, Conart JB, Nanteau C, et al. Characterization and transplantation of CD73-positive photoreceptors isolated from human iPSC-derived retinal Organoids. Stem Cell Reports. 2018;11:665–80.
Breitbach M, Kimura K, Luis TC, Fuegemann CJ, Woll PS, Hesse M, et al. In vivo labeling by CD73 Marks multipotent stromal cells and highlights endothelial heterogeneity in the bone marrow niche. Cell Stem Cell. 2018;22:262–76.
Samanta D, Park Y, Ni X, Li H, Zahnow CA, Gabrielson E, et al. Chemotherapy induces enrichment of CD47(+)/CD73(+)/PDL1(+) immune evasive triple-negative breast cancer cells. Proc Natl Acad Sci U S A. 2018;115:E1239–48.
Lupia M, Angiolini F, Bertalot G, Freddi S, Sachsenmeier KF, Chisci E, et al. CD73 regulates Stemness and epithelial-Mesenchymal transition in ovarian Cancer-initiating cells. Stem Cell Reports. 2018;10:1412–25.
Katsuta E, Tanaka S, Mogushi K, Shimada S, Akiyama Y, Aihara A, et al. CD73 as a therapeutic target for pancreatic neuroendocrine tumor stem cells. Int J Oncol. 2016;48:657–69.
Al-Salama ZT, Syed YY, Scott LJ. Lenvatinib: a review in hepatocellular carcinoma. Drugs. 2019;79:665–74.
Kudo M, Finn RS, Qin S, Han KH, Ikeda K, Piscaglia F, et al. Lenvatinib versus sorafenib in first-line treatment of patients with unresectable hepatocellular carcinoma: a randomised phase 3 non-inferiority trial. Lancet. 2018;391:1163–73.
Jang JW, Song Y, Kim SH, Kim JS, Kim KM, Choi EK, et al. CD133 confers cancer stem-like cell properties by stabilizing EGFR-AKT signaling in hepatocellular carcinoma. Cancer Lett. 2017;389:1–10.
Xiao Y, Sun Y, Liu G, Zhao J, Gao Y, Yeh S, et al. Androgen receptor (AR)/miR-520f-3p/SOX9 signaling is involved in altering hepatocellular carcinoma (HCC) cell sensitivity to the Sorafenib therapy under hypoxia via increasing cancer stem cells phenotype. Cancer Lett. 2019;444:175–87.
Kawai T, Yasuchika K, Ishii T, Miyauchi Y, Kojima H, Yamaoka R, et al. SOX9 is a novel cancer stem cell marker surrogated by osteopontin in human hepatocellular carcinoma. Sci Rep. 2016;6:30489.
Santos JC, Carrasco-Garcia E, Garcia-Puga M, Aldaz P, Montes M, Fernandez-Reyes M, et al. SOX9 elevation acts with canonical WNT signaling to drive gastric Cancer progression. Cancer Res. 2016;76:6735–46.
Ren-Quan Lu was supported by the National Natural Science Foundation of China (81572552 and 81772774), and the project from the Shanghai Science and Technology Commission (17411963500).
Lin Guo was supported by the National Natural Science Foundation of China (81772808), the project from Shanghai Science and Technology Commission (17411951000), and the construction project for important weak discipline from Shanghai Municipal Health Commission (2015ZB0204).
Ning Ren was supported by grants from the International Science and Technology Cooperation Project of Shanghai (18410721900),
Wei-Guo Tang was supported by the project from Minhang Hospital (2018MHJ05), the project from Minhang Science and Technology Commission (2019MHZ077), and the fusion grant of Minhang Hospital-College of Pharmacy, Fudan University (RO-MY201805),
Bo Hu was supported by the National Natural Science Foundation of China (8180101167 and 31771472).
Ethics approval and consent to participate
The present study was approved by the Zhongshan Hospital Research Ethics Committee, and all individuals provided their informed consent. Humane care of animals was objected to the “Guide for the Care and Use of Laboratory Animals” criteria of National Academy of Science (National Institute of Health publication 86-23, revised 1985).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Tables; Description: Table S1 and S2.
Supplementary figure legends.
Figure S1 Description: CD73 expression was increased in HCC spheres.
Figure S2 Description: Dynamic change pattern of CD73 according to serial differentiation assays.
Figure S3 Description: CD73 positive percentages in indicated HCC cell lines.
Figure S4 Description: CD73 triggers SOX9 transcription by c-Myc and enhances Sox9 protein stability via inhibiting GSK3β activity.
Figure S5 Description: CD73 mainly depended on its enzyme activity to promote CSC traits in HCC.
Figure S6 CD73 was critical for the resistance to sorafenib or Cabozantinib in HCC.
About this article
Cite this article
Ma, XL., Hu, B., Tang, WG. et al. CD73 sustained cancer-stem-cell traits by promoting SOX9 expression and stability in hepatocellular carcinoma. J Hematol Oncol 13, 11 (2020). https://doi.org/10.1186/s13045-020-0845-z
- Hepatocellular carcinoma
- Cancer stem cells
- AKT signaling
- Lenvatinib resistance