eISSN: 1897-4309
ISSN: 1428-2526
Contemporary Oncology/Współczesna Onkologia
Current issue Archive Manuscripts accepted About the journal Supplements Addendum Special Issues Editorial board Reviewers Abstracting and indexing Subscription Contact Instructions for authors Publication charge Ethical standards and procedures
Editorial System
Submit your Manuscript
SCImago Journal & Country Rank
3/2023
vol. 27
 
Share:
Share:
Original paper

NamiRNA-mediated high expression of KNSTRN correlates with poor prognosis and immune infiltration in hepatocellular carcinoma

Liang Jin
1
,
Xiaojing Zhang
1
,
Ming Fan
1
,
Weimin Li
1
,
Xuan Zhang
1

  1. Department of Hepatobiliary Surgery, Xijing Hospital, Air Force Medical University, Xi’an, Shaanxi, China
Contemp Oncol (Pozn) 2023; 27 (3): 163–175
Online publish date: 2023/12/06
Get citation
 
PlumX metrics:
 

Introduction

Hepatocellular carcinoma (HCC) is the fifth most common cancer and the second leading cause of cancer-related death worldwide [13]. The risk factors associated with the initiation and progression of HCC include viral infection [47], alcohol abuse [8], aflatoxin-contaminated foods [9], and immune system imbalances [10]. The overall 5-year survival rate for patients with HCC is less than 10% [11]. To improve patient outcomes, typical combination treatments, such as surgery, radiotherapy, transarterial chemoembolization, targeted therapy, and immunotherapy, have been investigated. Currently, surgical resection can only be performed in patients with early- stage HCC. The comprehensive treatment modality is rapidly becoming a standard treatment for liver cancer, which has been effective, but results remain unsatisfactory. Understanding the mechanisms and prognostic biomarkers of effective antitumour responses is essential for developing novel immunotherapeutic and targeted agents.

Kinetochore-localized astrin/sperm-associated antigen 5 binding protein (KNSTRN/SPAG5), a basic component of cell progression, plays an important role in cell division [12, 13]. Although most studies have focused on the association between KNSTRN and mitosis [14, 15], there are also some reports on the relationship between KNSTRN and tumourigenesis; however, it remains rarely studied [16, 17]. The oncogene, KNSTRN, appears to be mutated upon exposure to ultraviolet light [17]. Mutations of KNSTRN can interfere with chromatid cohesion, increase aneuploidy in tumours, and enhance tumourigenesis [17, 18]. Because few studies have been conducted on KNSTRN in cancer, we performed comprehensive and integrated bioinformatics analysis to determine its biological functions and potential mechanisms in HCC. Tumour immunotherapy, DNA methylation, and nuclear-activating miRNA (NamiRNA) regulatory networks are new directions in tumour gene therapy and are widely used to study the mechanism of liver hepatocellular carcinoma (LIHC). A recent study showed that KNSTRN expression was correlated with diverse immune infiltration levels in lung adenocarcinoma (LUAD) [19]. Therefore, it is possible that aberrant KNSTRN expression is related to immune escape and cancer development. In particular, the concept of NamiRNAs has been proposed in recent years, with encouraging results [2022]. However, the NamiRNA-enhancer-target gene activation network requires further attention.

Therefore, we comprehensively evaluated the association between KNSTRN expression and the prognosis of patients with LIHC using The Cancer Genome Atlas (TCGA) database. We investigated the correlation between KNSTRN expression and immune cell infiltration, immune cell biomarkers, and immune checkpoints in HCC. Finally, the NamiRNA-associated regulation of KNSTRN was explored, which would undoubtedly facilitate the identification of the underlying mechanisms and the development of new therapeutic methods for HCC. Collectively, our findings reveal the potential role of KNSTRN in patients with HCC and provide a clear understanding of the mechanisms underlying the association between KNSTRN and tumour immunity.

Material and methods

Data source

