Chenming Zhang
Sicheng Ma1†- 1The Second Clinical Medical School, Henan University of Chinese Medicine, Zhengzhou, China
- 2Andrology Department, The First Affiliated Hospital of Henan University of Chinese Medicine, Zhengzhou, China
- 3Reproductive Medicine Department, Henan Province Hospital of Traditional Chinese Medicine, Zhengzhou, China
TK1 is overexpressed in numerous cancers and is associated with to a poor prognosis. However, the relationship between methylation status of TK1 and Immune Infiltrates in Prostate Cancer (PCa) is unknown. The goal of this study was to use comprehensive bioinformatic analyses to elucidate the involvement relationship between methylation status of TK1 and Immune Infiltrates in PCa. TK1 mRNA expression and methylation data in PCa were investigated via GEPIA, TIMER, and UALCAN coupled with MEXPRESS data resources. We employed the LinkedOmics data resource to determine the signaling cascades linked to TK1 expression. Single-cell analysis was performed using the CancerSEA data resource. GeneMANIA and CancerSEA were used to analyze the correlation between TK1 and TK1 coexpressed genes. In addition, TIMER and TISIDB were adopted to assess tumor-invading immune cells and immunomodulators. CTD was utilized to detect the drugs acting on TK1. This study found that TK1 was overexpressed in PCa, and its contents were linked to tumor stage and prognosis. Genes co-expressed with TK1 were enriched in cascades involved in the ribosome, cell cycle, oxidative phosphorylation, DNA replication, oocyte meiosis, and the proteasome. The expression of TK1 along with its methylation status was found to be linked to tumor-invading immune cells, as well as PCa immunomodulators. We also examined the prospect of employing TK1 as a possible target for PCa therapy. This work provides the clinical value of TK1 hypermethylation in PCa and highlights new insights into its novel immunomodulatory functions.
Introduction
Prostate cancer (PCa) is one of the most frequent cancers in men, with approximately 1.1 million reported cases worldwide. (Luszczak et al., 2020). Furthermore, PCa has surpassed bladder and kidney cancers in terms of incidence and mortality (Mayer et al., 2001). PCa progresses in a multistep process and is thought to progress from precancerous or intraepithelial neoplasia of the prostate to locally aggressive, metastatic, and hormone-resistant disease (Khan et al., 2008). Several factors can increase the risk of PCa, such as family factors, baldness, gonorrhea, and smoking (Rizzardi et al., 2014; Rao et al., 2015; Ju-Kun et al., 2016), and some unchangeable risk factors include age, race, and family history (Facompre & El-Bayoumy, 2009). Crosstalk between tumor cells the tumor microenvironment (TME) provides new perspectives for understanding the molecular drivers of PCa onset, metastasis, and recurrence (Matos et al., 2021; Melo et al., 2021). However, the molecular pathology of PCa is complex (Xiao et al., 2018). In addition to somatic mutations and chromosomal abnormalities that result in the dysregulation of tumor suppressor genes and oncogenes, changes in the expression of growth factors and their receptors are also involved in the pathogenesis of PCa (Cazares et al., 2010).
Thymidine kinase 1 (TK1) is an enzyme in the DNA salvage pathway, that involved in DNA synthesis and damage through regenerating thymidine. Through facilitated diffusion, thymidine is transferred from the extracellular space to the cell membrane. Where it is converted into its monophosphate form (dTMP) by cytoplasmic TK1 (Bello, 1974; Johnson et al., 1982). As early as the 1960s, the activity of TK1 in tumors was found to be elevated (Weber & Lea, 1966). Numerous studies have further demonstrated that elevated TK1 levels have occur in the serum of many different types of cancers, including some of the most common cancers of the lung cancer, colon cancer, and breast cancer (Alegre et al., 2012; Jagarlamudi & Shaw, 2018; Topolcan & Holubec, 2008; Zhou et al., 2013). Recent findings suggest that intracellular TK1 is correlated with cancer cell invasion and progression, along with its signature role in cancer cell proliferation (Alshabi et al., 2019). The molecular mechanisms of intracellular TK1 overexpression have not been established, and its association with cancer progression has not been sufficiently explored. TK1 overexpression may not only be a by-product of cancer cell processes but rather, part of a selective process that aids cancer cell progression. Tumor growth supported by TK1 has been demonstrated in lung and breast cancer cell lines (Bitter et al., 2020). However, the function of TK1 in PCa remains unclear.
The control and adaptability of virtually all biological processes involve postsynthesis chemical modification of three classes of fundamental macromolecules: DNA, RNA and proteins. One of the most abundant modifications, methylation, is widespread throughout all kingdoms of life, and involves an alkylation reaction whereby a methyl group replaces a hydrogen atom. Methylation is catalysed by methyltransferases (‘writers’), which use S-adenosylmethionine (SAM) as the methyl donor; writers cooperate with dedicated ‘erasers’ (demethylases) and methyl ‘readers’. The importance of the various methylation pathways is highlighted by the fact that their deregulation is linked to many diseases.
Here, we used several available data resources to conduct a comprehensively analyze TK1 in PCa. We found that TK1 was remarkably overexpressed in PCa and that TK1 content was linked to tumor prognosis in individuals with PCa. TK1 coexpressed genes were grouped into cascades involving the cell cycle, oxidative phosphorylation, DNA replication, oocyte meiosis, and the proteasome. In PCa, the content of TK1 and its methylation status were linked to tumor-invading immune cells and immunomodulators. In addition, we explored the possibility of using TK1 as a treatment target for PCa. This findings of this study highlight a novel immunomodulatory function of TK1 hypermethylation in PCa.
Materials and methods
Analysis of gene expression differences
We used the UALCAN data resource (http://ualcan.path.uab.edu/) (Chandrashekar et al., 2017) to study the expression of TK1 mRNA in tumor tissues and in normal tissues and between different nodal metastasis statuses of PCa tissues, the GEPIA data resource (http://gepia2.cancer-pku.cn/) (Tang et al., 2019) was utilized to investigate the differential expression of TK1 in several cancer species with significantly different TKl expression, The HPA data resource (Human Protein Atlas, www.proteinatlas.org) (Thul & Lindskog, 2018) was used to determine the differential protein expression levels of TK1 between PCa and nontumorous tissues.
Prognostic value of TK1 and pathologic parameter correlation
We used the GEPIA data resource to determine the relationship between TK1 expression with overall survival (OS) and disease-free survival (DFS). Single- and multifactor Cox regression analyses for OS and DFS were performed using the survivor package in R (v3.2-10).
Methylation analysis
We analyzed the methylation levels of the TK1 promoter in normal and cancerous tissues at different stages using the UALCAN data resource. The methylation levels of TK1 CpG islands in TK1 in PCa were then analyzed using the MEXPRESS data resource (https://mexpress.be) (Koch et al., 2019).
Coexpressed gene analysis and functional enrichment analyses
We analyzed the list of genes coexpressed with TK1 in PCa using the LinkedOmics data resource (http://www.linkedomics.org/login.php) (Vasaikar et al., 2018). The gene and protein interaction networks of TK1 were constructed using GeneMANIA data resource (http://www.genemania.org) (Warde-Farley et al., 2010) and STRING data resource (https://string-db.org/) (Szklarczyk et al., 2021), respectively. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were also used for the LinkedOmics data. To elucidate the function of TK1 in PCa tissues, we carried out a single-cell analysis using the CancerSEA data resource (http://biocc.hrbmu.edu.cn/CancerSEA/) (Yuan et al., 2019).
Association of TK1 and E2F transcription factors in PCa
Human expression levels of E2Fs in PCa from The Cancer Genome Atlas (TCGA) data resource were estimated using the ggplot2 R package (v3.3.3).
Tumor immune estimation analysis
Immune cell richness was estimated via the GSA R package (v1.34.0) algorithm (Hanzelmann et al., 2013). The correlation module was utilized to establish the relationship between TK1 RNA-seq expression profile data in PCa and immune cells, which included T cells, NK cells, Tfh (T follicular helper) cells, aDC (activated DC), macrophages, Tcm (T central memory) cells, B cells, cytotoxic cells, Th17 cells, DC, iDC (immature DC), mast cells, NK CD56bright cells, eosinophils, pDC (plasmacytoid DC), T helper cells, NK CD56dim cells, Th1 cells, Tem (T effector memory) cells, CD8 T cells, Tgd (T gamma delta) cells, neutrophils, Th2 cells, and Treg cells. The gene markers for immune cells also use gene modules associated with TK1 expression as referenced in previous studies (Bindea et al., 2013).
Association of TK1 methylation with immune infiltrates in PCa
Our study analyzed the relationship between TK1 methylation and immunocytes and immunomodulators by the TISIDB data resource (http://cis.hku.hk/TISIDB) (Ru et al., 2019).
Drug response linked to TK1 expression in PCa
To assess whether TK1 could be a treatment target for individuals with cancer, we obtained a targeted drug for TK1 from the CTD data resource (http://ctdbase.org/) (Davis et al., 2021). We then visualized them via the ggalluvial R package (v3.3.3) algorithm.
Results
Expression of TK1 in PCa
The analysis of the UALCAN data resource further confirmed that expression of TK1 transcript in primary PCa was higher than that in paraneoplastic tissues (Figure 1A). In addition, TK1 presented strong positivity in PCa tissues with different nodal metastasis statuses, and this was in contrast TK1 levels to non-tumorous tissues (Figure 1B). Interactive analysis of the gene expression profiles also verified that TK1 levels were higher in PCa tumor tissues than in paraneoplastic tissues (Figure 1C). The GEPIA data resource confirmed these results (Figure 1D).
FIGURE 1. Overexpression of TK1 in PCa. (A) TK1 expression levels in various human malignancies were analyzed by TIMER. (B) UALCAN analyzed TK1 mRNA levels in tumor tissues and nontumorous tissues. (C) The expression of TK1 mRNA between different nodal metastasis statuses of PCa tissues was analyzed by UALCAN. (D) Differential protein expression of TK1 in tumor and nontumorous tissues in the HPA. (E) GEPIA analyzed high or low expression levels of TK1 in various malignant tumors.
The prognostic relevance of TK1 and TK1’s relationship with clinicopathological parameters in PCa
The baseline features of the patients are given in Table 1. Among the patients, 499 had both the clinical and genetic characteristics of PCa. The GEPIA assessment illustrated that high TK1 expression was significantly linked to poor OS (HR = 8.2, p = 0.046) and DFS (HR = 3, p = 3.2e-06) in individuals with PCa (Figures 2A,B). The univariate Cox data showed that OS was linked to the M stage, Gleason score, primary therapy outcome, and TK1 expression. (Table 2). Furthermore, the univariate Cox regression analysis showed that T stage, Gleason score, N stage, primary therapy outcome, and TK1 expression were associated with PFS. Lastly, the multivariate Cox regression revealed that the Gleason score, TK1 content, and primary treatment outcome were prognostic parameters in individuals with PCa (Table 3). These data indicate that there is an association between TK1 and a worsening prognosis.
TABLE 1. Relationship between TK1 expression and clinicopathological features in patients with PCa.
FIGURE 2. Prognostic value as a result of TK1 in PCa patients. (A,B) GEPIA profiling showed that PCa patients with high TK1 expression were associated with poorer OS. (A) and PFS (B) relative to those with low TK1 expression.
TABLE 2. Univariate and multivariable Cox regression of TK1 expression for overall survival in PCa patients.
TABLE 3. Univariate and multivariable Cox regression of TK1 expression for progression-free survival in PCa patients.
DNA methylation of TK1 and PCa prognosis
TK1 expression was negatively correlated with the methylation of CpG islands and promoter regions. In PCa cancer tissues, the DNA methylation level of TK1 was remarkably lower than that in non-tumorous samples. (Figure 3A). Furthermore the DNA methylation level of N1 was lower than that in N0 in PCa tissues with different nodal metastasis statuses, while the DNA methylation level of N0 was lower than that in nontumorous tissues (Figure 3B). The expression of TK1 was negatively correlated with methylation around CpG islands and promoter regions (Figure 3C). Thirteen aberrant methylation sites (cg10456035, r = -0.255; cg20740903, r = -0.334; cg25069807, r = -0.386; cg21519872, r = -0.414; cg08115732, r = -0.264; cg00715343, r = -0.167; cg18767057, r = -0.137; cg03291825, r = -0.361; cg26206461, r = -0.101; cg07314523, r = -0.209; cg22061523, r = -0.384; cg02441982, r = -0.335; cg21940220, r = -0.162), and the observed Pearson correlation coefficients are indicated on the right. Here, DNA methylation of the TK1 promoter may further support the TK1 overexpression in PCa.
FIGURE 3. DNA methylation level of TK1 in PCa. (A) UALCAN data resource showing promoter methylation levels of TK1 in normal and PCa tissues. (B) UALCAN data resource offering promoter methylation levels of TK1 in PCa tissues with different nodal metastasis statuses. (C) Relationship between TK1 expression and PCa and methylation around CpG islands and promoter regions.
Gene and functional enrichment analysis of TK1 coexpression
Figure 4 Ashows the coexpression trend of TK1. The heat map contains the top 50 positively correlated genes along with the negative genes linked to TK1 (Figures 4B,C). The genes along with the protein-protein interaction networks created via GeneMANIA and STRING revealed that 20 potential target genes (DTYMK, LRRK2, CRMP1, DDR1, CCND1, BIRC5, MAMDC2, CEBPA, PCLAF, E2F1/3/4/6, RELA, CRIP1 RRM2, CDK1, CDKN3, MELK, RAD51) and ten potential target proteins (PNP, TYMP, CDA, DCTD, UPP1, DTYMK, TYMS, DUT, BIRC5, E2F1) that interacted with TK1 (Figures 4D,E). To further characterize the function of TK1 in PCa, we conducted a single-cell analysis using CancerSEA. The results illustrated that TK1 was mainly involved in the cell cycle, DNA repair, and apoptosis (Figure 4F). The annotation of the GO terms showed that the co-expressed genes of TK1 mainly participated in chromosome segregation, organelle fission, DNA replication, mitotic cell cycle phase transition, meiotic cell cycle, DNA conformation change, spindle organization, chromosome localization, the cell cycle G2/M phase transition, harmful modulation of the cell cycle process, mitochondrial respiratory chain complex assembly, translational elongation, positive modulation of the cell cycle, the cell cycle G1/S phase transition, protein localization to chromosome, negative regulation of mitotic cell cycle, NADH dehydrogenase complex assembly, signal transduction in response to DNA damage, etc. (Figure 4G). The KEGG pathway analysis indicated enrichment in the ribosome, cell cycle, oxidative phosphorylation, DNA replication, oocyte meiosis, proteasome, pyrimidine metabolism, base excision repair, Huntington disease, Parkinson disease, homologous recombination, p53 signaling cascade, nonalcoholic fatty liver disease, progesterone-mediated oocyte maturation, spliceosome, mismatch repair, Fanconi anemia cascade, nucleotide excision repair, etc. (Figure 4H). Single-cell analysis of the functional characteristics of TK1 in prostate cancer showed that TK1 expression positively correlated with Cell cycle (p = 0.001), DNA repair (p = 0.003), Apoptosis (p = 0.005) Figure 5
FIGURE 4. Gene and functional enrichment analysis of TK1 coexpression. (A) Volcano map of TK1 coexpression in PCa according to the LinkedOmics data resource. (B,C) Heatmap of the top 50 genes positively (B) and negatively (C) linked to TK1 are shown. (D) Gene network associated with TK1, analysis by GeneMANIA. (E) Interaction network map between proteins encoded by TK1 genes plotted with STRING. (F) Single-cell analysis of the CancerSEA data resource indicates that TK1 is mainly involved in the cell cycle, DNA repair, and apoptosis. (G,H) TK1 coexpressed genes were annotated on LinkedOmics via Gene Ontology (GO) assessment (G) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis (H).
FIGURE 5. Single-cell analysis of the functional characteristics of TK1 in prostate cancer. (A) Correlation between TK1 and 14 cell functions. (B) The screened functional state with significant indigenous correlation with TK1.
Association of TK1 and E2F transcription factors in PCa
The association between TK1 and E2F transcription factors was investigated by using a ggplot2 R package (v3.3.3) algorithm analysis (Figure 6). E2F1 (r = 0.818, p = 0.001), E2F2 (r = 0.766, p = 0.001), E2F3 (r = 0.183, p = 0.05), E2F4 (r = 0.081, p = 0.072), E2F5 (r = 0.251, p = 0.001), E2F7 (r = 0.497, p = 0.001), and E2F8 (r = 0.526, p = 0.001) were all positively linked to TK1 expression. E2F6 (r = -0.121, p = 0.007) was negatively correlated with TK1 expression (Figure 6A). The results of the ggplot2 package (v3.3.3) algorithm analysis further confirmed that E2F1, E2F2, E2F3 and E2F5 expression were higher in primary PCa cancer tissues than in precancer tissues; however, E2F4 and E2F6 were lower in primary PCa tissues than in precancerous tissues (Figure 6B).
FIGURE 6. Association of TK1 and E2Fs in PCa. (A) Molecular correlation analysis of TK1 and E2Fs. (B) Human expression levels of E2Fs in PCa from The Cancer Genome Atlas (TCGA) data resource were estimated by the ggplot2 R package (v3.3.3) algorithm. (*<0.05, **<0.01, ***<0.001).
Relationship of TK1 with immune cell invasion in PCa
The infiltration of 24 types of immune cells into PCa was first determined using the ssGSEA method, followed by a Spearman analysis to assess the relationship between TK1 with immune cell invasion (Figure 7A). Th2 cells (r = 0.466, p = 0.001), TReg cells (r = 0.162, p = 0.001), and NK CD56bright cells (r = 0.148, p = 0.05) were all positively linked to TK1 contents. Nonetheless, Mast cells (r = -0.307, p = 0.001), eosinophils (r = -0.294, p = 0.001), NK cells (r = -0.286, p = 0.001), Neutrophils (r = -0.237, p = 0.001), Tcm (r = -0.202, p = 0.001), Tgd (r = -0.192, p = 0.001) (Figure 7B). In addition, TK1 was remarkably linked to the gene signatures of monocytes, M1 macrophages, dendritic cells, Th1 cells, Th2 cells, and Tregs (Table 4).
FIGURE 7. Effect of immune cell invasion on TK1 expression in individuals with PCa. (A) Relationship between the invasion contents of 24 immune cell types and the expression profile of TK1. (B) Comparison of the most associated levels of immune cell invasion in the high and low TK1 expression groups, including eosinophils, Tcm, neutrophils, mast cells, Tgd, NK cells, Th2 cells, TRed cells, and NK CD56bright cells. (*<0.05, **<0.01, ***<0.001).
TABLE 4. The correlations between TK1 and gene markers of immune cells in PCa and normal cells by TIMER.
Association of TK1 methylation with immune infiltrates in PCa
To assess the effect of TK1 methylation on PCa progression, we evaluated the relationship between TK1 methylation and immune cell invasion using the TISIDB platform. The results showed that TK1 methylation status correlated with Th2 cells (r = 0.418, p = 2.2e-16), TReg cells (r = 0.448, p = 2.2e-16), mast cells (r = 0.497, p = 2.25e-32), eosinophils (r = 0.458, p = 2.2e-16), NK cells (r = 0.567, p = 2.2e-16), Tcm cells (r = 0.563, p = 2.2e-16), and Tgd cells (r = 0.324, p = 1.67e-13) (Figure 8A). Similarly, the methylation status of TK1 was also linked to immunostimulants, such as C10orf54 (r = 0.605, p = 2.2e-16), CD40 (r = 0.552, p = 2.2e-16), TMEM173 (r = 0.557, p = 2.2e-16), and NT5E (r = 0.647, p = 2.2e-16) (Figure 7B).
FIGURE 8. The relationship of the methylation status of TK1 with immune infiltration in PCa. (A) Relationship of the methylation status of TK1 with Th2 cells, Tregs, NK eosinophils, NK cells, neutrophils, Tcm, and Tgd in PCa available from the TISIDB data resource. (B) Relationship of the methylation status of TK1 with immunostimulator PCa available from the TISIDB data resource.
Drug responses related to TK1 expression in PCa
As shown in Figure 9, there are 38 drugs that are helpful in the treatment of prostate cancer with high TK1 expression (two drugs can affect TK1 cotreatment, 24 drugs can increase TK1’s mutagenesis, two drugs affect the response to the substance of TK1, three drugs decrease TK1 reaction, one drug increase the response to TK1, four drugs can increase TK1 methylation, and one drug decrease metabolic processing).
FIGURE 9. Drug response related to TK1 expression in PCa.
Discussion
PCa is the most commen cancer threat in men, and biochemical failure in a short time after radical treatment is linked to a high risk of tumor progression (Hernandez et al., 2019). Thymidine kinase 1 (TK1) is a key enzyme involved in synthesis DNA precursors. It is upregulated in the S phase of the cell cycle, and its presence in cells is an indicator of active cell multiplication (Chen et al., 2010). TK1 levels are dramatically increased when DNA is damaged by radiation or chemotherapy drugs (Haveman et al., 2006; Chen et al., 2010). This may improve the restoration efficiency of thymidine released during DNA breakdown. TK1 is involved in DNA repair, and the impairment or low expression of TK1 during DNA restoration may lead to cell death (Chen et al., 2010). Collectively, these findings reveal that TK1 may be involved in the development of PCa in various ways and that repressing its functional site may be a novel approach for treating PCa.
Herein, the contents of TK1 in PCa and its relevance to the diagnosis, along with the prognosis of PCa, came to light. TK1 is dysregulated in several cancer types and plays a vital role in cancer onset and progression (Topolcan & Holubec, 2008). Based on TCGA data analysis, PCa shows remarkably high TK1 mRNA expression. Sun R. et al. Reported that TK1 expression and activity are affected by oxidative stress and nucleoside analog treatment (Sun et al., 2020). Researchers have also proven that the TK1 gene promoter’s abnormal CpG island methylation is involved in controlling TK1 expression in failing and non-failing human left ventricular myocardium (Koczor et al., 2013a).
Consequently, we then explored the promoter methylation levels of TK1 in PCa. We found that the DNA methylation level of TK1 in cancer tissues was markedly lower than that in non-tumorous tissues, which may suggests that low TK1 promoter methylation levels are due to the overexpression of TK1 in PCa. In addition, we found that TK1 expression was remarkably positively linked to the expression of the E2F gene family in PCa, with E2F1/2/3/5 being highly expressed in PCa. Several prior studies have shown that the transcription factor E2F1 can directly regulate TK1 and promote tumor proliferation (Zhu et al., 2018). We speculate that E2F2/3/5 have a similar effects.
Recently, it was shown that despite being a biomarker, TK1 also participates in cancer invasion and progression (Tilli et al., 2016; Alshabi et al., 2019; Song et al., 2019). TK1 has been previously documented to enhance cancer progression, and high TK1 levels are predictive of poor prognosis in individuals with cancer. To determine whether TK1 is an independent predictive factor for PCa, we studied the prognosis of patients with PCa with diverse TK1 levels. We found that high a TK1 content in PCa tissues was linked to poor prognosis in individuals with PCa. A multivariate Cox regression further supports that TK1 overexpression is a factor of poor prognosis in PCa patients, showing that TK1 is a prognostic biomarker for PCa. Genes in the same subpopulation tend to be coexpressed and synergistically coregulated.
To reveal the biological function of TK1, we performed co-expression and functional enrichment assessments. In addition, to investigate the function of TK1 in PCa, we performed a single-cell analysis. The results revealed that TK1 primarily participates in the cell cycle, DNA repair, and apoptosis. The dynamic crosstalk of tumor cells with the TME may trigger a chronic inflammatory environment that drives the onset and progression of cancer (Fridlender et al., 2009; Mehla & Singh, 2019).
To establish the role of TK1 in the TME, the relationship between TK1 and immune invasion in PCa was investigated. Our findings suggest that TK1 expression positively correlates with the immune invasion of immune cell populations consisting of Th2 cells, Treg cells, and NK CD56bright cells. However, TK1 expression was negatively linked to mast cells, eosinophils, NK cells, neutrophils, Tcm, and Tgd. The relationship between TK1 levels and genetic markers of immune cells further showed that TK1 interacts with Th2 and Treg cells. In the PCa microenvironment of PCa, high levels of Treg cells were associated with poor prognosis (Shen et al., 2021). These data suggest that TK1 may have a sophisticated modulatory role in immune-related processes.
The dysregulation of DNA methylation in the epigenome affects the immunogenicity of tumors and immune cells in the TME (Hogg et al., 2020). Our study found that the hypomethylation status of TK1 tended to be more common in high-grade and advanced PCa, which may illustrate that a change in TK1 methylation patterns inhibits PCa progression. The possible reason for this may be because the hypermethylation status of TK1 is detrimental to cell proliferation (Koczor et al., 2013b). To determine the mechanism by which TK1 methylation inhibits PCa progression, we explored the relationship between TK1 methylation and immune invasion. Our findings illustrated that TK1 methylation was positively linked to immune cells and immunostimulatory factors. The hypermethylation status of TK1 may contribute to the inhibition of tumor progression in PCa, which helps explain the hypomethylation status of TK1 in PCa. Methylation of TK1 may indicate cancer immune infiltration and response to immunotherapeutic agents in patients with PCa. It may be a potential predictor of PCa patients’ response to immunotherapeutic agents. The methylation level of TK1 was shown to act as a valid prognostic biomarker for PCa.
Based on clinical and pathophysiological data, TK1 is a promising therapeutic target for PCa. There are 38 drugs that are useful for treating PCa with high TK1 expression. Interestingly, we found that multiple drugs acted on TK1 methylation. Therefore, we hypothesize that the PCa development could be blocked by drugs that induce TK1 hypermethylation.
In summary, this study showed that TK1 is a prognostically relevant marker for PCa. TK1 expression and methylation status are linked to immune cell invasion and immune modulation. However, this study had certain limitations. First, gene expression analysis based on open source data resources may have limited accuracy. This possibility requires further investigation using in vitro and in vivo models to assess TK1 and the possible biological mechanisms of tumor immune crosstalk in PCa. Second, the impact of the E2F family targets on TK1 need to be further explored. Third, the TK1 gene methylation changes that occur during PCa development need to be further investigated. In summary, our study highlights a novel immunomodulatory function of TK1 hypermethylation in PCa.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
CZ and SM: Methodology, Project administration, Data curation, Writing-Original draft preparation, Validation, Visualization, Conceptualization, Funding acquisition; XH: Resources, Software; ZW: Validation; ZS: Reviewing and Editing, Funding acquisition. The manuscript was read and approved by all authors.
Funding
This work was supported by the National Natural Science Foundation of China (NSFC) (82174377, 81974573, and 81603632); the Natural Science Foundation of Henan Province (222300420485); the Special Scientific Research Project of National TCM Clinical Research Base (2021JDZX 2056); the University Student of Henan province Innovation and Entrepreneurship Training Program Project (202210471022 and 202210471057); Science and Technology Program of Henan Province (192102310159).
Acknowledgments
All the manuscript authors are immensely grateful to the foundations for their valuable support.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Alegre, M. M., Robison, R. A., and O'Neill, K. L. (2012). Thymidine kinase 1 upregulation is an early event in breast tumor formation. J. Oncol. 2012, 575647. doi:10.1155/2012/575647
Alshabi, A. M., Vastrad, B., Shaikh, I. A., and Vastrad, C. (2019). Identification of important invasion and proliferation related genes in adrenocortical carcinoma. Med. Oncol. 36 (9), 73. doi:10.1007/s12032-019-1296-7
Bello, L. J. (1974). Regulation of thymidine kinase synthesis in human cells. Exp. Cell Res. 89 (2), 263–274. doi:10.1016/0014-4827(74)90790-3
Bindea, G., Mlecnik, B., Tosolini, M., Kirilovsky, A., Waldner, M., Obenauf, A. C., et al. (2013). Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity 39 (4), 782–795. doi:10.1016/j.immuni.2013.10.003
Bitter, E. E., Townsend, M. H., Erickson, R., Allen, C., and O'Neill, K. L. (2020). Thymidine kinase 1 through the ages: A comprehensive review. Cell Biosci. 10 (1), 138. doi:10.1186/s13578-020-00493-1
Cazares, L. H., Drake, R. R., Esquela-Kirscher, A., Lance, R. S., Semmes, O. J., and Troyer, D. A. (2010). Molecular pathology of prostate cancer. Cancer Biomark. 9 (1-6), 441–459. doi:10.3233/CBM-2011-0181
Chandrashekar, D. S., Bashel, B., Balasubramanya, S. A. H., Creighton, C. J., Ponce-Rodriguez, I., Chakravarthi, B., et al. (2017). Ualcan: A portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia 19 (8), 649–658. doi:10.1016/j.neo.2017.05.002
Chen, Y. L., Eriksson, S., and Chang, Z. F. (2010). Regulation and functional contribution of thymidine kinase 1 in repair of DNA damage. J. Biol. Chem. 285 (35), 27327–27335. doi:10.1074/jbc.M110.137042
Davis, A. P., Grondin, C. J., Johnson, R. J., Sciaky, D., Wiegers, J., Wiegers, T. C., et al. (2021). Comparative toxicogenomics database (CTD): Update 2021. Nucleic Acids Res. 49 (D1), D1138–D1143. doi:10.1093/nar/gkaa891
Facompre, N., and El-Bayoumy, K. (2009). Potential stages for prostate cancer prevention with selenium: Implications for cancer survivors. Cancer Res. 69 (7), 2699–2703. doi:10.1158/0008-5472.CAN-08-4359
Fridlender, Z. G., Sun, J., Kim, S., Kapoor, V., Cheng, G., Ling, L., et al. (2009). Polarization of tumor-associated neutrophil phenotype by TGF-beta: "N1" versus "N2" TAN. Cancer Cell 16 (3), 183–194. doi:10.1016/j.ccr.2009.06.017
Hanzelmann, S., Castelo, R., and Guinney, J. (2013). Gsva: Gene set variation analysis for microarray and RNA-seq data. BMC Bioinforma. 14, 7. doi:10.1186/1471-2105-14-7
Haveman, J., Sigmond, J., van Bree, C., Franken, N. A., Koedooder, C., Peters, G. J., et al. (2006). Time course of enhanced activity of deoxycytidine kinase and thymidine kinase 1 and 2 in cultured human squamous lung carcinoma cells, SW-1573, induced by gamma-irradiation. Oncol. Rep. 16 (4), 901–905. doi:10.3892/or.16.4.901
Hernandez, G., Ramirez, J. L., Pedroza-Torres, A., Herrera, L. A., and Jimenez-Rios, M. A. (2019). The secret life of translation initiation in prostate cancer. Front. Genet. 10, 14. doi:10.3389/fgene.2019.00014
Hogg, S. J., Beavis, P. A., Dawson, M. A., and Johnstone, R. W. (2020). Targeting the epigenetic regulation of antitumour immunity. Nat. Rev. Drug Discov. 19 (11), 776–800. doi:10.1038/s41573-020-0077-5
Jagarlamudi, K. K., and Shaw, M. (2018). Thymidine kinase 1 as a tumor biomarker: Technical advances offer new potential to an old biomarker. Biomark. Med. 12 (9), 1035–1048. doi:10.2217/bmm-2018-0157
Johnson, L. F., Rao, L. G., and Muench, A. J. (1982). Regulation of thymidine kinase enzyme level in serum-stimulated mouse 3T6 fibroblasts. Exp. Cell Res. 138 (1), 79–85. doi:10.1016/0014-4827(82)90093-3
Ju-Kun, S., Yuan, D. B., Rao, H. F., Chen, T. F., Luan, B. S., Xu, X. M., et al. (2016). Association between Cd exposure and risk of prostate cancer: A PRISMA-compliant systematic review and meta-analysis. Med. Baltim. 95 (6), e2708. doi:10.1097/MD.0000000000002708
Khan, N., Afaq, F., and Mukhtar, H. (2008). Cancer chemoprevention through dietary antioxidants: Progress and promise. Antioxid. Redox Signal. 10 (3), 475–510. doi:10.1089/ars.2007.1740
Koch, A., Jeschke, J., Van Criekinge, W., van Engeland, M., and De Meyer, T. (2019). MEXPRESS update 2019. Nucleic Acids Res. 47 (W1), W561–W565. doi:10.1093/nar/gkz445
Koczor, C. A., Lee, E. K., Torres, R. A., Boyd, A., Vega, J. D., Uppal, K., et al. (2013a). Detection of differentially methylated gene promoters in failing and nonfailing human left ventricle myocardium using computation analysis. Physiol. Genomics 45 (14), 597–605. doi:10.1152/physiolgenomics.00013.2013
Koczor, C. A., Torres, R. A., Fields, E. J., Boyd, A., He, S., Patel, N., et al. (2013b). Thymidine kinase and mtDNA depletion in human cardiomyopathy: Epigenetic and translational evidence for energy starvation. Physiol. Genomics 45 (14), 590–596. doi:10.1152/physiolgenomics.00014.2013
Luszczak, S., Kumar, C., Sathyadevan, V. K., Simpson, B. S., Gately, K. A., Whitaker, H. C., et al. (2020). PIM kinase inhibition: Co-targeted therapeutic approaches in prostate cancer. Signal Transduct. Target. Ther. 5 (1), 7. doi:10.1038/s41392-020-0109-y
Matos, A., Carvalho, M., Bicho, M., and Ribeiro, R. (2021). Arginine and arginases modulate metabolism, tumor microenvironment and prostate cancer progression. Nutrients 13 (12), 4503. doi:10.3390/nu13124503
Mayer, R., Klemen, H., Quehenberger, F., Sankin, O., Mayer, E., Hackl, A., et al. (2001). Hyperbaric oxygen--an effective tool to treat radiation morbidity in prostate cancer. Radiother. Oncol. 61 (2), 151–156. doi:10.1016/s0167-8140(01)00430-3
Mehla, K., and Singh, P. K. (2019). Metabolic regulation of macrophage polarization in cancer. Trends Cancer 5 (12), 822–834. doi:10.1016/j.trecan.2019.10.007
Melo, C. M., Vidotto, T., Chaves, L. P., Lautert-Dutra, W., Reis, R. B. D., Squire, J. A., et al. (2021). The role of somatic mutations on the immune response of the tumor microenvironment in prostate cancer. Int. J. Mol. Sci. 22 (17), 9550. doi:10.3390/ijms22179550
Rao, D., Yu, H., Bai, Y., Zheng, X., and Xie, L. (2015). Does night-shift work increase the risk of prostate cancer? A systematic review and meta-analysis. Onco. Targets. Ther. 8, 2817–2826. doi:10.2147/OTT.S89769
Rizzardi, A. E., Vogel, R. I., Koopmeiners, J. S., Forster, C. L., Marston, L. O., Rosener, N. K., et al. (2014). Elevated hyaluronan and hyaluronan-mediated motility receptor are associated with biochemical failure in patients with intermediate-grade prostate tumors. Cancer 120 (12), 1800–1809. doi:10.1002/cncr.28646
Ru, B., Wong, C. N., Tong, Y., Zhong, J. Y., Zhong, S. S. W., Wu, W. C., et al. (2019). Tisidb: An integrated repository portal for tumor-immune system interactions. Bioinformatics 35 (20), 4200–4202. doi:10.1093/bioinformatics/btz210
Shen, Y., Xu, H., Long, M., Guo, M., Li, P., Zhan, M., et al. (2021). Screening to identify an immune landscape-based prognostic predictor and therapeutic target for prostate cancer. Front. Oncol. 11, 761643. doi:10.3389/fonc.2021.761643
Song, Z. Y., Chao, F., Zhuo, Z., Ma, Z., Li, W., Chen, G., et al. (2019). Identification of hub genes in prostate cancer using robust rank aggregation and weighted gene co-expression network analysis. Aging (Albany NY) 11 (13), 4736–4756. doi:10.18632/aging.102087
Sun, R., Eriksson, S., and Wang, L. (2020). The expression and activity of thymidine kinase 1 and deoxycytidine kinase are modulated by hydrogen peroxide and nucleoside analogs. Nucleosides Nucleotides Nucleic Acids 39 (10-12), 1347–1358. doi:10.1080/15257770.2020.1720234
Szklarczyk, D., Gable, A. L., Nastou, K. C., Lyon, D., Kirsch, R., Pyysalo, S., et al. (2021). Correction to 'the STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 49 (18), 10800. doi:10.1093/nar/gkab835
Tang, Z., Kang, B., Li, C., Chen, T., and Zhang, Z. (2019). GEPIA2: An enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 47 (W1), W556–W560. doi:10.1093/nar/gkz430
Thul, P. J., and Lindskog, C. (2018). The human protein atlas: A spatial map of the human proteome. Protein Sci. 27 (1), 233–244. doi:10.1002/pro.3307
Tilli, T. M., Carels, N., Tuszynski, J. A., and Pasdar, M. (2016). Validation of a network-based strategy for the optimization of combinatorial target selection in breast cancer therapy: siRNA knockdown of network targets in MDA-MB-231 cells as an in vitro model for inhibition of tumor development. Oncotarget 7 (39), 63189–63203. doi:10.18632/oncotarget.11055
Topolcan, O., and Holubec, L. (2008). The role of thymidine kinase in cancer diseases. Expert Opin. Med. diagn. 2 (2), 129–141. doi:10.1517/17530059.2.2.129
Vasaikar, S. V., Straub, P., Wang, J., and Zhang, B. (2018). LinkedOmics: Analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 46 (D1), D956–D963. doi:10.1093/nar/gkx1090
Warde-Farley, D., Donaldson, S. L., Comes, O., Zuberi, K., Badrawi, R., Chao, P., et al. (2010). The GeneMANIA prediction server: Biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 38, W214–W220. Web Server issue). doi:10.1093/nar/gkq537
Weber, G., and Lea, M. A. (1966). The molecular correlation concept of neoplasia. Adv. Enzyme Regul. 4, 115–145. doi:10.1016/0065-2571(66)90011-2
Xiao, J., Cohen, P., Stern, M. C., Odedina, F., Carpten, J., Reams, R., et al. (2018). Mitochondrial biology and prostate cancer ethnic disparity. Carcinogenesis 39 (11), 1311–1319. doi:10.1093/carcin/bgy133
Yuan, H., Yan, M., Zhang, G., Liu, W., Deng, C., Liao, G., et al. (2019). CancerSEA: A cancer single-cell state atlas. Nucleic Acids Res. 47 (D1), D900–D908. doi:10.1093/nar/gky939
Zhou, J., He, E., and Skog, S. (2013). The proliferation marker thymidine kinase 1 in clinical use. Mol. Clin. Oncol. 1 (1), 18–28. doi:10.3892/mco.2012.19
Keywords: prostate cancer, thymidine kinase 1, prognostic biomarker, immune infiltrates, DNA methylation
Citation: Zhang C, Ma S, Hao X, Wang Z and Sun Z (2022) Methylation status of TK1 correlated with immune infiltrates in prostate cancer. Front. Genet. 13:899384. doi: 10.3389/fgene.2022.899384
Received: 18 March 2022; Accepted: 07 July 2022;
Published: 11 August 2022.
Edited by:
Valentina Silvestri, Sapienza University of Rome, ItalyReviewed by:
Pradeep Reddy Cingaram, University of Miami, United StatesSohini Chakraborty, New York University, United States
Copyright © 2022 Zhang, Ma, Hao, Wang and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chenming Zhang, emNtaG56eUAxNjMuY29t; Zixue Sun, c3VuaGh6eEAxNjMuY29t
†These authors have contributed equally to this work and share first authorship