Research Paper Volume 14, Issue 23 pp 9550—9578

Comprehensive analysis of prognostic value, relationship to cell cycle, immune infiltration and m6A modification of ZSCAN20 in hepatocellular carcinoma

Fang Jiayu1,2, *, , Yike Jiang2, *, , Xuanrui Zhou2, , Minqin Zhou2, , Jingying Pan2, , Yun Ke2, , Jing Zhen2, , Da Huang3, , Weifan Jiang4, ,

  • 1 Second Affiliated Hospital of Nanchang University, Nanchang, China
  • 2 Second College of Clinical Medicine, Nanchang University, Nanchang, China
  • 3 Department of Thyroid Surgery, Second Affiliated Hospital of Nanchang University, Nanchang, China
  • 4 Department of Urology, Second Affiliated Hospital of Nanchang University, Nanchang, China
* Equal contribution

Received: May 25, 2022       Accepted: September 17, 2022       Published: December 3, 2022      

https://doi.org/10.18632/aging.204312
How to Cite

Copyright: © 2022 Jiayu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Hepatocellular carcinoma (HCC) is a common tumor across the globe with a high mortality rate. ZSCAN20 is a ZNF transcription factor, a key determinant of gene expression. Nonetheless, the mechanism of ZSCAN20 as a potential clinical biomarker and therapeutic target for HCC is not understood. Here, TIMER, TCGA, ICGC databases and immunohistochemical (IHC) and Western Blot found ZSCAN20 mRNA and protein levels were upregulated. Additionally, Kaplan-Meier Plotter, GEPIA and TCGA databases showed high ZSCAN20 expression was related to the short survival time of HCC patients. Multivariate Cox analysis exposed that ZSCAN20 can act as an independent prognostic factor. We observed methylation level of ZSCAN20 was associated with the clinicopathological characteristics and prognosis of HCC patients through UALCAN. Furthermore, enrichment examination exposed functional association between ZSCAN20 and cell cycle, immune infiltration. Functional experiments showed that interference with ZSCAN20 significantly reduced the invasion, migration and proliferation abilities of HCC cells. An immune infiltration analysis showed that ZSCAN20 was associated with immune cells, particularly T cells. The expression of ZSCAN20 was correlated with poor prognosis in the Regulatory T-cell. And Real-Time RT-PCR analysis found interference with ZSCAN20 significantly reduced the expression of some chemokines. Finally, the TCGA and ICGC data analysis suggested that the ZSCAN20 expression was greatly related to m6A modifier related genes. In conclusion, ZSCAN20 can serve as a prognostic biomarker for HCC and provide clues about cell cycle, immune infiltration, and m6A modification.

Introduction

Hepatocellular carcinoma (HCC) accounts for 90% of primary liver cancer [1]. HCC is the second reason to lead cancer-related demise all over the world, because of the accelerated progress of HCC [24]. Also, many patients have reached the advanced stage when they are diagnosed, bringing about a poor prognosis for HCC patients, with a 5-year survival rate of only 18% [5]. If patients get an early treatment, it will extremely better the survival rate of patients. In view of the discovery of HCC biomarkers, the diagnosis of HCC has made accelerated progress. For example, Alpha Fetoprotein (AFP) check is utilized in the early diagnosis of HCC, but it still has a certain misdiagnosis rate [6]. Early patient with HCC can be treated by surgical resection and ablation, but the tumor recurrence rate is still high, which is not so ideal [7]. Meanwhile, more and more evidence reveals that immunotherapy including immune checkpoint blockade can become a likely optional therapy for HCC [810]. However, useful immunotherapy targets for HCC are still deficient. Therefore, it is still crucial to find some more effective HCC biomarkers and therapeutic targets.

The ZSCAN family belongs to zinc finger transcription factors. As the family that contains the most transcription factors [11], zinc finger transcription factors are characterized by finger-like DNA binding domains, which require one or more zinc ions to stabilize the structure and are used in many biological processes and play an important role [12]. The human ZSCAN transcription factor family has 54 identified members. As a member of this family, ZSCAN20 has a chromosome located at 1p35.1. ZSCAN transcription factors may be involved in angiogenesis, cell apoptosis, cell differentiation, cell migration and invasion, cell proliferation, stem cell characteristics, and chemotherapy sensitivity. Increasingly, studies are detecting that abnormal ZSCAN transcription factor expression in liver cancer is linked to significant biological functions. In HCC, ZSCAN3 directly binds to the promoter of β-catenin to activate its transcription and promote cell proliferation [13], and ZSCAN33 can inhibit HCC cell growth, migration, and invasion [14]. In addition, ZSCAN35 is directly related to carcinoembryonic antigen and takes an active part in the growth of HCC metastasis; ZSCAN35 may promote liver metastasis [15]. Furthermore, recent studies have found that ZSCAN20 may become a new target for tumor treatment. The article reports that ZSCAN20 is highly expressed in low-grade gliomas and can be used as a part of the model to predict the prognosis of low-grade gliomas. Nevertheless, the effect and mechanism of ZSCAN20 in HCC haven’t been reported, and its relationship with the prognosis is still unclear.

There are few studies on ZSCAN20, and the mechanism of ZSCAN20’s role in cancer has not yet been elucidated. As a consequence, here we explore the specific mechanism of ZSCAN20 as a prognosis and potential therapeutic target for HCC patients.

Materials and Methods

Patients and tumor specimens

We gathered 40 pairs of paraffin-embedded HCC and corresponding normal tissues from 40 patient samples from the Second Affiliated Hospital of Nanchang University from January 2018 to January 2021. The above patients did not receive adjuvant chemotherapy. Entire patients endorsed a knowledgeable consent form. Also, this research was accepted by the Ethics Committee of the Second Affiliated Hospital of Nanchang University.

Data collection

