miR-10b-5p promotes tumor growth by regulating cell metabolism in liver cancer via targeting SLC38A2

Metabolic reprogramming plays a critical role in hepatocarcinogenesis. However, the mechanisms regulating metabolic reprogramming in primary liver cancer (PLC) are unknown. Differentially expressed miRNAs between PLC and normal tissues were identified using bioinformatic analysis. RT-qPCR was used to determine miR-10b-5p and SCL38A2 expression levels. IHC, WB, and TUNEL assays were used to assess the proliferation and apoptosis of the tissues. The proliferation, migration, invasion, and apoptosis of PLC cells were determined using the CCK-8 assay, Transwell assay, and flow cytometry. The interaction between miR-10b-5p and SLC38A2 was determined using dual-luciferase reporter assay. A PLC xenograft model in BALB/c nude mice was established, and tumorigenicity and SLC38A2 expression were estimated. Finally, liquid chromatography - mass spectrometry (LC-MS) untargeted metabolomics was used to analyze the metabolic profiles of xenograft PLC tissues in nude mice. miR-10b-5p was a key molecule in the regulation of PLC. Compared with para-carcinoma tissues, miR-10b-5p expression was increased in tumor tissues. miR-10b-5p facilitated proliferation, migration, and invasion of PLC cells. Mechanistically, miR-10b-5p targeted SLC38A2 to promote PLC tumor growth. Additionally, miR-10b-5p altered the metabolic features of PLC in vivo. Overexpression of miR-10b-5p resulted in remarkably higher amounts of lumichrome, folic acid, octanoylcarnitine, and Beta-Nicotinamide adenine dinucleotide, but lower levels of 2-methylpropanal, glycyl-leucine, and 2-hydroxycaproic acid. miR-10b-5p facilitates the metabolic reprogramming of PLC by targeting SLC38A2, which ultimately boosts the proliferation, migration, and invasion of PLC cells. Therefore, miR-10b-5p and SLC38A2 are potential targets for PLC diagnosis and treatment.

Plasma methylated GNB4 and Riplet as a novel dual-marker panel for the detection of hepatocellular carcinoma

Early detection of hepatocellular carcinoma (HCC) can greatly improve the survival rate of patients. We aimed to develop a novel marker panel based on cell-free DNA (cfDNA) methylation for the detection of HCC. The differentially methylated CpG sites (DMCs) specific for HCC blood diagnosis were selected from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, then validated by the whole genome bisulphite sequencing (WGBS) of 12 paired HCC and paracancerous tissues. The clinical performance of the panel was evaluated using tissue samples [32 HCC, chronic liver disease (CLD), and healthy individuals] and plasma cohorts (173 HCC, 199 CLD, and 98 healthy individuals). The combination of G protein subunit beta 4 (GNB4) and Riplet had the optimal area under the curve (AUC) in seven candidates through TCGA, GEO, and WGBS analyses. In tissue validation, the GNB4 and Riplet showed an AUC of 100% with a sensitivity and specificity of 100% for detecting any-stage HCC. In plasma, it demonstrated a high sensitivity of 84.39% at 91.92% specificity, with an AUC of 92.51% for detecting any-stage HCC. The dual-marker panel had a higher sensitivity of 78.26% for stage I HCC than alpha-fetoprotein (AFP) of 47.83%, and a high sensitivity of 70.27% for detecting a single tumour (size ≤3 cm). In conclusion, we developed a novel dual-marker panel that demonstrates high accuracy in detecting HCC, surpassing the performance of AFP testing.

Establishment and characterization of a highly metastatic hepatocellular carcinoma cell line

The prevalence of alcohol-related hepatocellular carcinoma (HCC) has been increasing during the last decade. Cancer research requires cell lines suitable for both in vitro and in vivo assays. However, there is a lack of cell lines with a high in vivo metastatic capacity for this HCC subtype. Herein, a new HCC cell line was established, named HCC-ZJ, using cells from a patient diagnosed with alcohol-related HCC. The karyotype of HCC-ZJ was 46, XY, del (p11.2). Whole-exome sequencing identified several genetic variations in HCC-Z that occur frequently in alcohol-associated HCC, such as mutations in TERT, CTNNB1, ARID1A, CDKN2A, SMARCA2, and HGF. Cell counting kit-8 assays, colony formation assays, and Transwell assays were performed to evaluate the proliferation, migration, and sensitivity to sorafenib and lenvatinib of HCC-Z in vitro. HCC-ZJ showed a robust proliferation rate, a weak foci-forming ability, a strong migration capacity, and a moderate invasion tendency in vitro. Finally, the tumorigenicity and metastatic capacity of HCC-Z were evaluated using a subcutaneous xenograft model, an orthotopic xenograft model, and a tail-veil injection model. HCCZJ exhibited strong tumorigenicity in the subcutaneous xenograft and orthotopic tumor models. Moreover, HCC-ZJ spontaneously formed pulmonary metastases in the orthotopic tumor model. In summary, a new HCC cell line derived from a patient with alcohol-related HCC was established, which showed a high metastatic capacity and could be applied for in vitro and in vivo experiments during pre-clinical research.Highlights• An alcohol-related HCC cell line, HCC-ZJ, was established• HCC-ZJ was applicable for in vitro functional experiment and gene editing• HCC-ZJ was applicable for in vivo tumor growth and spontaneous metastasis models.

A prognostic four-gene signature and a therapeutic strategy for hepatocellular carcinoma: Construction and analysis of a circRNA-mediated competing endogenous RNA network

BACKGROUND

Hepatocellular carcinoma (HCC) has a poor long-term prognosis. The competition of circular RNAs (circRNAs) with endogenous RNA is a novel tool for predicting HCC prognosis. Based on the alterations of circRNA regulatory networks, the analysis of gene modules related to HCC is feasible.

METHODS

Multiple expression datasets and RNA element targeting prediction tools were used to construct a circRNA-microRNA-mRNA network in HCC. Gene function, pathway, and protein interaction analyses were performed for the differentially expressed genes (DEGs) in this regulatory network. In the protein-protein interaction network, hub genes were identified and subjected to regression analysis, producing an optimized four-gene signature for prognostic risk stratification in HCC patients. Anti-HCC drugs were excavated by assessing the DEGs between the low- and high-risk groups. A circRNA-microRNA-hub gene subnetwork was constructed, in which three hallmark genes, KIF4A, CCNA2, and PBK, were subjected to functional enrichment analysis.

RESULTS

A four-gene signature (KIF4A, CCNA2, PBK, and ZWINT) that effectively estimated the overall survival and aided in prognostic risk assessment in the The Cancer Genome Atlas (TCGA) cohort and International Cancer Genome Consortium (ICGC) cohort was developed. CDK inhibitors, PI3K inhibitors, HDAC inhibitors, and EGFR inhibitors were predicted as four potential mechanisms of drug action (MOA) in high-risk HCC patients. Subsequent analysis has revealed that PBK, CCNA2, and KIF4A play a crucial role in regulating the tumor microenvironment by promoting immune cell invasion, regulating microsatellite instability (MSI), and exerting an impact on HCC progression.

CONCLUSIONS

The present study highlights the role of the circRNA-related regulatory network, identifies a four-gene prognostic signature and biomarkers, and further identifies novel therapy for HCC.

Comprehensive meta-analysis reveals distinct gene expression signatures of MASLD progression

Metabolic dysfunction-associated steatotic liver disease (MASLD) and its progressive form, metabolic dysfunction-associated steatohepatitis (MASH), pose significant risks of severe fibrosis, cirrhosis, and hepatocellular carcinoma. Despite their widespread prevalence, the molecular mechanisms underlying the development and progression of these common chronic hepatic conditions are not fully understood. Here, we conducted the most extensive meta-analysis of hepatic gene expression datasets from liver biopsy samples to date, integrating 10 RNA-sequencing and microarray datasets (1,058 samples). Using a random-effects meta-analysis model, we compared over 12,000 shared genes across datasets. We identified 685 genes differentially expressed in MASLD versus normal liver, 1,870 in MASH versus normal liver, and 3,284 in MASLD versus MASH. Integrating these results with genome-wide association studies and coexpression networks, we identified two functionally relevant, validated coexpression modules mainly driven by SMOC2, ITGBL1, LOXL1, MGP, SOD3, and TAT, HGD, SLC25A15, respectively, the latter not previously associated with MASLD and MASH. Our findings provide a comprehensive and robust analysis of hepatic gene expression alterations associated with MASLD and MASH and identify novel key drivers of MASLD progression.

Trim21 deficiency in mice increases HCC carcinogenesis in a NASH context and is associated with immune checkpoint upregulation

The global pandemic of metabolic diseases has increased the incidence of hepatocellular carcinoma (HCC) in the context of non-alcoholic steatohepatitis (NASH). The downregulation of the E3 ubiquitin ligase TRIM21 has been linked to poor prognosis in different cancers including HCC. In order to investigate the role of TRIM21 in liver cancer progression on NASH, Trim21+/+ and Trim21-/- male mice were injected with streptozotocin at the neonatal stage. The hypoinsulinemic mice were then fed with a high-fat high-cholesterol diet (HFHCD) for 4, 8 or 12 weeks. All mice developed NASH which systematically resulted in HCC progression. Interestingly, compared to the Trim21+/+ control mice, liver damage was worsened in Trim21-/- mice, with more HCC nodules found after 12 weeks on HFHCD. Immune population analysis in the spleen and liver revealed a higher proportion of CD4+PD-1+ and CD8+PD-1+ T cells in Trim21-/- mice. The liver and HCC tumors of Trim21-/- mice also exhibited an increase in the number of PD-L1+ and CD68+ PD-L1+ cells. Thus, TRIM21 limits the emergence of HCC nodules in mice with NASH by potentially restricting the expression of PD-1 in lymphocytes and PD-L1 in tumors.

