- 1Department of General Surgery, Tianjin Fifth Central Hospital, Tianjin, China
- 2Tianjin Key Laboratory of Epigenetics for Organ Development of Premature Infants, Tianjin Fifth Central Hospital, Tianjin, China
- 3High Altitude Characteristic Medical Research Institute, Huangnan Tibetan Autonomous Prefecture People’s Hospital, Huangnan Prefecture, Qinghai, China
- 4Department of Pathology, Tianjin Fifth Central Hospital, Tianjin, China
- 5North China University of Science and Technology, Tangshan, Hebei, China
- 6Department of Ophthalmology, Peking University First Hospital, Beijing, China
- 7Department of Pharmacy, Tianjin Fifth Central Hospital, Tianjin, China
- 8Department of General Surgery, Peking University First Hospital, Beijing, China
The Ten-Eleven Translocation (TET) family genes are implicated in a wide array of biological functions across various human cancers. Nonetheless, there is a scarcity of studies that comprehensively analyze the correlation between TET family members and the molecular phenotypes and clinical characteristics of different cancers. Leveraging updated public databases and employing several bioinformatics analysis methods, we assessed the expression levels, somatic variations, methylation levels, and prognostic values of TET family genes. Additionally, we explored the association between the expression of TET family genes and pathway activity, tumor microenvironment (TME), stemness score, immune subtype, clinical staging, and drug sensitivity in pan-cancer. Molecular biology and cytology experiments were conducted to validate the potential role of TET3 in tumor progression. Each TET family gene displayed distinct expression patterns across at least ten detected tumors. The frequency of Single Nucleotide Variant (SNV) in TET genes was found to be 91.24%, primarily comprising missense mutation types, with the main types of copy number variant (CNV) being heterozygous amplifications and deletions. TET1 gene exhibited high methylation levels, whereas TET2 and TET3 genes displayed hypomethylation in most cancers, which correlated closely with patient prognosis. Pathway activity analysis revealed the involvement of TET family genes in multiple signaling pathways, including cell cycle, apoptosis, DNA damage response, hormone AR, PI3K/AKT, and RTK. Furthermore, the expression levels of TET family genes were shown to impact the clinical staging of tumor patients, modulate the sensitivity of chemotherapy drugs, and thereby influence patient prognosis by participating in the regulation of the tumor microenvironment, cellular stemness potential, and immune subtype. Notably, TET3 was identified to promote cancer progression across various tumors, and its silencing was found to inhibit tumor malignancy and enhance chemotherapy sensitivity. These findings shed light on the role of TET family genes in cancer progression and offer insights for further research on TET3 as a potential therapeutic target for pan-cancer.
Introduction
Epigenetics, including DNA methylation, histone modifications, and non-coding RNA regulation, plays a crucial role in human disease, particularly cancer, by orchestrating gene expression dynamics without altering DNA sequence. Aberrant epigenetic changes, such as hypermethylation of tumor suppressor gene promoters and dysregulated histone modifications, disrupt normal cellular processes, fostering tumor initiation, progression, and metastasis. Despite the continuous advancements in surgical and targeted therapies in recent years, a considerable portion of tumor patients still present with advanced and widespread metastasis at the time of diagnosis, resulting in a bleak prognosis (Golemis et al., 2018). Consequently, conducting comprehensive research into the molecular mechanisms underlying cancer development holds paramount importance for tumor prevention, early screening and diagnosis, clinical treatment, as well as enhancing patient prognosis and survival rates (Alqahtani et al., 2019; Ban et al., 2021; Rodriguez et al., 2022).
Genetic variations and epigenetic changes play pivotal roles in the initiation and progression of tumors (Talib et al., 2021). Nearly all tumors involve single or tandem mutations in one or multiple genes (Kossenas and Constantinou, 2021). To devise precise treatment strategies capable of effectively targeting tumor cells while minimizing damage to normal tissues, from the myriad of small molecule chemical drugs available, a comprehensive understanding of the underlying mechanisms of the implicated genes is imperative (Srivastava and Goodwin, 2020). The advent and advancement of modern bioinformatics methodologies offer a streamlined approach for humanity to efficiently explore and analyze tumor occurrences within established gene information databases, thereby providing guidance for tumor treatment (Nooter and Stoter, 1996; Zaridze, 2008; Gonçalves et al., 2021).
Recent studies have revealed that epigenetic regulations, such as DNA methylation and demethylation modifications, are intricately linked to the onset and progression of human cancers (Kulis and Esteller, 2010; Klutstein, Nejman, Greenfield and Cedar, 2016; Nishiyama and Nakanishi, 2021). DNA methylation often occurs on the CpG islands of the promoter regions, resulting in the activation of oncogenes and the inactivation of tumor suppressor genes (Kulis and Esteller, 2010; Papanicolau-Sengos and Aldape, 2022). Aberrant DNA methylation can disrupt the apoptosis process of cells, render tumor cells insensitive to growth inhibition signals, and induce uncontrolled replication of tumor cells (Pan et al., 2018). Consequently, epigenetic research, particularly focusing on DNA methylation, holds significant value for tumor molecular diagnosis, prevention, and treatment, as well as predicting therapeutic outcomes and prognosis (C. Chen et al., 2022; Meng et al., 2015).
Studies indicate that both hypermethylation and hypomethylation play crucial roles in tumor development. The Ten-Eleven Translocation (TET) enzymes facilitate the oxidation of 5-methylcytosines (5mCs) and promote site-specific reversal of DNA methylation (Kohli and Zhang, 2013). TET family genes, including TET1, TET2, and TET3, exhibit diverse vital biological functions in mammals (Kinney and Pradhan, 2013; Zeng and Chen, 2019; Joshi, Liu, Breslin and Zhang, 2022). At their core, TET family genes consist of α-Ketoglutaric acid (α-KG) and Fe2+-dependent dioxygenases, which regulate DNA demethylation, modulate gene transcription, and influence various life processes such as embryonic development, as well as the onset of various diseases including tumors (Rasmussen and Helin, 2016; Ma et al., 2021). Additionally, studies have demonstrated that TET family genes not only facilitate active DNA demethylation but also impede methylation propagation by maintaining a low DNA methylation state, thus closely associating them with tumorigenesis (Shekhawat et al., 2021). Notably, TET family genes exhibit a dual role in different tumors, displaying both carcinogenic and anticancer effects.
While TET family genes are known to play an indispensable role in the progression of numerous tumors, there is currently scarce literature on pan-cancer analysis of these genes. Therefore, in this study, we conducted an analysis of the relationship between TET family genes and pan-cancer, exploring various aspects including mRNA expression levels, somatic cell variations, DNA methylation patterns, pathway activities, survival rates and prognoses, immune subtypes specific to individual tumors, tumor microenvironment (TME), stem cell indices, and drug sensitivities. Our aim is to provide valuable reference data and constructive suggestions based on a comprehensive understanding of the landscape.
Results
Expression and differential analysis of TET family genes in pan-cancer
Firstly, utilizing transcriptome data from 33 tumors available in the UCSC Xena database, we conducted an extraction and analysis of the expression profiles of the target genes belonging to the TET family (specifically, TET1, TET2, and TET3) across various types of cancer. The results depicted in the Boxplot revealed that the expression level of the TET3 gene ranked highest across pan-cancer samples, followed by TET2, with TET1 exhibiting the lowest expression (Figure 1A, Supplementary Table S1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418456/fphar-15-1418456-HTML-r1/image_m/fphar-15-1418456-g001.jpg)
Figure 1. Expression and differential analysis of TET family genes in pan-cancer. (A) TET1, TET2, and TET3 expression was assessed across pan-cancer samples. (B–D) Differences in TET1-3 expression between pan-cancer and normal tissues were analyzed. (E) A heatmap visualized expression disparities among tumors using Wilcox analysis. (F) The expression correlation among TET family genes was analyzed in pan-cancer. (G–I) TET3 protein expression was evaluated across various cancer cell types (G), lung tumor subtypes (H), and gastric cancer cell differentiation levels (I). Data are shown as mean ± SD, with statistical significance denoted as *p < 0.05, **p < 0.01, ***p < 0.001.
Next, we screened 18 types of tumors, each containing at least five corresponding normal samples, and analyzed the expression differences of TET family members in pan-cancer and its normal tissues. As shown in Figure 1B, TET1 exhibited abnormal expression in 10 types of tumors, with significantly downregulated expression observed in BRCA, KICH, KIRC, KIRP, and THCA, and upregulated expression noted in CHOL, HNSC, LIHC, LUAD, and LUSC tumors compared to their corresponding normal samples. The expression of TET2, in tumor samples such as BRCA, COAD, HNSC, KIRP, READ, and THCA, was lower, while in CHOL, GBM, KIRC, and UCEC tumors, it was significantly higher than that in corresponding normal control tissues (Figure 1C). The expression of the TET3 gene varied across 15 tumors and corresponding normal samples, including BLCA, BRCA, CHOL, ESCA, HNSC, KIRC, KIRP, LIHC, LUAD, LUSC, PRAD, STAD, THCA, UCEC, and KICH. Among these, except for KICH, TET3 was significantly overexpressed in all 14 other tumors (Figure 1D). These findings suggest that the high expression of TET3 may serve as a key factor in the occurrence, development, and poor prognosis of various clinical tumors.
To further elucidate the differential expression patterns of TET family members in various tumors, we conducted Wilcox analysis on the expression differences in each tumor and obtained the aforementioned differential expression cluster analysis graph (Figure 1E; Table 1). As previously mentioned, TET3 was significantly upregulated, while TET1 and TET2 were downregulated in most tumors. It remains unclear whether there is a correlation between the expression of the three members of the TET family in pan-cancer. Therefore, we conducted a detailed analysis of the correlation between TET family genes in pan-cancer. As illustrated in Figure 1F, the expression of TET1, TET2, and TET3 positively correlated in pan-cancer, with a stronger positive correlation observed between TET2 and TET3.
Given the high expression of TET3 in numerous tumors, we speculate that it may represent a potential therapeutic target. Consequently, we assessed the expression of TET3 protein in four tumor cell lines (BGC-823, HepG2, A549, TPC-1) and their corresponding normal cell lines (GES-1, LX2, 16HBE, Nthy-ori-3-1). The results demonstrated higher expression levels of TET3 in all four tumor cell lines compared to the corresponding control cell lines (Figures 1G–I). Subsequently, we examined the expression of TET3 protein in three different pathological subtypes of lung cancer cell lines (A549, SKMES1, NCIH460). The findings revealed significantly higher expression levels of TET3 protein in lung adenocarcinoma A549 and lung squamous cell carcinoma SKMES1 compared to normal controls and large cell lung cancer cell lines. This observation suggests differential expression of TET3 protein within the same type of tumors with varying pathological subtypes. Finally, we evaluated the expression of TET3 protein in gastric cancer cell lines (NCI-N87, BGC-823, HGC-27) with differing degrees of differentiation. The results indicated that as the degree of cell differentiation decreased, the expression of TET3 protein gradually increased. This implies that the expression intensity of TET3 protein may serve as an indicator of the malignant level of gastric cancer.
Somatic cell variation landscape of TET family genes
To comprehend the genomic alterations of TET1-3 in tumors, we scrutinized single nucleotide variant data (SNV, Table 2) and copy number variant data (CNV, Table 3) from a total of 10,289 pan-cancer samples spanning 33 types of tumors, and illustrated SNV and CNV landscapes.
Initially, we analyzed the SNP data pertaining to the TET gene to ascertain the frequency and variant types across each cancer subtype. As depicted in Figures 2A, B, among these cancers, the SNV frequency of UCEC, SKCM, and COAD ranges from 19% to 57%. The SNV frequency of the TETs gene collectively is 91.24% (1052/1153). Variant analysis unveiled missense mutation as the predominant SNP type. SNV percentage analysis revealed mutation rates of 23%, 17%, and 16% for TET1-3, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418456/fphar-15-1418456-HTML-r1/image_m/fphar-15-1418456-g002.jpg)
Figure 2. Somatic mutations in TET family genes were analyzed. (A) Single Nucleotide Variants (SNV) frequencies of TETs across cancers are presented in a table, indicating mutated gene counts per cancer type. ‘0’ indicates that there was no mutation in the gene coding region, and no number indicates there was no mutation in any region of the gene. (B) SNV variant types of TETs illustrate mutation distribution and categorization. (C) copy number variant (CNV) distribution across 33 cancers is depicted in pie charts, showing combined heterozygous/homozygous CNV proportions for each gene in each cancer. Hete Amp = heterozygous amplification; Hete Del = heterozygous deletion; Homo Amp = homozygous amplification; Homo Del = homozygous deletion; None = no CNV. (D, E) Heterozygous and homozygous CNV profiles display the percentage of CNV, including amplification and deletion, for each gene in each cancer. Only genes with >5% CNV in a given cancer are shown as a point on figure. (F) CNV correlation with mRNA expression is depicted through Person’s correlation analysis, with point size indicating statistical significance. FDR, false discovery rate.
To elucidate the mutation pattern of CNV, we examined the CNV data of TETs genes in the TCGA database. The distribution of CNV pie charts illustrated heterozygous amplification and heterozygous deletion as the primary types (Figure 2C). CNV percentage analysis unveiled heterozygous amplification rates of TET1 in ACC and UCS, TET2 in ACC and KICH, and TET3 in ESCA, LUSC, OV, TGCT, and UCS, all exceeding 20% (Figure 2D). However, homozygous amplification was weak, observed solely in CHOL for TET1 (Figure 2E). Heterozygous deletions of TET1 in GBM, KICH, SARC, SKCM, and TGCT, TET2 in CHOL, ESCA, LIHC, LUSC, MESO, OV, READ, TGCT, and UCS, and TET3 in KICH and SARC, all surpassed 40% (Figure 2D). No homozygous deletions were noted across all tumors (Figure 2E).
Correlation analysis revealed a positive association between the mRNA expression of TETs and CNV, notably with TET2 in CHOL and ACC. Conversely, TET2 in LAML and TET3 in THMY exhibited negative correlations with CNV (Figure 2F; Table 4). This suggests that CNV of the TETs genes mediate abnormal expression in pan-cancer, potentially playing a pivotal role in cancer progression.
Methylation variation in TET family genes across pan-cancer
We delved into the methylation levels of TET genes across various cancers to decipher their epigenetic regulation. Our findings revealed heterogeneity in the methylation levels of TET family genes among different tumors. Specifically, TET1 gene exhibited hypermethylation in most cancers, with hypomethylation observed solely in LIHC. Conversely, TET2 and TET3 genes displayed hypomethylation in the majority of tumors, except for instances of hypermethylation in THCA, COAD, and KIRP for TET2 (Figure 3A). Correlation analysis between methylation and mRNA expression unveiled a negative correlation between the expression of these genes and their methylation levels across pan-cancer (Figure 3B; Table 2). Furthermore, survival analysis indicated that methylation of TET1 was linked to a lower survival rate in ACC and a higher survival rate in UVM; methylation of TET2 was associated with increased survival rates in ACC, KIRP, and PCPG; while methylation of TET3 correlated with decreased survival rates in SARC and BLCA (Figure 3C, Supplementary Table S2). These observations suggest that the methylation levels of TET genes are primarily associated with the survival outcomes of a limited number of tumors.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418456/fphar-15-1418456-HTML-r1/image_m/fphar-15-1418456-g003.jpg)
Figure 3. Methylation analysis of TET family genes in pan-cancer. (A) Methylation levels of TETs genes were examined to assess epigenetic regulation. (B) Correlation between methylation and mRNA expression of TETs genes in pan-cancer was analyzed. Then, association of survival rate with methylation level of TETs genes was investigated (C).
Pathway activity analysis of TET family genes in pan-cancer
Pathway activity analysis revealed significant involvement of TETs in cancer-related signaling pathways, including Cell Cycle, Apoptosis, DNA Damage Response, Hormone AR, PI3K/AKT, and RTK. Notably, TET1 primarily activated the DNA Damage Response pathway (19% activation vs. 0% inhibition), while TET2 mainly participated in the Cell Cycle (4% activation vs. 28% inhibition) and DNA Damage Response pathways (10% activation vs. 9% inhibition). TET3 showed activation primarily in the PI3K/AKT (19%) and the RTK pathway (16% activation vs. 6% inhibition) (Figure 4A; Table 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418456/fphar-15-1418456-HTML-r1/image_m/fphar-15-1418456-g004.jpg)
Figure 4. Pathway activity analysis of TET family genes in pan-cancer. (A) Combined percentage of TET genes’ effect on pathway activity is depicted. Red indicates activation, while blue indicates inhibition. Numbers in the table represent the percentage of pathway activity on TETs. Red represents activate effect and blue represents inhibition. (B) Impact of TET3 on protein expression levels of key signaling pathway members was assessed via Western blot. (C, D) Therapeutic sensitivity of lung and gastric cancer cells to Sunitinib, a multi-target tyrosine kinase inhibitor, was tested to elucidate the role of TET3 in regulating the RTK pathway as shown in Figure 4A.
Given the significance of TET3, we modulated its expression in the NCI-N87 gastric cancer cell line. Overexpression of TET3 led to increased phosphorylation levels of AKT (Ser 308 and Ser 473) and mTOR (Ser2448), suppressed P21 expression, and enhanced HIF1α and c-Myc protein expression, whereas silencing TET3 exhibited the opposite effect (Figure 4B).
Sunitinib, a selective multi-target tyrosine kinase inhibitor, was employed to assess whether TET3 affects tumor behavior via the RTK pathway. We tested the therapeutic sensitivity of three different pathological subtypes of lung cancer cell lines (A549, SKMES1, NCIH460) to sunitinib. The results showed that A549 and SKMES1 cell lines with relatively high expression of TET3 exhibited higher sensitivity to sunitinib (Figure 4C). Similarly, we also tested the therapeutic sensitivity of gastric cancer cell lines (NCI-N87, BGC-823, HGC-27) with different degrees of differentiation to sunitinib. The results revealed that those with higher TET3 expression exhibited increased sensitivity to the drug (Figure 4D). These findings underscore the involvement of TET3 in regulating the PI3K/AKT and RTK pathways during tumor progression.
Influence of TET gene family on pan-cancer survival and prognosis
We conducted an analysis to assess the correlation between TET gene expression levels and the survival outcomes of patients across 33 types of tumors using data from the TCGA database. Kaplan-Meier survival curve analysis revealed that high expression of TET1 in ACC, KIRP, LIHC, SARC, and STAD was associated with shorter patient survival, while low expression in LGG and THYM indicated reduced survival. Similarly, high expression of TET2 in OV, PCPG, and UCS was linked to shorter survival, whereas low expression in KIRC correlated with decreased survival. Additionally, high expression of TET3 in ACC, MESO, and PCPG was associated with shorter survival, while low expression in THCA indicated reduced survival (Figure 5, Supplementary Table S3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418456/fphar-15-1418456-HTML-r1/image_m/fphar-15-1418456-g005.jpg)
Figure 5. Analysis of TET family genes expression and patient prognosis using univariate KM risk proportional regression model.
Furthermore, we evaluated the relationship between TET expression and cancer prognosis (Figure 6). COX regression analysis across 33 cancer types revealed that high expression of TET1 in ACC, BLCA, CESC, LIHC, PCPG, and SARC was indicative of poor prognosis, while in LGG, THYM, and UVM, it suggested a favorable prognosis. Conversely, in KIRC, high expression of TET2 was associated with a favorable prognosis. Notably, high expression of TET3 in ACC alone indicated a poor prognosis among tumor patients.
Correlation between TET family genes expression and immune subtypes in pan-cancer and single tumors
To elucidate the association between TETs and immune responses, we conducted a pan-cancer analysis correlating TETs expression with immune subtypes using the TCGA database. As depicted in Figure 7A, TETs expression exhibited a significant positive correlation with pan-cancer immune subtypes (p < 0.05): notably, TET1 and TET2 were markedly upregulated in immune subtype C5, whereas their expression remained stable in other subtypes; TET3 displayed variable expression levels across immune subtypes. To further delineate the specific immune subtype associations with TETs, we performed KS test correlation analysis in eight common tumors. Results indicated significant correlations between TETs expression and immune subtypes in most tumors (p < 0.05), with TET1 consistently exhibiting the lowest expression across all subtypes, while TET3 demonstrated high expression levels (Figure 7B). In STAD, LIHC, and LUAD, TETs expression showed significant positive correlations with various immune responses. Notably, TET1 displayed highest expression in STAD’s C3 subtype and minimal expression in C2, while TET2 was consistently upregulated across all subtypes. Conversely, in LIHC, TET1 and TET3 exhibited decreasing expression trends across subtypes, with TET1 almost absent in C6, and TET2 lowest in C2. In LUAD, TET1 expression remained stable, while TET2 peaked in C3 and declined in C4, and TET3 expression was highest in C1 and lowest in C4. However, COAD and LUSC immune subtypes showed significant correlations only with TET1 and TET3 expression. TET1 was highly expressed across all COAD immune subtypes, while TET3 was least expressed in C3. In LUSC, TET1 displayed minimal expression in C3 and C4, while TET3 was lowest in C6 and stable in other subtypes. In summary, TET1 exhibited consistently low expression across various tumors, while TET3 showed highest expression levels, both closely associated with pan-cancer immune subtypes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418456/fphar-15-1418456-HTML-r1/image_m/fphar-15-1418456-g007.jpg)
Figure 7. Correlation between expression of TET family members and immune subtypes. (A) Analysis of the relationship between TET family genes expression and immune subtypes across pan-cancer, based on the TCGA database. (B) Further analysis of the correlation between immune subtypes of eight common clinical tumors and TETs expression levels. Data are presented as the mean ± SD. *p < 0.05, **p < 0.01, ***p < 0.001.
Correlation analysis of TET family genes expression with tumor microenvironment (TME) and stem cell index across pan-cancer
The efficacy of radiotherapy, chemotherapy, and immunotherapy in treating tumors is greatly influenced by the Tumor Microenvironment (TME). Utilizing the ESTIMATE algorithm, we assessed TME characteristics across 33 tumor types by computing Stromal Score, Immune Score, and ESTIMATE Score. Analysis from Figures 8A–C reveals that TET1 exhibits higher Stromal Scores in COAD, ESCA, KIRP, MESO, PAAD, and STAD, with elevated Immune and ESTIMATE Scores observed specifically in PAAD. These findings suggest that upregulating TET1 expression may enhance TME in PAAD, potentially suppressing tumor invasion, metastasis, and augmenting immunotherapy sensitivity. Conversely, heightened TET1 expression may exacerbate TME conditions in GBM, LGG, and TGCT, potentially escalating invasion, metastasis, and drug resistance.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418456/fphar-15-1418456-HTML-r1/image_m/fphar-15-1418456-g008.jpg)
Figure 8. TET family members and tumor microenvironment in pan-cancer. (A–C) Association between TET family genes expression and tumor microenvironment parameters including tumor matrix score, immune score, and tumor purity score in pan-cancer. (D, E) Relationship between TET family genes expression and RNAs and DNAs in pan-cancer.
Similarly, TET2 displays elevated Stromal Scores in CHOL, LAML, and MESO, alongside increased Immune Scores in DLBC and LAML, and higher ESTIMATE Scores in DLBC and LAML. Conversely, reduced Stromal Scores are noted in BLCA and GBM, while diminished Immune Scores are observed in ACC, BLCA, GBM, PCPG, SARC, and THYM, collectively translating to overall lower scores in ACC, BLCA, GBM, and PCPG. This suggests that restoring TET2 expression may ameliorate TME in LAML, potentially suppressing invasion, metastasis, and bolstering immunotherapy sensitivity. However, it may exacerbate TME conditions in BLCA and GBM, fostering invasion, metastasis, and drug resistance.
Meanwhile, TET3 exhibits heightened Stromal, Immune, and ESTIMATE Scores in KICH, KIRC, and LAML, while displaying diminished Stromal Scores in BLCA, ESCA, GBM, LUSC, and STAD, alongside reduced Immune Scores in BLCA, CESC, ESCA, GBM, LUSC, and UCEC, culminating in overall lower scores in BLCA, CESC, ESCA, GBM, LUSC, and UCEC. This suggests that inducing TET3 expression may enhance TME in KICH, KIRC, and LAML, potentially suppressing invasion, metastasis, and augmenting immunotherapy sensitivity. Conversely, it may deteriorate TME conditions in BLCA, CESC, ESCA, GBM, LUSC, and UCEC, fostering invasion, metastasis, and drug resistance.
Moreover, the therapeutic efficacy of radiotherapy, chemotherapy, and immunotherapy is closely linked to the tumor stem cell index. Figures 8D, E indicates that TET1 expression is positively correlated with RNAss in LGG and TGCT, while negatively correlated with KIRC, PAAD, and STAD. TET2 levels are negatively correlated with RNAss in CHOL, KIRP, and THYM. Conversely, TET3 levels are negatively correlated in KICH, KIRP, THCA, and THYM, and positively correlated in BLCA, BRCA, and PRAD. Additionally, TET1 expression is positively correlated with DNAss in GBM and TGCT, but negatively correlated with LGG, OV, and THYM. TET2 shows a positive correlation in TGCT and a negative correlation in OV, while TET3 exhibits a negative correlation in OV.
Next, we conducted a detailed analysis of the tumor microenvironment and stem cell index of eight common clinical cancers (Figure 9). The results showed that for tumor microenvironment-related scores, TET1 was significantly positively correlated with the three scores of PAAD, and negatively in LUAC. Its stem cell index is negatively correlated in STAD and COAD, and positively in LUAD and LUAC. TET2 shows a significant positive correlation with the TME score in COAD, and a negative correlation with the stem cell index in PAAD. TET3 is significantly negatively correlated with TME scores in STAD and LUAC, and positively with the stem cell index in STAD, LUAD, and LUAC. In addition, several tumors had no significant correlation with their microenvironment score and stem cell index.
TET family genes expression was correlated with clinical stages across various cancers, indicating their potential as prognostic markers
These findings underscore the close relationship between TET gene expression and both the tumor microenvironment and tumor stemness potential. Further analysis revealed distinct correlations between TET gene expression and clinical stages across eight common tumor types. As shown in Figure 10, there was no significant difference between the expression level of TETs genes and clinical stage of tumors in COAD and READ. Notably, TET2 and TET3 exhibited significant correlations with clinical stages in LIHC, with elevated expression in Stage III and reduced expression in Stage IV. Similarly, in LUSC, TET1 and TET3 expression levels were notably associated with clinical stages, showing upregulation in Stage III and downregulation in Stage IV. Furthermore, TET2 expression levels were significantly linked to clinical staging in STAD and LUAD, with varying expressions across different stages. In PAAD, TET3 expression levels exhibited significant correlation with clinical characteristics, displaying the lowest expression in Stage III and the highest expression in Stage IV. These results maybe highlight the distinctive correlations between TET gene expression and clinical stages across various tumors.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418456/fphar-15-1418456-HTML-r1/image_m/fphar-15-1418456-g010.jpg)
Figure 10. Relationship between TET gene family expression and clinical characteristics in eight common tumors.
TET family genes levels were analyzed for their impact on drug sensitivity
Using the CellMiner ™ database, we conducted an analysis to explore potential correlations between gene expression levels of TET family members and drug sensitivity across various human tumors (Figure 11A, Supplementary Table S4). We identified the top 16 chemotherapy-sensitive drugs exhibiting the highest correlation coefficients with TETs genes. Notably, the expression level of TET1 displayed significant positive correlations with drug sensitivity to Arsenic trioxide, Fenretinide, Dimethylaminophenhenolide, Daunorubicin, Homoharringtonine, Imatinib, Testolactone, Pipappererone, and Lomustine. Conversely, TET2 expression was positively correlated with drug sensitivity to Fulvestrant and Raloxifene but significantly negatively correlated with Vemurafenib and Dabrafenib. Moreover, the expression level of TET3 exhibited significant positive correlations with drug sensitivity to Lapatinib, AZD-9291, and (+)-JQ1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418456/fphar-15-1418456-HTML-r1/image_m/fphar-15-1418456-g011.jpg)
Figure 11. Correlation analysis between TET gene family expression and drug sensitivity. (A) The potential correlation between gene expression levels and drug sensitivity of TET family members were analyzed based on the CellMiner™ database. (B) We tested the therapeutic sensitivity of three different pathological subtypes of lung cancer cell lines to Lapatinib, AZD-9291, and (+) - JQ1.
Furthermore, we assessed the therapeutic sensitivity of three distinct pathological subtypes of lung cancer cell lines (A549, SKMES1, NCIH460) to Lapatinib, AZD-9291, and (+)-JQ1. The findings indicated that the therapeutic efficacy of these drugs on tumor cells positively correlated with the expression level of TET3 (Figure 11B).
Silencing TET3 effectively suppresses the malignant behavior of tumor cells
Following our comprehensive analysis spanning various cancer types, TET3 has emerged as a prime candidate for targeted therapeutic intervention. In subsequent experiments, where we selectively silenced TET3 expression in four distinct tumor cell lines (BGC-823, HepG2, A549, TPC-1) (Figure 12A), we observed a profound attenuation of multiple malignant phenotypes. Specifically, the inhibition of TET3 led to a notable reduction in clone formation capacity (Figure 12B), as well as a significant impairment in cell migration and invasion abilities (Figures 12C, D). Furthermore, TET3 silencing induced cell cycle arrest (Figure 12E), resulting in decreased proliferation rates, while simultaneously promoting apoptosis (Figure 12F). These compelling findings underscore the pivotal oncogenic role of TET3 across various tumor contexts and highlight the potential of TET3-targeted therapies as a promising avenue for effective cancer treatment.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418456/fphar-15-1418456-HTML-r1/image_m/fphar-15-1418456-g012.jpg)
Figure 12. (Continued). Silencing TET3 can inhibit the malignant behavior of tumor cells. (A) The expression of TET3 in four tumor cell lines (BGC-823, HepG2, A549, TPC-1) was silenced. (B–F) Silencing TET3 significantly inhibited tumor cell clone formation ability (B), migration ability (C), and invasion ability (D); and also inhibited the cell proliferation cycle (E) while increased cell apoptosis (F).
Discussion
DNA methylation patterns impact gene expression and are often disrupted in diseases like inflammation, hypertension, diabetes, and tumors (Bowman and Levine, 2017; Ismail et al., 2020; Bray et al., 2021; Matuleviciute et al., 2021). TET enzymes catalyze the oxidation of 5-methylcytosine to stable epigenetic modifications, playing crucial roles in gene regulation (Delatte et al., 2014; Scott-Browne et al., 2017; Bray et al., 2021; Matuleviciute et al., 2021; Joshi et al., 2022). These oxidized 5-methylcytosine derivatives serve as stable epigenetic modification that exert distinctive regulatory roles (Kao et al., 2016). Dysregulated TET protein expression is common in various cancers (Pan et al., 2015; Smeets et al., 2018; Kunimoto and Nakajima, 2021).
In our analysis of transcriptome data from 33 tumors in the UCSC Xena database, we found distinct expression patterns among TET genes across different cancer types. Specifically, TET1 showed notable overexpression in epithelioid squamous cell carcinoma (HNSC and LUSC) and hepatobiliary duct tissue (CHOL and LIHC), while consistently exhibiting low expression in various kidney-derived tumors (KICH, KIRC, and KIRP). Conversely, TET2 was predominantly expressed in adenocarcinomas (BRCA, COAD, READ, and THCA). In contrast, TET3 demonstrated widespread high expression across most tumors, with the exception of KICH. These expression patterns may hold significance for clinical diagnostics and the development of targeted therapies.
Certain SNPs within genes can directly impact protein structure or expression levels, potentially influencing tumor genetic mechanisms (Kuhlen et al., 2019). In this study, we conducted a systematic analysis of somatic cell variations in TET gene family members. We observed a high frequency of Single Nucleotide Variants (SNV) and found that the main type of copy number variant (CNV) was heterozygous amplification and deletion. Correlation analysis revealed positive associations between TET2 expression and CNV in CHOL and ACC, while TET2 in LAML and TET3 in THMY showed negative correlations with CNV. Identifying pathogenic CNV and interpreting their clinical significance will be crucial for future research, despite the challenges associated with this endeavor.
In general, the extent of whole-genome hypomethylation in tumor cells correlates closely with disease progression, tumor size, and malignancy (Esteller and Herman, 2002). DNA methylation holds significant value in assessing tumor malignancy and guiding targeted drug selection. The TET family genes, pivotal in modulating the methylation levels of numerous genes, contribute to the complexity of gene expression regulation due to their own methylation status. Although these mechanisms enhance gene regulation accuracy, they pose challenges for researchers. Despite the negative correlation between TET family genes expression and methylation levels in pan-cancer, unique patterns were observed, such as high methylation of TET1 in most analyzed cancers (except LIHC) and low methylation of TET2 and TET3 in most tumors. Further investigation is required to elucidate whether this differential expression pattern suggests functional compensation among TET family genes. Notably, the state of DNA methylation evolves during tumor progression, necessitating dynamic and cautious interpretation of its clinical significance (Miyamoto and Ushijima, 2003).
Our analysis of TET family genes in pan-cancer revealed their involvement in various signaling pathways. Notably, TET1 was found to predominantly activate the PI3K/AKT pathway. Prior studies have highlighted TET1’s role in driving the proliferation of insulin-dependent endometrial cancer by enhancing G protein-coupled estrogen receptor expression and PI3K/AKT pathway activation (Xie et al., 2017). Moreover, research by Huang elucidated TET1’s impact on stem cell development through modulation of the Wnt and PI3K-Akt pathways, while TET2 deficiency in mice led to a progressive reduction in spermatogonia stem cells (Huang et al., 2020). TET2 plays a pivotal role in cell cycle regulation and DNA damage responses. Studies by Zhong have revealed that 5 mC oxidation is cell-cycle dependent, occurring primarily during the S and G2/M phases. Notably, TET2 depletion diminishes the observed elevation in 5hmC, indicating its dependence on TET2, particularly in response to idarubicin stimulation, a topoisomerase II inhibitor (Panigrahi et al., 2015). Additionally, Chen found that SMAD nuclear interacting protein one recruits TET2 to regulate c-MYC target genes and the cellular DNA damage response (Chen et al., 2018). TET3, a pivotal enzyme, showcases its versatility within cellular processes, including the cell cycle, apoptosis, hormone AR regulation, and the DNA damage response. Research led by Jiang and colleagues revealed that upon DNA damage, ATR kinase activation leads to TET3 phosphorylation in mammalian cells. This phosphorylation enhances DNA demethylation and the accumulation of 5-hydroxymethylcytosine, underlining TET3’s essential role in DNA repair and cell survival (Jiang et al., 2017).
In the landscape of cancer, the expression levels of the TET family offer prognostic insights. Elevated TET1 expression correlates with worse outcomes in solid organ tumors such as ACC, KIRP, LIHC, and SARC, yet it signifies better prognoses in LGG and THYM cancers. Similarly, increased TET2 expression suggests poorer outcomes in tumors of the female reproductive system, including OV and UCS, but predicts favorable outcomes in KIRC. Notably, high TET3 expression is linked to unfavorable prognoses in ACC but indicates better survival rates in THCA. These observations underscore the multifaceted roles that TET enzymes play in the vast expanse of cancer prognosis.
Diving deeper into the complex interplay between TET enzyme expression and immune subtypes, this pan-cancer study unveils their pivotal influence within the tumor immune microenvironment. The findings suggest that TET1 and TET2, especially upregulated in the C5 immune subtype known for its inflammatory profile, might influence immune evasion or activation. TET3, with its variable expression, appears to affect a broad range of immune responses across tumors. The distinct expression patterns of TET enzymes across different cancers and immune subtypes underscore their potential in modulating tumor immunity, progression, and therapy response. This research not only sheds light on the complex interactions between TET enzymes and the immune system but also emphasizes the importance of TET genes as biomarkers and therapeutic targets in cancer. It advances our understanding of the tumor immune microenvironment, paving the way for novel immunotherapeutic strategies.
Building on this understanding, the study further explores the relationship between TET expression and the TME, a critical factor in tumor survival, immune evasion, and drug resistance (Arneth, 2019). For instance, augmenting TET1 expression may improve the TME in PAAD while potentially exacerbating it in GBM, LGG, and TGCT. Insights from Li suggest that TET1 could inhibit epithelial-mesenchymal transition and increase PAAD cells’ sensitivity to chemotherapy agents like 5FU and gemcitabine (Li et al., 2020). However, the impact of TET1 overexpression on GBM, LGG, and TGCT TME remains unexplored. Secondly, TET2 expression restoration may improve LAML’s TME but worsen that of BLCA and GBM. Cimmino et al. found that TET2 restoration might reverse aberrant hematopoietic stem and progenitor cell self-renewal in vitro and in vivo (Cimmino et al., 2017). Similarly, the effect of inducing TET2 expression on BLCA and GBM TME remains unclear. Lastly, inducing TET3 expression could enhance KICH, KIRC, and LAML TME while deteriorating BLCA, CESC, ESCA, GBM, LUSC, and UCEC TME. This comprehensive analysis sheds further light on the nuanced roles of TET enzymes in the TME, highlighting their potential as biomarkers and therapeutic targets in the ongoing battle against cancer.
Further investigation into the tumor stem cell index revealed a positive correlation between heightened TET1 expression and increased levels of ribonucleic acid synthesis (RNAss) and deoxyribonucleic acid synthesis (DNAss) in testicular germ cell tumors (TGCT). Notably, Benešová et al. observed a marked increase in TET1 dioxygenase expression in most seminomas, suggesting its potential utility as a marker for seminomas and mixed germ cell tumors (Benešová et al., 2017). These findings underscore the potential of TET1 as a pivotal indicator for evaluating stemness maintenance and chemoradiotherapy resistance in TGCT tumor stem cells.
The examination of TET family genes expression concerning the clinical staging of common tumors reveals nuanced patterns. While TET1 demonstrates only limited correlation with tumor staging, TET2 exhibits significant associations in STAD, LIHC, and LUAD (p < 0.01). Similarly, TET3 shows correlations with staging in PAAD, LIHC, and LUSC (p < 0.05). These findings are in line with the research by Liu, indicating decreased genomic 5hmC and 5 fC contents in early LIHC stages, with further reductions in late stages. Moreover, the significantly positive correlations among the expression levels of TET2 in para-tumor tissues were generally attenuated or even disappeared in LIHC tumor tissues (Liu et al., 2019). Moreover, Sajadian underscore the impaired expression and activity of TET2 and TET3 in hepatocellular carcinoma, further validating our analysis (Sajadian et al., 2015). These conclusions are close to our analysis results.
Effective cancer treatment relies on understanding drug sensitivity, a cornerstone of personalized therapy and precision medicine advancement (Chaudhry and Asselin, 2009). Yet, due to individual variations and disparate drug responses, optimal resource utilization remains a challenge (Panczyk, 2014). Therefore, investigating molecules influencing drug reactions is essential for refining treatment strategies. Our study delves into the interplay between TET family genes expression and drug sensitivity, yielding significant insights. Notably, we find that TET genes expression levels correlate with the efficacy of specific drugs, highlighting the importance of assessing TET1, TET2, and TET3 expression for informed clinical drug selection.
Our study delved into the intricate role of TET family genes across various aspects of cancer biology. We identified TET1, TET2, and TET3 as pivotal players influencing tumor progression, prognosis, immune response, tumor microenvironment, and drug sensitivity. Notably, the analysis reveals that each TET family member, including TET3, displays unique expression patterns in at least ten detected tumor types. This heterogeneity suggests that TET3 may play distinct roles in different cancer subtypes. Furthermore, the finding that TET3 genes display hypomethylation in most cancers, which correlates closely with patient prognosis, highlights its potential involvement in cancer progression and metastasis. The association between TET3 expression and various cancer-related factors, such as pathway activity, tumor microenvironment, stemness score, immune subtype, clinical staging, and drug sensitivity, further underscores its relevance as a potential therapeutic target. The results from molecular biology and cytology experiments validating the potential role of TET3 in tumor progression strengthen this argument.
In summary, studying TET3 in different cancers is highly relevant due to its potential as a therapeutic target. The comprehensive pan-cancer analysis presented in the manuscript provides a foundation for future research aimed at developing targeted therapies that may improve cancer treatment outcomes.
Conclusion
The comprehensive analysis of TET family genes in pan-cancer unveiled their multifaceted roles. Through transcriptome data analysis, distinct expression patterns were observed across various tumor types, indicating their potential diagnostic significance. Moreover, correlation with prognosis highlighted their prognostic value, while associations with the tumor microenvironment and drug sensitivity underscored their therapeutic implications. These findings suggest that TET genes could serve as valuable targets for precision medicine approaches in cancer treatment.
Materials and methods
Data download preparation
Based on data obtained from the UCSC Xena database, RNA-Seq and clinical data for 33 tumor types prefixed with “GDC TCGA” were downloaded. This includes “HTSeq FPKM (n = 151) GDC Hub” data for gene expression RNAseq, as well as “Phenotype (n = 697) GDC Hub (Clinical Traits)" and “Survival Data (n = 626) GDC Hub” under the Phenotype category. TCGA pan-cancer (PANCAN) data, such as “Immune subtype” under Phenotype and “Stemness score (DNA methylation based) pan-cancer Atlas Hub” and “Stemness score (RNA based) pan-cancer Atlas Hub” under Signatures, was also retrieved. Additionally, drug sensitivity information was obtained from the CellMiner™ database (https://discover.nci.mih.gov/cellminer/home.do).
Differential analysis of gene expression
Utilize Perl software to extract, transform, and integrate the ‘HTSeq FPKM (n = 151) GDC Hub’ data obtained from the Gene expression RNAseq item. Generate boxplots to illustrate the expression profiles of TET family genes across diverse tumor samples. Next, filter samples with a minimum of five normal controls per tumor type and create gene expression boxplots accordingly. Conduct expression difference analysis of TET family genes in different cancer types using the ‘Wilcox. test’ method, with statistical significance levels indicated by '*', '**', and '***' for p < 0.05, <0.01, and <0.001, respectively. Utilize the R package ‘Pheatmap’ to generate a heatmap based on the resulting p-values. Finally, examine the correlation between genes within the TET family using the R package ‘Coreplot'.
Somatic mutation analysis
Data on single nucleotide variant (SNV) and copy number variant (CNV) for 33 tumor types were retrieved from the TCGA database through the Xena Functional Genomics Explorer (https://xenabrowser.net/datapages/). The SNV data encompasses various non-silent mutations, including Missense_Mutation, Nonsense_Mutation, Frame_Shift_Del, Splice_Site, Frame_Shift_Ins, In_Frame_Del, In_Frame_Ins, Translation_Start_Site, and Multi_Hit. The SNV mutation frequency (%) for each gene coding region is calculated as the number of mutation samples divided by the total number of cancer samples. Finally, SNV landscape maps were generated using Maftools.
For CNV analysis, we classified CNV into two types: homozygous and heterozygous, representing amplifications and deletions, respectively. We then computed the percentage statistics for each CNV subtype using GISTIC processed CNV data. Next, we investigated the correlation between CNV and mRNA expression levels. Genes with CNV exceeding 5% were identified, and their association with TET expression was explored. Utilizing the method described by Tyagi (Tyagi et al., 2024), we assessed the correlation between mRNA expression and CNV percentage samples using Pearson’s product-moment correlation coefficient, with p-values corrected for false discovery rate (FDR).
Methylation analysis
To perform methylation analysis, we initially obtained DNA methylation data from the UCSC Xena database. We focused on 14 tumor types with a minimum of 10 paired samples of tumor and normal tissues. Differential methylation between tumor and normal samples was assessed using Student’s t-test, with p-values adjusted for false discovery rate (FDR). Significance was determined at FDR <0.05. Subsequently, we integrated the methylation data with TET gene expression data. Spearman correlation coefficient was computed to evaluate the correlation between methylation levels and gene expression.
Further analysis involved categorizing the median methylation level of genes into two groups. This categorization was based on the threshold defined by the median methylation level. Cox regression analysis was then conducted to estimate the hazard ratio (HR) of gene methylation, considering covariates specific to each cancer type. A Cox coefficient >0 indicated worse survival rates for the high methylation group, hence defining it as high-risk, while a Cox coefficient ≤0 indicated low-risk. Additionally, we compared the distribution of the two methylation groups using the log-rank test to assess their association with clinical outcomes. A significance level of p < 0.05 was considered statistically significant in these comparisons.
Signal pathway activity analysis
The reverse protein array (RPPA) data from the UCSC Xena database within TCGA is leveraged to assess pathway activity across various tumor samples. Nine key pathways associated with cancer progression are scrutinized, including Apoptosis, Cell cycle, DNA Damage Response, Epithelial Mesenchyme Transition (EMT), Hormone androgen receptor (AR), Hormone estrogen receptor (ER), Phosphatidylinositol-4,5-bisphosphate-3-kinase (PI3K)/protein kinase B (AKT), Rasopathies (RAS)/mitogen-activated protein kinase (MAPK), and Receptor Tyrosine Kinase (RTK).
Pathway activity scores are calculated by summing the relative protein levels of positive regulatory components and subtracting those of negative regulatory components within each pathway. Following this, employing methodologies outlined in prior studies, such as those by Akbani and Ye. (Akbani et al., 2014), pathway activity scores (PAS) are derived. A higher PAS in one group compared to another suggests an activating effect of certain genes on the pathway, while a lower PAS indicates an inhibitory effect.
Survival and prognostic analysis
We performed expression survival analysis by integrating mRNA expression data of TETs genes with clinical survival data across 33 cancer types, utilizing the Survival Data (n = 626) GDC Hub from the UCSC Xena database. Employing the Kaplan-Meier method and log-rank test (p < 0.05), tumor samples were stratified into high and low expression groups based on median gene expression levels.
Subsequently, we applied a univariate Cox proportional hazards regression model to investigate the relationship between TET family genes expression and patient prognosis across various cancers. Finally, the results were visualized using forest plots generated with the “survival” and “Forestplot” R packages.
Immunological subtype analysis
Utilizing the “Immune subtype” data from the Phenotype category within the TCGA pan-cancer (PANCAN) dataset available on the UCSC Xena database, we conducted immune subtype analysis on the TETs genes. Employing R packages “Limma”, “Ggplot2”, and “Reshape2”, we applied the Kruskal–Wallis (KS) test method to detect expression differences among various immune subtypes. A significance threshold of p < 0.05 was set for statistical significance.
Analysis of tumor microenvironment and stem cell index
Using the “Stemness score (DNA methylation based) pan-cancer Atlas Hub” and “Stemness score (RNA based) pan-cancer Atlas Hub” data within the Signatures section of TCGA pan-cancer (PANCAN), we employed R packages “Estimate” and “Limma” to predict Stromal Score, Immune Score, and ESTIMATE Score for 33 tumor types, enabling an analysis of tumor purity. Following this, we conducted a correlation analysis between TETs gene expression levels and ESTIMATE Scores across the 33 tumor types using Spearman’s correlation coefficient. Additionally, we performed Spearman correlation tests based on transcriptional data and stemness scores (RNA expression and DNA methylation levels).
Drug sensitivity analysis
Download drug sensitivity data from the CellMiner™ database and preprocess it using R packages such as “Input”, “Limma”, “Ggplot2”, and “Ggpubr” to ensure accurate analysis and visualization. Utilize the Cor. Test function to perform correlation analysis on the data, considering a significance level of p < 0.05 as indicative of significant drug sensitivity.
Evaluate the sensitivity of tumor cells to drugs using the Cell Counting Kit-8 (CCK-8) assay following the manufacturer’s protocol (Yeasen, Shanghai, China). This involves treating cultured tumor cells with various drugs and assessing cell viability using the CCK-8 assay.
Cell culture, plasmids, and transfections
All cell lines utilized in this study, including those derived from gastric mucosa (GES-1), liver (LX2), lung epithelium (16HBE), thyroid (Nthy-ori-3-1), as well as various cancer cell lines like BGC-823 (poorly differentiated gastric cancer), HepG2 (hepatoblastoma), A549 (lung adenocarcinoma), TPC-1 (papillary thyroid cancer), SKMES1 (lung squamous cell carcinoma), NCIH460 (large cell lung cancer), NCI-N87 (moderately differentiated gastric cancer), and HGC-27 (undifferentiated gastric cancer), were procured from the American Type Culture Collection (ATCC; United States). These cell lines were cultured following ATCC’s recommended protocols.
TET3-targeting interfering RNAs (siRNA-TET3) and TET3 overexpression plasmid (pcDNA3.1-3*Flag-TET3) were synthesized by GenePharma Co., Ltd (Shanghai, China) and transfected into cells at a concentration of 20 nM using Lipofectamine RNAiMax (Invitrogen, United States). The siRNA sequence used was: siRNA-TET3: 5′-GGAAAGAGCUCCCGCGGUUTT-3'. The plasmid was constructed using the Fast MultiSite Mutagenesis System Kit (FM201-01, TransGen Biotech, China) according to the manufacturer’s instructions. Transfections of the mentioned plasmids were performed using Lipofectamine 2000 reagent (Thermo Fisher Scientific, United States), and the transfected cells were utilized in experiments 48 h post-transfection.
Western blot
Cell lysates were prepared using RIPA lysis buffer, and protein concentration was determined with a BCA kit (Beyotime, Shanghai, China, P0010) following the manufacturer’s instructions. Protein samples were separated on SDS-PAGE gels, transferred to PVDF membranes, and blocked with 5% skimmed milk powder. Primary antibodies were applied overnight at 4°C, followed by secondary antibodies at room temperature. Immunoreactive bands were visualized using ECL (Thermo Fisher Scientific, 34580) and imaged with a ChemiDoc™ MP Imaging System (BioRad, United States of America).
The antibodies used for Western blotting were as follows: TET3 (Abcam, ab153724, 1:1000), β-Actin (CST, #3700, 1:1000), p-AKT (CST, #4060, 1:2000), AKT (CST, #9272, 1:1000), p-mTOR (CST, #2974, 1:1000), mTOR (CST, #2983, 1:1000), p21 (Invitrogen, MA5-14949, 1:1000), Cyclin D1 (Invitrogen, MA5-16356, 1:200), HIF-1α (CST, #3716, 1:1000), c-Myc (Invitrogen, MA1-980, 1:1000).
Behavioral detection of tumor cells
According to the experimental requirements, the flat plate cloning experiment, wound healing assay, cell invasion assay, flow cytometry cell cycle and apoptosis detection involved in this project were all implemented according to standard methods as follows:
The colony formation assay, cells are first seeded at low density into culture dishes or plates and allowed to grow undisturbed for a period of time, typically several days to weeks, depending on the cell type and experimental requirements. During this time, individual cells proliferate and form colonies derived from a single progenitor cell. Once the colonies have reached a suitable size, the cells are fixed and stained to visualize them. Finally, the number of colonies and their size are quantified using microscopy or image analysis software, providing valuable information about cell proliferation, survival, and clonogenic potential. This assay is commonly used to assess the effects of various treatments or genetic manipulations on cell growth and survival.
The wound healing, cells are first cultured in a monolayer until they reach confluence. Then, a scratch or wound is created in the cell layer using a pipette tip or specialized tool. The cells are then washed to remove debris and allowed to incubate in fresh media. Images of the scratch are taken at regular intervals over a specified period, allowing researchers to monitor and measure cell migration into the scratch area. Finally, data analysis involves quantifying the extent of scratch closure or the rate of migration, typically by measuring the remaining scratch width using image analysis software. This assay provides insights into cell migration dynamics and can be used to assess the effects of various factors on cell motility.
The Transwell assay, cells are seeded into the upper chamber of a Transwell insert, while the lower chamber is filled with medium containing chemoattractant. The cells are allowed to migrate or invade through the porous membrane of the insert towards the lower chamber for a specified period of time. After incubation, non-migratory or non-invasive cells on the upper surface of the membrane are removed, while cells that have migrated or invaded to the lower surface are fixed, stained, and counted under a microscope. The number of migrated or invaded cells is quantified to assess the migratory or invasive capacity of the cells in response to different experimental conditions or treatments.
Flow cytometry analysis of the cell cycle and apoptosis, cells are typically harvested, fixed, and permeabilized to preserve their structural integrity. For cell cycle analysis, the fixed cells are treated with DNA intercalating dyes, such as propidium iodide (PI), to stain DNA. The stained cells are then subjected to flow cytometry analysis to measure the DNA content, allowing for the identification of cells in different phases of the cell cycle. Conversely, for apoptosis analysis, cells are stained with fluorescent dyes that selectively bind to apoptotic cells, such as Annexin V and PI. The stained cells are then analyzed by flow cytometry to quantify the percentage of apoptotic cells based on their fluorescence properties. By comparing treated samples to untreated controls, the effects of various treatments or experimental conditions on the cell cycle progression and apoptosis induction can be assessed.
Statistics methods
Each validation experiment included three replicates and was repeated thrice for reliability. Data analysis and graphing were performed using GraphPad software v.5.01, with results displayed as mean ± SEM. Student’s t-test compared two independent groups, while One-way ANOVA assessed multiple groups. Statistical significance was set at p < 0.05, denoted as *p < 0.05, **p < 0.01, ***p < 0.001, and ns for no significant difference.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
CZ: Methodology, Writing–original draft. JZ: Methodology, Writing–original draft. JL: Software, Validation, Writing–original draft. YL: Methodology, Validation, Writing–original draft. GX: Software, Writing–review and editing. SZ: Validation, Writing–review and editing. HC: Formal Analysis, Writing–original draft. JW: Formal Analysis, Writing–review and editing. ZS: Data curation, Writing–review and editing. YL: Methodology, Software, Writing–original draft. ZJ: Resources, Writing–original draft. YP: Data curation, Writing–review and editing. XL: Writing–review and editing. PX: Project administration, Writing–review and editing. WW: Project administration, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by grant from National Natural Science Foundation of China (Grant No. 82272880), Beijing Natural Science Foundation (Grant No. 7212112), Tianjin Municipal Science and Technology Plan Project (Grant No. 22ZYQYSY00030), Tianjin Municipal Natural Science Foundation (Grant No. 21JCZDJC01270), Tianjin Health Technology Project (TJWJ2022XK041, TJWJ2022XK042,TJWJ2022XK043, TJWJ2022MS051, TJWJ2022QN101, TJWJ2023MS053), Tianjin Binhai New Area Health Commission Science and Technology Project (2022BWKQ003), Bethune Charity Foundation Project (B-0307-H20200302), Tianjin Municipal Health Commission research project on the integration of traditional and Western medicine (2021061, 2023062), and also funded by Tianjin Key Medical Discipline (Specialty) Construction Project (TJYXZDXK-062B, TJYXZDXK-079D).
Acknowledgments
We thank the databases of TCGA, Kaplan-Meier Plotter, PrognoScan, and CellMiner™ for the availability of the data. We also would like to express my sincere gratitude to ChatGPT 3.5 for its invaluable assistance in correcting the language aspects of this article.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2024.1418456/full#supplementary-material
References
Akbani, R., Ng, P. K. S., Werner, H. M. J., Shahmoradgoli, M., Zhang, F., Ju, Z., et al. (2014). A pan-cancer proteomic perspective on the Cancer Genome Atlas. Nat. Commun. 5, 3887. doi:10.1038/ncomms4887
Alqahtani, A., Khan, Z., Alloghbi, A., Said Ahmed, T. S., Ashraf, M., and Hammouda, D. M. (2019). Hepatocellular carcinoma: molecular mechanisms and targeted therapies. Med. Kaunas. 55, 526. doi:10.3390/medicina55090526
Ban, J., Fock, V., Aryee, D. N. T., and Kovar, H. (2021). Mechanisms, diagnosis and treatment of bone metastases. Cells 10, 2944. doi:10.3390/cells10112944
Benešová, M., Trejbalová, K., Kučerová, D., Vernerová, Z., Hron, T., Szabó, A., et al. (2017). Overexpression of TET dioxygenases in seminomas associates with low levels of DNA methylation and hydroxymethylation. Mol. Carcinog. 56, 1837–1850. doi:10.1002/mc.22638
Bowman, R. L., and Levine, R. L. (2017). TET2 in normal and malignant hematopoiesis. Cold Spring Harb. Perspect. Med. 7, a026518. doi:10.1101/cshperspect.a026518
Bray, J. K., Dawlaty, M. M., Verma, A., and Maitra, A. (2021). Roles and regulations of TET enzymes in solid tumors. Trends Cancer 7, 635–646. doi:10.1016/j.trecan.2020.12.011
Chaudhry, P., and Asselin, E. (2009). Resistance to chemotherapy and hormone therapy in endometrial cancer. Endocr. Relat. Cancer 16, 363–380. doi:10.1677/erc-08-0266
Chen, C., Wang, Z., Ding, Y., Wang, L., Wang, S., Wang, H., et al. (2022). DNA methylation: from cancer biology to clinical perspectives. Front. Biosci. Landmark Ed. 27, 326. doi:10.31083/j.fbl2712326
Chen, L. L., Lin, H. P., Zhou, W. J., He, C. X., Zhang, Z. Y., Cheng, Z. L., et al. (2018). SNIP1 recruits TET2 to regulate c-MYC target genes and cellular DNA damage response. Cell Rep. 25, 1485–1500. doi:10.1016/j.celrep.2018.10.028
Cimmino, L., Dolgalev, I., Wang, Y., Yoshimi, A., Martin, G. H., Wang, J., et al. (2017). Restoration of TET2 function blocks aberrant self-renewal and leukemia progression. Cell 170, 1079–1095. doi:10.1016/j.cell.2017.07.032
Delatte, B., Deplus, R., and Fuks, F. (2014). Playing TETris with DNA modifications. Embo J. 33, 1198–1211. doi:10.15252/embj.201488290
Esteller, M., and Herman, J. G. (2002). Cancer as an epigenetic disease: DNA methylation and chromatin alterations in human tumours. J. Pathol. 196, 1–7. doi:10.1002/path.1024
Golemis, E. A., Scheet, P., Beck, T. N., Scolnick, E. M., Hunter, D. J., Hawk, E., et al. (2018). Molecular mechanisms of the preventable causes of cancer in the United States. Genes Dev. 32, 868–902. doi:10.1101/gad.314849.118
Gonçalves, A. C., Richiardone, E., Jorge, J., Polónia, B., Xavier, C. P. R., Salaroglio, I. C., et al. (2021). Impact of cancer metabolism on therapy resistance - clinical implications. Drug Resist Updat 59, 100797. doi:10.1016/j.drup.2021.100797
Huang, G., Liu, L., Wang, H., Gou, M., Gong, P., Tian, C., et al. (2020). Tet1 deficiency leads to premature reproductive aging by reducing spermatogonia stem cells and germ cell differentiation. iScience 23, 100908. doi:10.1016/j.isci.2020.100908
Ismail, J. N., Ghannam, M., Al Outa, A., Frey, F., and Shirinian, M. (2020). Ten-eleven translocation proteins and their role beyond DNA demethylation - what we can learn from the fly. Epigenetics 15, 1139–1150. doi:10.1080/15592294.2020.1767323
Jiang, D., Wei, S., Chen, F., Zhang, Y., and Li, J. (2017). TET3-mediated DNA oxidation promotes ATR-dependent DNA damage response. EMBO Rep. 18, 781–796. doi:10.15252/embr.201643179
Joshi, K., Liu, S., Breslin, S. J. P., and Zhang, J. (2022). Mechanisms that regulate the activities of TET proteins. Cell Mol. Life Sci. 79, 363. doi:10.1007/s00018-022-04396-x
Joshi, K., Zhang, L., Breslin, S. J. P., Kini, A. R., and Zhang, J. (2022). Role of TET dioxygenases in the regulation of both normal and pathological hematopoiesis. J. Exp. Clin. Cancer Res. 41, 294. doi:10.1186/s13046-022-02496-x
Kao, S. H., Wu, K. J., and Lee, W. H. (2016). Hypoxia, epithelial-mesenchymal transition, and TET-mediated epigenetic changes. J. Clin. Med. 5, 24. doi:10.3390/jcm5020024
Kinney, S. R., and Pradhan, S. (2013). Ten eleven translocation enzymes and 5-hydroxymethylation in mammalian development and cancer. Adv. Exp. Med. Biol. 754, 57–79. doi:10.1007/978-1-4419-9967-2_3
Klutstein, M., Nejman, D., Greenfield, R., and Cedar, H. (2016). DNA methylation in cancer and aging. Cancer Res. 76, 3446–3450. doi:10.1158/0008-5472.CAN-15-3278
Kohli, R. M., and Zhang, Y. (2013). TET enzymes, TDG and the dynamics of DNA demethylation. Nature 502, 472–479. doi:10.1038/nature12750
Kossenas, K., and Constantinou, C. (2021). Epidemiology, molecular mechanisms, and clinical trials: an update on research on the association between red meat consumption and colorectal cancer. Curr. Nutr. Rep. 10, 435–467. doi:10.1007/s13668-021-00377-x
Kuhlen, M., Taeubner, J., Brozou, T., Wieczorek, D., Siebert, R., and Borkhardt, A. (2019). Family-based germline sequencing in children with cancer. Oncogene 38, 1367–1380. doi:10.1038/s41388-018-0520-9
Kulis, M., and Esteller, M. (2010). DNA methylation and cancer. Adv. Genet. 70, 27–56. doi:10.1016/b978-0-12-380866-0.60002-2
Kunimoto, H., and Nakajima, H. (2021). TET2: a cornerstone in normal and malignant hematopoiesis. Cancer Sci. 112, 31–40. doi:10.1111/cas.14688
Li, H., Jiang, W., Liu, X. N., Yuan, L. Y., Li, T. J., Li, S., et al. (2020). TET1 downregulates epithelial-mesenchymal transition and chemoresistance in PDAC by demethylating CHL1 to inhibit the Hedgehog signaling pathway. Oncogene 39, 5825–5838. doi:10.1038/s41388-020-01407-8
Liu, J., Jiang, J., Mo, J., Liu, D., Cao, D., Wang, H., et al. (2019). Global DNA 5-hydroxymethylcytosine and 5-formylcytosine contents are decreased in the early stage of hepatocellular carcinoma. Hepatology 69, 196–208. doi:10.1002/hep.30146
Ma, C., Seong, H., Liu, Y., Yu, X., Xu, S., and Li, Y. (2021). Ten-eleven translocation proteins (TETs): tumor suppressors or tumor enhancers? Front. Biosci. Landmark Ed. 26, 895–915. doi:10.52586/4996
Matuleviciute, R., Cunha, P. P., Johnson, R. S., and Foskolou, I. P. (2021). Oxygen regulation of TET enzymes. Febs J. 288, 7143–7161. doi:10.1111/febs.15695
Meng, H., Cao, Y., Qin, J., Song, X., Zhang, Q., Shi, Y., et al. (2015). DNA methylation, its mediators and genome integrity. Int. J. Biol. Sci. 11, 604–617. doi:10.7150/ijbs.11218
Miyamoto, K., and Ushijima, T. (2003). DNA methylation and cancer--DNA methylation as a target of cancer chemotherapy. Gan Kagaku Ryoho 30, 2021–2029.
Nishiyama, A., and Nakanishi, M. (2021). Navigating the DNA methylation landscape of cancer. Trends Genet. 37, 1012–1027. doi:10.1016/j.tig.2021.05.002
Nooter, K., and Stoter, G. (1996). Molecular mechanisms of multidrug resistance in cancer chemotherapy. Pathol. Res. Pract. 192, 768–780. doi:10.1016/s0344-0338(96)80099-9
Pan, F., Weeks, O., Yang, F. C., and Xu, M. (2015). The TET2 interactors and their links to hematological malignancies. IUBMB Life 67, 438–445. doi:10.1002/iub.1389
Pan, Y., Liu, G., Zhou, F., Su, B., and Li, Y. (2018). DNA methylation profiles in cancer diagnosis and therapeutics. Clin. Exp. Med. 18, 1–14. doi:10.1007/s10238-017-0467-0
Panczyk, M. (2014). Pharmacogenetics research on chemotherapy resistance in colorectal cancer over the last 20 years. World J. Gastroenterol. 20, 9775–9827. doi:10.3748/wjg.v20.i29.9775
Panigrahi, G. K., Yadav, A., Srivastava, A., Tripathi, A., Raisuddin, S., and Das, M. (2015). Mechanism of rhein-induced apoptosis in rat primary hepatocytes: beneficial effect of cyclosporine A. Chem. Res. Toxicol. 28, 1133–1143. doi:10.1021/acs.chemrestox.5b00063
Papanicolau-Sengos, A., and Aldape, K. (2022). DNA methylation profiling: an emerging paradigm for cancer diagnosis. Annu. Rev. Pathol. 17, 295–321. doi:10.1146/annurev-pathol-042220-022304
Rasmussen, K. D., and Helin, K. (2016). Role of TET enzymes in DNA methylation, development, and cancer. Genes Dev. 30, 733–750. doi:10.1101/gad.276568.115
Rodriguez, R., Schreiber, S. L., and Conrad, M. (2022). Persister cancer cells: iron addiction and vulnerability to ferroptosis. Mol. Cell 82, 728–740. doi:10.1016/j.molcel.2021.12.001
Sajadian, S. O., Ehnert, S., Vakilian, H., Koutsouraki, E., Damm, G., Seehofer, D., et al. (2015). Induction of active demethylation and 5hmC formation by 5-azacytidine is TET2 dependent and suggests new treatment strategies against hepatocellular carcinoma. Clin. Epigenetics 7, 98. doi:10.1186/s13148-015-0133-x
Scott-Browne, J. P., Lio, C. J., and Rao, A. (2017). TET proteins in natural and induced differentiation. Curr. Opin. Genet. Dev. 46, 202–208. doi:10.1016/j.gde.2017.07.011
Shekhawat, J., Gauba, K., Gupta, S., Choudhury, B., Purohit, P., Sharma, P., et al. (2021). Ten-eleven translocase: key regulator of the methylation landscape in cancer. J. Cancer Res. Clin. Oncol. 147, 1869–1879. doi:10.1007/s00432-021-03641-3
Smeets, E., Lynch, A. G., Prekovic, S., Van den Broeck, T., Moris, L., Helsen, C., et al. (2018). The role of TET-mediated DNA hydroxymethylation in prostate cancer. Mol. Cell Endocrinol. 462, 41–55. doi:10.1016/j.mce.2017.08.021
Srivastava, S. P., and Goodwin, J. E. (2020). Cancer biology and prevention in diabetes. Cells 9, 1380. doi:10.3390/cells9061380
Talib, W. H., Mahmod, A. I., Kamal, A., Rashid, H. M., Alashqar, A. M. D., Khater, S., et al. (2021). Ketogenic diet in cancer prevention and therapy: molecular targets and therapeutic opportunities. Curr. Issues Mol. Biol. 43, 558–589. doi:10.3390/cimb43020042
Tyagi, N., Roy, S., Vengadesan, K., and Gupta, D. (2024). Multi-omics approach for identifying CNV-associated lncRNA signatures with prognostic value in prostate cancer. Noncoding RNA Res. 9, 66–75. doi:10.1016/j.ncrna.2023.10.001
Xie, B. Y., Lv, Q. Y., Ning, C. C., Yang, B. Y., Shan, W. W., Cheng, Y. L., et al. (2017). TET1-GPER-PI3K/AKT pathway is involved in insulin-driven endometrial cancer cell proliferation. Biochem. Biophys. Res. Commun. 482, 857–862. doi:10.1016/j.bbrc.2016.11.124
Zaridze, D. G. (2008). Molecular epidemiology of cancer. Biochem. (Mosc) 73, 532–542. doi:10.1134/s0006297908050064
Keywords: TET family genes, pan-cancer analysis, tumor microenvironment, drug sensitivity, therapeutic target
Citation: Zhang C, Zheng J, Liu J, Li Y, Xing G, Zhang S, Chen H, Wang J, Shao Z, Li Y, Jiang Z, Pan Y, Liu X, Xu P and Wu W (2024) Pan-cancer analyses reveal the molecular and clinical characteristics of TET family members and suggests that TET3 maybe a potential therapeutic target. Front. Pharmacol. 15:1418456. doi: 10.3389/fphar.2024.1418456
Received: 16 April 2024; Accepted: 28 May 2024;
Published: 22 July 2024.
Edited by:
Gregory Wiedman, Seton Hall University, United StatesReviewed by:
Fatemeh Saheb Sharif-Askari, University of Sharjah, United Arab EmiratesAdrian Bogdan Tigu, University of Medicine and Pharmacy Iuliu Hatieganu, Romania
Copyright © 2024 Zhang, Zheng, Liu, Li, Xing, Zhang, Chen, Wang, Shao, Li, Jiang, Pan, Liu, Xu and Wu. 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: Ping Xu, zishanshuxu@sina.com; Wenhan Wu, wuwenhan88@126.com