The TCGA database (https://cancergenome.nih.gov) [16] includes sequencing and clinicopathological data from more than 32 cancer patients. We evaluated the ZSCAN20 expression pattern according to various cancer research from the LIHC data of the Cancer Genome Atlas (TCGA), including 377 cancer samples and 50 normal samples. The main goal of ICGC is to comprehensively clarify the genomic changes in many cancers that contribute to the global burden of human disease [17]. ICGC collects tumor data of 50 different cancer types. The used file type was [LINC-JP] Liver Cancer-NCC, JP data set, which includes 243 tumor samples and 260 clinical samples.

TIMER database analysis

TIMER (https://cistrome.shinyapps.io/timer/) [1821] is a broad capability library for analyzing the immune infiltration of different cancer types systematically. It concludes 1,000 samples of 32 cancer types. The degree of immune cell infiltration in tumor tissues can be appraised from RNA-seq expression profile data, reflecting the correlation hidden in cancer and immune cells. We analyzed the expression of ZSCAN20 in various cancer types and the interaction between ZSCAN20 expression and immune infiltration, along with conventional immune cells. Then, we used the SCAN module of TIMER database to analyze the relationship between ZSCAN20 CNV and immune cell infiltration. In addition, the TIMER database was also utilized to study the interaction hidden in ZSCAN20 expression and various immune gene marker sets of immune cells.

Cell culture

Transfection Human HCC cell lines LM3 were maintained with 5% CO2 at 37° C in DMEM (HyClone Germany) with 10% fetal bovine serum (GIBCO USA).

Immunohistochemistry

The HCC tissues and paired adjacent tissues fixed with 10% formalin and implanted in paraffin were cut into 4um thick sections. After deparaffinization, rehydration, and microwave heating in sodium citrate buffer (10 mmol/L, pH 6.0) for 25 minutes to restore the antigen, the sections were sealed with goat serum for 30 minutes. Then, the sections were incubated overnight with anti-ZSCAN20 polyclonal antibodies (ZSCAN20-LifeSpan BioSciences-Catalog number: C660829) at 4° C. Next, the HRP conjugated second antibody (booster) was allowed to stand at room temperature for 2 hours. Subsequently, immunostaining was performed using a two-step method. Pathologists who were unaware of the clinical parameters assessed the staining intensity and the percentage of positive cells semi-quantitatively [22].

Western blot

For the preparation mean of Western Blot total protein extract, saw [23]. In general, protein extraction on ice using RIPA buffer (Beyotime, Shanghai, China) which was a mixture of protease and inhibitor (Thermo Fisher Scientific, New York, USA). After centrifugation, the protein concentration was observed using the BCA Protein Detection Kit (Thermo Scientific, Waltham, MA, USA), and an equal amount of protein was isolated on SDS-PAGE by electrophoresis, which was transferred to a polyvinylidene fluoride (PVDF) membrane. Then incubate overnight at 4° C applying the primary antibody, followed by washing the closed membrane 3 times in TBST, while incubating with the secondary antibody for 1h at room temperature. After the above process was accomplished, the protein expression was checked employing electro-chemiluminescence (ECL) assay.

Quantitative real-time PCR

Real-time PCR was performed on fresh frozen samples. Total RNA was isolated from tissues using the TRIzol reagent (Thermo Fisher Scientific). Then, according to standard methods, qRTPCR analysis was performed to observe the expression level of ZSCAN20. Primers for ZSCAN20 from 5’ to 3’: (F), (R). Primers for β-actin from 5’ to 3’: (F), (R).

HPA analysis

HPA (https://www.Proteinatlas.org/) uses antibody methods for immunostaining of tissues and cell lines, along with contrasting expression analysis of proteins in normal and tumor tissues [24]. In this research, the expression of ZSCAN20 was disclosed by using the protein expression module in the HPA database, and the immunohistochemical outcomes of ZSCAN20 in tumor tissues and normal tissues were concluded.

UALCAN analysis

UALCAN (http://ualcan.path.uab.edu/) is an interactive evaluation and excavation web on the Internet, which can be used to analyze the relationship between tumor and normal specimens and the relative expression of genes with diverse clinicopathological characteristics [25, 26]. UALCAN was used to evaluate the differential expression of ZSCAN20 in HCC tissues and normal tissues. The research further investigated various clinicopathological characteristics of ZSCAN20 and the methylation of ZSCAN20 promoter.

MethSurv database analysis

MethSurv (https://biit.cs.ut.ee/methsurv/) is an open online instrument, which can evaluate the prognostic value of CpG methylation data. The methylation data comes from The Cancer Genome Atlas (TCGA), which can further provide overall survival (OS) of DNA methylation levels. We could use this database to analyze the DNA methylation data of ZSCAN20 and evaluate the prognostic value of a single CpG site in HCC patients [27, 28].

GEPIA database analysis

Gene expression profiling transactional estimation (GEPIA) (http://gepia.cancer-pku.cn/) is an interactional website application consisting of thousands of tumor and normal tissue sample data, which can be used to visualize clinicopathological characteristics. The tumor data comes from the TCGA database [24]. We used GEPIA to parallel tumors with normal tissues and analyzed related prognosis in this research. In addition, the relationship between ZSCAN20 and CD274, CTLA4, PDCD1 was firmed in the “correlation analysis”. The overall survival of the m6A-related genes HNRNPA2B1, RBMX, YTHDC1 and YTHDF2 in LIHC were also generated through this locus.

TISIDB analysis

Tumor and immune system interaction database (TISIDB) (http://cis.hku.hk/TISIDB) is a portal site, which can be used to evaluate how the tumor interacts with the immune system, and it integrates various heterogeneous data types [29]. Using the TISIDB database, we investigated the features of ZSCAN20 and its function in tumor immune interaction. The “Chemokine” module was used to analyze spearman correlations between ZSCAN20 and chemokines across LIHC.

Kaplan-Meier plotter database analysis

The Kaplan-Meier database (http://kmplot.com/analysis/) can be used to appraise the effect of genes on the survival of cancer tissue samples [30, 31]. We used 364 LIHC samples to evaluate the interaction hidden in ZSCAN20 expression and overall survival (OS), relapse-free survival (RFS), and disease-specific survival (DSS). We also explored the difference in LIHC patients’ survival under different immune cell subtypes. A 95% confidence interval and log-rank p < 0.05 was studied analytically.

LinkedOmics database analysis

The website named LinkedOmics database has been utilized to analyze 32 TCGA cancer-related data sets (http://www.linkedomics.org/login.php) [3235]. Our team utilized the LinkFinder module of LinkedOmics to research the differently expressed genes connected with ZSCAN20 in TCGA LIHC (n = 515). Perform statistical analysis on the conclusions and display them in volcano maps, heat maps, and scatter plots. Gene set enrichment analysis (GSEA) has been utilized for LinkedOmics functional modules to conduct gene ontology (GO) analysis and KEGG pathway analysis. False discovery rate (FDR) less than 0.01 is a noteworthy expression, P-values less than 0.05 is a significantly related gene.

GSEA software analysis

GSEA is an analysis way that confirms if a previously confirmed set of genes reveals a statistically significant and consistent difference amongst two phenotypes [34, 36]. The R package clusterProfiler (4.1.0) [35] has been utilized for performing GSEA to clarify the noteworthy functional and pathway differences amongst ZSCAN20 high expression group and ZSCAN20 low expression group. We performed 1,000 analyses of the genome arrangement and used the expression level of ZSCAN20 mRNA as a phenotypic marker. After running GSEA, P < 0.05 and FDR < 0.25 are both meaningful.

Protein-protein interaction (PPI) network analysis

The online inquire implement (STRING, https://string-db.org) for finding interacting genes which are utilized for PPI network construction and pivot gene screening [37, 38]. Besides, with the purpose of exploring the correlation amongst the top 500, our team used STRING database for analysis, medium confidence=0.4 for screening, and visualization using Cytoscape software. To search cluster sub-networks, the Cytoscape Molecular Complex Detection (MCODE) plug-in has been utilized [39]. The acquiescent parameters are listed below: degree cutoff = 5, node score cutoff = 0.2, k-core = 9, maximum depth = 100.

Transwell assay

The stably transfected cells 48 hours after seeding were subjected to Transwell invasion assay to evaluate the invasive ability of ZSCAN20. For migration, 6 x 104 cells were plated in the upper chamber with serum-free medium. For invasion, 1 x 105 cells were placed in Matrigel-coated chambers (BD Biosciences). Subsequently, after incubation at 37° C for 24 hours (to check migration) or 48 hours (to check invasion), Transwell inserts and invasive cells were fixed in methanol for 30 minutes, stained with 0.1% crystal violet, using a CCD with DP70. The light microscope of the system (Olympus Corp.) counted the number of cells in five random fields.

Colony formation assay

For the colony formation assay, 0.6×103 cells were seeded in 6-well plates. Once the appropriate colony size had formed, cells were fixed with 4% paraformaldehyde for 30 minutes, and after that stained by 1.0% crystal violet for 30 mins until formed visible clones. Colonies were stained with Giemsa and counted in 10 different fields [40].

Statistical analysis

R software (version 3.6.3/4.1.2) was utilized to complete all statistical analyses amongst our research. The “limma” and “beeswarm” packages of “R” and the rank-sum test were utilized to detect the difference of ZSCAN20 expressions between LIHC samples and normal samples. We explored the association between ZSCAN20 and clinicopathological features through logistic regression. Then, we used the Kaplan-Meier method to generate survival curves, and performed univariate and multivariate Cox regression analysis in R. Our team established time-dependent receiver operating characteristic curves (ROC) for comparing disparate survival factors. The graphics part below the curve (AUC) was obtained to utilize the R package “pROC”. The p value < 0.05 was affirmed statistically significant. In order to study the discrepant expression of m6A-related genes amongst ZSCAN20 high expression group and ZSCAN20 low expression group, we utilized T test analysis.

Results

ZSCSN20 expression is increased in HCC patients

The difference between ZSCAN20 expression in tumor and normal tissues was detected. Initially, we used the TCGA, TIMER and ICGC databases to assess the difference between ZSCAN20 expression in specific tumor types. We discovered that the expression of ZSCAN20 was up-regulated in LIHC and many other tumors based on the TIMER database (Figure 1A). In addition, compared with 50 adjacent normal tissue samples, ZSCAN20 in 377 LIHC samples was clearly up-regulated (Figure 1B). Subsequently, compared with adjacent normal samples, the expression of ZSCAN20 in 58 pairs of LIHC tumor samples was significantly increased (Figure 1C). The data set LIRI-JP (liver cancer-RIKEN, JP) we downloaded from the ICGC database consistently showed that the expression of LIHC in tumor samples increased (Figure 1D). In addition, we also used HPA database to detect that ZSCAN20 protein expression in tumor tissues was remarkably higher than normal tissues in the liver of patients with HCC (Figure 1E). Further studying ZSCAN20 protein expression in HCC tissues by IHC staining, we discovered that ZSCAN20 protein expression in HCC tissues was considerably higher than adjacent tissues (Figure 1F). For the aim of examine the expression level of ZSCAN20 in HCC tissues, we performed IHC staining and examined 40 pairs tumor tissues of paraffin embedded archived HCC and corresponding 40 pairs paracancerous tissues. Our data indicated that ZSCAN20 was significantly overexpressed in HCC tissues (Figure 1G). Eventually, we did a set of Western Blot experiments. As presented in Figure 1H, we dug out that in comparison with neighboring normal tissues, ZSCAN20 was more expressed in HCC tissues. Consistently, profit from qRT-PCR, our result demonstrated that the mRNA expression of ZSCAN20 was dramatically upregulated in HCC specimens in comparation with the level in adjacent normal tissues (Figure 1I). These results indicate that ZSCAN20 expression is up-regulated in HCC, and ZSCAN20 may play a crucial regulatory part in the development of HCC.

Expression of ZSCAN20 in HCC. (A) The expression level of ZSCAN20 in different types of tumor tissues and normal tissues in the TIMER database (P B) Expression levels of ZSCAN20 were higher than corresponding normal tissues in LIHC samples (TCGA-LIHC) (p=1.476e-13). (C) ZSCAN20 expression in 50 paired LIHC tissues and corresponding adjacent non-tumor tissues (TCGA-LIHC) (p=1.37e-08). (D) Expression levels of ZSCAN20 were higher than corresponding normal tissues in LIHC samples by using ICGC-LIRI-JP liver datasets. (E) ZSCAN20 protein expression in normal and LIHC tissues (HPA). (F) Typical images of immunohistochemistry (IHC) in 35 pairs of LIHC tissues showing the protein expression of ZSCAN20 in LIHC and adjacent non-tumor tissues. (G) Diagram of ZSCAN20 staining score in IHC staining. (H) The protein expression of ZSCAN20 in tumor tissues and adjacent normal tissues was detected by Western Blot, and the coloring depth represented the level of protein expression. (I) Quantification of ZSCAN20 protein expression based on western blot.

Figure 1. Expression of ZSCAN20 in HCC. (A) The expression level of ZSCAN20 in different types of tumor tissues and normal tissues in the TIMER database (P < 0.05). (B) Expression levels of ZSCAN20 were higher than corresponding normal tissues in LIHC samples (TCGA-LIHC) (p=1.476e-13). (C) ZSCAN20 expression in 50 paired LIHC tissues and corresponding adjacent non-tumor tissues (TCGA-LIHC) (p=1.37e-08). (D) Expression levels of ZSCAN20 were higher than corresponding normal tissues in LIHC samples by using ICGC-LIRI-JP liver datasets. (E) ZSCAN20 protein expression in normal and LIHC tissues (HPA). (F) Typical images of immunohistochemistry (IHC) in 35 pairs of LIHC tissues showing the protein expression of ZSCAN20 in LIHC and adjacent non-tumor tissues. (G) Diagram of ZSCAN20 staining score in IHC staining. (H) The protein expression of ZSCAN20 in tumor tissues and adjacent normal tissues was detected by Western Blot, and the coloring depth represented the level of protein expression. (I) Quantification of ZSCAN20 protein expression based on western blot.

ZSCAN20 clinical parameters of HCC patients

Because we are still uncertain about the role of ZSCAN20 in HCC, studying the correlation between ZSCAN20 expression and clinicopathological features would help us to uncover ZSCAN20’s function in the development of HCC. Consequently, we used the UALCAN online tool to investigate the expression of ZSCAN20 in the patient group based on different clinical parameters. The results showed that ZSCAN20 expressed significantly different in various HCC samples, gender, age, cancer stages, tumor grade, nodal metastasis status, histological subtypes (Figure 2A2G). We used logistic regression analysis to further study ZSCAN20 expression connected with clinicopathological variables. The results indicated the high expression of ZSCAN20 was obviously associated with grade, stage and tumor size (Supplementary Table 1). All results show that ZSCAN20 expression has a close connection with clinicopathological factors such as tumor progression and metastasis.

Box-plots exploring the relationship between ZSCAN20 expression and clinicopathological characteristics (UALCAN). Increased ZSCAN20 expression was significantly with (A) sample type, (B) gender, (C) age, (D) cancer stage, (E) tumor grade, (F) nodal metastasis status. (G) histological subtype *P P P

Figure 2. Box-plots exploring the relationship between ZSCAN20 expression and clinicopathological characteristics (UALCAN). Increased ZSCAN20 expression was significantly with (A) sample type, (B) gender, (C) age, (D) cancer stage, (E) tumor grade, (F) nodal metastasis status. (G) histological subtype *P < 0.05; **P < 0.01; ***P < 0.001.

High ZSCAN20 expression is associated with poor prognosis of HCC

To understand the prognostic value of the expression of ZSCAN20 in HCC better, the Kaplan-Meier Plotter tool was used to detect the prognosis of ZSCAN20. Kaplan-Meier curve and log-rank test analysis showed that elevated ZSCAN20 mRNA levels correlated with overall survival (OS)(Figure 3A), relapse-free survival (RFS) (Figure 3B), and disease-specific survival (DSS) (Figure 3C). The consequence demonstrated that high ZSCAN20 expression corresponded to undesirable prognosis (OS, RFS, and DSS) in HCC patients (p < 0.05) (Figure 3A3C). We used the GEPIA and TCGA database to further verify the conjecture, and the results were consistent with the Kaplan-Meier Plotter tool (Figure 3D, 3E). Subsequently, the prediction accuracy of survival curve was determined by the construction of ROC curve. The graphics part below the ROC curve (AUC) for 1, 3, and 5 years were severally 0.664, 0.619, and 0.561. The result indicated favorable prediction accuracy (Figure 3F). In order to determine whether ZSCAN20 is an independent prognostic factor for HCC sufferers, we conducted the Cox regression analysis. Univariate Cox regression analysis showed that grade, T stage, and ZSCAN20 were noteworthy related to OS in patients with HCC. Multivariate Cox regression analysis revealed that ZSCAN20 was an independent factor affecting the prognosis in HCC (Table 1), also did forest plot (Figure 3G). The above results indicate that ZSCAN20 can be deemed to be an independent prognostic factor for the survival of HCC sufferers, and high expression of ZSCAN20 corresponds with the poor prognosis of HCC patients.

ZSCAN20 expression was associated with poor survival in HCC patients. (A–C) Survival curves using the Kaplan-Meier plotter are shown for OS, RFS, and DSS. (D) Survival curves using GEPIA is shown for OS. (E) Associations with overall survival and the expression of ZSCAN20 in TCGA patients. (F) ROC, receiver operating characteristic; ROC curve sensitivity and specificity analysis of LIHC, AUC, area under the curve. (G) A forest plot showed the correlation between ZSCAN20 expression and clinicopathological parameters in LIHC patients. *P P P

Figure 3. ZSCAN20 expression was associated with poor survival in HCC patients. (AC) Survival curves using the Kaplan-Meier plotter are shown for OS, RFS, and DSS. (D) Survival curves using GEPIA is shown for OS. (E) Associations with overall survival and the expression of ZSCAN20 in TCGA patients. (F) ROC, receiver operating characteristic; ROC curve sensitivity and specificity analysis of LIHC, AUC, area under the curve. (G) A forest plot showed the correlation between ZSCAN20 expression and clinicopathological parameters in LIHC patients. *P < 0.05; **P < 0.01; ***P < 0.001. HR, hazard ratio; CI, confidence interval; T, tumor; N, node, M, metastasis; OS, overall survival; AIC, Akaike information criterion.

Table 1. Univariate and multivariate COX regression analysis of factors associated with OS in HCC patients.

VariableUnivariate analysisMultivariate analysis
HR95%CIP-valueHR95%CIP-value
age1.0070.990-1.0240.4411.0090.991-1.0260.337
gender0.8390.536-1.3140.4431.1750.722-1.9120.516
grade1.0730.795-1.4490.6451.1190.814-1.5380.488
stage1.8091.426-2.294< 0.0011.2040.520-2.7900.665
T1.7671.415-2.207< 0.0011.3840.635-3.0180.414
M3.8921.223-12.3860.0212.1120.561-7.9510.269
ZSCAN2016.8834.940-57.699< 0.00110.7522.853-40.530< 0.001
OS, overall survival; HR, hazard ratio; CI, confidence interval; T, tumor; N, node; M, metastasis; Bold values indicate P-values<0.05.

DNA methylation level of ZSCAN20 in HCC patients

According to the report, DNA methylation is a chemical modification of DNA. For the purpose of studying the expression level of ZSCAN20 methylation in HCC, we based on the methylation characteristics of ZSCAN20 promoter in HCC patients and then UALCAN database was utilized to verify the methylation level of ZSCAN20 promoter in HCC. We found that the ZSCAN20 promoter in HCC was hypermethylated compared with normal tissues (Figure 4A). At the same time, it was worth noting that the ZSCAN20 promoter methylation profile based on gender, age, cancer stage, tumor grade and nodal metastasis status showed significant differences (all P < 0.01) (Figure 4B4F). In addition, we investigated the relationship between ZSCAN20 methylation level and OS. MethSurv analysis showed that the 6 CpG sites stood on CpG islands have an undesirable prognosis, containing cg02487499, cg14041778, cg18028289, cg18085998, cg19855402, and cg26805255 (Figure 4G4L). These results testify that the methylation of ZSCAN20 promoter is bound up with HCC’s occurrence and development.

Correlation between ZSCAN20 promoter methylation level and prognostic value of DNA methylation in HCC. (A) normal vs primary tumor, (B) gender, (C) age, (D) cancer stage, (E) tumor grade, (F) lymph node metastasis status; *P P P G), cg14041778 (H), cg18028289 (I), cg18085998 (J), cg19855402 (K) and cg26805255 (L) correlated with worse OS.

Figure 4. Correlation between ZSCAN20 promoter methylation level and prognostic value of DNA methylation in HCC. (A) normal vs primary tumor, (B) gender, (C) age, (D) cancer stage, (E) tumor grade, (F) lymph node metastasis status; *P < 0.05; **P < 0.01;***P < 0.001. High methylation level of cg02487499 (G), cg14041778 (H), cg18028289 (I), cg18085998 (J), cg19855402 (K) and cg26805255 (L) correlated with worse OS.

Co-expression networks, GO and KEGG enrichment analysis reveal pathways Of ZSCAN20 in HCC

For the sake of better exploring the biological functions of ZSCAN20 in HCC, “LinkFinder” in LinkedOmics was utilized to analyze 371 LIHC patients’ mRNA samples in TCGA. In accordance with the volcano map, there were 6713 truly positively correlated genes (red dots), and 2689 truly negatively correlated genes (green dots) (FDR < 0.01, p-value < 0.05) (Figure 5A). In addition, the first 50 negatively significant and positively significant related genes are shown in the heat maps (Figure 5B, 5C). We further explored the potential functional pathways through Linkedomics. The GO annotation showed that the co-expressed genes of ZSCAN20 mainly in chromosome separation, cell cycle G2/M phase transition, p53 binding etc. (Figure 5D5F). KEGG pathway analysis showed that it was rich in homologous recombination, cell cycle, Fanconi anemia pathway, DNA replication, etc. (Figure 5G).

Co-expression genes of ZSCAN20 in HCC. (A) Volcano plot of genes highly correlated with ZSCAN20 identified by the Pearson test in LIHC. Red and green dots represent genes significantly positively and negatively correlated with ZSCAN20, respectively. Heatmaps of the top 50 genes (B) positively and (C) negatively correlated with ZSCAN20. (D–G) Significantly enriched GO and KEGG pathways of ZSCAN20. GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes.

Figure 5. Co-expression genes of ZSCAN20 in HCC. (A) Volcano plot of genes highly correlated with ZSCAN20 identified by the Pearson test in LIHC. Red and green dots represent genes significantly positively and negatively correlated with ZSCAN20, respectively. Heatmaps of the top 50 genes (B) positively and (C) negatively correlated with ZSCAN20. (DG) Significantly enriched GO and KEGG pathways of ZSCAN20. GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes.

In contemplation of seeking the possible biological pathways of ZSCAN20 regulation in HCC, we selected the high and low ZSCAN20 expression groups to perform GSEA analysis. Cell cycle pathways, like cell cycle, Wnt signaling pathway; some popular cancer pathways, such as MAPK signaling pathway; immune-related pathways, like T cell receptor signaling pathways, are very abundant (Supplementary Table 2) (Nom p < 0.05; FDR < 0.05; NES > 1.7). These results indicate that the ZSCAN20 expression network acts a vital part on the immunity and cell cycle mechanism of HCC (Figure 6A6I).

GSEA results showed differential enrichment of genes with high ZSCAN20 expression. (A) Cell cycle. (B) DNA replication. (C) WNT signaling pathway. (D) JAK STAT signaling pathway. (E) MAPK signaling pathway. (F) P53 signaling pathway. (G) Neurotrophin signaling pathway. (H) T cell receptor signaling pathway. (I) Toll like receptor signaling pathway. NES, normalized enrichment score; ES, enrichment score; FDR, false discovery rate.

Figure 6. GSEA results showed differential enrichment of genes with high ZSCAN20 expression. (A) Cell cycle. (B) DNA replication. (C) WNT signaling pathway. (D) JAK STAT signaling pathway. (E) MAPK signaling pathway. (F) P53 signaling pathway. (G) Neurotrophin signaling pathway. (H) T cell receptor signaling pathway. (I) Toll like receptor signaling pathway. NES, normalized enrichment score; ES, enrichment score; FDR, false discovery rate.

Identification of key ZSCAN20-interacting genes related to cell cycle

To deeply study the hidden biological functions of ZSCAN20, we pricked off the first 500 genes in LinkedOmics through the STRING database, and then built up the PPI network of ZSCAN20 (Figure 7A). In addition, we used Cytoscape’s MCODE plug-in to build the most important modules. Then, according to k-core = 2, we calculated their connection degree and the highest average connection score (21.091) to identify the most relevant 23 central genes (Figure 7B). In addition, Supplementary Table 3 details central genes. This result initially screened 23 genes closely related to ZSCAN20. We subsequently studied these 23 genes and found that the above-mentioned genes are closely related to the cell cycle. Then, we further produced scatter plots and survival curves of ZSCAN20 and the highest-scoring genes (PRC1 and CEP55) through GEPIA. The outcomes revealed that ZSCAN20 was truly connected with PRC1 and CEP55, and the high expression of PRC1 and CEP55 led to a poor prognosis (Figure 7C7F). The above results suggest that ZSCAN20 may impact on adjusting the cell cycle.

Cell cycle related genes correlation analysis. (A) Top 500 genes in LinkedOmics were filtered into protein–protein interaction network complex. (B) The PPI network of ZSCAN20 was generated using STRING. (C, D) Draw a scatter plot to show the correlation between the ZSCAN20 and the highest score genes (PRC1 and CEP55). (E, F) Associations with overall survival and the expression of PRC1 and CEP55.

Figure 7. Cell cycle related genes correlation analysis. (A) Top 500 genes in LinkedOmics were filtered into protein–protein interaction network complex. (B) The PPI network of ZSCAN20 was generated using STRING. (C, D) Draw a scatter plot to show the correlation between the ZSCAN20 and the highest score genes (PRC1 and CEP55). (E, F) Associations with overall survival and the expression of PRC1 and CEP55.

ZSCAN20 knockdown inhibits HCC cell invasion, migration and proliferation

To explore whether ZSCAN20 has a role in regulating HCC cell invasion, migration and proliferation. Invasion and migration abilities of ZSCAN20 cells were analyzed by Transwell assay to determine the role of ZSCAN20 in HCC cells. As shown, knockdown of ZSCAN20 significantly reduced the invasive and migratory abilities of HCC cells (Figure 8A, 8B). Based on colony formation assay, knockdown of ZSCAN20 strongly lowered the proliferation ability of HCC cells (Figure 8C, 8D). Collectively, our data demonstrate that ZSCAN20 knockdown significantly inhibit invasion, migration and proliferation abilities of HCC cell.

Transwell and colony formation assays. (A, B) Migration and invasion capacity for HCC cells treated with shZSCAN20 or shNC was detected by Transwell separately. **p C, D) Colonies formed by HCCLM3 cells transfected with control shRNA or shRNA targeting ZSCAN20. The right panels are quantification of the results of the colony formation assay.

Figure 8. Transwell and colony formation assays. (A, B) Migration and invasion capacity for HCC cells treated with shZSCAN20 or shNC was detected by Transwell separately. **p < 0.01, ***p < 0.001. (C, D) Colonies formed by HCCLM3 cells transfected with control shRNA or shRNA targeting ZSCAN20. The right panels are quantification of the results of the colony formation assay.

The association of ZSCAN20 and immune infiltration in HCC

In order to understand the relationship between ZSCAN20 expression and tumor infiltrating immune cells, accordingly, we appraised the nexus between ZSCAN20 and a variety of immune cells. Via TIMER database, we knew that ZSCAN20 was related to some immune-infiltrating cells, like B cell, CD8+ T cell, CD4+ T cell, macrophage, neutrophil and dendritic cell (Figure 9A). Then, ZSCAN20 CNV was utilized to clarify the hidden mechanisms associated with different immune cell infiltration. Absence and deep loss at arm level seriously disturbed the infiltration level of B cell, CD8+ T cell, neutrophil and dendritic cell in ZSCAN20. These outcomes indicated that ZSCAN20 was closely related to the infiltration of LIHC immune cells, especially the infiltration of T cell and neutrophil (Figure 9B). Then, TISIDB was utilized to evaluate the association between immune subtypes and ZSCAN20 expression. We utilized it to explore the mRNA expression of ZSCAN20 and diverse immune subgroups in HCC, and we found that the expression of ZSCAN20 was higher in the C1 (wound healing) and C2 (IFN-γ dominant) subgroups, while the expression in the C6 (TGF-β dominant) subgroup was lower (Figure 9C). We also further studied the relationship between ZSCAN20 expression in the GEPIA database and T cell checkpoints (such as CD274, CTLA4, and PDCD1). The expression of ZSCAN20 had a significant relationship with CD274 expression (R = 0.25 P = 1.5E-06), CTLA4 (R = 0.19, P = 1.9E-04) and PDCD1 (R = 0.14, P = 6.6E-03) in HCC (Figure 9D). These outcomes prove that ZSCAN20 is intimately associated with the infiltration of immune cells.

Association between ZSCAN20 and immune cell infiltration in HCC. (A) ZSCAN20 showed a significant correlation with the infiltrating abundance of B cell, CD8+ T cell, CD4+ T cell, Macrophage, Neutrophil, and Dendritic cell using the TIMER database. (B) The relationship between ZSCAN20 CNV and immune infiltration.*p C) Distribution of ZSCAN20 expression across immune subtypes in LIHC (TISIDB). The different color plots represent the five immune subtypes (C1: wound healing; C2: IFN-gamma dominant; C3: inflammatory; C4: lymphocyte-depleted and C6: TGF-b dominant). (D) Scatter plots of the correlations between ZSCAN20 expression and CD274, CTLA4 and PDCD1 in LIHC using the GEPIA database.

Figure 9. Association between ZSCAN20 and immune cell infiltration in HCC. (A) ZSCAN20 showed a significant correlation with the infiltrating abundance of B cell, CD8+ T cell, CD4+ T cell, Macrophage, Neutrophil, and Dendritic cell using the TIMER database. (B) The relationship between ZSCAN20 CNV and immune infiltration.*p < 0.05, **p < 0.01, ***p < 0.00. (C) Distribution of ZSCAN20 expression across immune subtypes in LIHC (TISIDB). The different color plots represent the five immune subtypes (C1: wound healing; C2: IFN-gamma dominant; C3: inflammatory; C4: lymphocyte-depleted and C6: TGF-b dominant). (D) Scatter plots of the correlations between ZSCAN20 expression and CD274, CTLA4 and PDCD1 in LIHC using the GEPIA database.

The correlation among ZSCAN20 and immune markers, T cell immune markers in HCC

With the purpose of going deep on ZSCAN20 expression related to different markers of multiple immune cells. We made use of the TIMER database to further study the relation, and the genes listed in Table 2 and Figure 10A10F were used to characterize different immune cells. Tumor purity is a vital factor affecting clinical tumor biopsy immune infiltration anatomy. The expression of ZSCAN20 had a considerable positive correlation with the levels of most markers in distinguished kinds of LIHC immune cells after regulating tumor purity. In addition, we evaluated the link between ZSCAN20 and multiple T cell immune markers. We found that ZSCAN20 expression level was remarkably associated with 25 of the 26 T cell markers in LIHC after adjusting for tumor purity (Table 3). Moreover, ZSCAN20 was also closely linked to the LIHC-related chemokines including CXCL5, CXCL6 and CXCL8 (Figure 10G10I). Besides, we did the immune-related chemokines expression analysis in HCC cells. CXCL5, CXCL6 and CXCL8 were all down-regulated in ZSCAN20 shNC cells (Figure 10J10L). It could be known from the results that ZSCAN20 was positively associated with these chemokines.

Table 2. Relationship between ZSCAN20 and gene marker sets of different immune cells using the TIMER database.

DescriptionGene markersLIHC
NonePurity
CorpCorp
B cellCD190.185.65E-040.184.30E-04
CD79A0.131.55E-020.131.46E-02
T cell (general)CD3D0.214.43E-050.252.35E-06
CD3E0.201.37E-040.253.36E-06
CD2-0.161.51E-010.245.37E-06
CD8+ T cellCD8A0.221.70E-050.252.00E-06
CD8B0.171.13E-030.193.32E-04
MonocyteCD860.398.92E-150.462.89E-19
CSF1R0.341.09E-110.405.32E-15
TAMCCL20.292.00E-080.335.25E-10
CD680.291.52E-080.312.50E-09
IL100.352.34E-120.393.18E-14
M1IRF50.391.27E-140.385.35E-13
PTGS20.336.93E-110.393.26E-14
M2CD1630.303.40E-090.344.63E-11
VSIG40.336.25E-110.384.26E-13
MS4A4A0.318.00E-100.371.09E-12
NeutrophilsCEACAM80.053.01E-010.071.95E-01
ITGAM0.521.51E-270.571.13E-30
CCR70.201.36E-040.231.82E-05
Natural killer cellKIR2DL10.044.21E-010.035.80E-01
KIR2DL30.192.50E-040.201.26E-04
KIR2DL40.192.81E-040.193.08E-04
KIR3DL10.122.67E-020.122.63E-02
KIR3DL20.071.55E-010.098.17E-02
KIR3DL30.044.89E-010.036.19E-01
Dendritic cellHLA-DPB10.277.73E-080.316.42E-09
HLA-DQB10.209.16E-050.249.28E-06
HLA-DRA0.321.50E-100.363.64E-12
HLA-DPA10.317.28E-100.367.17E-12
CD1C0.146.97E-030.156.06E-03
NRP10.452.55E-200.459.17E-19
ITGAX0.406.13E-160.461.29E-19
ZSCAN20 was associated with immune cell infiltration in HCC. The scatter plots showed relationship between ZSCAN20 and various gene markers of (A) B cell, (B) T cell, (C) CD8+T cell, (D) TAM, (E) M1 macrophage and (F) M2 macrophage in LIHC. (G–I) The association between ZSCAN20 and LIHC-related chemokines; (J–L) Immune-related chemokines expression in sh-ZSCAN20 HCC cells. *p p p

Figure 10. ZSCAN20 was associated with immune cell infiltration in HCC. The scatter plots showed relationship between ZSCAN20 and various gene markers of (A) B cell, (B) T cell, (C) CD8+T cell, (D) TAM, (E) M1 macrophage and (F) M2 macrophage in LIHC. (GI) The association between ZSCAN20 and LIHC-related chemokines; (JL) Immune-related chemokines expression in sh-ZSCAN20 HCC cells. *p < 0.05, **p < 0.01, ***p < 0.001.

Table 3. Relationship between ZSCAN20 and gene marker sets of different T cells using the TIMER database.

DescriptionGene markersLIHC
NonePurity
CorpCorp
Th1 cellTBX210.152.90E-030.186.32E-04
STAT40.235.63E-060.269.94E-07
STAT10.412.76E-160.425.11E-16
TNF0.356.32E-120.395.29E-14
HAVCR20.424.77E-170.508.99E-23
Th1-like cellCXCR30.209.27E-050.232.14E-05
BHLHE400.459.79E-200.459.69E-19
CD40.303.31E-090.327.80E-10
Th2 cellSTAT60.311.13E-090.302.04E-08
STAT5A0.381.61E-140.399.88E-14
FOXP30.245.78E-160.262.22E-17
Treg cellCCR80.451.40E-380.497.97E-44
TGFB10.371.58E-130.401.31E-14
Resting Treg cellFOXP30.245.78E-160.262.22E-17
IL2RA0.381.98E-140.441.21E-17
FOXP30.245.78E-160.262.22E-17
Effector Treg cellCCR80.451.40E-380.497.97E-44
TNFRSF90.412.11E-160.452.19E-18
CX3CR10.406.37E-160.418.63E-16
Effector T cellFGFBP2-0.062.60E-01-0.054.02E-01
FCGR3A0.424.40E-170.462.33E-19
CCR70.206.47E-080.231.23E-09
Naive T cellSELL0.311.95E-180.352.29E-21
Effector memory T cellDUSP40.391.43E-140.432.94E-17
GZMK0.112.78E-020.148.94E-03
GZMA0.146.21E-030.171.34E-03
Resident memory T cellCD690.251.17E-060.301.37E-08
CXCR60.207.26E-050.261.33E-06
MYADM0.442.79E-190.452.86E-18
General memoryCCR70.206.47E-080.231.23E-09
SELL0.311.95E-180.352.29E-21

According to the results above, we have confirmed that the expression of ZSCAN20 was positively correlated with the infiltration of most immune cells in HCC, and the expression of ZSCAN20 was also associated with the poor prognosis of these tumors. Therefore, we speculated that part of the reason that ZSCAN20 expression affected the poor prognosis of HCC is immune infiltration. We performed a prognostic analysis of ZSCAN20 expression of patients in different immune cell subsets by Kaplan-Meier Plotter. We found that in enriched regulatory T cell cohorts, high expression of ZSCAN20 led to poor prognosis in HCC patients, but there was no significant association between low/high ZSCAN20 expression and HCC patient prognosis in decreased Regulatory T-cells cohorts (Figure 11A). Nonetheless, in other both enriched and decreased immune cell groups, it revealed that high expression of ZSCAN20 resulted in poor prognosis (Figure 11B11H). These results supported our speculation that high ZSCAN20 expression levels all contributed to the poor prognosis of HCC patients in part due to the level of regulatory T-cell immune infiltration. It is these consequences that show ZSCAN20 has a close association with the infiltration of HCC immune infiltrating cells, and may be associated with immune infiltration of regulatory T cells.

Kaplan-Meier survival curves according to high and low expression of ZSCAN20 in immune cell subgroups in HCC. (A–H) Correlations between ZSCAN20 expression and OS in different immune cell subgroups in HCC patients were estimated by Kaplan-Meier plotter.

Figure 11. Kaplan-Meier survival curves according to high and low expression of ZSCAN20 in immune cell subgroups in HCC. (AH) Correlations between ZSCAN20 expression and OS in different immune cell subgroups in HCC patients were estimated by Kaplan-Meier plotter.

Correlation between ZSCAN20 expression and m6A modification in HCC

Exploring whether ZSCAN20 can be modified by m6A played a significant role in the growth of HCC. We obtained the TCGA and ICGC data sets to research the relation between ZSCAN20 expression and the expression of 21 m6a-related genes in LIHC. The outcomes showed that the expression of ZSCAN20 had a positive relationship with many m6a-related genes (Figure 12A, 12B) (P < 0.05). Meanwhile, With the purpose of ascertaining whether there was a difference between high and low expression of m6A modification expressed by ZSCNA20 in LIHC, we tried to constrain the differential expression of m6A-related genes between high and low ZSCAN20 expression groups, and it was the outcomes that showed ZSCAN20 high expression groups were higher than those in the low expression groups (P < 0.05) (Figure 12C). We selected genes to draw a Venn diagram (correlation coefficient of TCGA and ICGC > 0.55), which showed the correlation and differential expression of genes, including METL3, HNRNPA2B1, RBM15B, RBMX and YTHDF1 (Figure 12D). Then, we produced scatter plots of ZSCAN20 and m6A related genes, which showed that ZSCAN20 was momentously positively associated with METL3, HNRNPA2B1, RBM15B, RBMX and YTHDF1 (Figure 12E). We further used GEPIA to study the prognosis of related genes, and the outcomes of GEPIA displayed the high expression of METL3, HNRNPA2B1, RBM15B, RBMX and YTHDF1 were closely related to the poor prognosis of LIHC (P < 0.05) (Figure 12F). These results indicate that ZSCAN20 may have an intimate relationship with the m6A modification of LIHC, may affecting the occurrence and development of tumors together with m6A-related genes.

Correlations of ZSCAN20 expression with m6A related genes in HCC. (A, B) TCGA LIHC and ICGC data sets were used to analyze the correlation between the ZSCAN20 and the m6A related genes expression in LIHC. (C) Draw a scatter plot to show the correlation between the ZSCAN20 and the m6A related genes expression, including METTL3, HNRNPA2B1, RBM15B, RBMX and YTHDF1. (D) The differential expression of m6A related genes between high and low ZSCNA20 expression groups in LIHC tumor samples. (E) Venn diagram showed both expression correlation and differential expression of genes. (F) Kaplan-Meier curve of METTL3, HNRNPA2B1, RBM15B, RBMX, YTHDF1. *p

Figure 12. Correlations of ZSCAN20 expression with m6A related genes in HCC. (A, B) TCGA LIHC and ICGC data sets were used to analyze the correlation between the ZSCAN20 and the m6A related genes expression in LIHC. (C) Draw a scatter plot to show the correlation between the ZSCAN20 and the m6A related genes expression, including METTL3, HNRNPA2B1, RBM15B, RBMX and YTHDF1. (D) The differential expression of m6A related genes between high and low ZSCNA20 expression groups in LIHC tumor samples. (E) Venn diagram showed both expression correlation and differential expression of genes. (F) Kaplan-Meier curve of METTL3, HNRNPA2B1, RBM15B, RBMX, YTHDF1. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001. ns, not significant.

Discussion

As the fifth most common cancer globally, HCC is a malignant tumor with high mortality. Patients with HCC often lack evident clinical symptoms in the first instance. Despite the discovery and treatment, patients may still have tumor progression or recurrence. In order to study new tumor treatment methods, such as the use of checkpoint inhibitors for immunotherapy [4143] except for AFP [44], new biomarkers have not yet entered daily practice. Therefore, research to find a new and meaningful biomarker to assist HCC detection and treatment is imminent. Here we applied various bioinformatics analysis for the purpose of investigating ZSCAN20’s expression and prognosis in HCC. At the same time, study its correlation with cell cycle, immune infiltration and m6A.

First, we found that ZSCAN20 was highly expressed in HCC through public databases and HCC clinical disease samples. The UALCAN results suggested that ZSCAN20 was closely connected to the clinicopathological features. The Kaplan-Meier Plotter survival analysis showed high ZSCAN20 expression was linked to low survival rate. And Cox analysis suggested that ZSCAN20 can be viewed as an independent predictor. Studies have shown that epigenetic changes may increase the risk of cancer by altering the expression of oncogenes [45]. DNA methylation is one type of epigenetic change that can be considered as a potential marker of cancer development and progression [46]. In HCC, according to our results, ZSCAN20 was hypermethylated in HCC. Meanwhile, the hypermethylation of ZSCAN20 was also associated with low survival rates. These results indicated that in HCC, ZSCAN20 was highly expressed and also hypermethylated. Both its high expression and hypermethylation indicated a poor prognosis.

In order to confirm our hypothesis, we further studied the biological function of ZSCAN20 in HCC. In order to analyze the biological pathways that ZSCAN20 can regulate in HCC, we used LinkedOmics to investigate the biological functions of ZSCAN20 through co-expressed genes. GO and KEGG results suggested that ZSCAN20 was mainly associated with cell cycle, DNA replication and immune infiltration. In addition, GSEA results showed that ZSCAN20 was related to cancer pathways such as the famous TP53 signaling pathway, MAPK signaling pathway. At the same time, GSEA results also showed that ZSCAN20 was related to cell cycle, Wnt signaling pathway and other famous cell cycles. Wnt signaling pathway vitally participates in cell cycle regulation, physiological behavior, disease and tumorigenesis, which is related to T cell factors and lymphatic enhancement. Factors combine to form complexes which affect cell proliferation, apoptosis, and migration [47, 48]. We still need to further study whether ZSCAN20 affected the occurrence and development of cancer over Wnt signaling pathway in the cell cycle mechanism.

In order to further study its potential biological functions, the PPI results showed that the first 23 genes associated with ZSCAN20 were correlated with the cell cycle. It was well known that the anomalous expression of cell cycle related proteins may affect cell growth and chromosomal instability, causing cancer progression. We selected the two co-expressed genes (PRC1 and CEP55) with the highest cluster scores and found that they had a strong correlation with ZSCAN20. Subsequently, we conducted survival analysis of these two genes, and the results indicated their high expression can cause poor prognosis. Studies have pointed out that the dysregulation of the protein regulator 1 (PRC1) of cytokinesis led to cytokinesis defects in the cell cycle, which in turn promoted chromosomal instability (CIN), thereby promoting tumor heterogeneity and cancer progression [49]. In addition, the absence of CEP55 leads to mitotic arrest and mitotic cell death in the cell cycle, and CEP55 has been proven to be an anti-mitotic drug targeting the mitotic mechanism of the cell cycle to eliminate cancer cells [50]. This implied that ZSCAN20 may get involved in HCC through its association with cell cycle-related signal pathways.

The above studies confirmed that ZSCAN20 was highly expressed in HCC and was closely related to poor prognosis. At the same time, the abilities of invasion, migration and proliferation are the main factors to evaluate the progression of malignant tumor [51, 52]. To investigate the effect of ZSCAN20 on malignant pathological processes in HCC. Through Transwell and colony formation assays, we found that interfering with the expression of ZSCAN20 could significantly inhibit the invasion, migration and proliferation abilities of HCC cells.

As one of the most major organs, the liver has a variety of immune functions. HCC is closely associated with inflammation, and immune tolerance and escape have an important impact on inflammation [53]. Most cancers develop immune tolerance and escape by impairing the function of immune cells or releasing inhibitory cytokines [54]. Recent studies have shown that tumor infiltrating immune cells (TIIC) have a key influence regulating and controlling tumor progression. Moreover, TIICs are preferentially enriched in HCC and may be cloned and expanded [55], and the accumulation of TIICs is associated with poor prognosis of HCC [56]. Functional network analysis results showed the biological function of ZSCAN20 in HCC was directly associated with the immune response. Thus, we explored whether the expression of ZSCAN20 can affect HCC immune cells in the TIMER database. We found ZSCAN20 expression was extremely related to tumor purity. Among the six immune cells, the level of ZSCAN20 was positively correlated with the infiltration of CD8+T cells, CD4+T cells, macrophages, neutrophils and dendritic cells. Significant associations between ZSCAN20 and diverse immune cell marker sets were discovered in HCC. Studies have reported that insufficient crosstalk between dendritic cells (DC) and T cells can cause HCC tumor tolerance [57]. In HCC patients, Treg binds to lots of activated CD8+ cytotoxic cells and is related to OS and DFS [58]. Furthermore, Tregs have been found to have profound roles in regulating immune responses and maintaining self-tolerance and immune homeostasis. While Tregs play an important role in preventing autoimmune diseases, Tregs can regulate the adaptive immune system by regulating CD4 Th cells, CD8 T cells and B cells. At the same time, Tregs also regulate innate immune system cells (including dendritic cells, macrophages, neutrophils, γδ T cells, natural killer and innate lymphoid cells) [59]. The study also found that Tregs also played an important role in HCC [60]. Here ZSCAN20 was closely related to various subtypes of T cells, including memory T, effector T, Th1-like, effector Treg and exhausted T cell. Meanwhile, the high expression of ZSCAN20 was correlated with poor prognosis in the enriched Regulatory T-cell. Recently, Immune Checkpoint Inhibitor (ICI) is a swiftly developing immunotherapy method for liver cancer [61]. ICI can activate tumor-specific T cells to fight tumors [62]. Here we found that ZSCAN20 was significantly related to T cell checkpoints (CD274, CTLA4 and PDCD1). Additionally, it was found through TISIDB that ZSCAN20 was positively correlated with CXCL5, CXCL6 and CXCL8, and interference with ZSCAN20 can significantly reduce the expression of these chemokines. Our findings suggested that in HCC ZSCAN20 may exert biological functions through immune cells, and ZSCAN20 may regulate the development and prognosis of HCC through Regulatory T-cell.

Methylation occurs at the N position of the sixth adenylate, that is, N6-methyladenosine, referred to as m6A [63], which is common and numerous post-transcriptional modification on mRNA and lncRNA [64]. M6A is used in various biological processes and diseases. Disease pathogenesis and tumorigenesis and development process play a pivotal role [65]. Studies have shown that HNRNPA2B1, as an m6A reader, mediates alternative splicing of target RNA and enhances primary miRNA processing [66], which is highly expressed in liver cancer [67]. RBMX is a m6A regulatory factor [68], that is upregulated in HCC, and RBMX expression is positively related to the viability and proliferation of HCC cells, promoting the development of HCC and sorafenib resistance [69]. The results of our study showed that there was a positive correlation between ZSCAN20 and the expression of various m6A-related genes in TCGA and ICGC, and it was found that METL3, HNRNPA2B1, RBM15B, RBMX and YTHDF1 have the most obvious relationship with ZSCAN20. Further research we also found that when ZSCAN20 was highly expressed and lowly expressed, and most of the m6A-related genes showed differential expression. Also, by constructing a Kaplan-Meier curve, it was found that the high expression group of METTL3, HNRNPA2B1, RBM15B, RBMX and YTHDF1 had a shorter survival time. In summary, our research results indicated that the mechanism of ZSCAN20’s role in HCC may be related to m6A modification. ZSCAN20 may affect HCC mRNA methylation level through its association with METTL3, HNRNPA2B1, RBM15B, RBMX, YTHDF1, and ultimately affected the progress of HCC.

This study has some limitations. Firstly, the data in the database is constantly updated and may affect the results obtained on the results website. We need to collect more clinical data for verification. Secondly, our study only focused on HCC and normal tissues. Jiang et al. found that Hsa_circ_0028502 and hsa_circ_0076251 distinguish HCC tissues, liver cirrhosis (LC) and chronic hepatitis (CH) [70]. We can collect more clinical data to further explore the development of HCC, analyzing the expression of ZSCAN20 in HCC tissues, LC tissues and CH tissues and whether ZSCAN20 can be used as a biomarker to distinguish HCC tissues, LC tissues and CH tissues. Finally, Tian et al. found the prognostic value of TPM1–4 as a risk model in HCC [71]. Although ZSCAN20 can be used as an independent prognostic factor for HCC, we can also study other genes in ZSCAN20’s family and combine them into a risk model in future studies. The accuracy of prediction may be greatly improved, and it may also be able to predict different types of HCC patients.

On a final note, using bioinformatics methods as well as experimental verification, this study has proved ZSCAN20 is upregulated in HCC. ZSCAN20 expression is correlated with the clinicopathological characteristics and prognosis of HCC patients. In addition, ZSCAN20 may have an influence on HCC through cell cycle and immune infiltration, m6A modification and other processes. Interference with ZSCAN20 could inhibit invasion, migration and proliferation abilities of HCC cell. Collectively, ZSCAN20 can be used as a biomarker for the diagnosis, treatment and prognosis of HCC.

Supplementary Materials

Supplementary Tables

Abbreviations

HCC: Hepatocellular carcinoma; ZSCAN20: Zinc Finger And SCAN Domain Containing 20; TIMER: Tumor Immune Estimation Resource; TCGA: The Cancer Genome Atlas; ICGC: International Cancer Genome Consortium; GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; GSEA: Gene Set Enrichment Analysis; PPI: protein-protein interaction; TIIC: Tumor infiltrating immune cells; AFP: Alpha Fetoprotein; CNV: copy number variation.

Author Contributions

Da Huang and Weifan Jiang contributed to conception and design of the study. Jiayu Fang and Yike jiang organized the database. Jingying Pan, Yun Ke, Jing Zhen performed statistical analysis. Jiayu Fang, Yike Jiang, Xuanrui Zhou and Minqin Zhou wrote sections of the manuscript. Jiayu Fang designed and performed the IHC, Western Blot and chemokines experiments. Yike Jiang designed and performed Transwell assay and colony formation assay. All authors contributed to manuscript revision, read, and approved the submitted version.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Ethical Statement and Consent

This study has been approved by the Ethics Committee of the Second Affiliated Hospital of Nanchang University. All 40 patients that provided samples for this study endorsed a acknowledgeable consent form

Funding

This study was supported by grants from the Project of Department of Education of Jiangxi Procince of China (no.GJJ200244).

References

  • 1. European Association for the Study of the Liver. Electronic address: easloffice@easloffice.eu, and European Association for the Study of the Liver. EASL Clinical Practice Guidelines: Management of hepatocellular carcinoma. J Hepatol. 2018; 69:182–236. https://doi.org/10.1016/j.jhep.2018.03.019 [PubMed]
  • 2. Asafo-Agyei KO, Samant H. Hepatocellular Carcinoma. In: StatPearls. Treasure Island (FL): StatPearls Publishing. 2022. [PubMed]
  • 3. Llovet JM, Zucman-Rossi J, Pikarsky E, Sangro B, Schwartz M, Sherman M, Gores G. Hepatocellular carcinoma. Nat Rev Dis Primers. 2016; 2:16018. https://doi.org/10.1038/nrdp.2016.18 [PubMed]
  • 4. Villanueva A. Hepatocellular Carcinoma. N Engl J Med. 2019; 380:1450–62. https://doi.org/10.1056/NEJMra1713263 [PubMed]
  • 5. Yang JD, Hainaut P, Gores GJ, Amadou A, Plymoth A, Roberts LR. A global view of hepatocellular carcinoma: trends, risk, prevention and management. Nat Rev Gastroenterol Hepatol. 2019; 16:589–604. https://doi.org/10.1038/s41575-019-0186-y [PubMed]
  • 6. Duan H, Zhao G, Xu B, Hu S, Li J. Maternal Serum PLGF, PAPPA, β-hCG and AFP Levels in Early Second Trimester as Predictors of Preeclampsia. Clin Lab. 2017; 63:921–5. https://doi.org/10.7754/Clin.Lab.2016.161103 [PubMed]
  • 7. Falzone L, Salomone S, Libra M. Evolution of Cancer Pharmacological Treatments at the Turn of the Third Millennium. Front Pharmacol. 2018; 9:1300. https://doi.org/10.3389/fphar.2018.01300 [PubMed]
  • 8. Sprinzl MF, Galle PR. Current progress in immunotherapy of hepatocellular carcinoma. J Hepatol. 2017; 66:482–4. https://doi.org/10.1016/j.jhep.2016.12.009 [PubMed]
  • 9. Topalian SL, Drake CG, Pardoll DM. Immune checkpoint blockade: a common denominator approach to cancer therapy. Cancer Cell. 2015; 27:450–61. https://doi.org/10.1016/j.ccell.2015.03.001 [PubMed]
  • 10. Rotte A, Jin JY, Lemaire V. Mechanistic overview of immune checkpoints to support the rational design of their combinations in cancer immunotherapy. Ann Oncol. 2018; 29:71–83. https://doi.org/10.1093/annonc/mdx686 [PubMed]
  • 11. Klug A, Schwabe JW. Protein motifs 5. Zinc fingers. FASEB J. 1995; 9:597–604. [PubMed]
  • 12. Huang M, Chen Y, Han D, Lei Z, Chu X. Role of the zinc finger and SCAN domain-containing transcription factors in cancer. Am J Cancer Res. 2019; 9:816–36. [PubMed]
  • 13. Liu G, Jiang S, Wang C, Jiang W, Liu Z, Liu C, Saiyin H, Yang X, Shen S, Jiang D, Zhou P, Han D, Hu X, et al. Zinc finger transcription factor 191, directly binding to β-catenin promoter, promotes cell proliferation of hepatocellular carcinoma. Hepatology. 2012; 55:1830–9. https://doi.org/10.1002/hep.25564 [PubMed]
  • 14. Yao Z, Luo J, Hu K, Lin J, Huang H, Wang Q, Zhang P, Xiong Z, He C, Huang Z, Liu B, Yang Y. ZKSCAN1 gene and its related circular RNA (circZKSCAN1) both inhibit hepatocellular carcinoma cell growth, migration, and invasion but through different signaling pathways. Mol Oncol. 2017; 11:422–37. https://doi.org/10.1002/1878-0261.12045 [PubMed]
  • 15. Kim CW, Roh SA, Tak KH, Koh BM, Ha YJ, Cho DH, Kim SY, Kim YS, Kim JC. ZKSCAN3 Facilitates Liver Metastasis of Colorectal Cancer Associated with CEA-expressing Tumor. Anticancer Res. 2016; 36:2397–406. [PubMed]
  • 16. Zhou J, Chen WR, Yang LC, Wang J, Sun JY, Zhang WW, He ZY, Wu SG. KIF11 Functions as an Oncogene and Is Associated with Poor Outcomes from Breast Cancer. Cancer Res Treat. 2019; 51:1207–21. https://doi.org/10.4143/crt.2018.460 [PubMed]
  • 17. Mosquera Orgueira A, Ferreiro Ferro R, Díaz Arias JÁ, Aliste Santos C, Antelo Rodríguez B, Bao Pérez L, Alonso Vence N, Bendaña López Á, Abuin Blanco A, Melero Valentín P, Peleteiro Raindo AR, Cid López M, Pérez Encinas MM, et al. Detection of new drivers of frequent B-cell lymphoid neoplasms using an integrated analysis of whole genomes. PLoS One. 2021; 16:e0248886. https://doi.org/10.1371/journal.pone.0248886 [PubMed]
  • 18. Liu J, Zhang S, Dai W, Xie C, Li JC. A Comprehensive Prognostic and Immune Analysis of SLC41A3 in Pan-Cancer. Front Oncol. 2021; 10:586414. https://doi.org/10.3389/fonc.2020.586414 [PubMed]
  • 19. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 2017; 77:e108–10. https://doi.org/10.1158/0008-5472.CAN-17-0307 [PubMed]
  • 20. Aran D, Sirota M, Butte AJ. Systematic pan-cancer analysis of tumour purity. Nat Commun. 2015; 6:8971. https://doi.org/10.1038/ncomms9971 [PubMed]
  • 21. Li B, Severson E, Pignon JC, Zhao H, Li T, Novak J, Jiang P, Shen H, Aster JC, Rodig S, Signoretti S, Liu JS, Liu XS. Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol. 2016; 17:174. https://doi.org/10.1186/s13059-016-1028-7 [PubMed]
  • 22. Zhang Y, Yuan Y, Liang P, Zhang Z, Guo X, Xia L, Zhao Y, Shu XS, Sun S, Ying Y, Cheng Y. Overexpression of a novel candidate oncogene KIF14 correlates with tumor progression and poor prognosis in prostate cancer. Oncotarget. 2017; 8:45459–69. https://doi.org/10.18632/oncotarget.17564 [PubMed]
  • 23. Kirby ED, Kuwahara AA, Messer RL, Wyss-Coray T. Adult hippocampal neural stem and progenitor cells regulate the neurogenic niche by secreting VEGF. Proc Natl Acad Sci USA. 2015; 112:4128–33. https://doi.org/10.1073/pnas.1422448112 [PubMed]
  • 24. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017; 45:W98–102. https://doi.org/10.1093/nar/gkx247 [PubMed]
  • 25. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi BVSK, Varambally S. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia. 2017; 19:649–58. https://doi.org/10.1016/j.neo.2017.05.002 [PubMed]
  • 26. Liao Y, Xiao H, Cheng M, Fan X. Bioinformatics Analysis Reveals Biomarkers With Cancer Stem Cell Characteristics in Lung Squamous Cell Carcinoma. Front Genet. 2020; 11:427. https://doi.org/10.3389/fgene.2020.00427 [PubMed]
  • 27. Kadomoto S, Izumi K, Mizokami A. The CCL20-CCR6 Axis in Cancer Progression. Int J Mol Sci. 2020; 21:5186. https://doi.org/10.3390/ijms21155186 [PubMed]
  • 28. 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:277–88. https://doi.org/10.2217/epi-2017-0118 [PubMed]
  • 29. Ru B, Wong CN, Tong Y, Zhong JY, Zhong SSW, Wu WC, Chu KC, Wong CY, Lau CY, Chen I, Chan NW, Zhang J. TISIDB: an integrated repository portal for tumor-immune system interactions. Bioinformatics. 2019; 35:4200–2. https://doi.org/10.1093/bioinformatics/btz210 [PubMed]
  • 30. Yuan Q, Sun N, Zheng J, Wang Y, Yan X, Mai W, Liao Y, Chen X. Prognostic and Immunological Role of FUN14 Domain Containing 1 in Pan-Cancer: Friend or Foe? Front Oncol. 2020; 9:1502. https://doi.org/10.3389/fonc.2019.01502 [PubMed]
  • 31. Peng WX, Huang JG, Yang L, Gong AH, Mo YY. Linc-RoR promotes MAPK/ERK signaling and confers estrogen-independent growth of breast cancer. Mol Cancer. 2017; 16:161. https://doi.org/10.1186/s12943-017-0727-3 [PubMed]
  • 32. Yan T, Tian D, Chen J, Tan Y, Cheng Y, Ye L, Deng G, Liu B, Yuan F, Zhang S, Cai L, Chen Q. FCGBP Is a Prognostic Biomarker and Associated With Immune Infiltration in Glioma. Front Oncol. 2022; 11:769033. https://doi.org/10.3389/fonc.2021.769033 [PubMed]
  • 33. 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–63. https://doi.org/10.1093/nar/gkx1090 [PubMed]
  • 34. Wang J, Vasaikar S, Shi Z, Greer M, Zhang B. WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 2017; 45:W130–7. https://doi.org/10.1093/nar/gkx356 [PubMed]
  • 35. Liberzon A, Subramanian A, Pinchback R, Thorvaldsdóttir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011; 27:1739–40. https://doi.org/10.1093/bioinformatics/btr260 [PubMed]
  • 36. Wang Z, Wang Y, Peng M, Yi L. UBASH3B Is a Novel Prognostic Biomarker and Correlated With Immune Infiltrates in Prostate Cancer. Front Oncol. 2020; 9:1517. https://doi.org/10.3389/fonc.2019.01517 [PubMed]
  • 37. Li Q, Pan Y, Cao Z, Zhao S. Comprehensive Analysis of Prognostic Value and Immune Infiltration of Chromobox Family Members in Colorectal Cancer. Front Oncol. 2020; 10:582667. https://doi.org/10.3389/fonc.2020.582667 [PubMed]
  • 38. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, Kuhn M, Bork P, Jensen LJ, von Mering C. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015; 43:D447–52. https://doi.org/10.1093/nar/gku1003 [PubMed]
  • 39. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003; 13:2498–504. https://doi.org/10.1101/gr.1239303 [PubMed]
  • 40. Ma X, Wang N, Chen K, Zhang C. Oncosuppressive role of MicroRNA-205-3p in gastric cancer through inhibition of proliferation and induction of senescence: Oncosuppressive role of MicroRNA-205 in gastric cancer. Transl Oncol. 2021; 14:101199. https://doi.org/10.1016/j.tranon.2021.101199 [PubMed]
  • 41. El-Khoueiry AB, Sangro B, Yau T, Crocenzi TS, Kudo M, Hsu C, Kim TY, Choo SP, Trojan J, Welling TH Rd, Meyer T, Kang YK, Yeo W, et al. Nivolumab in patients with advanced hepatocellular carcinoma (CheckMate 040): an open-label, non-comparative, phase 1/2 dose escalation and expansion trial. Lancet. 2017; 389:2492–502. https://doi.org/10.1016/S0140-6736(17)31046-2 [PubMed]
  • 42. Yau T, Hsu C, Kim TY, Choo SP, Kang YK, Hou MM, Numata K, Yeo W, Chopra A, Ikeda M, Kuromatsu R, Moriguchi M, Chao Y, et al. Nivolumab in advanced hepatocellular carcinoma: Sorafenib-experienced Asian cohort analysis. J Hepatol. 2019; 71:543–52. https://doi.org/10.1016/j.jhep.2019.05.014 [PubMed]
  • 43. Zhu AX, Finn RS, Edeline J, Cattan S, Ogasawara S, Palmer D, Verslype C, Zagonel V, Fartoux L, Vogel A, Sarker D, Verset G, Chan SL, et al, and KEYNOTE-224 investigators. Pembrolizumab in patients with advanced hepatocellular carcinoma previously treated with sorafenib (KEYNOTE-224): a non-randomised, open-label phase 2 trial. Lancet Oncol. 2018; 19:940–52. https://doi.org/10.1016/S1470-2045(18)30351-6 [PubMed]
  • 44. Zhu AX, Kang YK, Yen CJ, Finn RS, Galle PR, Llovet JM, Assenat E, Brandi G, Pracht M, Lim HY, Rau KM, Motomura K, Ohno I, et al, and REACH-2 study investigators. Ramucirumab after sorafenib in patients with advanced hepatocellular carcinoma and increased α-fetoprotein concentrations (REACH-2): a randomised, double-blind, placebo-controlled, phase 3 trial. Lancet Oncol. 2019; 20:282–96. https://doi.org/10.1016/S1470-2045(18)30937-9 [PubMed]
  • 45. Jia Y, Zeng Z, Li Y, Li Z, Jin L, Zhang Z, Wang L, Wang FS. Impaired function of CD4+ T follicular helper (Tfh) cells associated with hepatocellular carcinoma progression. PLoS One. 2015; 10:e0117458. https://doi.org/10.1371/journal.pone.0117458 [PubMed]
  • 46. Nebbioso A, Tambaro FP, Dell’Aversana C, Altucci L. Cancer epigenetics: Moving forward. PLoS Genet. 2018; 14:e1007362. https://doi.org/10.1371/journal.pgen.1007362 [PubMed]
  • 47. Peng YY, He YH, Chen C, Xu T, Li L, Ni MM, Meng XM, Huang C, Li J. NLRC5 regulates cell proliferation, migration and invasion in hepatocellular carcinoma by targeting the Wnt/β-catenin signaling pathway. Cancer Lett. 2016; 376:10–21. https://doi.org/10.1016/j.canlet.2016.03.006 [PubMed]
  • 48. Zhao ZQ, Wu XJ, Cheng YH, Zhou YF, Ma XM, Zhang J, Heng XY, Feng F. TROAP regulates cell cycle and promotes tumor progression through Wnt/β-Catenin signaling pathway in glioma cells. CNS Neurosci Ther. 2021; 27:1064–76. https://doi.org/10.1111/cns.13688 [PubMed]
  • 49. Chen J, Rajasekaran M, Xia H, Zhang X, Kong SN, Sekar K, Seshachalam VP, Deivasigamani A, Goh BK, Ooi LL, Hong W, Hui KM. The microtubule-associated protein PRC1 promotes early recurrence of hepatocellular carcinoma in association with the Wnt/β-catenin signalling pathway. Gut. 2016; 65:1522–34. https://doi.org/10.1136/gutjnl-2015-310625 [PubMed]
  • 50. Zhou W, Yang Y, Xia J, Wang H, Salama ME, Xiong W, Xu H, Shetty S, Chen T, Zeng Z, Shi L, Zangari M, Miles R, et al. NEK2 induces drug resistance mainly through activation of efflux drug pumps and is associated with poor prognosis in myeloma and other cancers. Cancer Cell. 2013; 23:48–62. https://doi.org/10.1016/j.ccr.2012.12.001 [PubMed]
  • 51. Wang Z, Su G, Dai Z, Meng M, Zhang H, Fan F, Liu Z, Zhang L, Weygant N, He F, Fang N, Zhang L, Cheng Q. Circadian clock genes promote glioma progression by affecting tumour immune infiltration and tumour cell proliferation. Cell Prolif. 2021; 54:e12988. https://doi.org/10.1111/cpr.12988 [PubMed]
  • 52. Lee EY, Kim M, Choi BK, Kim DH, Choi I, You HJ. TJP1 Contributes to Tumor Progression through Supporting Cell-Cell Aggregation and Communicating with Tumor Microenvironment in Leiomyosarcoma. Mol Cells. 2021; 44:784–94. https://doi.org/10.14348/molcells.2021.0130 [PubMed]
  • 53. Shlomai A, de Jong YP, Rice CM. Virus associated malignancies: the role of viral hepatitis in hepatocellular carcinoma. Semin Cancer Biol. 2014; 26:78–88. https://doi.org/10.1016/j.semcancer.2014.01.004 [PubMed]
  • 54. Dienstmann R, Vermeulen L, Guinney J, Kopetz S, Tejpar S, Tabernero J. Consensus molecular subtypes and the evolution of precision medicine in colorectal cancer. Nat Rev Cancer. 2017; 17:79–92. https://doi.org/10.1038/nrc.2016.126 [PubMed]
  • 55. Zheng C, Zheng L, Yoo JK, Guo H, Zhang Y, Guo X, Kang B, Hu R, Huang JY, Zhang Q, Liu Z, Dong M, Hu X, et al. Landscape of Infiltrating T Cells in Liver Cancer Revealed by Single-Cell Sequencing. Cell. 2017; 169:1342–56.e16. https://doi.org/10.1016/j.cell.2017.05.035 [PubMed]
  • 56. Sun H, Liu L, Huang Q, Liu H, Huang M, Wang J, Wen H, Lin R, Qu K, Li K, Wei H, Xiao W, Sun R, et al. Accumulation of Tumor-Infiltrating CD49a+ NK Cells Correlates with Poor Prognosis for Human Hepatocellular Carcinoma. Cancer Immunol Res. 2019; 7:1535–46. https://doi.org/10.1158/2326-6066.CIR-18-0757 [PubMed]
  • 57. Lurje I, Hammerich L, Tacke F. Dendritic Cell and T Cell Crosstalk in Liver Fibrogenesis and Hepatocarcinogenesis: Implications for Prevention and Therapy of Liver Cancer. Int J Mol Sci. 2020; 21:7378. https://doi.org/10.3390/ijms21197378 [PubMed]
  • 58. Gao Q, Qiu SJ, Fan J, Zhou J, Wang XY, Xiao YS, Xu Y, Li YW, Tang ZY. Intratumoral balance of regulatory and cytotoxic T cells is associated with prognosis of hepatocellular carcinoma after resection. J Clin Oncol. 2007; 25:2586–93. https://doi.org/10.1200/JCO.2006.09.4565 [PubMed]
  • 59. Okeke EB, Uzonna JE. The Pivotal Role of Regulatory T Cells in the Regulation of Innate Immune Cells. Front Immunol. 2019; 10:680. https://doi.org/10.3389/fimmu.2019.00680 [PubMed]
  • 60. Zhou YF, Song SS, Tian MX, Tang Z, Wang H, Fang Y, Qu WF, Jiang XF, Tao CY, Huang R, Zhou PY, Zhu SG, Zhou J, et al. Cystathionine β-synthase mediated PRRX2/IL-6/STAT3 inactivation suppresses Tregs infiltration and induces apoptosis to inhibit HCC carcinogenesis. J Immunother Cancer. 2021; 9:e003031. https://doi.org/10.1136/jitc-2021-003031 [PubMed]
  • 61. Cheng AL, Hsu C, Chan SL, Choo SP, Kudo M. Challenges of combination therapy with immune checkpoint inhibitors for hepatocellular carcinoma. J Hepatol. 2020; 72:307–19. https://doi.org/10.1016/j.jhep.2019.09.025 [PubMed]
  • 62. Presti M, Westergaard MCW, Draghi A, Chamberlain CA, Gokuldass A, Svane IM, Donia M. The effects of targeted immune-regulatory strategies on tumor-specific T-cell responses in vitro. Cancer Immunol Immunother. 2021; 70:1771–6. https://doi.org/10.1007/s00262-020-02760-z [PubMed]
  • 63. Chen X, Sun YZ, Liu H, Zhang L, Li JQ, Meng J. RNA methylation and diseases: experimental results, databases, Web servers and computational models. Brief Bioinform. 2019; 20:896–917. https://doi.org/10.1093/bib/bbx142 [PubMed]
  • 64. Meyer KD, Jaffrey SR. Rethinking m6A Readers, Writers, and Erasers. Annu Rev Cell Dev Biol. 2017; 33:319–42. https://doi.org/10.1146/annurev-cellbio-100616-060758 [PubMed]
  • 65. Winkler R, Gillis E, Lasman L, Safra M, Geula S, Soyris C, Nachshon A, Tai-Schmiedel J, Friedman N, Le-Trilling VTK, Trilling M, Mandelboim M, Hanna JH, et al. m6A modification controls the innate immune response to infection by targeting type I interferons. Nat Immunol. 2019; 20:173–82. https://doi.org/10.1038/s41590-018-0275-z [PubMed]
  • 66. Li Y, Sun N, Lu Z, Sun S, Huang J, Chen Z, He J. Prognostic alternative mRNA splicing signature in non-small cell lung cancer. Cancer Lett. 2017; 393:40–51. https://doi.org/10.1016/j.canlet.2017.02.016 [PubMed]
  • 67. Cui H, Wu F, Sun Y, Fan G, Wang Q. Up-regulation and subcellular localization of hnRNP A2/B1 in the development of hepatocellular carcinoma. BMC Cancer. 2010; 10:356. https://doi.org/10.1186/1471-2407-10-356 [PubMed]
  • 68. Elliott DJ, Dalgliesh C, Hysenaj G, Ehrmann I. RBMX family proteins connect the fields of nuclear RNA processing, disease and sex chromosome biology. Int J Biochem Cell Biol. 2019; 108:1–6. https://doi.org/10.1016/j.biocel.2018.12.014 [PubMed]
  • 69. Song Y, He S, Ma X, Zhang M, Zhuang J, Wang G, Ye Y, Xia W. RBMX contributes to hepatocellular carcinoma progression and sorafenib resistance by specifically binding and stabilizing BLACAT1. Am J Cancer Res. 2020; 10:3644–65. [PubMed]
  • 70. Jiang Z, Shen L, Wang S, Wu S, Hu Y, Guo J, Fu L. Hsa_circ_0028502 and hsa_circ_0076251 are potential novel biomarkers for hepatocellular carcinoma. Cancer Med. 2019; 8:7278–87. https://doi.org/10.1002/cam4.2584 [PubMed]
  • 71. Tian Z, Zhao J, Wang Y. The prognostic value of TPM1-4 in hepatocellular carcinoma. Cancer Med. 2022; 11:433–46. https://doi.org/10.1002/cam4.4453 [PubMed]