WDR4 promotes HCC pathogenesis through N7-methylguanosine by regulating and interacting with METTL1

BACKGROUND

The N7-methylguanosine (m7G), a modification at defined internal positions within tRNAs and rRNAs, is correlated with tumor progression. Methyltransferase like 1 (METTL1)/ WD repeat domain 4 (WDR4) mediated tRNA m7G modification, which could alter many oncogenic mRNAs translation to promote progress of multiple cancer types. However, whether and how the internal mRNA m7G modification is involved in tumorigenesis remains unclear.

METHODS

The immunohistochemistry assay was conducted to detect the expression of WDR4 and METTL1 in hepatocellular carcinoma (HCC) and the expression of both genes whether contributes to the prognosis of the survival rate of HCC patients. Then, CCK8, colony formation assays and tumor xenograft models were conducted to determine the effects of WDR4 on HCC cells in vitro and vivo. Besides, dot blot assay, m7G-MeRIP-seq and RNA-seq analysis were conducted to determine whether WDR4 contributes to m7G modification and underlying mechanism in HCC cells. Finally, rescue and CO-IP assay were conducted to explore whether WDR4 and METTL1 proteins form a complex in Huh7 cells.

RESULTS

WDR4 modulates m7G modification at the internal sites of tumor-promoting mRNAs by forming the WDR4-METTL1 complex. WDR4 knockdown downregulated the expression of mRNA and protein levels of METTL1 gene and thus further modulate the formation of WDR4-METTL1 complex indirectly. METTL1 expression was markedly correlated with WDR4 expression in HCC tissues. HCC patients with high expression of both genes had a poor prognosis.

CONCLUSIONS

WDR4 may contribute to HCC pathogenesis by interacting with and regulating the expression of METTL1 to synergistically modulate the m7G modification of target mRNAs in tumor cells.

Research progress on oncoprotein hepatitis B X‑interacting protein (Review)

Hepatitis B X‑interacting protein (HBXIP) is a membrane protein located on the lysosomal surface and encoded by the Lamtor gene. It is expressed by a wide range of tumor types, including breast cancer, esophageal squamous cell carcinoma and hepatocellular carcinoma, and its expression is associated with certain clinicopathological characteristics. In the past decade, research on the oncogenic mechanisms of HBXIP has increased and the function of HBXIP in normal cells has been gradually elucidated. In the present review, the following was discussed: The normal physiological role of the HBXIP carcinogenic mechanism; the clinical significance of high levels of HBXIP expression in different tumors; HBXIP regulation of transcription, post‑transcription and post‑translation processes in tumors; the role of HBXIP in improving the antioxidant capacity of tumor cells; the inhibition of ferroptosis of tumor cells and regulating the metabolic reprogramming of tumor cells; and the role of HBXIP in promoting the malignant progression of tumors. In conclusion, the present review summarized the existing knowledge of HBXIP, established its carcinogenic mechanism and discussed future related research on HBXIP.

Anoikis-related lncRNA signature predicts prognosis and is associated with immune infiltration in hepatocellular carcinoma

Anoikis is a programmed cell death process triggered when cells are dislodged from the extracellular matrix. Numerous long noncoding RNAs (lncRNAs) have been identified as significant factors associated with anoikis resistance in various tumor types, including glioma, breast cancer, and bladder cancer. However, the relationship between lncRNAs and the prognosis of hepatocellular carcinoma (HCC) has received limited research attention. Further research is needed to investigate this potential link and understand the role of lncRNAs in the progression of HCC. We developed a prognostic signature based on the differential expression of lncRNAs implicated in anoikis in HCC. A co-expression network of anoikis-related mRNAs and lncRNAs was established using data obtained from The Cancer Genome Atlas (TCGA) for HCC. Cox regression analyses were conducted to formulate an anoikis-related lncRNA signature (ARlncSig) in a training cohort, which was subsequently validated in both a testing cohort and a combined dataset comprising the two cohorts. Receiver operating characteristic curves, nomograms, and decision curve analyses based on the ARlncSig score and clinical characteristics demonstrated robust predictive ability. Moreover, gene set enrichment analysis revealed significant enrichment of several immune processes in the high-risk group compared to the low-risk group. Furthermore, significant differences were observed in immune cell subpopulations, expression of immune checkpoint genes, and response to chemotherapy and immunotherapy between the high- and low-risk groups. Lastly, we validated the expression levels of the five lncRNAs included in the signature using quantitative real-time PCR. In conclusion, our ARlncSig model holds substantial predictive value regarding the prognosis of HCC patients and has the potential to provide clinical guidance for individualized immunotherapy. In this study, we obtained 36 genes associated with anoikis from the Gene Ontology and Gene Set Enrichment Analysis databases. We also identified 22 differentially expressed lncRNAs that were correlated with these genes using data from TCGA. Using Cox regression analyses, we developed an ARlncSig in a training cohort, which was then validated in both a testing cohort and a combined cohort comprising data from both cohorts. Additionally, we collected eight pairs of liver cancer tissues and adjacent tissues from the Affiliated Tumor Hospital of Nantong University for further analysis. The aim of this study was to investigate the potential of ARlncSig as a biomarker for liver cancer prognosis. The study developed a risk stratification system called ARlncSig, which uses five lncRNAs to categorize liver cancer patients into low- and high-risk groups. Patients in the high-risk group exhibited significantly lower overall survival rates compared to those in the low-risk group. The model's predictive performance was supported by various analyses including the receiver operating characteristic curve, nomogram calibration, clinical correlation analysis, and clinical decision curve. Additionally, differential analysis of immune function, immune checkpoint, response to chemotherapy, and immune cell subpopulations revealed significant differences between the high- and low-risk groups. Finally, quantitative real-time PCR validated the expression levels of the five lncRNAs. In conclusion, the ARlncSig model demonstrates critical predictive value in the prognosis of HCC patients and may provide clinical guidance for personalized immunotherapy.

Nur77 inhibition of β-catenin expression mediates Hepatoblastoma progression and enhances cisplatin's therapeutic effect

Hepatoblastoma (HB) is the most common malignant tumor in children under 5 years old, but its pathogenesis remains unclear. Nur77 has been reported to be an important regulator for cancer progression in various cancer types. This study found that Nur77 was downregulated in HB tumors, compared with paracancer tissue. Knockout or overexpression of Nur77 in HB tumor cell line HepG2 and HuH6 could significantly enhance or inhibit the proliferation, migration and invasion of tumor cells both in vitro and in vivo. Further studies illustrated that Nur77 regulated the proliferation of tumor cells by affecting the expression of β-catenin. Nur77 agonist Csn-B effectively enhanced the therapeutic effect of cisplatin on HB tumors both in vitro and in vivo. This study confirms that Nur77 may act as an oncogene in HB tumors and mediate the progression of HB by inhibiting the expression of β-catenin, which provides a new targeted therapy for the clinical treatment of HB patients; meanwhile, the combination of Nur77 agonist and cisplatin treatment may improve the chemotherapeutic efficacy of HB patients, which provides a new idea for the improvement of the clinical prognosis of HB patients.

Systemic loss of CD36 aggravates NAFLD-related HCC through MEK1/2-ERK1/2 signaling pathway

BACKGROUND & AIMS

CD36, a membrane protein widely present in various tissues, is crucial role in regulating energy metabolism. The rise of HCC as a notable outcome of NAFLD is becoming more apparent. Patients with hereditary CD36 deficiency are at increased risk of NAFLD. However, the impact of CD36 deficiency on NAFLD-HCC remains unclear.

METHODS

Global CD36 knockout mice (CD36KO) and wild type mice (WT) were induced to establish NAFLD-HCC model by N-nitrosodiethylamine (DEN) plus high fat diet (HFD). Transcriptomics was employed to examine genes that were expressed differentially.

RESULTS

Compared to WT mice, CD36KO mice showed more severe HFD-induced liver issues and increased tumor malignancy. The MEK1/2-ERK1/2 pathway activation was detected in the liver tissues of CD36KO mice using RNA sequencing and Western blot analysis.

CONCLUSION

Systemic loss of CD36 leaded to the advancement of NAFLD to HCC by causing lipid disorders and metabolic inflammation, a process that involves the activation of MAPK signaling pathway. We found that CD36 contributes significantly to the maintenance of metabolic homeostasis in NAFLD-HCC.

Blocking of FGFR4 signaling by F30 inhibits hepatocellular carcinoma cell proliferation through HMOX1-dependent ferroptosis pathway

