DNMT/TET genes expression profiles in pan-cancer
We summarized the chromosomal positions of all DNMT/TET family members across the published literature (Table 1). We then analyzed the gene expression differences in TCGA. The results showed heterogeneous expression of DNTM and TET genes in different cancer types. DNMT3B was highly expressed in the majority of tumors; DNMT3L was expressed in low levels in CHOL, KICH, KIRP, and KIRC, and had increased expression in UCEC READ, ESCA, and LUSC. (Figure 1a). The Log2FC and adjusted P-value of the results are listed in Supplementary Table 1. We visualized the differential expression of DNMT3B in different cancers (Figure 1b). DNMT3B expression was up regulated in many cancer types, including BLCA, BRCA, CESC, CHOL, COAD, ESCA, GBM, HNSC, KIRC, KIRP, LIHC, LUAD, LUSC, PRAD, READ, STAD and UCEC.
Table 1 Basic characteristics of DNMT/TET genes
HGNC ID (gene)
|
GENE ID
|
Approved symbol
|
Synonym(s)
|
Exon
|
Chromosomal location
|
2976
|
1786
|
DNMT1
|
AIM
|
41
|
19p13.2
|
DNMT
|
MCMT
|
CXXC9
|
HSN1E
|
ADCADN
|
m.HsaI
|
2977
|
1787
|
TRDMT1
|
DMNT2
|
15
|
10p13
|
DNMT2
|
PUMET
|
RNMT1
|
MHSAIIP
|
2978
|
1788
|
DNMT3A
|
TBRS
|
34
|
2p23.3
|
HESJAS
|
DNMT3A2
|
M.HsaIIIA
|
2979
|
1789
|
DNMT3B
|
ICF
|
24
|
20q11.21
|
ICF1
|
M.HsaIIIB
|
2980
|
29947
|
DNMT3L
|
-
|
12
|
21q22.3
|
18291
|
55929
|
DMAP1
|
EAF2
|
11
|
1p34.1
|
SWC4
|
MEAF2
|
DNMAP1
|
DNMTAP1
|
29484
|
80312
|
TET1
|
LCX
|
20
|
10q21.3
|
CXXC6
|
bA119F7.1
|
25941
|
54790
|
TET2
|
MDS
|
15
|
4q24
|
KIAA1546
|
28313
|
200424
|
TET3
|
BEFAHRS
|
17
|
2p13.1
|
hCG_40738
|
Furthermore, we analyzed DNMT/TET protein expression levels in different cancers types. In COAD, STCA, LIHC and PACA, DNMT3A, TRDMT1, and TET2 were high expressed, while TET3 expression was significantly reduced (Figure2a). Detailed DNMT/TET protein expression information in 21 cancers is shown in Supplementary Table 2. DNMT3B protein expression in different human cancers is shown in Figure 2b.
Additionally, we elucidated the correlations between DNMT and TET genes expression in pan-cancer. Significant correlations were identified between DNMT3A-TET1 (r = 0.57, p<0.05), TET2-TET3 (r = 0.55, p<0.05), and TET1-TET3 (r = 0.50, p<0.05) (Figure3).
DNMT/TET genetic alterations in pan-cancer
We examined the DNMT/TET gene mutation frequency in 33 cancer types. Most DNMT/TET genes were frequently mutated in UCEC, and rare mutations were observed in CHOL, PCPG, and THCA. TET1 and TET2 had higher mutation frequencies than did other genes in pan-cancer, and the overall average mutation frequency of DNMT/TET genes in different cancer types ranged from 0% to 34.9% (Figure 4a). We also investigated DNMT/TET gene CNVs in pan-cancer (Figure 4b). DNMT3B exhibited inclusive copy number amplification in STAD, COAD, READ, BLCA, ESCA, and LUSC, but almost no copy number deletions in multiple cancer types. Most DNMT/TET genes displayed copy number deletions in KICH, but copy number variations were rarely observed in LIHC and LAML.
Next, we analyzed the relationship between DNMT/TET gene mutations and expression levels in human cancers and visualized the statistically significant results (Figure 5). DNMT1 expression level changes were related to its mutation in COAD, ESCA, STAD, UCEC, and PRAD. Similarly, mutations in DNMT3L and TET2 also affected their expression levels in COAD. The details were shown in Supplementary Table 3. The association between CNVs and DNMT/TET gene expression levels in multiple cancer types is shown in Figure 6. We found that almost all expression changes in DNMT/TET genes were associated with CNVs.
The correlation between DNMT/TET genes and cancer-related pathways
We investigated the correlation between DNMT/TET gene expression and cancer related pathways to clarify the molecular significance of DNMT/TET genes in carcinogenesis. DNMT/TET genes were significantly correlated with multiple carcinogenic pathways, playing both activation and suppression roles (Figure 7a). We concluded that DNMT/TET genes are mainly involved in the following cancer-related pathways: UV response DN, mitotic spindle, cholesterol homeostasis, TGF beta signaling, xenobiotic metabolism, G2/M checkpoint, and E2F targets (Supplementary Table4). We also found that DNMT3A, TET1, TET2, and TET3 were more likely to be related to the tumor occurrence and progression pathways (Figure 7b).
The correlation between DNMT/TET genes and immune infiltration
The correlation between DNMT/TET gene expression and immune cell infiltration in pan-cancer were assessed. DNMT/TET genes were significantly correlated with NK cells, CD4 positive T cells, and Tfh cells. Notably, the correlation between DNMT/TET genes and immune cell infiltration tended to be positive. In particular, TET2, TRDMT1, and TET3 showed more stronger correlation. (Figure 8)
Prognostic value of DNMT/TET genes in pan-cancer
The association between DNMT/TET gene expression and prognosis in patients with different cancers is shown in Figure 9a. Most DNMT/TET genes were significantly associated with lower survival rates in ACC, MESO, and LIHC, but played a protective role in THYM. DNMT3B was related to poor prognosis in many cancers, including ACC, KIRC, KIRP, LGG, LIHC, MESO, SARC, and UCEC. TRDMT1 and DNMT3A had clearly different effects on prognosis in various cancer types (Figure 9b). In ACC, KICH, LIHC, PRAD, and SARC, TRDMT1 was a predictor of poor prognosis but not in COAD, GBM, KIRC, PCPG, and READ. Meanwhile, DNMT3A expression was associated with poor prognosis in ACC, KICH, LGG, LIHC, and MESO, but not in LAML, PAAD, THYM, and UCS. Moreover, high or low DNMT gene expression was statistically significantly related to patient survival status (Figure 10). Increased DNMT gene expression, except for DMAP1, was associated worse prognoses in pan-cancer. Compared to other genes, DNMT3B had the worst prognosis in cancer patients. No significant correlation was observed between the expression levels of TET1 or TET2 and overall cancer survival.