Abstract
Background/Aim: Clear-cell renal cell carcinoma (ccRCC) is the most common and aggressive form of all urological cancers, with poor prognosis and high mortality. Despite growing evidence of involvement in carcinogenesis, the role of KRAB-ZFP in ccRCC has not been fully explored. KRAB Zinc finger proteins (KRAB-ZFPs) are the largest family of mammalian transcription regulators. They are differentially expressed in various tissues during cellular development and phenotypic differentiation. Materials and Methods: In this study, the levels of transcripts of ccRCC from The Cancer Genome Atlas (TCGA) dataset were used to identify prognostic biomarkers in this disease. Results: Using bioinformatics techniques, we demonstrate that approximately 60% of KRAB zinc finger proteins located on chromosome 19p13.2 are differentially expressed, with all but two being down-regulated in ccRCC. Moreover, ZNF844, a paralog of ZNF433, was the most down-regulated across all histological grades and pathological stages (p<0.001). In addition, the decrease in ZNF844 expression was associated with poor patient survival (HR=0.41; 95% CI=0.3-0.56; p<0.0001). Gene Set Enrichment Analysis of genes inversely co-expressed with ZNF844 revealed that enriched pathways were consistently related to immune and translation processes (p<0.05, FDR <0.05). Lastly, ZNF844 expression showed moderate, inverse correlation to Helper T-cell (CD4 or Th1) subtype 1 (R=–0.558, p=5.15×10–39) infiltration and with the exhausted T-cell phenotype (R=–0.37; p=4.1×10–21). Conclusion: Down-regulation of KRAB-ZFPs at 19p13.2 may represent a signature for ccRCC. Moreover, ZNF844 is a prognostic marker for ccRCC and may serve as a putative immune-related tumor suppressor gene.
- ZNF844
- ZNF433
- clear cell renal carcinoma
- Kruppel-associated box
- KRAB-ZFPs
- tumor microenvironment
- ccRCC
- Th1 cells
- biomarker
- Chr19p13.2
Clear cell renal carcinoma (ccRCC) is the most common neoplastic disease affecting the kidneys. It accounts for at least 70% of all kidney cancers (1, 2). In advanced stages, ccRCC becomes resistant to conventional therapies (1-3). Moreover, later stages are associated with poor prognosis with a 5-year overall survival rate of less than 10% (1-3). Despite growing evidence of involvement in carcinogenesis, the role of KRAB-ZFP in ccRCC has not been fully explored. Kruppel-associated box domain-containing zinc-finger proteins (KRAB-ZFP) represent the most prominent family of transcriptional repressors (4, 5). The human genome contains over 300 KRAB-ZFPs genes that produce more than 800 transcripts that encode proteins, pseudogenes, and splice variants (6-8). KRAB-ZFPs are differentially expressed in various tissues during cellular development. Several of these genes are involved in embryogenesis, genome imprinting, and biological processes that determine phenotypes in stem cells (9-11). Given their roles in cellular differentiation and their ability to alter gene expression via epigenetic processes and chromatin remodeling, it is conceivable that dysregulation of KRAB-ZFP transcription factors could facilitate tumorigenesis and tumor progression.
Interestingly, several large KRAB-ZFP clusters are located on chromosome 19 (chr19), and several studies have linked genetic alterations in chromosome 19p13 to cancers (12-18). Using data mining techniques, we recently identified ZNF433, a zinc finger that is differentially expressed in ccRCC (19) and that may be a biomarker for ccRCC. ZNF433 is part of a cluster of KRAB-ZFPs that are located on chromosome 19p13.2. Thus, we focused on this locus to ascertain if the expression of other members of this cluster was altered. We have identified several KRAB-ZFPs that are differentially expressed and are associated with changes in overall survival in patients. Since the gene ZNF844 is a paralog of ZNF433, it appears to be a strong candidate biomarker and prognostic factor for ccRCC.
Materials and Methods
UALCAN (gene expression) analysis. In this study, UALCAN (http://ualcan.path.uab.edu/analysis.html) (accessed on June 5, 2021 and on August 25, 2021) was used to ascertain whether chromosome 19p13.2 KRAB-ZFPs were differentially expressed in ccRCC. UALCAN is an online open-access web tool that provides graphical visualizations of pathological clinic features, gene expression data, promoter methylation data derived from TCGA KIRC RNA-Seq data, as well as protein expression data derived from the CPTAC dataset (20). UALCAN was used to examine the mRNA expression of ZNF844 in relationship to clinicopathological parameters (cancer stage and grade) of ccRCC. To illustrate gene expression data in graphs, gene symbols for various KRAB-ZFPs were inputted into the ‘scan by gene’ query box, and the kidney renal clear cell carcinoma (KIRC) was selected. Graphical representations of gene expression data as a function of clinicopathological and demographic data were obtained by selecting the search feature on the website.
UCSC Xena and GEO databases (gene expression validation). To validate gene expression results obtained from UALCAN, TCGA KIRC clinical information, and RNA sequencing dataset (n=538) were downloaded from the UCSC XENA portal (https://xena.ucsc.edu/public/) (21) on July 27, 2021. To obtain the data, the KIRC-TCGA study was selected. Under the querying system, genomic data of KRAB ZFP of interest was chosen as the first variable. For the second variable, phenotypic data for the sample type, age at initial pathologic diagnosis, sex, pathologic N, pathologic stage, neoplastic histological grade, and genomic information for ID TCGA-HiSeqV2 data were selected. Primary tumors were selected as a filter to remove duplications in the data, leaving 538 samples. To ascertain survival data, KM analysis feature was selected under the first variable results, and the raw data was downloaded. Data was utilized to generate Kaplan-Meier plot to ascertain the impact of gene expression on patient survival. Moreover, data was also used to perform ROC/AUC analysis as well as univariate and multivariate analyses to ascertain whether ZNF844 and other KRAB-ZFPs were predictive ccRCC biomarkers. In addition, we also downloaded patient data (GSE12606) from Gene Expression Omnibus GEO (http://www.ncbi.nlm.nih.gov/geo/) database (22). GEO2R was used to analyze ZNF844 gene expression between normal and metastatic ccRCC. The adjusted p-Value significance was set at 0.05.
GEPIA database analysis (gene expression and correlation analysis). GEPIA2 database is an enhanced version of GEPIA. This gene expression database includes gene and corresponding isoforms in normal and tumor samples of 84 cancer subtypes derived from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases (23-25). GEPIA2 was used to conduct correlation analysis between ZNF844 expression and T-cell signatures. GEPIA2 was accessed on December 12, 2021 (http://gepia2.cancer-pku.cn). To perform this function, the “Correlation Analysis” module was chosen, and ZNF844 for was entered into the Gene A query box, and the desired T-cell signature was chosen. These included naïve T-cell signature (genes: CCR7, LEF1, TCF7, SELL), Helper T-cell Type 1 signature (genes: CXCL13, HAVCR2, IFNG, CXCR3, CD4, BHLHE40), Helper T-cell Type 2 signature (genes: CCL26, IL5, IL13, GATA3, STAT6), Effector T-reg signature (genes: FOXP3, CTLA4, CCR8, TNFRSF9), Effector T-cell signature (genes: CX3CR1, FGFBP2, FCGR3A), and the Exhausted T-cell signature (genes: HAVCR2, TIGIT, LAG3, PDCD1, CXCL13, LAYN). To analyze the correlation between ZNF844 and various aforementioned gene signatures, Pearson correlation coefficient was utilized. p-Value significance was set at 0.05.
Human Protein Atlas Database (immunohistochemistry). Human Protein Atlas (HPA) was used to ascertain if the ZNF844 protein was differentially expressed in normal kidney and renal carcinomas (HPA, https://www.proteinatlas.org/) (accessed on June 5, 2021) (26). HPA provides qualitative analysis of the protein expression relative to normal tissues. In addition, HPA RNA-seq data were used to identify single cell types and gene clusters associated with expression of ZNF844 (accessed on February 2, 2022).
Timer 2.0. Database (tumor microenvironment). TIMER2.0 (http://timer.cistrome.org/) (accessed on June 26, 2021) was used to examine ZNF844 expression versus immune infiltration using the XCell algorithm. The TIMER 2.0 web tool provides a robust estimation of immune infiltration levels for The Cancer Genome Atlas (TCGA) or user-provided tumor profiles using six state-of-the-art algorithms (27, 28). In addition, TIMER2.0 provides modules for investigating the associations between immune infiltrates and genetic or clinical features, clinical outcomes and for exploring cancer-related associations in the TCGA cohorts (27, 28). All TCGA tumor data, including transcriptome profiles, somatic mutation calls, somatic copy number variations, and patient clinical outcomes, are collected from the GDAC firehose website (27, 28).
Linkedomics webtool (gene set enrichment analysis). The gene ontology (GO) enrichment analysis results of the genes co-expressed with ZNF844 were visualized using LinkedOmics (http://www.linkedomics.org) (accessed on July 10, 2021). The LinkedOmics database contains multi-omics data and clinical data for 32 cancer types comprising 11,158 patients from the TCGA project (28). The “LinkFinder” module was used to visualize the differentially expressed genes within the TCGA KIRC cohort. The “LinkInterpreter” module was used to perform GSEA and ORA enrichment analysis methods to identify pathways that are enriched among genes positively and inversely expressed with ZNF844 within the GO Biological Process KEGG and Reactome databases. Results were analyzed for significance using the Pearson’s correlation test with the p-value, and false discovery rate (FDR) set at 0.05.
Statistical analysis. For comparison gene expression in normal versus tumor samples, between male versus female between low grade and high-grade tumors, or between low and high-grade tumors, the Student t-test was used. Comparisons of ZNF844 gene expression across tumor stages and grades were conducted using one-way ANOVA with Tukey-Kramer post hoc analysis. Analysis of gene expression between the low histological grades or pathological stages versus higher grades and stages was conducted using the Student t-test. For the Kaplan-Meier (KM) survival outcome plots, Cox risk hazard ratio (HR), 95% confidence intervals (CI), and log-rank p-Values were generated by the KM plotter web tool. Receiver Operator Curve (ROC) analysis was performed to validate KM results using MedCalc statistical software and TCGA datasets downloaded via UCSC Xena webtool. An AUC greater than 0.6 is considered to mean that the gene has potential as a prognostic factor. The p-Value for all statistical sets was set at p<0.05 to indicate a statistical significance.
Results
Differential expression of KRAB-ZFPs at Chromosome 19p13.2. We analyzed 20 members of a cluster of KRAB-ZFPs located at the chromosome 19p13.2 locus. All (Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 17, Figure 18, Figure 19 and Figure 20) but one, ZNF562 (Figure 16), were differentially expressed. Four KRAB-ZFP transcripts [ZNF121 (Figure 4A), ZNF224 (Figure 6A), ZNF 266 (Figure 8A), and ZNF700 (Figure 18A)] were up-regulated in ccRCC and were associated with unfavorable survival outcomes. Out of the fifteen KRAB-ZFPs that were down-regulated (Figure 1A, Figure 2A, Figure 3A, Figure 5A, Figure 7A, Figure 9A, Figure 10A, Figure 11A, Figure 12A, Figure 13A, Figure 14A, Figure 15A, Figure 17A, Figure 19A and Figure 20A), ten proteins [ZNF69 (Figure 3B), ZNF426 (Figure 9B), ZNF440 (Figure 11B), ZNF441 (Figure 12B), ZNF442 (Figure 13B), ZNF443 (Figure 14B), ZNF561 (Figure 15B), ZNF564 (Figure 17B), ZNF709 (Figure 19B) and ZNF844 (Figure 20B)] were associated with significant favorable survival outcomes (p<0.05). Analysis of the chr19p13.2 KRAB-ZFPs revealed that ZNF844 had the most significant changes in both gene expression (p=1.624×10–12) and overall survival (p=9.4×10–8), suggesting that this transcription factor may be a possible biomarker and prognostic indicator (Figure 20). Indeed, survival outcomes for the low expression group (less than 7.512 transcripts/million) were significantly lower than the high expression group (greater than 7.512 transcripts/million). To ascertain whether ZNF844 expression could discriminate between high-expression and low-expression groups in relationship to survival outcome, we conducted ROC analysis. The AUC was 0.642, 95% CI=0.599 to 0.683, indicating that ZNF844 expression is a moderate predictor of survival outcomes. Therefore, we further analyzed ZNF844 relative to clinicopathological features of ccRCC.
ZNF844 gene and protein expression in patients with clear-cell renal cell carcinoma. Figure 21 depicts the expression of ZNF844 protein in clear cell renal carcinoma and normal tissues. We utilized Human Protein Atlas (https://www.proteinatlas.org/) to identify immunohistochemistry images representing human normal and renal carcinoma tissues stained with antibodies raised against the ZNF844 protein (Figure 21A). The renal carcinoma images showed qualitatively that there was lower ZNF844 protein expression, as documented by descriptions that accompanied each image. Tissues from renal adenocarcinoma had “low” or “not detected” staining with “weak” intensities, whereas all normal tissues were described as having “medium” staining with “moderate” intensities. This suggests that these representative, qualitative immunohistochemical images are consistent with reductions in transcript expression observed in Figure 20A. In addition, single cell map of the kidney (Figure 21B) revealed that ZNF844 expression within the kidney is associated with the cells proximal convoluted tubules and immune cells (macrophages, B-cells and T-cells). Lastly, analysis (Figure 21C) of expression data derived from metastatic and normal tissues (GSE12606) corroborated our TCGA-KIRC data by demonstrating that ZNF844 is significantly down-regulated (adjusted p-Value=0.025928).
ZNF844 gene expression in relationship to clinicopathological features. Table I provides the types of clinicopathological features, the corresponding numbers of each phenotype that were contained in the TCGA KIRC dataset. We further explored ZNF844 in relationship to age, pathological stage, histological grades, and nodal involvement. ZNF844 expression was significantly altered (p<0.001) across stages, grades, and nodal involvement but was not associated with sex (p=0.838) or age (p=0.442) (Table I).
To further characterize ZNF844, we examined the relationship between ZNF844 different clinicopathological features of clear cell renal carcinoma and survival outcome; we conducted multivariate and univariate Cox regression analysis (Table II). In the univariate model, age, stage, grade, and ZNF844 expression were predictive of survival outcome. In the multivariate model, age, stage, and ZNF844 expression retained predictive status, thereby indicating that ZNF844 may be a predictive biomarker for ccRCC.
In terms of ZNF844 expression across stages, ZNF844 transcripts were significantly lower (p<0.001) in all tumor stages relative to normal tissues (Figure 22A). Moreover, ZNF844 was differentially expressed across tumor stages [F(4, 596)=186.314, df=4, p<0.001], with the more advanced stages (3 and 4) having a lower expression of transcripts relative to stages 1 and 2 (Figure 22B). A similar pattern was also noted with regards to histological grades whereby ZNF844 expression was significantly down-regulated across all grades relative to normal tissues [F(4, 571)=216.528, df=4, p<0.001] (Figure 22C) and was significantly lower (p<0.05) in grades (3 and 4) relative to normal, grades 1 and 2 (Figure 22D).
ZNF844 pathway enrichment and gene ontology. Since there is limited information regarding ZNF844 target genes, we examined whether genes (Figure 23A-B) co-expressed with ZNF844 were components of overrepresented pathways. Specifically, we used Gene Set Enrichment Analysis (GSEA) to analyze KEGG,GO Biological Process and Reactome databases to identify enriched pathways (Figure 23C-E). Figure 23C provides the GSEA analysis of the genes co-expressed with ZNF844 using the GO biological processes (GOBP) database. For GSEA GOBP analysis, the most significantly enriched pathways [p<0.05; FDR (q value)<0.05] were found to be associated with inversely correlated genes. The top three enriched pathways were related to immune and translation processes. Immune-related pathways included humoral immune response (GO:0006959, 0.0091878, p<1×10–20, NES=–1.8716), response to interferon-gamma (GO:0034341 FDR=0.040801, p<1×10–20, NES=–1.6651), and adaptive immune response (GO:0002250, FDR=0.048321, p<1×10–20, NES=–1.6470). GSEA analysis of the KEGG database (Figure 23D) also revealed that immune-related and translation/transcription pathways were among the most significantly enriched pathways among inversely correlated genes (p<0.05; FDR<0.05). Immune-related pathways included Epstein-Barr virus infection (hsa05169), Cytokine-cytokine receptor interaction (hsa04060; p<1×10–20, FDR=0.012899, NES=–1.7121), Systemic lupus erythematosus (hsa05322; p<1×10–20, FDR=0.0024186, NES=–1.9218) and Pathogenic Escherichia coli infection (hsa05130; p=<1×10–20, FDR=0.016744, NES=–1.6932). Similar to the GOBP and KEGG databases, GSEA analysis of the Reactome database (Figure 23E) revealed the three most significantly enriched pathways among genes inversely expressed with ZNF844 were also related to transcription and antigen presentation (immune system pathway) with the most significant immune pathway in the Reactome database being antigen-processing cross-presentation (R-HSA-1236975; p<1×10–20, FDR<0.001, NES=–2.1904). GSEA analysis of the GO Biological Processes, KEGG, and Reactome databases did not reveal any significantly enriched pathways among the genes positively correlated with ZNF844 expression.
The ORA method was also used to identify over-represented pathways among negatively correlated genes. The ORA method also yielded similar results whereby 14 out of the 20 topmost enriched pathways are immune related pathways. These included GO:0006955-immune response (FDR<1×10–20, p<1×10–20, enrichment ratio=1.3643), GO:0006952-defense response (FDR£1×10–20, p<1×10–20, enrichment ratio=1.2939), GO:0034097-response to cytokine (FDR£0.001, p<1E-20, enrichment ratio=1.3397), GO:0019221-cytokine-mediated signaling pathway (FDR<1×10–20, p<1×10–20, enrichment ratio=1.4676), GO:0050776: regulation of immune response (FDR=2.22×10–16, p=2.82×10–13, enrichment ratio=1.3634), GO:0045087-innate immune response (FDR=4.32×10–15, p=3.67×10–12, enrichment ratio=1.3683), GO:0045321-leukocyte activation (FDR=4.21×10–15, p=3.67×10–12, enrichment ratio=1.2977), GO:0002252-immune effector process (FDR=2.14×10–14, p=1.61×10–11, enrichment ratio=1.29610), GO:0006954-inflammatory response (FDR=6.39×1014, p=4.28×10–11, enrichment ratio=1.3725) and GO:0006959-humoral immune response (FDR= 8.63×10–14, p=5.37×10–11, enrichment ratio=1.6696). Table III presents the top 20 most enriched pathways based on (FDR value) among inversely correlated genes.
ZNF844 expression and tumor microenvironment. The consistency of both GSEA and ORA methods with regards to the immune-related pathways being highly enriched among negatively correlated genes suggests that ZNF844 may be associated with immune function and tumor microenvironment. More specifically, data indicates that there is an inverse relationship between immune system function and ZNF844 expression, whereby down-regulation of ZNF844 is associated with activation of the immune system. Therefore, we analyzed ZNF844 transcript levels in relationship to immune cell infiltration within the context of the ccRCC tumor microenvironment (Figure 24A and B). We focused on T-cells, B-cells (Figure 24A), and antigen-presenting cells (macrophages and dendritic cells) (Figure 24B) infiltrates because most of the enriched pathways associated with ZNF844 expression were primarily related to adaptive immune system function. ZNF844 showed negligible to weak, negative correlations with CD4-Th2, CD4 memory, CD8 cells, B-cells, and naïve CD4 infiltrates (Figure 24A). Weak to negligible correlations were also noted for antigen-presenting cell infiltrates (Figure 24B). However, ZNF844 expression did show moderate negative correlations (R=–0.558, p=5.15×10–39) with CD4-Th1 infiltrates (Figure 24A). To further explore this observation, we conducted a correlation analysis between ZNF844 expression and several T-cell signatures (Figure 24C). The strongest correlations were seen between ZNF844 and TH1-like T-cells (R= -0.36, p=3.7×10–19) and between ZNF844 and exhausted T-cell signatures (R=–0.37, p=4.1×10–21). This data further suggests a role for ZNF844 in immune function.
Expression profile of ZNF844 in different cancers. Finally, we examined ZNF844 expression across cancer types (Figure 25 and Table IV). ZNF844 was significantly down-regulated (p<0.05) in breast, cervical, colon, esophageal, head and neck, clear cell renal carcinoma, renal papillary, lung, prostate, rectal, thyroid, and uterine cancers (Figure 25). The greatest degree of significance was seen in renal papillary and clear cell carcinomas. In terms of overall survival, changes in ZNF844 expression were only associated with survival in patients with ccRCC and head and neck squamous cell carcinoma (HNSCC) (Table IV).
Discussion
Though poorly understood, it is becoming increasingly clear that KRAB-ZFPs play a major role in the development of cancer and tumor progression (9, 30). The ubiquitous nature of the KRAB-ZFPs, which constitute the largest family of transcription factors, have prevented in-depth analysis of how these proteins contribute or prevent tumorigenesis. KRAB-ZFPs are found on several chromosomes. Clusters of closely related ZFPs are located on chromosomes 4, 7, 8, 10, 19, with the largest clusters located on the p and q arms of chromosome 19 (6, 31-35). We recently reported that ZNF433, a KRAB-ZFP located within a cluster found at chr19p13.2, is down-regulated in ccRCC, and its expression is closely associated with clinicopathological features of ccRCC (19). Given that there are several related KRAB-ZFPs in this region, we examined whether these proteins were also differentially expressed in clear cell renal carcinoma. Seventy percent (14 out of 20) were down-regulated similarly to ZNF433 (Table I). Moreover, decreases in the transcripts of ZNF20, ZNF426, ZNF440, ZNF441, ZNF442, ZNF443, ZNF561, ZNF564, ZNF709, and ZNF844 were associated with lower survival outcomes in ccRCC patients. This data suggests that downregulation of certain chr19.13.2 KRAB-ZFPs may represent a genomic signature for ccRCC.
Of the KRAB-ZFPs, ZNF844 was found to be exceptional because its expression, as well as its impact on survival outcome, had higher statistical relevance than other members of the chr19.13.2 cluster. Therefore, we examined ZNF844 expression in relation to the clinicopathological features of ccRCC. ZNF844 transcripts and proteins were found to be lower in renal carcinomas (See Figure 20, Figure 21A and C). Reductions in ZNF844 expression were closely associated with higher stage, grade, metastasis, as well as with unsatisfactory survival in ccRCC patients. Moreover, ZNF844 was shown to be an independent, moderate prognostic marker for overall survival. This data suggests that elevated or over-expression of ZNF844 has a protective effect against renal carcinomas and may act as a putative tumor suppressor gene.
Gene Set Enrichment Analysis of correlated genes indicated that immune-related pathways were over-represented in an inverse manner (Figure 23C-E), thereby suggesting that ZNF844 is associated with the regulation of the immune system. This supposition is corroborated by earlier reports that clusters of KRAB-ZFPs located on chromosome 19 are highly expressed in human T-lymphocytes (31). Indeed, analysis of immune infiltrates in ccRCC revealed that ZNF844 expression correlates moderately and inversely but in a highly significant manner to TH1 CD4+ T-lymphocytes activity (Figure 24A). Conversely, ZNF844 expression is only negligibly associated with other types of lymphocytes (B-cells) or myeloid lineage cells (Figure 24B). However, Figure 24D does suggest ZNF844 expression is highest in naïve B-cell cells in comparison to other immune cells. This discrepancy may be due to the markers utilized to identify the naïve B-cell subtype across different databases. An additional explanation is that differences in tissue (in this case non-cancerous renal tissue for Figure 24D) could dictate which genes are expressed. Nonetheless, Figure 24D further supports the observations that ZNF844 is associated with immune cell function.
It should be noted that renal carcinoma is highly enriched with immune cells; particularly T-lymphocytes such as the Th1 lineage (36-40). The Th1 lymphocyte phenotype is generally understood to promote anti-tumor immunity and thus is associated with more favorable outcomes in patients with renal carcinomas (40). Conversely, high levels of specific subsets of CD8+, T-regs, and tumor-associated macrophages are associated with unfavorable prognosis in ccRCC (41-45). The high level of T-cell infiltrates is utilized as the basis for the use of PDL1 and checkpoint inhibitors to attenuate; presumably, T-cell inhibition, exhaustion, and anergy, as well as tumor immune evasion in metastatic renal carcinoma (37, 46-49). Our observations indicate that down-regulation of ZNF844 is associated with more advanced disease with poor outcomes and increases in the anti-tumor (Th1) lymphocytic phenotype. These seemingly contradictory results probably suggest that ZNF844 downregulation may be associated with T-cell dysfunction or T-cell exhaustion. In fact, an examination of the relationship between ZNF844 with different T-cell signatures (Figure 24C) revealed that ZNF844 expression had the strongest inverse correlation as well as statistical relevance with the exhausted T-cell phenotype. Moreover, Figure 24D suggest that ZNF844 expression lower in activated lymphocytes, further suggesting that a loss of ZNF844 function is associated with immune cell activation.
Finally, we examined ZNF844 expression in relation to other cancers to determine whether this protein was a possible pan-cancer biomarker (Figure 25 and Table IV). Though the transcription factor was differentially expressed in numerous cancers, the most significant statistical relevance with regards to gene expression and overall survival occurred in ccRCC patients (Table IV). A similar pattern but less statistically robust effect was observed in head and neck squamous cell carcinomas (HNSCC). Hence, ZNF844 prognostic value regarding cancer is primarily limited to ccRCC and possibly HNSCC.
Conclusion
In summary, KRAB-ZFPs gene expression appears to be dysregulated in ccRCC and may represent a genomic signature of ccRCC. We have identified KRAB zinc finger protein ZNF844 as a potential predictive biomarker that may also act as a tumor suppressor in clear cell carcinoma. In silico analysis suggests that down-regulation of ZNF844 is associated with advanced clinicopathological features of ccRCC and poorer patient survival. ZNF844 is more closely associated with the inverse expression of several immune-related pathways, the helper T-cell subtype 1 signature, and the exhausted T-cell phenotype. Our observations therefore suggest that reduction in ZNF844 transcripts correlates with the activation of certain elements of the immune system as the renal carcinoma progresses to advanced stages. However, a weakness of this study is that these observations are made via in silico analysis. Moreover, protein expression data reflects a qualitative assessment of representative immunohistochemistry images. Further investigations in cellular and tumor models are warranted to validate these observations. In addition, functional studies are required to ascertain the role of ZNF844 in ccRCC tumor pathology and to identify the gene targets of this novel transcription factor. Nonetheless, our observations indicate that ZNF844 is indeed a biomarker for ccRCC. Moreover, our study suggests novel roles for KRAB-ZFP in the modulation of cellular functions in ccRCC.
Acknowledgements
This research was supported by the National Institute of Minority Health and Health Disparities of the National Institutes of Health through Grant Number U54 MD 007582 and by the National Cancer Institute of the NIH through Grant Number U54CA233396.
Footnotes
Authors’ Contributions
The Authors contributed to the preparation in the following ways: Conceptualization, SOH, MDS, and RRR; methodology, SOH; resources, KFAS; writing—original draft preparation, SOH; writing— review and editing, KFAS, RRR, MDS; funding acquisition, KFAS. All Authors have read and agreed to the published version of the manuscript.
Conflicts of Interest
The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.
- Received January 12, 2022.
- Revision received February 5, 2022.
- Accepted February 11, 2022.
- Copyright© 2022, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY-NC-ND) 4.0 international license (https://creativecommons.org/licenses/by-nc-nd/4.0).