Excessive activation of FGF19/fibroblast growth factor receptor 4 (FGFR4) signaling is associated with poor survival of patients with hepatocellular carcinoma (HCC). FGFR4 inhibitors show promise for HCC treatment. F30, an indazole derivative designed through computer-aided drug design targeting FGFR4, demonstrated anti-HCC activity as described in our previous studies. However, the precise molecular mechanisms underlying F30's anticancer effects remain largely unexplored. We report here that F30 could effectively induce ferroptosis in HCC cells. The concentrations of cellular ferrous iron, the peroxidation of cell membranes and the homeostasis of reduced glutathione (GSH)/oxidized glutathione disulfide (GSSG) were dysregulated by F30, thereby affecting cellular redox status. Induction of ferroptosis in HCC by F30 was inhibited by specific ferroptosis inhibitor ferrostatin-1. F30 upregulates various ferroptosis-related genes, including the heme oxygenase enzymes 1 (HMOX1), a key mediator of redox regulation. Surprisingly, F30-induced ferroptosis in HCC is dependent on HMOX1. The dysregulation of cellular ferrous iron concentrations and cell membrane peroxidation was rescued when knocking down HMOX1 with specific small interfering RNA. These findings shed light on the molecular mechanisms underlying FGFR4-targeting F30's anti-HCC effects and suggest that FGFR4 inactivation could be beneficial for HCC treatment involving ferroptosis.

Combining a glycolysis‑related prognostic model based on scRNA‑Seq with experimental verification identifies ZFP41 as a potential prognostic biomarker for HCC

Hepatocellular carcinoma (HCC) is a common malignancy with a poor prognosis, and its heterogeneity affects the response to clinical treatments. Glycolysis is highly associated with HCC therapy and prognosis. The present study aimed to identify a novel biomarker for HCC by exploring the heterogeneity of glycolysis in HCC. The intersection of both marker genes of glycolysis‑related cell clusters from single‑cell RNA sequencing analysis and mRNA data of liver HCC from The Cancer Genome Atlas were used to construct a prognostic model through Cox proportional hazard regression and the least absolute shrinkage and selection operator Cox regression. Data from the International Cancer Genome Consortium were used to validate the results of the analysis. Immune status analysis was then conducted. A significant gene in the prognostic model was identified as a potential biomarker and was verified through in vitro experiments. The results revealed that the glycolysis‑related prognostic model divided patients with HCC into high‑ and low‑risk groups. A nomogram combining the model and clinical features exhibited accurate predictive ability, with an area under the curve of 0.763 at 3 years. The high‑risk group exhibited a higher expression of checkpoint genes and lower tumor immune dysfunction and exclusion scores, suggesting that this group may be more likely to benefit from immunotherapy. The tumor tissues had a higher zinc finger protein (ZFP)41 mRNA and protein expression compared with the adjacent tissues. In vitro analyses revealed that ZFP41 played a crucial role in cell viability, proliferation, migration, invasion and glycolysis. On the whole, the present study demonstrates that the glycolysis‑related prognostic gene, ZFP41, is a potential prognostic biomarker and therapeutic target, and may play a crucial role in glycolysis and malignancy in HCC.

Featured lncRNA-based signature for discriminating prognosis and progression of hepatocellular carcinoma

Long non-coding RNAs (lncRNAs) have been implicated in carcinogenesis and progression of hepatocellular carcinoma (HCC). This study aimed to identify a robust lncRNA signature for predicting the survival of HCC patients. We performed an integrated analysis of the lncRNA expression profiling in The Cancer Genome Atlas (TCGA)-liver hepatocellular carcinoma database to identify the prognosis-related lncRNA for the HCC. The HCC cohort was randomly divided into a training set (n = 250) and a testing set (n = 113). Following a two-step screening, we identified an 18-lncRNA signature risk score. The high-risk subgroups had significantly shorter survival time than the low-risk group in both the training set (P < 0.0001) and the testing set (P = 0.005). Stratification analysis revealed that the prognostic value of the lncRNA-based signature was independent of the tumor stage and pathologic stage. The area under the receiver operating characteristic curve (AUROC) of the 18-lncRNA signature risk score was 0.826 (95%CI, 0.764-0.888), 0.817 (95%CI, 0.759-0.876), and 0.799 (95%CI, 0.731-0.867) for 1-year, 3-year, and 5-year follow-up, respectively. Bioinformatics analyses indicated that the 18 lncRNA might mediate cell cycle, DNA replication processes, and canonical cancer-related pathways, in which MCM3AP-AS1 was a potential target for HCC. In conclusion, the 18-lncRNA signature was a robust predictive biomarker for the prognosis and progression of HCC.

Crosstalk between m6A modification and non-coding RNAs in HCC

Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-related deaths worldwide, with high morbidity and occurrence. Although various therapeutic approaches have been rapidly developed in recent years, the underlying molecular mechanisms in the pathogenesis of HCC remain enigmatic. The N6-methyladenosine (m6A) RNA modification is believed to regulate RNA metabolism and further gene expression. This process is intricately regulated by multiple regulators, such as methylases and demethylases. Non-coding RNAs (ncRNAs) are involved in the regulation of the epigenetic modification, mRNA transcription and other biological processes, exhibiting crucial roles in tumor occurrence and development. The m6A-ncRNA interaction has been implicated in the malignant phenotypes of HCC and plays an important role in drug resistance. This review summarizes the effect of m6A-ncRNA crosstalk on HCC progression and their clinical implications as prognostic markers and therapeutic targets in this disease.

Albumin promoter-driven FlpO expression induces efficient genetic recombination in mouse liver

Tissue-specific gene manipulations are widely used in genetically engineered mouse models. A single recombinase system, such as the one using Alb-Cre, has been commonly used for liver-specific genetic manipulations. However, most diseases are complex, involving multiple genetic changes and various cell types. A dual recombinase system is required for conditionally modifying different genes sequentially in the same cell or inducing genetic changes in different cell types within the same organism. A FlpO cDNA was inserted between the last exon and 3'-UTR of the mouse albumin gene in a bacterial artificial chromosome (BAC-Alb-FlpO). The founders were crossed with various reporter mice to examine the efficiency of recombination. Liver cancer tumorigenesis was investigated by crossing the FlpO mice with FSF-KrasG12D mice and p53frt mice (KPF mice). BAC-Alb-FlpO mice exhibited highly efficient recombination capability in both hepatocytes and intrahepatic cholangiocytes. No recombination was observed in the duodenum and pancreatic cells. BAC-Alb-FlpO-mediated liver-specific expression of mutant KrasG12D and conditional deletion of p53 gene caused the development of liver cancer. Remarkably, liver cancer in these KPF mice manifested a distinctive mixed hepatocellular carcinoma and cholangiocarcinoma phenotype. A highly efficient and liver-specific BAC-Alb-FlpO mouse model was developed. In combination with other Cre lines, different genes can be manipulated sequentially in the same cell, or distinct genetic changes can be induced in different cell types of the same organism.NEW & NOTEWORTHY A liver-specific Alb-FlpO mouse line was generated. By coupling it with other existing CreERT or Cre lines, the dual recombinase approach can enable sequential gene modifications within the same cell or across various cell types in an organism for liver research through temporal and spatial gene manipulations.

Deciphering the immunological and prognostic features of hepatocellular carcinoma through ADP-ribosylation-related genes analysis and identify potential therapeutic target ARFIP2

BACKGROUND

Hepatocellular carcinoma is one of the most common malignancies, and its prognosis and treatment outcome cannot be accurately predicted. ADP-ribosylation (ADPR) is a post-translationa modification of proteins involved in protein trafficking and immune response. Therefore, it is necessary to explore the ADPR-related genes associated with the prognosis and therapeutic efficacy of hepatocellular carcinoma treatments.

METHODS

We downloaded the data of hepatocellular carcinoma samples to identify ADPR-related genes as prognostic markers, and established a novel ADPR-related index (ADPRI) based on univariate and multivariate COX regression analyses. Patients' prognosis, clinical features, somatic variant, tumor immune microenvironment, chemotherapeutic response and immunotherapeutic response were systematically analyzed. Finally, the role of ARFIP2 in hepatocellular carcinoma cells was preliminarily explored in vitro.

RESULTS

The ADPRI consisting of four ADPR related genes (ARL8B, ARFIP2, PARP12, ADPRHL1) was established to be a reliable predictor of survival in patients with hepatocellular carcinoma and was validated using external datasets. Compared with the low ADPRI group, the high ADPRI group presented higher levels of mutation frequency, immune infiltration and patients in high ADPRI group benefit more from immune checkpoint inhibitor treatment. In addition, we predicted some natural small molecule drugs as potential therapeutic targets for hepatocellular carcinoma. Finally, Knockdown of ARFIP2 inhibits the proliferation and migration of hepatocellular carcinoma cells by inducing the G1/S phase cell cycle arrest in HCC cells.

CONCLUSIONS

The ADPRI can be used to accurately predict the prognosis and immunotherapeutic response of hepatocellular carcinoma patients and providing valuable insights for future precision treatment of patients with hepatocellular carcinoma.

Detection of Naturally occurring abortive transcripts by Base-Stacking Hybridization Assisted Ligation and PCR amplification

Abortive transcripts (ATs) refer to nascent 2-10 nucleotides (nt) RNAs released by RNA polymerases before synthesizing productive RNAs. The quantitative detection of ATs is important for studying transcription initiation and the biological function of ATs; however, no method is available for the qualitative and quantitative assessment of such ultra-short oligonucleotides (typically shorter than 11 nt) in vivo at present, even with the LNA probes, the detection limit can only reach 11 nt. Here, we demonstrated the base stacking hybridization assisted ligation (BSHAL) technique, combined with TaqMan-MGB qPCR, can detect 4-10 nt ATs with a specificity of nucleotide resolution and a sensitivity of approximately 10 pM. By this technique, we detected endogenous ATs in cell lines, mice plasmas, and mice liver tissues, respectively, and proved that naturally occurring ATs do exist. We found that the 8 nt ATs of HMSB and Gapdh could be used as reference ATs for data normalization in Homo and mouse respectively, and 8 nt ATs of Afp and Gpc3 were suitable for use as plasma biomarkers of Hepatocellular carcinoma in mouse, indicate ATs are promising biomarkers. This study offers opportunities to study ATs and other ultra-short oligonucleotides in biological samples.

