Skip to main content

ORIGINAL RESEARCH article

Front. Surg., 23 September 2022
Sec. Surgical Oncology

CCDC25 may be a potential diagnostic and prognostic marker of hepatocellular carcinoma: Results from microarray analysis

Updated
\r\nHongyang Deng,&#x;Hongyang Deng1,†Jiaxing Zhang,&#x;Jiaxing Zhang2,†Yijun ZhengYijun Zheng1Jipin LiJipin Li1Qi XiaoQi Xiao1Fengxian WeiFengxian Wei1Wei HanWei Han1Xiaodong XuXiaodong Xu1Youcheng Zhang
\r\nYoucheng Zhang1*
  • 1Department of General Surgery, Hepatic-biliary-pancreatic Institute, Lanzhou University Second Hospital, Lanzhou, China
  • 2Key Laboratory of the Digestive System Tumors of Gansu Province, Lanzhou University Second Hospital, Lanzhou, China

Background: Hepatocellular carcinoma (HCC) is a tumor with a high recurrence rate, poor prognosis, and rapid progression. Therefore, it is necessary to find a novel biomarker for HCC. Coiled-coil domain containing 25 (CCDC25) has been identified as a target molecule that mediates liver metastasis in colon cancer. However, the molecular mechanisms of CCDC25 in HCC are unknown. This study aimed to explore the role of CCDC25 in HCC.

Methods: The expression of CCDC25 in HCC was identified through The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Receiver operating characteristic curve (ROC) curves were drawn to evaluate the diagnostic value of CCDC25 for HCC. The effect of CCDC25 on the prognosis of HCC was analyzed by using the Kaplan–Meier plotter. Co-expressed genes and Gene Set Enrichment Analysis (GSEA) were used to explore the related functions and regulatory signaling pathways of CCDC25. Moreover, we employed the Tumor Immune Estimation Resource (TIMER) database and CIBERSORT algorithm to investigate the relationship between CCDC25 and the tumor immune microenvironment (TME) in HCC. Meanwhile, the effect of CCDC25 on the sensitivity of HCC patients to chemotherapy drugs was evaluated. Finally, we explored the prognostic methylation sites of CCDC25 using the MethSurv database.

Results: CCDC25 expression was low in HCC. Low CCDC25 expression was significantly associated with poor overall survival of HCC and may be comparable to the ability of AFP to diagnose HCC. Dysregulation of glucose metabolism, fatty acid metabolism, amino acid metabolism, ubiquitination modification, and apoptosis inhibition caused by CCDC25 downregulation may be the causes and results of HCC. In addition, CCDC25 was positively correlated with the infiltration level of various adaptive antitumor immune cells. The levels of immune cell infiltration and immune checkpoint expression were lower in the samples with high CCDC25 expression. What is more, we found that downregulated CCDC25 may increase the sensitivity or resistance of HCC patients to multiple drugs, including sorafenib. We also identified a methylation site for CCDC25, which may be responsible for poor prognosis and low CCDC25 expression in HCC patients. Finally, CCDC25 may be associated with HCC ferroptosis.

Conclusions: CCDC25 may be a potential diagnostic and prognostic marker for HCC and is associated with immune infiltration and ferroptosis.

Introduction

As the sixth most common cancer, hepatocellular carcinoma (HCC) ranks fourth among cancer deaths (1). According to statistics, 90,600 new HCC patients were diagnosed in 2020, resulting in at least 83,000 deaths (2). At present, the incidence of HCC is still on the rise worldwide (2). On the surface, the 5-year survival rate after local resection or liver transplantation for early HCC patients is approximately 70% (3). However, studies have confirmed that the recurrence rate of HCC after hepatectomy is as high as 50%, and the prognosis of advanced HCC patients is poor (4). What is more, the complex molecular mechanisms of HCC are not yet understood. Therefore, it is necessary to explore new HCC biomarkers to improve the treatment outcome and long-term prognosis of HCC patients.

Coiled-coil domain containing 25 (CCDC25) is located on chromosome 8p, and a group of genes, including CCDC25, is deleted on chromosome 8P in patients with poor HCC prognosis (5). It has been found that CCDC25 overexpression in cholangiocarcinoma (CCA) can enhance the migration ability of CCA cells (6). Later, another study (7) found that the level of CCDC25 in the blood of CCA patients was significantly increased. In addition, CCDC25 is a target molecule (810) for the neutrophil extracellular traps (NETs) associated with cancer metastasis. However, the role of CCDC25 in HCC is still lacking.

In this study, the expression of CCDC25 in HCC and its related functional mechanism were explored using multiple online databases and R language, hoping to provide a theoretical basis for the diagnosis and treatment of HCC.

Materials and methods

Expression analysis of CCDC25 in pan cancer