Messenger RNA (mRNA) expression data (including 374 HCC and 50 adjacent non-tumour tissues) and clinical information were downloaded from the TCGA database (http://cancergenome.nih.gov). Samples with RNA-seq data or missing related clinical data were excluded from the analysis.

Survival analysis

To investigate the prognostic relevance of differential gene and NamiRNA expression, overall survival (OS) and progression-free interval (PFI) analyses were performed. We used the “survival” (v3.2–10) R package for statistical analysis of survival data and the “survminer” (v0.4.9) R package for visualization to conduct survival analysis. Forest plots were generated using the “ggplot2” (v3.3.3) R package. Meanwhile, we also used the “ggplot2” (v3.3.3) (visualization) and the “pROC” (v1.17.0.1) (statistical analysis) R packages to investigate the predictive value of KNSTRN expression levels in HCC diagnosis. Statistical significance was set at a p-value < 0.05 for the log-rank test.

Analysis of the DNA methylation status of the CpG islands of KNSTRN

MethSurv (https://biit.cs.ut.ee/methsurv/) is a website that performs survival analyses based on DNA methylation biomarkers using the TCGA data. KNSTRN DNA methylation at CpG sites and the prognostic value of these CpG sites in HCC were analysed using MethSurv.

Functional enrichment and protein-protein interaction network analysis

The LinkedOmics database (http://www.linkedomics.org/login.php) is a web-based platform that provides comprehensive multi-omics data analysis tools for the TCGA database [23]. The Pearson correlation coefficient was calculated, and a volcano map of KNSTRN-related genes was plotted using the LinkedOmics website.

Hepatocellular carcinoma C patients from the TCGA dataset were divided into 2 groups (KNSTRNlow and KNSTRNhigh) according to the median values of KNSTRN mRNA expression from the RNAseq data. The “DESeq2” (v1.36.0) and “edgeR” (v3.38.2) packages in R 4.2.1 were used to screen differentially expressed genes (DEGs) between the KNSTRNlow and KNSTRNhigh groups of patients with HCC and plot a volcano map. The “ggplot2” (v3.3.6) and “VennDiagram” (v1.7.3) packages in R 4.2.1 were applied to draw Venn diagrams to determine the overlapping genes between DEGs and KNSTRN-related genes for further enrichment analysis. The “ClusterProfiler” (v3.14.3) R package was employed to analyse the gene ontology (GO) functions and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment of the overlapping genes. Immediately afterward, the “ggplot2” (v3.3.3) R package was used to visualize the above analysis data.

Protein-protein interaction (PPI) network analysis was performed using the STRING database. STRING (www.string-db.org) is an online analysis platform that includes all the publicly available PPI data [24].

Analysis of immune characteristics

Single-sample Gene Set Enrichment Analysis (ssGSEA) was performed using the “GSVA” (v1.34.0) R package [25, 26] to analyse infiltration by 24 immune cells, including neutrophils, cytotoxic cells, dendritic cells (DCs), CD8+ T-cells, plasmacytoid DC (pDC), natural killer (NK) cells, mast cells, T-γ δ (Tgd), type 17 Th (Th17) cells, immature DCs (iDCs), eosinophils, NK CD56dim cells, regulatory T-cells (Tregs), T-effector memory, T-cells, T-central memory, B-cells, type 1 Th (Th1) cells, macrophages, NK CD56 bright cells, activated DC, T-follicular helper (TFH), T-helper cells, and type 2 Th (Th2) cells. In addition, we used the “ggplot2” (v3.3.3) R package to analyse the correlation between KNSTRN expression and T-cell exhaustion markers (including PDCD1, CD274, CTLA4, HAVCR2, LAG3, and TOX). Pearson rank correlation analysis was applied to evaluate partial correlations and p-values. Correlations with p < 0.05 were considered significant.

We also used the TIMER2.0 (http://timer.cistrome.org/) server to analyse the correlation between KNSTRN expression and cancer-associated fibroblasts (CAFs) and tumour infiltration of 3 immunosuppressive cell types [27], including myeloid-derived suppressor cells (MDSCs), the M2 subtype of tumour-associated macrophages (M2-TAMs), and Tregs.

NamiRNA-enhancer-KNSTRN activation network construction

The LinkedOmics (http://www.linkedomics.org/login.php) website was used to obtain KNSTRN-associated miRNAs. Additionally, the TCGA data were analysed using the “DSEeq2”(v1.36.0) and “edgeR” (v3.38.2) packages in R to obtain the differentially expressed miRNAs. Additionally, if the miRNA transcriptional locus overlapped with the enhancer regions, the miRNA was chosen as the NamiRNA candidate, as described previously [20]. Overlapping miRNAs, as shown in the Venn diagram, among the KNSTRN-associated miRNAs, differentially expressed miRNAs, and enhancer-associated miRNAs were regarded as potential NamiRNAs. Comprehensive analysis of positively correlated NamiRNA-mRNA expression levels was performed to establish a key NamiRNA-enhancer-KNSTRN network for LIHC.

Statistical analysis

R software version 4.2.1 (https://cran.r-project.org/) was the main tool used to conduct the statistical analysis, and p < 0.05 was considered statistically significant. Additionally, other statistical analyses were automatically performed using the online databases mentioned above.

Results

Expression levels of KNSTRN in normal and cancerous human tissues

KNSTRN expression was analysed in 33 cancer datasets from the TCGA database. KNSTRN expression was significantly upregulated in 26 of the 33 cancer types (Fig. 1 A). As shown in Figure 1 A, KNSTRN expression was significantly downregulated in testicular germ cell tumours and thyroid carcinoma. No KNSTRN expression was detected in the adjacent normal tissues of mesothelioma (MESO) and uveal melanoma. KNSTRN expression was significantly higher in HCC tissues than in normal liver tissues (p < 0.001) (Fig. 1 B). In addition, receiver operating characteristic (ROC) curves were used to analyse the diagnostic value of KNSTRN expression. The area under the curve (AUC) for KNSTRN was 0.935, suggesting that KNSTRN is a potential diagnostic biomarker (Fig. 1 C).

Fig. 1

KNSTRN expression significantly up-regulated in multiple cancers, including liver hepatocellular carcinoma. The Cancer Genome Atlas (TCGA) database analysis shows the KNSTRN expression levels in 33 types of cancer tissues and their corresponding adjacent normal liver tissues. ns, p ≥ 0.05; **, p < 0.01; ***, p < 0.001 (A), KNSTRN expression levels were significantly higher in the hepatocellular carcinoma tissues compared to the normal liver tissues in the TCGA liver hepatocellular carcinoma (LIHC) datasets (B), diagnostic value of KNSTRN expression in LIHC (C)

/f/fulltexts/WO/51979/WO-27-51979-g001_min.jpg

Correlation analysis between KNSTRN expression and prognosis and clinicopathological characteristics

To better understand the role of KNSTRN expression in cancer, we used clinical data and TCGA RNAseq data (downloaded from UCSC Xena) to examine the prognosis of all TCGA cancer types to determine whether high KNSTRN expression influences the prognosis of patients with cancer. Two prognostic parameters, OS and PFI, were included. For OS, high expression of KNSTRN in pancreatic adenocarcinoma (PAAD), lung adenocarcinoma (LUAD), Adrenocortical carcinoma (ACC), bladder urothelial carcinoma (BLCA), Cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), kidney renal papillary cell carcinoma (KIRP), Lower Grade Glioma (LGG), Mesothelioma (MESO), Uterine Corpus Endometrial Carcinoma (UCEC), Glioma (GBM/LGG), oral squamous cell carcinoma (OSCC), and Liver Hepatocellular Carcinoma (LIHC) had unfavorable prognosis, but patients with thymoma (THYM) and Esophageal Squamous Cell Carcinoma (ESCC) with higher expression of KNSTRN displayed better prognosis (Fig. 2 A). For PFI, increased expression of KNSTRN indicated a poor prognosis in all affected tumours (Fig. 2 B). By combining the OS and PFI, KNSTRN may be utilized as a prognostic biomarker for unfavourable outcomes in patients with HCC. Detailed correlation analysis results revealed that elevated expression of KNSTRN was considerably associated with poor OS (HR = 1.48, 95% CI: 1.05–2.09, p = 0.027) and PFI (HR = 1.41, 95% CI: 1.05–1.89, p = 0.021) in HCC, as shown in Figures 3 A, B.

Fig. 2

Correlation analysis between KNSTRN expression and the prognostic value. The link between the expression of KNSTRN and prognosis in patients (overall survival and progression-free interval) of various cancers in The Cancer Genome Atlas database (A, B)

/f/fulltexts/WO/51979/WO-27-51979-g002_min.jpg

Additionally, ROC curves were used to determine the predictive value of KNSTRN mRNA levels. The area under the ROC curve was evaluated to predict the 1-, 2-, and 5-year OS risks in patients with HCC (1-year, AUC = 0.687; 2-year, AUC = 0.621; 5-year, AUC = 0.634) (Fig. 3 C). For PFI, KNSTRN expression also predicted the risk of HCC (1-year, AUC = 0.655; 2-year, AUC = 0.638; 5-year, AUC = 0.659) (Fig. 3 D).

Fig. 3

Correlation Analysis between KNSTRN expression and the prognostic value and clinicopathological characteristics. Survival curves with statistical significance in The Cancer Genome Atlas for overall survival (OS) (A) and progression-free interval (PFI) (B), the predictive value of KNSTRN mRNA levels for OS (C), PFI (D), relationship of KNSTRN expression with AFP (E), histological grade (F), pathologic stage (G), T stage (H), and TP53 mutation (I)

* p < 0.05

** p < 0.01

*** p < 0.001

/f/fulltexts/WO/51979/WO-27-51979-g003_min.jpg

We further assessed the expression of KNSTRN in various clinical phases of HCC and found that KNSTRN expression levels were significantly related to AFP, histological grade, pathologic stage, and T stage, as shown in Figu- res 3 E–H. KNSTRN was more highly expressed in the TP53-mutant group than in the TP53-nonmutant group (p < 0.001) (Fig. 3 I). Therefore, KNSTRN may promote HCC progression.

KNSTRN methylation status and its prognostic value

DNA methylation of KNSTRN and the prognostic value of CpG islands in KNSTRN were evaluated using the MetSurv tool. Ten CpG islands, including cg08036289 and cg25710630, exhibited elevated levels of DNA methylation (Fig. 4 A). Correlation analysis indicated that KNSTRN expression was significantly negatively correlated with methylation status (cg08036289) (Fig. 4 B). Furthermore, the methylation levels of 4 CpG islands, cg01301233 (HR = 0.639), cg08036289 (HR = 0.693), cg17805752 (HR = 1.662), and cg25947121 (HR = 1.979), were associated with prognosis (p < 0.05) (Table 1). Elevated levels of KNSTRN methylation in these 4 CpG islands, particularly cg08036289, were associated with better OS in patients with HCC than that in those with low levels of KNSTRN CpG methylation.

Fig. 4

DNA methylation levels in the KNSTRN gene are associated with the prognosis of hepatocellular carcinoma patients. The heat map of DNA methylation at CpG sites in the KNSTRN gene by the MethSurv database (A), correlation analysis of KNSTRN mRNA expression with KNSTRN methylation status by The Cancer Genome Atlas liver hepatocellular carcinoma database (B)

/f/fulltexts/WO/51979/WO-27-51979-g004_min.jpg
Table 1

Effects of methylation levels in the CpG sites of the KNSTRN gene on the prognosis of hepatocellular carcinoma patients

CpG islandHRp-value
Body-cg08036289#0.6930.042462946*
Body;3’UTR-cg25710630#1.3440.109162787
Body-cg059933860.8140.274468019
1stExon;5’UTR-cg156745960.7610.125774446
1stExon;5’UTR-cg178057521.6620.01790666*
Body-cg013012330.6390.011333953*
Body-cg262315820.7790.172668362
Body-cg259471211.9790.005027299*
TSS200-cg047210870.8530.383155881
TSS200-cg199476110.8260.288088112

HCC – hepatocellular carcinoma, HR – hazard ratio

# indicating a high level of KNSTRN methylation

* p < 0.05

KNSTRN functional clustering and interaction network analysis of KNSTRN-related genes

To further determine the biological significance of KNSTRN in HCC, the LinkedOmics database was used to determine the co-expression pattern of KNSTRN in the TCGA-LIHC cohort. As shown in Figure 5 A, 5655 genes were positively correlated with KNSTRN, and 2661 genes were significantly negatively correlated with KNSTRN (false discovery rate, FDR < 0.05). Data from TCGA were analysed using the DSEeq2 package in R (|logFC| > 1.5, adjusted p-value < 0.05), and 1120 DEGs were identified between the groups with high and low KNSTRN expression, including 870 upregulated and 250 downregulated genes (Fig. 5 B). Detailed descriptions of the related genes and DEGs are provided in Supplementary 1, 2. The overlapping KNSTRN-related genes and DEGs are shown in Figure 5 C.

Fig. 5

Kinetochore-localized astrin/SPAG5-binding protein (KNSTRN) functional clustering and interaction network analysis of KNSTRN-related genes. Volcano map of KNSTRN-related genes in hepatocellular carcinoma (HCC) by the LinkedOmics database (A), volcano plot of differentially expressed genes (DEGs) between the high and low KNSTRN expression groups (B), Venn diagram of the KNSTRN-related genes and DEGs in HCC (C), gene ontology enrichment (D), and Kyoto Encyclopedia of Genes and Genomes pathway analysis (E) by the overlapping genes of the KNSTRN-related genes and DEGs, network of co-expressed genes of KNSTRN (F)

/f/fulltexts/WO/51979/WO-27-51979-g005_min.jpg

We used the R software package to perform GO and KEGG enrichment analyses of the overlapping genes (Suppl. 3) among the KNSTRN-related genes and DEGs. With adjusted p < 0.1, 298 biological processes, 59 cellular components (CC), 23 molecular functions, and 10 KEGG genes were identified. The bubble chart shows the first 3 GO items and all of the KEGG results. Gene ontology functional annotation showed that the overlapping genes were mainly involved in organelle fission, chromosomal regions, and tubulin binding (Fig. 5 D). Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the overlapping genes were mainly involved in the cell cycle, oocyte meiosis, and cellular senescence signalling pathways (Fig. 5 E). Supplementary 4 summarizes the GO and KEGG enrichment analysis results of the overlapping genes. To further understand the internal mechanism of KNSTRN, the STRING database was used to construct the KNSTRN PPI network. The analysis showed that KNSTRN was associated with SPAG5. Centromere protein E (CENPE), cyclin B1 (CCNB1), IQ motif containing GTPase-activating protein 1 (IQGAP1), cell division cycle 20 (CDC20), component of NDC80 kinetochore complex (NUF2), mitotic checkpoint serine/threonine kinase, cyclin B2 (CCNB2), charged multivesicular body protein 1B (CHMP1B), and aurora kinase B were 0.97, 0.922, 0.844, 0.838, 0.834, 0.833, 0.825, 0.802, 0.775, and 0.774, respectively (Fig. 5 F, Suppl. 5). In general, these related genes are highly expressed in several human cancers, and they promote tumourigenesis.

KNSTRN expression is associated with immune signatures in liver hepatocellular carcinoma

The infiltration status of 24 immune cell types in HCC tissues was evaluated using ssGSEA. As shown in Figure 6 A, KNSTRN expression showed a negative correlation with the levels of CD8+ T-cells (r = –0.199, p < 0.001), cytotoxic cells (r = –0.254, p < 0.001), dendritic cells (DCs) (r = –0.288, p < 0.001), mast cells (r = –0.155, p < 0.001), neutrophils (r = –0.288, p < 0.001), pDCs (r = –0.233, p < 0.001), γ δ T-cells (Tgd) (r = –0.183, p = 0.007), natural killer cells (NK) (r = –0.130, p = 0.012), and iDC (r = –0.119, p = 0.022). In contrast, KNSTRN expression was positively correlated with the infiltration levels of Th2-cells (r = 0.585, p < 0.001), T-helper cells (r = 0.287, p < 0.001), and TFH (r = 0.128, p = 0.013).

In addition, we divided the 374 tumour samples into 2 groups based on KNSTRN expression, with 187 samples in the high-expression group and 187 samples in the low-expression group. We used the TIMER database to analyse immune infiltration in LIHC with different KNSTRN expression levels (Figs. 6 B, C). The results showed that the infiltration levels of CD8+ T-cells, cytotoxic cells, DCs, iDC, mast cells, neutrophils, and pDCs were significantly lower in patients with LIHC with high KNSTRN expression than in those with low KNSTRN expression. In contrast, compared with that in the low expression group, infiltration of T-helper cells, TFH, and Th2-cells increased in the high KNSTRN expression group (p < 0.05).

Fig. 6

Correlation investigation of KNSTRN expression, immune infiltration, gene markers in exhausted T-cells in The Cancer Genome Atlas liver hepatocellular carcinoma dataset. The correlation analysis results between the expression levels of KNSTRN and the expression levels of 24 immune cell types (A–C), regulatory T-cells, macrophages-M2, MDSC, and CAF (D). The correlation analysis results between the expression levels of KNSTRN and the expression levels of PDCD1 (E), CD274 (F), CTLA-4 (G), HAVCR2 (H), LAG3 (I), TOX (J), and TP53 (K)

/f/fulltexts/WO/51979/WO-27-51979-g006_min.jpg

The TIMER tool was used to analyse the correlation between KNSTRN expression and the markers of various immune cells in HCC. The results showed that the expression of KNSTRN was significantly correlated with the immune markers CD19, CD79A, and KRT20 in B-cells in HCC (p < 0.001) (Table 2). We also analysed the levels of various T-cells with different functions, including CD8+ T, Tfh, Th1, Th2, Th17, and Treg cells. After adjusting for tumour purity, the results showed that the expression level of KNSTRN was significantly correlated with the immune markers of most T-cells in HCC, including CD8B, CD25, CD183, CD185, CD278, CD360, IL23R, and CD196 (Table 2). This indicates that KNSTRN may be involved in the T-cell immune response in HCC.

Table 2

Correlation analysis between KNSTRN and biomarkers of immune cells in hepatocellular carcinoma

DescriptionGene markersCorp-value
TAMCCL20.1726841150.001281692
CD680.2202950473.66E-05
IL100.2993830551.42E-08
MonocyteCD860.3690728411.42E-12
CSF1R(CD115)0.2197700683.83E-05
M1-macrophageNOS20.0646329580.231150943
M2-macrophageCD1630.1701700450.001511499
VSIG40.20005830.000183791
CSF1R(CD115)0.2197700683.83E-05
DCsITGAX(CD11C)0.4121176491.41E-15
CD1C0.1569910910.003460685
NRP10.2664517195.11E-07
NeutrophilsCCR70.1751685760.001086553
ITGAM(CD11b)0.3687047291.50E-12
CD590.1851405790.000547686
Th1STAT10.3885299127.07E-14
STAT40.3130038732.80E-09
TBX210.1665130870.00191368
CD40.3201098651.16E-09
IFNG (IFN-g)0.2652762965.76E-07
Th2GATA30.2791148421.36E-07
STAT60.2069071150.000108322
CXCR40.3984498241.41E-14
Th2CCR40.3262690955.34E-10
TregFOXP30.3053464337.05E-09
CCR80.46640584.90E-20
STAT5B0.3126177262.94E-09
TGFB10.2743604612.25E-07
NKsKIR2DL10.0250792440.642498226
KIR2DS40.1157803410.031558788
KIR3DL10.0517556860.33782438
B-cellCD190.2401849966.44E-06
CD79A0.1985385630.000206183
KRT200.2738164542.38E-07
CD8+ T-cellCD8B0.1904221060.00037555
CD25(IL2RA)0.322057979.11E-10
TFHCD183(CXCR3)0.2926718263.06E-08
  CD185(CXCR5)0.2526980151.99E-06
CD278(ICOS)0.3418862476.81E-11
Th17CD360(IL21R)0.3391126919.90E-11
IL23R0.1915274210.000346592
CD196(CCR6)0.2868802715.84E-08

* p < 0.05

** p < 0.01

*** p < 0.001

In addition, we used the TIMER database to investigate whether KNSTRN expression is related to the number of common immunosuppressive cells and CAFs in LIHC. As shown in Figure 6 D, KNSTRN expression was positively correlated with the levels of Tregs (r = 0.37, p < 0.001), macrophages-M2 (r = 0.339, p < 0.001), MDSC (r = 0.576, p < 0.001), and CAFs (r = 0.304, p < 0.001). Furthermore, positive correlations between the expression levels of KNSTRN and gene markers in exhausted T-cells, such as PDCD1 (PD-1, r = 0.282, p < 0.001), CD274(PDL-1, r = 0.207, p < 0.001), CTLA4 (r = 0.271, p < 0.001), HAVCR2 (TIM-3, r = 0.259, p < 0.001), LAG3 (r = 0.233, p < 0.001), and TOX (r = 0.161, p = 0.002) were observed in HCC (Figs. 6 E–J).

NamiRNA−enhancer−KNSTRN activation network construction in liver hepatocellular carcinoma

As shown in Figure 7 A, 798 KNSTRN-associated miRNAs were obtained from the LinkedOmics website. Additionally, 257 differentially expressed miRNAs were identified between the normal and tumour groups from the TCGA database, and 303 enhancer-associated miRNAs were obtained as potential NamiRNAs from a published study [20]. See Supplementary 6 for a detailed description of these miRNAs. Furthermore, correlation analysis revealed that KNSTRN was significantly and positively correlated with hsa-miR-107 (r = 0.119, p < 0.05), as shown in Figure 7 B. Next, we determined the expression and prognostic value of hsa-miR-107 in HCC cells. As presented in Figures 7 C, D, hsa-miR-107 expression was markedly upregulated in HCC, which was negatively correlated with patient prognosis (HR = 0.62, p < 0.05). These findings suggest that hsa-miR-107 may be the most potent regulatory NamiRNA of KNSTRN in HCC. Finally, a potential NamiRNA-enhancer network of hsa-miR-107 that activates KNSTRN expression in LIHC was constructed (Fig. 8).

Fig. 7

Identification of hsa-miR-107 as a potential NamiRNA of KNSTRN in liver hepatocellular carcinoma (LIHC). The overlapping miRNA of the KNSTRN-associated miRNAs, differentially expressed miRNAs, and enhancer-associated miRNAs was regarded as potential NamiRNA (A), correlation analysis of KNSTRN expression with hsa-miR-107 by The Cancer Genome Atlas (TCGA) LIHC database (B), the prognostic value of hsa-miR-107 in the TCGA-LIHC database (C), TCGA database analysis shows the hsa-miR-107 expression levels in 33 types of cancer tissues and their corresponding adjacent normal liver tissues (D)

/f/fulltexts/WO/51979/WO-27-51979-g007_min.jpg
Fig. 8

The model of NamiRNA (hsa-miR-107)-enhancer-KNSTRN network in carcinogenesis of hepatocellular carcinoma

/f/fulltexts/WO/51979/WO-27-51979-g008_min.jpg

Discussion

KNSTRN, also known as C15orf23, is a mitosis-related protein that binds directly to microtubules [17]. However, the role of KNSTRN in tumourigenesis has rarely been investigated. In recent years, studies have reported functional links between KNSTRN and endometrial cancer and LUAD [16, 19]. In this study, we sought to determine the expression and prognostic value of KNSTRN in HCC. We found that KNSTRN expression was upregulated in tumour tissues, and high KNSTRN expression was closely associated with poor OS and PFI in HCC. Furthermore, high KNSTRN expression was associated with worse histological-grade features in HCC cohorts. The area under the curve values for 1-, 3-, and 5-year predicted survival rates were all greater than 0.6. The area under the curve for diagnostic efficiency was 0.935. These findings strongly suggest that KNSTRN is an oncogenic molecule involved in tumour progression and is a potential prognostic and diagnostic biomarker for HCC.

The tumour microenvironment (TME) and tumour immune evasion correlate with cancer prognosis and therapeutic efficacy [28]. Tumour infiltration by immune cells leads to T-cell anergy or dysfunctional T-cell phenotypes [29, 30], which promote tumour escape from the host immune system, and tumour progression, invasion, metastasis, and therapeutic resistance [29]. In this study, given that KNSTRN was abnormally highly expressed in LIHC, we speculated that KNSTRN may be involved in regulating the tumour immune response. Therefore, we explored the possible correlation between the infiltration levels of various immune cells and expression of KNSTRN in LIHC (TCGA data). As shown in Figures 6 A–C, a significant positive correlation was observed between KNSTRN expression and the proportion of TFH, T-helper, and Th2-cells in HCC tissues. These results are consistent with those of previous studies showing that TFH-cells regulate tumour growth and progression via the chemokine receptor CXCR5 [26, 31]. Th2 induces M2 macrophage polarization to promote tumour progression via cytokines [32]. Moreover, Th2 plays an important role in the mechanisms of resistance against cancer immunotherapy [33].

One of the most important features of immunosuppression is T-cell exhaustion. In our study, upregulated KNSTRN expression was strongly associated with several T-cell exhaustion markers (Figs. 6 E–J), including PD-1(PDCD1), PD-L1(CD274), CTLA4, TIM-3(HAVCR2), LAG3, and TOX. PD1/PDL1 is an important immune checkpoint that regulates the function of tumour-infiltrating lymphocytes. Notably, although PD1/PDL1 checkpoint blockade therapy is widely used for various malignancies, including HCC, some studies have reported that PD-1 plays a critical role in tumour antigen tolerance, leading to poor therapeutic effects in some patients who receive PD1 therapy [3436]. Concurrently, combinations of inhibitory receptors, such as PD-1 and TIM-3, regulate exhausted T-cells and cooperate to induce CD8+ T-cell deficiency [37]. Furthermore, the activation of TIM-3 and CTLA-4 suppresses T-cell-mediated immune responses [38]. In addition, Khan et al. reported that TOX+ cells express multiple inhibitory receptors, such as LAG3 and PDCD1, suggesting that TOX is a major regulatory molecule of differentiation of TEX at the transcriptional and epigenetic levels [39]. Additionally, T-cell exclusion is dependent on the infiltration of immunosuppressive cells, including CAFs, Tregs, M2-TAMs, and MDSCs [40, 41]. Notably, as shown in Figure 6 D, we found that KNSTRN expression levels strongly correlated with the number of immunosuppressive cells in HCC. Thus, we speculate that T-cell exclusion is the major mechanism through which KNSTRN regulates tumour escape of immune cells, thereby contributing to HCC growth and progression.

The correlations between KNSTRN and the immune markers of various immune cells in HCC are shown in Table 2. In our study, we found that the expression level of KNSTRN was significantly correlated with immune markers of M2 macrophages in HCC, such as CD163, immunoglobulin domain containing 4, and colony stimulating factor 1 receptor (p < 0.001). This shows that KNSTRN has the potential to regulate macrophage polarization in HCC. Moreover, the expression of KNSTRN was positively correlated with immune markers of neutrophils and DCs in HCC, including ITGAX, CD1C, NRP1, CCR7, ITGAM, and CD59 (p < 0.05) (Table 2). After adjusting for tumour purity, the results revealed that the expression level of KNSTRN was significantly correlated with most of the immune markers of different T-cells in HCC. In summary, these findings indicate that KNSTRN plays a different crucial role in the regulation of the TME.

However, cancer development and therapeutic limitations cannot be attributed to a single gene. Therefore, we constructed a PPI network to screen hub genes involved in KNSTRN regulation. Notably, we found significant correlations (r > 0.8) between KNSTRN expression levels and the expression levels of SPAG5, CENPE, CCNB1, IQGAP1, CDC20, NUF2, BUB1B, and CCNB2 in HCC (Fig. 5 F, Suppl. 5). SPAG5 acts as an oncogene and is a promising prognostic factor in HCC [4245]. Human CENPE are potential biomarkers for clinical diagnosis [46]. Cyclin B1 and CCNB2 are independent risk factors for HCC [47]. Additionally, overexpression of IQGAP1 has been observed in HCC [48]. Zhao et al. and Lai et al. reported that CDC20 is a tumour promoter in HCC [49, 50]. Upregulated expression of BUB1B, CCNB1, and CDC20 in tumour tissues predicts worse OS and disease-free survival in patients with HCC [5153]. NUF2 is highly expressed in several human cancers, and it promotes tumourigenesis [5456]. Collectively, the expression of these related genes suggests that KNSTRN plays a key role in the progression of HCC.

DNA methylation during carcinogenesis affects gene expression and the prognosis of patients with cancer [57]. Notably, our data revealed that the methylation level of KNSTRN-cg08036289 was significantly correlated with the expression of KNSTRN. Moreover, hypermethylation of the CpG site KNSTRN-cg08036289 (p = 0.042) was associated with improved prognosis in HCC.

We performed GO and KEGG functional enrichment analyses to further clarify the underlying biological functions of KNSTRN. Cell cycle disorder is one of the most remarkable features of tumour cells. In this study, KEGG pathway analysis showed that KNSTRN was mainly involved in cell cycle, oocyte meiosis, and cellular senescence signalling pathways. Furthermore, using GO functional enrichment analysis, we found that KNSTRN co-expression was mainly involved in organelle fission, chromosomal regions, and tubulin binding. Thus, GO and KEGG analyses revealed that KNSTRN may play a role in influencing the cell cycle and regulating cellular senescence.

Most existing studies regarding miRNA are based on miRNA binding to the 3’UTR of target mRNA to degrade or repress expression of the corresponding gene [58, 59]. Notably, growing evidence has shown that miRNAs can also exhibit positive roles and promote gene expression under certain circumstances [20, 6062]. Liang et al.[22] defined a subset of NamiRNAs with an activation function in the nucleus. Several studies have proposed a novel regulatory network for NamiRNA-enhancer gene activation [63, 64], which is crucial for cell identification throughout the transformation of cells from normal to cancerous [21]. To ascertain whether KNSTRN expression was positively modulated by NamiRNAs, we first predicted a NamiRNA that could potentially bind to KNSTRN and identified hsa-miR-107. Based on the mechanism of action of NamiRNAs in the regulation of target gene expression, a positive correlation should exist between hsa-miR-107 and KNSTRN expression. In our study, correlation analysis revealed that KNSTRN expression was significantly and positively correlated with hsa-miR-107 (r = 0.119, p < 0.05), as shown in Figure 7 B. Next, we determined the expression and prognostic value of hsa-miR-107 in HCC cells. As presented in Figures 7 C, D, hsa-miR-107 expression was markedly upregulated in HCC, and its upregulation was negatively correlated with patient prognosis (HR = 0.62, p < 0.05). These findings suggest that hsa-miR-107 may be the most potent regulatory NamiRNA of KNSTRN in HCC. Finally, a possible NamiRNA-enhancer network of hsa-miR-107 that activates KNSTRN expression in LIHC cells was determined, as shown in Figure 8. Although the NamiRNA- enhancer-KNSTRN regulatory axis was identified through bioinformatics analysis, further experiments are required to confirm our prediction.

Conclusions

In summary, the current study revealed the diagnostic and prognostic value of KNSTRN in HCC. Simultaneously, KNSTRN expression significantly correlated with a variety of tumour-infiltrating immune cells in patients with HCC. Moreover, we constructed a possible NamiRNA-enhancer network through which hsa-miR-107 activates KNSTRN expression, indicating that KNSTRN may participate in various molecular regulatory mechanisms in HCC. Collectively, KNSTRN could serve as an effective biomarker for diagnosis and prognosis and support the development of new therapeutic strategies for HCC. However, further studies are required to validate these findings.

Notes

[7] Conflicts of interest The authors declare no conflict of interest.

References

1 

Luo Z, Lu L, Tang Q, et al. CircCAMSAP1 promotes hepatocellular carcinoma progression through miR-1294/GRAMD1A pathway. J Cell Mol Med 2021; 25: 3793-802.

2 

Sung H, Ferlay J, Siegel RL, et al. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 2021; 71: 209-249.

3 

Wen N, Cai Y, Li F, et al. The clinical management of hepatocellular carcinoma worldwide: a concise review and comparison of current guidelines: 2022 update. Biosci Trends 2022; 16: 20-30.

4 

Lou W, Liu J, Ding B, et al. Identification of potential miRNA-mRNA regulatory network contributing to pathogenesis of HBV-related HCC. J Transl Med 2019; 17: 7.

5 

Chang TE, Su CW, Huang YS, Huang YH, Hou MC, Wu JC. Hepatitis D virus dual infection increased the risk of hepatocellular carcinoma compared with hepatitis B virus mono infection: a meta-analysis. J Chinese Med Assoc 2022; 85: 30-41.

6 

Yu JH, Cho SG, Jin Y-J, Lee J-W. The best predictive model for hepatocellular carcinoma in patients with chronic hepatitis B infection. Clin Mol Hepatol 2022; 28: 351-361.

7 

Badami E, Busà R, Douradinha B, et al. Hepatocellular carcinoma, hepatitis C virus infection and miRNA involvement: perspectives for new therapeutic approaches. World J Gastroenterol 2022; 28: 2417-2428.

8 

Liu Z, Song C, Suo C, Fan H, Zhang T, Jin L, Chen X. Alcohol consumption and hepatocellular carcinoma: novel insights from a prospective cohort study and nonlinear Mendelian randomization analysis. BMC Med 2022; 20: 413.

9 

Turner PC SA, Diallo MS, Castegnaro JJ, Hall AJ, Wild CP. The role of aflatoxins and hepatitis viruses in the etiopathogenesis of hepatocellular carcinoma-A basis for primary prevention in Guinea-Conakry, West Africa. J Gastroenterol Hepatol 2002; 17: S441-8.

10 

Chen Y, Tian Z. HBV-induced immune imbalance in the development of HCC. Front Immunol 2019; 10: 2048.

11 

Wang H, Lu Z, Zhao X. Tumorigenesis, diagnosis, and therapeutic potential of exosomes in liver cancer. J Hematol Oncol 2019; 12: 133.

12 

Friese A, Faesen AC, Veld PJH, et al. Molecular requirements for the inter-subunit interaction and kinetochore recruitment of SKAP and Astrin. Nat Commun 2016; 7: 11407.

13 

Wong K, van der Weyden L, Schott CR, et al. Cross-species genomic landscape comparison of human mucosal melanoma with canine oral and equine melanoma. Nat Commun 2019; 10: 353.

14 

Kern DM, Monda JK, Su KC, Wilson-Kubalek EM, Cheeseman IM. Astrin-SKAP complex reconstitution reveals its kinetochore interaction with microtubule-bound Ndc80. eLife 2017; 6.

15 

Kern DM, Nicholls PK, Page DC, Cheeseman IM. A mitotic SKAP isoform regulates spindle positioning at astral microtubule plus ends. J Cell Biol 2016; 213: 315-328.

16 

Xiong Y, Ju L, Yuan L, et al. KNSTRN promotes tumorigenesis and gemcitabine resistance by activating AKT in bladder cancer. Oncogene 2021; 40: 1595-1608.

17 

Lee CS, Bhaduri A, Mah A, et al. Recurrent point mutations in the kinetochore gene KNSTRN in cutaneous squamous cell carcinoma. Nat Genet 2014; 46: 1060-1062.

18 

Schmitz L, Grinblat B, Novak B, et al. Somatic mutations in kinetochore gene KNSTRN are associated with basal proliferating actinic keratoses and cutaneous squamous cell carcinoma. J Eur Acad Dermatol Venereol 2019; 33: 1535-1540.

19 

Deng P, Zhou R, Zhang J, Cao L. Increased expression of KNSTRN in lung adenocarcinoma predicts poor prognosis: a bioinformatics analysis based on TCGA data. J Cancer 2021; 12: 3239-3248.

20 

Xiao M, Li J, Li W, et al. MicroRNAs activate gene transcription epigenetically as an enhancer trigger. RNA Biol 2017; 14: 1326-1334.

21 

Liang Y, Xu P, Zou Q, Luo H, Yu W. An epigenetic perspective on tumorigenesis: loss of cell identity, enhancer switching, and NamiRNA network. Semin Cancer Biol 2022; 83: 596-604.

22 

Liang Y, Zou Q, Yu W. Steering against wind: a new network of NamiRNAs and enhancers. Genomics Proteomics Bioinformatics 2017; 15: 331-337.

23 

Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res 2018; 46: D956-D63.

24 

Szklarczyk D, Gable AL, Lyon D, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res 2019; 47: D607-D13.

25 

Bindea G, Mlecnik B, Tosolini M, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity 2013; 39: 782-795.

26 

Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics 2013; 14: 7.

27 

Li T, Fu J, Zeng Z, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res 2020; 48: W509-W14.

28 

Lawal B, Lin LC, Lee JC, et al. Multi-omics data analysis of gene expressions and alterations, cancer-associated fibroblast and immune infiltrations, reveals the onco-immune prognostic relevance of STAT3/CDK2/4/6 in human malignancies. Cancers 2021; 13: 954.

29 

Yu GP, Chiang D, Song SJ, et al. Regulatory T cell dysfunction in subjects with common variable immunodeficiency complicated by autoimmune disease. Clin Immunol 2009; 131: 240-253.

30 

Taube JM, Klein A, Brahmer JR, et al. Association of PD-1, PD-1 ligands, and other features of the tumor immune microenvironment with response to anti-PD-1 therapy. Clin Cancer Res 2014; 20: 5064-5074.

31 

Breitfeld D, Ohl L, Kremmer E, et al. Follicular B helper T cells express CXC chemokine receptor 5, localize to B cell follicles, and support immunoglobulin production. J Exp Med 2000; 192: 1545-1552.

32 

Biswas SK, Mantovani A. Macrophage plasticity and interaction with lymphocyte subsets: cancer as a paradigm. Nat Immunol 2010; 11: 889-896.

33 

Li S, Liu M, Do MH, et al. Cancer immunotherapy via targeted TGF-beta signalling blockade in T(H) cells. Nature 2020; 587: 121-125.

34 

Waidmann O. Recent developments with immunotherapy for hepatocellular carcinoma. Expert Opin Biol Ther 2018; 18: 905-910.

35 

Fu C, Jiang A. Dendritic cells and CD8 T cell immunity in tumor microenvironment. Front Immunol 2018; 9: 3059.

36 

Zhao JJ, Yap DWT, Chan YH, et al. Low programmed death-ligand 1-expressing subgroup outcomes of first-line immune checkpoint inhibitors in gastric or esophageal adenocarcinoma. J Clin Oncol 2022; 40: 392-402.

37 

Liu Z, Ravindranathan R, Kalinski P, Guo ZS, Bartlett DL. Rational combination of oncolytic vaccinia virus and PD-L1 blockade works synergistically to enhance therapeutic efficacy. Nat Commun 2017; 8: 14754.

38 

Speranza MC, Passaro C, Ricklefs F, et al. Preclinical investigation of combined gene-mediated cytotoxic immunotherapy and immune checkpoint blockade in glioblastoma. Neuro Oncol 2018; 20: 225-235.

39 

Khan O, Giles JR, McDonald S, et al. TOX transcriptionally and epigenetically programs CD8+ T cell exhaustion. Nature 2019; 571: 211-218.

40 

Komohara Y, Fujiwara Y, Ohnishi K, Takeya M. Tumor-associated macrophages: potential therapeutic targets for anti-cancer therapy. Adv Drug Deliv Rev 2016; 99: 180-185.

41 

Zhou S, Zhao Z, Zhong H, et al. The role of myeloid-derived suppressor cells in liver cancer. Discov Oncol 2023; 14: 77.

42 

Yang YF, Zhang MF, Tian QH, et al. SPAG5 interacts with CEP55 and exerts oncogenic activities via PI3K/AKT pathway in hepatocellular carcinoma. Mol Cancer 2018; 17: 117.

43 

Zhou H, Wang SC, Ma JM, Yu LQ, Jing JS. Sperm-associated antigen 5 expression is increased in hepatocellular carcinoma and indicates poor prognosis. Med Sci Monit 2018; 24: 6021-6028.

44 

Liu H, Hu J, Wei R, et al. SPAG5 promotes hepatocellular carcinoma progression by downregulating SCARA5 through modifying beta-catenin degradation. J Exp Clin Cancer Res 2018; 37: 229.

45 

Chen W, Chen X, Li S, Ren B. Expression, immune infiltration and clinical significance of SPAG5 in hepatocellular carcinoma: a gene expression-based study. J Gene Med 2020; 22: e3155.

46 

Liu D, Wang Z, Zhong L, et al. Targets and potential mechanism of scutellaria baicalensis in treatment of primary hepatocellular carcinoma based on bioinformatics analysis. J Oncol 2022; 2022: 8762717.

47 

Su Q, Li W, Zhang X, et al. Integrated bioinformatics analysis for the screening of hub genes and therapeutic drugs in Hepatocellular carcinoma. Curr Pharm Biotechnol 2023; 24:1035-1058.

48 

Zoheir KMA, Abdelhafez MA, Darwish AM, Mahrous KF. New approach about the signaling crosstalk between IQGAPs/ NF-kappaB/IL-8 and PDCD5/p53/TRAIL pathways that modulate malignant transformation in hepatocellular carcinoma. Asian Pac J Cancer Prev 2022; 23: 271-279.

49 

Zhao S, Zhang Y, Lu X, et al. CDC20 regulates the cell proliferation and radiosensitivity of P53 mutant HCC cells through the Bcl-2/Bax pathway. Int J Biol Sci 2021; 17: 3608-3621.

50 

Lai X, Wu YK, Hong GQ, et al. A novel gene signature based on CDC20 and FCN3 for prediction of prognosis and immune features in patients with hepatocellular carcinoma. J Immunol Res 2022; 2022: 9117205.

51 

Rong MH, Li JD, Zhong LY, et al. CCNB1 promotes the development of hepatocellular carcinoma by mediating DNA replication in the cell cycle. Exp Biol Med (Maywood) 2022; 247: 395-408.

52 

Zhuang L, Yang Z, Meng Z. Upregulation of BUB1B, CCNB1, CDC7, CDC20, and MCM3 in tumor tissues predicted worse overall survival and disease-free survival in hepatocellular carcinoma patients. Biomed Res Int 2018; 2018: 7897346.

53 

Ando K, Kakeji Y, Kitao H, et al. High expression of BUBR1 is one of the factors for inducing DNA aneuploidy and progression in gastric cancer. Cancer Sci 2010; 101: 639-645.

54 

Kaneko N, Miura K, Gu Z, et al. siRNA-mediated knockdown against CDCA1 and KNTC2, both frequently overexpressed in colorectal and gastric cancers, suppresses cell proliferation and induces apoptosis. Biochem Biophysic Res Commun 2009; 390: 1235-1240.

55 

Sethi G, Pathak HB, Zhang H, et al. An RNA interference lethality screen of the human druggable genome to identify molecular vulnerabilities in epithelial ovarian cancer. PLoS One 2012; 7: e47086.

56 

Hu P, Shangguan J, Zhang L. Downregulation of NUF2 inhibits tumor growth and induces apoptosis by regulating lncRNA AF339813. Int J Clin Exp Pathol 2015; 8: 2638-2648.

57 

Lin SH, Wang J, Saintigny P, et al. Genes suppressed by DNA methylation in non-small cell lung cancer reveal the epigenetics of epithelial-mesenchymal transition. BMC Genomics 2014; 15: 1079.

58 

Krol J, Loedige I, Filipowicz W. The widespread regulation of microRNA biogenesis, function and decay. Nat Rev Genet 2010; 11: 597-610.

59 

Wust S, Drose S, Heidler J, et al. Metabolic maturation during muscle stem cell differentiation is achieved by miR-1/133a-mediated inhibition of the Dlk1-Dio3 mega gene cluster. Cell Metab 2018; 27: 1026-39e6.

60 

Vasudevan S, Tong Y, Steitz JA. Switching from repression to activation: microRNAs can up-regulate translation. Science 2007; 318: 1931-1934.

61 

Buchan JR, Parker R. Molecular biology. The two faces of miRNA. Science 2007; 318: 1877-1878.

62 

Huang V, Place RF, Portnoy V, et al. Upregulation of cyclin B1 by miRNA and its implications in cancer. Nucleic Acids Res 2012; 40: 1695-1707.

63 

Shen H, Xu W, Guo R, et al. Suppression of enhancer overactivation by a RACK7-histone demethylase complex. Cell 2016; 165: 331-342.

64 

Liu S, He X, Di Y, et al. NamiRNA-enhancer network of miR-492 activates the NR2C1-TGF-beta/Smad3 pathway to promote epithelial-mesenchymal transition of pancreatic cancer. Carcinogenesis 2023; 44:153-165.

Copyright: © 2023 Termedia Sp. z o. o. This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International (CC BY-NC-SA 4.0) License (http://creativecommons.org/licenses/by-nc-sa/4.0/), allowing third parties to copy and redistribute the material in any medium or format and to remix, transform, and build upon the material, provided the original work is properly cited and states its license.
 
Quick links
© 2024 Termedia Sp. z o.o.
Developed by Bentus.