Alpha-fetoprotein: Past, present, and future

Alpha-fetoprotein (AFP) is a glycoprotein that plays an important role in immune regulation with critical involvement in early human development and maintaining the immune balance during pregnancy. Postfetal development, the regulatory mechanisms controlling AFP undergo a shift and AFP gene transcription is suppressed. Instead, these enhancers refocus their activity to maintain albumin gene transcription throughout adulthood. During the postnatal period, AFP expression can increase in the setting of hepatocyte injury, regeneration, and malignant transformation. It is the first oncoprotein discovered and is routinely used as part of a screening strategy for HCC. AFP has been shown to be a powerful prognostic biomarker, and multiple HCC prognosis models confirmed the independent prognostic utility of AFP. AFP is also a useful predictive biomarker for monitoring the treatment response of HCC. In addition to its role as a biomarker, AFP plays important roles in immune modulation to promote tumorigenesis and thus has been investigated as a therapeutic target in HCC. In this review article, we aim to provide an overview of AFP, encompassing the discovery, biological role, and utility as an HCC biomarker in combination with other biomarkers and how it impacts clinical practice and future direction.

Cuproptosis-associated genes (CAGs) contribute to the prognosis prediction and potential therapeutic targets in hepatocellular carcinoma

BACKGROUND

Cuproptosis is a novel form of cell death that exhibits close association with mitochondrial respiration and occurs through distinct mechanisms compared to previously characterized forms of cell death. However, the precise impact of cuproptosis-associated genes (CAGs) on prognosis, immune profiles, and treatment efficacy in hepatocellular carcinomas (HCC) remains poorly understood.

METHODS

A comprehensive analysis of CAGs in hepatocellular carcinoma (HCC) prognosis was conducted using genomic data from HCC patients. Consensus clustering analysis was performed to determine molecular subtypes related to cuproptosis in HCC. The single-sample gene set enrichment analysis (ssGSEA) algorithm was applied to quantify the infiltration levels of immune cells, while the "ESTIMATE" package was employed to calculate tumor purity, stromal scores, and immune scores in the tumor microenvironment (TME). Principal component analysis (PCA) algorithm was utilized to construct a risk score related to CAGs. Finally, CCK8, wound healing, Transwell migration/invasion, EDU and xenograft model were employed to explore the potential oncogenic role of MTF1.

RESULTS

Three distinct patterns of cuproptosis modification were identified, each associated with unique functional enrichments, clinical characteristics, immune cell infiltration, immune checkpoints, tumor microenvironment (TME), and prognosis. A CAGs-related risk score (Cuscore) was developed to predict prognosis in TCGA and validated in GSE76427 and ICGC datasets. Notably, patients with a low Cuscore had better prognoses and were more likely to benefit from immunotherapy.Additionally, the high Cuscore group in HCC also revealed three potential therapeutic targets (TUBA1B, CDC25B, and CSNK2A1) as well as several therapeutic compounds. Moreover, the experiment measured the expression levels of six prognosis-related CAGs, wherein knockdown of MTF1 exhibited suppression of proliferation, invasion, and migration formation in HCC cell lines.

CONCLUSION

The findings have enhanced our comprehension of the cuproptosis characteristics in HCC, and stratification based on CuScore may potentially enhance the prediction of patients' prognosis and facilitate the development of effective and innovative treatment strategies.

The significance of m6A RNA methylation regulators in diagnosis and subtype classification of HBV-related hepatocellular carcinoma

In recent years, abnormal m6A alteration in hepatocellular carcinoma (HCC) has been a focus on investigating the biological implications. In this study, our objective is to determine whether m6A modification contributes to the progression of HBV-related HCC. To achieve this, we employed a random forest model to screen top 8 characteristic m6A regulators from 19 candidate genes. Subsequently, we developed a nomogram model that utilizes these 8 characteristic m6A regulators to predict the prevalence of HBV-related HCC. According to decision curve analysis, patients may benefit from the nomogram model. The clinical impact curves exhibited a robust predictive capability of the nomogram models. Additionally, consensus molecular subtyping was employed to identify m6A modification patterns and m6A-related gene signature. The quantification of immune cell subsets was accomplished through the implementation of ssGSEA algorithms. PCA algorithms were developed to compute the m6A score for individual tumors. Two distinct m6A modification patterns, namely cluster A and cluster B, exhibited significant correlations with distinct immune infiltration patterns and biological pathways. Notably, patients belonging to cluster B demonstrated higher m6A scores compared to those in cluster A, as determined by the m6A score metric. Furthermore, the expression of IGFBP3 proteins was validated through immunofluorescence, revealing their pronounced lower expression in tumor tissues. In summary, our study underscores the importance of m6A modification in the advancement of HBV-related HCC. This research has the potential to yield novel prognostic biomarkers and therapeutic targets for the identification of HBV-related HCC.

Prognostic value and potential mechanism of cellular senescence and tumor microenvironment in hepatocellular carcinoma: Insights from bulk transcriptomics and single-cell sequencing analysis

The high mortality rate and postoperative recurrence of hepatocellular carcinoma (HCC) contribute to the burden on society and healthcare. The prognostic value and underlying mechanisms of cellular senescence and tumor microenvironment (TME) in HCC remain unclear. Bulk transcriptomic data were obtained from 368 HCC samples in The Cancer Genome Atlas-liver hepatocellular carcinoma cohort and 64 samples from the GSE116174 dataset. Single-cell RNA sequencing (scRNA-seq) data of HCC were obtained from the GSE149614 dataset, including 18 tumor samples from 10 patients. Prognosis-related cellular senescence genes and immune cells were identified through univariate analysis. Least absolute shrinkage and selection operator regression analysis was performed to construct the CellAge score and TME score, both of which were identified as independent prognostic factors for HCC based on multivariate Cox analysis. The combined CellAge and TME scores showed improved prognostic stratification for HCC patients, as confirmed by multivariate Cox analysis (p < .001). The gene set enrichment analysis (GSEA) revealed enrichment of the extracellular matrix receptor interaction signaling pathway in the group with high CellAge scores and low TME scores, which exhibited a worse prognosis. Single-cell sequencing results revealed higher expression activity of the cAMP response element modulator (CREM) extended transcription factor in HCC cells and most immune cells, indicating its involvement in TME remodeling. Finally, the tumor immune dysfunction and exclusion (TIDE) analysis demonstrated that the combined scores could predict the outcomes of immune therapy in patients with HCC. In conclusion, cellular senescence contributes to TME remodeling in HCC, and the developed CellAge and TME scores serve as independent prognostic factors and predictors of immune therapy in HCC.

Stratifying ICIs-responsive tumor microenvironment in HCC: from parsing out immune-hypoxic crosstalk to clinically applicable MRI-radiomics models

BACKGROUND

We aimed to redefine Immune checkpoint inhibitors (ICIs)-responsive "hot" TME and develop a corresponding stratification model to maximize ICIs-efficacy in Hepatocellular Carcinoma (HCC).

METHODS

Hypoxic scores were designed, and the relevance to immunotherapy responses were validated in pan-cancers through single cell analysis. Multi-omics analysis using the hypoxic scores and immune infiltrate abundance was performed to redefine the ICIs-responsive TME subtype in HCC patients from TCGA (n = 363) and HCCDB database (n = 228). The immune hypoxic stress index (IHSI) was constructed to stratify the ICIs-responsive TME subtype, with exploring biological mechanism in vitro and in vivo. MRI-radiomics models were built for clinical applicability.

RESULTS

The hypoxic scores were lower in the dominant cell-subclusters of responders in pan-cancers. The higher immune infiltrate-normoxic (HIN) subtype was redefined as the ICIs-responsive TME. Stratification of the HIN subtype using IHSI effectively identified ICIs-responders in Melanoma (n = 122) and urological cancer (n = 22). TRAF3IP3, the constituent gene of IHSI, was implicated in ICIs-relevant "immune-hypoxic" crosstalk by stimulating MAVS/IFN-I pathway under normoxic condition. MRI-radiomics models assessing TRAF3IP3 with HIF1A expression (AUC > 0.80) screened ICIs-Responders in HCC cohort (n = 75).

CONCLUSION

The hypoxic-immune stratification redefined ICIs-responsive TME and provided MRI-Radiomics models for initial ICIs-responders screening, with IHSI facilitating further identification.

Characterization of the distinct immune microenvironments between hepatocellular carcinoma and intrahepatic cholangiocarcinoma

As two major types of primary liver cancers, the tumor immune microenvironment (TIME) of hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (ICC) have been well studied separately. However, a systemic assessment of the similarities and differences between the TIME of HCC and ICC is still lacking. In this study, we pictured a landscape of combined TIME of HCC and ICC by sequencing and integrating 41 single-cell RNA-seq samples from four different tissue types of both malignancies. We found that T cells in HCC tumors generally exhibit higher levels of immunosuppression and exhaustion than those in ICC tumors. Myeloid cells in HCC and ICC tumors also exhibit distinct phenotypes and may serve as a key factor driving the differences between their TIMEs. Besides, we identified a cluster of EGR1+ macrophages specifically enriched in HCC tumors. Together, our study provides new insights into cellular composition, states and interactions in the TIMEs of HCC and ICC, which could pave the way for the development of future therapeutic targets for liver cancers.

