Hepatocellular Carcinoma Stem Research Articles

Abstract 6065: ASF1B potentiates stem cell traits and tumor progression in hepatocellular carcinoma via histone H3.3-dependent transcriptional reprogramming

Abstract Hepatocellular carcinoma (HCC) is one of the frequently diagnosed cancers worldwide and has emerged as the main cause of cancer-related death with a very poor prognosis, which is closely related to cancer stem cell (CSC), a subset of cells known to be a root of tumor recurrence, metastasis, and chemoresistance in HCC. However, the characteristics of CSCs in hepatocellular carcinoma (HCC) are still elusive. Our previous in vitro hepatocyte differentiation model revealed that the expression of anti-silencing function 1B histone chaperone (ASF1B) was remarkably upregulated in the stem cell-like stages. In HCC patients, ASF1B was highly expressed in cancer tissues compared with normal counterparts and was positively correlated with poor survival. Functional assays indicated that ASF1B could considerably promote the self-renewal, growth, and metastasis of HCC both in vitro and in vivo. Mechanistically, ASF1B could directly bind to histone H3 Variant H3.3 and enhance its incorporation into the promoter regions of certain stemness- and metastasis-associated genes including SOX9, KLF4, and TBX3, contributing to their augmented transcriptional and protein expression levels. H3.3 knock-down decreased the expression of these stemness and metastasis-associated genes, eliminating the stemness and oncogenic properties of ASF1B. Overall, our study revealed that ASF1B could facilitate HCC stemness and progression via H3.3-dependent transcriptional reprogramming. ASF1B might be a novel diagnostic and prognostic biomarker for HCC patients. Citation Format: Xiaona Fang, Beilei Liu, Shan Liu, Yu Zhang, Lanqi Gong, Qian Yan, Xinyuan Guan. ASF1B potentiates stem cell traits and tumor progression in hepatocellular carcinoma via histone H3.3-dependent transcriptional reprogramming [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 6065.

Abstract 6064: ADAR1-dependent RNA editing of GLI1 drives hepatocellular carcinoma stem cell self-renewal by initiating mitophagy

Abstract Background: Hepatocellular carcinoma (HCC) is one of the most common and most difficult to cure malignancies worldwide. Identifying new effective therapeutic targets is of great significance for the management of HCC. Cumulative evidence suggests cancer stem cells (CSCs) are key drivers of tumor growth and heterogeneity. Post-transcriptional RNA editing of adenosine to inosine (A-to-I) catalyzed by ADAR1 dramatically alters cellular transcriptome in cancers. However, ADAR1 editase-dependent mechanisms governing hepatic CSC generation have not been elucidated. Therefore, we have systematically and comprehensively investigated ADAR1’s effect on CSC self-renewal in HCC. Methods: An in vitro human hepatocyte differentiation model along with relevant transcriptomic data form GEO and TCGA were analyzed to characterize the oncofetal role of ADAR1. Using genetic approaches, RNA sequencing was performed to identify putative ADAR1-mediated recoding edited genes in HCC cell lines. Clinical implication of GLI1 editing was studied in a cohort of 88 HCC patients. GLI1 editing-mediated change in its promoter activity and protein stability was investigated by dual reporter assay, CHIP, Co-IP, and PLA assays. Functional difference in stemness properties, including abilities of self-renewal, differentiation, tumorigenesis, chemo-resistance and metastasis between wild-type and edited GLI1 (GLI1wt vs. GLI1R701G) and the exact mechanisms were also studied in cell models and mice. Results: Lentiviral ADAR1 wild-type but not editing-defective ADAR1E912A mutant, editing at nt2101 of GLI1 transcript caused arginine-to-glycine substitution at the residue 701. Importantly, increased editing of GLI1 was implicated in the pathogenesis of HCC. Upon editing, C-terminal half of GLI1 harbored a lower susceptibility to the inhibition of SUFU, thus promoting its nuclear translocation and activation. Moreover, GLI1R701G appeared more stable than GLI1wt due to reduced formation of specific K63-GLI1 substrate and β-TrCP-GLI1 complex. Edited GLI1 was found to strongly enhance targeted activation of NANOG and SOX9, resulting in accumulation of hepatic CSCs population, hepatocarcinogenesis, sorafenib-resistance and metastasis. Additionally, GLI1 editing initiated mitophagy via PINK1/Parkin-dependent pathway. Impaired mitophagy effectively antagonized the functions of GLI1R701G on hepatic CSC self-renewal and overcome chemoresistance. Conclusion: The critical advance of this study is that ADAR1 editase activity drives GLI1-dependent maintenance of the hepatic CSC population. Our discovery of a pivotal ADAR1-GLI1-PINK1 self-renewal axis provides the first mechanistic link between RNA-editing-driven malignant progression and mitochondrial homeostasis. Thus, ADAR1 represents a unique therapeutic vulnerability in liver CSCs with active mitophagy mediated by edited GLI1. Citation Format: Jie Luo, Lanqi Gong, Yuma Yang, Jiao Huang, Xiaona Fang, Baifeng Zhang, Ying Tang, Beilei Liu, Ming Liu, Lu Bai, Victor Ho-Fun Lee, Xin-Yuan Guan. ADAR1-dependent RNA editing of GLI1 drives hepatocellular carcinoma stem cell self-renewal by initiating mitophagy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 6064.

Elucidation of the mechanism, target gene, and the most active parts of Fuzheng Yiliu formula in the treatment of hepatocellular carcinoma

Fuzheng Yiliu formula (FZYL) is a traditional Chinese medicine prescription in treating hepatocellular carcinoma (HCC) that can improve immunity and inhibit tumor, but its mechanism and molecular targets remain elusive. In this study, we seek to identify the active ingredients, the oral bioavailability, drug similarity, and target genes of FZYL. Then the possible targets of FZYL were selected by differential expression investigation and weighted co-expression framework investigation based on the data from the databases of Cancer Genome Atlas database and Gene Expression Omnibus. Ten hub genes of FZYL were identified according to the degree score, and assessed in terms of immune infiltration, tumor stemness, and tumor micro-scale environment of HCC. The most important active components of FZYL are quercetin, formononetin and (+) - medicarpin. The hub genes screened by bioinformatics were BUB1, RRM2, MYC, CHEK1, MCM2, TOP2A, CCNB2, CDKN3, CCNA2, and CCNB1. The expressing state achieved by most hub genes in the Wound healing immunization infiltration subdivided type reached over that in the rest subdivided types. Moreover, we identified that all ten hub genes affect the DNAss or RNAss. Moreover, all ten hub genes also affect StromalScore or ImmuneScore. According to the mentioned results, FZYL is likely to impact tumor inhibition by affecting immune infiltration, tumor stemness, and tumor micro-scale environment of HCC.

Novel Paired Cell Lines for the Study of Lipid Metabolism and Cancer Stemness of Hepatocellular Carcinoma.

There are few well-characterized syngeneic murine models for hepatocellular carcinoma (HCC), which limits immunological studies and the development of immunotherapies for HCC. We previously established an oncogene-induced spontaneous HCC mouse model based on transposon-mediated oncogene (AKT and NRASV12) insertion into the genome of hepatocytes to induce tumorigenesis. Two tumor clones with different levels of lipid droplets (LDs) showed similar in vitro growth but distinctive in vivo phenotypes, including divergent proliferative capability and varying induction of myeloid-derived suppressor cells (MDSCs). The two clones showed distinct gene expression related to lipid metabolism, glycolysis, and cancer stemness. Endogenous fatty acid (FA) synthesis and exogenous monounsaturated fatty acid (MUFA) consumption promoted both tumor proliferation and cancer stemness, and upregulated c-Myc in the HCC cell lines. Moreover, the LDhi HCC cell line expressed a higher level of type II IL-4 receptor, which promoted tumor proliferation through binding IL-4 or IL-13. The chromosomal DNA of two tumor clones, NHRI-8-B4 (LDhi) and NHRI-1-E4 (LDlo) showed five identical AKT insertion sites in chromosomes 9, 10, 13, 16 and 18 and two NRAS integration sites in chromosomes 2 and 3. Herein, we describe two novel HCC cell lines with distinct features of lipid metabolism related to cancer stemness and differential interplay with the immune system, and present this syngeneic HCC mouse model as a practical tool for the study of cancer stemness and discovery of new therapies targeting liver cancers.

Targeting TROY-mediated P85a/AKT/TBX3 signaling attenuates tumor stemness and elevates treatment response in hepatocellular carcinoma

BackgroundPrevious in vitro hepatocyte differentiation model showed that TROY was specifically expressed in liver progenitor cells and a small proportion of hepatocellular carcinoma cells, suggesting that TROY may participate in hepatocellular carcinoma (HCC) stemness regulation. Here, we aim to investigate the role and mechanism of TROY in HCC pathogenesis.MethodBioinformatics analysis of the TCGA dataset has been used to identify the function and mechanism of TROY. Spheroid, apoptosis, and ALDH assay were performed to evaluate the stemness functions. Validation of the downstream pathway was based on Western blot, co-immunoprecipitation, and double immunofluorescence.ResultsHCC tissue microarray study found that a high frequency of TROY-positive cells was detected in 53/130 (40.8%) of HCC cases, which was significantly associated with poor prognosis and tumor metastasis. Functional studies revealed that TROY could promote self-renewal, drug resistance, tumorigenicity, and metastasis of HCC cells. Mechanism study found that TROY could interact with PI3K subunit p85α, inducing its polyubiquitylation and degradation. The degradation of p85α subsequently activate PI3K/AKT/TBX3 signaling and upregulated pluripotent genes expression including SOX2, NANOG, and OCT4, and promoted EMT in HCC cells. Interestingly, immune cell infiltration analysis found that upregulation of TROY in HCC tissues was induced by TGF-β1 secreted from CAFs. PI3K inhibitor wortmannin could effectively impair tumor stemness to sorafenib.ConclusionWe demonstrated that TROY is an HCC CSC marker and plays an important role in HCC stemness regulation. Targeting TROY-positive CSCs with PI3K inhibitor wortmannin combined with chemo- or targeted drugs might be a novel therapeutic strategy for HCC patients.Graphical abstract

Ivermectin synergizes sorafenib in hepatocellular carcinoma via targeting multiple oncogenic pathways.

Advanced hepatocellular carcinoma (HCC) results in generally poor clinical outcomes and necessitates better therapeutic strategies. Ivermectin, which is an existing anti‐parasitic drug, has been recently identified as a novel anti‐cancer drug. In line with previous efforts, this work demonstrates the translational potential of ivermectin to treat advanced HCC. We demonstrated that ivermectin at clinically relevant concentrations was active against growth and survival in multiple HCC cell lines. We showed that ivermectin had the potential to inhibit metastasis and target HCC stem cell functions. Mechanism studies correlated well with cellular phenotypes observed in ivermectin‐treated cells, and demonstrated inhibition of mTOR/STAT3 pathway, suppression of epithelial mesenchymal transition (EMT) and reduced expression of stem cell markers. We further demonstrated that ivermectin inhibited tumor formation and growth in HCC xenograft mouse model, without causing significant toxicity in the mice. Using combination index (CI), we showed that ivermectin and sorafenib were synergistic in HCC in vitro, and this was further confirmed in vivo. Our work demonstrates the potent anti‐HCC activities of ivermectin and its multiple targets on essential oncogenic pathways. Our findings provide preclinical evidence to initialize clinical trial using ivermectin and sorafenib for treating advanced HCC.

LncRNA FGD5-AS1 enhances the proliferation and stemness of hepatocellular carcinoma cells through targeting miR-223 and regulating the expression of ECT2 and FAT1.

Hepatocellular carcinoma (HCC) is common and causes many deaths worldwide. The aim of this study is to explore the mechanism by which long non-coding RNA FGD5-AS1 regulates HCC cell proliferation and stemness. Tumor and normal adjacent tissues were harvested from HCC patients. Real-time quantitative reverse transcription-PCR was applied to examine the expression of FGD5-AS1, miR-223, Epithelial cell transforming sequence 2 (ECT2) and FAT1. The protein levels of ECT2, FAT1, proliferating cell nuclear antigen (PCNA), OCT4, CD133 and CD90 were analyzed by western blot. The localization of FGD5-AS1 was examined by Fluorescence in situ hybridization. Cell proliferation was analyzed with CCK-8 and colony formation assays. Spheroid formation was used for analyzing cell stemness. Gene interaction was examined by RNA immunoprecipitation and luciferase activity assays. A subcutaneous xenograft mouse model was established to analyze HCC growth and stemness in vivo. Immunohistochemistry staining was used to analyze the expression PCNA and OCT4 in subcutaneous tumors. FGD5-AS1 was upregulated in HCC and its high expression indicated poor prognosis of patients. High expression of FGD5-AS1 enhanced HCC cell proliferation and stemness. Knockdown of FGD5-AS1 restrained tumor growth and stemness in mice. FGD5-AS1 directly sponged miR-223 and promoted the expression of ECT2 and FAT1 in HCC. Both knockdown of miR-223 and overexpression of ECT2 and FAT1 reversed FGD5-AS1 silencing-mediated suppression of HCC cell proliferation and stemness. FGD5-AS1 directly sponged miR-223 and promoted the expression of ECT2 and FAT1 in HCC, thus enhancing HCC cell proliferation and stemness. Our study identifies potential prognostic biomarkers and therapeutic targets for HCC.

Uncovering N4-Acetylcytidine-Related mRNA Modification Pattern and Landscape of Stemness and Immunity in Hepatocellular Carcinoma.

N4-acetylcytidine (ac4C) is an ancient and conserved RNA modification. Previously, ac4C mRNA modification has been reported promoting proliferation and metastasis of tumor cells. However, it remains unclear whether and how ac4C-related mRNA modification patterns influencing the prognosis of hepatocellular carcinoma (HCC) patients. Hereby, we constructed an ac4Cscore model and classified patients into two groups and investigated the potential intrinsic and extrinsic characteristics of tumor. The ac4Cscore model, including COL15A1, G6PD and TP53I3, represented ac4C-related mRNA modification patterns in HCC. According to ac4Cscore, patients were stratified to high and low groups with distinct prognosis. Patients subject to high group was related to advanced tumor stage, higher TP53 mutation rate, higher tumor stemness, more activated pathways in DNA-repair system, lower stromal score, higher immune score and higher infiltrating of T cells regulatory. While patients attributed to low group were correlated with abundance of T cells CD4 memory, less aggressive immune subtype and durable therapy benefit. We also found ac4Cscore as a novel marker to predict patients’ prognosis with anti-PD1 immunotherapy and/or mTOR inhibitor treatment. Our study for the first time showed the association between ac4C-related mRNA modification patterns and tumor intrinsic and extrinsic characteristics, thus influencing the prognosis of patients.

Hepatitis B Virus Pregenomic RNA Reflecting Viral Replication in Distal Non-tumor Tissues as a Determinant of the Stemness and Recurrence of Hepatocellular Carcinoma.

BackgroundThe existence of hepatic cancer stem cells (CSCs) contributes to chemotherapy resistance and cancer recurrence after treatment or surgery. However, very little is known about the hepatitis B virus (HBV) replication and its relationship with the stemness of hepatocellular carcinoma (HCC) in HBV-related HCC patients.MethodsWe collected tumor tissues (T), matched adjacent non-tumor tissues (NT), and distal non-tumor tissues (FNT) from 55 HCC patients for analysis.ResultsWe found HBV DNA levels were higher in T samples than NT and FNT samples, but HBV pgRNA and total RNA expressed lower in T samples. HBV pgRNA and total RNA correlate to HBV DNA among the T, NT, and FNT samples. Further evidence for HBV replication in T samples was provided by HBV S, reverse transcriptase, and X genes sequencing, showing that HBV sequences and genotypes differed between T and matched NT and FNT samples. HBV pgRNA and total RNA showed more frequent significant correlations with CSC markers in NT samples in HBsAg-positive patients. The markers CD133 and OCT4 expressed higher in FNT samples, and HBV replication marker of pgRNA levels was significantly positively correlated to these two markers only in FNT samples. The detection of pgRNA and OCT4 in FNT was correlated to the recurrence of HCC in the resection of HCC patients. Analysis of HBV receptor, sodium taurocholate co-transporting polypeptide (NTCP), showed that NTCP was correlated negatively to CSC markers in T samples, except for the CD44.ConclusionHBV replication may present in HCC with a weak transcriptomic signature. Moreover, the expression level of HBV pgRNA in distal non-tumor tissues is a sensitive marker for HBV replication and prognosis, which is associated with CSC-related markers especially with OCT4 in distal non-tumor tissues and recurrence of HCC in HBV-related HCC patients.

From the Editor’s Desk...

From the Editor’s Desk...

SP1-induced lncRNA DUBR promotes stemness and oxaliplatin resistance of hepatocellular carcinoma via E2F1-CIP2A feedback.

Oxaliplatin-based chemotherapy is widely used to treat advanced hepatocellular carcinoma (HCC), but many patients develop drug resistance that leads to tumor recurrence. Cancer stem cells (CSCs) are known to contribute to chemoresistance, the underlying mechanism, however, remains largely unknown. In this study, we discovered a specificity protein 1 (SP1)-induced long noncoding RNA--DPPA2 upstream binding RNA (DUBR) and its high expression in HCC tissues and liver CSCs. DUBR was associated with HCC progression and poor chemotherapy response. Moreover, DUBR facilitated the stemness and oxaliplatin resistance of HCC in vitro and in vivo. Mechanistically, DUBR upregulated cancerous inhibitor of protein phosphatase 2A (CIP2A) expression through E2F1-mediated transcription regulation. DUBR also exerted function by binding microRNA (miR)-520d-5p as a competing endogenous RNA to upregulate CIP2A at mRNA level. CIP2A, in turn, stabilized E2F1 protein and activated the Notch1 signaling pathway, thereby increasing the stemness feature of HCC and leading to chemoresistance. In conclusion, we identified SP1/DUBR/E2F1-CIP2A as a critical axis to activate the Notch1 signaling pathway and promote stemness and chemoresistance of HCC. Therefore, DUBR could be a potential target in HCC treatment.

Ephrin-A3/EphA2 axis regulates cellular metabolic plasticity to enhance cancer stemness in hypoxic hepatocellular carcinoma

The highly proliferative nature of hepatocellular carcinoma (HCC) frequently results in a hypoxic intratumoural microenvironment, which creates a therapeutic challenge owing to a lack of mechanistic understanding of the phenomenon. We aimed to identify critical drivers of HCC development and progression in the hypoxic microenvironment. We performed integrative analysis of multiple transcriptomic and genomic profiles specific for HCC and hypoxia and identified the Ephrin-A3/Eph receptor A2 (EphA2) axis as a clinically relevant and hypoxia-inducible signalling axis in HCC. The functional significance and mechanistic consequences of the Ephrin-A3/EphA2 axis were examined in EFNA3- and EPHA2- knockdown/overexpressing HCC cells. The potential downstream pathways were investigated by transcriptome sequencing, quantitative reverse-transcription PCR, western blotting analysis and metabolomics. EFNA3 was frequently upregulated in HCC and its overexpression was associated with more aggressive tumour behaviours. HIF-1α directly and positively regulated EFNA3 expression under hypoxia. EFNA3 functionally contributed to self-renewal, proliferation and migration in HCC cells. EphA2 was identified as a key functional downstream mediator of EFNA3. Functional characterisation of the Ephrin-A3/EphA2 forward-signalling axis demonstrated a promotion of self-renewal ability and tumour initiation. Mechanistically, the Ephrin-A3/EphA2 axis promoted the maturation of SREBP1 and expression of its transcriptional target, ACLY, was significantly associated with the expression of EFNA3 and hypoxia markers in clinical cohorts. The metabolic signature of EPHA2 and ACLY stable knockdown HCC cells demonstrated significant overlap in fatty acid, cholesterol and tricarboxylic acid cycle metabolite profiles. ACLY was confirmed to mediate the self-renewal function of the Ephrin-A3/EphA2 axis. Our findings revealed the novel role of the Ephrin-A3/EphA2 axis as a hypoxia-sensitive modulator of HCC cell metabolism and a key contributor to HCC initiation and progression. Hepatocellular carcinoma (HCC) is a fast-growing tumour; hence, areas of the tumour often have insufficient vasculature and become hypoxic. The presence of hypoxia within tumours has been shown to negatively impact on the survival of patients with tumours, including HCC. Herein, we identified the Ephrin-A3/EphA2 axis as a key functional driver of tumour initiation and progression in response to hypoxia. Additionally, we showed that SREBP1-ACLY-mediated metabolic rewiring was an important downstream effector that induced cancer stemness in response to Ephrin-A3/EphA2 forward-signalling.

CHRNA5 Contributes to Hepatocellular Carcinoma Progression by Regulating YAP Activity.

Hepatocellular carcinoma (HCC) is a major health concern worldwide. A better understanding of the mechanisms underlying the malignant phenotype is necessary for developing novel therapeutic strategies for HCC. Signaling pathways initiated by neurotransmitter receptors, such as α5-nicotinic acetylcholine receptor (CHRNA5), have been reported to be implicated in tumor progression. However, the functional mechanism of CHRNA5 in HCC remains unclear. In this study, we explored the role of CHRNA5 in HCC and found that CHRNA5 expression was increased in human HCC tissues and positively correlated with the T stage (p < 0.05) and AJCC phase (p < 0.05). The KM plotter database showed that the high expression level of CHRNA5 was strongly associated with worse survival in HCC patients. Both in vitro and in vivo assays showed that CHRNA5 regulates the proliferation ability of HCC by regulating YAP activity. In addition, CHRNA5 promotes the stemness of HCC by regulating stemness-associated genes, such as Nanog, Sox2 and OCT4. Cell migration and invasion assays demonstrated that CHRNA5 significantly enhanced the metastasis of HCC by regulating epithelial–mesenchymal transition (EMT)-associated genes. Furthermore, we found that CHRNA5 regulates the sensitivity of sorafenib in HCC. Our findings suggest that CHRNA5 plays a key role in the progression and drug resistance of HCC, and targeting CHRNA5 may be a strategy for the treatment of HCC.

The Natural Compound Dehydrocrenatidine Attenuates Nicotine-Induced Stemness and Epithelial-Mesenchymal Transition in Hepatocellular Carcinoma by Regulating a7nAChR-Jak2 Signaling Pathways.

Background Exposure to nicotine has been observed associated with tumor progression, metastasis, and therapy resistance of many cancers. Hepatocellular carcinoma (HCC) is one major cancer related to the liver and the most difficult to treat malignancies worldwide. The underlying mechanism of nicotine in the stimulation of HCC tumorigenesis is still not studied well. Methods Classically, nicotine binds to nicotinic acetylcholine receptors (nAChRs) and induces many downstream cancer-associated signaling pathways. Big data analysis is used to explore the importance of a7nAChR-Jak2 axis in the progression of hepatocellular carcinoma. Bioinformatic analysis was performed to determine gene associated with a7nAChR-Jak2 axis of HCC patients. Biological importance of a7nAChR-Jak2 axis was investigated in vitro (Hun7 and HepG2 cell lines), and athymic nude mouse models bearing HepG2-HCC cells xenografts were established in vivo. Result We found that nicotine exposure stimulated the HCC tumorigenicity by inducing the expression of one of the key nAChRs subunit that is α7nAChR as well as the expression of Janus kinase (JAK)-2. In both the in vitro and in vivo studies, the reduced overexpression of α7nAChR and increased sensitization of HCC towards treatment is observed with dehydrocrenatidine (DHCT), a novel and potent JAK family kinase inhibitor. Interestingly, DHCT treatment results in the reduction of the epithelial-mesenchymal transition process which leads to a significant reduction of clonogenicity, migratory, and invasive ability of HCC cells. Moreover, DHCT treatment also inhibits the cancer stem cell phenotype by inhibiting the tumor-sphere formation and reducing the number of ALDH1+ cells population in nicotine-stimulated HCC cells. Conclusions Taken together, the presented results indicate the positive effect of inhibition of nicotine induced overexpression of α7nAChR and JAK2, unique to HCC. Thus, these findings suggest the nicotine effect on HCC progression via α7nAChR-mediated JAK2 signaling pathways, and DHCT treatment enhances the therapeutic potential of HCC patients via overcoming/reversing the effect of nicotine in HCC patients.

Implications of Stemness Features in 1059 Hepatocellular Carcinoma Patients from Five Cohorts: Prognosis, Treatment Response, and Identification of Potential Compounds.

Simple SummaryCancer stemness has been reported to drive hepatocellular carcinoma (HCC) tumorigenesis and treatment resistance. However, comprehensive interpretations of transcriptomic stemness features in HCC patients have not been conducted in multiple cohorts. Our aim was to interpret clinical and therapeutic implications of transcriptional stemness features and explore potential compounds for HCC treatment. We found that transcriptional stemness indexes (mRNAsi) were independently associated with worse HCC prognosis. The HCC stemness risk model (HSRM) developed in this study significantly predicted prognosis and treatment response in various HCC cohorts. Analysis of two stemness subtypes suggested several liver-specific metabolic pathways, and mutations of TP53 and RB1 were associated with HCC transcriptional stemness. Moreover, we also identified potential compounds that target HCC transcriptional stemness. Our findings comprehensively characterized transcriptional stemness as a risk factor in HCC progression and treatment.Cancer stemness has been reported to drive hepatocellular carcinoma (HCC) tumorigenesis and treatment resistance. In this study, five HCC cohorts with 1059 patients were collected to calculate transcriptional stemness indexes (mRNAsi) by the one-class logistic regression machine learning algorithm. In the TCGA-LIHC cohort, we found mRNAsi was an independent prognostic factor, and 626 mRNAsi-related genes were identified by Spearman correlation analysis. The HCC stemness risk model (HSRM) was trained in the TCGA-LIHC cohort and significantly discriminated overall survival in four independent cohorts. HSRM was also significantly associated with transarterial chemoembolization treatment response and rapid tumor growth in HCC patients. Consensus clustering was conducted based on mRNAsi-related genes to divide 1059 patients into two stemness subtypes. On gene set variation analysis, samples of subtype I were found enriched with pathways such as DNA replication and cell cycle, while several liver-specific metabolic pathways were inhibited in these samples. Somatic mutation analysis revealed more frequent mutations of TP53 and RB1 in the subtype I samples. In silico analysis suggested topoisomerase, cyclin-dependent kinase, and histone deacetylase as potential targets to inhibit HCC stemness. In vitro assay showed two predicted compounds, Aminopurvalanol-a and NCH-51, effectively suppressed oncosphere formation and impaired viability of HCC cell lines, which may shed new light on HCC treatment.

Single-Cell RNA Sequencing Reveals the Role of Phosphorylation-Related Genes in Hepatocellular Carcinoma Stem Cells.

Abnormal activation of protein kinases and phosphatases is implicated in various tumorigenesis, including hepatocellular carcinoma (HCC). Advanced HCC patients are treated with systemic therapy, including tyrosine kinase inhibitors, which extend overall survival. Investigation of the underlying mechanism of protein kinase signaling will help to improve the efficacy of HCC therapy. Combining single-cell RNA sequencing data and TCGA RNA-seq data, we profiled the protein kinases, phosphatases, and other phosphorylation-related genes (PRGs) of HCC patients in this study. We found nine protein kinases and PRGs with high expression levels that were mainly detected in HCC cancer stem cells, including POLR2G, PPP2R1A, POLR2L, PRC1, ITBG1BP1, MARCKSL1, EZH2, DTYMK, and AURKA. Survival analysis with the TCGA dataset showed that these genes were associated with poor prognosis of HCC patients. Further correlation analysis showed that these genes were involved in cell cycle-related pathways that may contribute to the development of HCC. Among them, AURKA and EZH2 were identified as two hub genes by Ingenuity Pathway Analysis. Treatment with an AURKA inhibitor (alisertib) and an EZH2 inhibitor (gambogenic) inhibited HCC cell proliferation, migration, and invasion. We also found that both AURKA and EZH2 were highly expressed in TP53-mutant HCC samples. Our comprehensive analysis of PRGs contributes to illustrating the mechanisms underlying HCC progression and identifying potential therapeutic targets for future clinical trials.

Identification of Biomarkers Associated with Hepatocellular Carcinoma Stem Cell Characteristics Based on Co-Expression Network Analysis of Transcriptome Data and Stemness Index.

Mounting evidence has revealed the key role of cancer stem cells in hepatocellular carcinoma (HCC) metastasis and therapy resistance, yet the genes maintaining HCC stem cell features remain to be explored. This study aimed to identify and validate the key biomarkers associated with HCC stemness. mRNA expression-based stemness index (mRNAsi) was calculating using one-class logistic regression algorithm. RNA-sequencing data and clinical information of HCC samples were downloaded from the cancer genome atlas (TCGA) and merged with the corresponding mRNAsi. We investigated the correlation between mRNAsi and HCC clinical characteristics, including tumor grades, pathologic stages, vascular invasion, and survival outcomes. Significant genes associated HCC stemness features were screened through weighted gene co-expression network analysis (WGCNA) and were functionally annotated using enrichment analysis. Protein-protein interaction network was constructed among significant genes and the key biomarkers were finally identified based on the maximal clique centrality (MCC) method. The expression of key biomarkers and its correlation with HCC clinical outcomes were validated using oncomine and gene expression omnibus (GEO) database. mR-NAsi was significantly higher in HCC tissues and gradually increased according to tumor grades and pathologic stages. Patients with vascular invasion or poor survival exhibited higher mRNAsi. Forty-four highly-correlated significant gens were screened through WGCNA and functionally related to cell cycle, cellular senescence, p53 signaling pathway, DNA replication, and mismatch repair. Four different GEO datasets confirmed that the expression levels of these 44 genes were notably higher in HCC tissues. We finally identified 15 key biomarkers (KIF4A, TTK, CCNB1, CDC20, NCAPG, CCNB2, CDC45, UBE2C, CENPA, AURKB, RRM2, CDCA8, BIRC5, TPX2, and KIF2C) through MCC method. The expression of these biomarkers was up-regulated in multiple types of cancers and showed a gradually increasing trend with HCC tumor grades. Furthermore, high expression levels of these biomarkers were also correlated with HCC metastasis, recurrence, sorafenib resistance, and poor overall survival. We identified 15 key biomarkers associated with HCC stemness features and these genes might serve as promising therapeutic targets for HCC.

MRNAsi-related genes can effectively distinguish hepatocellular carcinoma into new molecular subtypes

BackgroundRecent studies have shown that the mRNA expression-based stemness index (mRNAsi) can accurately quantify the similarity of cancer cells to stem cells, and mRNAsi-related genes are used as biomarkers for cancer. However, mRNAsi-driven tumor heterogeneity is rarely investigated, especially whether mRNAsi can distinguish hepatocellular carcinoma (HCC) into different molecular subtypes is still largely unknown. MethodsUsing OCLR machine learning algorithm, weighted gene co-expression network analysis, consistent unsupervised clustering, survival analysis and multivariate cox regression etc. to identify biomarkers and molecular subtypes related to tumor stemness in HCC. ResultsWe firstly demonstrate that the high mRNAsi is significantly associated with the poor survival and high disease grades in HCC. Secondly, we identify 212 mRNAsi-related genes that can divide HCC into three molecular subtypes: low cancer stemness cell phenotype (CSCP-L), moderate cancer stemness cell phenotype (CSCP-M) and high cancer stemness cell phenotype (CSCP-H), especially over-activated ribosomes, spliceosomes and nucleotide metabolism lead to the worst prognosis for the CSCP-H subtype patients, while activated amino acids, fatty acids and complement systems result in the best prognosis for the CSCP-L subtype. Thirdly, we find that three CSCP subtypes have different mutation characteristics, immune microenvironment and immune checkpoint expression, which may cause the differential prognosis for three subtypes. Finally, we identify 10 robust mRNAsi-related biomarkers that can effectively predict the survival of HCC patients. ConclusionsThese novel cancer stemness-related CSCP subtypes and biomarkers in this study will be of great clinical significance for the diagnosis, prognosis and targeted therapy of HCC patients.

SENP2 Reduces Hepatocellular Carcinoma Stemness and Improves Sorafenib Sensitivity Through Inactivating the AKT/GSK3β/CTNNB1 Pathway.

BackgroundSmall ubiquitin-like modifier specific peptidase 2 (SENP2) suppresses the progression and chemoresistance of several cancers, while few studies report its role in hepatocellular carcinoma (HCC). This study aimed to evaluate the effect of SENP2 on stemness, sorafenib sensitivity, and downstream pathway in HCC, with validation of its molecular mechanisms by compensation experiment.MethodsSENP2 was regulated by plasmid transfection; meanwhile, in a compensation experiment, protein kinase B (AKT) was activated by SC79 treatment and β-catenin (CTNNB1) was overexpressed by plasmid transfection. After modification, sorafenib sensitivity was detected by cell counting kit-8 assay; stemness was evaluated by CD133+ cell proportion and sphere formation assay.ResultsSENP2 was decreased in HCC cell lines (including Hep3B, Li7, and Huh7) compared with normal human liver epithelial cell lines, which was further reduced in HCC stem cells than in normal HCC cells. Subsequently, SENP2 overexpression inhibited CD133+ cell proportion, decreased sphere formation ability, promoted sorafenib sensitivity, suppressed AKT and glycogen synthase kinase-3β (GSK3β) phosphorylation, and reduced CTNNB1 expression in Huh7 and Hep3B cells, while SENP2 knockdown showed the reverse effects. The following compensation experiment revealed that activating AKT or overexpressing CTNNB1 promoted CD133+ cell proportion and sphere formation ability but suppressed sorafenib sensitivity in Huh7 and Hep3B cells. Moreover, activating AKT or overexpressing CTNNB1 attenuated the effect of SENP2 overexpression on stemness and sorafenib sensitivity in Huh7 and Hep3B cells.ConclusionSENP2 suppresses HCC stemness and increases sorafenib sensitivity through inactivating the AKT/GSK3β/CTNNB1 signaling pathway.

Expression, purification and characterization of anti-α2δ1/NanoLuc fusion protein

The existence of cancer stem cells is regarded as the major cause for therapeutic resistance and relapse of a variety of cancer types including hepatocellular carcinoma (HCC). However, the tracing of such a subpopulation in vivo has been challenging. We have previously demonstrated that the isoform 5 of the voltage-gated calcium channel α2δ1 subunit, which can be recognized specifically by a monoclonal antibody 1B50-1, is a bona fide surface marker for HCC stem cells. Here we developed a strategy for optical imaging of α2δ1-positive cells by using a fusion protein containing the single chain variable fragment (scFv) of Mab1B50-1 and the luciferase NanoLuc which was tagged with Flag in the C-terminal. The scFv of Mab1B50-1 was fused to the N-terminal of NanoLucFlag using overlap PCR, and the recombinant fragment, which was named as 1B50-1scFv-NanoLucFlag, was subsequently cloned into a eukaryotic expression vector. The resulting construct was transfected into FreeStyle 293F cells in suspension using PEI reagent. The expression of the fusion protein was identified as a protein with molecular weight about 50 kDa by Western blotting. After purification by ANTI-FLAG® M2 affinity chromatography, 1B50-1scFv-NanoLucFlag was demonstrated to bind to α2δ1 positive cells specifically with a Kd value of (18.62±1.84) nmol/L. Furthermore, a strong luciferase activity of 1B50-1scFv-NanoLucFlag was detected in α2δ1 positive cells following incubation with the fusion protein, indicating that the presence of α2δ1 could be quantified using this fusion protein. Hence, 1B50-1scFv-NanoLucFlag provides a potential tool for optical imaging of α2δ1 positive cancer stem cells both in vitro and in vivo.