Skip to main content

Comprehensive analysis of ZNF family genes in prognosis, immunity, and treatment of esophageal cancer

Abstract

Background

As a common malignant tumor, esophageal carcinoma (ESCA) has a low early diagnosis rate and poor prognosis. This study aimed to construct the prognostic features composed of ZNF family genes to effectively predict the prognosis of ESCA patients.

Methods

The mRNA expression matrix and clinical data were downloaded from TCGA and GEO database. Using univariate Cox analysis, lasso regression and multivariate Cox analysis, we screened six prognosis-related ZNF family genes to construct the prognostic model. We then used Kaplan-Meier plot, time-dependent receiver operating characteristic (ROC), multivariable Cox regression analysis of clinical information, and nomogram to evaluate the prognostic value within and across sets, separately and combined. We also validated the prognostic value of the six-gene signature using GSE53624 dataset. The different immune status was observed in the single sample Gene Set Enrichment Analysis (ssGSEA). Finally, real-time quantitative PCR was used to detect the expression of six prognostic ZNF genes in twelve pairs of ESCA and adjacent normal tissues.

Results

A six prognosis-related ZNF family genes model consisted of ZNF91, ZNF586, ZNF502, ZNF865, ZNF106 and ZNF225 was identified. Multivariable Cox regression analysis revealed that six prognosis-related ZNF family genes were independent prognostic factors for overall survival of ESCA patients in TCGA and GSE53624. Further, a prognostic nomogram including the riskScore, age, gender, T, stage was constructed, and TCGA/GSE53624-based calibration plots indicated its excellent predictive performance. Drug Sensitivity and ssGSEA analysis showed that the six genes model was closely related to immune cells infiltration and could be used as a potential predictor of chemotherapy sensitivity.

Conclusion

We identified six prognosis-related ZNF family genes model of ESCA, which provide evidence for individualized prevention and treatment.

Peer Review reports

Introduction

Esophageal cancer (ESCA) ranks seventh for incidence (604,000) and sixth for mortality (544,000), which means it is responsible for one in eighteen cancer deaths [1]. Traditional prognostic methods, such as histopathology and tumor staging systems, are of limited use and early detection remains a difficult goal [2]. Due to the lack of specific methods for early diagnosis and treatment, ESCA patients’ five-year survival rate is remains dismal [3]. The poor outcome urges to identify robust biomarkers for predicting the prognosis of ESCA patients.

Zinc finger protein encoded by nearly 5% of the human genome is the largest family of transcription factor proteins, which has finger-like DNA binding do-mains and plays an important role in many biological processes [4]. So far, zinc finger motifs have been classified into eight different classes based on their mainchain conformations and secondary structures around zinc-binding sites, including Cys2his2(C2H2)-like, ZN2/Cys6, Treble clef, zinc band, Gag joint, Taz2 domain-like, zinc-binding ring and metallothionein [5]. Due to the diversity of zinc finger motifs and these domains, ZFPs play different gene regulatory roles in different cellular environments and stimuli. ZNF306 promotes the development of colorectal cancer by transcriptionally activating integrin β4 and vascular endothelial growth factor [6]. ZNF384 promotes the proliferation of Hepatocellular carcinoma by directly up-regulating the expression of cyclin D1 [7]. Upregulation of ZNF554 is a potential tumor suppressor and its decreased expression may lead to the loss of oncogene suppression, activation of tumor pathways, and shorter survival of patients with malignant glioma [8]. The overexpression of ZNF655 promoted the progression of glioma by binding to the promoter of AURKA [9]. ZNF410 represents a special class of gene regulators, a conserved TF, which has a unique regulatory role on chromatin subcomplexes [10].Taken together, these studies indicate that ZNF genes may function as oncogenes involved in the occurrence and progression of cancer.

With the development of large-scale genome sequencing technologies, the integration of prognostic-related genetic markers has improved the level of early diagnosis of cancer compared with traditional clinical parameters. In the current study, we screened prognostic-associated ZNF family genes from TCGA dataset and validated the prognostic value of the six-gene signature using GEO dataset. We also constructed a nomogram based on the riskScore and clinical characteristics to predict individual overall survival (OS). In conclusion, our work may contribute to the early diagnosis of ESCA patients.

Materials and methods

Collection of data