Activation of CTU2 expression by LXR promotes the development of hepatocellular carcinoma

Cytosolic thiouridylase 2 (CTU2) is an enzyme modifying transfer RNAs post-transcriptionally, which has been implicated in breast cancer and melanoma development. And we found CTU2 participated in hepatocellular carcinoma (HCC) progression here. HepG2 cells as well as xenograft nude mice model were employed to investigate the role of CTU2 in HCC development in vitro and in vivo respectively. Further, we defined CTU2 as a Liver X receptor (LXR) targeted gene, with a typical LXR element in the CTU2 promoter. CTU2 expression was activated by LXR agonist and depressed by LXR knockout. Interestingly, we also found CTU2 took part in lipogenesis by directly enhancing the synthesis of lipogenic proteins, which provided a novel mechanism for LXR regulating lipid synthesis. Meanwhile, lipogenesis was active during cell proliferation, particularly in tumor cells. Reduction of CTU2 expression was related to reduced tumor burden and synergized anti-tumor effect of LXR ligands by inducing tumor cell apoptosis and inhibiting cell proliferation. Taken together, our study identified CTU2 as an LXR target gene. Inhibition of CTU2 expression could enhance the anti-tumor effect of LXR ligand in HCC, identifying CTU2 as a promising target for HCC treatment and providing a novel strategy for the application of LXR agonists in anti-tumor effect.

DNA-based molecular classifiers for the profiling of gene expression signatures

Although gene expression signatures offer tremendous potential in diseases diagnostic and prognostic, but massive gene expression signatures caused challenges for experimental detection and computational analysis in clinical setting. Here, we introduce a universal DNA-based molecular classifier for profiling gene expression signatures and generating immediate diagnostic outcomes. The molecular classifier begins with feature transformation, a modular and programmable strategy was used to capture relative relationships of low-concentration RNAs and convert them to general coding inputs. Then, competitive inhibition of the DNA catalytic reaction enables strict weight assignment for different inputs according to their importance, followed by summation, annihilation and reporting to accurately implement the mathematical model of the classifier. We validated the entire workflow by utilizing miRNA expression levels for the diagnosis of hepatocellular carcinoma (HCC) in clinical samples with an accuracy 85.7%. The results demonstrate the molecular classifier provides a universal solution to explore the correlation between gene expression patterns and disease diagnostics, monitoring, and prognosis, and supports personalized healthcare in primary care.

A hypoxia-glycolysis-lactate-related gene signature for prognosis prediction in hepatocellular carcinoma

BACKGROUND

Liver cancer ranks sixth in incidence and third in mortality globally and hepatocellular carcinoma (HCC) accounts for 90% of it. Hypoxia, glycolysis, and lactate metabolism have been found to regulate the progression of HCC separately. However, there is a lack of studies linking the above three to predict the prognosis of HCC. The present study aimed to identify a hypoxia-glycolysis-lactate-related gene signature for assessing the prognosis of HCC.

METHODS

This study collected 510 hypoxia-glycolysis-lactate genes from Molecular Signatures Database (MSigDB) and then classified HCC patients from TCGA-LIHC by analyzing their hypoxia-glycolysis-lactate genes expression. Differentially expressed genes (DEGs) were screened out to construct a gene signature by LASSO-Cox analysis. Univariate and multivariate regression analyses were used to evaluate the independent prognostic value of the gene signature. Analyses of immune infiltration, somatic cell mutations, and correlation heatmap were conducted by "GSVA" R package. Single-cell analysis conducted by "SingleR", "celldex", "Seurat", and "CellCha" R packages revealed how signature genes participated in hypoxia/glycolysis/lactate metabolism and PPI network identified hub genes.

RESULTS

We classified HCC patients from TCGA-LIHC into two clusters and screened out DEGs. An 18-genes prognostic signature including CDCA8, CBX2, PDE6A, MED8, DYNC1LI1, PSMD1, EIF5B, GNL2, SEPHS1, CCNJL, SOCS2, LDHA, G6PD, YBX1, RTN3, ADAMTS5, CLEC3B, and UCK2 was built to stratify the risk of HCC. The risk score of the hypoxia-glycolysis-lactate gene signature was further identified as a valuable independent factor for estimating the prognosis of HCC. Then we found that the features of clinical characteristics, immune infiltration, somatic cell mutations, and correlation analysis differed between the high-risk and low-risk groups. Furthermore, single-cell analysis indicated that the signature genes could interact with the ligand-receptors of hepatocytes/fibroblasts/plasma cells to participate in hypoxia/glycolysis/lactate metabolism and PPI network identified potential hub genes in this process: CDCA8, LDHA, YBX1.

CONCLUSION

The hypoxia-glycolysis-lactate-related gene signature we built could provide prognostic value for HCC and suggest several hub genes for future HCC studies.

New treatment alternatives for primary and metastatic colorectal cancer by an integrated transcriptome and network analyses

Metastatic colorectal cancer (CRC) is still in need of effective treatments. This study applies a holistic approach to propose new targets for treatment of primary and liver metastatic CRC and investigates their therapeutic potential in-vitro. An integrative analysis of primary and metastatic CRC samples was implemented for alternative target and treatment proposals. Integrated microarray samples were grouped based on a co-expression network analysis. Significant gene modules correlated with primary CRC and metastatic phenotypes were identified. Network clustering and pathway enrichments were applied to gene modules to prioritize potential targets, which were shortlisted by independent validation. Finally, drug-target interaction search led to three agents for primary and liver metastatic CRC phenotypes. Hesperadin and BAY-1217389 suppress colony formation over a 14-day period, with Hesperadin showing additional efficacy in reducing cell viability within 48 h. As both candidates target the G2/M phase proteins NEK2 or TTK, we confirmed their anti-proliferative properties by Ki-67 staining. Hesperadinin particular arrested the cell cycle at the G2/M phase. IL-29A treatment reduced migration and invasion capacities of TGF-β induced metastatic cell lines. In addition, this anti-metastatic treatment attenuated TGF-β dependent mesenchymal transition. Network analysis suggests IL-29A induces the JAK/STAT pathway in a preventive manner.

An exosome mRNA-related gene risk model to evaluate the tumor microenvironment and predict prognosis in hepatocellular carcinoma

BACKGROUND

The interplay between exosomes and the tumor microenvironment (TME) remains unclear. We investigated the influence of exosomes on the TME in hepatocellular carcinoma (HCC), focusing on their mRNA expression profile.

METHODS

mRNA expression profiles of exosomes were obtained from exoRBase. RNA sequencing data from HCC patients' tumors were acquired from The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC). An exosome mRNA-related risk score model of prognostic value was established. The patients in the two databases were divided into high- and low-risk groups based on the median risk score value, and used to validate one another. Functional enrichment analysis was performed based on a differential gene prognosis model (DGPM). CIBERSORT was used to assess the abundance of immune cells in the TME. The correlation between the expression levels of immune checkpoint-related genes and DGPM was analyzed alongside the prediction value to drug sensitivity.

RESULTS

A prognostic exosome mRNA-related 4-gene signature (DYNC1H1, PRKDC, CCDC88A, and ADAMTS5) was constructed and validated. A prognostic nomogram had prognostic ability for HCC. The genes for this model are involved in extracellular matrix, extracellular matrix (ECM)-receptor interaction, and the PI3K-Akt signaling pathway. Expression of genes here had a positive correlation with immune cell infiltration in the TME.

CONCLUSIONS

Our study results demonstrate that an exosome mRNA-related risk model can be established in HCC, highlighting the functional significance of the molecules in prognosis and risk stratification.

Th17 Cells Secrete TWEAK to Trigger Epithelial-Mesenchymal Transition and Promote Colorectal Cancer Liver Metastasis

Liver metastasis is the leading cause of mortality in patients with colorectal cancer. Given the significance of both epithelial-mesenchymal transition (EMT) of tumor cells and the immune microenvironment in colorectal cancer liver metastasis (CRLM), the interplay between them could hold the key for developing improved treatment options. We employed multiomics analysis of 130 samples from 18 patients with synchronous CRLM integrated with external datasets to comprehensively evaluate the interaction between immune cells and EMT of tumor cells in liver metastasis. Single-cell RNA sequencing analysis revealed distinct distributions of nonmalignant cells between primary tumors from patients with metastatic colorectal cancer (mCRC) and non-metastatic colorectal cancer, showing that Th17 cells were predominantly enriched in the primary lesion of mCRC. TWEAK, a cytokine secreted by Th17 cells, promoted EMT by binding to receptor Fn14 on tumor cells, and the TWEAK-Fn14 interaction enhanced tumor migration and invasion. In mouse models, targeting Fn14 using CRISPR-induced knockout or lipid nanoparticle-encapsulated siRNA alleviated metastasis and prolonged survival. Mice lacking Il17a or Tnfsf12 (encoding TWEAK) exhibited fewer metastases compared with wild-type mice, while cotransfer of Th17 with tumor cells promoted liver metastasis. Higher TWEAK expression was associated with a worse prognosis in patients with colorectal cancer. In addition, CD163L1+ macrophages interacted with Th17 cells, recruiting Th17 via the CCL4-CCR5 axis. Collectively, this study unveils the role of immune cells in the EMT process and identifies TWEAK secreted by Th17 as a driver of CRLM.

SIGNIFICANCE

TWEAK secreted by Th17 cells promotes EMT by binding to Fn14 on colorectal cancer cells, suggesting that blocking the TWEAK-Fn14 interaction may be a promising therapeutic approach to inhibit liver metastasis.

