Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 10 November 2022
Sec. Gastrointestinal Cancers: Gastric and Esophageal Cancers
This article is part of the Research Topic Biomarkers, Functional Mechanisms, and Therapeutic Potentials in Gastrointestinal Cancers View all 40 articles

SKA1/2/3 is a biomarker of poor prognosis in human hepatocellular carcinoma

Guo-Qiang Song&#x;Guo-Qiang Song1†Tian-Li He&#x;Tian-Li He2†Ke-Jie JiKe-Jie Ji1Yi-Meng DuanYi-Meng Duan2Jia-Wen ZhangJia-Wen Zhang1Guo-Qiang Hu,*Guo-Qiang Hu1,3*
  • 1Department of Respiratory, Changxing Hospital of Traditional Chinese Medicine, Huzhou, China
  • 2Department of Radiotherapy, Changxing People’s Hospital, Huzhou, China
  • 3Department of Cancer Center, Changxing Hospital of Traditional Chinese Medicine, Huzhou, China

Background: Spindle and kinetochore-associated complex subunits 1–3 (SKA1–3) stabilize the kinetochore-attached spindle microtubules in metaphase. Due to the dysregulation in multiple cancers, SKA1–3 is considered a predictor for the prognosis of the patients. However, the potential clinical applications of SKA1–3, particularly in hepatocellular carcinoma (HCC) prognosis and progression, have completely unknown yet.

Methods: For the analysis of SKA1–3 expression and applications in clinics in HCC patients, several databases, such as STRING, UALCAN, GEO, and TCGA, were searched. In addition, the underlying mechanisms of SKA for the regulation of HCC occurrence, development, and progression were also explored.

Results: Compared to the normal controls, HCC patients showed dramatically elevated SKA1–3 expression at the mRNA level, and the values of the area under the curve (AUC) were 0.982, 0.887, and 0.973, respectively. Increased SKA1–3 expression levels were associated with the clinical stage, age, body mass index, tumor grade, tissue subtype, and Tp53 mutation status in HCC patients. The analyses of Kyoto Encyclopedia of Genes and Genome (KEGG) and Gene ontology (GO) demonstrated that SKA1–3 are enriched mainly in the Fanconi anemia, homologous recombination, spliceosome, DNA replication, and cell cycle signaling pathways. The hub genes, such as CDK1, CCNB1, CCNA2, TOP2A, BUB1, AURKB, CCNB2, BUB1B, NCAPG, and KIF11, were identified in protein–protein interactions (PPIs). The expression levels of hub genes were increased in HCC patients and predictive of a poor prognosis. Finally, the expression levels of SKA1–3 were determined using the GEO database.

Conclusions: SKA1–3 are potential prognostic biomarkers of and targets for HCC. In addition, SKA1–3 may affect HCC prognosis via the Fanconi anemia pathway, homologous recombination, spliceosome, DNA replication, and cell cycle signaling pathway.

Introduction

Liver cancer is the most common cancer worldwide and has high morbidity and mortality rates (1). The patients who were diagnosed as liver cancer worldwide increased from 471,000 in 1990 to 1,007,800 in 2016 (1). Hepatocellular carcinoma (HCC) is considered the most frequently identified primary liver cancer. Patients with similar pathological types may have different underlying molecular etiologies, leading to different responses to treatment (2). The commonly used cancer treatments, such as surgical resection, radiation therapy, chemotherapy, and immune checkpoint therapy, have limited efficacy against HCC (36). Therefore, new therapies are required for advanced HCC to improve patient outcomes; in particular, multimodal therapies are required that can improve the tumor response in HCC patients.

Spindle and kinetochore-associated complex subunits (SKAs) stabilize the centromere-attached spindle microtubules during mid-mitosis, which promotes the completion of mitosis. SKAs include SKA1–3 and play important roles in non-neoplastic diseases, such as obsessive-compulsive disorder (4); they also predict the development of post-traumatic stress disorder (5). SKA1–3 affect the occurrence and development of several cancers and is a poor prognosis marker in lung cancer, esophageal carcinoma, and kidney renal papillary cell carcinoma (KIRP) (6).

The SKA1 expression level is increased in non-small cell lung cancer and is associated with cancer progression by regulating cell proliferation, migration, and invasion (7). The SKA2 expression level is upregulated in esophageal carcinoma tissues compared with the adjacent normal tissues. Through enhancing the AKT signaling pathway activity and elevating the expression of certain epithelial–mesenchymal transition-associated markers, such as N-cadherin and Snail, SKA2 could significantly participated in the regulation of esophageal carcinoma cells proliferation, migration, and invasion (8). SKA3 is a target of KIRP and predicts a poor prognosis due to its effects on the RAS/MAPK, PI3K/Akt, hormone estrogen receptor, hormone androgen receptor, DNA damage, and cell cycle (9). Therefore, SKA1–3 are useful prognostic biomarkers and potential therapeutic targets for several cancers.