RNA sequencing data of ESCA patients were extracted from the Cancer Genome Atlas (TCGA, https://cancergenome.nih.gov) and Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/). Relevant clinical information, including age, gender grade, survival status, and TNM stage, were also acquired. In this study, we collected 163 ESCA samples and 11 adjacent normal samples from TCGA and 119 ESCA samples from GSE53624 dataset. It has become evident that normalization of RNA-Seq data is necessary for reliable inferences and replication of results [11] [12]. GEO data is in TPM format, whereas TCGA data is in FPKM format. In this study, we first transform TCGA data into TPM format, then normalize TCGA and GEO data using the SVA package to minimize the batch effect.The expression levels of ZNF family genes were extracted for survival analysis, and the prognosis-related ZNF genes were identified. The prognostic model was constructed by using these genes, then verify the accuracy of the model.

Establishment and testing of the risk score model

We randomly separated the TCGA-ESCA patient samples into training and testing groups. The training cohort was used to establish construct the prognosis model and the testing cohort was used to verify the model [13]. By applying differential expression analysis, univariate Cox regression analysis, least absolute shrinkage and selection operator (LASSO) regression analysis and multivariate Cox regression analysis, we identified six prognostic related ZNF family genes which can be viewed as a signature to predict the disease outcomes of patients with ESCA. Subsequently, we divided the training group into high- and low-risk groups using the median risk score as the cut-off point. Kaplan–Meier survival curve analysis was performed to show the difference in OS between the two groups. We also plotted the time-dependent receiver operating characteristic (ROC) curve and risk nomogram to evaluate the prediction accuracy of the model.

Furthermore, univariate and multivariate Cox analyses were performed to confirm that the six prognostic-related genes signature was an independent prognostic factor for ESCA compared to other clinical characteristics—such as age, sex, clinical stage, TNM, and risk score. We analyzed the correlation between the risk scores and clinical parameters to investigate whether there was any difference in risk scores among different clinical parameter stratification.

Gene-set enrichment analysis

Gene ontology (GO) term enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to investigate potential signaling pathways and functions related to the six ZNF family genes included in the model [14].

Immune infiltration analysis and prediction of the sensitivity toward chemotherapeutic agents

We used the R package “GSEABase” to investigate the differential expression and function of 23 infiltrating immune cells between the high- and low-risk groups. We also used the R package “pRRophetic” to analyze drug sensitivity in the high- and low-risk groups.

Quantitative real-time PCR

Total RNA was extracted using Trizol (Cat#9109; TaKaRa, Japan) and reversed into cDNA by PrimeScript RT Master Mix kit (RR047A; Takara, Japan). All mRNA levels were assessed using the SYBR Green PCR Mix (RR420A; Takara, Japan) and the CFX Connect (BIO-RAD, USA). All experiments were performed in triplicate and analyzed with the 2−ΔΔCT method. Supplementary Table 1 lists the primer sequences.

Statistical analysis

All data were expressed as the mean and standard deviation of at least three independent experiments and analyzed by GraphPad Prism 7.0 software (La Jolla, CA, USA). R software (version x64 4.0.5) was used for statistical analysis, including Cox regression analysis, ROC curve analysis, gene enrichment analysis, and immune infiltration analysis.

Results

Features of patients with ESCA enrolled in this study

Sequencing data and corresponding clinical data of ESCA, including 163 tumors and 11 paracancerous tissues, were obtained from the TCGA database. Clinical information of 163 ESCA patients including survival status, grade of tumor, American Joint Committee on Cancer (AJCC) stage, classification of tumor, lymph node, and metastasis was summarized in Supplementary Table 2.

Construction and evaluation of the prognostic risk signature

This study was designed to investigate the prognostic significance of ZNF family in ESCA. The process of constructing and validating the six prognosis-related ZNF family genes signature is shown in Fig. 1. First, we performed univariate Cox regression and identified 12 remarkably prognosis-associated gens (ZNF641, ZNF484, ZNF629, ZNF91, ZNF45, ZNF502, ZNF25, ZNF586, ZNF510, ZNF865, ZNF106, and ZNF225) (Fig. 2A). Then, we used LASSO regression to build the prognostic signature and determine the coefficients (Fig. 2B -C). Finally, six ZNF family genes were enrolled in the signature, and each coefficient represented the weight of the expression of the corresponding coefficients (P < 0.05, Supplementary Table 3). The unpaired t-test was used to analyze the expression differences of six prognosis-related ZNF family genes in 163 ESCA patients and 11 normal esophageal samples. These findings suggested that compared with normal esophageal tissues, the expression of ZNF91 and ZNF586 were upregulated, while ZNF502, ZNF865, ZNF106, ZNF225 were downregulated in ESCA patients (Fig. 3A-B).

Fig. 1
figure 1

Flowchart for generating and validating the six prognosis-related ZNF family genes signature

Fig. 2
figure 2

Selection of prognostic ZNF family genes with prognostic value. (A) Risk ratio forest plot shows that twelve prognosis-related ZNF family genes, were significantly related to OS of ESCA patients. (B) Adjusted parameters of LASSO regression model. (C) Figure for LASSO coefficient spectrum of prognostic ZNF family genes

Fig. 3
figure 3

ZNF family genes signature predicts overall survival in patients with ESCA. (A-B) A heatmap and box plot showed the differential expression of six prognosis-related ZNF family genes between high-and low-risk subgroups. The gene expression was scaled by log2 (original expression of gene + 1). (C-D) The distribution of risk scores for each patient. With the median risk score as the cutoff, ESCA patients were divided into high- and low-risk subgroups. (E) Relationship between survival time (years) and survival status for each patient. (F) Kaplan-Meier curve of patients in the high- and low-risk subgroups to validate the predictive value of ZNF genes signature. The difference between the high- and low-risk subgroups was measured by the log-rank test, with a P-value < 0.001. OS, Overall Survival

Next, we used multivariate Cox regression to calculate their respective coefficients (βi) to establish a risk score model. We set the median risk.

score as the cutoff value and divided 173 patients into high-risk and low-risk groups (Fig. 3C). Then, we performed PCA to assess the distinct distribution between the high- and low-risk groups. Patients tended to separate into two clusters, which clearly indicated that the status of ESCA patients in the two risk score groups was different (Fig. 3D). As shown in the scatter plot, The higher the risk the more death patients (Fig. 3E). Additionally, the negative correlation between risk score and prognosis was affirmed by Kaplan-Meier survival curve (P < 0.001, Fig. 3F).

Testing of the risk score model

After calculating the risk scores of all patients in TCGA, we divided the training set and testing set samples into high- and low-risk groups according to the median value of risk score, as shown in Fig. 4A and B. We found that in both the training and testing sets, the proportion of patients with ESCA who died in the high-risk group was higher than that in the low-risk group (Fig. 4C-D). Moreover, six prognosis-related ZNF family genes in the risk model showed the same expression pattern in the training and testing sets (Fig. 4E-F). The Kaplan-Meier survival curve showed that the clinical outcomes of patients in the low-risk group were better than those in the high-risk group (P = 0.001; Fig. 4G), both in the testing set (P = 0.02; Fig. 4H).

Fig. 4
figure 4

Development and validation of the risk model for patients with ESCA. (A, B) Distribution of the ESCA patients with different risk scores in the training set and testing set. According to the median of the patient’s risk score, the ESCA patients were divided into high- (red) and low-risk(blue) groups. (C, D) The distribution of survival status of ESCA patients. The blue blots represent the patients who are alive, and red represents the patients who are dead. (E, F) Heat map depicting the expression patterns in the six prognosis-related ZNF family genes between high- and low-risk groups. (G, H) Overall survival (OS) of high-risk and low-risk patients in the testing group and training group

Correlation between the model and the clinical parameters

We then analyzed the correlation between the risk scores and other clinical parameters. According to the Kaplan-Meier analysis, the model based on six prognosis-related ZNF family genes had significantly distinct risk stratification ability in ESCA. The results presented that the high-risk patients in ESCA exhibited obviously worse prognosis (Fig. 5A-G, I). Whereas, in the patients with N2-3 (Fig. 5H), M1 (Fig. 5J) subgroups, this conclusion did not hold.

Fig. 5
figure 5

Kaplan-Meier curves showing the differences in prognosis between the high- and low- risk groups in different clinical subgroups, including male (A), female (B), stage I-II (C), stage III-IV (D), T0-2 (E), T3-4 (F), N0-1 (G), N2-3 (H), M0(I), M1(J)

Independent prognostic analysis and construction of a nomogram

Both univariate and multivariate Cox analyses showed gender, stage, T grade and the prognostic risk model could be used independently to predict the prognosis of ESCA (Fig. 6A-B). We then further compared these variables and found that the risk score was more accurate than the pathological stage and age in predicting OS at one years. The AUCs at one years for the risk score, gender, stage and the T grade were 0.848, 0.500, 0.625 and 0.536, respectively (Fig. 6C). We also drew a time-dependent ROC curve for the patients in the two groups (Fig. 6D). The AUC values for the risk score at 1, 3 and 5 years were 0.848, 0.872, 0.952. We also constructed a nomogram to estimate the probability of survival at 1, 3 and 5 years. The predictive factors including gender, stage, T stage, and the ZNF family genes prognostic signature, were used to construct the nomogram for OS (Fig. 6E). The C-index value of the nomogram was 0.826. The calibration curves depicting the actual and nomogram-predicted survival at 1, 3 and 5 years were relatively in accord with the reference lines (Fig. 6F). These results suggest that the nomogram including our prognostic signature is precise and reliable.

Fig. 6
figure 6

Assessing risk factors and constructing nomogram of prognosis. Univariate analysis (A) and multivariate analysis (B) were performed for screening of risk factors. (C) The ROC curves of clinicopathological characteristics and risk score for 1-year OS. (D) The ROC curves for 1-, 3-, and 5-year OS (E) An established nomogram model incorporated with the six prognosis-related ZNF family genes and clinicopathological parameters for prediction of OS in the TCGA dataset. (F) Calibration curves showed the concordances between predicted and observed 1-, 3-, and 5-year survival rates of ESCA patients based on the nomogram after bias corrections

External validation of the prognostic gene signature

To confirm the prognostic model had similar predictive values in different populations, the GEO cohort was used for external confirmation. Supplementary Table 2 shows the demographics and clinicopathologic characteristics of ESCA patients in the GEO validation cohort. Similarly, we performed univariate (riskScore (P < 0.001)) and multivariate Cox regression analyses (riskScore (P < 0.001)) to evaluate the prognostic significance of the model combined with various clinicopathologic parameters (Fig. 7A-B). In addition, risk score was more accurate than the clinicopathologic parameters in predicting OS at one years. The AUCs at one years for the riskScore, age, gender, T, N, and stage were 0.815, 0.616, 0.511, 0.605, 0.505 and 0.557, respectively (Fig. 7C). Therefore, the prognostic riskScore model constructed by the GEO validation cohorts was an independent prognostic factor for ESCA. The riskScore model also showed a favorable predictive ability for the 1-, 3- and 5-year OS rates, with AUC values of 0.815, 0.798 and 0.783, respectively (Fig. 7D). Furthermore, total of 119 ESCA patients in the GEO set were classified into low- and high-risk group (Fig. 7E-F), and the OS of the ESCA patients in the high-risk group was significantly lower than that of the patients in the low-risk group (P < 0.001; Fig. 7G). Finally, we established a prognostic nomogram to predict the survival probability based on the GEO validation cohort (Fig. 8A). The calibration curves depicting the actual and nomogram-predicted survival at 1, 3 and 5 years were relatively in accord with the reference lines (Fig. 8B).

Fig. 7
figure 7

External validation of the prognostic gene signature. Univariate analysis (A) and multivariate analysis (B) were performed for screening of risk factors in GEO dataset. (C) The ROC curves of clinicopathological characteristics and risk score for 1-year OS. (D) The ROC curves for 1-, 3-, and 5-year OS. (E) Distribution of the ESCA patients with different risk scores in high- and low riskScore groups. (F) The distribution of survival status of ESCA patients. (G) Overall survival (OS) of high-risk and low-risk patients

Fig. 8
figure 8

Established nomogram model. (A) An established nomogram model incorporated with the six prognosis-related ZNF family genes and clinicopathological parameters for prediction of OS in the GEO dataset. (B) Calibration curves showed the concordances between predicted and observed 1-, 3-, and 5-year survival rates of ESCA patients based on the nomogram after bias corrections

Functional enrichment analysis and evaluation of immune cell infiltration

To evaluate the infiltration scores of immune cells and immune-related functions, we performed ssGSEA analysis to quantify the scores of immune cell infiltration and immunity-related functions. In the high-risk patients, activated dendritic cell, CD56dim natural killer cell, macrophage, neutrophil and type 17 T helper cell were significantly higher than in low-risk patients (P < 0.01; Fig. 9A). In the high-risk patients, the functions were at higher levels, including ABC-co-inhibition, CCR, cytolytic activity, parainflammation, and TIL (P < 0.05; Fig. 9B). Our investigation indicated that the high-risk group had elevated immune activity, which might contribute to the occurrence of cancer. We performed GO analysis of the mRNAs co-expressed with the six prognosis-related ZNF family genes. GO and KEGG enrichment analysis showed that the target genes were mainly enriched in cell cycle, DNA replication and Fanconi anemia pathway (Fig. 9C-D).

Fig. 9
figure 9

Functional enrichment analysis and immune analysis. (A) Comparison of the infiltration of 21 immune cells between the different risk-groups. (B) Comparison of 14 immune-related functions between the different risk-groups. (C) GO analysis of highly related mRNAs. (D) KEGG analysis of highly related mRNAs. BP: biological process. CC: cellular component. MF: molecular function. *P < 0.05, **P < 0.01, ***P < 0.001

Analysis of drug sensitivity

We analyzed the sensitivity difference of chemotherapy drugs for ESCA in the current stage of clinical trials between the high- and low-risk groups, with the drug sensitivity expressed by IC50. We showed that patients in the high-risk group were more sensitive to AP-24,534 (P = 0.019, Fig. 10A), BMS-509,744 (P = 0.011, Fig. 10B), CGP-082996 (P = 0.04, Fig. 10C), HG-6-64-1 (P = 0.05, Fig. 10E), MG-132 (P = 0.015, Fig. 10F), Midostaurin (P = 0.037 Fig. 10G,), Ruxolitinib (P = 0.015 Fig. 10H,) Sunitinib (P = 0.022, Fig. 10I), TAE684 (P = 0.041, Fig. 10J), and Thapsigargin (P = 0.0084, Fig. 10K), whereas patients in low-risk group were more sensitive to Erlotinib (P = 0.017, Fig. 10D), which indicated that the model could be used as a potential predictor of chemotherapy sensitivity.

Fig. 10
figure 10

Drug sensitivity analysis to drugs of high- and low-risk subgroups. Differential chemotherapeutic responses in high- and low-risk patients (A-K).

External validation in expression

RT-qPCR assay was conducted to quantify the mRNA level of the six ZNF family genes in 12 pairs of esophageal cancer tissue samples (Fig. 11). Results indicate that ZNF91 (Fig. 11A), and ZNF586 (Fig. 11C) were upregulated in esophageal cancer tissue than adjacent tissue, whereas ZNF502 (Fig. 11B), ZNF865 (Fig. 11D), ZNF106 (Fig. 11E), and ZNF225 (Fig. 11F) was downregulated.

Fig. 11
figure 11

Expression of six prognosis-related ZNF family genes in twelve pairs of esophageal cancer tissue samples. *P < 0.05; **P < 0.01; ns > 0.05

Discussion

Esophageal cancer is one of the most common cancers with high morbidity and mortality [15]. Current therapeutic strategies for ESCA include surgery, chemotherapy, radiotherapy, molecular targeted therapy and their combination [16, 17]. In addition, immunotherapy is also playing an increasingly significant role [2, 18]. However, the prognosis remains poor and the overall Five-year survival rate is very low [19]. Therefore, achieving early diagnosis and effective treatment remains challenging. Identification of novel biomarkers will help assess prognosis, screen out patients in need of Immune intervention and drug therapy. Zinc finger proteins (ZFPs) primarily function as transcription factors in tumorigenesis and tumor progression involved in various tumor, such as esophageal squamous cell carcinoma cells [20],lung cancer [21], hepatocellular carcinoma [7], kidney renal clear cell carcinoma [22], oral squamous cell carcinoma [23]. Transcription factors (TFS) are proteins that play important roles in complex biological processes, such as metabolism, autophagy, apoptosis, immune response, stem cell maintenance and differentiation [24].

The TCGA database of 163 cases of esophageal cancer has improved our ability to diagnose, treat, and prevent cancer [25]. Based on the mRNA expression matrix and clinical data from the TCGA-ESCA cohort, we identified six prognostic-associated ZNF genes that may be clinically valuable biomarkers. Patients with ESCA were divided into two subgroups with different survival outcomes based on a prognostic model of six ZNF family genes. We also established a risk score model to predict the prognoses of patients with ESCA based on these prognostic genes. Importantly, the ability of the prognostic model to distinguish high-and low-risk patients, and to estimate OS, was similarly validated in the GSE53624 dataset. Moreover, we combined risk score with other clinical variables to conduct a nomogram to establish a quantitative prognostic evaluation method for patients with ESCA.

Our ZNF family gene-based signature included six genes, i.e., ZNF91, ZNF586, ZNF502, ZNF865, ZNF106 and ZNF225. ZNF91 is likely to play an important role in cell proliferation and/or anti-apoptosis, and may serve as a molecular marker for AML [26]. Upregulation of ZNF91 could promote irradiation resistance by regulating the stem cell-like properties of NSCLC cells. Abnormal expression of ZNF91 is related to the occurrence and development of bladder cancer [26, 27], colorectal cancer (CRC) [28]and ovarian cancer [29, 30]. Our results also showed that ZNF91 is upregulated and plays an oncogenic role in ESCA. Genome-wide differential gene/microRNA signatures show that ZNF502 might be a prognostic biomarker in cytogenetically normal acute myeloid leukemia [31]. In our study, the expression of ZNF502 is low in esophageal carcinoma, which is associated with poor prognosis. ZNF106 is a RNA-binding protein that binds to the core splicing factor RNA-binding motif protein 39 and localizes to nuclear spots near the spliceosome [32]. We found that mRNA level of 106 was significantly reduced in ESCA tissues, this was similar to the reported results that ZNF106 expression was downregulated and associated with a good predictive value in Bladder Cancer [33, 34]. Little research has been done on the role of ZNF225 in ESCA, and a few evidences suggest ZNF225 inhibits autophagy and promotes apoptosis of hepatocellular carcinoma cells. Our findings on ZNF225 are supported by evidence that this ZNF protein serves as prognostic genes. In contrast, the roles of ZNF586, and ZNF865 in ESCA onset and development had not, to our knowledge, been as yet explored. Based on current knowledge, our findings suggest that the six prognosis-related ZNF family genes may exert important roles in the tumorigenesis and progression of ESCA.

By univariate Cox analysis, lasso regression and multivariate Cox analysis, we screened six prognosis-related ZNF family genes to construct the prognostic model. Survival and ROC curve analyses showed that these six genes had good diagnostic ability and could be used to screen out ESCA patients who had poor prognoses.when compared to previously reported models [35, 36], our model(At 1, 3, and 5 years, the risk score’s AUC values were 0.848, 0.872, and 0.952, respectively) has greater predictive power. However, the specific molecular mechanisms of these six prognosis-related ZNF family genes in ESCA remain unclear, and the underlying molecular mechanisms should be explored. Subsequently, we assessed the relationship between risk score model and clinical variables and found that the risk score model had significantly distinct risk stratification ability in ESCA. Nomogram has long been used in oncology to calculate the prognosis of patients with esophageal cancer based on the relevant clinical parameters [37,38,39]. We then established a nomogram to more intuitively predict 1-year, 3-year, and 5-year survival estimates in patients with ESCA and found that the risk score was more accurate than the pathological stage and age in predicting OS from TCGA and GEO dataset.

Next, we analyzed the differences in immune cell infiltration and response rates to chemotherapy sensitivity among different groups of patients with ESCA based on the model. ESCA was enriched in immune-suppressive cell populations, including Tregs, exhausted CD8 T, CD4 T and NK cells, M2 macrophages, and tDCs [40]. Accumulating evidence regards the tumor immune microenvironment can potentially influence the patient’s response to immune checkpoint inhibitors, tumor immunity, such as PD-L1 expression on tumors, tumor-infiltrating lymphocytes and tumor-associated macrophages [41]. Thus, our model can be used as an indicator to predict immune cell infiltration and immune response in patients with ESCA. At present, surgical resection, radiotherapy and chemotherapy are the main clinical treatment methods for ESCA, However, due to the limited efficacy and serious adverse effects of conventional treatment, the result is still unsatisfactory. As a new treatment method, target therapy has a good application prospect [41, 42]. In our study, patients in the high-risk group were more sensitive to AP-24,534, BMS-509,744, CGP-082996, HG-6-64-1, MG-132, Midostaurin, Ruxolitinib, Sunitinib, TAE684and Thapsigargin. From what has been discussed above, our results revealed differences in immune cell infiltration and immune response between the groups. ZNF-gene signature for ESCA was able to predict chemotherapy sensitivity and may thus help guide treatment selection.

We further conducted GO and KEGG analyses to evaluate biological functions. Enrichment analysis of biological functions and pathways of the ZNF family gene indicated that these prognosis gene were significantly enriched in Cell cycle, DNA replication and Fanconi anemia pathway.

Summarily, this study found that ZNF genes were differentially expressed in ESCA tissues and the reason may be different from the mechanism in the process of tumor formation. As the largest transcription regulator family in mammals, zinc finger (ZNF) protein expression regulation mechanism is very complex, including Genetic variation [43], Epigenetic modifications [44] and Posttranslational regulation [45]. We used the Kaplan-Meier analysis to study the prognostic significance of the six prognosis-related ZNF family genes and found that ZNF502, ZNF865, ZNF106 and ZNF225 gene expressions were related to good prognoses in patients with ESCA, while high ZNF586 and ZNF91 gene expressions were related to poor prognoses. We further confirmed the expressions of these genes at the tissue level. The results suggested that the signatures of these six genes may assess treatment outcomes and predict patient survival.

However, the current study has multiple limitations. Firstly, there are few normal tissues in TCGA database, which need to be verified by expanded samples. Second, the functional relationship between the ZNF gene signature members and non-tumor cells in the tumor microenvironment, especially infiltrating immune cells, could not be elucidated and requires future in vitro and in vivo studies. The effect on proliferation, invasion and migration of ZNF family genes in ESCA requires further be verified in vitro and in vivo.

Conclusion

In summary, we first constructed a prognostic model of ESCA based on features of ZNF family genes that divides ESCA patients into two subgroups with different survival outcomes and constructed a nomogram to help clinical decision-makers provide optimal treatment. The prognostic signature is related to different immune cell and predicts sensitivity to chemotherapeutic agents which might be novel targets for developing immunotherapies in low-risk and high-risk ESCA patients. These prognosis-related ZNF family genes may play vital roles in ESCA occurrence, progression, invasion and metastasis. Additionally, these findings have led to the development of new clinical therapeutic targets or prognostic marker.

Data Availability

(ADM)

The Supplementary Material for this article can be found online at this study, further inquiries can be directed to the corresponding authors. The datasets analyzed during the current study are available online from the GSE53624 dataset (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE53624) and The Cancer Genome Atlas Program ().

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. 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. https://doi.org/10.3322/caac.21660.

    Article  PubMed  Google Scholar 

  2. Kelly RJ. Emerging multimodality approaches to treat localized esophageal Cancer. J Natl Compr Canc Netw. 2019;17(8):1009–14. https://doi.org/10.6004/jnccn.2019.7337.

    Article  CAS  PubMed  Google Scholar 

  3. Pennathur A, Gibson MK, Jobe BA, Luketich JD. Oesophageal carcinoma. Lancet. 2013;381(9864):400–12. https://doi.org/10.1016/S0140-6736(12)60643-6.

    Article  PubMed  Google Scholar 

  4. Li X, Han M, Zhang H, Liu F, Pan Y, Zhu J, Liao Z, Chen X, Zhang B. Structures and biological functions of zinc finger proteins and their roles in hepatocellular carcinoma. Biomark Res. 2022;10(1):2. https://doi.org/10.1186/s40364-021-00345-1.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Klug A. The discovery of zinc fingers and their development for practical applications in gene regulation and genome manipulation. Q Rev Biophys. 2010;43(1):1–21. https://doi.org/10.1017/S0033583510000089.

    Article  CAS  PubMed  Google Scholar 

  6. Yang L, Hamilton SR, Sood A, Kuwai T, Ellis L, Sanguino A, Lopez-Berestein G, Boyd DD. The previously undescribed ZKSCAN3 (ZNF306) is a novel “driver” of colorectal cancer progression. Cancer Res. 2008;68(11):4321–30. https://doi.org/10.1158/0008-5472.CAN-08-0407.

    Article  CAS  PubMed  Google Scholar 

  7. He L, Fan X, Li Y, Chen M, Cui B, Chen G, Dai Y, Zhou D, Hu X, Lin H. Overexpression of zinc finger protein 384 (ZNF 384), a poor prognostic predictor, promotes cell growth by upregulating the expression of cyclin D1 in hepatocellular carcinoma. Cell Death Dis. 2019;10(6):444. https://doi.org/10.1038/s41419-019-1681-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Balogh A, Reiniger L, Hetey S, Kiraly P, Toth E, Karaszi K, Juhasz K, Gelencser Z, Zvara A, Szilagyi A, et al. Decreased expression of ZNF554 in Gliomas is Associated with the activation of Tumor Pathways and shorter patient survival. Int J Mol Sci. 2020;21(16). https://doi.org/10.3390/ijms21165762.

  9. Chen X, Liu C, Zhang Z, Wang M, Guo S, Li T, Sun H, Zhang P. ZNF655 promotes the progression of Glioma through Transcriptional Regulation of AURKA. Front Oncol. 2022;12:770013. https://doi.org/10.3389/fonc.2022.770013.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Vinjamur DS, Yao Q, Cole MA, McGuckin C, Ren C, Zeng J, Hossain M, Luk K, Wolfe SA, Pinello L, et al. ZNF410 represses fetal globin by singular control of CHD4. Nat Genet. 2021;53(5):719–28. https://doi.org/10.1038/s41588-021-00843-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Zhao Y, Li MC, Konate MM, Chen L, Das B, Karlovich C, Williams PM, Evrard YA, Doroshow JH, McShane LM. TPM, FPKM, or normalized Counts? A comparative study of quantification measures for the analysis of RNA-seq data from the NCI patient-derived models repository. J Transl Med. 2021;19(1):269. https://doi.org/10.1186/s12967-021-02936-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Zhao S, Ye Z, Stanton R. Misuse of RPKM or TPM normalization when comparing across samples and sequencing protocols. RNA. 2020;26(8):903–9. https://doi.org/10.1261/rna.074922.120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Lu Y, Wang W, Liu Z, Ma J, Zhou X, Fu W. Long non-coding RNA profile study identifies a metabolism-related signature for colorectal cancer. Mol Med. 2021;27(1):83. https://doi.org/10.1186/s10020-021-00343-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–92. https://doi.org/10.1093/nar/gkac963.

    Article  PubMed  Google Scholar 

  15. Fan J, Liu Z, Mao X, Tong X, Zhang T, Suo C, Chen X. Global trends in the incidence and mortality of esophageal cancer from 1990 to 2017. Cancer Med. 2020;9(18):6875–87. https://doi.org/10.1002/cam4.3338.

    Article  PubMed  Google Scholar 

  16. He S, Xu J, Liu X, Zhen Y. Advances and challenges in the treatment of esophageal cancer. Acta Pharm Sin B. 2021;11(11):3379–92. https://doi.org/10.1016/j.apsb.2021.03.008.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Yang YM, Hong P, Xu WW, He QY, Li B. Advances in targeted therapy for esophageal cancer. Signal Transduct Target Ther. 2020;5(1):229. https://doi.org/10.1038/s41392-020-00323-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. van den Ende T, van den Boorn HG, Hoonhout NM, van Etten-Jamaludin FS, Meijer SL, Derks S, de Gruijl TD, Bijlsma MF, van Oijen MGH, van Laarhoven HWM. Priming the tumor immune microenvironment with chemo(radio)therapy: a systematic review across tumor types. Biochim Biophys Acta Rev Cancer. 2020;1874(1):188386. https://doi.org/10.1016/j.bbcan.2020.188386.

    Article  CAS  PubMed  Google Scholar 

  19. Thrift AP. Global burden and epidemiology of Barrett oesophagus and oesophageal cancer. Nat Rev Gastroenterol Hepatol. 2021;18(6):432–43. https://doi.org/10.1038/s41575-021-00419-3.

    Article  PubMed  Google Scholar 

  20. Gong W, Xu J, Wang G, Li D, Zhan Q. ZNF292 suppresses proliferation of ESCC cells through ZNF292/SKP2/P27 signaling axis. Chin J Cancer Res. 2021;33(6):637–48. https://doi.org/10.21147/j.issn.1000-9604.2021.06.01.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Liao SY, Kuo IY, Chen YT, Liao PC, Liu YF, Wu HY, Lai WW, Wang YC. AKT-mediated phosphorylation enhances protein stability and transcription activity of ZNF322A to promote lung cancer progression. Oncogene. 2019;38(41):6723–36. https://doi.org/10.1038/s41388-019-0928-x.

    Article  CAS  PubMed  Google Scholar 

  22. Czerwinska P, Mackiewicz AA. Low levels of TRIM28-Interacting KRAB-ZNF genes associate with Cancer Stemness and Predict Poor prognosis of kidney renal clear cell carcinoma patients. Cancers (Basel). 2021;13(19). https://doi.org/10.3390/cancers13194835.

  23. Starzynska A, Sobocki BK, Sejda A, Sakowicz-Burkiewicz M, Szot O, Jereczek-Fossa BA. ZNF-281 as the potential diagnostic marker of oral squamous cell carcinoma. Cancers (Basel). 2021;13(11). https://doi.org/10.3390/cancers13112661.

  24. Jen J, Wang YC. Zinc finger proteins in cancer progression. J Biomed Sci. 2016;23(1):53. https://doi.org/10.1186/s12929-016-0269-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Liu J, Lichtenberg T, Hoadley KA, Poisson LM, Lazar AJ, Cherniack AD, Kovatich AJ, Benz CC, Levine DA, Lee AV, et al. An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics. Cell. 2018;173(2):400–416e411. https://doi.org/10.1016/j.cell.2018.02.052.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Unoki M, Okutsu J, Nakamura Y. Identification of a novel human gene, ZFP91, involved in acute myelogenous leukemia. Int J Oncol. 2003;22(6):1217–23.

    CAS  PubMed  Google Scholar 

  27. Han Y, Zheng Q, Tian Y, Ji Z, Ye H. Identification of a nine-gene panel as a prognostic indicator for recurrence with muscle-invasive bladder cancer. J Surg Oncol. 2019;119(8):1145–54. https://doi.org/10.1002/jso.25446.

    Article  CAS  PubMed  Google Scholar 

  28. Zhao Q, Wang F, Chen YX, Chen S, Yao YC, Zeng ZL, Jiang TJ, Wang YN, Wu CY, Jing Y, et al. Comprehensive profiling of 1015 patients’ exomes reveals genomic-clinical associations in colorectal cancer. Nat Commun. 2022;13(1):2342. https://doi.org/10.1038/s41467-022-30062-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Micci F, Skotheim RI, Haugom L, Weimer J, Eibak AM, Abeler VM, Trope CG, Arnold N, Lothe RA, Heim S. Array-CGH analysis of microdissected chromosome 19 markers in ovarian carcinoma identifies candidate target genes. Genes Chromosomes Cancer. 2010;49(11):1046–53. https://doi.org/10.1002/gcc.20813.

    Article  CAS  PubMed  Google Scholar 

  30. Dufresne J, Bowden P, Thavarajah T, Florentinus-Mefailoski A, Chen ZZ, Tucholska M, Norzin T, Ho MT, Phan M, Mohamed N, et al. The plasma peptides of ovarian cancer. Clin Proteom. 2018;15:41. https://doi.org/10.1186/s12014-018-9215-z.

    Article  CAS  Google Scholar 

  31. Jinlong S, Lin F, Yonghui L, Li Y, Weidong W. Identification of let-7a-2-3p or/and mir-188-5p as prognostic biomarkers in cytogenetically normal acute myeloid leukemia. PLoS ONE. 2015;10(2):e0118099. https://doi.org/10.1371/journal.pone.0118099.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Anderson DM, Cannavino J, Li H, Anderson KM, Nelson BR, McAnally J, Bezprozvannaya S, Liu Y, Lin W, Liu N, et al. Severe muscle wasting and denervation in mice lacking the RNA-binding protein ZFP106. Proc Natl Acad Sci U S A. 2016;113(31):E4494–4503. https://doi.org/10.1073/pnas.1608423113.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Guo C, Shao T, Jiang X, Wei D, Wang Z, Li M, Bao G. Comprehensive analysis of the functions and prognostic significance of RNA-binding proteins in bladder urothelial carcinoma. Am J Transl Res. 2020;12(11):7160–73.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Wu Y, Liu Z, Wei X, Feng H, Hu B, Liu B, Luan Y, Ruan Y, Liu X, Liu Z, et al. Identification of the functions and prognostic values of RNA binding proteins in bladder Cancer. Front Genet. 2021;12:574196. https://doi.org/10.3389/fgene.2021.574196.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Zhang Z, Chen C, Fang Y, Li S, Wang X, Sun L, Zhou G, Ye J. Development of a prognostic signature for esophageal cancer based on nine immune related genes. BMC Cancer. 2021;21(1):113. https://doi.org/10.1186/s12885-021-07813-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Du H, Xie S, Guo W, Che J, Zhu L, Hang J, Li H. Development and validation of an autophagy-related prognostic signature in esophageal cancer. Ann Transl Med. 2021;9(4):317. https://doi.org/10.21037/atm-20-4541.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Semenkovich TR, Yan Y, Subramanian M, Meyers BF, Kozower BD, Nava R, Patterson GA, Kreisel D, Puri V. A clinical Nomogram for Predicting Node-positive disease in Esophageal Cancer. Ann Surg. 2021;273(6):e214–21. https://doi.org/10.1097/SLA.0000000000003450.

    Article  PubMed  Google Scholar 

  38. Sudo K, Wang X, Xiao L, Wadhwa R, Shiozaki H, Elimova E, Rice DC, Lee JH, Weston B, Bhutani MS, et al. A Nomogram to predict distant metastases after multimodality therapy for patients with localized esophageal Cancer. J Natl Compr Canc Netw. 2016;14(2):173–9. https://doi.org/10.6004/jnccn.2016.0020.

    Article  PubMed  Google Scholar 

  39. Gaur P, Sepesi B, Hofstetter WL, Correa AM, Bhutani MS, Vaporciyan AA, Watson TJ, Swisher SG et al. Members of the MDAECG, the University of Rochester School of M. A clinical nomogram predicting pathologic lymph node involvement in esophageal cancer patients. Ann Surg 2010; 252(4):611–617. doi:https://doi.org/10.1097/SLA.0b013e3181f56419

  40. Zheng Y, Chen Z, Han Y, Han L, Zou X, Zhou B, Hu R, Hao J, Bai S, Xiao H, et al. Immune suppressive landscape in the human esophageal squamous cell carcinoma microenvironment. Nat Commun. 2020;11(1):6268. https://doi.org/10.1038/s41467-020-20019-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Baba Y, Nomoto D, Okadome K, Ishimoto T, Iwatsuki M, Miyamoto Y, Yoshida N, Baba H. Tumor immune microenvironment and immune checkpoint inhibitors in esophageal squamous cell carcinoma. Cancer Sci. 2020;111(9):3132–41. https://doi.org/10.1111/cas.14541.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Wu X, Hu W, Lu L, Zhao Y, Zhou Y, Xiao Z, Zhang L, Zhang H, Li X, Li W, et al. Repurposing vitamin D for treatment of human malignancies via targeting tumor microenvironment. Acta Pharm Sin B. 2019;9(2):203–19. https://doi.org/10.1016/j.apsb.2018.09.002.

    Article  PubMed  Google Scholar 

  43. Hawe JS, Wilson R, Schmid KT, Zhou L, Lakshmanan LN, Lehne BC, Kuhnel B, Scott WR, Wielscher M, Yew YW, et al. Genetic variation influencing DNA methylation provides insights into molecular mechanisms regulating genomic function. Nat Genet. 2022;54(1):18–29. https://doi.org/10.1038/s41588-021-00969-x.

    Article  CAS  PubMed  Google Scholar 

  44. Longo M, Raciti GA, Zatterale F, Parrillo L, Desiderio A, Spinelli R, Hammarstedt A, Hedjazifar S, Hoffmann JM, Nigro C, et al. Epigenetic modifications of the Zfp/ZNF423 gene control murine adipogenic commitment and are dysregulated in human hypertrophic obesity. Diabetologia. 2018;61(2):369–80. https://doi.org/10.1007/s00125-017-4471-4.

    Article  CAS  PubMed  Google Scholar 

  45. Liu P, Xiao J, Wang Y, Song X, Huang L, Ren Z, Kitazato K, Wang Y. Posttranslational modification and beyond: interplay between histone deacetylase 6 and heat-shock protein 90. Mol Med. 2021;27(1):110. https://doi.org/10.1186/s10020-021-00375-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank AiMi Academic Services (www.aimieditor.com) for English language editing and review services.

Funding

This research was supported by the National Natural Science Foundation of China (81770638), the Science and Technology Fund of the Guizhou Provincial Health Commission (gzwkj2022-049), and Guizhou Provincial Science and Technology Projects (QKH JC-ZK [2022]-264).

Author information

Authors and Affiliations

Authors

Contributions

KQH designed the project and wrote the manuscript. QY processed the datasets and revised the manuscript. HSY and NW analyzed the data. WW and BPY oversaw project design. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Wei Wang or Baoping Yu.

Ethics declarations

Ethics approval and consent to participate

This study involving human participants was reviewed and approved by the Institutional Review Board of the Renmin Hospital of Wuhan University (No. WDRY2021-K062). All the experimental protocols were performed in accordance with the Declaration of Helsinki. Written informed consent was obtained from all the patients who participated in this study.

Consent for publication

Not applicable.

Competing interests

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.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Table 1

Primer sequence

Supplementary Table 3

Coefficients of six prognosis-related ZNF family genes in LASSO regression

Supplementary Table 2

The clinical data of esophageal cancer patients in TCGA and GEO

Supplementary Material 4

Supplementary Material 5

Supplementary Material 6

Supplementary Material 7

Supplementary Material 8

Supplementary Material 9

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hong, K., Yang, Q., Yin, H. et al. Comprehensive analysis of ZNF family genes in prognosis, immunity, and treatment of esophageal cancer. BMC Cancer 23, 301 (2023). https://doi.org/10.1186/s12885-023-10779-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-10779-5

Keywords