Scutellarin activates IDH1 to exert antitumor effects in hepatocellular carcinoma progression

Isochlorate dehydrogenase 1 (IDH1) is an important metabolic enzyme for the production of α-ketoglutarate (α-KG), which has antitumor effects and is considered to have potential antitumor effects. The activation of IDH1 as a pathway for the development of anticancer drugs has not been attempted. We demonstrated that IDH1 can limit glycolysis in hepatocellular carcinoma (HCC) cells to activate the tumor immune microenvironment. In addition, through proteomic microarray analysis, we identified a natural small molecule, scutellarin (Scu), which activates IDH1 and inhibits the growth of HCC cells. By selectively modifying Cys297, Scu promotes IDH1 active dimer formation and increases α-KG production, leading to ubiquitination and degradation of HIF1a. The loss of HIF1a further leads to the inhibition of glycolysis in HCC cells. The activation of IDH1 by Scu can significantly increase the level of α-KG in tumor tissue, downregulate the HIF1a signaling pathway, and activate the tumor immune microenvironment in vivo. This study demonstrated the inhibitory effect of IDH1-α-KG-HIF1a on the growth of HCC cells and evaluated the inhibitory effect of Scu, the first IDH1 small molecule agonist, which provides a reference for cancer immunotherapy involving activated IDH1.

Tyrosine hydroxylase inhibits HCC progression by downregulating TGFβ/Smad signaling

The alteration of metabolic processes has been found to have significant impacts on the development of hepatocellular carcinoma (HCC). Nevertheless, the effects of dysfunction of tyrosine metabolism on the development of HCC remains to be discovered. This research demonstrated that tyrosine hydroxylase (TH), which responsible for the initial and limiting step in the bio-generation of the neuro-transmitters dopamine and adrenaline, et al. was shown to be reduced in HCC. Increased expression of TH was found facilitates the survival of HCC patients. In addition, decreased TH indicated larger tumor size, much more numbers of tumor, higher level of AFP, and the presence of cirrhosis. TH effectively impairs the growth and metastasis of HCC cells, a process dependent on the phosphorylation of serine residues (S19/S40). TH directly binds to Smad2 and hinders the cascade activation of TGFβ/Smad signaling with the treatment of TGFβ1. In summary, our study uncovered the non-metabolic functions of TH in the development of HCC and proposes that TH might be a promising biomarker for diagnosis as well as an innovative target for metastatic HCC.

Lifestyle factors, glycemic traits, and lipoprotein traits and risk of liver cancer: a Mendelian randomization analysis

The current state of knowledge on the relationship between lifestyle factors, glycemic traits, lipoprotein traits with liver cancer risk is still uncertain despite some attempts made by observational studies. This study aims to investigate the causal genetic relationship between factors highly associated with liver cancer incidence by using Mendelian randomization (MR) analysis. Employing MR analysis, this study utilized previously published GWAS datasets to investigate whether lifestyle factors, glycemic traits, and lipoprotein traits would affect the risk of liver cancer. The study utilized three MR methods, including inverse variance-weighted model (IVW), MR Egger, and weighted median. Furthermore, MR-Egger analyses were performed to detect heterogeneity in the MR results. The study also conducted a leave-one-out analysis to assess the potential influence of individual SNPs on the MR analysis results. MR-PRESSO was used to identify and remove SNP outliers associated with liver cancer. MR analyses revealed that 2-h glucose (odds ratio, OR 2.33, 95% confidence interval, CI 1.28-4.21), type 2 diabetes mellitus (T2DM, OR 1.67, 95% CI 1.18-2.37), body mass index (BMI, OR 1.67, 95% CI 1.18-2.37), waist circumference (OR 1.78, 95% CI 1.18-2.37) were associated with increased risk of liver cancer. On the contrary, apolipoproteins B (APOB, OR 0.67, 95% CI 0.47-0.97), and low-density lipoprotein (LDL, OR 0.62, 95% CI 0.42-0.92) were negatively related to liver cancer risk. Additionally, after adjusting for BMI, apolipoproteins A-I (APOA-I, OR 0.56, 95% CI, 0.38-0.81), total cholesterol (TC, OR 0.72, 95% CI, 0.54-0.94), and total triglycerides (TG, OR 0.57, 95% CI, 0.40-0.78) exhibited a significant inverse correlation with the risk of liver cancer. This study supports a causal relationship between 2-h glucose, T2DM, BMI, and waist circumference with the increased risk of liver cancer. Conversely, the study reveals a cause-effect relationship between TC, TG, LDL, APOA-I, and APOB with a decreased risk of liver cancer.

Tumor phylogeography reveals block-shaped spatial heterogeneity and the mode of evolution in Hepatocellular Carcinoma

Solid tumors are complex ecosystems with heterogeneous 3D structures, but the spatial intra-tumor heterogeneity (sITH) at the macroscopic (i.e., whole tumor) level is under-explored. Using a phylogeographic approach, we sequence genomes and transcriptomes from 235 spatially informed sectors across 13 hepatocellular carcinomas (HCC), generating one of the largest datasets for studying sITH. We find that tumor heterogeneity in HCC segregates into spatially variegated blocks with large genotypic and phenotypic differences. By dissecting the transcriptomic heterogeneity, we discover that 30% of patients had a "spatially competing distribution" (SCD), where different spatial blocks have distinct transcriptomic subtypes co-existing within a tumor, capturing the critical transition period in disease progression. Interestingly, the tumor regions with more advanced transcriptomic subtypes (e.g., higher cell cycle) often take clonal dominance with a wider geographic range, rejecting neutral evolution for SCD patients. Extending the statistical tests for detecting natural selection to many non-SCD patients reveal varying levels of selective signal across different tumors, implying that many evolutionary forces including natural selection and geographic isolation can influence the overall pattern of sITH. Taken together, tumor phylogeography unravels a dynamic landscape of sITH, pinpointing important evolutionary and clinical consequences of spatial heterogeneity in cancer.

Identification of differentially expressed mRNA/lncRNA modules in acutely regorafenib-treated sorafenib-resistant Huh7 hepatocellular carcinoma cells

The multikinase inhibitor sorafenib is the standard first-line treatment for advanced hepatocellular carcinoma (HCC), but many patients become sorafenib-resistant (SR). This study investigated the efficacy of another kinase inhibitor, regorafenib (Rego), as a second-line treatment. We produced SR HCC cells, wherein the PI3K-Akt, TNF, cAMP, and TGF-beta signaling pathways were affected. Acute Rego treatment of these cells reversed the expression of genes involved in TGF-beta signaling but further increased the expression of genes involved in PI3K-Akt signaling. Additionally, Rego reversed the expression of genes involved in nucleosome assembly and epigenetic gene expression. Weighted gene co-expression network analysis (WGCNA) revealed four differentially expressed long non-coding RNA (DElncRNA) modules that were associated with the effectiveness of Rego on SR cells. Eleven putative DElncRNAs with distinct expression patterns were identified. We associated each module with DEmRNAs of the same pattern, thus obtaining DElncRNA/DEmRNA co-expression modules. We discuss the potential significance of each module. These findings provide insights and resources for further investigation into the potential mechanisms underlying the response of SR HCC cells to Rego.

Psilostachyin C reduces malignant properties of hepatocellular carcinoma cells by blocking CREBBP-mediated transcription of GATAD2B

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related mortality globally. Many herbal medicines and their bioactive compounds have shown anti-tumor properties. This study was conducted to examine the effect of psilostachyin C (PSC), a sesquiterpenoid lactone isolated from Artemisia vulgaris L., in the malignant properties of HCC cells. CCK-8, flow cytometry, wound healing, and Transwell assays revealed that 25 μM PSC treatment significantly suppressed proliferation, cell cycle progression, migration, and invasion of two HCC cell lines (Hep 3B and Huh7) while promoting cell apoptosis. Bioinformatics prediction suggests CREB binding protein (CREBBP) as a promising target of PSC. CREBBP activated transcription of GATA zinc finger domain containing 2B (GATAD2B) by binding to its promoter. CREBBP and GATAD2B were highly expressed in clinical HCC tissues and the acquired HCC cell lines, but their expression was reduced by PSC. Either upregulation of CREBBP or GATAD2B restored the malignant properties of HCC cells blocked by PSC. Collectively, this evidence demonstrates that PSC pocessess anti-tumor functions in HCC cells by blocking CREBBP-mediated transcription of GATAD2B.

S100A6 drives lymphatic metastasis of liver cancer via activation of the RAGE/NF-kB/VEGF-D pathway

Patients diagnosed with lymph node (LN) metastatic liver cancer face an exceedingly grim prognosis. In-depth analysis of LN metastatic patients' characteristics and tumor cells' interactions with human lymphatic endothelial cells (HLECs), can provide important biological and therapeutic insights. Here we identify at the single-cell level that S100A6 expression differs between primary tumor and their LN metastasis. Of particular significance, we uncovered the disparity in S100A6 expression between tumors and normal tissues is greater in intrahepatic cholangiocarcinoma (ICC) patients, frequently accompanied by LN metastases, than that in hepatocellular carcinoma (HCC), with rare occurrence of LN metastasis. Furthermore, in the infrequent instances of LN metastasis in HCC, heightened S100A6 expression was observed, suggesting a critical role of S100A6 in the process of LN metastasis. Subsequent experiments further uncovered that S100A6 secreted from tumor cells promotes lymphangiogenesis by upregulating the expression and secretion of vascular endothelial growth factor-D (VEGF-D) in HLECs through the RAGE/NF-kB/VEGF-D pathway while overexpression of S100A6 in tumor cells also augmented their migration and invasion. Taken together, these data reveal the dual effects of S100A6 in promoting LN metastasis in liver cancer, thus highlighting its potential as a promising therapeutic target.