We used the Tumor Immune Estimation Resource (TIMER) database (11) (https://cistrome.shinyapps.io/timer/) and the HCC (HCCDB) database (12) (http://lifeome.net/database/hccdb/home.html) to detect the CCDC25 expression in pan cancer. The gene expression profile in the TIMER database is based on The Cancer Genome Atlas (TCGA) database, and the expression levels are calculated in log2TPM format. HCCDB (12) is a database integrating 15 HCC public data sets, including Gene Expression Omnibus (GEO), the Liver Hepatocellular Carcinoma Project of TCGA (TCGA-LIHC), and the Liver Cancer-RIKEN JP Project from the International Cancer Genome Consortium (ICGC LIRI-JP).

Data acquisition and preprocessing

The TCGA-LIHC RNA sequencing data (RNA-seq), consisting of 371 HCC samples and 53 adjacent samples, were downloaded in TPM format from the TCGA database (https://portal.gdc.cancer.gov/), and the batch effects were removed using limma R packages. In addition, we used GSE14520, GSE22058, GSE25097, GSE36376, GSE54236, GSE63898, and GSE64041, which were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/), to verify the CCDC25 expression level in HCC. The probe ID of each gene was converted to the gene symbol, and the average expression of the same gene symbol was taken. The count format expression profile was converted to TPM format, and log2(TPM+1) standardization was performed.

Clinical correlation analysis

Clinical characteristics and survival information corresponding to the TCGA-LIHC data set were downloaded from the Xena (https://xenabrowser.net/datapages/) database. A total of 365 HCC patients had RNA expression profiles and clinical information. The clinical characteristics included age, sex, race, TNM stage, clinical stage, pathological grade, survival status, and duration of survival.

Kaplan–Meier plotter and the diagnostic ROC

The Kaplan–Meier plotter (13, 14) (https://kmplot.com/analysis/) database, the data of which emanate from GEO, EGA (European Genome-Phenome Archive), and TCGA, was applied to analyze the relationship between CCDC25 expression and the prognosis of HCC. Log-rank p-values <0.05 represented statistical significance. The diagnostic value of CCDC25 was evaluated through diagnostic receiver operating characteristic curve (ROC) of TCGA-LIHC data and six GEO data sets (GSE101685, GSE121248, GSE36376, GSE59259, GSE62232, and GSE76427). We calculated AUC values and plotted ROC curves using the pROC package in R-4.1.2. The diagnostic values of alpha fetoprotein (AFP) and CCDC25 for HCC were compared by using the Delong test.

Analysis of CCDC25 co-expressed genes and their functional enrichment in HCC

LinkedOmics (15) (http://www.linkedomics.org/login.php) is a public database that contains multiomics data from all 32 TCGA cancer types and 10 Clinical Proteomics Tumor Analysis Consortium (CPTAC) cancer cohorts. The Pearson correlation coefficient was calculated through this database to analyze CCDC25 co-expressed genes, and the results were visualized by using a volcano map and heat map. Genes with a co-expression coefficient greater than 0.2 were analyzed for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment.

Gene set enrichment analysis

Gene Set Enrichment Analysis (16) (GSEA) is a computational method that can be used to determine the effect of co-variation of a set of predefined genes on phenotypic changes. The level of CCDC25 expression was recognized as a phenotypic label. The gene set was permuted 1,000 times for each analysis. KEGG pathways with p < 0.05 and false discovery rate (FDR) < 0.25 were considered to be significantly enriched.

Immune cell infiltration and immune checkpoints

CIBERSORT (17) is an analysis tool for estimating the infiltration degree of 22 kinds of human immune cells by using gene transcriptome data. The TIMER database and CIBERSORT algorithm were used to predict the correlation between CCDC25 expression and tumor-infiltrating immune cells. TISIDB (17) (http://cis.hku.hk/TISIDB/index.php), integrating data from various databases such as Pubmed, TCGA, and DrugBank, is a web portal that analyzes the relationship between tumor and immune system. So, the relationships between CCDC25 expression and the immune checkpoints in HCC were analyzed through TISIDB.

Drug sensitivity and the methylation site analysis

To investigate the relationship between the sensitivity of HCC patients to chemotherapy drugs and CCDC25 expression, a semi-inhibitory concentration (IC50) of commonly used drugs in the LIHC-STAD cohort was predicted using the “pRRophetic” package. MethSurv (18) (https://biit.cs.ut.ee/methsurv/) is an online database containing information on 25 different cancer types and 7,358 patients in TCGA that can be used for multivariate survival analyses using DNA methylation data. Therefore, MethSurv was used to explore CCDC25 DNA methylation sites associated with survival.

Statistical method

Data from TCGA and GEO were collated by R-4.1.2. OS, PFS, and DSS were calculated by using the log-rank test. Correlation analysis was performed by using Pearson's correlation test. The Wilcox test was used to compare the differences of continuous variables between the two groups. A value of p < 0.05 was considered statistically significant.

Result

The expression of CCDC25 in pan cancer and HCC

An analysis of CCDC25 expression in pan cancer based on the TIMER database showed that CCDC25 expression levels in bladder urothelial carcinoma (BLCA), breast invasive carcinoma (BRCA), colon adenocarcinoma (COAD), kidney chromophobe (KICH), kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), prostate adenocarcinoma (PRAD), rectum adenocarcinoma (READ), thyroid carcinoma (THCA), and uterine corpus endometrial carcinoma (UCEC) were significantly lower than those in normal tissues. As for head and neck squamous cell carcinoma (HNSC), CCDC25 expression in samples with HPV+ was significantly higher than that in samples with HPV. In addition, the skin cutaneous melanoma (SKCM) with metastasis had a higher CCDC25 expression compared with the SKCM group (Figure 1A). Interestingly, CCDC25 was also significantly underexpressed in HCC in the HCCDB pan-cancer analysis (Figure 1B). This result was confirmed in GSE14520, GSE22058, GSE25097, GSE36376, GSE54236, GSE63898, GSE64041, the TCGA-LIHC data set, and the TCGA-LIHC paired data set (Figure 2). The above results suggested that CCDC25 may be involved in the regulation of a variety of cancers. Particularly in HCC, downregulated CCDC25 may promote the progression of HCC.

FIGURE 1
www.frontiersin.org

Figure 1. The expression of CCDC25 in pan cancer. (A) TIMER, (B) HCCDB. Red: HCC tissues, Blue: adjacent tissues. (*p < 0.05; **p < 0.01; ***p < 0.001; ns, not significant, p > 0.05.)

FIGURE 2
www.frontiersin.org

Figure 2. CCDC25 expression in different HCC data sets. (A) GSE14520, (B)GSE22058, (C) GSE25097, (D) GSE36376, (E) GSE54236, (F) GSE63898, (G) GSE64041, (H) ALL TCGA-LIHC samples, and (I) TCGA-LIHC paired samples.

Correlation between CCDC25 expression and clinical characteristics in HCC

As we can see in Figures 3A–H, downregulated CCDC25 was associated with patients in age (Age < 60 vs. Age ≥ 60, p = 0.030), gender (female vs. male, p = 0.046), and M stage (M0 vs. M1, p = 0.013).

FIGURE 3
www.frontiersin.org

Figure 3. The difference of CCDC25 expression in different clinical subgroups: (A) age, (B) gender, and (C) M stage.

Diagnostic and prognostic values of CCDC25 in HHC patients

Figures 4A–C show that a lower CCDC25 expression was significantly related to poor overall survival (OS) (p = 0.013), progression-free survival (PFS) (p = 0.014), and disease-specific survival (DSS) (p = 0.013). As shown in Figure 4D, CCDC25 expression had a strong diagnostic value for HCC (the AUC values of TCGA-LIHC, GSE101685, GSE121248, GSE36376, GSE59259, GSE62232, and GSE76427 were 0.732, 0.823, 0.798, 0.756, 0.875, 0.895, and 0.796). In addition, the diagnostic value of CCDC25 was similar to that of AFP (AUC: 0.723 vs. 0.706, p > 0.05) (Figure 4E), and CCDC25 was significantly negatively correlated with AFP (Figure 4F). These results suggested that CCDC25 expression may be a potential predictor of prognosis and diagnosis in HCC.

FIGURE 4
www.frontiersin.org

Figure 4. Diagnostic and prognostic values of CCDC25 in HHC patients. (A) OS, (B) PFS, (C) DDS, and (D) the ROC curve of CCDC25 diagnostic value for HCC in the TCGA-LIHC and GEO data sets. (E) The ROC curve of CCDC25 and AFP diagnostic value of HCC in TCGA-LIHC. (F) The correlation between CCDC25 and AFP in TCGA-LIHC.

Analysis of the CCDC25 co-expressed gene network

With a correlation coefficient greater than 0.2, we identified 1,557 genes that were significantly positively correlated with CCDC25. The heat map shows the top 50 genes with the most significant positive and negative associations with CCDC25 (Figures 5A–C). GO and KEGG enrichment analysis of 1,557 positively related genes showed that CCDC25 expression was mainly involved in metabolism-related functions such as the citrate cycle, fatty acid metabolism, amino acid metabolism, and PPAR signaling pathway (Figures 5D,E). These results suggested that CCDC25 was mainly involved in liver energy metabolism, and its downregulated expression may cause metabolic disorders, thereby promoting the occurrence of HCC.

FIGURE 5
www.frontiersin.org

Figure 5. Function enrichment analysis of CCDC25 co-expressed genes. (A) The volcanic map shows the co-expression genes of CCDC25 in HCC, (B) the top 50 genes positively correlate with CCDC25, (C) the top 50 genes negatively correlate with CCDC25, (D) GO enrichment analysis for CCDC25 co-expressed genes, and (E) KEGG enrichment analysis for CCDC25 co-expressed genes.

Gene set enrichment analysis

GSEA results showed (Figure 6) that CCDC25 was related to the nuclear-transcribed mRNA catabolic process, deadenylation-dependent decay, RNA polyadenylation, helicase activity, ubiquitin-like protein-specific protease activity, damaged DNA binding, adherens junction, apoptosis, lysine degradation, ubiquitin-mediated proteolysis, and valine, leucine, and isoleucine degradation (Figure 6), indicating that CCDC25 is involved in a variety of metabolisms, nucleotide metabolisms, and the promotion of apoptosis processes in the body to maintain normal functions. When CCDC25 expression level is downregulated, abnormalities of these processes may lead to the occurrence of HCC.

FIGURE 6
www.frontiersin.org

Figure 6. Enrichment plots of GSEA. (A) The results of GO enrichment. (B) The top 10 results of KEGG enrichment.

Correlation of CCDC25 expression with immune cell infiltration and immune checkpoints

We comprehensively evaluated the biological role of CCDC25 in TME through TIMER (Figure 7A) and found that the CD8 + T cell, macrophage, neutrophil, and dendritic cells were positively correlated with CCDC25 expression in HCC. The gene expression profile of TCGA-LIHC was then processed by using the CIBERSORT algorithm to explore the difference in 22 immune cell infiltration between the high CCDC25 expression group and the low CCDC25 expression group in HCC. The results show that the naïve B cell, M2 macrophage, activated mast cell, and neutrophil were elevated in the high CCDC25 expression cohort compared with the low one. Contrastively, the proportion of memory B cells and regulatory T cells (Treg) in the group with high CCDC25 expression was significantly lower than that in the group with low CCDC25 expression (Figure 7B). In addition, CCDC25 expression was found to be negatively correlated with the immune checkpoints PDCD1, CTLA4, TIGIT, and TIM3 (HAVCR2) by TISIDB (Figure 8). These findings suggested that CCDC25 expression may promote antitumor immunity by upregulating CD8 + T cells and dendritic cells and downregulating Treg cells, and it may inhibit tumor immune escape by inhibiting immune checkpoint expression.

FIGURE 7
www.frontiersin.org

Figure 7. CCDC25 affects immune cell infiltration in HCC patients. (A) CCDC25 expression is significantly correlated with the TIMER immune cell infiltration level. (B) The boxplot shows the relationship between CCDC25 expression and the abundance of 22 kinds of immune cells in HCC.

FIGURE 8
www.frontiersin.org

Figure 8. The relationship between CCDC25 expression and the immune checkpoints of HCC. (A) CTLA4, (B) HAVCR2, (C) PDCD1, and (D) TIGIT.

Drug sensitivity analysis

For patients with unresectable HCC, chemotherapy and targeted therapy may improve patient outcomes (19). Sorafenib is currently the first-line chemotherapy for HCC patients, but the tolerability of sorafenib shown by many patients remains a major challenge (20). To identify the efficacy of CCDC25 as a biomarker for predicting treatment response in HCC patients, we evaluated the IC50 of 138 drugs in TCGA-LIHC patients. We discovered that patients with low CCDC25 expression may be more sensitive to Sorafenib, Lapatinib, Gefitinib, and Selumetinib, while patients with high CCDC25 expression may respond better to Ponatinib, Axitinib, Bleomycin, Bosutinib, and Methotrexate (Figure 9). Taken together, these results suggested that CCDC25 was associated with a drug sensitivity of HCC, and patients with low CCDC25 expression may have high sensitivity to Sorafenib and may be expected to help in the sensitization treatment of Sorafenib.

FIGURE 9
www.frontiersin.org

Figure 9. Relationship between chemotherapeutic drug sensitivity and CCDC25 expression in HCC. (A) Sorafenib, (B) Lapatinib, (C) Gefitinib, (D) Selumetinib, (E) Ponatinib, (F) Axitinib, (G) Bleomycin, (H) Bosutinib, and (I) Methotrexate.

The methylation site of CCDC25

We identified eight DNA methylation sites of CCDC25. Among them, CCDC25-TSS200-Island-cg16735490 was significantly positively related to a shorter survival time (HR = 1.679, p = 0.013) (Figure 10). These results supported previous findings that HCC patients with low CCDC25 expression had a poorer prognosis, and the low CCDC25 expression may be associated with hypermethylation at site CCDC25-TSS200-Island-cg16735490.

FIGURE 10
www.frontiersin.org

Figure 10. DNA methylation sites of CCDC25. (A) Blue to red indicates low to high levels of methylation. The different colored boxes represent ethnicity, race, Body Mass Index (BMI), age, event, relation to UCSC_CpG_Island, UCSC_RefGene_Group. (B) The HCC survival analysis of CCDC25-TSS200-Island-cg16735490.

CCDC25 expression may be associated with ferroptosis in HCC

Apoptosis is a classical form of cell death and plays an important role in maintaining metabolism. In the above analysis, CCDC25 was found to be involved in apoptosis. Ferroptosis, as a new type of regulatory cell death, has been reported to have negative effects on cancer (21). Therefore, to investigate the role of CCDC25 in the ferroptosis of HCC cells, we analyzed the correlation between CCDC25 and the expression levels of various key enzymes of ferroptosis through the TIMER database. The results suggested that CCDC25 was significantly negatively correlated with GPX4 (cor = −0.208, p = 5.26 × 10−5) and ACSL4 (cor = −0.152, p = 3.26 × 10−3) and was significantly positively correlated with ATP5G3 (cor = 0.116, p = 2.5 × 10−2), GLS2 (cor = 0.188, p = 2.64 × 10−4), NCOA4 (cor = 0.461, p = 0.589 × 10−21), PTGS2 (cor = 0.198, p = 1.2 × 10−4), NFE2L2 (cor = 0.441, p = 3.99 × 10−19), SLC7A11 (cor = 0.104, p = 4.55 × 10−2), TFRC (cor = 0.208, p = 5.24 × 10−5), FANCD2 (cor = 0.175, p = 7.23 × 10−4), and SAT1 (cor = 0.143, p = 5.72 × 10−3) (Figure 11 and Table 1). These results suggested that CCDC25 may not only participate in apoptosis, but also play an antitumor role by regulating ferroptosis.

FIGURE 11
www.frontiersin.org

Figure 11. The relationship between CCDC25 expression and ferroptosis key enzymes (A) GPX4, (B) ACSL4, (C) ATP5G3, (D) GLS2, (E) NCOA4, (F) PTGS2, (G) NFE2L2, (H) SLC7A11, (I) TFRC, (J) FANCD2, and (Q) SAT1.

TABLE 1
www.frontiersin.org

Table 1. The relationship between CCDC25 and ferroptosis key enzymes

Discussion

HCC is a serious global health problem, and Asia is particularly vulnerable to HCC as a region with a high incidence of HBV infection (22). The incidence of HCC will continue to increase in the coming decades (23). Therefore, it is necessary to find biomarkers that can be used for early diagnosis or as therapeutic targets for HCC, which may throw new light into the diagnosis and treatment of the disease.

Compared with normal tissues, CCDC25 expression was significantly reduced in HCC tissues of seven GEO data sets, which was confirmed by the TCGA-LIHC data. At the same time, its expression was reduced in many tumors. We suspected that it may play an inhibitory role in HCC or may play a protective role in normal liver tissue. Subsequently, we analyzed the relationship between CCDC25 expression and prognosis of HCC patients and found that HCC patients with low CCDC25 expression had a poor prognosis, which was consistent with our conjecture. Meanwhile, the calculation results of multiple data sets indicated that CCDC25 had a good diagnostic ability for HCC, and its diagnostic value was comparable to that of APF. Based on the above analysis, we have reason to believe that CCDC25 may be a promising prognostic indicator in HCC.

Co-expressed genes are generally involved in similar biological functions, so we explored the co-expressed gene network of CCDC25 in HCC through LinkedOmics. The enrichment analysis suggested that CCDC25 was mainly involved in the citrate cycle, fatty acid metabolism, amino acid metabolism, and PPAR signaling pathway. The citric acid cycle is the body’s main way of obtaining ATP, and cancer cells seem to prefer aerobic glycolysis, which is known as the Warburg effect (24, 25). If CCDC25 is involved in the citrate cycle, to which cancer cells are not prone, it may explain its downregulation in HCC. Many studies have shown that lipid (2628) and amino acid (29) metabolism disorders are closely related to the occurrence and progression of HCC. Low CCDC25 expression may promote metabolic reprogramming and thus contribute to the development of HCC. The PPAR signaling pathway (30) mainly includes members of the peroxisome proliferator-activated receptor family, which is the main regulator of glucose and lipid metabolism. In conclusion, CCDC25 with a high expression in normal liver tissue may be involved in three metabolic pathways to maintain body function, and its downregulation may promote metabolic reprogramming in favor of HCC.

In addition, the GSEA results showed that CCDC25 was associated with nucleotide metabolism, apoptosis, ubiquitination modification, and DNA damage binding. Controlled cell growth and programmed death are keys to homeostasis. Apoptosis is important for the elimination of excess, virus-infected, and damaged cells (31). CCDC25 may be a tumor suppressor gene that mediates apoptosis. Protein ubiquitination is a form of post-translational modification that regulates protein activity. It involves protein degradation, DNA damage repair, gene transcription, and signal transduction, and the dysregulation of this process is closely related to cancer (32). In conclusion, we found that CCDC25 may mainly play a role in metabolism, epigenetic modification, and regulation of programmed cell death. In conclusion, this study is the first to report that downregulated CCDC25 may inhibit the above-mentioned pathways and promote the occurrence and development of HCC.

The relationship between the TME and the cancer cells is very complex (33), and TME is closely related to tumor cell proliferation, apoptosis, and metastasis (34). It has been argued (35) that there are fundamental differences in the composition of immune cells in healthy tissue, adjacent tissue, and cancer tissue. A number of studies (36) have shown that B cells play an antitumor role by directly interacting with tumor cells or assisting other immune functions. CD8 + T cells are the main antitumor effector cells, while Treg cells often suppress antitumor immunity (37). We found that CCDC25 seems to increase antitumor immunity by increasing CD8 + T cells and decreasing Treg cells. Yang et al. reported (8) that the DNA of NETs, which are released by neutrophils, promotes cancer metastasis. Our findings suggested that patients with a high CCDC25 expression indeed had more neutrophil infiltration. Therefore, CCDC25 may negatively affect HCC by changing the degree of immune cell infiltration in TME.

It has been demonstrated that PDCD-1 (38), TIGIT (39), CTLA-4 (40), and TIM-3 (41) help tumor cell immune evasion through different mechanisms. Our study found that CCDC25 expression was negatively correlated with these four immune checkpoints in HCC, suggesting that CCDC25 may negatively regulate these immune checkpoints to play an antitumor role.

Sorafenib (42) is the first first-line chemotherapy agent for advanced HCC approved by the US Food and Drug Administration (FDA). However, advanced HCC patients will fail treatment due to sorafenib tolerance (22). In this study, we found that HCC patients with high CCDC25 expression may have higher sorafenib tolerance, and patients with low CCDC25 expression may have better sensitivity to sorafenib treatment. In addition, we found that patients with low CCDC25 expression had better sensitivity to Lapatinib, Gefitinib, and Selumetinib, while Ponatinib, Axitinib, Bleomycin, Bosutinib, and Methotrexate may be potentially sensitive drugs for HCC patients with a high expression of CCDC25.

Ferroptosis (43) is a newly discovered type of programmed cell death characterized by iron-dependence lipid peroxidation. Emerging evidence shows that ferroptosis is involved in various liver diseases, especially HCC (44). Glutathione peroxidase 4 (GPX4) (45) has been shown to inhibit ferroptosis in cancer cells. The decreased activity of GPX4 (46) can hinder the metabolism of lipid peroxides, which leads to a breakdown of cell membrane integrity and ultimately ferroptosis. Interestingly, GPX4 played a carcinogenic role in HCC and was negatively correlated with CCDC25 in this study. Subsequently, we found that CCDC25 was associated with many ferroptosis key enzymes, which suggested that CCDC25 might be a key molecule in the development of ferroptosis.

This study is the first to reveal the relationship between CCDC25 and HCC; it suggested that CCDC25 might have a complex impact on the prognosis, diagnosis, and TME in HCC. In conclusion, CCDC25 is expected to be a promising diagnostic and prognostic biomarker for HCC patients, and its function in HCC deserves further exploration.

Data Sharing Statement

All data generated or analyzed during this study are available in the TCGA database (https://portal.gdc.cancer.gov) and GEO database (https://www.ncbi.nlm.nih.gov/geo/).

Data availability statement

Publicly available data sets were analyzed in this study. These data sets can be found in https://portal.gdc.cancer.gov https://www.ncbi.nlm.nih.gov/geo/.

Ethics statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

HD and JZ analyzed and interpreted the data and were major contributors to the writing of the manuscript. YZ refined the language of this article. JL and QX collected data. WH, FW, and XX provided technical support. YZ, who provided ideas and financial support, was our corresponding author. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (82060800), Doctoral Students Training Research Fund of Lanzhou University Second Hospital (YJS-BD-32), Science and Technology Projects in Chengguan District of Lanzhou City (2014-4-4), Gansu Provincial Science Foundation for Young Scholars of China (21JR1RA161, 21JR1RA149, 21JR11RA114), and Cuiying Scientific and Technological Innovation Program of Lanzhou University Second Hospital (CY2019-BJ02).

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. Lee IC, Huang J-Y, Chen T-C, Yen CH, Chiu NC, Hwang HE, et al. Evolutionary learning-derived clinical-radiomic models for predicting early recurrence of hepatocellular carcinoma after resection. Liver Cancer. (2021) 10(6):572–82. doi: 10.1159/000518728

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, 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(3):209–49. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Kulik L, El-Serag HB. Epidemiology and management of hepatocellular carcinoma. Gastroenterology. (2019) 156(2):477–91.e1. doi: 10.1053/j.gastro.2018.08.065

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Javan H, Dayyani F, Abi-Jaoudeh N. Therapy in advanced hepatocellular carcinoma. Semin Intervent Radiol. (2020) 37(5):466–74. doi: 10.1055/s-0040-1719187

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Roessler S, Long EL, Budhu A, Chen Y, Zhao X, Ji J, et al. Integrative genomic identification of genes on 8p associated with hepatocellular carcinoma progression and patient survival. Gastroenterology. (2012) 142(4):957–66.e12. doi: 10.1053/j.gastro.2011.12.039

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Proungvitaya S, Klinthong W, Proungvitaya T, Limpaiboon T, Jearanaikoon P, Roytrakul S, et al. High expression of CCDC25 in cholangiocarcinoma tissue samples. Oncol Lett. (2017) 14(2):2566–72. doi: 10.3892/ol.2017.6446

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Chanakankun R, Proungvitaya T, Chua-On D, Limpaiboon T, Roytrakul S, Jusakul A, et al. Serum coiled-coil domain containing 25 protein as a potential screening/diagnostic biomarker for cholangiocarcinoma. Oncol Lett. (2020) 19(1):930–42. doi: 10.3892/ol.2019.11162

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Yang L, Liu Q, Zhang X, Liu X, Zhou B, Chen J, et al. DNA of neutrophil extracellular traps promotes cancer metastasis via CCDC25. Nature. (2020) 583(7814):133–8. doi: 10.1038/s41586-020-2394-6

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Dickson I. NETs promote liver metastasis via CCDC25. Nat Rev Gastroenterol Hepatol. (2020) 17(8):451. doi: 10.1038/s41575-020-0345-1

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Liu R, Zhao E, Wang F, Cui H. CCDC25: precise navigator for neutrophil extracellular traps on the prometastatic road. Signal Transduct and Targeted Ther. (2020) 5(1):162. doi: 10.1038/s41392-020-00285-6

CrossRef Full Text | Google Scholar

11. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. (2017) 77(21):e108–10. doi: 10.1158/0008-5472.CAN-17-0307

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Lian Q, Wang S, Zhang G, Wang D, Luo G, Tang J, et al. HCCDB: a database of hepatocellular carcinoma expression atlas. Genomics Proteomics Bioinformatics. (2018) 16(4):269–75. doi: 10.1016/j.gpb.2018.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Lánczky A, Győrffy B. Web-based survival analysis tool tailored for medical research (KMplot): development and implementation. J Med Internet Res. (2021) 23(7):e27633. doi: 10.2196/27633

CrossRef Full Text | Google Scholar

14. Menyhárt O, Nagy Á, Győrffy B. Determining consistent prognostic biomarkers of overall survival and vascular invasion in hepatocellular carcinoma. R Soc Open Sci. (2018) 5(12):181006. doi: 10.1098/rsos.181006

CrossRef Full Text | Google Scholar

15. Vasaikar SV, Straub P, Wang J, Zhang B. Linkedomics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. (2018) 46(D1):D956–63. doi: 10.1093/nar/gkx1090

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. (2005) 102(43):15545–50. doi: 10.1073/pnas.0506580102

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Ru B, Wong CN, Tong Y, Zhong JY, Zhong SSW, Wu WC, et al. TISIDB: an integrated repository portal for tumor-immune system interactions. Bioinformatics (Oxford, England). (2019) 35(20):4200–2. doi: 10.1093/bioinformatics/btz210

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Modhukur V, Iljasenko T, Metsalu T, Lokk K, Laisk-Podar T, Vilo J. Methsurv: a web tool to perform multivariable survival analysis using DNA methylation data. Epigenomics. (2018) 10(3):277–88. doi: 10.2217/epi-2017-0118

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Tella SH, Kommalapati A, Mahipal A, Jin Z. First-line targeted therapy for hepatocellular carcinoma: role of atezolizumab/bevacizumab combination. Biomedicines. (2022) 10(6):1304. doi: 10.3390/biomedicines10061304

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Chakraborty E, Sarkar D. Emerging therapies for hepatocellular carcinoma (HCC). Cancers (Basel). (2022) 14(11):2798. doi: 10.3390/cancers14112798

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Bekric D, Ocker M, Mayr C, Stintzing S, Ritter M, Kiesslich T, et al. Ferroptosis in hepatocellular carcinoma: mechanisms, drug targets and approaches to clinical translation. Cancers (Basel). (2022) 14(7):1826. doi: 10.3390/cancers14071826

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Chen X, Yang G, Guo X, Zhang J, Sun W, Liu D, et al. DJ-1/FGFR-1 signaling pathway contributes to sorafenib resistance in hepatocellular carcinoma . Oxid Med Cell Longev. (2022) 2022:2543220. doi: 10.1155/2022/2543220

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Petrick JL, Kelly SP, Altekruse SF, McGlynn KA, Rosenberg PS. Future of hepatocellular carcinoma incidence in the United States forecast through 2030 . J Clin Oncol. (2016) 34(15):1787–94. doi: 10.1200/JCO.2015.64.7412

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Todisco S, Convertini P, Iacobazzi V, Infantino V. TCA cycle rewiring as emerging metabolic signature of hepatocellular carcinoma. Cancers (Basel). (2019) 12(1):68. doi: 10.3390/cancers12010068

CrossRef Full Text | Google Scholar

25. Du D, Liu C, Qin M, Zhang X, Xi T, Yuan S, et al. Metabolic dysregulation and emerging therapeutical targets for hepatocellular carcinoma. Acta Pharm Sin B. (2022) 12(2):558–80. doi: 10.1016/j.apsb.2021.09.019

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Paul B, Lewinska M, Andersen JB. Lipid alterations in chronic liver disease and liver cancer. JHEP Rep. (2022) 4(6):100479. doi: 10.1016/j.jhepr.2022

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Satriano L, Lewinska M, Rodrigues PM, Banales JM, Andersen JB. Metabolic rearrangements in primary liver cancers: cause and consequences. Nat Rev Gastroenterol Hepatol. (2019) 16(12):748–66. doi: 10.1038/s41575-019-0217-8

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Calvisi DF, Wang C, Ho C, Ladu S, Lee SA, Mattu S. Increased lipogenesis, induced by AKT-mTORC1-RPS6 signaling, promotes development of human hepatocellular carcinoma. Gastroenterology. (2011) 140(3):1071–83. doi: 10.1053/j.gastro.2010.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Chisari A, Golán I, Campisano S, Gélabert C, Moustakas A, Sancho P, et al. Glucose and amino acid metabolic dependencies linked to stemness and metastasis in different aggressive cancer types. Front Pharmacol. (2021) 12:723798. doi: 10.3389/fphar.2021.723798

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Mello T, Simeone I, Galli A. Mito-nuclear communication in hepatocellular carcinoma metabolic rewiring. Cells. (2019) 8(5):417. doi: 10.3390/cells8050417

CrossRef Full Text | Google Scholar

31. Arbuthnot P, Capovilla A, Kew M. Putative role of hepatitis B virus X protein in hepatocarcinogenesis: effects on apoptosis, DNA repair, mitogen-activated protein kinase and JAK/STAT pathways. J Gastroenterol Hepatol. (2000) 15(4):357–68. doi: 10.1046/j.1440-1746.2000.02069.x

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zhang X, Meng T, Cui S, Liu D, Pang Q, Wang P. Roles of ubiquitination in the crosstalk between tumors and the tumor microenvironment (review). Int J Oncol. (2022) 61(1):84. doi: 10.3892/ijo.2022.5374

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Wu F, Wei H, Liu G, Zhang Y. Bioinformatics profiling of five immune-related lncRNAs for a prognostic model of hepatocellular carcinoma. Front Oncol. (2021) 11:667904. doi: 10.3389/fonc.2021.667904

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Liu R, Li R, Yu H, Liu J, Zheng S, Li Y, et al. NTF3 Correlates with prognosis and immune infiltration in hepatocellular carcinoma. Front Med (Lausanne). (2021) 8:795849. doi: 10.3389/fmed.2021.795849

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Woller N, Engelskircher SA, Wirth T, Wedemeyer H. Prospects and challenges for T cell-based therapies of HCC. Cells. (2021) 10(7):1651. doi: 10.3390/cells10071651

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Qin M, Wang D, Fang Y, Zheng Z, Liu X, Wu F, et al. Current perspectives on B lymphocytes in the immunobiology of hepatocellular carcinoma. Front Oncol. (2021) 11:647854. doi: 10.3389/fonc.2021.647854

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Ohue Y, Nishikawa H. Regulatory T (Treg) cells in cancer: can Treg cells be a new therapeutic target?. Cancer Sci. (2019) 110(7):2080–9. doi: 10.1111/cas.14069

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Ishida Y, Agata Y, Shibahara K, Honjo T. Induced expression of PD-1, a novel member of the immunoglobulin gene superfamily, upon programmed cell death. EMBO J. (1992) 11(11):3887–95. doi: 10.1002/j.1460-2075.1992.tb05481.x

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Johnston RJ, Comps-Agrar L, Hackney J, Yu X, Huseni M, Yang Y, et al. The immunoreceptor TIGIT regulates antitumor and antiviral CD8(+) T cell effector function. Cancer Cell. (2014) 26(6):923–37. doi: 10.1016/j.ccell.2014.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Ozer M, George A, Goksu SY, George TJ, Sahin I. The role of immune checkpoint blockade in the hepatocellular carcinoma: a review of clinical trials. Front Oncol. (2021) 11:801379. doi: 10.3389/fonc.2021.801379

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Makaremi S, Asadzadeh Z, Hemmat N, Baghbanzadeh A, Sgambato A, Ghorbaninezhad F, et al. Immune checkpoint inhibitors in colorectal cancer: challenges and future prospects. Biomedicines. (2021) 9(9):1075. doi: 10.3390/biomedicines9091075

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Lu Y, Chan Y-T, Tan H-Y, Zhang C, Guo W, Xu Y, et al. Epigenetic regulation of ferroptosis via ETS1/miR-23a-3p/ACSL4 axis mediates sorafenib resistance in human hepatocellular carcinoma. J Exp Clin Cancer Res. (2022) 41(1):3. doi: 10.1186/s13046-021-02208-x

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE, et al. Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell. (2012) 149(5):1060–72. doi: 10.1016/j.cell.2012.03.042

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Chen S, Zhu J-Y, Zang X, Zhai YZ. The emerging role of ferroptosis in liver diseases. Front Cell Dev Biol. (2021) 9:801365. doi: 10.3389/fcell.2021.801365

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Zhang X, Sui S, Wang L, Li H, Zhang L, Xu S, et al. Inhibition of tumor propellant glutathione peroxidase 4 induces ferroptosis in cancer cells and enhances anticancer effect of cisplatin. J Cell Physiol. (2020) 235(4):3425–37. doi: 10.1002/jcp.29232

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Jiang X, Stockwell BR, Conrad M. Ferroptosis: mechanisms, biology and role in disease. Nat Rev Mol Cell Biol. (2021) 22(4):266–82. doi: 10.1038/s41580-020-00324-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: coiled-coil domain containing 25, hepatocellular carcinoma, prognosis, diagnosis, immune infiltrates, ferroptosis

Citation: Deng H, Zhang J, Zheng Y, Li J, Xiao Q, Wei F, Han W, Xu X and Zhang Y (2022) CCDC25 may be a potential diagnostic and prognostic marker of hepatocellular carcinoma: Results from microarray analysis. Front. Surg. 9:878648. doi: 10.3389/fsurg.2022.878648

Received: 18 February 2022; Accepted: 31 August 2022;
Published: 23 September 2022.

Edited by:

Pasquale Cianci, Azienda Sanitaria Localedella Provincia di Barletta Andri Trani (ASL BT), Italy

Reviewed by:

Qi Zhang, Huazhong University of Science and Technology, China
Shijie Li, Zhejiang University, China

© 2022 Deng, Zhang, Zheng, Li, Xiao, Wei, Han, Xu and Zhang. 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: Youcheng Zhang zhangychmd@126.com

These authors have contributed equally to this work

Specialty Section: This article was submitted to Surgical Oncology, a section of the journal Frontiers in Surgery

Abbreviations HCC, hepatocellular carcinoma; CCDC25, coiled-coil domain containing 25; CCA, cholangiocarcinoma; NETs, neutrophil extracellular traps; TIMER, Tumor IMmune Estimation Resource; TCGA, The Cancer Genome Atlas; GEO, Gene Expression Omnibus; ROC, receiver operating characteristic curve; AFP, alpha fetoprotein; KEGG, Kyoto Encyclopedia of Genes and Genomes; GSEA, Gene Set Enrichment Analysis; BLCA, bladder urothelial carcinoma; BRCA, breast invasive carcinoma; COAD, colon adenocarcinoma; KICH, kidney chromophobe; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; PRAD, prostate adenocarcinoma; READ, rectum adenocarcinoma; THCA, thyroid carcinoma; UCEC, uterine corpus endometrial carcinoma; HNSC, head and neck squamous cell carcinoma; SKCM, the skin cutaneous melanoma; EGA, European Genome-Phenome Archive; OS, overall survival; PFS, progression-free survival; DSS, disease-specific survival; TAM, tumor-associated macrophages; FDA, Food and Drug Administration

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.