Published online Jan 24, 2024. doi: 10.5306/wjco.v15.i1.62
Peer-review started: October 23, 2023
First decision: November 23, 2023
Revised: December 12, 2023
Accepted: December 28, 2023
Article in press: December 28, 2023
Published online: January 24, 2024
Processing time: 91 Days and 12.6 Hours
Transcatheter arterial embolisation (TACE) is the primary treatment for intermediate-stage hepatocellular carcinoma (HCC) patients while some HCC cases have shown resistance to TACE.
To investigate the key genes and potential mechanisms correlated with TACE refractoriness in HCC.
The microarray datasets of TACE-treated HCC tissues, HCC and non-HCC tissues were collected by searching multiple public databases. The respective differentially expressed genes (DEGs) were attained via limma R package. Weighted gene co-expression network analysis was employed for identifying the significant modules related to TACE non-response. TACE refractoriness-related genes were obtained by intersecting up-regulated TACE-associated and HCC-associated DEGs together with the genes in significant modules related to TACE non-response. The key genes expression in the above two pairs of samples was compared respectively via Wilcoxon tests and standard mean differences model. The prognostic value of the key genes was evaluated by Kaplan-Meier curve. Multivariate analysis was utilised to investigate the independent prognostic factor in key genes. Single-cell RNA (scRNA) sequencing analysis was conducted to explore the cell types in HCC. TACE refractoriness-related genes activity was calculated via AUCell packages. The CellChat R package was used for the investigation of the cell–cell communication between the identified cell types.
HCC tissues of TACE non-responders (n = 66) and TACE responders (n = 81), HCC (n = 3941) and non-HCC (n = 3443) tissues were obtained. The five key genes, DLG associated protein 5 (DLGAP5), Kinesin family member 20A (KIF20A), Assembly factor for spindle microtubules (ASPM), Kinesin family member 11 (KIF11) and TPX2 microtubule nucleation factor (TPX2) in TACE refractoriness-related genes, were identified. The five key genes were all up-regulated in the TACE non-responders group and the HCC group. High expression of the five key genes predicted poor prognosis in HCC. Among the key genes, TPX2 was an independent prognostic factor. Four cell types, hepatocytes, embryonic stem cells, T cells and B cells, were identified in the HCC tissues. The TACE refractoriness-related genes expressed primarily in hepatocytes and embryonic stem cells. Hepatocytes, as the providers of ligands, had the strongest interaction with embryonic stem cells that provided receptors.
Five key genes (DLGAP5, KIF20A, ASPM, KIF11 and TPX2) were identified as promoting refractory TACE. Hepatocytes and embryonic stem cells were likely to boost TACE refractoriness.
Core Tip: This is a study that explored the key genes and mechanisms related to transcatheter arterial chemoembolisation (TACE) refractoriness in hepatocellular carcinoma. Through the combination of tissue microarrays and RNA-seq with single-cell RNA sequencing, the TACE refractoriness-related genes were identified and five key genes (DLGAP5, KIF20A, ASPM, KIF11 and TPX2) associated with TACE refractoriness were revealed. The TACE refractoriness-related genes were found to mainly express in hepatocytes and embryonic stem cells. Hepatocytes providing ligands had the strongest interaction with embryonic stem cells as receptors providers.