Novel insights into small open reading frame-encoded micropeptides in hepatocellular carcinoma: A potential breakthrough

Traditionally, non-coding RNAs (ncRNAs) are regarded as a class of RNA transcripts that lack encoding capability; however, advancements in technology have revealed that some ncRNAs contain small open reading frames (sORFs) that are capable of encoding micropeptides of approximately 150 amino acids in length. sORF-encoded micropeptides (SEPs) have emerged as intriguing entities in hepatocellular carcinoma (HCC) research, shedding light on this previously unexplored realm. Recent studies have highlighted the regulatory functions of SEPs in the occurrence and progression of HCC. Some SEPs exhibit inhibitory effects on HCC, but others facilitate its development. This discovery has revolutionized the landscape of HCC research and clinical management. Here, we introduce the concept and characteristics of SEPs, summarize their associations with HCC, and elucidate their carcinogenic mechanisms in HCC metabolism, signaling pathways, cell proliferation, and metastasis. In addition, we propose a step-by-step workflow for the investigation of HCC-associated SEPs. Lastly, we discuss the challenges and prospects of applying SEPs in the diagnosis and treatment of HCC. This review aims to facilitate the discovery, optimization, and clinical application of HCC-related SEPs, inspiring the development of early diagnostic, individualized, and precision therapeutic strategies for HCC.

Single-cell multiomics reveals the interplay of clonal evolution and cellular plasticity in hepatoblastoma

Hepatoblastomas (HB) display heterogeneous cellular phenotypes that influence the clinical outcome, but the underlying mechanisms are poorly understood. Here, we use a single-cell multiomic strategy to unravel the molecular determinants of this plasticity. We identify a continuum of HB cell states between hepatocytic (scH), liver progenitor (scLP) and mesenchymal (scM) differentiation poles, with an intermediate scH/LP population bordering scLP and scH areas in spatial transcriptomics. Chromatin accessibility landscapes reveal the gene regulatory networks of each differentiation pole, and the sequence of transcription factor activations underlying cell state transitions. Single-cell mapping of somatic alterations reveals the clonal architecture of each tumor, showing that each genetic subclone displays its own range of cellular plasticity across differentiation states. The most scLP subclones, overexpressing stem cell and DNA repair genes, proliferate faster after neo-adjuvant chemotherapy. These results highlight how the interplay of clonal evolution and epigenetic plasticity shapes the potential of HB subclones to respond to chemotherapy.

Pharmacogenomic profiling of intra-tumor heterogeneity using a large organoid biobank of liver cancer

Inter- and intra-tumor heterogeneity is a major hurdle in primary liver cancer (PLC) precision therapy. Here, we establish a PLC biobank, consisting of 399 tumor organoids derived from 144 patients, which recapitulates histopathology and genomic landscape of parental tumors, and is reliable for drug sensitivity screening, as evidenced by both in vivo models and patient response. Integrative analysis dissects PLC heterogeneity, regarding genomic/transcriptomic characteristics and sensitivity to seven clinically relevant drugs, as well as clinical associations. Pharmacogenomic analysis identifies and validates multi-gene expression signatures predicting drug response for better patient stratification. Furthermore, we reveal c-Jun as a major mediator of lenvatinib resistance through JNK and β-catenin signaling. A compound (PKUF-01) comprising moieties of lenvatinib and veratramine (c-Jun inhibitor) is synthesized and screened, exhibiting a marked synergistic effect. Together, our study characterizes the landscape of PLC heterogeneity, develops predictive biomarker panels, and identifies a lenvatinib-resistant mechanism for combination therapy.

Loss of lncRNA LINC01056 leads to sorafenib resistance in HCC

BACKGROUND AND AIMS

Sorafenib is a major nonsurgical option for patients with advanced hepatocellular carcinoma (HCC); however, its clinical efficacy is largely undermined by the acquisition of resistance. The aim of this study was to identify the key lncRNA involved in the regulation of the sorafenib response in HCC.

MATERIALS AND METHODS

A clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) single-guide RNA (sgRNA) synergistic activation mediator (SAM)-pooled lncRNA library was applied to screen for the key lncRNA regulated by sorafenib treatment. The role of the identified lncRNA in mediating the sorafenib response in HCC was examined in vitro and in vivo. The underlying mechanism was delineated by proteomic analysis. The clinical significance of the expression of the identified lncRNA was evaluated by multiplex immunostaining on a human HCC microtissue array.

RESULTS

CRISPR/Cas9 lncRNA library screening revealed that Linc01056 was among the most downregulated lncRNAs in sorafenib-resistant HCC cells. Knockdown of Linc01056 reduced the sensitivity of HCC cells to sorafenib, suppressing apoptosis in vitro and promoting tumour growth in mice in vivo. Proteomic analysis revealed that Linc01056 knockdown in sorafenib-treated HCC cells induced genes related to fatty acid oxidation (FAO) while repressing glycolysis-associated genes, leading to a metabolic switch favouring higher intracellular energy production. FAO inhibition in HCC cells with Linc01056 knockdown significantly restored sensitivity to sorafenib. Mechanistically, we determined that PPARα is the critical molecule governing the metabolic switch upon Linc01056 knockdown in HCC cells and indeed, PPARα inhibition restored the sorafenib response in HCC cells in vitro and HCC tumours in vivo. Clinically, Linc01056 expression predicted optimal overall and progression-free survival outcomes in HCC patients and predicted a better sorafenib response. Linc01056 expression indicated a low FAO level in HCC.

CONCLUSION

Our study identified Linc01056 as a critical epigenetic regulator and potential therapeutic target in the regulation of the sorafenib response in HCC.

Proteomic analysis of DEN and CCl4-induced hepatocellular carcinoma mouse model

Hepatocellular carcinoma (HCC) seriously threatens human health, mostly developed from liver fibrosis or cirrhosis. Since diethylnitrosamine (DEN) and carbon tetrachloride (CCl4)-induced HCC mouse model almost recapitulates the characteristic of HCC with fibrosis and inflammation, it is taken as an essential tool to investigate the pathogenesis of HCC. However, a comprehensive understanding of the protein expression profile of this model is little. In this study, we performed proteomic analysis of this model to elucidate its proteomic characteristics. Compared with normal liver tissues, 432 differentially expressed proteins (DEPs) were identified in tumor tissues, among which 365 were up-regulated and 67 were down-regulated. Through Gene Ontology (GO) analysis, Ingenuity Pathway Analysis (IPA), protein-protein interaction networks (PPI) analysis and Gene-set enrichment analysis (GSEA) analysis of DEPs, we identified two distinguishing features of DEN and CCl4-induced HCC mouse model in protein expression, the upregulation of actin cytoskeleton and branched-chain amino acids metabolic reprogramming. In addition, matching DEPs from the mouse model to homologous proteins in the human HCC cohort revealed that the DEN and CCl4-induced HCC mouse model was relatively similar to the subtype of HCC with poor prognosis. Finally, combining clinical information from the HCC cohort, we screened seven proteins with prognostic significance, SMAD2, PTPN1, PCNA, MTHFD1L, MBOAT7, FABP5, and AGRN. Overall, we provided proteomic data of the DEN and CCl4-induced HCC mouse model and highlighted the important proteins and pathways in it, contributing to the rational application of this model in HCC research.

Depletion of proteasome subunit PSMD1 induces cancer cell death via protein ubiquitination and DNA damage, irrespective of p53 status

Hepatocellular carcinoma (HCC) is characterized by high incidence and fatality rates worldwide. In our exploration of prognostic factors in HCC, the 26s proteasome subunit, non-ATPase 1 (PSMD1) protein emerged as a significant contributor, demonstrating its potential as a therapeutic target in this aggressive cancer. PSMD1 is a subunit of the 19S regulatory particle in the 26S proteasome complex; the 19S particle controls the deubiquitination of ubiquitinated proteins, which are then degraded by the proteolytic activity of the complex. Proteasome-targeting in cancer therapy has received significant attention because of its practical application as an established anticancer agent. We investigated whether PSMD1 plays a critical role in cancer owing to its prognostic significance. PSMD1 depletion induced cell cycle arrest in G2/M phase, DNA damage and apoptosis of cancer cells, irrespective of the p53 status. PSMD1 depletion-mediated cell death was accompanied by an increase in overall protein ubiquitination. These phenotypes occurred exclusively in cancer cells, with no effects observed in normal cells. These findings indicate that PSMD1 depletion-mediated ubiquitination of cellular proteins induces cell cycle arrest and eventual death in cancer cells, emphasizing PSMD1 as a potential therapeutic target in HCC.

Molecular Insights of Cholestasis in MDR2 Knockout Murine Liver Organoids