Previous studies have found that SKA2 accelerates HCC progression via Wnt/β-catenin signaling upregulation (10). However, the potential clinical applications of SKA1–3, particularly for HCC prognosis and development, are yet to be fully elucidated. Thus, we performed a comprehensive analysis using several databases and web tools to investigate SKA expression and its relationships with clinical outcomes and prognosis in HCC patients.

Methods

The analysis using the cancer genome atlas

As a landmark cancer genomics program for the molecular characterization of more than 33 cancer types, 20,000 matched normal samples and primary cancer were contained in TCGA (http://cancergenome.nih.gov/).

The analysis using GEPIA and UALCAN

As two opened databases, GEPIA(http://gepia.cancer-pku.cn/) and UALCAN (http://ualcan.path.uab.edu/) are frequently used to analyze the gene expression of the data from TCGA. For the analysis of SKA1–3 and their hub genes expression in HCC and normal samples, the evaluation of the associations of the SKA1–3 with HCC patients’ clinicopathological characteristics (tumor grade, cancer stage, body mass index, age, and stage), and the determination of the SKA1–3’s prognostic values for HCC prognosis, the GEPIA and UALCAN databases were applied.

The analysis using cBioPortal

To analyze the clinical and genomic data comprehensively, we employed the cBioPortal database (www.cbioportal.org/) to integrates the data sourced from other databases, such as the International Cancer Genome Consortium and TCGA. Meanwhile, we also analyzed and evaluated the SKA1–3 expression at the mRNA level with the cBioPortal database (RNA Seq v2 RSEM).

The analysis of biological functions

The data for gene expression from 535 HCC patients in HTSeq-fragments per kilobase per million (FPKM) were obtained from TCGA website. The Pearson correlation coefficients (p < 0.001, |r| > 0.4) were applied to screen the co-expressed of genes with SKA1–3. For exploring the potential signal pathways and biological functions mediated by SKA1-3, the package “clusterProfiler” was applied to perform the Kyoto Encyclopedia of Genes and Genome (KEGG) and Gene ontology (GO) analyses of co-expressed genes. The categories of the molecular function (MF), cell composition (CC), and biological process (BP) were contained in the analysis of GO. The data of gene expression from TCGA was analyzed by GSEA. The value of P less than 0.05 was considered as significant difference.

Based on the string database, we constructed a PPI network of the gens co-expressed with SKA. The significant difference was set as the combined score more than 0.7. Then, the Cytoscape 3.6.1 software was applied to import the PPI network. The CytoHubba plug-in was used to screen the hub genes, and the first 10 genes were defined as hub genes.

The analysis using GEO

The GEO database (https://www.ncbi.nlm.nih.gov/geo/) was created by the National Center for Biotechnology Information in 2000 and contains the data for high-throughput gene expression. GEO is an international public database that collects and organizes microarray data, next-generation sequencing data, and other forms of high-throughput genomic data uploaded by researchers worldwide, including those related to tumors and non-tumor diseases. We searched the GEO database to obtain the relevant data.

Statistical analysis

The data with normal distribution were showed as means ± standard deviation. The intergroup comparisons were conducted by Student’s t-test. The values at different time points from same individual or group were compared using the paired samples t-test. The data with non-normal distribution were compared by the rank sum test. The P value less than 0.05 was set as the significant difference. All analyses were performed using R 3.6.3(http://www.R-project.org) and R online site(www.xiantao.love/products).

Results

SKA1–3 expression in human cancers

The RNAseq data with level 3 HTSeq-FPKM format was downloaded from the ALL (Pan-Cancer) project of TCGA (http://cancergenome.nih.gov/). Then, the log2-transformation was performed to the RNAseq data in the FPKM format. SKA1–3 expression between the adjacent normal and tumor tissues is shown in Figure 1. Compared with the normal tissues, SKA1 expression was significantly higher in tumor tissues, except in the case of kidney chromophobes. Except for kidney renal clear cell carcinoma, KIRP, and prostate adenocarcinoma, we observed that tumor tissues exhibited significantly high expression of SKA2 compared to that of normal tissues. Additionally, except for pheochromocytoma and paraganglioma, dramatically high expression of SKA3 in tumor tissues was also discovered compared to that in normal tissues. In HCC, SKA1–3 expression levels were all higher in cancer tissues than normal tissues (Figure 2).

FIGURE 1
www.frontiersin.org

Figure 1 The SKA1/2/3 expression levels in different types of human cancers. (A) SKA1, (B) SKA2 (C) SKA3 (ns, p≥0.05 *P<0.05, **P<0.01, ***P<0.001).

FIGURE 2
www.frontiersin.org

Figure 2 The expression of SKA1/2/3 in HCC in TGCA. SKA, spindle and kinetochore-associated; HCC, Hepatocellular carcinoma; the cancer genome atlas, TGCA (ns, p≥0.05, ***P<0.001).

SKA1–3 mRNA expression in UALCAN

Analysis of the cBioPortal and TCGA databases showed increased frequencies of SKA1–3 mRNA expression in HCC of 4%, 11%, and 9%, respectively (Figures 3 and 5), with AUC values of 0.982, 0.887, and 0.973, respectively (p < 0.05; Figure 4).

FIGURE 3
www.frontiersin.org

Figure 3 Expression of SKA1/2/3 in cBioPortal database in HCC patients; SKA, spindle and kinetochore-associated; HCC, Hepatocellular carcinoma.

FIGURE 4
www.frontiersin.org

Figure 4 Expression of SKA1/2/3 in TGCA database. (A) is SKA1, (B) is SKA2, (C) is SKA3 SKA1/2/3, spindle and kinetochore associated complex subunit 1/2/3; AUC, area under curve; Cl, confidence level.

High SKA1–3 levels were closely related to the cancer stage, weight, tumor grade, Tp53 mutation status, tumor histology, and HCC prognosis. The significant correlations of the elevated SKA1–3 expression with the HCC patient’s prognosis and clinicopathological features were observed from the UALCAN database. Additionally, a association of the high SKA1 expression with the tumor histology, Tp53 mutation status, tumor grade, weight, and cancer stage, was also discovered. Meanwhile, we also observed the close relation between the elevated SKA2 expression and the tumor histology, Tp53 mutation status, tumor grade, weight, age, and cancer stage, as well as the correlation between the increased SKA3 expression with the tumor histology, Tp53 mutation status, tumor grade, weight, and cancer stage (p < 0.05; Table 1).

TABLE 1
www.frontiersin.org

Table 1 The relationships between SKA1/2/3 expression and clinicopathological features in HCC patients in the UALCAN database.

In addition, compared to those low SKA1- and SKA3-expressed HCC patients, the SKA1- and SKA3-expressed patients exhibited a worse prognosis (p < 0.01; Figures 5A, C). However, no significant association of the SKA2 expression with the HCC prognosis was observed (p > 0.05; Figure 5B).

FIGURE 5
www.frontiersin.org

Figure 5 The relationship of SKA1 (A), SKA2 (B), SKA3 (C) with OS in HCC patients in TGCA database. SKA, spindle and kinetochore associated; OS, overall survive.

And we testify in TGCA database, we found the higher expression of SKA1-3 in tumor group, pathological stage III and IV group, tumor status and AFP more than 400ng/ml group. But there was no significant difference in gender. The expression of SKA2-3 in age less than 60 years was higher than older than 60 years. The expression of SKA1-2 in BMI less than 25Kg/m2 was higher than heavier than 25Kg/m2. The expression of SKA1 and SKA3 in tumor free group were lower than the group with tumor, as well as the patient with vascular invasion (Table 2).

TABLE 2
www.frontiersin.org

Table 2 The relationships between SKA1/2/3 with clinicopathological parameters in HCC patients.

Genes co-expressed with SKA1–3

In TCGA transcriptome data, the numbers of the gene positively correlated with SKA1, SKA2, and SKA3 were 2,757, 1,381, and 2,815, and the numbers of the gene negatively correlated with SKA1, SKA2, and SKA3 were 89, 12, and 82 respectively. Figures 6A–C showed the top 10 co-expressed genes correlated with SKA1–3 positively and negatively. Moreover, the 1,113 intersections among SKA1–3 co-expressing genes are shown in a Venn diagram (Figure 6D).

FIGURE 6
www.frontiersin.org

Figure 6 The top 10 genes with positive and negative co-expression of SKA1 (A), SKA2 (B), SKA3 (C) in TCGA database according to Heat map and Venn map. (D) Intersection co-expression genes of SKA1/2/3. Note: |r| >0.4 and P<0.001. SKA, spindle and kinetochore associated. ***P<0.001.

GO and KEGG analyses

Next, for exploration of the SKA1–3 functions in HCC, we conducted the analyses of GO and KEGG in the genes co-expressed with SKA1–3. The main BPs affected by the SKA1–3 co-expressed genes included DNA replication, chromosome segregation, nuclear division, and organelle fission. The main CCs in associated with the SKA1–3 co-expressed genes included centromeric region, chromosome, condensed chromosome, spindle, and chromosomal region. The main MFs influenced by the SKA1–3 co-expressed genes included DNA-dependent ATPase activity, DNA helicase activity, helicase activity, acting on DNA, and catalytic activity (Figure 7). The results from KEGG analysis indicated that the SKA1–3 co-expressed genes regulated the Fanconi anemia pathway, homologous recombination, cell cycle signaling pathway, spliceosome, and DNA replication (Figure 8). Our results suggest that in HCCs, SKA1–3 participate in the Fanconi anemia pathway, homologous recombination, spliceosome, DNA replication, and cell cycle signaling pathway.

FIGURE 7
www.frontiersin.org

Figure 7 GO analysis of SKA1 (A), SKA2 (B), SKA3 (C) co-expression genes, and intersection co-expression genes among with SKA1/2/3 (D). GO, Gene ontology.

FIGURE 8
www.frontiersin.org

Figure 8 KEGG pathway enrichment analysis of SKA1 (A), SKA2 (B), SKA3 (C) co-expression genes, and intersection co-expression genes among with SKA1/2/3 (D). SKA, spindle and kinetochore associated complex subunit. KEGG, Kyoto Encyclopedia of Genes and Genome.

Analysis of hub genes in the PPI network

Ten representative hub genes, such as CDK1, CCNB1, CCNA2, TOP2A, BUB1, AURKB, CCNB2, BUB1B, NCAPG, and KIF11, were identified using the analysis of PPI network (Table 3). TCGA analysis demonstrated that the expression levels of CDK1, CCNB1, CCNA2, TOP2A, BUB1, AURKB, CCNB2, BUB1B, NCAPG, and KIF11 were significantly elevated in the patients with HCC (p < 0.05; Figure 9) and associated with the disease-free progression and overall survival (OS) of HCC patients (Figures 10 and 11).

TABLE 3
www.frontiersin.org

Table 3 The correlations between the representative 10 Hub genes with SKA1/2/3 in HCC patients.

FIGURE 9
www.frontiersin.org

Figure 9 Hub gene expression was increased in HCC in TGCA database. The expressions of CDK1, CCNB1, CCAN2, TOP2A, BUB1, AURKB, CCNB2, BUB1B, NCAPG, and KIF11 were shown. ***P<0.001.

FIGURE 10
www.frontiersin.org

Figure 10 The hub gene related to the overall survival (OS) of HCC patients in TGCA database. The expressions of CDK1 (A). CCNB1 (B), CCAN2 (C), TOP2A (D), BUB1 (E), AURKB (F), CCNB2 (G), BUB1B (H), NCAPG (I), and KIF11 (J) were shown.

FIGURE 11
www.frontiersin.org

Figure 11 The hub gene related to the disease specific survival (DSF) of HCC patients in TGCA database. The expressions of CDK1 (A). CCNB1 (B), CCAN2 (C). TOP2A (D), BUB1 (E), AURKB (F), CCNB2 (G), BUBIB (H), NCAPG (I), and KIF11 (J) were shown.

External validation

SKA1–3 mRNA expression levels in the GEO database

For the validation of the results of SKA1–3 expression, we downloaded GSE84402 (5) datasets from GEO (http://www.ncbi.nlm.nih.gov/geo/) and generated the expression profiling arrays using GPL570 (HG-U133_Plus_2) Affymetrix Human Genome U133 Plus 2.0 Array (Affymetrix, Santa Clara, CA, USA). In addition, the GSE84402 dataset was obtained, which included 28 specimens from 14 paired HCC and corresponding non-cancerous tissues. Similar to TCGA analysis, SKA1–3 expression levels were higher in tumor than normal tissues (Figure 12).

FIGURE 12
www.frontiersin.org

Figure 12 The expression of SKA1 (A), 2 (B), 3 (C) in HCC patients in GEO database (GSE 884402). (*P<0.05, **P<0.01, ***P<0.001).

Discussion

As a key component in the motor-binding microtubules (11), SKA aberrantly expression results in spindle checkpoint defects, which play a key role in cell cycle regulation and tumorigenesis (12). Several studies have demonstrated that SKA1–3 are involved in apoptosis and tumorigenesis, and their abnormal expression or activation is common in malignant tumors. As a result, SKAs are significantly correlated with malignant tumors (13).

High SKA1 expression was significantly correlated with tumor size, cellular differentiation, and a poor prognosis in pancreatic ductal adenocarcinoma. In addition, SKA1 in pancreatic ductal adenocarcinoma acts as a tumor promoter, and SKA1 overexpression promoted cell proliferation, migration, and invasion in vitro and in vivo (14). SKA1 knockdown inhibited cell proliferation, invasion, migration, and cell cycle arrest in human adenoid cystic carcinoma (15). SKA1 was found to be significantly overexpressed in esophageal squamous cell carcinoma cells, and suppression of the SKA1 expression could significantly inhibit the cell cycle progress, cell migration and proliferation, and promote the cell apoptosis (16).

A recent study reported that the mRNA level and immunohistochemical staining of SKA2 were significantly increased in HCC tissues compared with normal tissues (17). Due to the direct targeting of miR-140-3p to SKA2, overexpression of SK2 could partially suppress the inhibitory effect of miR-140-3p restoration in breast cancer cells (18). Another study demonstrated that miR-301, the locus of which is in the intron of the SKA2 gene, is responsible for centromere assembly, and that the two genes are co-expressed in primary breast cancer samples (19). The association of SKA2 with the metastasis of breast cancer was verified by the inhibition of the metastasis and migration of breast cancer observed after the blocking the SKA2 expression, and the underlying mechanism of this association is the E-cadherin translocation from the cytoplasm to the nucleus (20).

SKA3 is significantly particapted in the processes of chromosome cohesion maintenance and silencing of spindle checkpoint during mitosis (21). In the female patients with early breast cancer, the close association of the elevated SKA3 expression with a poor prognosis was observed (22). The inhibition of cell proliferation, induction of cell cycle arrest in the G1/S phase, as well as in vitro and in vivo tumorigenesis, were significantly connected to the downregulation of SKA3. In addition, SKA3 downregulation led to decreased cyclin D1 expression and retinoblastoma phosphorylation, and increased p21 level, suggesting that SKA3 mediates HCC cell cycle and progression (23). In lung adenocarcinoma, SKA3 promoted metastasis by binding to EGFR and activating the PI3K/AKT signaling pathway (24). The analyses of GSEA and RNA-Seq showed the significant participation of SKA3 in the regulation of the PI3K/Akt signaling pathway and the progression of cell cycle in cervical cancer. Western blot analysis showed that elevated SKA3 expression is significantly and positively associated with the increased expression of p-Rb, E2F1, CDK4, cyclin D1, CDK2, cyclin E2, and p-Akt in HeLa cells (25). In addition, SKA3 knockdown inhibited epithelial–mesenchymal transition (26).

In summary, SKA1–3 expression is correlated with cell apoptosis, proliferation, and cell cycle progression and connected to cancer patients’ poor prognosis. However, the expression patterns and roles of SKA1–3 in HCC are unclear. In the present study, the expression, correlations with clinicopathological features, and underlying mechanisms of SKAs in HCC were explored.

In the present study, SKA1–3 expression was upregulated in HCC patients, which was associated with clinical stage, age, body mass index, tumor grade, tissue subtype, and Tp53 mutation status. Furthermore, increased SKA1–3 mRNA levels were associated with OS (SKA2: hazard ratio = 1.22; 95% confidence interval = 0.86–1.72; p = 0.263). Therefore, SKA1–3 could be used as novel targets of anticancer therapy. The main BPs associated with the genes co-expressed with SKA1–3 included DNA replication, chromosome segregation, nuclear division, and organelle fission. The main CCs in associated with the SKA1–3 co-expressed genes included centromeric region, chromosome, condensed chromosome, spindle, and chromosomal region. The main MFs influenced by the SKA1–3 co-expressed genes included DNA-dependent ATPase activity, DNA helicase activity, helicase activity, acting on DNA, and catalytic activity. The analysis of KEGG demonstrated that the genes co-expressed with SKA1–3 were involved in regulating the Fanconi anemia pathway, homologous recombination, spliceosome, DNA replication, and cell cycle signaling pathway.

Ten representative hub genes, such as CDK1, CCNB1, CCNA2, TOP2A, BUB1, AURKB, CCNB2, BUB1B, NCAPG, and KIF11, were identified using the analysis of PPI network. The elevated expression of CDK1, CCNB1, CCNA2, TOP2A, BUB1, AURKB, CCNB2, BUB1B, NCAPG, and KIF11 was observed by TCGA analysis in HCC (p < 0.05), and these elevated expressions were significantly connected to the disease-free progression and OS of HCC patients (27). The spindle rotation can be blocked by CDK-1 through the inhibition of the interaction between dynein with LIN-5-ASPM-1 and microtubules at the meiotic spindle poles, while through suppression of the CDK-1, APC could promote the spindle rotation (28). Due to the participation of CCNB1 in the cell cycle of HCC through the regulation of DNA replication, CCNB1 can be used for the diagnosis of early-stage HCC (29). Immunohistochemical analysis showed significantly high level of CCNA2 was observed in the tissues from breast cancer patients compared to that from normal controls (29). By regulating the p-ERK1/2/p-SMAD2/Snail pathway, TOP2A could enhanced the process of epithelial–mesenchymal transition and subsequently promote the metastasis of HCC (27). In addition, the HCC patients with the invasion into vascular, dramatically elevated BUB1B and CCNB1 expressions were observed (p < 0.05 for both) (30). Additionally, in comparison to the normal controls, the sarcoma tissues and cells exhibited obviously high BUB1, BUB1B, and BUB3 expression (31). High AURKB expression was significantly related to poor prognosis in neuroblastoma patients (32). AURKB knockdown promotes apoptosis in carboplatin-resistant cells in vitro (32). Knockdown of CCNB2 mRNA could inhibit cell proliferation, reduce breast cancer cell migration, block the G2/M cell cycle transition, and increase cell apoptosis (33). NCAPG silencing inhibited cell proliferation and induced apoptosis in endometrial cancer cells via inhibiting the Wnt/β-catenin pathway (34). Because of the depletion of KIF11 could suppress the growth of cells and tumors in vitro and in vivo, it’s believed that KIF1I could predict the prognosis of HCC (35)

We examined tissue samples from patients using the GEO database and found that SKA1–3 mRNA levels were higher in HCC tissues than in adjacent tissues. However, we did not perform in vitro or in vivo experiments. The potential application of SKA1-3 for the development and pathogenesis of HCC, as well as the related oncogenic signaling pathways were analyzed in our present study to develop multi-targeted and SKA1–3-targeted therapies.

Our study had some limitations. First, we used only the GEO database to validate our findings. Second, additional HCC samples are required to determine the interactions between SKAs and the mechanisms regulating HCC development and progression. Third, further research is needed to identify SKA functions at the cellular level.

SKA1–3 mRNA levels were dramatically increased in HCC samples and were significantly correlated with the clinical stage, age, body mass index, tumor grade, tissue subtype, Tp53 mutation status, and prognosis of HCC patients. SKA1–3 can be used to predict the prognosis of HCC patients and treat HCC. Our results will provide novel insight underlying the molecular mechanism of HCC and new therapeutic strategies related to SKA regulation. However, the verification of our present findings and development of clinical applications of SKA as promising treatment targets and prognostic biomarkers should be done in further using experimental studies.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author contributions

G-QS write this paper. T-LH polish this paper. K-JJ analysis all data. Y-MD check grammar. J-WZ check the tables and figures. G-QH design this paper. All authors contributed to the article and approved the submitted version.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Liu Z, Jiang Y, Yuan H, Fang Q, Cai N, Suo C, et al. The trends in incidence of primary liver cancer caused by specific etiologies: Results from the global burden of disease study 2016 and implications for liver cancer prevention. J Hepatol (2019) 70:674–83. doi: 10.1016/j.jhep.2018.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Chidambaranathan-Reghupaty S, Fisher PB, Sarkar D. Hepatocellular carcinoma (HCC): Epidemiology, etiology and molecular classification. Adv Cancer Res (2021) 149:1–61. doi: 10.1016/bs.acr.2020.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Makuuchi M, Sano K. The surgical approach to HCC: our progress and results in Japan. Liver Transplant Off Publ Am Assoc Study Liver Dis Int Liver Transplant Soc (2004) 10:S46–52. doi: 10.1002/lt.20044

CrossRef Full Text | Google Scholar

4. Lisoway AJ, Zai CC, Zai G, Nair A, Ebrahimi S, Tiwari AK, et al. PT613. the role of SKA2 genetic and epigenetic variation in obsessive-compulsive disorder. Int J Neuropsychopharmacol (2016) 19:25. doi: 10.1093/ijnp/pyw044.613

CrossRef Full Text | Google Scholar

5. Boks MP, Rutten BPF, Geuze E, Houtepen LC, Vermetten E, Kaminsky Z, et al. SKA2 methylation is involved in cortisol stress reactivity and predicts the development of post-traumatic stress disorder (PTSD) after military deployment. Neuropsychopharmacol Off Publ Am Coll Neuropsychopharmacol (2016) 41:1350–6. doi: 10.1038/npp.2015.286

CrossRef Full Text | Google Scholar

6. Yu S, Ma J. Spindle and kinetochore-associated complex is associated with poor prognosis in adrenocortical carcinoma. J Surg Res (2022) 277:50–9. doi: 10.1016/j.jss.2022.03.022

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Shen L, Yang M, Lin Q, Zhang Z, Miao C, Zhu B. SKA1 regulates the metastasis and cisplatin resistance of non-small cell lung cancer. Oncol Rep (2016) 35:2561–8. doi: 10.3892/or.2016.4670

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Chen J, Yang H-M, Zhou H-C, Peng R-R, Niu Z-X, Kang C-Y. PRR11 and SKA2 promote the proliferation, migration and invasion of esophageal carcinoma cells. Oncol Lett (2020) 20:639–46. doi: 10.3892/ol.2020.11615

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Feng D, Zhang F, Liu L, Xiong Q, Xu H, Wei W, et al. SKA3 serves as a biomarker for poor prognosis in kidney renal papillary cell carcinoma. Int J Gen Med (2021) 14:8591–602. doi: 10.2147/IJGM.S336799

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Jiang J, Xu B, Zheng Y, Guo X, Chen F. Spindle and kinetochore-associated protein 2 facilitates the proliferation and invasion of hepatocellular carcinoma via the regulation of wnt/β-catenin signaling. Exp Cell Res (2020) 395:112181. doi: 10.1016/j.yexcr.2020.112181

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Hanisch A, Silljé HHW, Nigg EA. Timely anaphase onset requires a novel spindle and kinetochore complex comprising Ska1 and Ska2. EMBO J (2006) 25:5504–15. doi: 10.1038/sj.emboj.7601426

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Monda JK, Whitney IP, Tarasovetc EV, Wilson-Kubalek E, Milligan RA, Grishchuk EL, et al. Microtubule tip tracking by the spindle and kinetochore protein Ska1 requires diverse tubulin-interacting surfaces. Curr Biol CB (2017) 27:3666–3675.e6. doi: 10.1016/j.cub.2017.10.018

CrossRef Full Text | Google Scholar

13. Chen C, Guo Q, Song Y, Xu G, Liu L. SKA1/2/3 serves as a biomarker for poor prognosis in human lung adenocarcinoma. Transl Lung Cancer Res (2020) 9:20. doi: 10.21037/tlcr.2020.01.20

CrossRef Full Text | Google Scholar

14. Li T, Liu X, Xu B, Wu W, Zang Y, Li J, et al. SKA1 regulates actin cytoskeleton remodelling via activating Cdc42 and influences the migration of pancreatic ductal adenocarcinoma cells. Cell Prolif (2020) 53:e12799. doi: 10.1111/cpr.12799

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Zhao L-J, Yang H-L, Li K-Y, Gao Y-H, Dong K, Liu Z-H, et al. Knockdown of SKA1 gene inhibits cell proliferation and metastasis in human adenoid cystic carcinoma. BioMed Pharmacother Biomed Pharmacother (2017) 90:8–14. doi: 10.1016/j.biopha.2017.03.029

CrossRef Full Text | Google Scholar

16. Hu D, Li Z, Li X, Fu H, Zhang M. SKA1 overexpression is associated with the prognosis of esophageal squamous cell carcinoma and regulates cell proliferation and migration. Int J Mol Med (2019) 44:1971–8. doi: 10.3892/ijmm.2019.4343

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wang D, Suo YJ, Gong L, Lv SQ. SKA2 promotes proliferation and invasion of hepatocellular carcinoma cells via activating the β-catenin signaling pathway. Neoplasma (2020) 67:743–50. doi: 10.4149/neo_2020_190709N609

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Dou D, Ren X, Han M, Xu X, Ge X, Gu Y, et al. Circ_0008039 supports breast cancer cell proliferation, migration, invasion, and glycolysis by regulating the miR-140-3p/SKA2 axis. Mol Oncol (2021) 15:697–709. doi: 10.1002/1878-0261.12862

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Shi W, Gerster K, Alajez NM, Tsang J, Waldron L, Pintilie M, et al. MicroRNA-301 mediates proliferation and invasion in human breast cancer. Cancer Res (2011) 71:2926–37. doi: 10.1158/0008-5472.CAN-10-3369

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Ren Z, Yang T, Zhang P, Liu K, Liu W, Wang P. SKA2 mediates invasion and metastasis in human breast cancer via EMT. Mol Med Rep (2019) 19:515–23. doi: 10.3892/mmr.2018.9623

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Daum JR, Wren JD, Daniel JJ, Sivakumar S, McAvoy JN, Potapova TA, et al. Ska3 is required for spindle checkpoint silencing and the maintenance of chromosome cohesion in mitosis. Curr Biol CB (2009) 19:1467–72. doi: 10.1016/j.cub.2009.07.017

CrossRef Full Text | Google Scholar

22. Zhong Y, Zhuang Z, Mo P, Lin M, Gong J, Huang J, et al. Overexpression of SKA3 correlates with poor prognosis in female early breast cancer. PeerJ (2021) 9:e12506. doi: 10.7717/peerj.12506

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Tang J, Liu J, Li J, Liang Z, Zeng K, Li H, et al. Upregulation of SKA3 enhances cell proliferation and correlates with poor prognosis in hepatocellular carcinoma. Oncol Rep (2021) 45:48. doi: 10.3892/or.2021.7999

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Hu D-D, Chen H-L, Lou L-M, Zhang H, Yang G-L. SKA3 promotes lung adenocarcinoma metastasis through the EGFR-PI3K-Akt axis. Biosci Rep (2020) 40:BSR20194335. doi: 10.1042/BSR20194335

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Hu R, Wang M-Q, Niu W-B, Wang Y-J, Liu Y-Y, Liu L-Y, et al. SKA3 promotes cell proliferation and migration in cervical cancer by activating the PI3K/Akt signaling pathway. Cancer Cell Int (2018) 18:183. doi: 10.1186/s12935-018-0670-4

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Wu Y, Dai F, Zhang Y, Zheng X, Li L, Zhang Y, et al. miR-1207-5p suppresses laryngeal squamous cell carcinoma progression by downregulating SKA3 and inhibiting epithelial-mesenchymal transition. Mol Ther Oncolytics (2021) 22:152–65. doi: 10.1016/j.omto.2021.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Dong Y, Sun X, Zhang K, He X, Zhang Q, Song H, et al. Type IIA topoisomerase (TOP2A) triggers epithelial-mesenchymal transition and facilitates HCC progression by regulating snail expression. Bioengineered (2021) 12:12967–79. doi: 10.1080/21655979.2021.2012069

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Ellefson ML, McNally FJ. CDK-1 inhibits meiotic spindle shortening and dynein-dependent spindle rotation in c. elegans. J Cell Biol (2011) 193:1229–44. doi: 10.1083/jcb.201104008

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Rong M-H, Li J-D, Zhong L-Y, Huang Y-Z, Chen J, Xie L-Y, et al. CCNB1 promotes the development of hepatocellular carcinoma by mediating DNA replication in the cell cycle. Exp Biol Med Maywood NJ (2022) 247:395–408. doi: 10.1177/15353702211049149

CrossRef Full Text | Google Scholar

30. 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. doi: 10.1155/2018/7897346

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Long Z, Wu T, Tian Q, Carlson LA, Wang W, Wu G. Expression and prognosis analyses of BUB1, BUB1B and BUB3 in human sarcoma. Aging (2021) 13:12395–409. doi: 10.18632/aging.202944

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Yang Y, Sheng Y, Sun D, Sun J, Li L, Sun L. AURKB promotes tumorigenesis and carboplatin resistance by regulating the ERK pathway in neuroblastoma cells. Int J Neurosci (2021), 1–11. doi: 10.1080/00207454.2021.1914610

CrossRef Full Text | Google Scholar

33. Aljohani AI, Toss MS, El-Sharawy KA, Mirza S, Ball GR, Green AR, et al. Upregulation of cyclin B2 (CCNB2) in breast cancer contributes to the development of lymphovascular invasion. Am J Cancer Res (2022) 12:469–89.

PubMed Abstract | Google Scholar

34. Liu C, Yan Y, Di F, Li W, Yin X, Dong L. Inhibition of NCAPG expression inactivates the wnt/β-catenin signal to suppresses endometrial cancer cell growth in vitro. Environ Toxicol (2021) 36:2512–20. doi: 10.1002/tox.23364

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Hu Z-D, Jiang Y, Sun H-M, Wang J-W, Zhai L-L, Yin Z-Q, et al. KIF11 promotes proliferation of hepatocellular carcinoma among patients with liver cancers. BioMed Res Int (2021) 2021:2676745. doi: 10.1155/2021/2676745

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: spindle and kinetochore-associated complex subunit (SKA), liver hepatocellular carcinoma, enrichment analysis, biomarker, bioinformatics analysis

Citation: Song G-Q, He T-L, Ji K-J, Duan Y-M, Zhang J-W and Hu G-Q (2022) SKA1/2/3 is a biomarker of poor prognosis in human hepatocellular carcinoma. Front. Oncol. 12:1038925. doi: 10.3389/fonc.2022.1038925

Received: 07 September 2022; Accepted: 25 October 2022;
Published: 10 November 2022.

Edited by:

Zequn Li, The Affiliated Hospital of Qingdao University, China

Reviewed by:

Jiang Chen, Zhejiang University, China
Xiaochen Wang, University of Texas Southwestern Medical Center, United States
Zhigang Ren, First Affiliated Hospital of Zhengzhou University, China

Copyright © 2022 Song, He, Ji, Duan, Zhang and Hu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Guo-Qiang Hu, Y2hhbmd4aW5naGdxQDE2My5jb20=

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.