- Open Access
Sulfarotene, a synthetic retinoid, overcomes stemness and sorafenib resistance of hepatocellular carcinoma via suppressing SOS2-RAS pathway
Journal of Experimental & Clinical Cancer Research volume 40, Article number: 280 (2021)
Recurrent hepatocellular carcinoma (HCC) shows strong resistance to sorafenib, and the tumor-repopulating cells (TRCs) with cancer stem cell-like properties are considered a driver for its high recurrent rate and drug resistance.
Suppression of TRCs may thus be an effective therapeutic strategy for treating this fatal disease. We evaluated the pharmacology and mechanism of sulfarotene, a new type of synthetic retinoid, on the cancer stem cell-like properties of HCC TRCs, and assessed its preclinical efficacy in models of HCC patient-derived xenografts (PDXs).
Sulfarotene selectively inhibited the growth of HCC TRCs in vitro and significantly deterred TRC-mediated tumor formation and lung metastasis in vivo without apparent toxicity, with an IC50 superior to that of acyclic retinoid and sorafenib, to which the recurrent HCC exhibits significant resistance at advanced stage. Sulfarotene promoted the expression and activation of RARα, which down-regulated SOS2, a key signal mediator associated with RAS activation and signal transduction involved in multiple downstream pathways. Moreover, sulfarotene selectively inhibited tumorigenesis of HCC PDXs with high expression for SOS2.
Our study identified sulfarotene as a selective inhibitor for the TRCs of HCC, which targets a novel RARα-SOS2-RAS signal nexus, shedding light on a new, promising strategy of target therapy for advanced liver cancer.
Hepatocellular carcinoma (HCC) is the main type of liver cancer and ranks fifth among the most common malignancies worldwide. The incidence of HCC is expected to continue rising to become the second leading cause of cancer-related deaths . Although several treatment strategies for HCC are currently used in clinical practice , the high recurrence rate and emerging drug resistance contribute significantly to this dire situation .
Cancer stem cells (CSCs) or tumor-initiating cells (TICs) are a small, specialized subpopulation of cells that arise from the precancerous lesion, which play a predominant role in the progression, recurrence and resistance to therapy of HCC . Instead of the complex approaches used for the identification and isolation of CSCs or TICs, e.g. by sets of stem cell surface markers , that hamper study applicability, we adapted and developed a simple mechanical approach to select cancer stem-like cells while they were cultured in three-dimensional (3D) soft fibrin gels  rather than two-dimensional (2D) rigid dishes. Such selectively survived cancer stem-like cells were shown to be highly tumorigenic, and thus were defined as tumor-repopulating cells (TRCs) [6, 7]. In a previous study, our group identified a new synthetic retinoid, named sulfarotene (WYC-209), that exerted potent, selective activity in suppressing the growth and tumor-initiating ability of TRCs derived from various types of cancer with negligible toxicity . The results suggested that the combination of sulfarotene with the efficient TRC selection for cancer stem-like cells could be used to better tackle the current HCC problems as aforementioned.
In the present study, we have found that sulfarotene selectively inhibited TRCs of HCC origins and metastatic tumor formation in multiple preclinical models, including TRCs-based xenografts and patient derived xenografts (PDXs). Of note, sulfarotene effectively suppressed tumor formation and lung metastasis of the HCC TRCs that otherwise were resistance to sorafenib and ACR. Mechanistically, we demonstrated that sulfarotene upregulated RARα in HCC TRCs, which downregulated SOS2, an important mediator of oncogenic RAS activation that is critical for multiple upstream as well as downstream signaling pathways. Thus, our study has identified sulfarotene as a potential therapeutic agent for treating the TRCs of HCC by targeting a novel RARα-SOS2-RAS signal axis, which plays critical roles in tumorigenicity of TRCs and mediation of the therapeutic effects of sulfarotene.
Materials and methods
TRCs proliferation, migration and invasion assay
TRCs cultured in 3D fibrin gels were treated with agents, after 24 h, cell viability was assessed by using the Cell Counting Kit-8 (CCK-8, Dojindo, Japan) according to the manufacturer’s protocol in which the absorbance of produced formazan dye was measured at 450 nm. The cell inhibition rate was calculated according to the formula [(OD control cells - OD treated cells) / (OD control cells - OD blanks)] × 100. The capacity of cell migration was determined by transwell migration and invasion assay as previously described .
PDX, orthotopic transplant mouse models and lung metastasis mouse models
Tumor tissues freshly isolated from patients in the operating room were dissected into small blocks of 1 mm3 under aseptic conditions. In a sterile environment, NOG mice were anesthetized and the HCC tissue blocks were implanted subcutaneously into the top right flank of the mice. About two months later, after reaching 1 cm in diameter, the subcutaneous PDX tumor nodes were removed, dissected into approximately 2 × 2 × 2 mm3 pieces, and re-transplanted into the flanks of nude mice for 30 days to permit growth as previously described [10, 11].
The mice were euthanized at day 30 or when a tumor node reached 15 mm in diameter. In metastasis experiments, the TRCs and treatment agents were injected i.v. through the tail vein into BALB/C nude mice. The lung was dissected and subjected to bioluminescence imaging in IVIS Lumina II with image radiance values normalized by the Living Image program (PerkinElmer, Boston, MA). The weight of the metastatic lung was measured and the metastasis rate was determined on paraffin embedded and HE stained lung sections by counting detectable tumor foci. PLC/PRF/5-TRCs stably overexpressing SOS2 or shRNA targeting SOS2 (see details in the Supplemental Materials) were established by puromycin antibiotic selection. The inoculation of these SOS2-manipulated TRCs into nude mice and subsequent drug treatment schemes essentially followed the methods described above.
Identification of critical responsive genes by ranking and network analyses
Critical genes responsive to sulfarotene treatment were selected and ranked based on an integrated analysis of RNA-Seq and ChIP-Seq data according to the priority criteria as follows: (1) degree of neighbor nodes. The correlation of the integrated genes obtained from RNA-Seq and ChIP-Seq data (|PCCs| > 0.9 and P < 0.05) was determined by Pearson’s correlation coefficients (PCCs) analysis, from which a molecular network was constructed with differential degrees determined by the number of neighboring genes linked to each critical gene that was considered to play an important role in the network; (2) total number of involved pathways. After selection of a top list of associated genes according to the degree value (≥ 18), critical genes were differentially ranked based on the number of the involved biological pathways, especially those closely related to tumorigenesis, e.g., the RAS signaling pathway and pathways associated with cell cycle progression, pluripotency of stem cells and ubiquitin mediated proteolysis in the KEGG annotated pathways ; (3) ChIP-Seq peak values, which indicate expression values of genes in association with RARα with statistic significance (M-value ≥ 0.5, P-value < 0.05); (4) the membership in 3 clusters. Assuming genes downregulated in response to sulfarotene treatment as having priority roles as oncogenes, genes enriched in clusters 1, 3 and 5 based on the change patterns were considered as targets for sulfarotene treatment with a membership score ≥ 0.2.
HCC patients recruitment and follow-up
A total of 282 patients with HCC who underwent curative resection between January 2009 and January 2010 were enrolled in 2 independent cohorts at the First Affiliated Hospital of Wenzhou Medical University (Zhejiang, China). All patients were identified by the pathologic diagnosis of HCC who had not yet received drug treatment. The clinicopathological data of all patients included the grade, stage and tumor location. However, among all the patients, 45 in cohort 1 had no follow-up while 237 were monitored until May 1, 2019, of whom 127 experienced postoperative pulmonary metastasis, as previously described . Each patient provided informed consent before they participated in the study. The study protocol was approved by the Ethics Committee of The First Affiliated Hospital of Wenzhou Medical University (Zhejiang, China). The overall survival (OS) was defined as the duration from surgery to death or final follow-up, while the time to recurrence (TTR) was defined as the interval between the surgery and recurrence.
Sulfarotene selectively targets tumor-repopulating cells of HCC in vitro and in vivo
First, sulfarotene was used to screen HCC cell lines, and effectively inhibited the proliferation of Hep3B and PLC/PRF/5 cells among the less invasive cells (Figure S1). Then, we established the TRC model with two human HCC cell lines, Hep3B and PLC/PRF/5, as previously described . In view of these two HCC cell lines are weaker than other HCC cell lines in terms of proliferation and invasion, the selected TRCs from Hep3B and PLC/PRF/5 cells were more stem-like compared with their parental cancer cells. FACS analyses revealed that 89.2 and 82.1 % of the resulting Hep3B-TRCs and PLC/PRF/5-TRCs, respectively, were enriched for the expression of EpCAM+ and CD133+, two surface antigens that mark the stemness of HCC stem cells  (Figure S2a, b). Transwell assays revealed significantly increased survival and migration of these TRCs over their parental cancer cells in vitro (Figure S2c). Furthermore, all inoculums of both types of TRCs, sampled at 5 × 105 cells selected from initial parental cancer cells, formed significant subcutaneous xenograft tumor nodes in BALB/c nude mice compared to none of their parental cancer cells at 30 days (Figure S2d, e), confirming the notable self-renewal and tumorigenic properties of the isolated TRCs.
Using these HCC TRC models, we evaluated the therapeutic effects of the recently discovered RA analog sulfarotene (Figure S3)  in comparison with ACR, the open-ring RA analog developed almost three decades ago, sorafenib, and the solvent carrier DMSO as the negative control. After treatment for 5 days, we found that the colony spheroids from HCC TRCs were markedly suppressed by sulfarotene at concentrations ranging from 1.0 to 10 µM, compared to 10 µM for ACR or sorafenib (Fig. 1a). Similarly, sulfarotene caused a significantly higher percent of apoptotic cell deaths in both types of HCC-related TRCs compared to the other two drugs tested (Fig. 1b and Figure S4a).
Sulfarotene exhibited the lowest IC50 values of 2.83 and 1.08 µM of the selected Hep3B-TRCs and PLC/PRF/5-TRCs after 48 h treatment, which increased to 10.20 and 6.04 µM for Hep3B and PLC/PRF/5 cancer cells (Fig. 1c and Figure S4b). However, the IC50 values of sorafenib increased 2.5 and 10 fold for these two types of TRCs, respectively, which is in accordance with the reported sorafenib resistance of CSCs and in the clinic . Analyses of cell proliferation and apoptosis biomarkers revealed that sulfarotene significantly reduced the expression of Ki67 while increasing Caspase-3 in TRCs compared with ACR, sorafenib and the carrier DMSO (Fig. 1d).
Next the effects of sulfarotene on the growth and formation of subcutaneous xenograft nodes derived from the inoculated Hep3B-TRCs in vivo in BALB/c nude mice were investigated. Both the volumes and weights of the formed Hep3B-TRCs tumor nodes were significantly reduced by sulfarotene administration every other day for a total of 25 days at a concentration of 0.22 mg/kg or 2.2 mg/kg, compared to 1.8 mg/kg for ACR and 0.1 % DMSO carrier (P < 0.001) (Fig. 1e, f). Five days after injection of 1 × 105 Hep3B-TRCs through the tail vein, treatment with 0.22 mg/kg sulfarotene for 25 days reduced the lung metastasis foci by 50 % and 2.2 mg/kg completely blocked lung metastasis. By contrast, a dose as high as 1.8 mg/kg of ACR only reduced lung metastasis by up to 25 %, and the control of 0.1 % DMSO had no effect on lung metastasis (Fig. 1 g, h). There were no detectable toxic or necrotic effects on the heart, liver, spleen, lung or kidney tissues after sulfarotene treatment, based on structural morphology of the HE stained sections (Figure S5). Taken together, these data suggest that sulfarotene exhibits significantly better selective therapeutic activity against tumorigenesis from the HCC tumor-repopulating cells than other HCC targeted drugs.
Sulfarotene modulates RARα to confer therapeutic sensitivity to HCC TRCs
RA is an important physiological ligand of the nuclear receptors RAR and RXR that when activated exert growth inhibitory effects on various tumors . However, whether sulfarotene that exerts high therapeutic selectivity and efficacy on the TRCs of HCC as described above utilizes a similar RAR- or RXR-dependent mechanism remains unclear. We first performed transcriptome sequencing and found that the mRNA levels of RARα increased under the action of sulfarotene at 5.0 µM, but not its homologs RARβ and RARγ (Fig. 2a). qPCR and western blot analyses revealed that both the mRNA and protein levels of RARα significantly increased in a dose-dependent manner in both Hep3B-TRCs and PLC/PRF/5-TRCs in response to the treatment with sulfarotene (Fig. 2b, c). Although sulfarotene treatment significantly inhibited colony formation in 3D soft fibrin gels, as well as migration and invasion in transwell culture of TRCs in vitro, pre-treatment with a selective RARα activation antagonist BMS195614 could significantly reverse these inhibitory effects (Fig. 2d, e). Similar results were obtained with RARα-specific siRNAs, which markedly reduced the mRNA and protein levels of RARα (Figure S6), accompanied by a reversal of the inhibitory effects imposed by sulfarotene on colony spheroid formation as well as the migration and invasion of TRCs (Figure S7a-c). As expected, immunofluorescence staining and nuclear fractionation studies revealed that sulfarotene treatment promoted the translocation of the activated RARα from the cytosol into the nucleus (Fig. 2f, g). These results strongly indicate that RARα is a potential target of sulfarotene in mediating the selective suppression of HCC TRCs.
SOS2 potentially serves as an oncogenic factor in HCC TRCs
RARα has been demonstrated to be a key transcriptional regulator of various oncogenes [16, 17]. To understand how RARα mediates the elevated therapeutic sensitivity of TRCs to sulfarotene (Fig. 1c), we analyzed the downstream target genes associated with RARα activation and sulfarotene effects by multiple integrative analyses. Multiple comparisons of clusters of the most differentially expressed genes (DEGs) were carried out with data obtained by RNA-Seq in Hep3B, Hep3B-TRCs and L02 cells treated with 0, 1.0 and 5.0 µM sulfarotene as well as by ChIP-Seq signals in Hep3B-TRCs treated with 5.0 µM sulfarotene compared to the control 0.1 % DMSO (Fig. 3a, b). The syngeneic DEGs analyses among the 3 groups classified 15,319 union genes of dynamic changes into Venn diagrams and 8 patterns (clusters 1 to 8) with Mfuzz (version 2.50.0)  (Fig. 3c, d).
A subsequent integrative analysis on the RNA-Seq clusters (genes from clusters 1, 3 and 5) and the ChIP-Seq signals (M-value > 0.5, P < 0.05) mapped 402 relevant peaks from which 367 RARα-associated genes were identified (Fig. 3e). Further ranking analyses on the 367 RARα-associated genes according to criteria with 4 priorities, including the degree of neighbor nodes, the total number of involved pathways, the ChIP-Seq peak values and the membership in clusters, identified a list of the top 10 transcriptional target genes of RARα, including SOS2, PTPN11, ATG2B, DCBLD2, SUPT16H, BIRC6, ARFGEF2, MYCBP2, TTBK2, and NDST1. Among them, SOS2 appeared to be one of the most relevant genes in determining the highly differential sensitivity of Hep3B-TRCs to sulfarotene in association with RARα relative to Hep3B cells and L02 cells (Fig. 3f and Supplementary Data 1).
Analyses of single-cell sequencing data of 21 HCC samples obtained from the Gene Expression Omnibus (GEO, GSE149614) (Figure S8a) allowed us to identify a total of 43 clusters that could be classified into 9 cell subsets, including B cells, endothelial cells, hepatocytes, macrophage, monocyte, neurons, NK cells, smooth muscle cells and T cells (Figure S8b). However, none of these cell subsets exhibited clustered distributions with SOS2 gene expression (Figure S8c). Consistent with the above data, we found that SOS2 protein was co-localized with tumor-specific antigen AFP in hepatoma cells but not in other cell types in human HCC tissues (Fig. 3 g).
Furthermore, analyses of the data gathered from human HCC patients by the Cancer Genome Atlas (TCGA) showed that the expression of SOS2 in tumor foci was higher than in the paired adjacent non-tumor tissues (Figure S8d). Analyses of SOS2 expression by qPCR from a cohort of 45 patients, as well as by IHC from a cohort of 237 HCC patients in which 127 had metastasis, revealed that the levels of SOS2 in both the primary and the metastatic tumor foci were higher than in the paired peritumor tissues (Figure S8e-h). Correlation analyses on clinical characteristics and multivariate parameters of 237 HCC patients indicated that SOS2 was an independent prognostic factor associated with both overall survival (OS) (HR = 1.442, P < 0.001) and time to recurrence (TTR) (HR = 1.485, P = 0.029) (Tables S1 and S2). Survival analyses indicated that patients who were high SOS2 expressors (comprehensive positive score (CPS) > 4) had markedly shorter median overall survival times and times to recurrence than those who were low SOS2 expressors (CPS ≤ 4) (Figure S8i, j) . Functionally, reduction of SOS2 expression by shRNAs significantly impaired the abilities of HCC TRCs to form colony spheroids and of PLC/PRF/5-TRCs to grow subcutaneous xenograft tumor nodes in nude mice (Fig. 3 h, i). Taken together, these results indicated that SOS2 is potentially oncogenic in HCC, and in particular in HCC TRCs.
Sulfarotene overcomes stemness of HCC via suppressing SOS2
SOS2 (Son of Sevenless Homolog 2) is an intracellular RAS/Rho guanine nucleotide exchange factor, which promotes the exchange of RAS-GDP to RAS-GTP thereby activating RAS to allow signal transduction to multiple downstream pathways . To understand the mechanistic relationship between SOS2 and the therapeutic effects of sulfarotene in association with the tumorigenic properties of HCC TRCs, we first found that the protein levels of SOS2 in TRCs were significantly higher than in their parental HCC cells and normal liver cells (Fig. 4a). Then, a significant reduction of SOS2 levels was detected in both types of TRCs in response to sulfarotene treatment in a concentration-dependent manner (Fig. 4b, c). Stable overexpression of SOS2 significantly enhanced colony spheroid formation and subsequent formation of subcutaneous xenograft tumor nodes of these TRCs. By contrast, treatment with sulfarotene at 0.22 and 2.2 mg/kg markedly repressed both effect parameters promoted by SOS2 overexpression (Fig. 4d-g). Such suppressive effects elicited by sulfarotene were associated with the loss of both SOS2 and the cell proliferation marker Ki67 and with the increase of apoptotic Caspase-3 (Fig. 4 h). These results indicated that sulfarotene potentially reduces SOS2 expression levels to elicit the observed anti-TRC and anti-tumor activity.
Oncogenic SOS2 is transcriptionally targeted by RARα in TRCs of HCC
To understand the mechanism(s) underlying the anti-tumor effect of sulfarotene through SOS2, we set about determining if RARα and SOS2 were closely associated in response to and in mediation of sulfarotene treatment as suggested by the data shown above. Treatment with sulfarotene inhibited the expression of SOS2 while the RARα antagonist BMS195614 rescued the reduction both in Hep3B-TRCs and PLC/PRF/5-TRCs (Figure S9a). Overexpression of RARα directly lowered the expression levels of SOS2 in these TRCs (Figure S9b). Moreover, IHC analyses on the previous xenograft tumors treated with sulfarotene (Fig. 1f) revealed that SOS2 was significantly suppressed while RARα was substantially elevated upon sulfarotene treatment in comparison with ACR and the negative control DMSO (Figure S10). Nevertheless, upon sulfarotene treatment, RARα and SOS2 were observed to co-localize in the nucleuses of HCC TRCs, even though the expression of SOS2 decreased (Fig. 5a), indicating an inverse association between SOS2 and RARα.
Based on the Eukaryotic Promoter Database, analysis of the SOS2 gene promoter region revealed 2 typical RARα binding elements located at -861 and − 641 bp, respectively (Fig. 5b), which were confirmed by a luciferase reporter assay showing that RARα overexpression decreased the luciferase activity driven by the wild-type SOS2 promoter. Interestingly, disruption of these two elements not only failed to rescue but even further reduced the SOS2 expression driven by RARα overexpression (Fig. 5b), indicating the existence of other RARα response transcriptional regulatory region(s) in the SOS2 gene locus. Indeed, ChIP-Seq data showed a strong increase in RARα binding due to sulfarotene treatment in Hep3B-TRCs at + 1154 to + 1264 bp of the exon 7 locus of the SOS2 gene, while sulfarotene treatment decreased the transcriptional levels and the signals of the surrounding enhancer-like elements marked by H3K27ac of the SOS2 gene (Fig. 5c, d). Subsequent ChIP-PCR analysis indicated that the RARα-binding sites were enriched in both the promoter and exon 7 regions on the SOS2 gene in both TRCs (Fig. 5e, f, Figure S11 and Table S3). 3D structure simulation using the MOE program illustrated a functional binding contact between the critical amino acids (Gly-94, Gly-97, His-112, Arg-126, Leu-132, Cys-135 and Glu-160) in the DNA-binding domain of human RARα and the RARα response element in exon 7 of the SOS2 gene (Figure S12).
Furthermore, we found that the ChIP-Seq signals for H3K27ac decreased at the SOS2 locus in sulfarotene-treated TRCs. We speculated that the sulfarotene-induced RARα might also inhibit the super enhancers (SEs) of SOS2 or other critical gene binding sites. By ranking according to increasing H3K27ac enrichment, we identified 1,028 SEs, among which 767 were associated with RARα (Fig. 5 g). Compared to the DMSO control, RARα binding was dramatically reduced at SEs by sulfarotene (Fig. 5 h). Although no SEs could be found at the SOS2 gene locus, the ChIP-Seq signals for H3K27ac were found to be decreased at the typical enhancers (TEs) in the SOS2 gene and at the SEs of HSPB1 as a control, which has been shown to promote oncogene addiction in many types of cancer including HCC, whereas RARα increased the binding to such loci in SOS2 and HSPB1 in sulfarotene-treated TRCs (Fig. 5i, j). Taken together, these results suggest that SOS2 is a direct transcriptional target of RARα that is targeted by sulfarotene.
Sulfarotene inhibits HCC TRCs by suppressing SOS2-RAS associated signaling pathways
As SOS2 serves as a signaling center in association with RAS , we analyzed the possible KEGG pathways associated with SOS2-RAS activation in these HCC TRCs. Analyses of the dynamic changes of SOS2 and its 220 neighboring genes identified by RNA-Seq in response to sulfarotene treatment, according to significant Pearson correlation coefficients (PCCs), allowed us to abstract and construct a SOS2-centered participating pathway network (Fig. 6a and Supplementary Data 2). Most of the identified pathways are involved in RTK-mediated RAS signaling and downstream pathways, including PI3K/AKT and MEK/ERK (Fig. 6b). Using a human phospho-RTK array, we showed that the phosphorylation levels of several RTKs were significantly elevated compared with those of the parental cancer cells (Fig. 6c). Among the commonly associated pathway mediators downstream of these RTKs, the levels of GTP-RAS, p-MEK1/2, p-ERK1/2 and p-AKT were significantly elevated (Figure S13a, b). Sulfarotene treatment for 48 h markedly reduced the levels of both SOS2 and GTP-RAS while significantly increasing RARα, which were companied by decreases in p-MEK1/2, p-ERK1/2 and p-AKT in a concentration-dependent manner (Fig. 6d e). Of note, overexpression of SOS2 significantly increased the activation of these SOS2-RAS associated pathways while sulfarotene treatment or SOS2 knockdown effectively reversed these changes (Figure S13c-f).
Sulfarotene overcomes sorafenib resistance via SOS2-RAS pathway suppression
As mentioned previously, both types of TRCs exhibited strong sorafenib resistance. We found that the phosphorylation levels of several known targets of sorafenib including PDGFRα, PDGFRβ and VEGFR3 were downregulated in PLC/PRF/5-TRCs, indicating a possible mechanism of sorafenib resistance, while the increased phosphorylation levels in the increased number of other RTKs might underlie sorafenib resistance in Hep3B-TRCs (Figure S14). Furthermore, a combination of sorafenib and sulfarotene effectively abolished the inability of sorafenib to suppress the active levels of the SOS2-RAS associated mediators in these HCC TRCs (Fig. 7a, b). In addition, we calculated the combination index (CI) of sulfarotene and sorafenib, and found that the two drugs have a synergistic effect (CI = 0.51 of Hep3B-TRC and CI = 0.73 of PLC/PRF/5-TRC). Further, SOS2, as the target of sulfarotene, knockdown combined with sorafenib exhibited the lower IC50 values of 7.68 µM and 14.39 µM of the selected Hep3B-TRCs and PLC/PRF/5-TRCs after 48 h treatment compared with sorafenib alone treated for HCC TRCs, which is similar to the effect of sorafenib treated for HCC cell lines (Figure S15). Consistent with the inhibition of such growth-promoting and apoptosis-suppressing pathways, the orthotopic xenograft tumors from PLC/PRF/5-TRCs inoculums succumbed to the combination regimen in vivo accompanied by the loss of sorafenib resistance, the decreases in p-ERK1/2 and p-AKT, and the increase in Caspase-3 (Fig. 7c-g). Collectively, the results suggested that sulfarotene effectively abolished the activation and oncogenicity of SOS2-RAS associated signaling pathways thus inhibiting the tumorigenicity and drug resistance of HCC TRCs.
Sulfarotene inhibits tumor progression of HCC PDXs with high SOS2 expression
To determine potential therapeutic efficacy of sulfarotene on HCC, we established PDX models in NOG mice (nonobese diabetic mice with severe combined immunodeficiency) with freshly resected specimens from HCC patients. The HCC specimens were divided into one group with high expression levels (CPS > 4) and the other with low expression levels for SOS2 based on the normalized intensity of IHC stains (Fig. 8a). We found that the volume and weight of tumor nodes derived from PDXs of high SOS2 expressors were larger than that of the low SOS2 expressors, suggesting that the tumor nodes of high SOS2 expressors retained the characteristics of the HCC TRCs (Fig. 8b, c). After treatment with sulfarotene at 0, 0.22 and 2.2 mg/kg by i.v. injection once every 2 days for 21 days, both the volume and weight of tumor nodes of high SOS2 expressors were markedly suppressed by both the low and high dose treatments compared to that of the low SOS2 expressors and the DMSO control group (Fig. 8d-f). Sulfarotene at 2.2 mg/kg imposed a as high as 90 % reduction in the volumes of tumor nodes that originated from PDXs with high SOS2 expression compared to only a 52 % reduction in those with low SOS2 expression, demonstrating a predominant role of SOS2 in determining the malignant behavior and the high sensitivity to sulfarotene of HCC TRCs. Consistent with the aforementioned changes in key SOS2-RAS associated pathway mediators, IHC staining of tumor node sections revealed that the levels of SOS2, Ki67 and SOS2-RAS associated downstream mediators (p-MEK1/2, p-ERK1/2 and p-AKT) were remarkably inhibited by sulfarotene in the high SOS2 expressors (Fig. 8 g, h). These results suggested that sulfarotene selectively and significantly inhibited human HCC PDX tumors predominantly by blocking the SOS2-RAS nexus and associated signaling pathways.
Despite significant progress being made in the development of novel therapies, HCC remains among the most recurrent, metastatic, and thus lethal malignancies worldwide . The notable self-renewal ability and high tumorigenicity of resident CSCs have been thought to be responsible for the observed high rate of recurrence, acquired drug resistance, poor prognosis, and treatment failure [6, 7, 22]. Historically, both the natural and synthetic trans-RAs have been used to treat cancers of the hemopoietic system, in particular leukemia, by targeting the nuclear receptor RAR-RXR system, however, the efficacy remains poor for solid tumors due to drug resistance, poor aqueous solubility, poor accessibility, and a short half-life . Only one RA derivative, the acyclic retinoid with a prolonged half-life, has been used to treat recurrent HCV-related solid liver cancer after curative therapy by targeting oncogenic MYCN, which transcriptionally regulates EPCAM, a biomarker of HCC CSCs , in addition to RAR-RXR [13, 25, 26]. Much effort has led to the recent development of a new class of synthetic RA-like compound, sulfarotene (also called WYC-209), which in one study exhibited high therapeutic activity against malignant melanoma with little toxicity, which is another type of drug-resistant cancer . Furthermore, a substantial portion of patients with advanced HCC also exhibits resistance to another widely used multi-kinase inhibitor sorafenib. In this endeavor, we first adapted a simple, recently developed mechanical approach to isolate CSC-like tumor-repopulating cells in 3D soft fibrin gels  from human HCC cell lines. We showed that the isolated TRCs recapitulated the tumorigenic and drug-resistant characteristics of CSCs while sulfarotene effectively and selectively suppressed these features in vitro and aborted the growth and formation of xenograft and PDX tumors and associated lung metastasis derived from TRCs inoculums in vivo, without notable side effects compared to ACR and sorafenib.
As a new class of structural analog of RA that serves as an agonist of RAR or RXR, sulfarotene indeed activates RARα, which is in agreement with a previous report . Analyses of the dynamic, differential changes in expression of genes in response to sulfarotene treatment and of the neighboring gene network functionally associated with RARα activation, revealed SOS2 as a critical player in mediating the therapeutic effects of sulfarotene. Furthermore, we found that SOS2 is a direct target of RARα as a transcription factor, with several typical RAR response elements to which RARα binds to repress expression. Like the family member SOS1, SOS2 belongs to the GEF family that promotes the exchange of RAS-GDP to RAS-GTP, thereby activating RAS, leading to activation of an array of downstream signaling pathways, notably the PI3K/AKT and MEK/ERK pathways, that promote cell growth, survival and migration . Aberrant signaling of these pathways is associated with many types of cancers [28,29,30]. SOS1 was identified as a predictive gene for HCC prognosis in association with RAS , and a small molecular inhibitor BAY-293 that disrupts SOS1-KRAS interaction blocked RAS activation, leading to potent anti-proliferative effects . Deficiency of SOS1 and SOS2 prevented malignant progression of chemical-induced skin papilloma . In the present study, we found that both human HCC tumor foci and the selected TRCs were among the high expressors for SOS2, but not for SOS1, which remained low and unchanged.
Of importance, we observed that the elevated SOS2 levels or RAS activation in the selected TRCs was in accordance with increased activation of several upstream RTKs as well as downstream MEK/ERK and PI3K/AKT pathways while loss of SOS2 after sulfarotene treatment blocked the accumulation of RAS-GTP and activation of associated pathways, highlighting SOS2-RAS as a central oncogenic nexus. Apart from the aforementioned possible mechanisms, this may also explain why HCC TRCs are resistant to sorafenib while sulfarotene or a combination of sorafenib with sulfarotene effectively overcomes the resistance. Although mutations in RAS were not frequently detected in HCC, overexpression or overactivation of RTKs upstream of RAS can be prominent . In addition to CSCs evolution resistance to sorafenib, there are also studies that have found that epithelial-mesenchymal transition (EMT) transition in tumor microenvironment led to sorafenib resistance. As the core factor of EMT, Snail, miR-182-5p/lncRNA-POIR and ZNF703 are highly expressed in sorafenib resistance cells, while Snail , miR-182-5p/lncRNA-POIR  and ZNF703  knockdown could restore the sensitivity to sorafenib. Therefore, Regorafenib reversed sorafenib resistance by inhibiting ERK and STAT3, and subsequently downregulating Snail and EMT . Moreover, dysregulated metabolism could be associated with increased sorafenib resistance in HCC . Sorafenib also promoted to export glucose absorption and lactic acid . PFKFB3 , HK2 , and PKM2 , as the key enzymes of glycolysis, have been shown to be overexpressed in sorafenib resistant HCC cell lines to increase glycolysis flux, which were silenced for contributing to synergetic effect with sorafenib. Compared to normal HCC cells, sorafenib resistant CSCs have higher rates of glucose consumption and lactate production, which are highly dependent on glycolysis [44,45,46]. In our study, the rates of glucose consumption and lactate production were found to be much higher in HCC TRCs, compared to the corresponding cell lines, which is consistent with other studies . However, the application of sulfarotene has not resulted in the decrease of glucose consumption and lactate production in HCC TRCs, suggesting that sulfarotene could not reverse sorafenib resistance through the glycolysis pathway (Figure S16). Overall, it would be interesting to reveal the plasticity of metabolism in liver CSCs and its contribution to sorafenib resistance.
Taken together, our findings have identified sulfarotene as a potentially effective agent that selectively targets HCC tumor-repopulating cells, which in many aspects resemble cancer stem cells and contribute significantly to the recurrence and drug resistance of HCC, and highlights SOS2 as a critical new oncogenic factor in association with RARα and RAS that together form a novel signal nexus in HCC.
Availability of data and materials
The obtained RNA-Seq and ChIP-Seq raw data were uploaded to the Sequence Read Archive (SRA) database of the National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/), with an accession number of PRJNA673935. All data are available in the main text or in the supplementary materials.
Cancer stem cells
Small hairpin RNA
Time to recurrence
Quantitative reverse transcription-polymerase chain reaction
Fitzmaurice C, Abate D, Abbasi N, Abbastabar H, Abd-Allah F, Abdel-Rahman O, et al. Global, Regional, and National Cancer Incidence, Mortality, Years of Life Lost, Years Lived With Disability, and Disability-Adjusted Life-Years for 29 Cancer Groups, 1990 to 2017: A Systematic Analysis for the Global Burden of Disease Study. JAMA Oncol. 2019;5(12):1749–68.
Gyawali B, Prasad V. Health policy: Me-too drugs with limited benefits - the tale of regorafenib for HCC. Nat Rev Clin Oncol. 2017;14(11):653–4.
Llovet JM, Montal R, Sia D, Finn RS. Molecular therapies and precision medicine for hepatocellular carcinoma. Nat Rev Clin Oncol. 2018;15(10):599–616.
Magee JA, Piskounova E, Morrison SJ. Cancer stem cells: impact, heterogeneity, and uncertainty. Cancer Cell. 2012;21(3):283–96.
Ma S, Chan KW, Hu L, Lee TK, Wo JY, Ng IO, et al. Identification and characterization of tumorigenic liver cancer stem/progenitor cells. Gastroenterology. 2007;132(7):2542–56.
Liu J, Tan Y, Zhang H, Zhang Y, Xu P, Chen J, et al. Soft fibrin gels promote selection and growth of tumorigenic cells. Nat Mater. 2012;11(8):734–41.
Huang W, Hu H, Zhang Q, Wu X, Wei F, Yang F, et al. Regulatory networks in mechanotransduction reveal key genes in promoting cancer cell stemness and proliferation. Oncogene. 2019;38(42):6818–34.
Chen J, Cao X, An Q, Zhang Y, Li K, Yao W, et al. Inhibition of cancer stem cell like cells by a synthetic retinoid. Nat Commun. 2018;9(1):1406.
Tang W, Lv B, Yang B, Chen Y, Yuan F, Ma L, et al. TREM2 acts as a tumor suppressor in hepatocellular carcinoma by targeting the PI3K/Akt/β-catenin pathway. Oncogenesis. 2019;8(2):9.
Nazzal M, Sur S, Steele R, Khatun M, Patra T, Phillips N, et al. Establishment of a Patient-Derived Xenograft Tumor From Hepatitis C-Associated Liver Cancer and Evaluation of Imatinib Treatment Efficacy. Hepatology. 2020;72(2):379–88.
Yang B, Li M, Tang W, Liu W, Zhang S, Chen L, et al. Dynamic network biomarker indicates pulmonary metastasis at the tipping point of hepatocellular carcinoma. Nat Commun. 2018;9(1):678.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Qin XY, Suzuki H, Honda M, Okada H, Kaneko S, Inoue I, et al. Prevention of hepatocellular carcinoma by targeting MYCN-positive liver cancer stem cells with acyclic retinoid. Proc Natl Acad Sci U S A. 2018;115(19):4969–74.
Tang W, Chen Z, Zhang W, Cheng Y, Zhang B, Wu F, et al. The mechanisms of sorafenib resistance in hepatocellular carcinoma: theoretical basis and therapeutic aspects. Signal Transduct Target Ther. 2020;5(1):87.
Wołoszynowska-Fraser MU, Kouchmeshky A, McCaffery P. Vitamin A and Retinoic Acid in Cognition and Cognitive Disease. Annu Rev Nutr. 2020;40:247–72.
Schwartz DM, Farley TK, Richoz N, Yao C, Shih HY, Petermann F, et al. Retinoic Acid Receptor Alpha Represses a Th9 Transcriptional and Epigenomic Program to Reduce Allergic Pathology. Immunity. 2019;50(1):106 – 20.e10.
Tan Y, Wang X, Song H, et al. A PML/RARα direct target atlas redefines transcriptional deregulation in acute promyelocytic leukemia. Blood. 2021;137(11):1503-16.
Zhou G, Soufan O, Ewald J, Hancock REW, Basu N, Xia J. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 2019;47(W1):W234-w41.
Otali D, Fredenburgh J, Oelschlager DK, Grizzle WE. A standard tissue as a control for histochemical and immunohistochemical staining. Biotech Histochem. 2016;91(5):309–26.
Sheffels E, Sealover NE, Wang C, Kim DH, Vazirani IA, Lee E, et al. Oncogenic RAS isoforms show a hierarchical requirement for the guanine nucleotide exchange factor SOS2 to mediate cell transformation. Sci Signal. 2018;11:546.
Llovet JM, Kelley RK, Villanueva A, Singal AG, Pikarsky E, Roayaie S, et al. Hepatocellular carcinoma. Nat Rev Dis Primers. 2021;7(1):6.
Jiang MJ, Chen YY, Dai JJ, Gu DN, Mei Z, Liu FR, et al. Dying tumor cell-derived exosomal miR-194-5p potentiates survival and repopulation of tumor repopulating cells upon radiotherapy in pancreatic cancer. Mol Cancer. 2020;19(1):68.
Masetti R, Vendemini F, Zama D, Biagi C, Gasperini P, Pession A. All-trans retinoic acid in the treatment of pediatric acute promyelocytic leukemia. Expert Rev Anticancer Ther. 2012;12(9):1191–204.
Yamashita T, Ji J, Budhu A, Forgues M, Yang W, Wang HY, et al. EpCAM-positive hepatocellular carcinoma cells are tumor-initiating cells with stem/progenitor cell features. Gastroenterology. 2009;136(3):1012–24.
Matsushima-Nishiwaki R, Okuno M, Takano Y, Kojima S, Friedman SL, Moriwaki H. Molecular mechanism for growth suppression of human hepatocellular carcinoma cells by acyclic retinoid. Carcinogenesis. 2003;24(8):1353–9.
Okita K, Izumi N, Ikeda K, Osaki Y, Numata K, Ikeda M, et al. Survey of survival among patients with hepatitis C virus-related hepatocellular carcinoma treated with peretinoin, an acyclic retinoid, after the completion of a randomized, placebo-controlled trial. J Gastroenterol. 2015;50(6):667–74.
Pierre S, Bats AS, Chevallier A, Bui LC, Ambolet-Camoit A, Garlatti M, et al. Induction of the Ras activator Son of Sevenless 1 by environmental pollutants mediates their effects on cellular proliferation. Biochem Pharmacol. 2011;81(2):304–13.
Lu S, Jang H, Gu S, Zhang J, Nussinov R. Drugging Ras GTPase: a comprehensive mechanistic and signaling structural view. Chem Soc Rev. 2016;45(18):4929–52.
Frau M, Feo F, Pascale RM. Pleiotropic effects of methionine adenosyltransferases deregulation as determinants of liver cancer progression and prognosis. J Hepatol. 2013;59(4):830–41.
Singh AM, Reynolds D, Cliff T, Ohtsuka S, Mattheyses AL, Sun Y, et al. Signaling network crosstalk in human pluripotent cells: a Smad2/3-regulated switch that controls the balance between self-renewal and differentiation. Cell Stem Cell. 2012;10(3):312–26.
Zhai X, Xue Q, Liu Q, Guo Y, Chen Z. Classifier of cross talk genes predicts the prognosis of hepatocellular carcinoma. Mol Med Rep. 2017;16(3):3253–61.
Hillig RC, Sautier B, Schroeder J, Moosmayer D, Hilpmann A, Stegmann CM, et al. Discovery of potent SOS1 inhibitors that block RAS activation via disruption of the RAS-SOS1 interaction. Proc Natl Acad Sci U S A. 2019;116(7):2551–60.
Liceras-Boillos P, Jimeno D, García-Navas R, et al. Differential Role of the RasGEFs Sos1 and Sos2 in Mouse Skin Homeostasis and Carcinogenesis. Mol Cell Biol. 2018;38(16):e00049-18.
Leung CON, Tong M, Chung KPS, Zhou L, Che N, Tang KH, et al. Overriding Adaptive Resistance to Sorafenib Through Combination Therapy With Src Homology 2 Domain-Containing Phosphatase 2 Blockade in Hepatocellular Carcinoma. Hepatology. 2020;72(1):155–68.
Zhao H, Cheng X, Yu J, Li Y. Stabilization of snail maintains the sorafenib resistance of hepatocellular carcinoma cells. Arch Biochem Biophys. 2021;699:108754.
Chen BW, Zhou Y, Wei T, Wen L, Zhang YB, Shen SC, et al. lncRNA-POIR promotes epithelial-mesenchymal transition and suppresses sorafenib sensitivity simultaneously in hepatocellular carcinoma by sponging miR-182-5p. J Cell Biochem. 2021;122(1):130–42.
Wang H, Xu H, Ma F, Zhan M, Yang X, Hua S, et al. Zinc finger protein 703 induces EMT and sorafenib resistance in hepatocellular carcinoma by transactivating CLDN4 expression. Cell Death Dis. 2020;11(4):225.
Chen W, Yang J, Zhang Y, Cai H, Chen X, Sun D. Regorafenib reverses HGF-induced sorafenib resistance by inhibiting epithelial-mesenchymal transition in hepatocellular carcinoma. FEBS Open Bio. 2019;9(2):335–47.
Xia S, Pan Y, Liang Y, Xu J, Cai X. The microenvironmental and metabolic aspects of sorafenib resistance in hepatocellular carcinoma. EBioMedicine. 2020;51:102610.
Huang Q, Li J, Xing J, Li W, Li H, Ke X, et al. CD147 promotes reprogramming of glucose metabolism and cell proliferation in HCC cells by inhibiting the p53-dependent signaling pathway. J Hepatol. 2014;61(4):859–66.
Li S, Dai W, Mo W, Li J, Feng J, Wu L, et al. By inhibiting PFKFB3, aspirin overcomes sorafenib resistance in hepatocellular carcinoma. Int J Cancer. 2017;141(12):2571–84.
Yoo JJ, Yu SJ, Na J, et al. Hexokinase-II Inhibition Synergistically Augments the Anti-tumor Efficacy of Sorafenib in Hepatocellular Carcinoma. Int J Mol Sci. 2019;20(6):1292.
Feng J, Wu L, Ji J, Chen K, Yu Q, Zhang J, et al. PKM2 is the target of proanthocyanidin B2 during the inhibition of hepatocellular carcinoma. J Exp Clin Cancer Res. 2019;38(1):204.
Zhang HL, Wang MD, Zhou X, Qin CJ, Fu GB, Tang L, et al. Blocking preferential glucose uptake sensitizes liver tumor-initiating cells to glucose restriction and sorafenib treatment. Cancer Lett. 2017;388:1–11.
Bort A, Sánchez BG, Mateos-Gómez PA, Vara-Ciruelos D, Rodríguez-Henche N, Díaz-Laviada I. Targeting AMP-activated kinase impacts hepatocellular cancer stem cells induced by long-term treatment with sorafenib. Mol Oncol. 2019;13(5):1311–31.
Wei Z, Jia J, Heng G, Xu H, Shan J, Wang G, et al. Sirtuin-1/Mitochondrial Ribosomal Protein S5 Axis Enhances the Metabolic Flexibility of Liver Cancer Stem Cells. Hepatology. 2019;70(4):1197–213.
Li Y, Luo S, Ma R, Liu J, Xu P, Zhang H, et al. Upregulation of cytosolic phosphoenolpyruvate carboxykinase is a critical metabolic event in melanoma cells that repopulate tumors. Cancer Res. 2015;75(7):1191–6.
We thank LC Biotech for aiding in bioinformatics analyses of RNA-Seq and ChIP-seq data.
This study was funded by the Natural Science Foundation of Shanghai (No. 20ZR1410400 to X. Cao), National Natural Science Foundation of China (No. 81772590 and 81572395 to JL. Xia, No. 82173662 to X. Cao, No. 11902121 and 32071306 to JW. Chen). The Open Funding of Key Laboratory of Diagnosis and Treatment of Severe Hepato-pancreatic Diseases of Zhejiang Province (No. 2018E10008 to X. Cao), and CAS Interdisciplinary Innovation Team (No. JCTD-2019-07 to X. Cao).
Ethics approval and consent to participate
The animal study protocols were performed in accordance with the Guide for the Care and Use of Laboratory Animals stipulated by the National Academy of Sciences and the National Institutes of Health (NIH publication 86 − 23, revised 1985) and approved by the Animal Care and Use Committee of Zhongshan Hospital, Fudan University, Shanghai, China. Each patient provided informed consent before they participated in the study. The study protocol was approved by the Ethics Committee of The First Affiliated Hospital of Wenzhou Medical University (Zhejiang, China, No. 123/2020).
Consent for publication
The authors declare that they agree to submit the article for publication.
All authors declare that no competing interests are associated with this study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Supplementary Information.
Additional file 2: Supplementary Data 1.
The critical genes clustered in response to sulfarotene and in association with RARα.
Additional file 3: Supplementary Data 2.
The association of SOS2 and associated neighboring genes.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Qi, F., Qin, W., Zhang, Y. et al. Sulfarotene, a synthetic retinoid, overcomes stemness and sorafenib resistance of hepatocellular carcinoma via suppressing SOS2-RAS pathway. J Exp Clin Cancer Res 40, 280 (2021). https://doi.org/10.1186/s13046-021-02085-4
- Hepatocellular carcinoma
- Tumor-repopulating cells