MDR3 (multidrug resistance 3) deficiency in humans (MDR2 in mice) causes progressive familial intrahepatic cholestasis type 3 (PFIC3). PFIC3 is a lethal disease characterized by an early onset of intrahepatic cholestasis progressing to liver cirrhosis, a preneoplastic condition, putting individuals at risk of hepatocellular carcinoma (HCC). Hepatocyte-like organoids from MDR2-deficient mice (MDR2KO) were used in this work to study the molecular alterations caused by the deficiency of this transporter. Proteomic analysis by mass spectrometry allowed characterization of 279 proteins that were differentially expressed in MDR2KO compared with wild-type organoids. Functional enrichment analysis indicated alterations in three main cellular functions: (1) interaction with the extracellular matrix, (2) remodeling intermediary metabolism, and (3) cell proliferation and differentiation. The affected cellular processes were validated by orthogonal molecular biology techniques. Our results point to molecular mechanisms associated with PFIC3 that may drive the progression to liver cirrhosis and HCC and suggest proteins and cellular processes that could be targeted for the development of early detection strategies for these severe liver diseases.

Cyclophilin D knockout significantly prevents HCC development in a streptozotocin-induced mouse model of diabetes-linked NASH

A family of Peptidyl-prolyl isomerases (PPIases), called Cyclophilins, localize to numerous intracellular and extracellular locations where they contribute to a variety of essential functions. We previously reported that non-immunosuppressive pan-cyclophilin inhibitor drugs like reconfilstat (CRV431) or NV556 decreased multiple aspects of non-alcoholic fatty liver disease (NAFLD) in mice under two different non-alcoholic steatohepatitis (NASH) mouse models. Both CRV431 and NV556 inhibit several cyclophilin isoforms, among which cyclophilin D (CypD) has not been previously investigated in this context. It is unknown whether it is necessary to simultaneously inhibit multiple cyclophilin family members to achieve therapeutic benefits or if loss-of-function of one is sufficient. Furthermore, narrowing down the isoform most responsible for a particular aspect of NAFLD/NASH, such as hepatocellular carcinoma (HCC), would allow for more precise future therapies. Features of human diabetes-linked NAFLD/NASH can be reliably replicated in mice by administering a single high dose of streptozotocin to disrupt pancreatic beta cells, in conjunction with a high sugar, high fat, high cholesterol western diet over the course of 30 weeks. Here we show that while both wild-type (WT) and Ppif-/- CypD KO mice develop multipe severe NASH disease features under this model, the formation of HCC nodules was significantly blunted only in the CypD KO mice. Furthermore, of differentially expressed transcripts in a qPCR panel of select HCC-related genes, nearly all were downregulated in the CypD KO background. Cyclophilin inhibition is a promising and novel avenue of treatment for diet-induced NAFLD/NASH. This study highlights the impact of CypD loss-of-function on the development of HCC, one of the most severe disease outcomes.

DLK1/DIO3 locus upregulation by a β-catenin-dependent enhancer drives cell proliferation and liver tumorigenesis

The CTNNB1 gene, encoding β-catenin, is frequently mutated in hepatocellular carcinoma (HCC, ∼30%) and in hepatoblastoma (HB, >80%), in which DLK1/DIO3 locus induction is correlated with CTNNB1 mutations. Here, we aim to decipher how sustained β-catenin activation regulates DLK1/DIO3 locus expression and the role this locus plays in HB and HCC development in mouse models deleted for Apc (ApcΔhep) or Ctnnb1-exon 3 (β-cateninΔExon3) and in human CTNNB1-mutated hepatic cancer cells. We identified an enhancer site bound by TCF-4/β-catenin complexes in an open conformation upon sustained β-catenin activation (DLK1-Wnt responsive element [WRE]) and increasing DLK1/DIO3 locus transcription in β-catenin-mutated human HB and mouse models. DLK1-WRE editing by CRISPR-Cas9 approach impaired DLK1/DIO3 locus expression and slowed tumor growth in subcutaneous CTNNB1-mutated tumor cell grafts, ApcΔhep HB and β-cateninΔExon3 HCC. Tumor growth inhibition resulted either from increased FADD expression and subsequent caspase-3 cleavage in the first case or from decreased expression of cell cycle actors regulated by FoxM1 in the others. Therefore, the DLK1/DIO3 locus is an essential determinant of FoxM1-dependent cell proliferation during β-catenin-driven liver tumorigenesis. Targeting the DLK1-WRE enhancer to silence the DLK1/DIO3 locus might thus represent an interesting therapeutic strategy to restrict tumor growth in primary liver cancers with CTNNB1 mutations.

Pre-enrichment-free detection of hepatocellular carcinoma-specific ctDNA via PDMS and MEMS-based microfluidic sensor

The growing interest in microfluidic biosensors has led to improvements in the analytical performance of various sensing mechanisms. Although various sensors can be integrated with microfluidics, electrochemical ones have been most commonly employed due to their ease of miniaturization, integration ability, and low cost, making them an established point-of-care diagnostic method. This concept can be easily adapted to the detection of biomarkers specific to certain cancer types. Pathological profiling of hepatocellular carcinoma (HCC) is heterogeneous and rather complex, and biopsy samples contain limited information regarding the tumor and do not reflect its heterogeneity. Circulating tumor DNAs (ctDNAs), which can contain information regarding cancer characteristics, have been studied tremendously since liquid biopsy emerged as a new diagnostic method. Recent improvements in the accuracy and sensitivity of ctDNA determination also paved the way for genotyping of somatic genomic alterations. In this study, three-electrode (Au-Pt-Ag) glass chips were fabricated and combined with polydimethylsiloxane (PDMS) microchannels to establish an electrochemical microfluidic sensor for detecting c.747G > T hotspot mutations in the TP53 gene of ctDNAs from HCC. The preparation and analysis times of the constructed sensor were as short as 2 h in total, and a relatively high flow rate of 30 µl/min was used during immobilization and hybridization steps. To the best of our knowledge, this is the first time a PDMS-based microfluidic electrochemical sensor has been developed to target HCC ctDNAs. The system exhibited a limit of detection (LOD) of 24.1 fM within the tested range of 2-200 fM. The sensor demonstrated high specificity in tests conducted with fully noncomplementary and one-base mismatched target sequences. The developed platform is promising for detecting HCC-specific ctDNA at very low concentrations without requiring pre-enrichment steps.

Ropivacaine combined with sorafenib attenuates hepatocellular carcinoma cell proliferation and metastasis by inhibiting the miR-224/HOXD10 axis

OBJECTIVE

The development of drug resistance in hepatocellular carcinoma (HCC) cells limits the effectiveness of sorafenib (Sor). However, the regulatory mechanisms underlying the effects of the combination Sor and ropivacaine (Rop) on HCC cells remain unclear.

METHODS

miR-224 and HOXD10 mRNA expression in HCC cells was analyzed using qRT-PCR. CCK-8, Transwell assays and tumor formation experiments in nude mice were used to assess HCC cell proliferation, migration, and invasion. Migration of HCC cells was also analyzed using a cell scratch assay. Hematoxylin and eosin staining was used to detect tumor area.

RESULTS

miR-224 expression profoundly increased in HepG2 and Huh7 cells. Treatment with Rop and/or Sor blocked miR-244 expression, especially the combination treatment. Transfection of miR-224 mimic increased HCC cell proliferation and tumor size in nude mice, and migration and invasion in vitro in the presence of Rop and Sor compared to the negative control mimic. Dual-luciferase reporter assays showed that HOXD10 was targeted by miR-224. HOXD10 protein expression and was markedly reduced in HepG2 and Huh7 cells. Rop and/or Sor treatment increased HOXD10 protein expression, particularly the combination treatment. miR-224 negatively regulated HOXD10 expression in HCC cells treated with Rop and Sor. Transfection-mediated silencing of HOXD10 increased HCC cell proliferation, migration, and invasion in the presence of Rop and Sor compared with negative control transfection.

CONCLUSION

The combination of Rop and Sor attenuates HCC cell proliferation and metastasis via the miR-224/HOXD10 axis.

Local Therapy Improves Survival for Early Recurrence After Resection of Colorectal Liver Metastases

BACKGROUND

Early recurrence following hepatectomy for colorectal liver metastases (CLM) is associated with worse survival; yet, impact of further local therapy is unclear. We sought to evaluate whether local therapy benefits patients with early recurrence following hepatectomy for CLM.

METHODS

Clinicopathologic and survival outcomes of patients managed with hepatectomy for CLM (1/2001-12/2020) were queried from a prospectively maintained database. Timing of recurrence was stratified as early (recurrence-free survival [RFS] < 6 months), intermediate (RFS 6-12 months), and later (RFS > 12 months). Local therapy was defined as ablation, resection, or radiation.

RESULTS

Of 671 patients, 541 (81%) recurred with 189 (28%) early, 180 (27%) intermediate, and 172 (26%) later recurrences. Local therapy for recurrence resulted in improved survival, regardless of recurrence timing (early 78 vs. 32 months, intermediate 72 vs. 39 months, later 132 vs. 65 months, all p < 0.001). Following recurrence, treatment with local therapy (hazard ratio [HR] = 0.24), liver and extrahepatic recurrence (HR = 1.81), RAS + TP53 co-mutation (HR = 1.52), and SMAD4 mutation (HR = 1.92) were independently associated with overall survival (all p ≤ 0.002). Among patients with recurrence treated by local therapy, patients older than 65 years (HR 1.79), liver and extrahepatic recurrence (HR 2.05), primary site or other recurrence (HR 1.90), RAS-TP53 co-mutation (HR 1.63), and SMAD4 mutation (HR 2.06) had shorter post-local therapy survival (all p ≤ 0.04).

CONCLUSIONS

While most patients recur after hepatectomy for CLM, local therapy may result in long-term survival despite early recurrence. Somatic mutational profiling may help to guide the multidisciplinary consideration of local therapy after recurrence.